SQUID: transcriptomic structural variation detection from RNA-seq

Epigenetic alterations affecting hematopoietic regulatory networks as drivers of mixed myeloid/lymphoid leukemia

Leukemias with ambiguous lineage comprise several loosely defined entities, often without a clear mechanistic basis. Here, we extensively profile the epigenome and transcriptome of a subgroup of such leukemias with CpG Island Methylator Phenotype. These leukemias exhibit comparable hybrid myeloid/lymphoid epigenetic landscapes, yet heterogeneous genetic alterations, suggesting they are defined by their shared epigenetic profile rather than common genetic lesions. Gene expression enrichment reveals similarity with early T-cell precursor acute lymphoblastic leukemia and a lymphoid progenitor cell of origin. In line with this, integration of differential DNA methylation and gene expression shows widespread silencing of myeloid transcription factors. Moreover, binding sites for hematopoietic transcription factors, including CEBPA, SPI1 and LEF1, are uniquely inaccessible in these leukemias. Hypermethylation also results in loss of CTCF binding, accompanied by changes in chromatin interactions involving key transcription factors. In conclusion, epigenetic dysregulation, and not genetic lesions, explains the mixed phenotype of this group of leukemias with ambiguous lineage. The data collected here constitute a useful and comprehensive epigenomic reference for subsequent studies of acute myeloid leukemias, T-cell acute lymphoblastic leukemias and mixed-phenotype leukemias.

SFyNCS detects oncogenic fusions involving non-coding sequences in cancer

Fusion genes are well-known cancer drivers. However, most known oncogenic fusions are protein-coding, and very few involve non-coding sequences due to lack of suitable detection tools. We develop SFyNCS to detect fusions of both protein-coding genes and non-coding sequences from transcriptomic sequencing data. The main advantage of this study is that we use somatic structural variations detected from genomic data to validate fusions detected from transcriptomic data. This allows us to comprehensively evaluate various fusion detection and filtering strategies and parameters. We show that SFyNCS has superior sensitivity and specificity over existing algorithms through extensive benchmarking in cancer cell lines and patient samples. We then apply SFyNCS to 9565 tumor samples across 33 tumor types in The Cancer Genome Atlas cohort and detect a total of 165,139 fusions. Among them, 72% of the fusions involve non-coding sequences. We find a long non-coding RNA to recurrently fuse with various oncogenes in 3% of prostate cancers. In addition, we discover fusions involving two non-coding RNAs in 32% of dedifferentiated liposarcomas and experimentally validated the oncogenic functions in mouse model.

On de novo Bridging Paired-end RNA-seq Data

The high-throughput short-reads RNA-seq protocols often produce paired-end reads, with the middle portion of the fragments being unsequenced. We explore if the full-length fragments can be computationally reconstructed from the sequenced two ends in the absence of the reference genome-a problem here we refer to as de novo bridging. Solving this problem provides longer, more informative RNA-seq reads, and benefits downstream RNA-seq analysis such as transcript assembly, expression quantification, and splicing differential analysis. However, de novo bridging is a challenging and complicated task owing to alternative splicing, transcript noises, and sequencing errors. It remains unclear if the data provides sufficient information for accurate bridging, let alone efficient algorithms that determine the true bridges. Methods have been proposed to bridge paired-end reads in the presence of reference genome (called reference-based bridging), but the algorithms are far away from scaling for de novo bridging as the underlying compacted de Bruijn graph (cdBG) used in the latter task often contains millions of vertices and edges. We designed a new truncated Dijkstra's algorithm for this problem, and proposed a novel algorithm that reuses the shortest path tree to avoid running the truncated Dijkstra's algorithm from scratch for all vertices for further speeding up. These innovative techniques result in scalable algorithms that can bridge all paired-end reads in a cdBG with millions of vertices. Our experiments showed that paired-end RNA-seq reads can be accurately bridged to a large extent. The resulting tool is freely available at https://github.com/Shao-Group/rnabridge-denovo.

Longitudinal characterization of primary osteosarcoma and derived subcutaneous and orthotopic relapsed patient-derived xenograft models

