Identifying homologous recombination deficiency in breast cancer: genomic instability score distributions differ among breast cancer subtypes.

PURPOSE: A 3-biomarker homologous recombination deficiency (HRD) score is a key component of a currently FDA-approved companion diagnostic assay to identify HRD in patients with ovarian cancer using a threshold score of ≥ 42, though recent studies have explored the utility of a lower threshold (GIS ≥ 33). The present study evaluated whether the ovarian cancer thresholds may also be appropriate for major breast cancer subtypes by comparing the genomic instability score (GIS) distributions of BRCA1/2-deficient estrogen receptor-positive breast cancer (ER + BC) and triple-negative breast cancer (TNBC) to the GIS distribution of BRCA1/2-deficient ovarian cancer. METHODS: Ovarian cancer and breast cancer (ER + BC and TNBC) tumors from ten study cohorts were sequenced to identify pathogenic BRCA1/2 mutations, and GIS was calculated using a previously described algorithm. Pathologic complete response (pCR) to platinum therapy was evaluated in a subset of TNBC samples. For TNBC, a threshold was set and threshold validity was assessed relative to clinical outcomes. RESULTS: A total of 560 ovarian cancer, 805 ER + BC, and 443 TNBC tumors were included. Compared to ovarian cancer, the GIS distribution of BRCA1/2-deficient samples was shifted lower for ER + BC (p = 0.015), but not TNBC (p = 0.35). In the subset of TNBC samples, univariable logistic regression models revealed that GIS status using thresholds of ≥ 42 and ≥ 33 were significant predictors of response to platinum therapy. CONCLUSIONS: This study demonstrated that the GIS thresholds used for ovarian cancer may also be appropriate for TNBC, but not ER + BC. GIS thresholds in TNBC were validated using clinical response data to platinum therapy.

Clinical significance of homologous recombination deficiency score testing in endometrial Cancer.

BACKGROUND: Homologous recombination deficiency (HRD) score is related to chemotherapy response in some cancers, but its role in endometrial cancer in not known. We determined frequency and clinical significance of alterations in the HR pathway in endometrial cancer. METHODS: 253 endometrioid endometrial adenocarcinoma (EEA) samples from two independent cohorts (discovery and replication) were tested for HRD score using the Myriad HRD assay, microsatellite instability (MSI) and tumor mutation burden (TMB) using a next generation sequencing assay. HRD scores were also generated on endometrial cancer cell lines and in vivo response to olaparib was assessed. RESULTS: ROC curves were employed to determine optimal cutoffs of HRD in relation to survival impact in endometrial cancer and a cutoff of HRD ≥ 4 was suggested for DFS using the discovery cohort. Patients from two independent cohorts with HRD score ≥ 4 trended toward worse survival as compared to those with HRD score < 4. Both cohorts were further separated into four groups according to molecular subtypes (TMB positive; MSI positive; HRD positive; all others). When grouped by molecular subtype, there was a significant difference between groups using an HRD ≥4 cutoff in the initial (p = 0.0024) and replication (p = 0.042) cohorts. The Hec1a model (HRD score = 19) was highly sensitive to olaparib in in vitro and in vivo experiments. CONCLUSIONS: High HRD score was associated with worse DFS in our patient cohort. These findings suggest that HRD score may have clinical utility in patients with advanced or recurrent endometrial cancer.

Homologous recombination deficiency (HRD) status predicts response to standard neoadjuvant chemotherapy in patients with triple-negative or BRCA1/2 mutation-associated breast cancer.

