A correction to Statistical modeling for selecting housekeeper genes by Aniko Szabo, Charles M Perou, Mehmet Karaca, Laurent Perreard, John F Quackenbush, and Philip S Bernard. Genome Biology 2004, 5:R59
Mammary gland morphology and physiology are supported by an underlying cellular differentiation hierarchy. Molecular features associated with particular cell types along this hierarchy may contribute to the biological and clinical heterogeneity observed in human breast carcinomas. Investigating the normal cellular developmental phenotypes in breast tumors may provide new prognostic paradigms, identify new targetable pathways, and explain breast cancer subtype etiology. We used transcriptomic profiles coming from fluorescence-activated cell sorted (FACS) normal mammary epithelial cell types from several independent human and murine studies. Using a meta-analysis approach, we derived consensus gene signatures for both species and used these to relate tumors to normal mammary epithelial cell phenotypes. We then compiled a dataset of breast cancer patients treated with neoadjuvant anthracycline and taxane chemotherapy regimens to determine if normal cellular traits predict the likelihood of a pathological complete response (pCR) in a multivariate logistic regression analysis with clinical markers and genomic features such as cell proliferation. Most human and murine tumor subtypes shared some, but not all, features with a specific FACS-purified normal cell type; thus for most tumors a potential distinct cell type of ‘origin’ could be assigned. We found that both human luminal progenitor and mouse fetal mammary stem cell features predicted pCR sensitivity across all breast cancer patients even after controlling for intrinsic subtype, proliferation, and clinical variables. This work identifies new clinically relevant gene signatures and highlights the value of a developmental biology perspective for uncovering relationships between tumor subtypes and their potential normal cellular counterparts.
Electronic supplementary material
The online version of this article (doi:10.1007/s10549-014-3262-6) contains supplementary material, which is available to authorized users.
Breast cancer; Comparative genomics; Genetically engineered mouse models; Genomic signatures; Neoadjuvant chemotherapy; Normal mammary tissue
Both abundant epidermal growth factor receptor (EGFR or ErbB1) and high activity of the phosphatidyl-inositol 3-kinase (PI3K)–Akt pathway are common and therapeutically targeted in triple-negative breast cancer (TNBC). However, activation of another EGFR family member [human epidermal growth factor receptor 3 (HER3) (or ErbB3)] may limit the antitumor effects of these drugs. We found that TNBC cell lines cultured with the EGFR or HER3 ligand EGF or heregulin, respectively, and treated with either an Akt inhibitor (GDC-0068) or a PI3K inhibitor (GDC-0941) had increased abundance and phosphorylation of HER3. The phosphorylation of HER3 and EGFR in response to these treatments was reduced by the addition of a dual EGFR and HER3 inhibitor (MEHD7945A). MEHD7945A also decreased the phosphorylation (and activation) of EGFR and HER3 and the phosphorylation of downstream targets that occurred in response to the combination of EGFR ligands and PI3K-Akt pathway inhibitors. In culture, inhibition of the PI3K-Akt pathway combined with either MEHD7945A or knockdown of HER3 decreased cell proliferation compared with inhibition of the PI3K-Akt pathway alone. Combining either GDC-0068 or GDC-0941 with MEHD7945A inhibited the growth of xenografts derived from TNBC cell lines or from TNBC patient tumors, and this combination treatment was also more effective than combining either GDC-0068 or GDC-0941 with cetuximab, an EGFR-targeted antibody. After therapy with EGFR-targeted antibodies, some patients had residual tumors with increased HER3 abundance and EGFR/HER3 dimerization (an activating interaction). Thus, we propose that concomitant blockade of EGFR, HER3, and the PI3K-Akt pathway in TNBC should be investigated in the clinical setting.
To use genetically engineered mouse models (GEMMs) and orthotopic syngeneic murine transplants (OSTs) to develop gene-expression based predictors of response to anti-cancer drugs in human tumors. These mouse models offer advantages including precise genetics and an intact microenvironment/immune system.
We examined the efficacy of four chemotherapeutic or targeted anti-cancer drugs, alone and in combination, using mouse models representing three distinct breast cancer subtypes: Basal-like (C3(1)-T-antigen GEMM), Luminal B (MMTV-Neu GEMM), and Claudin-low (T11/TP53−/− OST). We expression-profiled tumors to develop signatures that corresponded to treatment and response, then tested their predictive potential using human patient data.
Although a single agent exhibited exceptional efficacy (i.e. lapatinib in the Neu-driven model), generally single-agent activity was modest, while some combination therapies were more active and life-prolonging. Through analysis of RNA expression in this large set of chemotherapy-treated murine tumors, we identified a pair of gene expression signatures that predicted pathological complete response to neoadjuvant anthracycline/taxane therapy in human patients with breast cancer.
These results show that murine-derived gene signatures can predict response even after accounting for common clinical variables and other predictive genomic signatures, suggesting that mice can be used to identify new biomarkers for human cancer patients.
Brain metastases (BM) arising from Triple-negative breast cancer (TNBC) portend poor prognosis. TNBC is more common in premenopausal and African-American (AA) patients; both also confer poor prognosis. In a single institution cohort study, we sought to determine if inferior outcome of TN BCBM is more reflective of a higher-risk population or subtype itself.
The UNC Breast Cancer Database identified pts with BCBM diagnosed 1988 – 2008. BC subtype was assigned by IHC: HR+ (hormone receptor, ER+ and/or PR+)/HER2−, HR+/HER2+, HR−/HER2+ and TN (ER−/PR−/HER2−). Survival and recurrence patterns were evaluated by subtype, age (< vs ≥ 40 years) and race (AA vs non-AA) using the Kaplan-Meier method and Cox regression.
Among 119 patients with BCBM, 33% were AA and 31% aged < 40 yrs. BC subtype was confirmed in 98 patients: 30% HR+/HER2−, 21% HR+/HER2+, 18% HR−/HER2+, 31% TN. Survival after BM was impacted by subtype (p=0.002), shortest for TNBC (0.24 yrs, CI 0.17 – 0.48). There were no age- (p=0.84) or race-specific (p=0.09) differences in survival after BM; stratification of subtypes by age and race revealed no difference (all, p > 0.1). Receipt of systemic therapy after BCBM was an important predictor of survival following BCBM (HR = 0.29, p=0.002) when adjusted for race, age, number of CNS lesions and BC subtype.
TNBC confers a high risk of death following BM regardless of race and age supporting the need for novel agents capable of controlling both intra- and extracranial TNBC across all races and ages.
Obesity is associated with basal-like breast cancer (BBC), an aggressive breast cancer subtype. The objective of this study was to determine whether obesity promotes BBC onset in adulthood and to evaluate the role of stromal-epithelial interactions in obesity-associated tumorigenesis. We hypothesized that hepatocyte growth factor (HGF) plays a promoting role in BBC, which express the HGF receptor, c-Met. In C3(1)-Tag mice, a murine model of BBC, we demonstrated that obesity leads to a significant increase in HGF secretion and an associated decrease in tumor latency. By immunohistochemical analysis, normal mammary gland exhibited obesity-induced HGF, c-Met and phospho-c-Met, indicating that activation of the cascade was obesity-driven. HGF secretion was also increased from primary mammary fibroblasts isolated from normal mammary glands and tumors of obese mice compared to lean. These results demonstrate that obesity-induced elevation of HGF expression is a stable phenotype, maintained after several passages, and after removal of dietary stimulation. Conditioned media from primary tumor fibroblasts from obese mice drove tumor cell proliferation. In co-culture, neutralization of secreted HGF blunted tumor cell migration, further linking obesity-mediated HGF-dependent effects to in vitro measures of tumor aggressiveness. In sum, these results demonstrate that HGF/c-Met plays an important role in obesity-associated carcinogenesis. Understanding the effects of obesity on risk and progression is important given that epidemiologic studies imply a portion of BBC could be eliminated by reducing obesity.
Basal-like breast cancer; tumor latency; microenvironment; normal mammary gland; high fat diet-induced obesity; fibroblast
We describe the landscape of somatic genomic alterations based on multi-dimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.
Lysophosphatidic acid (LPA) acts through high affinity G protein-coupled receptors to mediate a plethora of physiological and pathological activities associated with tumorigenesis. LPA receptors and autotaxin (ATX/LysoPLD), the primary enzyme producing LPA, are aberrantly expressed in multiple cancer lineages. However, the role of ATX and LPA receptors in the initiation and progression of breast cancer has not been evaluated. We demonstrate that expression of ATX or each Edg-family LPA receptor in mammary epithelium of transgenic mice is sufficient to induce a high frequency of late-onset, estrogen receptor (ER) positive, invasive and metastatic mammary cancer. Thus ATX and LPA receptors can contribute to the initiation and progression of breast cancer.
LPA; ATX; Transgenic mouse model; Breast cancer; Metastasis
Triple-negative breast cancer is a broad and diverse category for which additional subclassifications are needed. Therefore, future clinical trials should stratify patients based on a tumor′s basal-like versus non-basal-like gene expression profile, which appears to be the main biological difference seen within triple-negative breast cancer.
Contrast the definitions of TN and basal-like.Describe the undistinguishable global gene expression patterns of non-basal-like TN tumors versus non-TN tumors that are non-basal-like.Describe the relationship between TN heterogeneity and tumor heterogeneity plus microenvironmental heterogeneity.
Triple-negative (TN) and basal-like (BL) breast cancer definitions have been used interchangeably to identify breast cancers that lack expression of the hormone receptors and overexpression and/or amplification of HER2. However, both classifications show substantial discordance rates when compared to each other. Here, we molecularly characterize TN tumors and BL tumors, comparing and contrasting the results in terms of common patterns and distinct patterns for each. In total, when testing 412 TN and 473 BL tumors, 21.4% and 31.5% were identified as non-BL and non-TN, respectively. TN tumors identified as luminal or HER2-enriched (HER2E) showed undistinguishable overall gene expression profiles when compared versus luminal or HER2E tumors that were not TN. Similar findings were observed within BL tumors regardless of their TN status, which suggests that molecular subtype is preserved regardless of individual marker results. Interestingly, most TN tumors identified as HER2E showed low HER2 expression and lacked HER2 amplification, despite the similar overall gene expression profiles to HER2E tumors that were clinically HER2-positive. Lastly, additional genomic classifications were examined within TN and BL cancers, most of which were highly concordant with tumor intrinsic subtype. These results suggest that future clinical trials focused on TN disease should consider stratifying patients based upon BL versus non-BL gene expression profiles, which appears to be the main biological difference seen in patients with TN breast cancer.
Breast cancer; Subtype; Gene expression; Triple-negative; basal-like
The clinical impact of the biological heterogeneity within HER2-positive (HER2+) breast cancer is not fully understood. Here, we evaluated the molecular features and survival outcomes of the intrinsic subtypes within HER2+ breast cancer.
We interrogated The Cancer Genome Atlas (n = 495) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) datasets (n = 1730) of primary breast cancers for molecular data derived from DNA, RNA and protein, and determined intrinsic subtype. Clinical HER2 status was defined according to American Society of Clinical Oncology (ASCO)/College of American Pathologists (CAP) guidelines or DNA copy-number aberration by single nucleotide polymorphism arrays. Cox models tested the prognostic significance of each variable in patients not treated with trastuzumab (n = 1711).
Compared with clinically HER2 (cHER2)-negative breast cancer, cHER2+ breast cancer had a higher frequency of the HER2-enriched (HER2E) subtype (47.0% vs 7.1%) and a lower frequency of Luminal A (10.7% vs 39.0%) and Basal-like (14.1% vs 23.4%) subtypes. The likelihood of cHER2-positivity in HER2E, Luminal B, Basal-like and Luminal A subtypes was 64.6%, 20.0%, 14.4% and 7.3%, respectively. Within each subtype, only 0.3% to 3.9% of genes were found differentially expressed between cHER2+ and cHER2-negative tumors. Within cHER2+ tumors, HER2 gene and protein expression was statistically significantly higher in the HER2E and Basal-like subtypes than either luminal subtype. Neither cHER2 status nor the new 10-subtype copy number-based classification system (IntClust) added independent prognostic value to intrinsic subtype.
When the intrinsic subtypes are taken into account, cHER2-positivity does not translate into large changes in the expression of downstream signaling pathways, nor does it affect patient survival in the absence of HER2 targeting.
Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype, and fewer than half of patients are cured with standard front-line therapy. To improve therapeutic options, better animal models that accurately mimic human DLBCL (hDLBCL) are needed. Canine DLBCL (cDLBCL), one of the most common cancers in veterinary oncology, is morphologically similar to hDLBCL and is treated using similar chemotherapeutic protocols. With genomic technologies, it is now possible to molecularly evaluate dogs as a potential large-animal model for hDLBCL. We evaluated canine B-cell lymphomas (cBCLs) using immunohistochemistry and gene expression profiling. Canine B-cell lymphoma expression profiles were similar in many ways to hDLBCLs. For instance, a subset had increased expression of NF-κB pathway genes, mirroring human activated B-cell (ABC)-type DLBCL. Furthermore, immunoglobulin heavy chain (IGH) ongoing mutation status, which is correlated with ABC/germinal center B-cell (GCB) cell of origin in hDLBCL, separated cBCL into two groups with statistically different progression-free and overall survival times. In contrast with hDLBCL, cBCL rarely expressed BCL6 and MUM1/IRF4 by immunohistochemistry. Collectively, these studies identify molecular similarities to hDLBCL that introduce pet dogs as a representative model of hDLBCL for future studies, including therapeutic clinical trials.
Current immunohistochemical (IHC)-based definitions of luminal A and B breast cancers are imperfect when compared with multigene expression-based assays. In this study, we sought to improve the IHC subtyping by examining the pathologic and gene expression characteristics of genomically defined luminal A and B subtypes.
Patients and Methods
Gene expression and pathologic features were collected from primary tumors across five independent cohorts: British Columbia Cancer Agency (BCCA) tamoxifen-treated only, Grupo Español de Investigación en Cáncer de Mama 9906 trial, BCCA no systemic treatment cohort, PAM50 microarray training data set, and a combined publicly available microarray data set. Optimal cutoffs of percentage of progesterone receptor (PR) –positive tumor cells to predict survival were derived and independently tested. Multivariable Cox models were used to test the prognostic significance.
Clinicopathologic comparisons among luminal A and B subtypes consistently identified higher rates of PR positivity, human epidermal growth factor receptor 2 (HER2) negativity, and histologic grade 1 in luminal A tumors. Quantitative PR gene and protein expression were also found to be significantly higher in luminal A tumors. An empiric cutoff of more than 20% of PR-positive tumor cells was statistically chosen and proved significant for predicting survival differences within IHC-defined luminal A tumors independently of endocrine therapy administration. Finally, no additional prognostic value within hormonal receptor (HR) –positive/HER2-negative disease was observed with the use of the IHC4 score when intrinsic IHC-based subtypes were used that included the more than 20% PR-positive tumor cells and vice versa.
Semiquantitative IHC expression of PR adds prognostic value within the current IHC-based luminal A definition by improving the identification of good outcome breast cancers. The new proposed IHC-based definition of luminal A tumors is HR positive/HER2 negative/Ki-67 less than 14%, and PR more than 20%.
Breast cancer research is hampered by difficulties in obtaining and studying primary human breast tissue, and by the lack of in vivo preclinical models that reflect patient tumor biology accurately. To overcome these limitations, we propagated a cohort of human breast tumors grown in the epithelium-free mammary fat pad of SCID/Beige and NOD/SCID/IL2γ-receptor null (NSG) mice, under a series of transplant conditions. Both models yielded stably transplantable xenografts at comparably high rates (~21% and ~19%, respectively). Of the conditions tested, xenograft take rate was highest in the presence of a low-dose estradiol pellet. Overall, 32 stably transplantable xenograft lines were established, representing 25 unique patients. Most tumors yielding xenografts were “triple-negative” (ER-PR-HER2+) (n=19). However, we established lines from three ER-PR-HER2+ tumors, one ER+PR-HER2−, one ER+PR+HER2− and one “triple-positive” (ER+PR+HER2+) tumor. Serially passaged xenografts show biological consistency with the tumor of origin, are phenotypically stable across multiple transplant generations at the histologic, transcriptomic, proteomic, and genomic levels, and show comparable treatment responses as those observed clinically. Xenografts representing 12 patients, including two ER+ lines, showed metastasis to the mouse lung. These models thus serve as a renewable, quality-controlled tissue resource for preclinical studies investigating treatment response and metastasis.
BRCA1-mutated breast cancer is associated with basal-like disease; however, it is currently unclear if the presence of a BRCA1 mutation depicts a different entity within this subgroup. In this study, we compared the molecular features among basal-like tumors with and without BRCA1 mutations. Fourteen patients with BRCA1-mutated (nine germline and five somatic) tumors and basal-like disease, and 79 patients with BRCA1 non-mutated tumors and basal-like disease, were identified from the cancer genome atlas dataset. The following molecular data types were evaluated: global gene expression, selected protein and phospho-protein expression, global miRNA expression, global DNA methylation, total number of somatic mutations, TP53 and PIK3CA somatic mutations, and global DNA copy-number aberrations. For intrinsic subtype identification, we used the PAM50 subtype predictor. Within the basal-like disease, we observed minor molecular differences in terms of gene, protein, and miRNA expression, and DNA methylation variation, according to BRCA1 status (either germinal or somatic). However, there were significant differences according to average number of mutations and DNA copy-number aberrations, and four amplified regions (2q32.2, 3q29, 6p22.3, and 22q12.2), which are characteristic in high-grade serous ovarian carcinomas, were observed in both germline and somatic BRCA1-mutated breast tumors. These results suggest that minor, but potentially relevant, baseline molecular features exist among basal-like tumors according to BRCA1 status. Additional studies are needed to better clarify if BRCA1 genetic status is an independent prognostic feature, and more importantly, if BRCA1 mutation status is a predictive biomarker of benefit from DNA-damaging agents among basal-like disease.
Electronic supplementary material
The online version of this article (doi:10.1007/s10549-014-3056-x) contains supplementary material, which is available to authorized users.
Basal-like; BRCA1; Intrinsic subtype; Breast cancer
High-throughput sequencing technologies, including RNA-seq, have made it possible to move beyond gene expression analysis to study transcriptional events including alternative splicing and gene fusions. Furthermore, recent studies in cancer have suggested the importance of identifying transcriptionally altered loci as biomarkers for improved prognosis and therapy. While many statistical methods have been proposed for identifying novel transcriptional events with RNA-seq, nearly all rely on contrasting known classes of samples, such as tumor and normal. Few tools exist for the unsupervised discovery of such events without class labels. In this paper, we present SigFuge for identifying genomic loci exhibiting differential transcription patterns across many RNA-seq samples. SigFuge combines clustering with hypothesis testing to identify genes exhibiting alternative splicing, or differences in isoform expression. We apply SigFuge to RNA-seq cohorts of 177 lung and 279 head and neck squamous cell carcinoma samples from the Cancer Genome Atlas, and identify several cases of differential isoform usage including CDKN2A, a tumor suppressor gene known to be inactivated in a majority of lung squamous cell tumors. By not restricting attention to known sample stratifications, SigFuge offers a novel approach to unsupervised screening of genetic loci across RNA-seq cohorts. SigFuge is available as an R package through Bioconductor.
Nearly half of patients with advanced triple negative breast cancer (TNBC) develop brain metastases (BM) and most will also have uncontrolled extracranial disease. This study evaluated the safety and efficacy of iniparib, a small molecule anti-cancer agent that alters reactive oxygen species tumor metabolism and penetrates the blood brain barrier, with the topoisomerase I inhibitor irinotecan in patients with TNBC-BM. Eligible patients had TNBC with new or progressive BM and received irinotecan and iniparib every 3 weeks. Time to progression (TTP) was the primary end point; secondary endpoints were response rate (RR), clinical benefit rate (CBR), overall survival (OS), toxicity, and health-related quality of life. Correlative endpoints included molecular subtyping and gene expression studies on pre-treatment archival tissues, and determination of germline BRCA1/2 status. Thirty-seven patients began treatment; 34 were evaluable for efficacy. Five of 24 patients were known to carry a BRCA germline mutation (4 BRCA1, 1 BRCA2). Median TTP was 2.14 months and median OS was 7.8 months. Intracranial RR was 12 %, while intracranial CBR was 27 %. Treatment was well-tolerated; the most common grade 3/4 adverse events were neutropenia and fatigue. Grade 3/4 diarrhea was rare (3 %). Intrinsic subtyping revealed 19 of 21 tumors (79 %) were basal-like, and intracranial response was associated with high expression of proliferation-related genes. This study suggests a modest benefit of irinotecan plus iniparib in progressive TNBC-BM. More importantly, this trial design is feasible and lays the foundation for additional studies for this treatment-refractory disease.
Electronic supplementary material
The online version of this article (doi:10.1007/s10549-014-3039-y) contains supplementary material, which is available to authorized users.
Irinotecan; Iniparib; Brain metastases; Breast cancer; Phase II; Triple negative
Mutations in PIK3CA, the gene encoding the p110α catalytic subunit of phosphatidylinositol-3 kinase (PI3K), have been shown to transform mammary epithelial cells (MECs). Studies suggest this transforming activity requires binding of mutant p110α via p85 to phosphorylated YXXM motifs in activated receptor tyrosine kinases (RTKs) or adaptors. Using transgenic mice, we examined if ErbB3, a potent activator of PI3K, is required for mutant PIK3CA-mediated transformation of MECs. Conditional loss of ErbB3 in mammary epithelium resulted in a delay of PIK3CAH1047R-dependent mammary gland hyperplasia, but tumor latency, gene expression and PI3K signaling were unaffected. In ErbB3-deficient tumors, mutant PI3K remained associated with several tyrosyl phosphoproteins, potentially explaining the dispensability of ErbB3 for tumorigenicity and PI3K activity. Similarly, inhibition of ErbB RTKs with lapatinib did not affect PI3K signaling in PIK3CAH1047R-expressing tumors. However, the p110α-specific inhibitor BYL719, in combination with lapatinib impaired mammary tumor growth and PI3K signaling more potently than BYL719 alone. Further, co-inhibition of p110α and ErbB3 potently suppressed proliferation and PI3K signaling in human breast cancer cells harboring PIK3CAH1047R. These data suggest that PIK3CAH1047R-driven tumor growth and PI3K signaling can occur independently of ErbB RTKs. However, simultaneous blockade of p110α and ErbB RTKs results in superior inhibition of PI3K and mammary tumor growth, suggesting a rational therapeutic combination against breast cancers harboring PIK3CA activating mutations.
ErbB3; PIK3CA; PI3K; Breast Cancer; transgenic mice
Identifying somatic mutations is critical for cancer genome characterization and for prioritizing patient treatment. DNA whole exome sequencing (DNA-WES) is currently the most popular technology; however, this yields low sensitivity in low purity tumors. RNA sequencing (RNA-seq) covers the expressed exome with depth proportional to expression. We hypothesized that integrating DNA-WES and RNA-seq would enable superior mutation detection versus DNA-WES alone. We developed a first-of-its-kind method, called UNCeqR, that detects somatic mutations by integrating patient-matched RNA-seq and DNA-WES. In simulation, the integrated DNA and RNA model outperformed the DNA-WES only model. Validation by patient-matched whole genome sequencing demonstrated superior performance of the integrated model over DNA-WES only models, including a published method and published mutation profiles. Genome-wide mutational analysis of breast and lung cancer cohorts (n = 871) revealed remarkable tumor genomics properties. Low purity tumors experienced the largest gains in mutation detection by integrating RNA-seq and DNA-WES. RNA provided greater mutation signal than DNA in expressed mutations. Compared to earlier studies on this cohort, UNCeqR increased mutation rates of driver and therapeutically targeted genes (e.g. PIK3CA, ERBB2 and FGFR2). In summary, integrating RNA-seq with DNA-WES increases mutation detection performance, especially for low purity tumors.
Motivation: Variant detection from next-generation sequencing (NGS) data is an increasingly vital aspect of disease diagnosis, treatment and research. Commonly used NGS-variant analysis tools generally rely on accurately mapped short reads to identify somatic variants and germ-line genotypes. Existing NGS read mappers have difficulty accurately mapping short reads containing complex variation (i.e. more than a single base change), thus making identification of such variants difficult or impossible. Insertions and deletions (indels) in particular have been an area of great difficulty. Indels are frequent and can have substantial impact on function, which makes their detection all the more imperative.
Results: We present ABRA, an assembly-based realigner, which uses an efficient and flexible localized de novo assembly followed by global realignment to more accurately remap reads. This results in enhanced performance for indel detection as well as improved accuracy in variant allele frequency estimation.
Availability and implementation: ABRA is implemented in a combination of Java and C/C++ and is freely available for download at https://github.com/mozack/abra.
Supplementary data are available at Bioinformatics online.
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.
Electronic supplementary material
The online version of this article (doi: 10.1186/1471-2164-15-419) contains supplementary material, which is available to authorized users.
RNA sequencing; FFPE; RNA depletion; Ribo-zero; Gene expression; Microarray
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.
Patients and Methods
In this randomized phase II trial, patients with metastatic TNBC received anti-EGFR antibody cetuximab (400 mg/m2 load then 250 mg/m2 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.
Anticancer drug development is inefficient, but genetically engineered murine models (GEMM) and orthotopic, syngeneic transplants (OST) of cancer may offer advantages to in vitro and xenograft systems.
We assessed the activity of 16 treatment regimens in a RAS-driven, Ink4a/Arf-deficient melanoma GEMM. In addition, we tested a subset of treatment regimens in three breast cancer models representing distinct breast cancer subtypes: claudin-low (T11 OST), basal-like (C3-TAg GEMM), and luminal B (MMTV-Neu GEMM).
Like human RAS-mutant melanoma, the melanoma GEMM was refractory to chemotherapy and single-agent small molecule therapies. Combined treatment with AZD6244 [mitogen-activated protein–extracellular signal-regulated kinase kinase (MEK) inhibitor] and BEZ235 [dual phosphoinositide-3 kinase (PI3K)/mammalian target of rapamycin (mTOR) inhibitor] was the only treatment regimen to exhibit significant antitumor activity, showed by marked tumor regression and improved survival. Given the surprising activity of the "AZD/BEZ" combination in the melanoma GEMM, we next tested this regimen in the "claudin-low" breast cancer model that shares gene expression features with melanoma. The AZD/BEZ regimen also exhibited significant activity in this model, leading us to testing in even more diverse GEMMs of basal-like and luminal breast cancer. The AZD/BEZ combination was highly active in these distinct breast cancer models, showing equal or greater efficacy compared with any other regimen tested in studies of over 700 tumor-bearing mice. This regimen even exhibited activity in lapatinib-resistant HER2+ tumors.
These results show the use of credentialed murine models for large-scale efficacy testing of diverse anticancer regimens and predict that combinations of PI3K/mTOR and MEK inhibitors will show antitumor activity in a wide range of human malignancies.
Perturbations in stem cell activity and differentiation can lead to developmental defects and cancer. We use an approach involving a quantitative model of cell-state transitions in vitro to gain insights into how SLUG/SNAI2, a key developmental transcription factor, modulates mammary epithelial stem cell activity and differentiation in vivo. In the absence of SLUG, stem cells fail to transition into basal progenitor cells, while existing basal progenitor cells undergo luminal differentiation; together, these changes result in abnormal mammary architecture and defects in tissue function. Furthermore, we show that in the absence of SLUG, mammary stem cell activity necessary for tissue regeneration and cancer initiation is lost. Mechanistically, SLUG regulates differentiation and cellular plasticity by recruiting the chromatin modifier lysine-specific demethylase 1 (LSD1) to promoters of lineage-specific genes to repress transcription. Together, these results demonstrate that SLUG plays a dual role in repressing luminal epithelial differentiation while unlocking stem cell transitions necessary for tumorigenesis.
•SLUG promotes stem cell-state transitions necessary for tumorigenesis•Loss of SLUG in mice leads to luminal differentiation of basal/myoepithelial cells•In vitro quantitative modeling can accurately predict epithelial plasticity in vivo•A SLUG/LSD1 complex epigenetically regulates mammary epithelial cell differentiation
Perturbations in stem cell activity and differentiation can lead to developmental defects and cancer. Here, Kuperwasser, Gupta, and colleagues show that the transcription factor SLUG plays a dual role in repressing luminal epithelial differentiation while unlocking stem cell transitions necessary for tissue regeneration and tumorigenesis.
Ionizing radiation is a well established carcinogen in rodent models and a risk factor associated with human cancer. We developed a mouse model that captures radiation effects on host biology by transplanting unirradiated Trp53 null mammary tissue to sham or irradiated hosts. Gene expression profiles of tumors that arose in irradiated mice are distinct from those that arose in naïve hosts. We asked whether expression metaprofiles could discern radiation-preceded human cancer or be informative in sporadic breast cancers.
Affymetrix microarray gene expression data from 56 Trp53 null mammary tumors were used to define gene profiles and a centroid that discriminate tumors arising in irradiated hosts. These were applied to publicly available human cancer data sets.
Host irradiation induces a metaprofile consisting of gene modules representing stem cells, cell motility, macrophages and autophagy. Human orthologs of the host irradiation metaprofile discriminated between radiation-preceded and sporadic human thyroid cancers. An irradiated host centroid was strongly associated with estrogen receptor negative breast cancer. When applied to sporadic human breast cancers, the irradiated host metaprofile strongly associated with basal-like and claudin-low breast cancer intrinsic subtypes. Comparing host irradiation in the context of TGFβ levels showed that inflammation was robustly associated with claudin-low tumors.
Detection of radiation-preceded human cancer by the irradiated host metaprofile raises possibilities of assessing human cancer etiology. Moreover, the association of the irradiated host metaprofiles with estrogen receptor negative status and claudin-low subtype suggests that host processes similar to those induced by radiation underlie sporadic cancers.
ionizing radiation; breast cancer; intrinsic subtypes; etiology
Cancer cells have altered metabolism, with increased glucose uptake, glycolysis, and biomass production. This study performed genomic and metabolomic analyses to elucidate how tumor and stromal genomic characteristics influence tumor metabolism.
Thirty-three breast tumors and six normal breast tissues were analyzed by gene expression microarray and by mass spectrometry for metabolites. Gene expression data and clinical characteristics were evaluated in association with metabolic phenotype. To evaluate the role of stromal interactions in altered metabolism, cocultures were performed using breast cancer cells and primary cancer-associated fibroblasts (CAFs).
Across all metabolites, unsupervised clustering resulted in two main sample clusters. Normal breast tissue and a subset of tumors with less aggressive clinical characteristics had lower levels of nucleic and amino acids and glycolysis byproducts, while more aggressive tumors had higher levels of these Warburg-associated metabolites. While tumor intrinsic subtype did not predict metabolic phenotype, metabolic cluster was significantly associated with expression of a wound response signature. In cocultures, CAFs from basal-like breast cancers increased glucose up take and basal-like epithelial cells increased glucose oxidation and glycogen synthesis, suggesting interplay of stromal and epithelial phenotypes on metabolism. Cytokine arrays identified hepatocyte growth factor (HGF) as a potential mediator of stromal-epithelial interaction, and antibody neutralization of HGF resulted in reduced expression of glucose transporter 1 (GLUT1) and decreased glucose uptake by epithelium.
Both tumor/epithelial and stromal characteristics play important roles in metabolism. Warburg-like metabolism is influenced by changes in stromal-epithelial interactions, including altered expression of HGF/Met pathway and GLUT1 expression.
tumor-stromal cell interactions; metabolomics; breast cancer; genomics; microenvironment