Osteosarcoma is a rare bone cancer in adolescents and young adults with a dismal prognosis because of metastatic disease and chemoresistance. Despite multiple clinical trials, no improvement in outcome has occurred in decades. There is an urgent need to better understand resistant and metastatic disease and to generate in vivo models from relapsed tumors. We developed eight new patient-derived xenograft (PDX) subcutaneous and orthotopic/paratibial models derived from patients with recurrent osteosarcoma and compared the genetic and transcriptomic landscapes of the disease progression at diagnosis and relapse with the matching PDX. Whole exome sequencing showed that driver and copy-number alterations are conserved from diagnosis to relapse, with the emergence of somatic alterations of genes mostly involved in DNA repair, cell cycle checkpoints, and chromosome organization. All PDX patients conserve most of the genetic alterations identified at relapse. At the transcriptomic level, tumor cells maintain their ossification, chondrocytic, and trans-differentiation programs during progression and implantation in PDX models, as identified at the radiological and histological levels. A more complex phenotype, like the interaction with immune cells and osteoclasts or cancer testis antigen expression, seemed conserved and was hardly identifiable by histology. Despite NSG mouse immunodeficiency, four of the PDX models partially reconstructed the vascular and immune-microenvironment observed in patients, among which the macrophagic TREM2/TYROBP axis expression, recently linked to immunosuppression. Our multimodal analysis of osteosarcoma progression and PDX models is a valuable resource to understand resistance and metastatic spread mechanisms, as well as for the exploration of novel therapeutic strategies for advanced osteosarcoma.

SFyNCS detects oncogenic fusions involving non-coding sequences in cancer

Fusion genes are well-known cancer drivers. However, very few known oncogenic fusions involve non-coding sequences. We develop SFyNCS with superior performance to detect fusions of both protein-coding genes and non-coding sequences from transcriptomic sequencing data. We validate fusions using somatic structural variations detected from the genomes. This allows us to comprehensively evaluate various fusion detection and filtering strategies and parameters. We detect 165,139 fusions in 9,565 tumor samples across 33 tumor types in the Cancer Genome Atlas cohort. Among them, 72% of the fusions involve non-coding sequences and many are recurrent. We discover two long non-coding RNAs recurrently fused with various partner genes in 32% of dedifferentiated liposarcomas and experimentally validated the oncogenic functions in mouse model.

Chimeric RNAs Discovered by RNA Sequencing and Their Roles in Cancer and Rare Genetic Diseases

Chimeric RNAs are transcripts that are generated by gene fusion and intergenic splicing events, thus comprising nucleotide sequences from different parental genes. In the past, Northern blot analysis and RT-PCR were used to detect chimeric RNAs. However, they are low-throughput and can be time-consuming, labor-intensive, and cost-prohibitive. With the development of RNA-seq and transcriptome analyses over the past decade, the number of chimeric RNAs in cancer as well as in rare inherited diseases has dramatically increased. Chimeric RNAs may be potential diagnostic biomarkers when they are specifically expressed in cancerous cells and/or tissues. Some chimeric RNAs can also play a role in cell proliferation and cancer development, acting as tools for cancer prognosis, and revealing new insights into the cell origin of tumors. Due to their abilities to characterize a whole transcriptome with a high sequencing depth and intergenically identify spliced chimeric RNAs produced with the absence of chromosomal rearrangement, RNA sequencing has not only enhanced our ability to diagnose genetic diseases, but also provided us with a deeper understanding of these diseases. Here, we reviewed the mechanisms of chimeric RNA formation and the utility of RNA sequencing for discovering chimeric RNAs in several types of cancer and rare inherited diseases. We also discussed the diagnostic, prognostic, and therapeutic values of chimeric RNAs.

Frequency and Genomic Aspects of Intrinsic Resistance to Vismodegib in Locally Advanced Basal Cell Carcinoma

PURPOSE

Vismodegib is approved for the treatment of locally advanced basal cell carcinoma (laBCC), but some cases demonstrate intrinsic resistance (IR) to the drug. We sought to assess the frequency of IR to vismodegib in laBCC and its underlying genomic mechanisms.

EXPERIMENTAL DESIGN