PURPOSE: Defects in the homologous recombination (HR) DNA repair pathway sensitize tumors to therapeutics that target this pathway. A significant proportion of triple-negative breast cancers (TNBC) carry HR defects. The HRD assay is highly associated with sensitivity to neoadjuvant platinum-based chemotherapy in TNBC. Standard chemotherapy consists of some combination of an anthracycline, cyclophosphamide, and taxane. This study assesses the association of HR deficiency status with response to standard neoadjuvant chemotherapy in TNBC or BRCA1/2 mutation-associated breast cancer. METHODS: Tumor samples were retrospectively obtained from 45 TNBC patients and 2 BRCA1/2 mutant, hormone receptor-positive/HER2-negative breast cancer patients who received anthracycline- and/or taxane-based neoadjuvant chemotherapy at Stanford University or Cedars-Sinai Medical Centers. The HRD score and tumor BRCA1/2 mutation status were determined from baseline tumor biopsies. HR deficient tumors were those with a HRD score of ≥ 42 or a tumor BRCA1/2 mutation. Response was categorized by the residual cancer burden (RCB) index. RESULTS: HR deficient patients were more likely to achieve a pathologic complete response (pCR) compared with non-deficient patients (OR 13.06, CI 1.52-11.241, p = 0.0028). Among BRCA1/2 mutation wild-type patients, HR deficient patients were more likely to achieve a pCR (OR 16, 95% CI 1.65-160.41, p = 0.0041) compared with HR non-deficient patients. Further, HRD scores were highly concordant pre- and post-therapy (Spearman correlation > 99%). CONCLUSIONS: HR deficiency status is significantly associated with response to standard neoadjuvant chemotherapy in TNBC. This observation is consistent with the mechanisms of action of doxorubicin and cyclophosphamide as DNA damaging agents.

Association of BRCA1/2 defects with genomic scores predictive of DNA damage repair deficiency among breast cancer subtypes.

INTRODUCTION: Homologous recombination (HR) DNA repair is of clinical relevance in breast cancer. Three DNA-based homologous recombination deficiency (HRD) scores (HRD-loss of heterozygosity score (LOH), HRD-telomeric allelic imbalance score (TAI), and HRD-large-scale state transition score (LST)) have been developed that are highly correlated with defects in BRCA1/2, and are associated with response to platinum therapy in triple negative breast and ovarian cancer. This study examines the frequency of BRCA1/2 defects among different breast cancer subtypes, and the ability of the HRD scores to identify breast tumors with defects in the homologous recombination DNA repair pathway. METHODS: 215 breast tumors representing all ER/HER2 subtypes were obtained from commercial vendors. Next-generation sequencing based assays were used to generate genome wide SNP profiles, BRCA1/2 mutation screening, and BRCA1 promoter methylation data. RESULTS: BRCA1/2 deleterious mutations were observed in all breast cancer subtypes. BRCA1 promoter methylation was observed almost exclusively in triple negative breast cancer. BRCA1/2 deficient tumors were identified with BRCA1/2 mutations, or BRCA1 promoter methylation, and loss of the second allele of the affected gene. All three HRD scores were highly associated with BRCA1/2 deficiency (HRD-LOH: P = 1.3 × 10(-17); HRD-TAI: P = 1.5 × 10(-19); HRD-LST: P = 3.5 × 10(-18)). A combined score (HRD-mean) was calculated using the arithmetic mean of the three scores. In multivariable analyses the HRD-mean score captured significant BRCA1/2 deficiency information not captured by the three individual scores, or by clinical variables (P values for HRD-Mean adjusted for HRD-LOH: P = 1.4 × 10(-8); HRD-TAI: P = 2.9 × 10(-7); HRD-LST: P = 2.8 × 10(-8); clinical variables: P = 1.2 × 10(-16)). CONCLUSIONS: The HRD scores showed strong correlation with BRCA1/2 deficiency regardless of breast cancer subtype. The frequency of elevated scores suggests that a significant proportion of all breast tumor subtypes may carry defects in the homologous recombination DNA repair pathway. The HRD scores can be combined to produce a more robust predictor of HRD. The combination of a robust score, and the FFPE compatible assay described in this study, may facilitate use of agents targeting homologous recombination DNA repair in the clinical setting.

Proteomic screening of variola virus reveals a unique NF-kappaB inhibitor that is highly conserved among pathogenic orthopoxviruses.

Identification of the binary interactions between viral and host proteins has become a valuable tool for investigating viral tropism and pathogenesis. Here, we present the first systematic protein interaction screening of the unique variola virus proteome by using yeast 2-hybrid screening against a variety of human cDNA libraries. Several protein-protein interactions were identified, including an interaction between variola G1R, an ankryin/F-box containing protein, and human nuclear factor kappa-B1 (NF-kappaB1)/p105. This represents the first direct interaction between a pathogen-encoded protein and NF-kappaB1/p105. Orthologs of G1R are present in a variety of pathogenic orthopoxviruses, but not in vaccinia virus, and expression of any one of these viral proteins blocks NF-kappaB signaling in human cells. Thus, proteomic screening of variola virus has the potential to uncover modulators of the human innate antiviral responses.

