Summary We analyzed primary breast cancers by genomic DNA copy number arrays, DNA methylation, exome sequencing, mRNA arrays, microRNA sequencing and reverse phase protein arrays. Our ability to integrate information across platforms provided key insights into previously-defined gene expression subtypes and demonstrated the existence of four main breast cancer classes when combining data from five platforms, each of which shows significant molecular heterogeneity. Somatic mutations in only three genes (TP53, PIK3CA and GATA3) occurred at > 10% incidence across all breast cancers; however, there were numerous subtype-associated and novel gene mutations including the enrichment of specific mutations in GATA3, PIK3CA and MAP3K1 with the Luminal A subtype. We identified two novel protein expression-defined subgroups, possibly contributed by stromal/microenvironmental elements, and integrated analyses identified specific signaling pathways dominant in each molecular subtype including a HER2/p-HER2/HER1/p-HER1 signature within the HER2-Enriched expression subtype. Comparison of Basal-like breast tumors with high-grade Serous Ovarian tumors showed many molecular commonalities, suggesting a related etiology and similar therapeutic opportunities. The biologic finding of the four main breast cancer subtypes caused by different subsets of genetic and epigenetic abnormalities raises the hypothesis that much of the clinically observable plasticity and heterogeneity occurs within, and not across, these major biologic subtypes of breast cancer.

ARID1A (the AT-rich interaction domain 1A, also known as BAF250a) is one of the most commonly mutated genes in cancer1,2. The majority of ARID1A mutations are inactivating mutations and lead to loss of ARID1A expression3, which makes ARID1A a poor therapeutic target. Therefore, it is of clinical importance to identify molecular consequences of ARID1A deficiency that create therapeutic vulnerabilities in ARIDIA-mutant tumors. In a proteomic screen, we found that ARID1A interacts with mismatch repair (MMR) protein MSH2. ARID1A recruited MSH2 to chromatin during DNA replication and promoted MMR. Conversely, ARID1A inactivation compromised MMR and increased mutagenesis. ARID1A deficiency correlated with microsatellite instability genomic signature and a predominant C>T mutation pattern and increased mutation load across multiple human cancer types. Tumors formed by an ARID1A-deficient ovarian cancer cell line in syngeneic mice displayed increased mutation load, elevated numbers of tumor-infiltrating lymphocytes, and PD-L1 expression. Notably, treatment with anti-PD-L1 antibody reduced tumor burden and prolonged survival of mice bearing ARIDIA-deficient but not ARID1A-wild-type ovarian tumors. Together, these results suggest ARID1A deficiency contributes to impaired MMR and mutator phenotype in cancer, and may cooperate with immune checkpoint blockade therapy.

The claudin-low subtype is a recently identified rare molecular subtype of human breast cancer that expresses low levels of tight and adherens junction genes and shows high expression of epithelial-to-mesenchymal transition (EMT) genes. These tumors are enriched in gene expression signatures derived from human tumorinitiating cells (TICs) and human mammary stem cells. Through crossspecies analysis, we discovered mouse mammary tumors that have similar gene expression characteristics as human claudin-low tumors and were also enriched for the human TIC signature. Such claudinlow tumors were similarly rare but came from a number of distinct mouse models, including the p53 null transplant model. Here we present a molecular characterization of 50 p53 null mammary tumors compared with other mouse models and human breast tumor subtypes. Similar to human tumors, the murine p53 null tumors fell into multiple molecular subtypes, including two basallike, a luminal, a claudin-low, and a subtype unique to this model. The claudin-low tumors also showed high gene expression of EMT inducers, low expression of the miR-200 family, and low to absent expression of both claudin 3 and E-cadherin. These murine subtypes also contained distinct genomic DNA copy number changes, some of which are similarly altered in their cognate human subtype counterpart. Finally, limiting dilution transplantation revealed that p53 null claudin-low tumors are highly enriched for TICs compared with the more common adenocarcinomas arising in the same model, thus providing a unique preclinical mouse model to investigate the therapeutic response of TICs. genetically engineered mouse model | gene profiling | array comparative genomic hybridization B reast cancer (BC) is the second leading cause of cancerrelated deaths among women in the United States (1). The large compendium of underlying genetic alterations and the resulting histological and molecular subtypes illustrate the heterogeneous nature of this disease. Both this intertumor heterogeneity and the cellular heterogeneity found within a breast tumor (intratumor heterogeneity) are major obstacles to effective treatments. One common feature of BC (and most cancers) is the loss of the tumor suppressor p53 function. p53 has been shown to be mutated in ≈40% of BCs, associated with poor clinical outcomes, and a higher frequency of mutations occurs in more-aggressive molecular subtypes, including the basal-like subtype of human BC (2).Mice homozygous for p53 loss have been shown to develop lymphomas and sarcomas with a short latency (3, 4). When crossed into the BALB/c background, mammary tumors were observed in p53 +/− mice (5). To circumvent the appearance of other tumor types that occurred with short latency, the model was further modified (6); namely, 6-wk-old p53 −/− glands were removed and transplanted into 3-wk-old wild-type BALB/c recipients. These mice develop mammary tumors stochastically with an average latency of approximately 12 mo. Interestingly, the p53 null epithelium initially forms normal du...