Response to vismodegib was evaluated in a cohort of 148 laBCC patients. Comprehensive genomic and transcriptomic profiling was performed in a subset of five intrinsically resistant BCC (IR-BCC).

RESULTS

We identified that IR-BCC represents 6.1% of laBCC in the studied cohort. Prior treatment with chemotherapy was associated with IR. Genetic events that were previously associated with acquired resistance (AR) in BCC or medulloblastoma were observed in three out of five IR-BCC. However, IR-BCCs were distinct by highly rearranged polyploid genomes. Functional analyses identified hyperactivation of the HIPPO-YAP and WNT pathways at RNA and protein levels in IR-BCC. In vitro assay on the BCC cell line further confirmed that YAP1 overexpression increases the cell proliferation rate.

CONCLUSIONS

IR to vismodegib is a rare event in laBCC. IR-BCCs frequently harbor resistance mutations in the Hh pathway, but also are characterized by hyperactivation of the HIPPO-YAP and WNT pathways.

A guide for the diagnosis of rare and undiagnosed disease: beyond the exome

Rare diseases affect 30 million people in the USA and more than 300-400 million worldwide, often causing chronic illness, disability, and premature death. Traditional diagnostic techniques rely heavily on heuristic approaches, coupling clinical experience from prior rare disease presentations with the medical literature. A large number of rare disease patients remain undiagnosed for years and many even die without an accurate diagnosis. In recent years, gene panels, microarrays, and exome sequencing have helped to identify the molecular cause of such rare and undiagnosed diseases. These technologies have allowed diagnoses for a sizable proportion (25-35%) of undiagnosed patients, often with actionable findings. However, a large proportion of these patients remain undiagnosed. In this review, we focus on technologies that can be adopted if exome sequencing is unrevealing. We discuss the benefits of sequencing the whole genome and the additional benefit that may be offered by long-read technology, pan-genome reference, transcriptomics, metabolomics, proteomics, and methyl profiling. We highlight computational methods to help identify regionally distant patients with similar phenotypes or similar genetic mutations. Finally, we describe approaches to automate and accelerate genomic analysis. The strategies discussed here are intended to serve as a guide for clinicians and researchers in the next steps when encountering patients with non-diagnostic exomes.

MINTIE: identifying novel structural and splice variants in transcriptomes using RNA-seq data

Calling fusion genes from RNA-seq data is well established, but other transcriptional variants are difficult to detect using existing approaches. To identify all types of variants in transcriptomes we developed MINTIE, an integrated pipeline for RNA-seq data. We take a reference-free approach, combining de novo assembly of transcripts with differential expression analysis to identify up-regulated novel variants in a case sample. We compare MINTIE with eight other approaches, detecting > 85% of variants while no other method is able to achieve this. We posit that MINTIE will be able to identify new disease variants across a range of disease types.

Identification and Characterization of Non-Coding RNAs in Thymoma

BACKGROUND Thymoma is the most common tumor of the anterior mediastinum, and can be caused by infrequent malignancies arising from the epithelial cells of the thymus. Unfortunately, blood-based diagnostic markers are not currently available. High-throughput sequencing technologies, such as RNA-seq with next-generation sequencing, have facilitated the detection and characterization of both coding and non-coding RNAs (ncRNAs), which play significant roles in genomic regulation, transcriptional and post-transcriptional regulation, and imprinting and epigenetic modification. The knowledge about fusion genes and ncRNAs in thymomas is scarce. MATERIAL AND METHODS For this study, we gathered large-scale RNA-seq data belonging to samples from 25 thymomas and 25 healthy thymus specimens and analyzed them to identify fusion genes, lncRNAs, and miRNAs. RESULTS We found 21 fusion genes, including KMT2A-MAML2, HADHB-REEP1, COQ3-CGA, MCM4-SNTB1, and IFT140-ACTN4, as the most frequent and significant in thymomas. We also detected 65 differentially-expressed lncRNAs in thymomas, including AFAP1-AS1, LINC00324, ADAMTS9-AS1, VLDLR-AS1, LINC00968, and NEAT1, that have been validated with the TCGA database. Moreover, we identified 1695 miRNAs from small RNA-seq data that were overexpressed in thymomas. Our network analysis of the lncRNA-mRNA-miRNA regulation axes identified a cluster of miRNAs upregulated in thymomas, that can trigger the expression of target protein-coding genes, and lead to the disruption of several biological pathways, including the PI3K-Akt signaling pathway, FoxO signaling pathway, and HIF-1 signaling pathway. CONCLUSIONS Our results show that overexpression of this miRNA cluster activates PI3K-Akt, FoxO, HIF-1, and Rap-1 signaling pathways, suggesting pathway inhibitors may be therapeutic candidates against thymoma.