Homologous Recombination Deficiency as an Ovarian Cancer Biomarker in a Real-World Cohort: Validation of Decentralized Genomic Profiling.

The diagnostic evaluation of homologous recombination deficiency (HRD) is central to define targeted therapy strategies for patients with ovarian carcinoma. We evaluated HRD in 514 ovarian carcinoma samples by next-generation sequencing of DNA libraries, including BRCA1/BRCA2 and 26,523 single-nucleotide polymorphisms using the standardized Myriad HRD assay, with the predefined cut point of ≥42 for a positive genomic instability score (GIS). All samples were measured in the central Myriad laboratory and in an academic molecular pathology laboratory. A positive GIS was detected in 196 (38.1%) of tumors, whereas 318 (61.9%) were GIS negative. Combining GIS and BRCA mutations, a total of 200 (38.9%) of the 514 tumors were HRD positive. A positive GIS was significantly associated with high-grade serous histology (P < 0.000001), grade 3 tumors (P = 0.001), and patient age <60 years (P = 0.0003). The concordance between both laboratories for the GIS status was 96.9% (P < 0.000001), with a sensitivity of 94.6% and a specificity of 98.4%. Concordance for HRD status was 97.1% (499 of 514 tumors). The percentage of HRD-positive tumors in our real-life cohort was similar to the proportion observed in the recently published PAOLA-1 trial, with high concordance between central and local laboratories. Our results support introduction of the standardized HRD assay in academic molecular pathology laboratories, thus broadening access to personalized oncology strategies for patients with ovarian cancer worldwide.

The myxoma virus m-t5 ankyrin repeat host range protein is a novel adaptor that coordinately links the cellular signaling pathways mediated by Akt and Skp1 in virus-infected cells.

Most poxviruses express multiple proteins containing ankyrin (ANK) repeats accounting for a large superfamily of related but unique determinants of poxviral tropism. Recently, select members of this novel family of poxvirus proteins have drawn considerable attention for their potential roles in modulating intracellular signaling networks during viral infection. The rabbit-specific poxvirus, myxoma virus (MYXV), encodes four unique ANK repeat proteins, termed M-T5, M148, M149, and M150, all of which include a carboxy-terminal PRANC domain which closely resembles a cellular protein motif called the F-box domain. Here, we show that each MYXV-encoded ANK repeat protein, including M-T5, interacts directly with the Skp1 component of the host SCF ubiquitin ligase complex, and that the binding of M-T5 to cullin 1 is indirect via binding to Skp1 in the host SCF complex. To understand the significance of these virus-host protein interactions, the various binding domains of M-T5 were mapped. The N-terminal ANK repeats I and II were identified as being important for interaction with Akt, whereas the C-terminal PRANC/F-box-like domain was essential for binding to Skp1. We also report that M-T5 can bind Akt and the host SCF complex (via Skp1) simultaneously in MYXV-infected cells. Finally, we report that M-T5 specifically mediates the relocalization of Akt from the nucleus to the cytoplasm during infection with the wild-type MYXV, but not the M-T5 knockout version of the virus. These results indicate that ANK/PRANC proteins play a critical role in reprogramming disparate cellular signaling cascades to establish a new cellular environment more favorable for virus replication.

Evidence for epistasis between SLC6A4 and a chromosome 4 gene as risk factors in major depression.

Linkage analysis on Utah pedigrees with strong family histories of major depression including only cases with the SLC6A4 HTTLPR short allele revealed a linkage peak on chromosome 4 (maximum HLOD = 3.5). This evidence suggests epistasis between SLC6A4 and an unknown gene as risk factors for major depression.

Association of Tumor-Infiltrating Lymphocytes with Homologous Recombination Deficiency and BRCA1/2 Status in Patients with Early Triple-Negative Breast Cancer: A Pooled Analysis.

