Evolution of synchronous female bilateral breast cancers and response to treatment.

Synchronous bilateral breast cancer (sBBC) occurs after both breasts have been affected by the same germline genetics and environmental exposures. Little evidence exists regarding immune infiltration and response to treatment in sBBCs. Here we show that the impact of the subtype of breast cancer on levels of tumor infiltrating lymphocytes (TILs, n = 277) and on pathologic complete response (pCR) rates (n = 140) differed according to the concordant or discordant subtype of breast cancer of the contralateral tumor: luminal breast tumors with a discordant contralateral tumor had higher TIL levels and higher pCR rates than those with a concordant contralateral tumor. Tumor sequencing revealed that left and right tumors (n = 20) were independent regarding somatic mutations, copy number alterations and clonal phylogeny, whereas primary tumor and residual disease were closely related both from the somatic mutation and from the transcriptomic point of view. Our study indicates that tumor-intrinsic characteristics may have a role in the association of tumor immunity and pCR and demonstrates that the characteristics of the contralateral tumor are also associated with immune infiltration and response to treatment.

Noncanonical splicing junctions between exons and transposable elements represent a source of immunogenic recurrent neo-antigens in patients with lung cancer.

Although most characterized tumor antigens are encoded by canonical transcripts (such as differentiation or tumor-testis antigens) or mutations (both driver and passenger mutations), recent results have shown that noncanonical transcripts including long noncoding RNAs and transposable elements (TEs) can also encode tumor-specific neo-antigens. Here, we investigate the presentation and immunogenicity of tumor antigens derived from noncanonical mRNA splicing events between coding exons and TEs. Comparing human non-small cell lung cancer (NSCLC) and diverse healthy tissues, we identified a subset of splicing junctions that is both tumor specific and shared across patients. We used HLA-I peptidomics to identify peptides encoded by tumor-specific junctions in primary NSCLC samples and lung tumor cell lines. Recurrent junction-encoded peptides were immunogenic in vitro, and CD8+ T cells specific for junction-encoded epitopes were present in tumors and tumor-draining lymph nodes from patients with NSCLC. We conclude that noncanonical splicing junctions between exons and TEs represent a source of recurrent, immunogenic tumor-specific antigens in patients with NSCLC.

Epigenetically controlled tumor antigens derived from splice junctions between exons and transposable elements.

Oncogenesis often implicates epigenetic alterations, including derepression of transposable elements (TEs) and defects in alternative splicing. Here, we explore the possibility that noncanonical splice junctions between exons and TEs represent a source of tumor-specific antigens. We show that mouse normal tissues and tumor cell lines express wide but distinct ranges of mRNA junctions between exons and TEs, some of which are tumor specific. Immunopeptidome analyses in tumor cell lines identified peptides derived from exon-TE splicing junctions associated to MHC-I molecules. Exon-TE junction-derived peptides were immunogenic in tumor-bearing mice. Both prophylactic and therapeutic vaccinations with junction-derived peptides delayed tumor growth in vivo. Inactivation of the TE-silencing histone 3-lysine 9 methyltransferase Setdb1 caused overexpression of new immunogenic junctions in tumor cells. Our results identify exon-TE splicing junctions as epigenetically controlled, immunogenic, and protective tumor antigens in mice, opening possibilities for tumor targeting and vaccination in patients with cancer.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion.

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

Assessing reliability of intra-tumor heterogeneity estimates from single sample whole exome sequencing data.

Tumors are made of evolving and heterogeneous populations of cells which arise from successive appearance and expansion of subclonal populations, following acquisition of mutations conferring them a selective advantage. Those subclonal populations can be sensitive or resistant to different treatments, and provide information about tumor aetiology and future evolution. Hence, it is important to be able to assess the level of heterogeneity of tumors with high reliability for clinical applications. In the past few years, a large number of methods have been proposed to estimate intra-tumor heterogeneity from whole exome sequencing (WES) data, but the accuracy and robustness of these methods on real data remains elusive. Here we systematically apply and compare 6 computational methods to estimate tumor heterogeneity on 1,697 WES samples from the cancer genome atlas (TCGA) covering 3 cancer types (breast invasive carcinoma, bladder urothelial carcinoma, and head and neck squamous cell carcinoma), and two distinct input mutation sets. We observe significant differences between the estimates produced by different methods, and identify several likely confounding factors in heterogeneity assessment for the different methods. We further show that the prognostic value of tumor heterogeneity for survival prediction is limited in those datasets, and find no evidence that it improves over prognosis based on other clinical variables. In conclusion, heterogeneity inference from WES data on a single sample, and its use in cancer prognosis, should be considered with caution. Other approaches to assess intra-tumoral heterogeneity such as those based on multiple samples may be preferable for clinical applications.

Author Correction: New insight for pharmacogenomics studies from the transcriptional analysis of two large-scale cancer cell line panels.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Lymphovascular invasion after neoadjuvant chemotherapy is strongly associated with poor prognosis in breast carcinoma.

PURPOSE: Few studies evaluated the prognostic value of the presence of lymphovascular invasion (LVI) after neoadjuvant chemotherapy (NAC) for breast cancer (BC). METHODS: The association between LVI and survival was evaluated in a cohort of BC patients treated by NAC between 2002 and 2011. Five post-NAC prognostic scores (ypAJCC, RCB, CPS, CPS + EG and Neo-Bioscore) were evaluated and compared with or without the addition of LVI. RESULTS: Out of 1033 tumors, LVI was present on surgical specimens in 29.2% and absent in 70.8% of the cases. Post-NAC LVI was associated with impaired disease-free survival (DFS) (HR 2.54; 95% CI 1.96-3.31; P < 0.001), and the magnitude of this effect depended on BC subtype (Pinteraction = 0.003), (luminal BC: HR 1.83; P = 0.003; triple negative BC: HR 3.73; P < 0.001; HER2-positive BC: HR 6.21; P < 0.001). Post-NAC LVI was an independent predictor of local relapse, distant metastasis, and overall survival; and increased the accuracy of all five post-NAC prognostic scoring systems. CONCLUSIONS: Post-NAC LVI is a strong independent prognostic factor that: (i) should be systematically reported in pathology reports; (ii) should be used as stratification factor after NAC to propose inclusion in second-line trials or adjuvant treatment; (iii) should be included in post-NAC scoring systems.

