Transcription-mediated gene fusion in the human genome

RTCpredictor: identification of read-through chimeric RNAs from RNA sequencing data

Read-through chimeric RNAs are being recognized as a means to expand the functional transcriptome and contribute to cancer tumorigenesis when mis-regulated. However, current software tools often fail to predict them. We have developed RTCpredictor, utilizing a fast ripgrep tool to search for all possible exon-exon combinations of parental gene pairs. We also added exonic variants allowing searches containing common SNPs. To our knowledge, it is the first read-through chimeric RNA specific prediction method that also provides breakpoint coordinates. Compared with 10 other popular tools, RTCpredictor achieved high sensitivity on a simulated and three real datasets. In addition, RTCpredictor has less memory requirements and faster execution time, making it ideal for applying on large datasets.

Sex-Biased Transcription Expression of Vitellogenins Reveals Fusion Gene and MicroRNA Regulation in the Sea Louse Caligus rogercresseyi

The caligid ectoparasite, Caligus rogercresseyi, is one of the main concerns in the Chilean salmon industry. The molecular mechanisms displayed by the parasite during the reproductive process represent an opportunity for developing novel control strategies. Vitellogenin is a multifunctional protein recognized as a critical player in several crustaceans' biological processes, including reproduction, embryonic development, and immune response. This study aimed to characterize the C. rogercresseyi vitellogenins, including discovering novel transcripts and regulatory mechanisms associated with microRNAs. Herein, vitellogenin genes were identified by homology analysis using the reference sea louse genome, transcriptome database, and arthropods vitellogenin-protein database. The validation of expression transcripts was conducted by RNA nanopore sequencing technology. Moreover, fusion gene profiling, miRNA target analysis, and functional validation were performed using luciferase assay. Six putative vitellogenin genes were identified in the C. rogercresseyi genome with high homology with other copepods vitellogenins. Furthermore, miR-996 showed a putative role in regulating the Cr_Vitellogenin1 gene, which is highly expressed in females. Moreover, vitellogenin-fusion genes were identified in adult stages and highly regulated in males, demonstrating sex-related expression patterns. In females, the identified fusion genes merged with several non-vitellogenin genes involved in biological processes of ribosome assembly, BMP signaling pathway, and biosynthetic processes. This study reports the genome array of vitellogenins in C. rogercresseyi for the first time, revealing the putative role of fusion genes and miRNA regulation in sea lice biology.

Uncovering the Localization and Function of a Novel Read-Through Transcript 'TOMM40-APOE'

Recent advancements in genome analysis technology have revealed the presence of read-through transcripts in which transcription continues by skipping the polyA signal. We here identified and characterized a new read-through transcript, TOMM40-APOE. With cDNA amplification from THP-1 cells, the TOMM40-APOE3 product was successfully generated. We also generated TOMM40-APOE4, another isoform, by introducing point mutations. Notably, while APOE3 and APOE4 exhibited extracellular secretion, both TOMM40-APOE3 and TOMM40-APOE4 were localized exclusively to the mitochondria. But functionally, they did not affect mitochondrial membrane potential. Cell death induction studies illustrated increased cell death with TOMM40-APOE3 and TOMM40-APOE4, and we did not find any difference in cellular function between the two isoforms. These findings indicated that the new mitochondrial protein TOMM40-APOE has cell toxic ability.

The fusion gene hsf5-rnf43 in Nile tilapia: A potential regulator in the maintenance of testis function and sexual differentiation

We identified that the genes heat shock transcription factor 5 (hsf5) and ring finger protein 43 (rnf43) happened fusion in Nile tilapia (Oreochromis niloticus), called hsf5-rnf43, and provided the characteristic and functional analysis of hsf5-rnf43 gene in fish for the first time. Analysis of spatiotemporal expression showed that hsf5-rnf43 was specifically expressed in the testis and located in primary spermatocytes of adult Nile tilapia and gradually increased during testis development from 5 to 180 days after hatching. We also found DNA methylation regulated sex-biased expression of hsf5-rnf43 in the early development of Nile tilapia, and was affected by high temperature during the thermosensitive period of Nile tilapia sex differentiation. Therefore, we first reported that the fusion gene hsf5-rnf43 was sex-biased expressed in the testis regulated by DNA methylation and affected by high temperature, which may be involved in the maintenance of testis function and sex differentiation of Nile tilapia.

Pattern recognition in the landscape of seemingly random chimeric transcripts

The molecular and functional diversity generated by chimeric transcripts (CTs) that are derived from two genes is indicated to contribute to tumor cell survival. Several gaps yet exist. The present research is a systematic study of the spectrum of CTs identified in RNA sequencing datasets of 160 ovarian cancer samples in the The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov). Structural annotation revealed complexities emerging from chromosomal localization of partner genes, differential splicing and inclusion of regulatory, untranslated regions. Identification of phenotype-specific associations further resolved a dynamically modulated mesenchymal signature during transformation. On an evolutionary background, protein-coding CTs were indicated to be highly conserved, while non-coding CTs may have evolved more recently. We also realized that the current premise postulating structural alterations or neighbouring gene readthrough generating CTs is not valid in instances wherein the parental genes are genomically distanced. In addressing this lacuna, we identified the essentiality of specific spatiotemporal arrangements mediated gene proximities in 3D space for the generation of CTs. All these features together suggest non-random mechanisms towards increasing the molecular diversity in a cell through chimera formation either in parallel or with cross-talks with the indigenous regulatory network.

The power of "controllers": Transposon-mediated duplicated genes evolve towards neofunctionalization

Since the discovery of the first transposon by Dr. Barbara McClintock, the prevalence and diversity of transposable elements (TEs) have been gradually recognized. As fundamental genetic components, TEs drive organismal evolution not only by contributing functional sequences (e.g., regulatory elements or "controllers" as phrased by Dr. McClintock) but also by shuffling genomic sequences. In the latter respect, TE-mediated gene duplications have contributed to the origination of new genes and attracted extensive interest. In response to the development of this field, we herein attempt to provide an overview of TE-mediated duplication by focusing on common rules emerging across duplications generated by different TE types. Specifically, despite the huge divergence of transposition machinery across TEs, we identify three common features of various TE-mediated duplication mechanisms, including end bypass, template switching, and recurrent transposition. These three features lead to one common functional outcome, namely, TE-mediated duplicates tend to be subjected to exon shuffling and neofunctionalization. Therefore, the intrinsic properties of the mutational mechanism constrain the evolutionary trajectories of these duplicates. We finally discuss the future of this field including an in-depth characterization of both the duplication mechanisms and functions of TE-mediated duplicates.