PURPOSE: Patients with triple-negative breast cancer (TNBC) with homologous recombination deficient tumors achieve significantly higher pathologic complete response (pCR) rates when treated with neoadjuvant platinum-based therapy. Tumor-infiltrating lymphocytes (TIL) are prognostic and predictive of chemotherapy benefit in early stage TNBC. The relationship between TILs, BRCA1/2 mutation status, and homologous recombination deficiency (HRD) status in TNBC remains unclear. EXPERIMENTAL DESIGN: We performed a pooled analysis of five phase II studies that included patients with TNBC treated with neoadjuvant platinum-based chemotherapy to evaluate the association of TILs with HRD status (Myriad Genetics) and tumor BRCA1/2 mutation status. Furthermore, the relationship between pathologic response assessed using the residual cancer burden (RCB) index and HRD status with adjustment for TILs was evaluated. RESULTS: Among 161 patients, stromal TIL (sTIL) density was not significantly associated with HRD status (P = 0.107) or tumor BRCA1/2 mutation status (P = 0.391). In multivariate analyses, sTIL density [OR, 1.23; 95% confidence interval (CI), 0.94-1.61; P = 0.139] was not associated with pCR, but was associated with RCB 0/I status (OR 1.62; 95% CI, 1.20-2.28; P = 0.001). HRD was significantly associated with both pCR (OR 12.09; 95% CI, 4.11-44.29; P = 7.82 × 10-7) and RCB 0/I (OR 10.22; 95% CI, 4.11-28.75; P = 1.09 × 10-7) in these models. CONCLUSIONS: In patients with TNBC treated with neoadjuvant platinum-based therapy, TIL density was not significantly associated with either tumor BRCA1/2 mutation status or HRD status. In this pooled analysis, HRD and sTIL density were independently associated with treatment response, with HRD status being the strongest predictor.

Intratumor Heterogeneity of Homologous Recombination Deficiency in Primary Breast Cancer.

Purpose: The 3-biomarker homologous recombination deficiency (HRD) assay measures the number of telomeric allelic imbalances, loss of heterozygosity, and large-scale state transitions in tumor DNA and combines these metrics into a single score that reflects DNA repair deficiency. The goal of this study is to assess the consistency of these HRD measures in different biopsies from distinct areas of the same cancer.Experimental Design: HRD scores, BRCA mutation status, and BRCA1 promoter methylation were assessed in 99 samples from 33 surgically resected, stage I-III breast cancers; each cancer was biopsied in three distinct areas. Homologous recombination repair (HR) deficiency was defined as either high HRD score (≥42) or tumor BRCA mutation.Results: Eighty-one biopsies from 32 cancers were analyzed. Tumor BRCA status was available for all samples, HRD scores for 70, and BRCA1 methylation values for 76 samples. The BRCA1/2 mutation and promoter methylation status and HR category showed perfect concordance across all biopsies from the same cancer. All tumors with BRCA1/2 mutations or promoter methylation had high HRD scores, as did 17% (4/24) of the BRCA1/2 wild-type and nonmethylated tumors. The HRD scores were also highly consistent between different biopsies from the same tumor with an intraclass correlation coefficient of 0.977, indicating that only 2.3% of the variance is attributed to within-tumor biopsy-to-biopsy variation.Conclusions: These results indicate that within-tumor spatial heterogeneity for HRD metrics and the technical noise in the assay are small and do not influence HRD scores and HR status. Clin Cancer Res; 23(5); 1193-9. ©2016 AACR.

Homologous Recombination Deficiency (HRD) Score Predicts Response to Platinum-Containing Neoadjuvant Chemotherapy in Patients with Triple-Negative Breast Cancer.

