In silico validation of RNA-Seq results can identify gene fusions with oncogenic potential in glioblastoma.

RNA-Sequencing (RNA-Seq) can identify gene fusions in tumors, but not all these fusions have functional consequences. Using multiple data bases, we have performed an in silico analysis of fusions detected by RNA-Seq in tumor samples from 139 newly diagnosed glioblastoma patients to identify in-frame fusions with predictable oncogenic potential. Among 61 samples with fusions, there were 103 different fusions, involving 167 different genes, including 20 known oncogenes or tumor suppressor genes (TSGs), 16 associated with cancer but not oncogenes or TSGs, and 32 not associated with cancer but previously shown to be involved in fusions in gliomas. After selecting in-frame fusions able to produce a protein product and running Oncofuse, we identified 30 fusions with predictable oncogenic potential and classified them into four non-overlapping categories: six previously described in cancer; six involving an oncogene or TSG; four predicted by Oncofuse to have oncogenic potential; and 14 other in-frame fusions. Only 24 patients harbored one or more of these 30 fusions, and only two fusions were present in more than one patient: FGFR3::TACC3 and EGFR::SEPTIN14. This in silico study provides a good starting point for the identification of gene fusions with functional consequences in the pathogenesis or treatment of glioblastoma.

RNA sequencing and Immunohistochemistry Reveal ZFN7 as a Stronger Marker of Survival than Molecular Subtypes in G-CIMP-negative Glioblastoma.

PURPOSE: Glioblastoma is the most aggressive brain tumor in adults and has few therapeutic options. The study of molecular subtype classifications may lead to improved prognostic classification and identification of new therapeutic targets. The Cancer Genome Atlas (TCGA) subtype classification has mainly been applied in U.S. clinical trials, while the intrinsic glioma subtype (IGS) has mainly been applied in European trials. EXPERIMENTAL DESIGN: From paraffin-embedded tumor samples of 432 patients with uniformly treated, newly diagnosed glioblastoma, we built tissue microarrays for IHC analysis and applied RNA sequencing to the best samples to classify them according to TCGA and IGS subtypes. RESULTS: We obtained transcriptomic results from 124 patients. There was a lack of agreement among the three TCGA classificatory algorithms employed, which was not solely attributable to intratumoral heterogeneity. There was overlapping of TCGA mesenchymal subtype with IGS cluster 23 and of TCGA classical subtype with IGS cluster 18. Molecular subtypes were not associated with prognosis, but levels of expression of 13 novel genes were identified as independent prognostic markers in glioma-CpG island methylator phenotype-negative patients, independently of clinical factors and MGMT methylation. These findings were validated in at least one external database. Three of the 13 genes were selected for IHC validation. In particular, high ZNF7 RNA expression and low ZNF7 protein expression were strongly associated with longer survival, independently of molecular subtypes. CONCLUSIONS: TCGA and IGS molecular classifications of glioblastoma have no higher prognostic value than individual genes and should be refined before being applied to clinical trials.

Differential expression of long non-coding RNAs are related to proliferation and histological diversity in follicular lymphomas.

Long non-coding RNAs (lncRNAs) comprise a family of non-coding transcripts that are emerging as relevant gene expression regulators of different processes, including tumour development. To determine the possible contribution of lncRNA to the pathogenesis of follicular lymphoma (FL) we performed RNA-sequencing at high depth sequencing in primary FL samples ranging from grade 1-3A to aggressive grade 3B variants using unpurified (n = 16) and purified (n = 12) tumour cell suspensions from nodal samples. FL grade 3B had a significantly higher number of differentially expressed lncRNAs (dif-lncRNAs) with potential target coding genes related to cell cycle regulation. Nine out of the 18 selected dif-lncRNAs were validated by quantitative real time polymerase chain reaction in an independent series (n = 43) of FL. RP4-694A7.2 was identified as the top deregulated lncRNA potentially involved in cell proliferation. RP4-694A7.2 silencing in the WSU-FSCCL FL cell line reduced cell proliferation due to a block in the G1/S phase. The relationship between RP4-694A7.2 and proliferation was confirmed in primary samples as its expression levels positively related to the Ki-67 proliferation index. In summary, lncRNAs are differentially expressed across the clinico-biological spectrum of FL and a subset of them, related to cell cycle, may participate in cell proliferation regulation in these tumours.