Targeting pre-mRNA splicing in cancers: roles, inhibitors, and therapeutic opportunities

Accumulating evidence has indicated that pre-mRNA splicing plays critical roles in a variety of physiological processes, including development of multiple diseases. In particular, alternative splicing is profoundly involved in cancer progression through abnormal expression or mutation of splicing factors. Small-molecule splicing modulators have recently attracted considerable attention as a novel class of cancer therapeutics, and several splicing modulators are currently being developed for the treatment of patients with various cancers and are in the clinical trial stage. Novel molecular mechanisms modulating alternative splicing have proven to be effective for treating cancer cells resistant to conventional anticancer drugs. Furthermore, molecular mechanism-based combination strategies and patient stratification strategies for cancer treatment targeting pre-mRNA splicing must be considered for cancer therapy in the future. This review summarizes recent progress in the relationship between druggable splicing-related molecules and cancer, highlights small-molecule splicing modulators, and discusses future perspectives of splicing modulation for personalized and combination therapies in cancer treatment.

Preleukemic Fusion Genes Induced via Ionizing Radiation

Although the prevalence of leukemia is increasing, the agents responsible for this increase are not definitely known. While ionizing radiation (IR) was classified as a group one carcinogen by the IARC, the IR-induced cancers, including leukemia, are indistinguishable from those that are caused by other factors, so the risk estimation relies on epidemiological data. Several epidemiological studies on atomic bomb survivors and persons undergoing IR exposure during medical investigations or radiotherapy showed an association between radiation and leukemia. IR is also known to induce chromosomal translocations. Specific chromosomal translocations resulting in preleukemic fusion genes (PFGs) are generally accepted to be the first hit in the onset of many leukemias. Several studies indicated that incidence of PFGs in healthy newborns is up to 100-times higher than childhood leukemia with the same chromosomal aberrations. Because of this fact, it has been suggested that PFGs are not able to induce leukemia alone, but secondary mutations are necessary. PFGs also have to occur in specific cell populations of hematopoetic stem cells with higher leukemogenic potential. In this review, we describe the connection between IR, PFGs, and cancer, focusing on recurrent PFGs where an association with IR has been established.

RTCpredictor: Identification of Read-Through Chimeric RNAs from RNA Sequencing Data

Read-through chimeric RNAs are gaining attention in cancer and other research fields, yet current tools often fail in predicting them. We have thus developed the first read-through chimeric RNA specific prediction method, RTCpredictor, utilizing a fast ripgrep algorithm to search for all possible exon-exon combinations of parental gene pairs. Compared with other ten popular tools, RTCpredictor achieved top performance on both simulated and real datasets. We randomly selected up to 30 candidate read-through chimeras predicted from each software method and experimentally validated a total of 109 read-throughs and on this set, RTCpredictor outperformed all the other methods. In addition, RTCpredictor ( https://github.com/sandybioteck/RTCpredictor ) has less memory requirements and faster execution time.

Knowing when to stop: Transcription termination on protein-coding genes by eukaryotic RNAPII

Gene expression is controlled in a dynamic and regulated manner to allow for the consistent and steady expression of some proteins as well as the rapidly changing production of other proteins. Transcription initiation has been a major focus of study because it is highly regulated. However, termination of transcription also plays an important role in controlling gene expression. Transcription termination on protein-coding genes is intimately linked with 3' end cleavage and polyadenylation of transcripts, and it generally results in the production of a mature mRNA that is exported from the nucleus. Termination on many non-coding genes can also result in the production of a mature transcript. Termination is dynamically regulated-premature termination and transcription readthrough occur in response to a number of cellular signals, and these can have varied consequences on gene expression. Here, we review eukaryotic transcription termination by RNA polymerase II (RNAPII), focusing on protein-coding genes.

Comprehensive analysis of full-length transcripts reveals novel splicing abnormalities and oncogenic transcripts in liver cancer

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. Although long-reads sequencing technologies would enable the sequencing of full-length transcripts, the data analysis is difficult. In this study, we developed an analysis pipeline named SPLICE and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC) and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 46,663 transcripts from the protein-coding genes in the HCCs and the matched non-cancerous livers, of which 5,366 (11.5%) were novel. A comparison of expression levels identified 9,933 differentially expressed transcripts (DETs) in 4,744 genes. Interestingly, 746 genes with DETs, including the LINE1-MET transcript, were not found by a gene-level analysis. We also found that fusion transcripts of transposable elements and hepatitis B virus (HBV) were overexpressed in HCCs. In vitro experiments on DETs showed that LINE1-MET and HBV-human transposable elements promoted cell growth. Furthermore, fusion gene detection showed novel recurrent fusion events that were not detected in the short-reads. These results suggest the efficiency of full-length transcriptome studies and the importance of splicing variants in carcinogenesis.

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. In this study, we developed an analysis pipeline named SPLICE for long-read transcriptome sequencing and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC), and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 5,366 novel transcripts and 9,933 differentially expressed transcripts in 4,744 genes between HCCs and non-cancerous livers. An analysis of hepatitis B virus (HBV) transcripts showed that fusion transcripts of the HBV gene and human transposable elements were overexpressed in HBV-infected HCCs. We also identified fusion genes that were not found in the short-reads. These results suggest that long-reads sequencing technologies provide a fuller understanding of cancer transcripts and that our method contributes to the analysis of transcriptome sequences by such technologies.

Long-Read Sequencing Annotation of the Transcriptome in DNA-PK Inactivated Cells

The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) with a Ku70/Ku80 heterodimer constitutes the intact DNA-PK kinase, which is an upstream component of the DNA repair machinery that signals the DNA damage, orchestrates the DNA repair, and serves to maintain genome integrity. Beyond its role in DNA damage repair, the DNA-PK kinase is also implicated in transcriptional regulation and RNA metabolism, with an illuminated impact on tumor progression and therapeutic responses. However, the efforts to identify DNA-PK regulated transcriptomes are limited by short-read sequencing to resolve the full complexity of the transcriptome. Therefore, we leveraged the PacBio Single Molecule, Real-Time (SMRT) Sequencing platform to study the transcriptome after DNA-PK inactivation to further underscore the importance of its role in diseases. Our analysis revealed additional novel transcriptome and complex gene structures in the DNA-PK inactivated cells, identifying 8,355 high-confidence new isoforms from 3,197 annotated genes and 523 novel genes. Among them, 380 lncRNAs were identified. We validated these findings using computational approaches and confirmatory transcript quantification with short-read sequencing. Several novel isoforms representing distinct splicing events have been validated through PCR experiments. Our analyses provide novel insights into DNA-PK function in transcriptome regulation and RNA metabolism.

Novel read-through fusion transcript Bcl2l2-Pabpn1 in glioblastoma cells

