Germline mutations in BRCA1/2 predispose individuals to breast cancer (termed germline-mutated BRCA1/2 breast cancer, gBRCA-BC) by impairing homologous recombination (HR) and causing genomic ...instability. HR also repairs DNA lesions caused by platinum agents and PARP inhibitors. Triple-negative breast cancers (TNBCs) harbor subpopulations with BRCA1/2 mutations, hypothesized to be especially platinum-sensitive. Cancers in putative 'BRCAness' subgroups-tumors with BRCA1 methylation; low levels of BRCA1 mRNA (BRCA1 mRNA-low); or mutational signatures for HR deficiency and those with basal phenotypes-may also be sensitive to platinum. We assessed the efficacy of carboplatin and another mechanistically distinct therapy, docetaxel, in a phase 3 trial in subjects with unselected advanced TNBC. A prespecified protocol enabled biomarker-treatment interaction analyses in gBRCA-BC and BRCAness subgroups. The primary endpoint was objective response rate (ORR). In the unselected population (376 subjects; 188 carboplatin, 188 docetaxel), carboplatin was not more active than docetaxel (ORR, 31.4% versus 34.0%, respectively; P = 0.66). In contrast, in subjects with gBRCA-BC, carboplatin had double the ORR of docetaxel (68% versus 33%, respectively; biomarker, treatment interaction P = 0.01). Such benefit was not observed for subjects with BRCA1 methylation, BRCA1 mRNA-low tumors or a high score in a Myriad HRD assay. Significant interaction between treatment and the basal-like subtype was driven by high docetaxel response in the nonbasal subgroup. We conclude that patients with advanced TNBC benefit from characterization of BRCA1/2 mutations, but not BRCA1 methylation or Myriad HRD analyses, to inform choices on platinum-based chemotherapy. Additionally, gene expression analysis of basal-like cancers may also influence treatment selection.
Purpose To investigate relationships between computer-extracted breast magnetic resonance (MR) imaging phenotypes with multigene assays of MammaPrint, Oncotype DX, and PAM50 to assess the role of ...radiomics in evaluating the risk of breast cancer recurrence. Materials and Methods Analysis was conducted on an institutional review board-approved retrospective data set of 84 deidentified, multi-institutional breast MR examinations from the National Cancer Institute Cancer Imaging Archive, along with clinical, histopathologic, and genomic data from The Cancer Genome Atlas. The data set of biopsy-proven invasive breast cancers included 74 (88%) ductal, eight (10%) lobular, and two (2%) mixed cancers. Of these, 73 (87%) were estrogen receptor positive, 67 (80%) were progesterone receptor positive, and 19 (23%) were human epidermal growth factor receptor 2 positive. For each case, computerized radiomics of the MR images yielded computer-extracted tumor phenotypes of size, shape, margin morphology, enhancement texture, and kinetic assessment. Regression and receiver operating characteristic analysis were conducted to assess the predictive ability of the MR radiomics features relative to the multigene assay classifications. Results Multiple linear regression analyses demonstrated significant associations (R
= 0.25-0.32, r = 0.5-0.56, P < .0001) between radiomics signatures and multigene assay recurrence scores. Important radiomics features included tumor size and enhancement texture, which indicated tumor heterogeneity. Use of radiomics in the task of distinguishing between good and poor prognosis yielded area under the receiver operating characteristic curve values of 0.88 (standard error, 0.05), 0.76 (standard error, 0.06), 0.68 (standard error, 0.08), and 0.55 (standard error, 0.09) for MammaPrint, Oncotype DX, PAM50 risk of relapse based on subtype, and PAM50 risk of relapse based on subtype and proliferation, respectively, with all but the latter showing statistical difference from chance. Conclusion Quantitative breast MR imaging radiomics shows promise for image-based phenotyping in assessing the risk of breast cancer recurrence.
RSNA, 2016 Online supplemental material is available for this article.
Chemotherapy of a combination of DNA alkylating agents, procarbazine and lomustine (CCNU), and a microtubule poison, vincristine, offers a significant benefit to a subset of glioma patients. The ...benefit of this regimen, known as PCV, was recently linked to IDH mutation that occurs frequently in glioma and produces D-2-hydroxyglutarate (D-2-HG), a competitive inhibitor of α-ketoglutarate (α-KG). We report here that D-2-HG inhibits the α-KG-dependent alkB homolog (ALKBH) DNA repair enzymes. Cells expressing mutant IDH display reduced repair kinetics, accumulate more DNA damages, and are sensitized to alkylating agents. The observed sensitization to alkylating agents requires the catalytic activity of mutant IDH to produce D-2-HG and can be reversed by the deletion of mutant IDH allele or overexpression of ALKBH2 or AKLBH3. Our results suggest that impairment of DNA repair may contribute to tumorigenesis driven by IDH mutations and that alkylating agents may merit exploration for treating IDH-mutated cancer patients.
Display omitted
•IDH mutations and 2HG inhibit the DNA repair enzyme ALKBH•IDH mutant cells are sensitized to alkylating agents•Alkylating chemotherapy agents should be explored for treating IDH-mutated tumors
Wang et al. demonstrate that D-2-HG produced by mutant IDH inhibits alkylated DNA repair enzymes, leading to DNA damage and sensitizing IDH mutant cells to alkylating agents. These results suggest that impairment of DNA repair may contribute to tumorigenesis driven by IDH mutations and that alkylating agents should be explored as a therapeutic option for IDH-mutated cancer patients.
We sought to define whether there are intrinsic molecular subtypes of high-grade bladder cancer. Consensus clustering performed on gene expression data from a meta-dataset of high-grade, ...muscle-invasive bladder tumors identified two intrinsic, molecular subsets of high-grade bladder cancer, termed "luminal" and "basal-like," which have characteristics of different stages of urothelial differentiation, reflect the luminal and basal-like molecular subtypes of breast cancer, and have clinically meaningful differences in outcome. A gene set predictor, bladder cancer analysis of subtypes by gene expression (BASE47) was defined by prediction analysis of microarrays (PAM) and accurately classifies the subtypes. Our data demonstrate that there are at least two molecularly and clinically distinct subtypes of high-grade bladder cancer and validate the BASE47 as a subtype predictor. Future studies exploring the predictive value of the BASE47 subtypes for standard of care bladder cancer therapies, as well as novel subtype-specific therapies, will be of interest.
RNA sequencing (RNA-Seq) is often used for transcriptome profiling as well as the identification of novel transcripts and alternative splicing events. Typically, RNA-Seq libraries are prepared from ...total RNA using poly(A) enrichment of the mRNA (mRNA-Seq) to remove ribosomal RNA (rRNA), however, this method fails to capture non-poly(A) transcripts or partially degraded mRNAs. Hence, a mRNA-Seq protocol will not be compatible for use with RNAs coming from Formalin-Fixed and Paraffin-Embedded (FFPE) samples.
To address the desire to perform RNA-Seq on FFPE materials, we evaluated two different library preparation protocols that could be compatible for use with small RNA fragments. We obtained paired Fresh Frozen (FF) and FFPE RNAs from multiple tumors and subjected these to different gene expression profiling methods. We tested 11 human breast tumor samples using: (a) FF RNAs by microarray, mRNA-Seq, Ribo-Zero-Seq and DSN-Seq (Duplex-Specific Nuclease) and (b) FFPE RNAs by Ribo-Zero-Seq and DSN-Seq. We also performed these different RNA-Seq protocols using 10 TCGA tumors as a validation set.The data from paired RNA samples showed high concordance in transcript quantification across all protocols and between FF and FFPE RNAs. In both FF and FFPE, Ribo-Zero-Seq removed rRNA with comparable efficiency as mRNA-Seq, and it provided an equivalent or less biased coverage on gene 3' ends. Compared to mRNA-Seq where 69% of bases were mapped to the transcriptome, DSN-Seq and Ribo-Zero-Seq contained significantly fewer reads mapping to the transcriptome (20-30%); in these RNA-Seq protocols, many if not most reads mapped to intronic regions. Approximately 14 million reads in mRNA-Seq and 45-65 million reads in Ribo-Zero-Seq or DSN-Seq were required to achieve the same gene detection levels as a standard Agilent DNA microarray.
Our results demonstrate that compared to mRNA-Seq and microarrays, Ribo-Zero-Seq provides equivalent rRNA removal efficiency, coverage uniformity, genome-based mapped reads, and consistently high quality quantification of transcripts. Moreover, Ribo-Zero-Seq and DSN-Seq have consistent transcript quantification using FFPE RNAs, suggesting that RNA-Seq can be used with FFPE-derived RNAs for gene expression profiling.
Research in several fields now requires the analysis of data sets in which multiple high-dimensional types of data are available for a common set of objects. In particular, The Cancer Genome Atlas ...(TCGA) includes data from several diverse genomic technologies on the same cancerous tumor samples. In this paper we introduce Joint and Individual Variation Explained (JIVE), a general decomposition of variation for the integrated analysis of such data sets. The decomposition consists of three terms: a low-rank approximation capturing joint variation across data types, low-rank approximations for structured variation individual to each data type, and residual noise. JIVE quantifies the amount of joint variation between data types, reduces the dimensionality of the data and provides new directions for the visual exploration of joint and individual structures. The proposed method represents an extension of Principal Component Analysis and has clear advantages over popular two-block methods such as Canonical Correlation Analysis and Partial Least Squares. A JIVE analysis of gene expression and miRNA data on Glioblastoma Multiforme tumor samples reveals gene—miRNA associations and provides better characterization of tumor types. Data and software are available at https://genome.unc.edu/jive/.
Immune infiltration of the tumor microenvironment has been associated with improved survival for some patients with solid tumors. The precise makeup and prognostic relevance of immune infiltrates ...across a broad spectrum of tumors remain unclear.
Using mRNA sequencing data from The Cancer Genome Atlas (TCGA) from 11 tumor types representing 3485 tumors, we evaluated lymphocyte and macrophage gene expression by tissue type and by genomic subtypes defined within and across tumor tissue of origin (Cox proportional hazards, Pearson correlation). We investigated clonal diversity of B-cell infiltrates through calculating B-cell receptor (BCR) repertoire sequence diversity. All statistical tests were two-sided.
High expression of T-cell and B-cell signatures predicted improved overall survival across many tumor types including breast, lung, and melanoma (breast CD8_T_Cells hazard ratio HR = 0.36, 95% confidence interval CI = 0.16 to 0.81, P = .01; lung adenocarcinoma B_Cell_60gene HR = 0.71, 95% CI = 0.58 to 0.87, P = 7.80E-04; melanoma LCK HR = 0.86, 95% CI = 0.79 to 0.94, P = 6.75E-04). Macrophage signatures predicted worse survival in GBM, as did B-cell signatures in renal tumors (Glioblastoma Multiforme GBM: macrophages HR = 1.62, 95% CI = 1.17 to 2.26, P = .004; renal: B_Cell_60gene HR = 1.17, 95% CI = 1.04 to 1.32, P = .009). BCR diversity was associated with survival beyond gene segment expression in melanoma (HR = 2.67, 95% CI = 1.32 to 5.40, P = .02) and renal cell carcinoma (HR = 0.36, 95% CI = 0.15 to 0.87, P = .006).
These data support existing studies suggesting that in diverse tissue types, heterogeneous immune infiltrates are present and typically portend an improved prognosis. In some tumor types, BCR diversity was also associated with survival. Quantitative genomic signatures of immune cells warrant further testing as prognostic markers and potential biomarkers of response to cancer immunotherapy.
Lymphocytic infiltration of tumors predicts improved survival in patients with breast cancer. Previous studies have suggested that this survival benefit is confined predominantly to the basal-like ...subtype. Immune infiltration in ovarian tumors is also associated with improved prognosis. Currently, it is unclear what aspects of the immune response mediate this improved outcome.
Using The Cancer Genome Atlas mRNA-seq data and a large microarray dataset, we evaluated adaptive immune gene expression by genomic subtype in breast and ovarian cancer. To investigate B-cells observed to be prognostic within specific subtypes, we developed methods to analyze B-cell population diversity and degree of somatic hypermutation (SHM) from B-cell receptor (BCR) sequences in mRNA-seq data.
Improved metastasis-free/progression-free survival was correlated with B-cell gene expression signatures, which were restricted mainly to the basal-like and HER2-enriched breast cancer subtypes and the immunoreactive ovarian cancer subtype. Consistent with a restricted epitope-driven response, a subset of basal-like and HER2-enriched breast tumors and immunoreactive ovarian tumors showed high expression of a low-diversity population of BCR gene segments. More BCR segments showed improved prognosis with increased expression in basal-like breast tumors and immunoreactive ovarian tumors compared with other subtypes. Basal-like and HER2-enriched tumors exhibited more BCR sequence variants in regions consistent with SHM.
Taken together, these data suggest the presence of a productive and potentially restricted antitumor B-cell response in basal-like breast and immunoreactive ovarian cancers. Immunomodulatory therapies that support B-cell responses may be a promising therapeutic approach to targeting these B-cell infiltrated tumors.
Pan-omics, pan-cancer analysis has advanced our understanding of the molecular heterogeneity of cancer. However, such analyses have been limited in their ability to use information from multiple ...sources of data (e.g., omics platforms) and multiple sample sets (e.g., cancer types) to predict clinical outcomes. We address the issue of prediction across multiple high-dimensional sources of data and sample sets by using molecular patterns identified by BIDIFAC+, a method for integrative dimension reduction of bidimensionally-linked matrices, in a Bayesian hierarchical model. Our model performs variable selection through spike-and-slab priors that borrow information across clustered data. We use this model to predict overall patient survival from the Cancer Genome Atlas with data from 29 cancer types and 4 omics sources and use simulations to characterize the performance of the hierarchical spike-and-slab prior. We found that molecular patterns shared across all or most cancers were largely not predictive of survival. However, our model selected patterns unique to subsets of cancers that differentiate clinical tumor subtypes with markedly different survival outcomes. Some of these subtypes were previously established, such as subtypes of uterine corpus endometrial carcinoma, while others may be novel, such as subtypes within a set of kidney carcinomas. Through simulations, we found that the hierarchical spike-and-slab prior performs best in terms of variable selection accuracy and predictive power when borrowing information is advantageous, but also offers competitive performance when it is not. We address the issue of prediction across multiple sources of data by using results from BIDIFAC+ in a Bayesian hierarchical model for overall patient survival. By incorporating spike-and-slab priors that borrow information across cancers, we identified molecular patterns that distinguish clinical tumor subtypes within a single cancer and within a group of cancers. We also corroborate the flexibility and performance of using spike-and-slab priors as a Bayesian variable selection approach.
Epidermal growth factor receptor (EGFR) is a targetable receptor frequently overexpressed in basal-like breast cancer, which comprises most triple-negative breast cancers (TNBCs), the only subtype ...without established targeted therapy.
In this randomized phase II trial, patients with metastatic TNBC received anti-EGFR antibody cetuximab (400 mg/m(2) load then 250 mg/m(2) per week intravenously IV) alone, with carboplatin (area under the curve of 2, once per week IV) added after progression or as concomitant therapy from the beginning. Response rate (RR) was the primary end point; others included time to progression (TTP), overall survival (OS), and toxicity. Embedded correlative studies included molecular subtyping on archival tissue. Fresh tumor tissue before and after 7 to 14 days of therapy was used for microarray analyses exploring EGFR pathway activity and inhibition.
In 102 patients with TNBC, RRs were 6% (two of 31) to cetuximab and 16% (four of 25) to cetuximab plus carboplatin after progression. RR to those treated from the beginning with cetuximab plus carboplatin was 17% (12 of 71); 31% of patients responded or had prolonged disease stabilization. The cetuximab plus carboplatin regimen was well tolerated, but both TTP and OS were short at 2.1 months (95% CI, 1.8 to 5.5 months) and 10.4 months (95% CI, 7.7 to 13.1 months), respectively. Of 73 patients with archival tissue for analysis, 74% had basal-like molecular subtype. Sixteen patients had tumor biopsies before and 1 week after therapy; genomic patterns of the EGFR pathway showed activation in 13 and inhibition by therapy in five.
Despite strong preclinical data, combination cetuximab plus carboplatin in metastatic TNBC produced responses in fewer than 20% of patients. EGFR pathway analysis showed that most TNBCs involved activation. However, cetuximab blocked expression of the EGFR pathway in only a minority, suggesting that most had alternate mechanisms for pathway activation.