Genetic Abnormalities in Large to Giant Congenital Nevi: Beyond NRAS Mutations.

Large and giant congenital melanocytic nevi (CMN) are rare melanocytic lesions mostly caused by postzygotic NRAS alteration. Molecular characterization is usually focused on NRAS and BRAF genes in a unique biopsy sample of the CMN. However, large/giant CMN may exhibit phenotypic differences among distinct areas, and patients differ in features such as presence of multiple CMN or spilus-like lesions. Herein, we have characterized a series of 21 large/giant CMN including patients with spilus-type nevi (9/21 patients, 42.8%). Overall, 53 fresh frozen biopsy samples corresponding to 40 phenotypically characterized areas of large/giant CMNs and 13 satellite lesions were analyzed with a multigene panel and RNA sequencing. Mutational screening showed mutations in 76.2% (16/21) of large/giant CMNs. A NRAS mutation was found in 57.1% (12/21) of patients, and mutations in other genes such as BRAF, KRAS, APC, and MET were detected in 14.3% (3/21) of patients. RNA sequencing showed the fusion transcript ZEB2-ALK and SOX5-RAF1 in large/giant CMN from two patients without missense mutations. Both alterations were not detected in unaffected skin and were detected in different areas of affected skin. These findings suggest that large/giant CMN may result from distinct molecular events in addition to NRAS mutations, including point mutations and fusion transcripts.

Deep analysis of acquired resistance to FGFR1 inhibitor identifies MET and AKT activation and an expansion of AKT1 mutant cells.

The development of acquired resistance (AR) to tyrosine kinase inhibitors (TKIs) of FGFR1 activation is currently not well understood. To gain a deeper insight into this matter in lung cancer, we used the FGFR1-amplified DMS114 cell line and generated multiple clones with AR to an FGFR1-TKI. We molecularly scrutinized the resistant cells, using whole-exome sequencing, RNA sequencing and global DNA methylation analysis. Our results show a de novo activation of AKT and ERK, and a reactivation of mTOR. Furthermore, the resistant cells exhibited strong upregulation and activation of MET, indicating crosstalk between the FGFR1 and MET axes. The resistant cells also underwent a global decrease in promoter hypermethylation of the CpG islands. Finally, we observed clonal expansion of a pre-existing change in AKT1, leading to S266L substitution, within the kinase domain of AKT. Our results demonstrate that AR to FGFR1-TKI involves deep molecular changes that promote the activation of MET and AKT, coupled with common gene expression and DNA methylation profiles. The expansion of a substitution at AKT1 was the only shared genetic change, and this may have contributed to the AR.

Histone H1 depletion triggers an interferon response in cancer cells via activation of heterochromatic repeats.

Histone H1 has seven variants in human somatic cells and contributes to chromatin compaction and transcriptional regulation. Knock-down (KD) of each H1 variant in breast cancer cells results in altered gene expression and proliferation differently in a variant specific manner with H1.2 and H1.4 KDs being most deleterious. Here we show combined depletion of H1.2 and H1.4 has a strong deleterious effect resulting in a strong interferon (IFN) response, as evidenced by an up-regulation of many IFN-stimulated genes (ISGs) not seen in individual nor in other combinations of H1 variant KDs. Although H1 participates to repress ISG promoters, IFN activation upon H1.2 and H1.4 KD is mainly generated through the activation of the IFN response by cytosolic nucleic acid receptors and IFN synthesis, and without changes in histone modifications at induced ISG promoters. H1.2 and H1.4 co-KD also promotes the appearance of accessibility sites genome wide and, particularly, at satellites and other repeats. The IFN response may be triggered by the expression of noncoding RNA generated from heterochromatic repeats or endogenous retroviruses upon H1 KD. In conclusion, redundant H1-mediated silencing of heterochromatin is important to maintain cell homeostasis and to avoid an unspecific IFN response.

