David Norman Brown

Principles Governing A-to-I RNA Editing in the Breast Cancer Transcriptome.

Summary Little is known about how RNA editing operates in cancer. Transcriptome analysis of 68 normal and cancerous breast tissues revealed that the editing enzyme ADAR acts uniformly, on the same loci, across tissues. In controlled ADAR expression experiments, the editing frequency increased at all loci with ADAR expression levels according to the logistic model. Loci-specific “editabilities,” i.e., propensities to be edited by ADAR, were quantifiable by fitting the logistic function to dose-response data. The editing frequency was increased in tumor cells in comparison to normal controls. Type I interferon response and ADAR DNA copy number together explained 53% of ADAR expression variance in breast cancers. ADAR silencing using small hairpin RNA lentivirus transduction in breast cancer cell lines led to less cell proliferation and more apoptosis. A-to-I editing is a pervasive, yet reproducible, source of variation that is globally controlled by 1q amplification and inflammation, both of which are highly prevalent among human cancers.

Genomic Characterization of Primary Invasive Lobular Breast Cancer

PURPOSE Invasive lobular breast cancer (ILBC) is the second most common histologic subtype after invasive ductal breast cancer (IDBC). Despite clinical and pathologic differences, ILBC is still treated as IDBC. We aimed to identify genomic alterations in ILBC with potential clinical implications. METHODS From an initial 630 ILBC primary tumors, we interrogated oncogenic substitutions and insertions and deletions of 360 cancer genes and genome-wide copy number aberrations in 413 and 170 ILBC samples, respectively, and correlated those findings with clinicopathologic and outcome features. RESULTS Besides the high mutation frequency of CDH1 in 65% of tumors, alterations in one of the three key genes of the phosphatidylinositol 3-kinase pathway, PIK3CA, PTEN, and AKT1, were present in more than one-half of the cases. HER2 and HER3 were mutated in 5.1% and 3.6% of the tumors, with most of these mutations having a proven role in activating the human epidermal growth factor receptor/ERBB pathway. Mutations in FOXA1 and ESR1 copy number gains were detected in 9% and 25% of the samples. All these alterations were more frequent in ILBC than in IDBC. The histologic diversity of ILBC was associated with specific alterations, such as enrichment for HER2 mutations in the mixed, nonclassic, and ESR1 gains in the solid subtype. Survival analyses revealed that chromosome 1q and 11p gains showed independent prognostic value in ILBC and that HER2 and AKT1 mutations were associated with increased risk of early relapse. CONCLUSION This study demonstrates that we can now begin to individualize the treatment of ILBC, with HER2, HER3, and AKT1 mutations representing high-prevalence therapeutic targets and FOXA1 mutations and ESR1 gains deserving urgent dedicated clinical investigation, especially in the context of endocrine treatment.

Uncovering the genomic heterogeneity of multifocal breast cancer

Multifocal breast cancer (MFBC), defined as multiple synchronous unilateral lesions of invasive breast cancer, is relatively frequent and has been associated with more aggressive features than unifocal cancer. Here, we aimed to investigate the genomic heterogeneity between MFBC lesions sharing similar histopathological parameters. Characterization of different lesions from 36 patients with ductal MFBC involved the identification of non‐silent coding mutations in 360 protein‐coding genes (171 tumour and 36 matched normal samples). We selected only patients with lesions presenting the same grade, ER, and HER2 status. Mutations were classified as ‘oncogenic’ in the case of recurrent substitutions reported in COSMIC or truncating mutations affecting tumour suppressor genes. All mutations identified in a given patient were further interrogated in all samples from that patient through deep resequencing using an orthogonal platform. Whole‐genome rearrangement screen was further conducted in 8/36 patients. Twenty‐four patients (67%) had substitutions/indels shared by all their lesions, of which 11 carried the same mutations in all lesions, and 13 had lesions with both common and private mutations. Three‐quarters of those 24 patients shared oncogenic variants. The remaining 12 patients (33%) did not share any substitution/indels, with inter‐lesion heterogeneity observed for oncogenic mutation(s) in genes such as PIK3CA, TP53, GATA3, and PTEN. Genomically heterogeneous lesions tended to be further apart in the mammary gland than homogeneous lesions. Genome‐wide analyses of a limited number of patients identified a common somatic background in all studied MFBCs, including those with no mutation in common between the lesions. To conclude, as the number of molecular targeted therapies increases and trials driven by genomic screening are ongoing, our findings highlight the presence of genomic inter‐lesion heterogeneity in one‐third, despite similar pathological features. This implies that deeper molecular characterization of all MFBC lesions is warranted for the adequate management of those cancers.