Read‐through fusion transcripts have recently been identified as chimeric RNAs and have since been linked to tumour growth in some cases. Many fusion genes generated by chromosomal rearrangements have been described in glioblastoma. However, read‐through fusion transcripts between neighbouring genes in glioblastoma remain unexplored. We performed paired‐end RNA‐seq of rat C6 glioma cells and normal cells and discovered a read‐through fusion transcript Bcl2l2‐Pabpn1 in which exon 3 of Bcl‐2‐like protein 2 (Bcl2l2) fused to exon 2 of Polyadenylate‐binding protein 1 (Pabpn1). This fusion transcript was found in both human glioblastoma and normal cells. Unlike other fusions reported in glioblastoma, Bcl2l2‐Pabpn1 appeared to result from RNA processing rather than genomic rearrangement. Bcl2l2‐Pabpn1 fusion transcript encoded a fusion protein with BH4, BCL and RRM domains. Functionally, Bcl2l2‐Pabpn1 knockdown by targeting its fusion junction decreased its expression, and suppressed cell proliferation, migration and invasion in vitro. Mechanistically, Bcl2l2‐Pabpn1 blocked Bax activity and activated PI3K/AKT pathway to promote glioblastoma progression. Together, our work characterized a glioblastoma‐associated Bcl2l2‐Pabpn1 fusion transcript shared by humans and rats.

CircRNA: An emerging star in the progression of glioma

Circular RNAs (circRNAs), a class of single-stranded noncoding RNAs with a covalently closed loop structure, are recognized as promising biomarkers and targets for diagnosing and treating dozens of diseases, especially cancers. CircRNAs are extremely stable, abundant and conserved and have tissue- or developmental stage-specific expression. Currently, the biogenesis and biological functions of circRNAs have been increasingly revealed with deep sequencing and bioinformatics. Studies have indicated that circRNAs are frequently expressed in brain tissues and that their expression levels change in different stages of neural development, suggesting that circRNAs may play an important role in diseases of the nervous system, such as glioma. However, because the biogenesis and functions of circRNAs do not depend on a single mechanism but are coregulated by multiple factors, it is necessary to further explore the underlying mechanisms. In this review, we summarized the classification, mechanisms of biogenesis and biological functions of circRNAs. Meanwhile, we emphatically expounded on the process of abnormal expression of circRNAs, methods used in circRNA research, and their effects on the malignant biological capabilities of glioma.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02696-w.

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.

Single-Cell Transcriptomics Reveals Splicing Features of Adult Neural Stem Cells in the Subventricular Zone

The central nervous system has enormously complex cellular diversity with hundreds of distinct cell types, yet alternative splicing features in single cells of important cell types at neurogenic regions are not well understood. By employing in silico analysis, we systematically identified 3,611 alternative splicing events from 1,908 genes in 28 single-cell transcriptomic data of adult mouse ependymal and subependymal regions, and found that single-cell RNA-seq has the advantage in uncovering rare splicing isoforms compared to bulk RNA-seq at the population level. We uncovered that the simultaneous presence of multiple isoforms from the same gene in a single cell is prevalent, and quiescent stem cells, activated stem cells, and neuroblast cells exhibit high heterogeneity of splicing variants. Furthermore, we also demonstrated the existence of novel bicistronic transcripts in quiescent stem cells.

Characterization of a Read-through Fusion Transcript, BCL2L2-PABPN1, Involved in Porcine Adipogenesis

cis-Splicing of adjacent genes (cis-SAGe) has been involved in multiple physiological and pathological processes in humans. However, to the best of our knowledge, there is no report of cis-SAGe in adipogenic regulation. In this study, a cis-SAGe product, BCL2L2–PABPN1 (BP), was characterized in fat tissue of pigs with RT-PCR and RACE method. BP is an in-frame fusion product composed of 333 aa and all the functional domains of both parents. BP is highly conserved among species and rich in splicing variants. BP was found to promote proliferation and inhibit differentiation of primary porcine preadipocytes. A total of 3074/44 differentially expressed mRNAs (DEmRs)/known miRNAs (DEmiRs) were identified in porcine preadipocytes overexpressing BP through RNA-Seq analysis. Both DEmRs and target genes of DEmiRs were involved in various fat-related pathways with MAPK and PI3K-Akt being the top enriched. PPP2CB, EGFR, Wnt5A and EHHADH were hub genes among the fat-related pathways identified. Moreover, ssc-miR-339-3p was found to be critical for BP regulating adipogenesis through integrated analysis of mRNA and miRNA data. The results highlight the role of cis-SAGe in adipogenesis and contribute to further revealing the mechanisms underlying fat deposition, which will be conductive to human obesity control.

Aquaporin (AQP) channels in the spiny dogfish, Squalus acanthias I: Characterization of AQP3 and AQP15 function and expression, and localization of the proteins in gill and spiral valve intestine

Complementary DNAs (cDNAs) for two aquaporin water channel genes (AQP3 and AQP15) were amplified cloned and sequenced to initiate this study. Northern blot analysis was carried out to confirm the mRNA sizes of these AQP genes with AQP3 mRNA bands exhibiting sizes of 1.2 and 1.6 k bases and AQP15 had a mRNA band of 2.1 k bases. Northern blot analysis was also performed on kidney and esophagus total RNA samples from fish acclimated to 75%, 100% or 120% seawater (SW). The level of AQP15 mRNA expression was shown to significantly decrease following salinity acclimation from 100 to 120% SW. An opposite but non-significantly different trend was observed for AQP3 mRNA levels. Full length cDNAs were then used to generate AQP3 and AQP15 mRNAs for microinjection into Xenopus oocytes. Both AQP3- and AQP15- microinjected oocytes exhibited significantly elevated apparent water permeability compared to control oocytes at neutral pH. The apparent water permeability was mercury-inhibitable, significantly so in the case of AQP3. AQP3 microinjected oocytes showed pH sensitivity in their apparent water permeability, showing a lack of permeability at acidic pH values. The Carboxyl-terminal derived amino acid sequences of AQP3 and AQP15 were used to generate rabbit affinity-purified polyclonal antibodies. Western blots with the antibodies showed a band of 31.3 kDa for AQP3 in the kidney, with minor bands at 26, 24 and 21 kDa. For AQP15 a band of 26 kDa was seen in gill and kidney. Fainter bands at 28 and 24 kDa were also seen in the kidney. There was also some higher molecular weight banding. None of the bands were seen when the antibodies were pre- blocked with their peptide antigens. Immunohistochemical localization studies were also performed in the gill and spiral valve intestine. In the gill, AQP15 antibody staining was seen sporadically in the membranes of surface epithelial cells of the secondary lamellae. Tyramide amplification of signals was employed in the spiral valve intestine. Tyramide-amplified AQP3 antibody staining was observed in the basal membrane of the invaginated epithelial cell layer of secondary intestinal folds in luminal surface of either the side wall of the spiral valve intestine or in internal valve tissue 'flaps'. For the AQP15 antibody, tyramide-amplified staining was instead found on the apical and to a lesser extent the lateral membranes of the same invaginated epithelial cell layer. The localization of AQP3 and AQP15 in the spiral valve intestine suggests that a trans-cellular water absorption pathway may exist in this tissue.