Wolf-Hirschhorn Syndrome Candidate 1 Is Necessary for Correct Hematopoietic and B Cell Development.

Immunodeficiency is one of the most important causes of mortality associated with Wolf-Hirschhorn syndrome (WHS), a severe rare disease originated by a deletion in chromosome 4p. The WHS candidate 1 (WHSC1) gene has been proposed as one of the main genes responsible for many of the alterations in WHS, but its mechanism of action is still unknown. Here, we present in vivo genetic evidence showing that Whsc1 plays an important role at several points of hematopoietic development. Particularly, our results demonstrate that both differentiation and function of Whsc1-deficient B cells are impaired at several key developmental stages due to profound molecular defects affecting B cell lineage specification, commitment, fitness, and proliferation, demonstrating a causal role for WHSC1 in the immunodeficiency of WHS patients.

Genomic and Molecular Screenings Identify Different Mechanisms for Acquired Resistance to MET Inhibitors in Lung Cancer Cells.

The development of resistance to tyrosine kinase inhibitors (TKI) limits the long-term efficacy of cancer treatments involving them. We aimed to understand the mechanisms that underlie acquired resistance (AR) to MET inhibitors in lung cancer. EBC1 cells, which have MET amplification and are sensitive to TKIs against MET, were used to generate multiple clones with AR to a MET-TKI. Whole-exome sequencing, RNA sequencing, and global DNA methylation analysis were used to scrutinize the genetic and molecular characteristics of the resistant cells. AR to the MET-TKI involved changes common to all resistant cells, that is, phenotypic modifications, specific changes in gene expression, and reactivation of AKT, ERK, and mTOR. The gene expression, global DNA methylation, and mutational profiles distinguished at least two groups of resistant cells. In one of these, the cells have acquired sensitivity to erlotinib, concomitantly with mutations of the KIRREL, HDAC11, HIATL1, and MAPK1IP1L genes, among others. In the other group, some cells have acquired inactivation of neurofibromatosis type 2 (NF2) concomitantly with strong overexpression of NRG1 and a mutational profile that includes changes in LMLN and TOMM34 Multiple independent and simultaneous strategies lead to AR to the MET-TKIs in lung cancer cells. The acquired sensitivity to erlotinib supports the known crosstalk between MET and the HER family of receptors. For the first time, we show inactivation of NF2 during acquisition of resistance to MET-TKI that may explain the refractoriness to erlotinib in these cells. Mol Cancer Ther; 16(7); 1366-76. ©2017 AACR.

A Comparison of RNA-Seq Results from Paired Formalin-Fixed Paraffin-Embedded and Fresh-Frozen Glioblastoma Tissue Samples.

The molecular classification of glioblastoma (GBM) based on gene expression might better explain outcome and response to treatment than clinical factors. Whole transcriptome sequencing using next-generation sequencing platforms is rapidly becoming accepted as a tool for measuring gene expression for both research and clinical use. Fresh frozen (FF) tissue specimens of GBM are difficult to obtain since tumor tissue obtained at surgery is often scarce and necrotic and diagnosis is prioritized over freezing. After diagnosis, leftover tissue is usually stored as formalin-fixed paraffin-embedded (FFPE) tissue. However, RNA from FFPE tissues is usually degraded, which could hamper gene expression analysis. We compared RNA-Seq data obtained from matched pairs of FF and FFPE GBM specimens. Only three FFPE out of eleven FFPE-FF matched samples yielded informative results. Several quality-control measurements showed that RNA from FFPE samples was highly degraded but maintained transcriptomic similarities to RNA from FF samples. Certain issues regarding mutation analysis and subtype prediction were detected. Nevertheless, our results suggest that RNA-Seq of FFPE GBM specimens provides reliable gene expression data that can be used in molecular studies of GBM if the RNA is sufficiently preserved.

A comprehensive assessment of somatic mutation detection in cancer using whole-genome sequencing.

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.