A mouse model featuring tissue-specific deletion of p53 and Brca1 gives rise to mammary tumors with genomic and transcriptomic similarities to human basal-like breast cancer.

PURPOSE AND METHODS: In human basal-like breast cancer, mutations and deletions in TP53 and BRCA1 are frequent oncogenic events. Thus, we interbred mice expressing the CRE-recombinase with mice harboring loxP sites at TP53 and BRCA1 (K14-Cre; p53f/f Brca1f/f) to test the hypothesis that tissue-specific deletion of TP53 and BRCA1 would give rise to tumors reflective of human basal-like breast cancer. RESULTS: In support of our hypothesis, these transgenic mice developed tumors that express basal-like cytokeratins and demonstrated intrinsic gene expression features similar to human basal-like tumors. Array comparative genomic hybridization revealed a striking conservation of copy number alterations between the K14-Cre; p53f/f Brca1f/f mouse model and human basal-like breast cancer. Conserved events included MYC amplification, KRAS amplification, and RB1 loss. Microarray analysis demonstrated that these DNA copy number events also led to corresponding changes in signatures of pathway activation including high proliferation due to RB1 loss. K14-Cre; p53f/f Brca1f/f also matched human basal-like breast cancer for a propensity to have immune cell infiltrates. Given the long latency of K14-Cre; p53f/f Brca1f/f tumors (~ 250 days), we created tumor syngeneic transplant lines, as well as in vitro cell lines, which were tested for sensitivity to carboplatin and paclitaxel. These therapies invoked acute regression, extended overall survival, and resulted in gene expression signatures of an anti-tumor immune response. CONCLUSION: These findings demonstrate that this model is a valuable preclinical resource for the study of human basal-like breast cancer.

Amplification of SOX4 promotes PI3K/Akt signaling in human breast cancer.

PURPOSE: The PI3K/Akt signaling axis contributes to the dysregulation of many dominant features in breast cancer including cell proliferation, survival, metabolism, motility, and genomic instability. While multiple studies have demonstrated that basal-like or triple-negative breast tumors have uniformly high PI3K/Akt activity, genomic alterations that mediate dysregulation of this pathway in this subset of highly aggressive breast tumors remain to be determined. METHODS: In this study, we present an integrated genomic analysis based on the use of a PI3K gene expression signature as a framework to analyze orthogonal genomic data from human breast tumors, including RNA expression, DNA copy number alterations, and protein expression. In combination with data from a genome-wide RNA-mediated interference screen in human breast cancer cell lines, we identified essential genetic drivers of PI3K/Akt signaling. RESULTS: Our in silico analyses identified SOX4 amplification as a novel modulator of PI3K/Akt signaling in breast cancers and in vitro studies confirmed its role in regulating Akt phosphorylation. CONCLUSIONS: Taken together, these data establish a role for SOX4-mediated PI3K/Akt signaling in breast cancer and suggest that SOX4 may represent a novel therapeutic target and/or biomarker for current PI3K family therapies.

Activating PIK3CA Mutations Induce an Epidermal Growth Factor Receptor (EGFR)/Extracellular Signal-regulated Kinase (ERK) Paracrine Signaling Axis in Basal-like Breast Cancer.

Mutations in PIK3CA, the gene encoding the p110α catalytic subunit of phosphoinositide 3-kinase (PI3K) have been shown to transform human mammary epithelial cells (MECs). These mutations are present in all breast cancer subtypes, including basal-like breast cancer (BLBC). Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we identified 72 protein expression changes in human basal-like MECs with knock-in E545K or H1047R PIK3CA mutations versus isogenic MECs with wild-type PIK3CA. Several of these were secreted proteins, cell surface receptors or ECM interacting molecules and were required for growth of PIK3CA mutant cells as well as adjacent cells with wild-type PIK3CA. The proteins identified by MS were enriched among human BLBC cell lines and pointed to a PI3K-dependent amphiregulin/EGFR/ERK signaling axis that is activated in BLBC. Proteins induced by PIK3CA mutations correlated with EGFR signaling and reduced relapse-free survival in BLBC. Treatment with EGFR inhibitors reduced growth of PIK3CA mutant BLBC cell lines and murine mammary tumors driven by a PIK3CA mutant transgene, all together suggesting that PIK3CA mutations promote tumor growth in part by inducing protein changes that activate EGFR.