Transfer of clinically relevant gene expression signatures in breast cancer: from Affymetrix microarray to Illumina RNA-Sequencing technology

BackgroundMicroarrays have revolutionized breast cancer (BC) research by enabling studies of gene expression on a transcriptome-wide scale. Recently, RNA-Sequencing (RNA-Seq) has emerged as an alternative for precise readouts of the transcriptome. To date, no study has compared the ability of the two technologies to quantify clinically relevant individual genes and microarray-derived gene expression signatures (GES) in a set of BC samples encompassing the known molecular BC’s subtypes. To accomplish this, the RNA from 57 BCs representing the four main molecular subtypes (triple negative, HER2 positive, luminal A, luminal B), was profiled with Affymetrix HG-U133 Plus 2.0 chips and sequenced using the Illumina HiSeq 2000 platform. The correlations of three clinically relevant BC genes, six molecular subtype classifiers, and a selection of 21 GES were evaluated.Results16,097 genes common to the two platforms were retained for downstream analysis. Gene-wise comparison of microarray and RNA-Seq data revealed that 52% had a Spearman’s correlation coefficient greater than 0.7 with highly correlated genes displaying significantly higher expression levels. We found excellent correlation between microarray and RNA-Seq for the estrogen receptor (ER; rs = 0.973; 95% CI: 0.971-0.975), progesterone receptor (PgR; rs = 0.95; 0.947-0.954), and human epidermal growth factor receptor 2 (HER2; rs = 0.918; 0.912-0.923), while a few discordances between ER and PgR quantified by immunohistochemistry and RNA-Seq/microarray were observed. All the subtype classifiers evaluated agreed well (Cohen’s kappa coefficients >0.8) and all the proliferation-based GES showed excellent Spearman correlations between microarray and RNA-Seq (all rs >0.965). Immune-, stroma- and pathway-based GES showed a lower correlation relative to prognostic signatures (all rs >0.6).ConclusionsTo our knowledge, this is the first study to report a systematic comparison of RNA-Seq to microarray for the evaluation of single genes and GES clinically relevant to BC. According to our results, the vast majority of single gene biomarkers and well-established GES can be reliably evaluated using the RNA-Seq technology.

Squalene epoxidase is a bona fide oncogene by amplification with clinical relevance in breast cancer.

SQLE encodes squalene epoxidase, a key enzyme in cholesterol synthesis. SQLE has sporadically been reported among copy-number driven transcripts in multi-omics cancer projects. Yet, its functional relevance has never been subjected to systematic analyses. Here, we assessed the correlation of SQLE copy number (CN) and gene expression (GE) across multiple cancer types, focusing on the clinico-pathological associations in breast cancer (BC). We then investigated whether any biological effect of SQLE inhibition could be observed in BC cell line models. Breast, ovarian, and colorectal cancers showed the highest CN driven GE among 8,783 cases from 22 cancer types, with BC presenting the strongest one. SQLE overexpression was more prevalent in aggressive BC, and was an independent prognostic factor of unfavorable outcome. Through SQLE pharmacological inhibition and silencing in a panel of BC cell lines portraying the diversity of SQLE CN and GE, we demonstrated that SQLE inhibition resulted in a copy-dosage correlated decrease in cell viability, and in a noticeable increase in replication time, only in lines with detectable SQLE transcript. Altogether, our results pinpoint SQLE as a bona fide metabolic oncogene by amplification, and as a therapeutic target in BC. These findings could have implications in other cancer types.