PURPOSE: BRCA1/2-mutated and some sporadic triple-negative breast cancers (TNBC) have DNA repair defects and are sensitive to DNA-damaging therapeutics. Recently, three independent DNA-based measures of genomic instability were developed on the basis of loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale state transitions (LST). EXPERIMENTAL DESIGN: We assessed a combined homologous recombination deficiency (HRD) score, an unweighted sum of LOH, TAI, and LST scores, in three neoadjuvant TNBC trials of platinum-containing therapy. We then tested the association of HR deficiency, defined as HRD score ≥42 or BRCA1/2 mutation, with response to platinum-based therapy. RESULTS: In a trial of neoadjuvant platinum, gemcitabine, and iniparib, HR deficiency predicted residual cancer burden score of 0 or I (RCB 0/I) and pathologic complete response (pCR; OR = 4.96, P = 0.0036; OR = 6.52, P = 0.0058). HR deficiency remained a significant predictor of RCB 0/I when adjusted for clinical variables (OR = 5.86, P = 0.012). In two other trials of neoadjuvant cisplatin therapy, HR deficiency predicted RCB 0/I and pCR (OR = 10.18, P = 0.0011; OR = 17.00, P = 0.0066). In a multivariable model of RCB 0/I, HR deficiency retained significance when clinical variables were included (OR = 12.08, P = 0.0017). When restricted to BRCA1/2 nonmutated tumors, response was higher in patients with high HRD scores: RCB 0/I P = 0.062, pCR P = 0.063 in the neoadjuvant platinum, gemcitabine, and iniparib trial; RCB 0/I P = 0.0039, pCR P = 0.018 in the neoadjuvant cisplatin trials. CONCLUSIONS: HR deficiency identifies TNBC tumors, including BRCA1/2 nonmutated tumors more likely to respond to platinum-containing therapy. Clin Cancer Res; 22(15); 3764-73. ©2016 AACR.

A virtual reality system for neurobehavioral and functional MRI studies.

We are developing a VR system of integrated software and hardware for scientific research and clinical application. The system is sufficiently flexible and broad-based in appeal that neurobehavioral researchers from a variety of disciplines might be interested in using it for basic research and clinical studies. The system runs on a standard Windows-based personal computer with a high-performance graphics card. Options allow a head-mounted display, dataglove, simultaneous physiological monitoring or use within neuroimaging machines such as magnetic resonance imaging (MRI) scanners. Currently, the software consists of a virtual world of nearly a dozen interconnected environments that the subject can freely navigate. Additional environments can be built and easily added to the application. A startup interface provides menus for selecting characters and objects that a researcher might want to put at specific locations within the simulation. Interactivity is provided for many typical objects such as doors, chairs and money. There are more than 50 characters in the world, most of them animated or interactive. All movements and actions of the subject within the world are tracked and recorded to an Excel spreadsheet for data analysis. Overlay maps are available as navigational aids. Concurrent physiological data can be acquired on up to 16 channels. The system provides synchronization of the VR simulation with physiological recordings and functional MR images. A spatial navigation memory task was performed with the integrated VR/fMRI system, and some pilot data is presented that shows robust activation in multiple cortical areas appropriate to the task.

The human-bacterial pathogen protein interaction networks of Bacillus anthracis, Francisella tularensis, and Yersinia pestis.

BACKGROUND: Bacillus anthracis, Francisella tularensis, and Yersinia pestis are bacterial pathogens that can cause anthrax, lethal acute pneumonic disease, and bubonic plague, respectively, and are listed as NIAID Category A priority pathogens for possible use as biological weapons. However, the interactions between human proteins and proteins in these bacteria remain poorly characterized leading to an incomplete understanding of their pathogenesis and mechanisms of immune evasion. METHODOLOGY: In this study, we used a high-throughput yeast two-hybrid assay to identify physical interactions between human proteins and proteins from each of these three pathogens. From more than 250,000 screens performed, we identified 3,073 human-B. anthracis, 1,383 human-F. tularensis, and 4,059 human-Y. pestis protein-protein interactions including interactions involving 304 B. anthracis, 52 F. tularensis, and 330 Y. pestis proteins that are uncharacterized. Computational analysis revealed that pathogen proteins preferentially interact with human proteins that are hubs and bottlenecks in the human PPI network. In addition, we computed modules of human-pathogen PPIs that are conserved amongst the three networks. Functionally, such conserved modules reveal commonalities between how the different pathogens interact with crucial host pathways involved in inflammation and immunity. SIGNIFICANCE: These data constitute the first extensive protein interaction networks constructed for bacterial pathogens and their human hosts. This study provides novel insights into host-pathogen interactions.