Cross-species DNA copy number analyses identifies multiple 1q21-q23 subtype-specific driver genes for breast cancer.

A large number of DNA copy number alterations (CNAs) exist in human breast cancers, and thus characterizing the most frequent CNAs is key to advancing therapeutics because it is likely that these regions contain breast tumor 'drivers' (i.e., cancer causal genes). This study aims to characterize the genomic landscape of breast cancer CNAs and identify potential subtype-specific drivers using a large set of human breast tumors and genetically engineered mouse (GEM) mammary tumors. Using a novel method called SWITCHplus, we identified subtype-specific DNA CNAs occurring at a 15% or greater frequency, which excluded many well-known breast cancer-related drivers such as amplification of ERBB2, and deletions of TP53 and RB1. A comparison of CNAs between mouse and human breast tumors identified regions with shared subtype-specific CNAs. Additional criteria that included gene expression-to-copy number correlation, a DawnRank network analysis, and RNA interference functional studies highlighted candidate driver genes that fulfilled these multiple criteria. Numerous regions of shared CNAs were observed between human breast tumors and GEM mammary tumor models that shared similar gene expression features. Specifically, we identified chromosome 1q21-23 as a Basal-like subtype-enriched region with multiple potential driver genes including PI4KB, SHC1, and NCSTN. This step-wise computational approach based on a cross-species comparison is applicable to any tumor type for which sufficient human and model system DNA copy number data exist, and in this instance, highlights that a single region of amplification may in fact harbor multiple driver genes.

SynthEx: a synthetic-normal-based DNA sequencing tool for copy number alteration detection and tumor heterogeneity profiling.

Changes in the quantity of genetic material, known as somatic copy number alterations (CNAs), can drive tumorigenesis. Many methods exist for assessing CNAs using microarrays, but considerable technical issues limit current CNA calling based upon DNA sequencing. We present SynthEx, a novel tool for detecting CNAs from whole exome and genome sequencing. SynthEx utilizes a "synthetic-normal" strategy to overcome technical and financial issues. In terms of accuracy and precision, SynthEx is highly comparable to array-based methods and outperforms sequencing-based CNA detection tools. SynthEx robustly identifies CNAs using sequencing data without the additional costs associated with matched normal specimens.

Comprehensive analysis of The Cancer Genome Atlas reveals a unique gene and non-coding RNA signature of fibrolamellar carcinoma.

Fibrolamellar carcinoma (FLC) is a unique liver cancer primarily affecting young adults and characterized by a fusion event between DNAJB1 and PRKACA. By analyzing RNA-sequencing data from The Cancer Genome Atlas (TCGA) for >9,100 tumors across ~30 cancer types, we show that the DNAJB1-PRKACA fusion is specific to FLCs. We demonstrate that FLC tumors (n = 6) exhibit distinct messenger RNA (mRNA) and long intergenic non-coding RNA (lincRNA) profiles compared to hepatocellular carcinoma (n = 263) and cholangiocarcinoma (n = 36), the two most common liver cancers. We also identify a set of mRNAs (n = 16) and lincRNAs (n = 4), including LINC00473, that distinguish FLC from ~25 other liver and non-liver cancer types. We confirm this unique FLC signature by analysis of two independent FLC cohorts (n = 20 and 34). Lastly, we validate the overexpression of one specific gene in the FLC signature, carbonic anhydrase XII (CA12), at the protein level by western blot and immunohistochemistry. Both the mRNA and lincRNA signatures support a major role for protein kinase A (PKA) signaling in shaping the FLC gene expression landscape, and present novel candidate FLC oncogenes that merit further investigation.

Characterization of copy number alterations in a mouse model of fibrosis-associated hepatocellular carcinoma reveals concordance with human disease.