Phylogenetic analysis of metastatic progression in breast cancer using somatic mutations and copy number aberrations

Several studies using genome-wide molecular techniques have reported various degrees of genetic heterogeneity between primary tumours and their distant metastases. However, it has been difficult to discern patterns of dissemination owing to the limited number of patients and available metastases. Here, we use phylogenetic techniques on data generated using whole-exome sequencing and copy number profiling of primary and multiple-matched metastatic tumours from ten autopsied patients to infer the evolutionary history of breast cancer progression. We observed two modes of disease progression. In some patients, all distant metastases cluster on a branch separate from their primary lesion. Clonal frequency analyses of somatic mutations show that the metastases have a monoclonal origin and descend from a common ‘metastatic precursor’. Alternatively, multiple metastatic lesions are seeded from different clones present within the primary tumour. We further show that a metastasis can be horizontally cross-seeded. These findings provide insights into breast cancer dissemination.

Somatic mutation, copy number and transcriptomic profiles of primary and matched metastatic estrogen receptor-positive breast cancers

BACKGROUND Estrogen receptor-positive (ER+) breast cancers (BCs) constitute the most frequent BC subtype. The molecular landscape of ER+ relapsed disease is not well characterized. In this study, we aimed to describe the genomic evolution between primary (P) and matched metastatic (M) ER+ BCs after failure of adjuvant therapy. MATERIALS AND METHODS A total of 182 ER+ metastatic BC patients with long-term follow-up were identified from a single institution. P tumor tissue was available for all patients, with 88 having matched M material. According to the availability of tumor material, samples were characterized using a 120 mutational hotspot qPCR, a 29 gene copy number aberrations (CNA) and a 400 gene expression panels. ESR1 mutations were assayed by droplet digital PCR. Molecular alterations were correlated with overall survival (OS) using the Cox proportional hazards regression models. RESULTS The median follow-up was 6.4 years (range 0.5-26.6 years). Genomic analysis of P tumors revealed somatic mutations in PIK3CA, KRAS, AKT1, FGFR3, HRAS and BRAF at frequencies of 41%, 6%, 5%, 2%, 1% and 2%, respectively, and CN amplification of CCND1, ZNF703, FGFR1, RSF1 and PAK1 at 23%, 19%, 17%, 12% and 11%, respectively. Mutations and CN amplifications were largely concordant between P and matched M (>84%). ESR1 mutations were found in 10.8% of the M but none of the P. Thirteen genes, among which ESR1, FOXA1, and HIF1A, showed significant differential expression between P and M. In P, the differential expression of 18 genes, among which IDO1, was significantly associated with OS (FDR < 0.1). CONCLUSIONS Despite the large concordance between P and matched M for the evaluated molecular alterations, potential actionable targets such as ESR1 mutations were found only in M. This supports the importance of characterizing the M disease. Other targets we identified, such as HIF1A and IDO1, warrant further investigation in this patient population.

HER2-overexpressing breast cancers amplify FGFR signaling upon acquisition of resistance to dual therapeutic blockade of HER2

Purpose: Dual blockade of HER2 with trastuzumab and lapatinib or pertuzumab has been shown to be superior to single-agent trastuzumab. However, a significant fraction of HER2-overexpressing (HER2+) breast cancers escape from these drug combinations. In this study, we sought to discover the mechanisms of acquired resistance to the combination of lapatinib + trastuzumab. Experimental Design: HER2+ BT474 xenografts were treated with lapatinib + trastuzumab long-term until resistance developed. Potential mechanisms of acquired resistance were evaluated in lapatinib + trastuzumab-resistant (LTR) tumors by targeted capture next-generation sequencing. In vitro experiments were performed to corroborate these findings, and a novel drug combination was tested against LTR xenografts. Gene expression and copy-number analyses were performed to corroborate our findings in clinical samples. Results: LTR tumors exhibited an increase in FGF3/4/19 copy number, together with an increase in FGFR phosphorylation, marked stromal changes in the tumor microenvironment, and reduced tumor uptake of lapatinib. Stimulation of BT474 cells with FGF4 promoted resistance to lapatinib + trastuzumab in vitro. Treatment with FGFR tyrosine kinase inhibitors reversed these changes and overcame resistance to lapatinib + trastuzumab. High expression of FGFR1 correlated with a statistically shorter progression-free survival in patients with HER2+ early breast cancer treated with adjuvant trastuzumab. Finally, FGFR1 and/or FGF3 gene amplification correlated with a lower pathologic complete response in patients with HER2+ early breast cancer treated with neoadjuvant anti-HER2 therapy. Conclusions: Amplification of FGFR signaling promotes resistance to HER2 inhibition, which can be diminished by the combination of HER2 and FGFR inhibitors. Clin Cancer Res; 23(15); 4323–34. ©2017 AACR.