Gene regulatory network analysis defines transcriptome landscape with alternative splicing of human umbilical vein endothelial cells during replicative senescence

Background

Endothelial cell senescence is the state of permanent cell cycle arrest and plays a critical role in the pathogenesis of age-related diseases. However, a comprehensive understanding of the gene regulatory network, including genome-wide alternative splicing machinery, involved in endothelial cell senescence is lacking.

Results

We thoroughly described the transcriptome landscape of replicative senescent human umbilical vein endothelial cells. Genes with high connectivity showing a monotonic expression increase or decrease with the culture period were defined as hub genes in the co-expression network. Computational network analysis of these genes led to the identification of canonical and non-canonical senescence pathways, such as E2F and SIRT2 signaling, which were down-regulated in lipid metabolism, and chromosome organization processes pathways. Additionally, we showed that endothelial cell senescence involves alternative splicing. Importantly, the first and last exon types of splicing, as observed in FLT1 and ACACA, were preferentially altered among the alternatively spliced genes during endothelial senescence. We further identified novel microexons in PRUNE2 and PSAP, each containing 9 nt, which were altered within the specific domain during endothelial senescence.

Conclusions

These findings unveil the comprehensive transcriptome pathway and novel signaling regulated by RNA processing, including gene expression and splicing, in replicative endothelial senescence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08185-x.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

•

There is widespread transcript diversity in the cortex and many novel transcripts

•

Some genes display big differences in isoform number between human and mouse cortex

•

There is evidence of differential transcript usage between human fetal and adult cortex

•

There are many novel isoforms of genes associated with human brain disease

Computational analysis of sense-antisense chimeric transcripts reveals their potential regulatory features and the landscape of expression in human cells

Many human genes are transcribed from both strands and produce sense-antisense gene pairs. Sense-antisense (SAS) chimeric transcripts are produced upon the coalescing of exons/introns from both sense and antisense transcripts of the same gene. SAS chimera was first reported in prostate cancer cells. Subsequently, numerous SAS chimeras have been reported in the ChiTaRS-2.1 database. However, the landscape of their expression in human cells and functional aspects are still unknown. We found that longer palindromic sequences are a unique feature of SAS chimeras. Structural analysis indicates that a long hairpin-like structure formed by many consecutive Watson-Crick base pairs appears because of these long palindromic sequences, which possibly play a similar role as double-stranded RNA (dsRNA), interfering with gene expression. RNA-RNA interaction analysis suggested that SAS chimeras could significantly interact with their parental mRNAs, indicating their potential regulatory features. Here, 267 SAS chimeras were mapped in RNA-seq data from 16 healthy human tissues, revealing their expression in normal cells. Evolutionary analysis suggested the positive selection favoring sense-antisense fusions that significantly impacted the evolution of their function and structure. Overall, our study provides detailed insight into the expression landscape of SAS chimeras in human cells and identifies potential regulatory features.

Regulatory variation within 3’UTR of STAT5A correlates with sudden cardiac death in Chinese populations

Definitive diagnosis to sudden cardiac death (SCD) is often challenging since the postmortem examination on SCD victims could hardly demonstrate an adequate cause of death. It is therefore important to uncover the inherited risk component to SCD. Signal transducer and activators of transcription 5 A (STAT5A) is a member of the STAT family and a transcription factor that is activated by many cell ligands and associated with various cardiovascular processes. In this study, we performed a systematic variant screening on the STAT5A to filter potential functional genetic variations. Based on the screening results, an insertion/deletion polymorphism (rs3833144) in 3'UTR of STAT5A was selected as the candidate variant. A total of 159 SCD cases and 668 SCD matched healthy controls was enrolled to perform a case-control study and evaluate the association between rs3833144 and SCD susceptibility in Chinese populations. Logistic regression analysis showed that the deletion allele of rs3833144 had significantly increased the SCD risk (odds ratio (OR) = 1.54; 95% confidence interval (CI) = 1.18-2.01; P = 0.000955). Further genotype-expression eQTL analysis showed that samples with deletion allele appeared to lower expression of STAT5A, and in silico prediction suggested the local 3 D structure changes of STAT5A mRNA caused by the variant. On the other hand, the bioinformatic analysis presented that promoters of RARA and PTGES3L-AARSD1 could interact with rs3833144, and eQTL analysis showed the higher expression of both genes in samples with deletion allele. Dual-luciferase activity assays also suggested the significant regulatory role of rs3833144 in gene transcription. Our current data thus suggested a possible involvement of rs3833144 to SCD predisposition in Chinese populations and rs3833144 with potential function roles may become a candidate marker for SCD diagnosis and prevention.

Widespread premature transcription termination of Arabidopsis thaliana NLR genes by the spen protein FPA

Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants.

Interstitial Deletions Generating Fusion Genes

A fusion gene is the physical juxtaposition of two different genes resulting in a structure consisting of the head of one gene and the tail of the other. Gene fusion is often a primary neoplasia-inducing event in leukemias, lymphomas, solid malignancies as well as benign tumors. Knowledge about fusion genes is crucial not only for our understanding of tumorigenesis, but also for the diagnosis, prognostication, and treatment of cancer. Balanced chromosomal rearrangements, in particular translocations and inversions, are the most frequent genetic events leading to the generation of fusion genes. In the present review, we summarize the existing knowledge on chromosome deletions as a mechanism for fusion gene formation. Such deletions are mostly submicroscopic and, hence, not detected by cytogenetic analyses but by array comparative genome hybridization (aCGH) and/or high throughput sequencing (HTS). They are found across the genome in a variety of neoplasias. As tumors are increasingly analyzed using aCGH and HTS, it is likely that more interstitial deletions giving rise to fusion genes will be found, significantly impacting our understanding and treatment of cancer.

Detection of novel and recurrent conjoined genes in non-Hodgkin B-cell lymphoma

