An Early Neoplasia Index (ENI10), Based on Molecular Identity of CD10 Cells and Associated Stemness Biomarkers, is a Predictor of Patient Outcome in Many Cancers.

An accurate estimate of patient survival at diagnosis is critical to plan efficient therapeutic options. A simple and multiapplication tool is needed to move forward the precision medicine era. Taking advantage of the broad and high CD10 expression in stem and cancers cells, we evaluated the molecular identity of aggressive cancer cells. We used epithelial primary cells and developed a breast cancer stem cell­based progressive model. The superiority of the early-transformed isolated molecular index was evaluated by large-scale analysis in solid cancers. BMP2-driven cell transformation increases CD10 expression which preserves stemness properties. Our model identified a unique set of 159 genes enriched in G2­M cell-cycle phases and spindle assembly complex. Using samples predisposed to transformation, we confirmed the value of an early neoplasia index associated to CD10 (ENI10) to discriminate premalignant status of a human tissue. Using a stratified Cox model, a large-scale analysis (>10,000 samples, The Cancer Genome Atlas Pan-Cancer) validated a strong risk gradient (HRs reaching HR = 5.15; 95% confidence interval: 4.00­6.64) for high ENI10 levels. Through different databases, Cox regression model analyses highlighted an association between ENI10 and poor progression-free intervals for more than 50% of cancer subtypes tested, and the potential of ENI10 to predict drug efficacy. The ENI10 index constitutes a robust tool to detect pretransformed tissues and identify high-risk patients at diagnosis. Owing to its biological link with refractory cancer stem cells, the ENI10 index constitutes a unique way of identifying effective treatments to improve clinical care. SIGNIFICANCE: We identified a molecular signature called ENI10 which, owing to its biological link with stem cell properties, predicts patient outcome and drugs efficiency in breast and several other cancers. ENI10 should allow early and optimized clinical management of a broad number of cancers, regardless of the stage of tumor progression.

Isocitrate dehydrogenase wt and IDHmut adult-type diffuse gliomas display distinct alterations in ribosome biogenesis and 2'O-methylation of ribosomal RNA.