New insight for pharmacogenomics studies from the transcriptional analysis of two large-scale cancer cell line panels.

One of the most challenging problems in the development of new anticancer drugs is the very high attrition rate. The so-called "drug repositioning process" propose to find new therapeutic indications to already approved drugs. For this, new analytic methods are required to optimize the information present in large-scale pharmacogenomics datasets. We analyzed data from the Genomics of Drug Sensitivity in Cancer and Cancer Cell Line Encyclopedia studies. We focused on common cell lines (n = 471), considering the molecular information, and the drug sensitivity for common drugs screened (n = 15). We propose a novel classification based on transcriptomic profiles of cell lines, according to a biological network-driven gene selection process. Our robust molecular classification displays greater homogeneity of drug sensitivity than cancer cell line grouped based on tissue of origin. We then identified significant associations between cell line cluster and drug response robustly found between both datasets. We further demonstrate the relevance of our method using two additional external datasets and distinct sensitivity metrics. Some associations were still found robust, despite cell lines and drug responses' variations. This study defines a robust molecular classification of cancer cell lines that could be used to find new therapeutic indications to known compounds.

A Stromal Immune Module Correlated with the Response to Neoadjuvant Chemotherapy, Prognosis and Lymphocyte Infiltration in HER2-Positive Breast Carcinoma Is Inversely Correlated with Hormonal Pathways.

INTRODUCTION: HER2-positive breast cancer (BC) is a heterogeneous group of aggressive breast cancers, the prognosis of which has greatly improved since the introduction of treatments targeting HER2. However, these tumors may display intrinsic or acquired resistance to treatment, and classifiers of HER2-positive tumors are required to improve the prediction of prognosis and to develop novel therapeutic interventions. METHODS: We analyzed 2893 primary human breast cancer samples from 21 publicly available datasets and developed a six-metagene signature on a training set of 448 HER2-positive BC. We then used external public datasets to assess the ability of these metagenes to predict the response to chemotherapy (Ignatiadis dataset), and prognosis (METABRIC dataset). RESULTS: We identified a six-metagene signature (138 genes) containing metagenes enriched in different gene ontologies. The gene clusters were named as follows: Immunity, Tumor suppressors/proliferation, Interferon, Signal transduction, Hormone/survival and Matrix clusters. In all datasets, the Immunity metagene was less strongly expressed in ER-positive than in ER-negative tumors, and was inversely correlated with the Hormonal/survival metagene. Within the signature, multivariate analyses showed that strong expression of the "Immunity" metagene was associated with higher pCR rates after NAC (OR = 3.71[1.28-11.91], p = 0.019) than weak expression, and with a better prognosis in HER2-positive/ER-negative breast cancers (HR = 0.58 [0.36-0.94], p = 0.026). Immunity metagene expression was associated with the presence of tumor-infiltrating lymphocytes (TILs). CONCLUSION: The identification of a predictive and prognostic immune module in HER2-positive BC confirms the need for clinical testing for immune checkpoint modulators and vaccines for this specific subtype. The inverse correlation between Immunity and hormone pathways opens research perspectives and deserves further investigation.

No evidence for TSLP pathway activity in human breast cancer.

Thymic stromal lymphopoietin (TSLP) is an epithelial cell-derived cytokine that primes dendritic cells for Th2 induction. It has been implicated in different types of allergic diseases. Recent work suggested that TSLP could play an important role in the tumor microenvironment and influence tumor progression, in particular in breast cancer. In this study we systematically assessed the production of TSLP at the mRNA and protein levels in several human breast cancer cell lines, large-scale public transcriptomics data sets, and primary human breast tumors. We found that TSLP production was marginal, and concerned less than 10% of the tumors, with very low mRNA and protein levels. In most cases TSLP was undetectable and found to be expressed at lower levels in breast cancer as compared to normal breast tissue. Last, we could not detect any functional TSLP receptor (TSLPR) expression neither on hematopoietic cells nor on stromal cells within the primary tumor microenvironment. We conclude that TSLP-TSLPR pathway activity is not significantly detected within human breast cancer. Taken together, these observations do not support TSLP targeting in breast cancer.

Impaired PRC2 activity promotes transcriptional instability and favors breast tumorigenesis.

Alterations of chromatin modifiers are frequent in cancer, but their functional consequences often remain unclear. Focusing on the Polycomb protein EZH2 that deposits the H3K27me3 (trimethylation of Lys27 of histone H3) mark, we showed that its high expression in solid tumors is a consequence, not a cause, of tumorigenesis. In mouse and human models, EZH2 is dispensable for prostate cancer development and restrains breast tumorigenesis. High EZH2 expression in tumors results from a tight coupling to proliferation to ensure H3K27me3 homeostasis. However, this process malfunctions in breast cancer. Low EZH2 expression relative to proliferation and mutations in Polycomb genes actually indicate poor prognosis and occur in metastases. We show that while altered EZH2 activity consistently modulates a subset of its target genes, it promotes a wider transcriptional instability. Importantly, transcriptional changes that are consequences of EZH2 loss are predominantly irreversible. Our study provides an unexpected understanding of EZH2's contribution to solid tumors with important therapeutic implications.