For this study, we investigated comprehensive expression of conjoined genes (CGs) in non-Hodgkin B-cell lymphoma (B-NHL) cell line KPUM-UH1 by using paired-end RNA sequencing. Furthermore, we analyzed the expression of these transcripts in an additional 21 cell lines, 37 primary samples of various malignancies and peripheral blood mononuclear cells of four normal individuals. Seventeen CGs were detected in KPUM-UH1: CTBS-GNG5, SRP9-EPHX1, RMND5A-ANAPC, OTX1-EHBP1, ATF2-CHN1, PRKAA1-TTC33, LARP1-MRPL22, LOC105379697-BAK1, TIAM2-SCAF8, SPAG1-VPS13B, WBP1L-CNNM2, NARS2-GAB2, CTSC-RAB38, VAMP1-CD27-AS1, LRRC37A2-NSF, UBA2-WTIP and ZNF600-ZNF611. To our knowledge, 10 of these genes have not been previously reported. The various characteristics of the CGs included in- and out-of-frame fusions, chimeras involving non-coding RNA and transcript variants. A finding of note was that LARP1-MRPL2 was characterized as in-frame fusion and was recurrently expressed in B-NHL samples. In this study, variety of CGs was expressed both in malignant and normal cells, some of which might be specific to lymphoma.

Aberrant splicing in neuroblastoma generates RNA-fusion transcripts and provides vulnerability to spliceosome inhibitors

The paucity of recurrent mutations has hampered efforts to understand and treat neuroblastoma. Alternative splicing and splicing-dependent RNA-fusions represent mechanisms able to increase the gene product repertoire but their role in neuroblastoma remains largely unexplored. Here we investigate the presence and possible roles of aberrant splicing and splicing-dependent RNA-fusion transcripts in neuroblastoma. In addition, we attend to establish whether the spliceosome can be targeted to treat neuroblastoma. Through analysis of RNA-sequenced neuroblastoma we show that elevated expression of splicing factors is a strong predictor of poor clinical outcome. Furthermore, we identified >900 primarily intrachromosomal fusions containing canonical splicing sites. Fusions included transcripts from well-known oncogenes, were enriched for proximal genes and in chromosomal regions commonly gained or lost in neuroblastoma. As a proof-of-principle that these fusions can generate altered gene products, we characterized a ZNF451-BAG2 fusion, producing a truncated BAG2-protein which inhibited retinoic acid induced differentiation. Spliceosome inhibition impeded neuroblastoma fusion expression, induced apoptosis and inhibited xenograft tumor growth. Our findings elucidate a splicing-dependent mechanism generating altered gene products in neuroblastoma and show that the spliceosome is a potential target for clinical intervention.

Read-through transcripts in lung: germline genetic regulation and correlation with the expression of other genes

Transcripts originating from the transcriptional read through of two adjacent, similarly oriented genes have been identified in normal and neoplastic tissues, but their functional role and the mechanisms that regulate their expression are mostly unknown. Here, we investigated whether the expression of read-through transcripts previously identified in the non-involved lung tissue of lung adenocarcinoma patients was genetically regulated. Data on genome-wide single nucleotide variant genotypes and expression levels of 10 read-through transcripts in 201 samples of lung tissue were combined to identify expression quantitative trait loci (eQTLs). Then, to identify genes whose expression levels correlated with the 10 read-through transcripts, we used whole transcriptome profiles available for 154 patients. For 8 read-though transcripts, we identified 60 eQTLs (false discovery rate <0.05), including 17 cis-eQTLs and 43 trans-eQTLs. These eQTLs did not maintain their behavior on the 'parental' genes involved in the read-through transcriptional event. The expression levels of 7 read-through transcripts were found to correlate with the expression of other genes: CHIA-PIFO and CTSC-RAB38 correlated with CHIA and RAB38, respectively, while 5 other read-through transcripts correlated with 43 unique non-parental transcripts; thus offering indications about the molecular processes in which these chimeric transcripts may be involved. We confirmed 9 eQTLs (for 4 transcripts) in the non-involved lung tissue from an independent series of 188 lung adenocarcinoma patients. Therefore, this study indicates that the expression of four read-through transcripts in normal lung tissue is under germline genetic regulation, and that this regulation is independent of that of the genes involved in the read-through event.

Compression of quantification uncertainty for scRNA-seq counts

Motivation

Quantification estimates of gene expression from single-cell RNA-seq (scRNA-seq) data have inherent uncertainty due to reads that map to multiple genes. Many existing scRNA-seq quantification pipelines ignore multi-mapping reads and therefore underestimate expected read counts for many genes. alevin accounts for multi-mapping reads and allows for the generation of ‘inferential replicates’, which reflect quantification uncertainty. Previous methods have shown improved performance when incorporating these replicates into statistical analyses, but storage and use of these replicates increases computation time and memory requirements.

Results

We demonstrate that storing only the mean and variance from a set of inferential replicates (‘compression’) is sufficient to capture gene-level quantification uncertainty, while reducing disk storage to as low as 9% of original storage, and memory usage when loading data to as low as 6%. Using these values, we generate ‘pseudo-inferential’ replicates from a negative binomial distribution and propose a general procedure for incorporating these replicates into a proposed statistical testing framework. When applying this procedure to trajectory-based differential expression analyses, we show false positives are reduced by more than a third for genes with high levels of quantification uncertainty. We additionally extend the Swish method to incorporate pseudo-inferential replicates and demonstrate improvements in computation time and memory usage without any loss in performance. Lastly, we show that discarding multi-mapping reads can result in significant underestimation of counts for functionally important genes in a real dataset.

Availability and implementation

makeInfReps and splitSwish are implemented in the R/Bioconductor fishpond package available at https://bioconductor.org/packages/fishpond. Analyses and simulated datasets can be found in the paper’s GitHub repo at https://github.com/skvanburen/scUncertaintyPaperCode.

Supplementary information

Supplementary data are available at Bioinformatics online.

Read-through circular RNAs reveal the plasticity of RNA processing mechanisms in human cells

In the human genome, there are several genes whose primary transcripts are both canonically and non-canonically spliced to generate mRNAs and RNA circles, respectively. These RNA circles are a novel class of long non-coding RNAs that became known as circular RNAs (circRNAs). Recently, a new type of circRNA was discovered and called read-through circRNAs (rt-circRNAs). They are hybrid circles that include coding exons from two adjacent and similarly oriented genes. The function of rt-circRNAs, as well as the impact of read-through transcription in our transcriptome, remains to be elucidated. Although we have just begun to scratch it, here I discussed some insights that these fascinating circRNAs are already giving us about the plasticity of RNA processing in our cells.

TIF-Seq2 disentangles overlapping isoforms in complex human transcriptomes