BACKGROUND: High-grade adult-type diffuse gliomas (HGGs) constitute a heterogeneous group of aggressive tumors that are mostly incurable. Recent advances highlighting the contribution of ribosomes to cancer development have offered new clinical perspectives. Here, we uncovered that isocitrate dehydrogenase (IDH)wt and IDHmut HGGs display distinct alterations of ribosome biology, in terms of rRNA epitranscriptomics and ribosome biogenesis, which could constitute novel hallmarks that can be exploited for the management of these pathologies. METHODS: We analyzed (1) the ribosomal RNA 2'O-ribose methylation (rRNA 2'Ome) using RiboMethSeq and in-house developed bioinformatics tools (https://github.com/RibosomeCRCL/ribomethseq-nfandrRMSAnalyzer) on 3 independent cohorts compiling 71 HGGs (IDHwt n = 30, IDHmut n = 41) and 9 non-neoplastic samples, (2) the expression of ribosome biogenesis factors using medium throughput RT-qPCR as a readout of ribosome biogenesis, and (3) the sensitivity of 5 HGG cell lines to RNA Pol I inhibitors (CX5461, BMH-21). RESULTS: Unsupervised analysis demonstrated that HGGs could be distinguished based on their rRNA 2'Ome epitranscriptomic profile, with IDHwt glioblastomas displaying the most significant alterations of rRNA 2'Ome at specific sites. In contrast, IDHmut HGGs are largely characterized by an overexpression of ribosome biogenesis factors compared to non-neoplastic tissues or IDHwt glioblastomas. Finally, IDHmut HGG-derived spheroids display higher cytotoxicity to CX5461 than IDHwt glioblastoma, while all HGG spheroids display a similar cytotoxicity to BMH-21. CONCLUSIONS: In HGGs, IDH mutational status is associated with specific alterations of the ribosome biology and with distinct sensitivities to RNA Pol I inhibitors.

Dissecting the Origin of Heterogeneity in Uterine and Ovarian Carcinosarcomas.

Gynecologic carcinosarcomas (CS) are biphasic neoplasms composed of carcinomatous (C) and sarcomatous (S) malignant components. Because of their rarity and histologic complexity, genetic and functional studies on CS are scarce and the mechanisms of initiation and development remain largely unknown. Whole-genome analysis of the C and S components reveals shared genomic alterations, thus emphasizing the clonal evolution of CS. Reconstructions of the evolutionary history of each tumor further reveal that C and S samples are composed of both ancestral cell populations and component-specific subclones, supporting a common origin followed by distinct evolutionary trajectories. However, while we do not find any recurrent genomic features associated with phenotypic divergence, transcriptomic and methylome analyses identify a common mechanism across the cohort, the epithelial-to-mesenchymal transition (EMT), suggesting a role for nongenetic factors in inflicting changes to cellular fate. Altogether, these data accredit the hypothesis that CS tumors are driven by both clonal evolution and transcriptomic reprogramming, essential for susceptibility to transdifferentiation upon encountering environmental cues, thus linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences. Significance: We have provided a detailed characterization of the genomic landscape of CS and identified EMT as a common mechanism associated with phenotypic divergence, linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences.

Erratum: Ribosomal RNA 2' O-methylation as a novel layer of inter-tumour heterogeneity in breast cancer.

[This corrects the article DOI: 10.1093/narcan/zcaa036.].

Comprehensive characterization of claudin-low breast tumors reflects the impact of the cell-of-origin on cancer evolution.

Claudin-low breast cancers are aggressive tumors defined by the low expression of key components of cellular junctions, associated with mesenchymal and stemness features. Although they are generally considered as the most primitive breast malignancies, their histogenesis remains elusive. Here we show that this molecular subtype of breast cancers exhibits a significant diversity, comprising three main subgroups that emerge from unique evolutionary processes. Genetic, gene methylation and gene expression analyses reveal that two of the subgroups relate, respectively, to luminal breast cancers and basal-like breast cancers through the activation of an EMT process over the course of tumor progression. The third subgroup is closely related to normal human mammary stem cells. This unique subgroup of breast cancers shows a paucity of genomic aberrations and a low frequency of TP53 mutations, supporting the emerging notion that the intrinsic properties of the cell-of-origin constitute a major determinant of the genetic history of tumorigenesis.

IFN-III is selectively produced by cDC1 and predicts good clinical outcome in breast cancer.

Dendritic cells play a key role in the orchestration of antitumor immune responses. The cDC1 (conventional dendritic cell 1) subset has been shown to be essential for antitumor responses and response to immunotherapy, but its precise role in humans is largely unexplored. Using a multidisciplinary approach, we demonstrate that human cDC1 play an important role in the antitumor immune response through their capacity to produce type III interferon (IFN-λ). By analyzing a large cohort of breast primary tumors and public transcriptomic datasets, we observed specific production of IFN-λ1 by cDC1. In addition, both IFN-λ1 and its receptor were associated with favorable patient outcomes. We show that IFN-III promotes a TH1 microenvironment through increased production of IL-12p70, IFN-γ, and cytotoxic lymphocyte-recruiting chemokines. Last, we showed that engagement of TLR3 is a therapeutic strategy to induce IFN-III production by tumor-associated cDC1. These data provide insight into potential IFN- or cDC1-targeting antitumor therapies.

Ribosomal RNA 2'O-methylation as a novel layer of inter-tumour heterogeneity in breast cancer.

Recent epitranscriptomics studies unravelled that ribosomal RNA (rRNA) 2'O-methylation is an additional layer of gene expression regulation highlighting the ribosome as a novel actor of translation control. However, this major finding lies on evidences coming mainly, if not exclusively, from cellular models. Using the innovative next-generation RiboMeth-seq technology, we established the first rRNA 2'O-methylation landscape in 195 primary human breast tumours. We uncovered the existence of compulsory/stable sites, which show limited inter-patient variability in their 2'O-methylation level, which map on functionally important sites of the human ribosome structure and which are surrounded by variable sites found from the second nucleotide layers. Our data demonstrate that some positions within the rRNA molecules can tolerate absence of 2'O-methylation in tumoral and healthy tissues. We also reveal that rRNA 2'O-methylation exhibits intra- and inter-patient variability in breast tumours. Its level is indeed differentially associated with breast cancer subtype and tumour grade. Altogether, our rRNA 2'O-methylation profiling of a large-scale human sample collection provides the first compelling evidence that ribosome variability occurs in humans and suggests that rRNA 2'O-methylation might represent a relevant element of tumour biology useful in clinic. This novel variability at molecular level offers an additional layer to capture the cancer heterogeneity and associates with specific features of tumour biology thus offering a novel targetable molecular signature in cancer.

Repurposing rotavirus vaccines for intratumoral immunotherapy can overcome resistance to immune checkpoint blockade.

Although immune checkpoint-targeted therapies are currently revolutionizing cancer care, only a minority of patients develop durable objective responses to anti-PD-1, PD-L1, and CTLA-4 therapy. Therefore, new therapeutic interventions are needed to increase the immunogenicity of tumors and overcome the resistance to these immunotherapies. Oncolytic properties of common viruses can be exploited for the priming of antitumor immunity, and such oncolytic viruses are currently in active clinical development in combination with immune checkpoint-targeted therapies. However, the routine implementation of these therapies is limited by their manufacturing constraints, the risk of exposure of clinical staff, and the ongoing regulations on genetically modified organisms. We sought to determine whether anti-infectious disease vaccines could be used as a commercially available source of immunostimulatory agents for cancer immunotherapy. We found that rotavirus vaccines have both immunostimulatory and oncolytic properties. In vitro, they can directly kill cancer cells with features of immunogenic cell death. In vivo, intratumoral rotavirus therapy has antitumor effects that are dependent on the immune system. In several immunocompetent murine tumor models, intratumoral rotavirus overcomes resistance to and synergizes with immune checkpoint-targeted therapy. Heat- and UV-inactivated rotavirus lost their oncolytic activity but kept their synergy with immune checkpoint-targeted antibodies through the up-regulation of the double-stranded RNA receptor retinoic acid-induced gene 1 (RIG-I). Rotavirus vaccines are clinical-grade products used in pediatric and adult populations. Therefore, in situ immunization strategies with intratumoral-attenuated rotavirus could be implemented quickly in the clinic.

Mutational Profile of Aggressive, Localised Prostate Cancer from African Caribbean Men Versus European Ancestry Men.

Causes of high mortality of prostate cancer in men of African ancestry living in the French West Indies are still debated, between suspicions of environmental factors and genetic susceptibility. We report an integrated genomic study of 25 tumour tissues from radical prostatectomy of aggressive (defined by International Society of Urological Pathology ≥3) prostate cancer patients (10 African Caribbean and 15 French Caucasian) using single nucleotide polymorphism arrays, whole-genome sequencing, and RNA sequencing. The results show that African Caribbean tumours are characterised by a more frequent deletion at 1q41-43 encompassing the DNA repair gene PARP1, and a higher proportion of intrachromosomal rearrangements including duplications associated with CDK12 truncating mutations. Transcriptome analyses show an overexpression of genes related to androgen receptor activity in African Caribbean tumours, and of PVT1, a long non-coding RNA located at 8q24 that confirms the strong involvement of this region in prostate tumours from men of African ancestry. Patient summary: Mortality of prostate cancer is higher in African Caribbean men than in French Caucasian men. Specificities of the former could be explained by genomic events linked with key genes such as DNA damage pathway genes PARP1, CDK12, and the oncogenic long non-coding RNA gene PVT1 at the 8q24 prostate cancer susceptibility locus.

Immunological and classical subtypes of oral premalignant lesions.

Oral squamous cell carcinoma (OSCC) is a major cause of cancer-associated morbidity and mortality and may develop from oral premalignant lesions (OPL). An improved molecular classification of OPL may help refining prevention strategies. We identified two main OPL gene-expression subtypes, named immunological and classical, in 86 OPL (discovery dataset). A gene expression-based score was then developed to classify OPL samples from three independent datasets, including 17 (GSE30784),13 (GSE10174) and 15 (GSE85195) OPLs, into either one of the two gene-expression subtypes. Using the single sample gene set enrichment analysis, enrichment scores for immune-related pathways were different between the two OPL subtypes. In OPL from the discovery set, loss of heterozygosities (LOH) at 3p14, 17p13, TP53, 9p21 and 8p22 and miRNA gene expression profiles were analyzed. Deconvolution of the immune infiltrate was performed using the Microenvironment Cell Populations-counter tool. A multivariate analysis revealed that decreased miRNA-142-5p expression (P = 0.0484) and lower T-cell, monocytic and myeloid dendritic cells (MDC) immune infiltration (T-cells, P = 0.0196; CD8 T cells, P = 0.0129; MDC, P = 0.0481; and monocytes, P = 0.0212) were associated with oral cancer development in the immunological subtype only. In contrast, LOH at 3p14 (P = 0.0241), 17p13 (P = 0.0348) and TP53 (P = 0.004) were associated with oral cancer development in the classical subtype only. In conclusion, we identified 2 subtypes of OPLs, namely immune and classical, which may benefit from different and specific personalized prevention interventions.

A 13-gene expression-based radioresistance score highlights the heterogeneity in the response to radiation therapy across HPV-negative HNSCC molecular subtypes.

BACKGROUND: Radiotherapy for head and neck squamous cell carcinomas (HNSCC) is associated with a substantial morbidity and inconsistent efficacy. Human papillomavirus (HPV)-positive status is recognized as a marker of increased radiosensitivity. Our goal was to identify molecular markers associated with benefit to radiotherapy in patients with HPV-negative disease. METHODS: Gene expression profiles from public repositories were downloaded for data mining. Training sets included 421 HPV-negative HNSCC tumors from The Cancer Genome Atlas (TCGA) and 32 HNSCC cell lines with available radiosensitivity data (GSE79368). A radioresistance (RadR) score was computed using the single sample Gene Set Enrichment Analysis tool. The validation sets included two panels of cell lines (NCI-60 and GSE21644) and HPV-negative HNSCC tumor datasets, including 44 (GSE6631), 82 (GSE39366), and 179 (GSE65858) patients, respectively. We finally performed an integrated analysis of the RadR score with known recurrent genomic alterations in HNSCC, patterns of protein expression, biological hallmarks, and patterns of drug sensitivity using TCGA and the E-MTAB-3610 dataset (659 pancancer cell lines, 140 drugs). RESULTS: We identified 13 genes differentially expressed between tumor and normal head and neck mucosa that were associated with radioresistance in vitro and in patients. The 13-gene expression-based RadR score was associated with recurrence in patients treated with surgery and adjuvant radiotherapy but not with surgery alone. It was significantly different among different molecular subtypes of HPV-negative HNSCC and was significantly lower in the "atypical" molecular subtype. An integrated analysis of RadR score with genomic alterations, protein expression, biological hallmarks and patterns of drug sensitivity showed a significant association with CCND1 amplification, fibronectin expression, seven hallmarks (including epithelial-to-mesenchymal transition and unfolded protein response), and increased sensitivity to elesclomol, an HSP90 inhibitor. CONCLUSIONS: Our study highlights the clinical relevance of the molecular classification of HNSCC and the RadR score to refine radiation strategies in HPV-negative disease.

SNP calling from RNA-seq data without a reference genome: identification, quantification, differential analysis and impact on the protein sequence.

SNPs (Single Nucleotide Polymorphisms) are genetic markers whose precise identification is a prerequisite for association studies. Methods to identify them are currently well developed for model species, but rely on the availability of a (good) reference genome, and therefore cannot be applied to non-model species. They are also mostly tailored for whole genome (re-)sequencing experiments, whereas in many cases, transcriptome sequencing can be used as a cheaper alternative which already enables to identify SNPs located in transcribed regions. In this paper, we propose a method that identifies, quantifies and annotates SNPs without any reference genome, using RNA-seq data only. Individuals can be pooled prior to sequencing, if not enough material is available from one individual. Using pooled human RNA-seq data, we clarify the precision and recall of our method and discuss them with respect to other methods which use a reference genome or an assembled transcriptome. We then validate experimentally the predictions of our method using RNA-seq data from two non-model species. The method can be used for any species to annotate SNPs and predict their impact on the protein sequence. We further enable to test for the association of the identified SNPs with a phenotype of interest.

A whole-genome sequence and transcriptome perspective on HER2-positive breast cancers.

HER2-positive breast cancer has long proven to be a clinically distinct class of breast cancers for which several targeted therapies are now available. However, resistance to the treatment associated with specific gene expressions or mutations has been observed, revealing the underlying diversity of these cancers. Therefore, understanding the full extent of the HER2-positive disease heterogeneity still remains challenging. Here we carry out an in-depth genomic characterization of 64 HER2-positive breast tumour genomes that exhibit four subgroups, based on the expression data, with distinctive genomic features in terms of somatic mutations, copy-number changes or structural variations. The results suggest that, despite being clinically defined by a specific gene amplification, HER2-positive tumours melt into the whole luminal-basal breast cancer spectrum rather than standing apart. The results also lead to a refined ERBB2 amplicon of 106 kb and show that several cases of amplifications are compatible with a breakage-fusion-bridge mechanism.

Exploration of the core metabolism of symbiotic bacteria.

BACKGROUND: A large number of genome-scale metabolic networks is now available for many organisms, mostly bacteria. Previous works on minimal gene sets, when analysing host-dependent bacteria, found small common sets of metabolic genes. When such analyses are restricted to bacteria with similar lifestyles, larger portions of metabolism are expected to be shared and their composition is worth investigating. Here we report a comparative analysis of the small molecule metabolism of symbiotic bacteria, exploring common and variable portions as well as the contribution of different lifestyle groups to the reduction of a common set of metabolic capabilities. RESULTS: We found no reaction shared by all the bacteria analysed. Disregarding those with the smallest genomes, we still do not find a reaction core, however we did find a core of biochemical capabilities. While obligate intracellular symbionts have no core of reactions within their group, extracellular and cell-associated symbionts do have a small core composed of disconnected fragments. In agreement with previous findings in Escherichia coli, their cores are enriched in biosynthetic processes whereas the variable metabolisms have similar ratios of biosynthetic and degradation reactions. Conversely, the variable metabolism of obligate intracellular symbionts is enriched in anabolism. CONCLUSION: Even when removing the symbionts with the most reduced genomes, there is no core of reactions common to the analysed symbiotic bacteria. The main reason is the very high specialisation of obligate intracellular symbionts, however, host-dependence alone is not an explanation for such absence. The composition of the metabolism of cell-associated and extracellular bacteria shows that while they have similar needs in terms of the building blocks of their cells, they have to adapt to very distinct environments. On the other hand, in obligate intracellular bacteria, catabolism has largely disappeared, whereas synthetic routes appear to have been selected for depending on the nature of the symbiosis. As more genomes are added, we expect, based on our simulations, that the core of cell-associated and extracellular bacteria continues to diminish, converging to approximately 60 reactions.

KISSPLICE: de-novo calling alternative splicing events from RNA-seq data.

BACKGROUND: In this paper, we address the problem of identifying and quantifying polymorphisms in RNA-seq data when no reference genome is available, without assembling the full transcripts. Based on the fundamental idea that each polymorphism corresponds to a recognisable pattern in a De Bruijn graph constructed from the RNA-seq reads, we propose a general model for all polymorphisms in such graphs. We then introduce an exact algorithm, called KISSPLICE, to extract alternative splicing events. RESULTS: We show that KISSPLICE enables to identify more correct events than general purpose transcriptome assemblers. Additionally, on a 71 M reads dataset from human brain and liver tissues, KISSPLICE identified 3497 alternative splicing events, out of which 56% are not present in the annotations, which confirms recent estimates showing that the complexity of alternative splicing has been largely underestimated so far. CONCLUSIONS: We propose new models and algorithms for the detection of polymorphism in RNA-seq data. This opens the way to a new kind of studies on large HTS RNA-seq datasets, where the focus is not the global reconstruction of full-length transcripts, but local assembly of polymorphic regions. KISSPLICE is available for download at http://alcovna.genouest.org/kissplice/.