Hepatocellular carcinoma (HCC) is a prevalent human cancer with rising incidence worldwide. Human HCC is frequently associated with chronic liver inflammation and cirrhosis, pathophysiological processes that are a consequence of chronic viral infection, disturbances in metabolism, or exposure to chemical toxicants. To better characterize the pathogenesis of HCC, we used a human disease-relevant mouse model of fibrosis-associated hepatocarcinogenesis. In this model, marked liver tumor response caused by the promutagenic chemical N-nitrosodiethylamine in the presence of liver fibrosis was associated with epigenetic events indicative of genomic instability. Therefore, we hypothesized that DNA copy number alterations (CNAs), a feature of genomic instability and a common characteristic of cancer, are concordant between human HCC and mouse models of fibrosis-associated hepatocarcinogenesis. We evaluated DNA CNAs and changes in gene expression in the mouse liver (normal, tumor, and nontumor fibrotic tissues). Additionally, we compared our findings to DNA CNAs in human HCC cases (tumor and nontumor cirrhotic/fibrotic tissues) using publicly available data from The Cancer Genome Atlas (TCGA). We observed that while fibrotic liver tissue is largely devoid of DNA CNAs, highly frequently occurring DNA CNAs are found in mouse tumors, which is indicative of a profound increase in chromosomal instability in HCC. The cross-species gene-level comparison of CNAs identified shared regions of CNAs between human fibrosis- and cirrhosis-associated liver tumors and mouse fibrosis-associated HCC. Our results suggest that CNAs most commonly arise in neoplastic tissue rather than in fibrotic or cirrhotic liver, and demonstrate the utility of this mouse model in replicating the molecular features of human HCC.

An integrated genomics approach identifies drivers of proliferation in luminal-subtype human breast cancer.

Elucidating the molecular drivers of human breast cancers requires a strategy that is capable of integrating multiple forms of data and an ability to interpret the functional consequences of a given genetic aberration. Here we present an integrated genomic strategy based on the use of gene expression signatures of oncogenic pathway activity (n = 52) as a framework to analyze DNA copy number alterations in combination with data from a genome-wide RNA-mediated interference screen. We identify specific DNA amplifications and essential genes within these amplicons representing key genetic drivers, including known and new regulators of oncogenesis. The genes identified include eight that are essential for cell proliferation (FGD5, METTL6, CPT1A, DTX3, MRPS23, EIF2S2, EIF6 and SLC2A10) and are uniquely amplified in patients with highly proliferative luminal breast tumors, a clinical subset of patients for which few therapeutic options are effective. This general strategy has the potential to identify therapeutic targets within amplicons through an integrated use of genomic data sets.

Loss of heterozygosity at the CYP2D6 locus in breast cancer: implications for germline pharmacogenetic studies.

BACKGROUND: Controversy exists regarding the impact of CYP2D6 genotype on tamoxifen responsiveness. We examined loss of heterozygosity (LOH) at the CYP2D6 locus and determined its impact on genotyping error when tumor tissue is used as a DNA source. METHODS: Genomic tumor data from the adjuvant and metastatic settings (The Cancer Genome Atlas [TCGA] and Foundation Medicine [FM]) were analyzed to characterize the impact of CYP2D6 copy number alterations (CNAs) and LOH on Hardy Weinberg equilibrium (HWE). Additionally, we analyzed CYP2D6 *4 genotype from formalin-fixed paraffin-embedded (FFPE) tumor blocks containing nonmalignant tissue and buccal (germline) samples from patients on the North Central Cancer Treatment Group (NCCTG) 89-30-52 tamoxifen trial. All statistical tests were two-sided. RESULTS: In TCGA samples (n =627), the CYP2D6 LOH rate was similar in estrogen receptor (ER)-positive (41.2%) and ER-negative (35.2%) but lower in HER2-positive tumors (15.1%) (P < .001). In FM ER+ samples (n = 290), similar LOH rates were observed (40.8%). In 190 NCCTG samples, the agreement between CYP2D6 genotypes derived from FFPE tumors and FFPE tumors containing nonmalignant tissue was moderate (weighted Kappa = 0.74; 95% CI = 0.63 to 0.84). Comparing CYP2D6 genotypes derived from buccal cells to FFPE tumor DNA, CYP2D6*4 genotype was discordant in six of 31(19.4%). In contrast, there was no disagreement between CYP2D6 genotypes derived from buccal cells with FFPE tumors containing nonmalignant tissue. CONCLUSIONS: LOH at the CYP2D6 locus is common in breast cancer, resulting in potential misclassification of germline CYP2D6 genotypes. Tumor DNA should not be used to determine germline CYP2D6 genotype without sensitive techniques to detect low frequency alleles and quality control procedures appropriate for somatic DNA.