Eukaryotic transcriptomes are complex, involving thousands of overlapping transcripts. The interleaved nature of the transcriptomes limits our ability to identify regulatory regions, and in some cases can lead to misinterpretation of gene expression. To improve the understanding of the overlapping transcriptomes, we have developed an optimized method, TIF-Seq2, able to sequence simultaneously the 5' and 3' ends of individual RNA molecules at single-nucleotide resolution. We investigated the transcriptome of a well characterized human cell line (K562) and identified thousands of unannotated transcript isoforms. By focusing on transcripts which are challenging to be investigated with RNA-Seq, we accurately defined boundaries of lowly expressed unannotated and read-through transcripts putatively encoding fusion genes. We validated our results by targeted long-read sequencing and standard RNA-Seq for chronic myeloid leukaemia patient samples. Taking the advantage of TIF-Seq2, we explored transcription regulation among overlapping units and investigated their crosstalk. We show that most overlapping upstream transcripts use poly(A) sites within the first 2 kb of the downstream transcription units. Our work shows that, by paring the 5' and 3' end of each RNA, TIF-Seq2 can improve the annotation of complex genomes, facilitate accurate assignment of promoters to genes and easily identify transcriptionally fused genes.

DNA hypomethylation drives changes in MAGE-A gene expression resulting in alteration of proliferative status of cells

Melanoma Antigen Genes (MAGEs) are a family of genes that have piqued the interest of scientists for their unique expression pattern. A subset of MAGEs (Type I) are expressed in spermatogonial cells and in no other somatic tissue, and then re-expressed in many cancers. Type I MAGEs are often referred to as cancer-testis antigens due to this expression pattern, while Type II MAGEs are more ubiquitous in expression. This study determines the cause and consequence of the aberrant expression of the MAGE-A subfamily of cancer-testis antigens. We have discovered that MAGE-A genes are regulated by DNA methylation, as revealed by treatment with 5-azacytidine, an inhibitor of DNA methyltransferases. Furthermore, bioinformatics analysis of existing methylome sequencing data also corroborates our findings. The consequence of expressing certain MAGE-A genes is an increase in cell proliferation and colony formation and resistance to chemo-therapeutic agent 5-fluorouracil and DNA damaging agent sodium arsenite. Taken together, these data indicate that DNA methylation plays a crucial role in regulating the expression of MAGE-A genes which then act as drivers of cell proliferation, anchorage-independent growth and chemo-resistance that is critical for cancer-cell survival.

Chimeric RNAs in cancer

Chimeric RNAs are hybrid transcripts containing exons from two separate genes. Chimeric RNAs are traditionally considered to be transcribed from fusion genes caused by chromosomal rearrangement. These canonical chimeric RNAs are well characterized to be expressed in a cancer-unique pattern and/or act as oncogene products. However, benefited by the development of advanced deep sequencing technologies, novel types of non-canonical chimeric RNAs have been discovered to be generated from intergenic splicing without genomic aberrations. They can be formed through trans-splicing or cis-splicing between adjacent genes (cis-SAGe) mechanisms. Non-canonical chimeric RNAs are widely detected in normal physiology, although several have been shown to have a cancer-specific expression pattern. Further studies have indicated that some of them play fundamental roles in controlling cell growth and motility, and may have functions independent of the parental genes. These discoveries are unveiling a new layer of the functional transcriptome and are also raising the possibility of utilizing non-canonical chimeric RNAs as cancer diagnostic markers and therapeutic targets. In this chapter, we will overview different categories of chimeric RNAs and their expression in various types of cancerous and normal samples. Acknowledging that chimeric RNAs are not unique to cancer, we will discuss both bioinformatic and biological methods to identify credible cancer-specific chimeric RNAs. Furthermore, we will describe downstream methods to explore their molecular processing mechanisms and potential functions. A better understanding of the biogenesis mechanisms and functional products of cancer-specific chimeric RNAs will pave ways for the development of novel cancer biomarkers and therapeutic targets.

Gene Fusions Derived by Transcriptional Readthrough are Driven by Segmental Duplication in Human

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products.

Computational Detection of Known Pathogenic Gene Fusions in a Normal Tissue Database and Implications for Genetic Disease Research

Several recent studies have demonstrated the utility of RNA-Seq in the diagnosis of rare inherited disease. Diagnostic rates 35% higher than those previously achievable with DNA-Seq alone have been attained. These studies have primarily profiled gene expression and splicing defects, however, some have also shown that fusion transcripts are diagnostic or phenotypically relevant in patients with constitutional disorders. Fusion transcripts have traditionally been studied as oncogenic phenomena, with relevance only to cancer testing. Consequently, fusion detection algorithms were biased toward the detection of well-known oncogenic fusions, hindering their application to rare Mendelian genetic disease studies. A recent methodology published by the authors successfully tailored a traditional algorithm to the detection of pathogenic fusion events in inherited disease. A key mechanism of decreasing false positive or biologically benign events was comparison to a database of events detected in normal tissues. This approach is akin to population frequency-based filtering of genetic variants. It is predicated on the idea that pathogenic fusion transcripts are absent from normal tissue. We report on an analysis of RNA-Seq data from the genotype-tissue expression (GTEx) project in which known pathogenic fusions are computationally detected at low levels in normal tissues unassociated with the disease phenotype. Examples include archetypal cancer fusion transcripts, as well as fusions responsible for rare inherited disease. We consider potential explanations for the detectability of such transcripts and discuss the bearing such results have on the future profiling of genetic disease patients for pathogenic gene fusions.

Fusion transcripts in normal human cortex increase with age and show distinct genomic features for single cells and tissues

Fusion transcripts can contribute to diversity of molecular networks in the human cortex. In this study, we explored the occurrence of fusion transcripts in normal human cortex along with single neurons and astrocytes. We identified 1305 non-redundant fusion events from 388 transcriptomes representing 59 human cortices and 329 single cells. Our results indicate while the majority of fusion transcripts in human cortex are intra-chromosomal (85%), events found in single neurons and astrocytes were primarily inter-chromosomal (80%). The number of fusions in single neurons was significantly higher than that in single astrocytes (p < 0.05), indicating fusion as a possible contributor towards transcriptome diversity in neuronal cells. The identified fusions were largely private and 4 specific recurring events were found both in cortex and in single neurons but not in astrocytes. We found a significant increase in the number of fusion transcripts in human brain with increasing age both in single cells and whole cortex (p < 0.0005 and < 0.005, respectively). This is likely one of the many possible contributors for the inherent plasticity of the adult brain. The fusion transcripts in fetal brain were enriched for genes for long-term depression; while those in adult brain involved genes enriched for long-term potentiation pathways. Our findings demonstrate fusion transcripts are naturally occurring phenomenon spanning across the health-disease continuum, and likely contribute to the diverse molecular network of human brain.