Full-Length Transcriptome Analysis of Plasmodium falciparum by Single-Molecule Long-Read Sequencing

Malaria, an infectious disease caused by Plasmodium parasites, still accounts for amounts of deaths annually in last decades. Despite the significance of Plasmodium falciparum as a model organism of malaria parasites, our understanding of gene expression of this parasite remains largely elusive since lots of progress on its genome and transcriptome are based on assembly with short sequencing reads. Herein, we report the new version of transcriptome dataset containing all full-length transcripts over the whole asexual blood stages by adopting a full-length sequencing approach with optimized experimental conditions of cDNA library preparation. We have identified a total of 393 alternative splicing (AS) events, 3,623 long non-coding RNAs (lncRNAs), 1,555 alternative polyadenylation (APA) events, 57 transcription factors (TF), 1,721 fusion transcripts in P. falciparum. Furthermore, the shotgun proteome was performed to validate the full-length transcriptome of P. falciparum. More importantly, integration of full-length transcriptomic and proteomic data identified 160 novel small proteins in lncRNA regions. Collectively, this full-length transcriptome dataset with high quality and accuracy and the shotgun proteome analyses shed light on the complex gene expression in malaria parasites and provide a valuable resource for related functional and mechanistic researches on P. falciparum genes.

Loss of Fbxw7 Impairs Development of and Induces Heterogeneous Tumor Formation in the Mouse Mammary Gland

Fbxw7 is an F-box protein that contributes to regulation of cell proliferation and cell fate determination as well as to tumor suppression in various tissues. In this study, we generated mice with mammary gland-specific ablation of Fbxw7 (Blg-Cre/Fbxw7 ^F/F mice) and found that most neonates born to mutant dams die soon after birth as a result of defective maternal lactation. The mammary gland of mutant dams was markedly atrophic and manifested both excessive cell proliferation and apoptosis in association with the accumulation of Notch1 and p63. Despite the hypoplastic nature of the mutant mammary gland, Blg-Cre/Fbxw7 ^F/F mice spontaneously developed mammary tumors that resembled basal-like carcinoma with marked intratumoral heterogeneity. Additional inactivation of Trp53 in Blg-Cre/Fbxw7 ^F/F mice further promoted onset and development of mammary tumors, suggesting that spontaneous mutation of Trp53 may facilitate transition of hypoplastic mammary lesions to aggressive cancer in mice lacking Fbxw7. RNA-sequencing analysis of epithelial- and mesenchymal-like cell lines from a Blg-Cre/Fbxw7 ^F/F mouse tumor revealed an increased mutation rate and structural alterations in the tumor and differential expression of upstream transcription factors including known targets of Fbxw7. Together, our results implicate Fbxw7 in the regulation of cell differentiation and in tumor suppression in the mammary gland. Loss of Fbxw7 increases mutation rate and chromosome instability, activates signaling pathways governed by transcription factors regulated by Fbxw7, and triggers the development of mammary tumors with prominent heterogeneity. SIGNIFICANCE: Mammary gland-specific ablation of Fbxw7 in mice results in defective gland development and spontaneous mammary tumor formation reminiscent of human basal-like carcinoma with intratumoral heterogeneity. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/24/5515/F1.large.jpg.

The mutational landscape of the SCAN-B real-world primary breast cancer transcriptome