BackgroundRNA 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.ResultsTo 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.ConclusionsOur 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 materialThe online version of this article (doi: 10.1186/1471-2164-15-419) contains supplementary material, which is available to authorized users.

Although the MYC oncogene has been implicated in cancer, a systematic assessment of alterations of MYC, related transcription factors, and co-regulatory proteins, forming the proximal MYC network (PMN), across human cancers is lacking. Using computational approaches, we define genomic and proteomic features associated with MYC and the PMN across the 33 cancers of The Cancer Genome Atlas. Pan-cancer, 28% of all samples had at least one of the MYC paralogs amplified. In contrast, the MYC antagonists MGA and MNT were the most frequently mutated or deleted members, proposing a role as tumor suppressors. MYC alterations were mutually exclusive with PIK3CA, PTEN, APC, or BRAF alterations, suggesting that MYC is a distinct oncogenic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such as immune response and growth factor signaling; chromatin, translation, and DNA replication/repair were conserved pan-cancer. This analysis reveals insights into MYC biology and is a reference for biomarkers and therapeutics for cancers with alterations of MYC or the PMN.

PARP inhibitors (PARPi) have shown remarkable therapeutic efficacy against BRCA1/2-mutant cancers through a synthetic lethal interaction. PARPi exert their therapeutic effects mainly through the blockade of ssDNA damage repair, which leads to the accumulation of toxic DNA double-strand breaks specifically in cancer cells with DNA repair deficiency (BCRAness), including those harboring BRCA1/2 mutations. Here we show that PARPi-mediated modulation of the immune response contributes to their therapeutic effects independently of BRCA1/2 mutations. PARPi promoted accumulation of cytosolic DNA fragments because of unresolved DNA lesions, which in turn activated the DNA-sensing cGAS-STING pathway and stimulated production of type I IFNs to induce antitumor immunity independent of BRCAness. These effects of PARPi were further enhanced by immune checkpoint blockade. Overall, these results provide a mechanistic rationale for using PARPi as immunomodulatory agents to harness the therapeutic efficacy of immune checkpoint blockade. Significance: This work uncovers the mechanism behind the clinical efficacy of PARPi in patients with both BRCAwild-type and BRCA-mutant tumors and provides a rationale for combining PARPi with immunotherapy in patients with cancer.

Poly(ADP-ribose) polymerase inhibitors (PARPi) are selectively active in cells with homologous recombination (HR) deficiency (HRD) caused by mutations in BRCA1, BRCA2, and other pathway members. We sought small molecules that induce HRD in HR-competent cells to induce synthetic lethality with PARPi and extend the utility of PARPi. We demonstrated that inhibition of bromodomain containing 4 (BRD4) induced HRD and sensitized cells across multiple tumor lineages to PARPi regardless of BRCA1/2, TP53, RAS, or BRAF mutation status through depletion of the DNA double-stand break resection protein CtIP (C-terminal binding protein interacting protein). Importantly, BRD4 inhibitor (BRD4i) treatment reversed multiple mechanisms of resistance to PARPi. Furthermore, PARPi and BRD4i are synergistic in multiple in vivo models.

Summary Cancer cell lines are major model systems for mechanistic investigation and drug development. However, protein expression data linked to high-quality DNA, RNA and drug screening data have not been available across a large number of cancer cell lines. Using reverse-phase protein arrays, we measured expression levels of ∼230 key cancer-related proteins in >650 independent cell lines, many of which have publically available genomic, transcriptomic and drug screening data. Our dataset recapitulates the effects of mutated pathways on protein expression observed in patient samples, and demonstrates that proteins and particularly phosphoproteins provide information for predicting drug sensitivity that is not available from the corresponding mRNAs. We also developed a user-friendly bioinformatic resource, MCLP, to help serve the biomedical research community.