Analysis of vaccinia virus-host protein-protein interactions: validations of yeast two-hybrid screenings.

Vaccinia virus, a large double-stranded DNA virus, is the prototype of the Orthopoxvirus genus, which includes several pathogenic poxviruses of humans, such as monkeypox virus and variola virus. Here, we report a comprehensive yeast two-hybrid (Y2H) screening for the protein-protein interactions between vaccinia and human proteins. A total of 109 novel vaccinia-human protein interactions were detected among 33 viral proteins. To validate subsets of those interactions, we constructed an ORFeome library of vaccinia virus strain WR using the Gateway plasmid cloning system. By co-expressing selected vaccinia and host proteins in a variety of expression systems, we found that at least 17 of the Y2H hits identified between vaccinia and human proteins can be verified by independent methods using GST pull-down assays, representing a 63% validation rate for the Y2H hits examined (17/27). Because the cloned ORFs are conveniently transferable from the entry vectors to various destination expression vectors, the vaccinia ORFeome library will be a useful resource for future high-throughput functional proteomic experiments.

TBC1D1 is a candidate for a severe obesity gene and evidence for a gene/gene interaction in obesity predisposition.

The molecular etiology of obesity predisposition is largely unknown. Here, we present evidence that genetic variation in TBC1D1 confers risk for severe obesity in females. We identified a coding variant (R125W) in TBC1D1 that segregated with the disease in 4p15-14-linked obesity pedigrees. In cases derived from pedigrees with the strongest linkage evidence, the variant was significantly associated with obesity (P=0.000007) and chromosomes carrying R125W accounted for the majority of the evidence that originally linked 4p15-14 with the disease. In addition, by selecting families that segregated R125W with obesity, we were able to generate highly significant linkage evidence for an obesity predisposition locus at 4q34-35. This result provides additional and confirming evidence that R125W affects obesity susceptibility, delimits the location of an obesity gene at 4q34-35 and identifies a gene/gene interaction that influences the risk for obesity predisposition. Finally, although the function of TBC1D1 is unknown, the protein is structurally similar to a known regulator of insulin-mediated Glut4 translocation.

A major predisposition locus for severe obesity, at 4p15-p14.

Although the predisposition to morbid obesity is heritable, the identities of the disease-causing genes are largely unknown. Therefore, we have conducted a genomewide search with 628 markers, using multigenerational Utah pedigrees to identify genes involved in predisposition to obesity. In the genomewide search, we identified a highly significant linkage to high body-mass index in female patients, at D4S2632, with a multipoint heterogeneity LOD (HLOD) score of 6.1 and a nonparametric linkage (NPL) score of 5.3. To further delineate the linkage, we increased both the marker density around D4S2632 and the size of our pedigree data set. As a result, the linkage evidence increased to a multipoint HLOD score of 9.2 (at D4S3350) and an NPL score of 11.3. Evidence from almost half of the families in this analysis support this linkage, and therefore the gene in this region might account for a significant percentage of the genetic predisposition to severe obesity in females. However, further studies are necessary to clarify the effect that this gene has in males and in the general population.

Predisposition locus for major depression at chromosome 12q22-12q23.2.

Major depression disorder is a common psychiatric disease with a major economic impact on society. In many cases, no effective treatment is available. The etiology of major depression is complex, but it is clear that the disease is, to a large extent, determined genetically, especially among individuals with a familial history of major depression, presumably through the involvement of multiple predisposition genes in addition to an environmental component. As a first step toward identification of chromosomal loci contributing to genetic predisposition to major depression, we have conducted a genomewide scan by using 628 microsatellite markers on 1,890 individuals from 110 Utah pedigrees with a strong family history of major depression. We identified significant linkage to major depression in males at marker D12S1300 (multipoint heterogeneity LOD score 4.6; P=.00003 after adjustment for multiple testing). With additional markers, the linkage evidence became highly significant, with the multipoint heterogeneity LOD score at marker D12S1706 increasing to 6.1 (P=.0000007 after adjustment for multiple testing). This study confirms the presence of one or more genes involved in psychiatric diseases on the q arm of chromosome 12 and provides strong evidence for the existence of a sex-specific predisposition gene to major depression at 12q22-q23.2.