ChiTaRS 5.0: the comprehensive database of chimeric transcripts matched with druggable fusions and 3D chromatin maps

Chimeric RNA transcripts are formed when exons from two genes fuse together, often due to chromosomal translocations, transcriptional errors or trans-splicing effect. While these chimeric RNAs produce functional proteins only in certain cases, they play a significant role in disease phenotyping and progression. ChiTaRS 5.0 (http://chitars.md.biu.ac.il/) is the latest and most comprehensive chimeric transcript repository, with 111 582 annotated entries from eight species, including 23 167 known human cancer breakpoints. The database includes unique information correlating chimeric breakpoints with 3D chromatin contact maps, generated from public datasets of chromosome conformation capture techniques (Hi-C). In this update, we have added curated information on druggable fusion targets matched with chimeric breakpoints, which are applicable to precision medicine in cancers. The introduction of a new section that lists chimeric RNAs in various cell-lines is another salient feature. Finally, using text-mining techniques, novel chimeras in Alzheimer's disease, schizophrenia, dyslexia and other diseases were collected in ChiTaRS. Thus, this improved version is an extensive catalogue of chimeras from multiple species. It extends our understanding of the evolution of chimeric transcripts in eukaryotes and contributes to the analysis of 3D genome conformational changes and the functional role of chimeras in the etiopathogenesis of cancers and other complex diseases.

Identification of Chimeric RNAs Using RNA-Seq Reads and Protein-Protein Interactions of Translated Chimeras

Chimeric RNA moieties typically consist of exons from two genes expressed from different genomic locations and produced by chromosomal translocations, trans-splicing or transcription errors. Recent advances in next-generation sequencing procedures have opened new horizons for identification of novel chimeric transcripts in various diseases in a personalized manner. Here we describe the detailed computational procedures to identify chimeric transcripts using RNA-seq reads. Moreover, we elaborate on the domain-domain co-occurrence method to detect alterations in chimeric protein-protein interaction (ChiPPI) networks produced by chimeric RNA that are translated to chimeric proteins.

Confirmation of Transcriptional Read-Through Events by RT-PCR

Chimeric RNAs can be formed by trans-splicing from different transcripts or cis-splicing of adjacent genes (cis-SAGe). Cis-SAGe results from read-through transcription of two neighbor genes. To investigate the mechanisms underlying intergenic splicing of adjacent genes, it is important to develop an assay to detect transcriptional read-through. Here, we describe a general RT-PCR based method to confirm the process for cis-SAGe candidates. In this method, we use PCR to amplify cDNA that is reverse transcribed from the read-through precursor mRNA. The result provides a foundation for further downstream mechanistic studies.

Fusion Transcripts of Adjacent Genes: New Insights into the World of Human Complex Transcripts in Cancer

The awareness of genome complexity brought a radical approach to the study of transcriptome, opening eyes to single RNAs generated from two or more adjacent genes according to the present consensus. This kind of transcript was thought to originate only from chromosomal rearrangements, but the discovery of readthrough transcription opens the doors to a new world of fusion RNAs. In the last years many possible intergenic cis-splicing mechanisms have been proposed, unveiling the origins of transcripts that contain some exons of both the upstream and downstream genes. In some cases, alternative mechanisms, such as trans-splicing and transcriptional slippage, have been proposed. Five databases, containing validated and predicted Fusion Transcripts of Adjacent Genes (FuTAGs), are available for the scientific community. A comparative analysis revealed that two of them contain the majority of the results. A complete analysis of the more widely characterized FuTAGs is provided in this review, including their expression pattern in normal tissues and in cancer. Gene structure, intergenic splicing patterns and exon junction sequences have been determined and here reported for well-characterized FuTAGs. The available functional data and the possible roles in cancer progression are discussed.

A tailored approach to fusion transcript identification increases diagnosis of rare inherited disease

BACKGROUND

RNA sequencing has been proposed as a means of increasing diagnostic rates in studies of undiagnosed rare inherited disease. Recent studies have reported diagnostic improvements in the range of 7.5-35% by profiling splicing, gene expression quantification and allele specific expression. To-date however, no study has systematically assessed the presence of gene-fusion transcripts in cases of germline disease. Fusion transcripts are routinely identified in cancer studies and are increasingly recognized as having diagnostic, prognostic or therapeutic relevance. Isolated reports exist of fusion transcripts being detected in cases of developmental and neurological phenotypes, and thus, systematic application of fusion detection to germline conditions may further increase diagnostic rates. However, current fusion detection methods are unsuited to the investigation of germline disease due to performance biases arising from their development using tumor, cell-line or in-silico data.

METHODS

We describe a tailored approach to fusion candidate identification and prioritization in a cohort of 47 undiagnosed, suspected inherited disease patients. We modify an existing fusion transcript detection algorithm by eliminating its cell line-derived filtering steps, and instead, prioritize candidates using a custom workflow that integrates genomic and transcriptomic sequence alignment, biological and technical annotations, customized categorization logic, and phenotypic prioritization.

RESULTS

We demonstrate that our approach to fusion transcript identification and prioritization detects genuine fusion events excluded by standard analyses and efficiently removes phenotypically unimportant candidates and false positive events, resulting in a reduced candidate list enriched for events with potential phenotypic relevance. We describe the successful genetic resolution of two previously undiagnosed disease cases through the detection of pathogenic fusion transcripts. Furthermore, we report the experimental validation of five additional cases of fusion transcripts with potential phenotypic relevance.

CONCLUSIONS

The approach we describe can be implemented to enable the detection of phenotypically relevant fusion transcripts in studies of rare inherited disease. Fusion transcript detection has the potential to increase diagnostic rates in rare inherited disease and should be included in RNA-based analytical pipelines aimed at genetic diagnosis.

Functional heritage: the evolution of chimeric RNA into a gene

Once believed to be unique features of neoplasia, chimeric RNAs are now being discovered in normal physiology. We speculated that some chimeric RNAs may be functional precursors of genes, and that forming chimeric RNA at the transcriptional level may be a 'trial' mechanism before the functional element is fixed into the genome. Supporting this idea, we identified a chimeric RNA, HNRNPA1L2-SUGT1 (H-S), whose sequence is highly similar to that of a 'pseudogene' MRPS31P5. Sequence analysis revealed that MRPS31P5 transcript is more similar to H-S chimeric RNA than its 'parent' gene, MRPS31. Evolutionarily, H-S precedes MRPS31P5, as it can be detected bioinformatically and experimentally in marmosets, which do not yet possess MRPS31P5 in their genome. Conversely, H-S is minimally expressed in humans, while instead, MRPS31P5 is abundantly expressed. Silencing H-S in marmoset cells resulted in similar phenotype as silencing MRPS31P5 in human cells. In addition, whole transcriptome analysis and candidate downstream target validation revealed common signalling pathways shared by the two transcripts. Interestingly, H-S failed to rescue the phenotype caused by silencing MPRS31P5 in human and rhesus cells, whereas MRPS31P5 can at least partially rescue the phenotype caused by silencing H-S in marmoset cells, suggesting that MRPS31P5 may have further evolved into a distinct entity. Thus, multiple lines of evidence support that MRPS31P5 is not truly a pseudogene of MRPS31, but a likely functional descendent of H-S chimera. Instead being a gene fusion product, H-S is a product of cis-splicing between adjacent genes, while MRPS31P5 is likely produced by genome rearrangement.

A cell-based splicing reporter system to identify regulators of cis-splicing between adjacent genes

Chimeric RNAs generated by cis-splicing between adjacent genes (cis-SAGe) are increasingly recognized as a widespread phenomenon. These chimeric messenger RNAs are present in normal human cells, and are also detected in various cancers. The mechanisms for how this group of chimeras is formed are not yet clear, in part due to the lack of a tractable system for their experimental investigation. Here we developed a fast, easy and versatile cell-based reporter system to identify regulators of cis-SAGe. The reporter, consisting of four main cassettes, simultaneously measures the effects of a candidate regulator on cis-SAGe and canonical splicing. Using this cell-based assay, we screened 102 candidate factors involved in RNA pol II cleavage and termination, elongation, splicing, alternative splicing and R-loop formation. We discovered that two factors, SRRM1 and SF3B1, affect not only cis-SAGe chimeras, but also other types of chimeric RNAs in a genome-wide fashion. This system can be used for studying trans-acting factors and cis-acting sequence elements and factors, as well as for screening small molecule inhibitors.

Alternative Transcription at Venom Genes and Its Role as a Complementary Mechanism for the Generation of Venom Complexity in the Common House Spider

The complex composition of venom, a proteinaceous secretion used by diverse animal groups for predation or defense, is typically viewed as being driven by gene duplication in conjunction with positive selection, leading to large families of diversified toxins with selective venom gland expression. Yet, the production of alternative transcripts at venom genes is often overlooked as another potentially important process that could contribute proteins to venom, and requires comprehensive datasets integrating genome and transcriptome sequences together with proteomic characterization of venom to be fully documented. In the common house spider, Parasteatoda tepidariorum, we used RNA sequencing of four tissue types in conjunction with the sequenced genome to provide a comprehensive transcriptome annotation. We also used mass spectrometry to identify a minimum of 99 distinct proteins in P tepidariorum venom, including at least 33 latrotoxins, pore-forming neurotoxins shared with the confamilial black widow. We found that venom proteins are much more likely to come from multiple transcript genes, whose transcripts produced distinct protein sequences. The presence of multiple distinct proteins in venom from transcripts at individual genes was confirmed for eight loci by mass spectrometry, and is possible at 21 others. Alternative transcripts from the same gene, whether encoding or not encoding a protein found in venom, showed a range of expression patterns, but were not necessarily restricted to the venom gland. However, approximately half of venom protein encoding transcripts were found among the 1,318 transcripts with strongly venom gland biased expression. Our findings revealed an important role for alternative transcription in generating venom protein complexity and expanded the traditional model of venom evolution.

CINSARC signature as a prognostic marker for clinical outcome in sarcomas and beyond

Prognostication is a key issue for sarcoma patients' care as it triggers the therapeutic approach including chemotherapy, which is still not standard for localized patients. Current prognostic evaluation, based on the FNCLCC grading system, has recently been improved by the CINSARC signature outperforming histology-based grading system by identifying high-risk patients in every grade, even in those considered as low. CINSARC is an expression-based signature related to mitosis and chromosome integrity with prognostic value in a wide range of cancers additional to sarcoma. First developed with frozen material, CINSARC is now coupled with NanoString technology allowing evaluation from FFPE blocks used in clinical practice. Consequently, CINSARC is currently evaluated in clinical trials with a dual objective of demonstrating the benefit of chemotherapy in sarcoma patients and testing its response prediction. Considering its overarching value in oncology, its development is welcome in any cancers where the prognostication needs to be improved.

The Landscape of Circular RNA in Cancer

Circular RNAs (circRNAs) are an intriguing class of RNA due to their covalently closed structure, high stability, and implicated roles in gene regulation. Here, we used an exome capture RNA sequencing protocol to detect and characterize circRNAs across >2,000 cancer samples. When compared against Ribo-Zero and RNase R, capture sequencing significantly enhanced the enrichment of circRNAs and preserved accurate circular-to-linear ratios. Using capture sequencing, we built the most comprehensive catalog of circRNA species to date: MiOncoCirc, the first database to be composed primarily of circRNAs directly detected in tumor tissues. Using MiOncoCirc, we identified candidate circRNAs to serve as biomarkers for prostate cancer and were able to detect circRNAs in urine. We further detected a novel class of circular transcripts, termed read-through circRNAs, that involved exons originating from different genes. MiOncoCirc will serve as a valuable resource for the development of circRNAs as diagnostic or therapeutic targets across cancer types.

Transcriptional activation of a chimeric retrogene PIPSL in a hominoid ancestor

Retrogenes are a class of functional genes derived from the mRNA of various intron-containing genes. PIPSL was created through a unique mechanism, whereby distinct genes were assembled at the RNA level, and the resulting chimera was then reverse transcribed and integrated into the genome by the L1 retrotransposon. Expression of PIPSL RNA via its transcription start sites (TSSs) has been confirmed in the testes of humans and chimpanzee. Here, we demonstrated that PIPSL RNA is expressed in the testis of the white-handed gibbon. The 5'-end positions of gibbon RNAs were confined to a narrow range upstream of the PIPSL start codon and overlapped with those of orangutan and human, suggesting that PIPSL TSSs are similar among hominoid species. Reporter assays using a luciferase gene and the flanking sequences of human PIPSL showed that an upstream sequence exhibits weak promoter activity in human cells. Our findings suggest that PIPSL might have acquired a promoter at an early stage of hominoid evolution before the divergence of gibbons and ultimately retained similar TSSs in all of the lineages. Moreover, the upstream sequence derived from the phosphatidylinositol-4-phosphate 5-kinase, type I, alpha 5' untranslated region and/or neighboring repetitive sequences in the genome possibly exhibits promoter activity. Furthermore, we observed that a TATA-box-like sequence has emerged by nucleotide substitution in a lineage leading to humans, with this possibly responsible for a broader distribution of the human PIPSL TSSs.