Breast cancer is a disease of genomic alterations, of which the panorama of somatic mutations and how these relate to subtypes and therapy response is incompletely understood. Within SCAN-B (ClinicalTrials.gov: NCT02306096), a prospective study elucidating the transcriptomic profiles for thousands of breast cancers, we developed a RNA-seq pipeline for detection of SNVs/indels and profiled a real-world cohort of 3,217 breast tumors. We describe the mutational landscape of primary breast cancer viewed through the transcriptome of a large population-based cohort and relate it to patient survival. We demonstrate that RNA-seq can be used to call mutations in genes such as PIK3CA, TP53, and ERBB2, as well as the status of molecular pathways and mutational burden, and identify potentially druggable mutations in 86.8% of tumors. To make this rich dataset available for the research community, we developed an open source web application, the SCAN-B MutationExplorer (http://oncogenomics.bmc.lu.se/MutationExplorer). These results add another dimension to the use of RNA-seq as a clinical tool, where both gene expression- and mutation-based biomarkers can be interrogated in real-time within 1 week of tumor sampling.

Fusion-Bloom: fusion detection in assembled transcriptomes

Summary

Presence or absence of gene fusions is one of the most important diagnostic markers in many cancer types. Consequently, fusion detection methods using various genomics data types, such as RNA sequencing (RNA-seq) are valuable tools for research and clinical applications. While information-rich RNA-seq data have proven to be instrumental in discovery of a number of hallmark fusion events, bioinformatics tools to detect fusions still have room for improvement. Here, we present Fusion-Bloom, a fusion detection method that leverages recent developments in de novo transcriptome assembly and assembly-based structural variant calling technologies (RNA-Bloom and PAVFinder, respectively). We benchmarked Fusion-Bloom against the performance of five other state-of-the-art fusion detection tools using multiple datasets. Overall, we observed Fusion-Bloom to display a good balance between detection sensitivity and specificity. We expect the tool to find applications in translational research and clinical genomics pipelines.

Availability and implementation

Fusion-Bloom is implemented as a UNIX Make utility, available at https://github.com/bcgsc/pavfinder and released under a Creative Commons License (Attribution 4.0 International), as described at http://creativecommons.org/licenses/by/4.0/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Chromosomal translocation-derived aberrant Rab22a drives metastasis of osteosarcoma

Osteosarcoma is a type of aggressive malignant bone tumour that frequently metastasizes to lungs, resulting in poor prognosis. However, the molecular mechanisms of lung metastasis of osteosarcoma remain poorly understood. Here we identify exon-intron fusion genes in osteosarcoma cell lines and tissues. These fusion genes are derived from chromosomal translocations that juxtapose the coding region for amino acids 1-38 of Rab22a (Rab22a1-38) with multiple inverted introns and untranslated regions of chromosome 20. The resulting translation products, designated Rab22a-NeoFs, acquire the ability to drive lung metastasis of osteosarcoma. The Rab22a1-38 moiety governs the function of Rab22a-NeoFs by binding to SmgGDS-607, a GTP-GDP exchange factor of RhoA. This association facilitates the release of GTP-bound RhoA from SmgGDS-607, which induces increased activity of RhoA and promotes metastasis. Disrupting the interaction between Rab22a-NeoF1 and SmgGDS-607 with a synthetic peptide prevents lung metastasis in an orthotopic model of osteosarcoma. Our findings may provide a promising strategy for a subset of osteosarcoma patients with lung metastases.

Detecting transcriptomic structural variants in heterogeneous contexts via the Multiple Compatible Arrangements Problem

BACKGROUND

Transcriptomic structural variants (TSVs)-large-scale transcriptome sequence change due to structural variation - are common in cancer. TSV detection from high-throughput sequencing data is a computationally challenging problem. Among all the confounding factors, sample heterogeneity, where each sample contains multiple distinct alleles, poses a critical obstacle to accurate TSV prediction.

RESULTS

To improve TSV detection in heterogeneous RNA-seq samples, we introduce the Multiple Compatible Arrangements Problem (MCAP), which seeks k genome arrangements that maximize the number of reads that are concordant with at least one arrangement. This models a heterogeneous or diploid sample. We prove that MCAP is NP-complete and provide a 1 4 -approximation algorithm for k = 1 and a 3 4 -approximation algorithm for the diploid case ( k = 2 ) assuming an oracle for k = 1 . Combining these, we obtain a 3 16 -approximation algorithm for MCAP when k = 2 (without an oracle). We also present an integer linear programming formulation for general k. We characterize the conflict structures in the graph that require k > 1 alleles to satisfy read concordancy and show that such structures are prevalent.

CONCLUSIONS

We show that the solution to MCAP accurately addresses sample heterogeneity during TSV detection. Our algorithms have improved performance on TCGA cancer samples and cancer cell line samples compared to a TSV calling tool, SQUID. The software is available at https://github.com/Kingsford-Group/diploidsquid.

Structural variant identification and characterization

Structural variant (SV) differences between human genomes can cause germline and mosaic disease as well as inter-individual variation. De-regulation of accurate DNA repair and genomic surveillance mechanisms results in a large number of SVs in cancer. Analysis of the DNA sequences at SV breakpoints can help identify pathways of mutagenesis and regions of the genome that are more susceptible to rearrangement. Large-scale SV analyses have been enabled by high-throughput genome-level sequencing on humans in the past decade. These studies have shed light on the mechanisms and prevalence of complex genomic rearrangements. Recent advancements in both sequencing and other mapping technologies as well as calling algorithms for detection of genomic rearrangements have helped propel SV detection into population-scale studies, and have begun to elucidate previously inaccessible regions of the genome. Here, we discuss the genomic organization of simple and complex SVs, the molecular mechanisms of their formation, and various ways to detect them. We also introduce methods for characterizing SVs and their consequences on human genomes.

Structural variant calling: the long and the short of it

Recent research into structural variants (SVs) has established their importance to medicine and molecular biology, elucidating their role in various diseases, regulation of gene expression, ethnic diversity, and large-scale chromosome evolution-giving rise to the differences within populations and among species. Nevertheless, characterizing SVs and determining the optimal approach for a given experimental design remains a computational and scientific challenge. Multiple approaches have emerged to target various SV classes, zygosities, and size ranges. Here, we review these approaches with respect to their ability to infer SVs across the full spectrum of large, complex variations and present computational methods for each approach.

NCLcomparator: systematically post-screening non-co-linear transcripts (circular, trans-spliced, or fusion RNAs) identified from various detectors

BACKGROUND

Non-co-linear (NCL) transcripts consist of exonic sequences that are topologically inconsistent with the reference genome in an intragenic fashion (circular or intragenic trans-spliced RNAs) or in an intergenic fashion (fusion or intergenic trans-spliced RNAs). On the basis of RNA-seq data, numerous NCL event detectors have been developed and detected thousands of NCL events in diverse species. However, there are great discrepancies in the identification results among detectors, indicating a considerable proportion of false positives in the detected NCL events. Although several helpful guidelines for evaluating the performance of NCL event detectors have been provided, a systematic guideline for measurement of NCL events identified by existing tools has not been available.

RESULTS

We develop a software, NCLcomparator, for systematically post-screening the intragenic or intergenic NCL events identified by various NCL detectors. NCLcomparator first examine whether the input NCL events are potentially false positives derived from ambiguous alignments (i.e., the NCL events have an alternative co-linear explanation or multiple matches against the reference genome). To evaluate the reliability of the identified NCL events, we define the NCL score (NCLscore) based on the variation in the number of supporting NCL junction reads identified by the tools examined. Of the input NCL events, we show that the ambiguous alignment-derived events have relatively lower NCLscore values than the other events, indicating that an NCL event with a higher NCLscore has a higher level of reliability. To help selecting highly expressed NCL events, NCLcomparator also provides a series of useful measurements such as the expression levels of the detected NCL events and their corresponding host genes and the junction usage of the co-linear splice junctions at both NCL donor and acceptor sites.

CONCLUSION

NCLcomparator provides useful guidelines, with the input of identified NCL events from various detectors and the corresponding paired-end RNA-seq data only, to help users selecting potentially high-confidence NCL events for further functional investigation. The software thus helps to facilitate future studies into NCL events, shedding light on the fundamental biology of this important but understudied class of transcripts. NCLcomparator is freely accessible at https://github.com/TreesLab/NCLcomparator .