Genomic DNA Copy Number Aberrations, Histological Diagnosis, Oral Subsite and Aneuploidy in OPMDs/OSCCs.

Oral potentially malignant disorders (OPMDs) characterized by the presence of dysplasia and DNA copy number aberrations (CNAs), may reflect chromosomal instability (CIN) and predispose to oral squamous cell carcinoma (OSCC). Early detection of OPMDs with such characteristics may play a crucial role in OSCC prevention. The aim of this study was to explore the relationship between CNAs, histological diagnosis, oral subsite and aneuploidy in OPMDs/OSCCs. Samples from OPMDs and OSCCs were processed by high-resolution DNA flow cytometry (hr DNA-FCM) to determine the relative nuclear DNA content. Additionally, CNAs were obtained for a subset of these samples by genome-wide array comparative genomic hybridization (aCGH) using DNA extracted from either diploid or aneuploid nuclei suspension sorted by FCM. Our study shows that: i) aneuploidy, global genomic imbalance (measured as the total number of CNAs) and specific focal CNAs occur early in the development of oral cancer and become more frequent at later stages; ii) OPMDs limited to tongue (TNG) mucosa display a higher frequency of aneuploidy compared to OPMDs confined to buccal mucosa (BM) as measured by DNA-FCM; iii) TNG OPMDs/OSCCs show peculiar features of CIN compared to BM OPMDs/OSCCs given the preferential association with total broad and specific focal CNA gains. Follow-up studies are warranted to establish whether the presence of DNA aneuploidy and specific focal or broad CNAs may predict cancer development in non-dysplastic OPMDs.

Immune Infiltration in Invasive Lobular Breast Cancer

Abstract Background Invasive lobular breast cancer (ILC) is the second most common histological subtype of breast cancer after invasive ductal cancer (IDC). Here, we aimed at evaluating the prevalence, levels, and composition of tumor-infiltrating lymphocytes (TILs) and their association with clinico-pathological and outcome variables in ILC, and to compare them with IDC. Methods We considered two patient series with TIL data: a multicentric retrospective series (n = 614) and the BIG 02-98 study (n = 149 ILC and 807 IDC). We compared immune subsets identified by immuno-histochemistry in the ILC (n = 159) and IDC (n = 468) patients from the Nottingham series, as well as the CIBERSORT immune profiling of the ILC (n = 98) and IDC (n = 388) METABRIC and The Cancer Genome Atlas patients. All ILC/IDC comparisons were done in estrogen receptor (ER)–positive/human epidermal growth factor receptor 2 (HER2)–negative tumors. All statistical tests were two-sided. Results TIL levels were statistically significantly lower in ILC compared with IDC (fold-change = 0.79, 95% confidence interval = 0.70 to 0.88, P < .001). In ILC, high TIL levels were associated with young age, lymph node involvement, and high proliferative tumors. In the univariate analysis, high TIL levels were associated with worse prognosis in the retrospective and BIG 02-98 lobular series, although they did not reach statistical significance in the latter. The Nottingham series revealed that the levels of intratumoral but not total CD8+ were statistically significantly lower in ILC compared with IDC. Comparison of the CIBERSORT profiles highlighted statistically significant differences in terms of immune composition. Conclusions This study shows differences between the immune infiltrates of ER-positive/HER2-negative ILC and IDC in terms of prevalence, levels, localization, composition, and clinical associations.