Tandem chimerism as a means to increase protein complexity in the human genome

The landscape of fusion transcripts in plants: a new insight into genome complexity

BACKGROUND

Fusion transcripts (FTs), generated by the fusion of genes at the DNA level or RNA-level splicing events significantly contribute to transcriptome diversity. FTs are usually considered unique features of neoplasia and serve as biomarkers and therapeutic targets for multiple cancers. The latest findings show the presence of FTs in normal human physiology. Several discrete reports mentioned the presence of fusion transcripts in planta, has important roles in stress responses, morphological alterations, or traits (e.g. seed size, etc.).

RESULTS

In this study, we identified 169,197 fusion transcripts in 2795 transcriptome datasets of Arabidopsis thaliana, Cicer arietinum, and Oryza sativa by using a combination of tools, and confirmed the translational activity of 150 fusion transcripts through proteomic datasets. Analysis of the FT junction sequences and their association with epigenetic factors, as revealed by ChIP-Seq datasets, demonstrated an organised process of fusion formation at the DNA level. We investigated the possible impact of three-dimensional chromatin conformation on intra-chromosomal fusion events by leveraging the Hi-C datasets with the incidence of fusion transcripts. We further utilised the long-read RNA-Seq datasets to validate the most reoccurring fusion transcripts in each plant species followed by further authentication through RT-PCR and Sanger sequencing.

CONCLUSIONS

Our findings suggest that a significant portion of fusion events may be attributed to alternative splicing during transcription, accounting for numerous fusion events without a proportional increase in the number of RNA pairs. Even non-nuclear DNA transcripts from mitochondria and chloroplasts can participate in intra- and inter-chromosomal fusion formation. Genes in close spatial proximity are more prone to undergoing fusion formation, especially in intra-chromosomal FTs. Most of the fusion transcripts may not undergo translation and serve as long non-coding RNAs. The low validation rate of FTs in plants indicated that the fusion transcripts are expressed at very low levels, like in the case of humans. FTs often originate from parental genes involved in essential biological processes, suggesting their relevance across diverse tissues and stress conditions. This study presents a comprehensive repository of fusion transcripts, offering valuable insights into their roles in vital physiological processes and stress responses.

Full-length RNA-Seq of the RHOH gene in human B cells reveals new exons and splicing patterns

The RhoH protein is a member of the Ras superfamily of guanosine triphosphate-binding proteins. RhoH is an atypical Rho family member that is always GTP-bound and thus always activated. It is restrictively expressed in normal hematopoietic cells, where it is a negative regulator of cell growth and survival. We previously analyzed the RHOH gene structure and demonstrated that this gene is composed of 7 exons, one single encoding exon located at the 3' extremity of the gene, preceded by 6 noncoding exons. To further understand the transcription events associated with this gene, we performed full-length RNA-Seq on 12 B-cell lines. We identified new exons, new splice events and new splice sites, leading to the discovery of 38 RHOH mRNA molecules, 27 of which have never been described before. Here, we also describe new fusion transcripts. Moreover, our method allowed quantitative measurements of the different mRNA species relative to each other in relation to B-cell differentiation.

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.

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.

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.

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.

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.

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.

Novel Gene Fusions in Glioblastoma Tumor Tissue and Matched Patient Plasma

Sequencing studies have provided novel insights into the heterogeneous molecular landscape of glioblastoma (GBM), unveiling a subset of patients with gene fusions. Tissue biopsy is highly invasive, limited by sampling frequency and incompletely representative of intra-tumor heterogeneity. Extracellular vesicle-based liquid biopsy provides a minimally invasive alternative to diagnose and monitor tumor-specific molecular aberrations in patient biofluids. Here, we used targeted RNA sequencing to screen GBM tissue and the matched plasma of patients (n = 9) for RNA fusion transcripts. We identified two novel fusion transcripts in GBM tissue and five novel fusions in the matched plasma of GBM patients. The fusion transcripts FGFR3-TACC3 and VTI1A-TCF7L2 were detected in both tissue and matched plasma. A longitudinal follow-up of a GBM patient with a FGFR3-TACC3 positive glioma revealed the potential of monitoring RNA fusions in plasma. In summary, we report a sensitive RNA-seq-based liquid biopsy strategy to detect RNA level fusion status in the plasma of GBM patients.

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.

Pan-Cancer Analysis Reveals the Diverse Landscape of Novel Sense and Antisense Fusion Transcripts

Gene fusions that contribute to oncogenicity can be explored for identifying cancer biomarkers and potential drug targets. To investigate the nature and distribution of fusion transcripts in cancer, we examined the transcriptome data of about 9,000 primary tumors from 33 different cancers in TCGA (The Cancer Genome Atlas) along with cell line data from CCLE (Cancer Cell Line Encyclopedia) using ChimeRScope, a novel fusion detection algorithm. We identified several fusions with sense (canonical, 39%) or antisense (non-canonical, 61%) transcripts recurrent across cancers. The majority of the recurrent non-canonical fusions found in our study are novel, unexplored, and exhibited highly variable profiles across cancers, with breast cancer and glioblastoma having the highest and lowest rates, respectively. Overall, 4,344 recurrent fusions were identified from TCGA in this study, of which 70% were novel. Additional analysis of 802 tumor-derived cell line transcriptome data across 20 cancers revealed significant variability in recurrent fusion profiles between primary tumors and corresponding cell lines. A subset of canonical and non-canonical fusions was validated by examining the structural variation evidence in whole-genome sequencing (WGS) data or by Sanger sequencing of fusion junctions. Several recurrent fusion genes identified in our study show promise for drug repurposing in basket trials and present opportunities for mechanistic studies.

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.

The UGTome: The expanding diversity of UDP glycosyltransferases and its impact on small molecule metabolism

The UDP glycosyltransferase (UGT) superfamily of enzymes is responsible for the metabolism and clearance of thousands of lipophilic chemicals including drugs, toxins and endogenous signaling molecules. They provide a protective interface between the organism and its chemical-rich environment, as well as controlling critical signaling pathways to maintain healthy tissue function. UGTs are associated with drug responses and interactions, as well as a wide range of diseases including cancer. The human genome contains 22 UGT genes; however as befitting their exceptionally diverse substrate ranges and biological activities, the output of these UGT genes is functionally diversified by multiple processes including alternative splicing, post-translational modification, homo- and hetero-oligomerization, and interactions with other proteins. All UGT genes are subject to extensive alternative splicing generating variant/truncated UGT proteins with altered functions including the capacity to dominantly modulate/inhibit cognate full-length forms. Heterotypic oligomerization of different UGTs can alter kinetic properties relative to monotypic complexes, and potentially produce novel substrate specificities. Moreover, the recently profiled interactions of UGTs with non-UGT proteins may facilitate coordination between different metabolic processes, as well as providing opportunities for UGTs to engage in novel 'moonlighting' functions. Herein we provide a detailed and comprehensive review of all known modes of UGT functional diversification and propose a UGTome model to describe the resulting expansion of metabolic capacity and its potential to modulate drug/xenobiotic responses and cell behaviours in normal and disease contexts.

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.

Identification of candidate neoantigens produced by fusion transcripts in human osteosarcomas

Osteosarcomas are characterized by highly disrupted genomes. Although osteosarcomas lack common fusions, we find evidence of many tumour specific gene-gene fusion transcripts, likely due to chromosomal rearrangements and expression of transcription-induced chimeras. Most of the fusions result in out-of-frame transcripts, potentially capable of producing long novel protein sequences and a plethora of neoantigens. To identify fusions, we explored RNA-sequencing data to obtain detailed knowledge of transcribed fusions, by creating a novel program to compare fusions identified by deFuse to de novo transcripts generated by Trinity. This allowed us to confirm the deFuse results and identify unusual splicing patterns associated with fusion events. Using various existing tools combined with this custom program, we developed a pipeline for the identification of fusion transcripts applicable as targets for immunotherapy. In addition to identifying candidate neoantigens associated with fusions, we were able to use the pipeline to establish a method for measuring the frequency of fusion events, which correlated to patient outcome, as well as highlight some similarities between canine and human osteosarcomas. The results of this study of osteosarcomas underscores the numerous benefits associated with conducting a thorough analysis of fusion events within cancer samples.

Intergenic Splicing between Four Adjacent UGT Genes (2B15, 2B29P2, 2B17, 2B29P1) Gives Rise to Variant UGT Proteins That Inhibit Glucuronidation via Protein-Protein Interactions

Recent studies have investigated alternative splicing profiles of UDP-glucuronosyltransferase (UGT) genes and identified over 130 different alternatively spliced UGT transcripts. Although UGT genes are highly clustered, the formation of chimeric transcripts by intergenic splicing between two or more UGT genes has not yet been reported. This study identified 12 chimeric transcripts (chimeras A-L) containing exons from two or three genes of the four neighboring UGT genes (UGT2B15, UGT2B29P2, UGT2B17, and UGT2B29P1) in human liver and prostate cancer cells. These chimeras typically contain the first five exons of UGT2B15 or UGT2B17 (exons 1-5) spliced to a terminal exon (exon 6) from a downstream UGT gene. Hence they encode truncated UGTs with novel C-terminal peptides. Functional assays of representative chimeric UGT proteins (termed chimeric UGT2B15 and chimeric UGT2B17) showed that they are inactive and can repress the activity of wild-type UGTs. Coimmunoprecipitation assays demonstrated heterotypic interactions between chimeric UGT2B15 (or chimeric UGT2B17) and the UGT2B7 protein. Thus oligomerization of the chimeric UGTs with wild-type UGTs may explain their inhibitory activity. Studies in breast and prostate cancer cells showed that both wild-type and chimeric UGT2B15 and UGT2B17 transcripts are regulated in a similar way at the transcriptional level by sex hormones through their canonical promoters but are differentially regulated at the post-transcriptional level by micro-RNA 376c via their unique 3'-untranslated regions. In conclusion, the formation of chimeric transcripts by intergenic splicing among UGT genes represents a novel mechanism contributing to the diversity of the human UGT transcriptome and proteome. The differential post-transcriptional regulation of wild-type and variant transcripts by micro-RNAs may contribute to their deregulated expression in cancer.

Identification of novel enriched recurrent chimeric COL7A1-UCN2 in human laryngeal cancer samples using deep sequencing

BACKGROUND

As hybrid RNAs, transcription-induced chimeras (TICs) may have tumor-promoting properties, and some specific chimeras have become important diagnostic markers and therapeutic targets for cancer.

METHODS

We examined 23 paired laryngeal cancer (LC) tissues and adjacent normal mucous membrane tissue samples (ANMMTs). Three of these pairs were used for comparative transcriptomic analysis using high-throughput sequencing. Furthermore, we used real-time polymerase chain reaction (RT-PCR) for further validation in 20 samples. The Kaplan-Meier method and Cox regression model were used for the survival analysis.

RESULTS

We identified 87 tumor-related TICs and found that COL7A1-UCN2 had the highest frequency in LC tissues (13/23; 56.5%), whereas none of the ANMMTs were positive (0/23; p < 0.0001). COL7A1-UCN2, generated via alternative splicing in LC tissue cancer cells, had disrupted coding regions, but it down-regulated the mRNA expression of COL7A1 and UCN2. Both COL7A1 and UCN2 were down-expressed in LC tissues as compared to their paired ANMMTs. The COL7A1:β-actin ratio in COL7A1-UCN2-positive LC samples was significantly lower than that in COL7A1-UCN2-negative samples (p = 0.019). Likewise, the UCN2:β-actin ratio was also decreased (p = 0.21). Furthermore, COL7A1-UCN2 positivity was significantly associated with the overall survival of LC patients (p = 0.032; HR, 13.2 [95%CI, 1.2-149.5]).

CONCLUSION

LC cells were enriched in the recurrent chimera COL7A1-UCN2, which potentially affected cancer stem cell transition, promoted epithelial-mesenchymal transition in LC, and resulted in poorer prognoses.

Novel Intergenically Spliced Chimera, NFATC3-PLA2G15, Is Associated with Aggressive T-ALL Biology and Outcome

Leukemias are frequently characterized by the expression of oncogenic fusion chimeras that normally arise due to chromosomal rearrangements. Intergenically spliced chimeric RNAs (ISC) are transcribed in the absence of structural genomic changes, and aberrant ISC expression is now recognized as a potential driver of cancer. To better understand these potential oncogenic drivers, high-throughput RNA sequencing was performed on T-acute lymphoblastic leukemia (T-ALL) patient specimens (n = 24), and candidate T-ALL-related ISCs were identified (n = 55; a median of 4/patient). In-depth characterization of the NFATC3-PLA2G15 chimera, which was variably expressed in primary T-ALL, was performed. Functional assessment revealed that the fusion had lower activity than wild-type NFATC3 in vitro, and T-ALLs with elevated NFATC3-PLA2G15 levels had reduced transcription of canonical NFAT pathway genes in vivo Strikingly, high expression of the NFATC3-PLA2G15 chimera correlated with aggressive disease biology in murine patient-derived T-ALL xenografts, and poor prognosis in human T-ALL patients. Mol Cancer Res; 16(3); 470-5. ©2018 AACR.

The human GCOM1 complex gene interacts with the NMDA receptor and internexin-alpha

The known functions of the human GCOM1 complex hub gene include transcription elongation and the intercalated disk of cardiac myocytes. However, in all likelihood, the gene's most interesting, and thus far least understood, roles will be found in the central nervous system. To investigate the functions of the GCOM1 gene in the CNS, we have cloned human and rat brain cDNAs encoding novel, 105 kDa GCOM1 combined (Gcom) proteins, designated Gcom15, and identified a new group of GCOM1 interacting genes, termed Gints, from yeast two-hybrid (Y2H) screens. We showed that Gcom15 interacts with the NR1 subunit of the NMDA receptor by co-expression in heterologous cells, in which we observed bi-directional co-immunoprecipitation of human Gcom15 and murine NR1. Our Y2H screens revealed 27 novel GCOM1 interacting genes, many of which are synaptic proteins and/or play roles in neurologic diseases. Finally, we showed, using rat brain protein preparations, that the Gint internexin-alpha (INA), a known interactor of the NMDAR, co-IPs with GCOM1 proteins, suggesting a GCOM1-GRIN1-INA interaction and a novel pathway that may be relevant to neuroprotection.

Read-through transcripts in normal human lung parenchyma are down-regulated in lung adenocarcinoma

Read-through transcripts result from the continuous transcription of adjacent, similarly oriented genes, with the splicing out of the intergenic region. They have been found in several neoplastic and normal tissues, but their pathophysiological significance is unclear. We used high-throughput sequencing of cDNA fragments (RNA-Seq) to identify read-through transcripts in the non-involved lung tissue of 64 surgically treated lung adenocarcinoma patients. A total of 52 distinct read-through species was identified, with 24 patients having at least one read-through event, up to a maximum of 17 such transcripts in one patient. Sanger sequencing validated 28 of these transcripts and identified an additional 15, for a total of 43 distinct read-through events involving 35 gene pairs. Expression levels of 10 validated read-through transcripts were measured by quantitative PCR in pairs of matched non-involved lung tissue and lung adenocarcinoma tissue from 45 patients. Higher expression levels were observed in normal lung tissue than in the tumor counterpart, with median relative quantification ratios between normal and tumor varying from 1.90 to 7.78; the difference was statistically significant (P < 0.001, Wilcoxon's signed-rank test for paired samples) for eight transcripts: ELAVL1–TIMM44, FAM162B–ZUFSP, IFNAR2–IL10RB, INMT–FAM188B, KIAA1841–C2orf74, NFATC3–PLA2G15, SIRPB1–SIRPD, and SHANK3–ACR. This report documents the presence of read-through transcripts in apparently normal lung tissue, with inter-individual differences in patterns and abundance. It also shows their down-regulation in tumors, suggesting that these chimeric transcripts may function as tumor suppressors in lung tissue.

Connections between Transcription Downstream of Genes and cis-SAGe Chimeric RNA

cis-Splicing between adjacent genes (cis-SAGe) is being recognized as one way to produce chimeric fusion RNAs. However, its detail mechanism is not clear. Recent study revealed induction of transcriptions downstream of genes (DoGs) under osmotic stress. Here, we investigated the influence of osmotic stress on cis-SAGe chimeric RNAs and their connection to DoGs. We found, the absence of induction of at least some cis-SAGe fusions and/or their corresponding DoGs at early time point(s). In fact, these DoGs and their cis-SAGe fusions are inversely correlated. This negative correlation was changed to positive at a later time point. These results suggest a direct competition between the two categories of transcripts when total pool of readthrough transcripts is limited at an early time point. At a later time point, DoGs and corresponding cis-SAGe fusions are both induced, indicating that total readthrough transcripts become more abundant. Finally, we observed overall enhancement of cis-SAGe chimeric RNAs in KCl-treated samples by RNA-Seq analysis.

The Output of Protein-Coding Genes Shifts to Circular RNAs When the Pre-mRNA Processing Machinery Is Limiting

Many eukaryotic genes generate linear mRNAs and circular RNAs, but it is largely unknown how the ratio of linear to circular RNA is controlled or modulated. Using RNAi screening in Drosophila cells, we identify many core spliceosome and transcription termination factors that control the RNA outputs of reporter and endogenous genes. When spliceosome components were depleted or inhibited pharmacologically, the steady-state levels of circular RNAs increased while expression of their associated linear mRNAs concomitantly decreased. Upon inhibiting RNA polymerase II termination via depletion of the cleavage/polyadenylation machinery, circular RNA levels were similarly increased. This is because readthrough transcripts now extend into downstream genes and are subjected to backsplicing. In total, these results demonstrate that inhibition or slowing of canonical pre-mRNA processing events shifts the steady-state output of protein-coding genes toward circular RNAs. This is in part because nascent RNAs become directed into alternative pathways that lead to circular RNA production.

Chimeric RNAs and their implications in cancer

Chimeric RNAs have been believed to be solely produced by gene fusions resulting from chromosomal rearrangement, thus unique features of cancer. Detected chimeric RNAs have also been viewed as surrogates for the presence of gene fusions. However, more and more research has demonstrated that chimeric RNAs in general are not a hallmark of cancer, but rather widely present in non-cancerous cells and tissues. At the same time, they may be produced by other mechanisms other than chromosomal rearrangement. The field of non-canonical chimeric RNAs is still in its infancy, with many challenges ahead, including the lack of a unified terminology. However, we believe that these non-canonical chimeric RNAs will have significant impacts in cancer detection and treatment.

RNA-seq Based Transcription Characterization of Fusion Breakpoints as a Potential Estimator for Its Oncogenic Potential

Based on high-throughput sequencing technology, the detection of gene fusions is no longer a big challenge but estimating the oncogenic potential of fusion genes remains challenging. Recent studies successfully applied machine learning methods and gene structural and functional features of fusion mutation to predict their oncogenic potentials. However, the transcription characterizations features of fusion genes have not yet been studied. In this study, based on the clonal evolution theory, we hypothesized that a fusion gene is more likely to be an oncogenic genomic alteration, if the neoplastic cells harboring this fusion mutation have larger clonal size than other neoplastic cells in a tumor. We proposed a novel method, called iFCR (internal Fusion Clone Ratio), given an estimation of oncogenic potential for fusion mutations. We have evaluated the iFCR method in three public cancer transcriptome sequencing datasets; the results demonstrated that the fusion mutations occurring in tumor samples have higher internal fusion clone ratio than normal samples. And the most frequent prostate cancer fusion mutation, TMPRSS2-ERG, appears to have a remarkably higher iFCR value in all three independent patients. The preliminary results suggest that the internal fusion clone ratio might potentially advantage current fusion mutation oncogenic potential prediction methods.

De novo assembly and characterization of breast cancer transcriptomes identifies large numbers of novel fusion-gene transcripts of potential functional significance

Background

Gene-fusion or chimeric transcripts have been implicated in the onset and progression of a variety of cancers. Massively parallel RNA sequencing (RNA-Seq) of the cellular transcriptome is a promising approach for the identification of chimeric transcripts of potential functional significance. We report here the development and use of an integrated computational pipeline for the de novo assembly and characterization of chimeric transcripts in 55 primary breast cancer and normal tissue samples.

Methods

An integrated computational pipeline was employed to screen the transcriptome of breast cancer and control tissues for high-quality RNA-sequencing reads. Reads were de novo assembled into contigs followed by reference genome mapping. Chimeric transcripts were detected, filtered and characterized using our R-SAP algorithm. The relative abundance of reads was used to estimate levels of gene expression.

Results

De novo assembly allowed for the accurate detection of 1959 chimeric transcripts to nucleotide level resolution and facilitated detailed molecular characterization and quantitative analysis. A number of the chimeric transcripts are of potential functional significance including 79 novel fusion-protein transcripts and many chimeric transcripts with alterations in their un-translated leader regions. A number of chimeric transcripts in the cancer samples mapped to genomic regions devoid of any known genes. Several ‘pro-neoplastic’ fusions comprised of genes previously implicated in cancer are expressed at low levels in normal tissues but at high levels in cancer tissues.

Conclusions

Collectively, our results underscore the utility of deep sequencing technologies and improved bioinformatics workflows to uncover novel and potentially significant chimeric transcripts in cancer and normal somatic tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0289-7) contains supplementary material, which is available to authorized users.

A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing

BACKGROUND

Despite the economic importance of sugarcane in sugar and bioenergy production, there is not yet a reference genome available. Most of the sugarcane transcriptomic studies have been based on Saccharum officinarum gene indices (SoGI), expressed sequence tags (ESTs) and de novo assembled transcript contigs from short-reads; hence knowledge of the sugarcane transcriptome is limited in relation to transcript length and number of transcript isoforms.

RESULTS

The sugarcane transcriptome was sequenced using PacBio isoform sequencing (Iso-Seq) of a pooled RNA sample derived from leaf, internode and root tissues, of different developmental stages, from 22 varieties, to explore the potential for capturing full-length transcript isoforms. A total of 107,598 unique transcript isoforms were obtained, representing about 71% of the total number of predicted sugarcane genes. The majority of this dataset (92%) matched the plant protein database, while just over 2% was novel transcripts, and over 2% was putative long non-coding RNAs. About 56% and 23% of total sequences were annotated against the gene ontology and KEGG pathway databases, respectively. Comparison with de novo contigs from Illumina RNA-Sequencing (RNA-Seq) of the internode samples from the same experiment and public databases showed that the Iso-Seq method recovered more full-length transcript isoforms, had a higher N50 and average length of largest 1,000 proteins; whereas a greater representation of the gene content and RNA diversity was captured in RNA-Seq. Only 62% of PacBio transcript isoforms matched 67% of de novo contigs, while the non-matched proportions were attributed to the inclusion of leaf/root tissues and the normalization in PacBio, and the representation of more gene content and RNA classes in the de novo assembly, respectively. About 69% of PacBio transcript isoforms and 41% of de novo contigs aligned with the sorghum genome, indicating the high conservation of orthologs in the genic regions of the two genomes.

CONCLUSIONS

The transcriptome dataset should contribute to improved sugarcane gene models and sugarcane protein predictions; and will serve as a reference database for analysis of transcript expression in sugarcane.

Genome and transcriptome profiling of fibrolamellar hepatocellular carcinoma demonstrates p53 and IGF2BP1 dysregulation

Fibrolamellar hepatocellular carcinoma (FL-HCC) is a rare variant of HCC that most frequently affects young adults. Because of its rarity and an absence of preclinical models, our understanding of FL-HCC is limited. Our objective was to analyze chromosomal alterations and dysregulated gene expression in tumor specimens collected at a single center during two decades of experience with FL-HCC. We analyzed 38 specimens from 26 patients by array comparative genomic hybridiziation (aCGH) and 35 specimens from 15 patients by transcriptome sequencing (RNA-seq). All tumor specimens exhibited genomic instability, with a higher frequency of genomic amplifications or deletions in metastatic tumors. The regions encoding 71 microRNAs (miRs) were deleted in at least 25% of tumor specimens. Five of these recurrently deleted miRs targeted the insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) gene product, and a correlating 100-fold upregulation of IGF2BP1 mRNA was seen in tumor specimens. Transcriptome analysis demonstrated intrapatient tumor similarity, independent of recurrence site or time. The p53 tumor suppressor pathway was downregulated as demonstrated by both aCGH and RNA-seq analysis. Notch, EGFR, NRAS, and RB1 pathways were also significantly dysregulated in tumors compared with normal liver tissue. The findings illuminate the genomic and transcriptomic landscape of this rare disease and provide insight into dysregulated oncogenic pathways and potential therapeutic targets in FL-HCC.

The Evolving Definition of the Term "Gene"

This paper presents a history of the changing meanings of the term "gene," over more than a century, and a discussion of why this word, so crucial to genetics, needs redefinition today. In this account, the first two phases of 20th century genetics are designated the "classical" and the "neoclassical" periods, and the current molecular-genetic era the "modern period." While the first two stages generated increasing clarity about the nature of the gene, the present period features complexity and confusion. Initially, the term "gene" was coined to denote an abstract "unit of inheritance," to which no specific material attributes were assigned. As the classical and neoclassical periods unfolded, the term became more concrete, first as a dimensionless point on a chromosome, then as a linear segment within a chromosome, and finally as a linear segment in the DNA molecule that encodes a polypeptide chain. This last definition, from the early 1960s, remains the one employed today, but developments since the 1970s have undermined its generality. Indeed, they raise questions about both the utility of the concept of a basic "unit of inheritance" and the long implicit belief that genes are autonomous agents. Here, we review findings that have made the classic molecular definition obsolete and propose a new one based on contemporary knowledge.

A Comprehensive Analysis of Cell Type-Specific Nuclear RNA From Neurons and Glia of the Brain

BACKGROUND

Studies in psychiatric genetics have identified >100 loci associated with disease risk, yet many of these loci are distant from protein coding genes. Recent characterization of the transcriptional landscape of cell lines and whole tissues has suggested widespread transcription in both coding and noncoding regions of the genome, including differential expression from loci that produce regulatory noncoding RNAs that function within the nucleus; however, the nuclear transcriptome of specific cell types in the brain has not been previously investigated.

METHODS

We defined the nuclear transcriptional landscape of the three major cellular divisions of the nervous system using flow sorting of genetically labeled nuclei from bacTRAP mouse lines. Next, we characterized the unique expression of coding, noncoding, and intergenic RNAs in the mature mouse brain with RNA-Seq and validation with independent methods.

RESULTS

We found diverse expression across the cell types of all classes of RNAs, including long noncoding RNAs, several of which were confirmed as highly enriched in the nuclei of specific cell types using anatomic methods. We also discovered several examples of cell type-specific expression of tandem gene fusions, and we report the first cell type-specific expression of circular RNAs-a neuron-specific and nuclear-enriched RNA arising from the gene Hnrnpu.

CONCLUSIONS

These data provide an important resource for studies evaluating the function of various noncoding RNAs in the brain, including noncoding RNAs that may play a role in psychiatric disease.

ChimPipe: accurate detection of fusion genes and transcription-induced chimeras from RNA-seq data

Background

Chimeric transcripts are commonly defined as transcripts linking two or more different genes in the genome, and can be explained by various biological mechanisms such as genomic rearrangement, read-through or trans-splicing, but also by technical or biological artefacts. Several studies have shown their importance in cancer, cell pluripotency and motility. Many programs have recently been developed to identify chimeras from Illumina RNA-seq data (mostly fusion genes in cancer). However outputs of different programs on the same dataset can be widely inconsistent, and tend to include many false positives. Other issues relate to simulated datasets restricted to fusion genes, real datasets with limited numbers of validated cases, result inconsistencies between simulated and real datasets, and gene rather than junction level assessment.

Results

Here we present ChimPipe, a modular and easy-to-use method to reliably identify fusion genes and transcription-induced chimeras from paired-end Illumina RNA-seq data. We have also produced realistic simulated datasets for three different read lengths, and enhanced two gold-standard cancer datasets by associating exact junction points to validated gene fusions. Benchmarking ChimPipe together with four other state-of-the-art tools on this data showed ChimPipe to be the top program at identifying exact junction coordinates for both kinds of datasets, and the one showing the best trade-off between sensitivity and precision. Applied to 106 ENCODE human RNA-seq datasets, ChimPipe identified 137 high confidence chimeras connecting the protein coding sequence of their parent genes. In subsequent experiments, three out of four predicted chimeras, two of which recurrently expressed in a large majority of the samples, could be validated. Cloning and sequencing of the three cases revealed several new chimeric transcript structures, 3 of which with the potential to encode a chimeric protein for which we hypothesized a new role. Applying ChimPipe to human and mouse ENCODE RNA-seq data led to the identification of 131 recurrent chimeras common to both species, and therefore potentially conserved.

Conclusions

ChimPipe combines discordant paired-end reads and split-reads to detect any kind of chimeras, including those originating from polymerase read-through, and shows an excellent trade-off between sensitivity and precision. The chimeras found by ChimPipe can be validated in-vitro with high accuracy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3404-9) contains supplementary material, which is available to authorized users.

InFusion: Advancing Discovery of Fusion Genes and Chimeric Transcripts from Deep RNA-Sequencing Data

Analysis of fusion transcripts has become increasingly important due to their link with cancer development. Since high-throughput sequencing approaches survey fusion events exhaustively, several computational methods for the detection of gene fusions from RNA-seq data have been developed. This kind of analysis, however, is complicated by native trans-splicing events, the splicing-induced complexity of the transcriptome and biases and artefacts introduced in experiments and data analysis. There are a number of tools available for the detection of fusions from RNA-seq data; however, certain differences in specificity and sensitivity between commonly used approaches have been found. The ability to detect gene fusions of different types, including isoform fusions and fusions involving non-coding regions, has not been thoroughly studied yet. Here, we propose a novel computational toolkit called InFusion for fusion gene detection from RNA-seq data. InFusion introduces several unique features, such as discovery of fusions involving intergenic regions, and detection of anti-sense transcription in chimeric RNAs based on strand-specificity. Our approach demonstrates superior detection accuracy on simulated data and several public RNA-seq datasets. This improved performance was also evident when evaluating data from RNA deep-sequencing of two well-established prostate cancer cell lines. InFusion identified 26 novel fusion events that were validated in vitro, including alternatively spliced gene fusion isoforms and chimeric transcripts that include intergenic regions. The toolkit is freely available to download from http:/bitbucket.org/kokonech/infusion.

A computational analysis of the genetic and transcript diversity at the kallikrein locus

The kallikrein related peptidase gene family (KLKs) comprises 15 genes located between 19q13.3-13.4. KLKs have chymotrypsin and/or trypsin like activity, but the tissue/organ expression profile of each KLK varies considerably. Thus, the role of KLKs in human biology is also very diverse, and the deregulation of their function results in a wide-range of diseases. Here, we have cataloged the transcript (variants and fusions) and genetic (single nucleotide polymorphisms, small insertions/deletions, copy number variations (CNVs), and short tandem repeats) diversity at the KLK locus, providing a data set for researchers to explore the mechanisms through which KLK function may be deregulated. We reveal that the KLK locus hosts 85 fusion transcripts, and 80 variant transcripts. Interestingly, some fusion transcripts comprise up to 6 KLK genes. Our analysis of genetic variations of 2504 individuals from the 1000 Genome Project indicated that the KLK locus is rich in genetic diversity, with some fusion transcripts harboring over 1000 single nucleotide variations. We also found evidence from the literature linking 2387 KLK genetic variants with many types of diseases. Finally, genotyping data from the 131 KLK genetic variants in the NCI-60 cancer cell lines is provided as a resource for the cancer and KLK field.

ChiTaRS-3.1-the enhanced chimeric transcripts and RNA-seq database matched with protein-protein interactions

Discovery of chimeric RNAs, which are produced by chromosomal translocations as well as the joining of exons from different genes by trans-splicing, has added a new level of complexity to our study and understanding of the transcriptome. The enhanced ChiTaRS-3.1 database (http://chitars.md.biu.ac.il) is designed to make widely accessible a wealth of mined data on chimeric RNAs, with easy-to-use analytical tools built-in. The database comprises 34 922 chimeric transcripts along with 11 714 cancer breakpoints. In this latest version, we have included multiple cross-references to GeneCards, iHop, PubMed, NCBI, Ensembl, OMIM, RefSeq and the Mitelman collection for every entry in the ‘Full Collection’. In addition, for every chimera, we have added a predicted chimeric protein–protein interaction (ChiPPI) network, which allows for easy visualization of protein partners of both parental and fusion proteins for all human chimeras. The database contains a comprehensive annotation for 34 922 chimeric transcripts from eight organisms, and includes the manual annotation of 200 sense-antiSense (SaS) chimeras. The current improvements in the content and functionality to the ChiTaRS database make it a central resource for the study of chimeric transcripts and fusion proteins.

Overexpression of KLHL23 protein from read-through transcription of PHOSPHO2-KLHL23 in gastric cancer increases cell proliferation

Gene fusion, as a prototypical pathognomonic mutation, contributes to genome complexity, and the cis‐transcription‐induced gene fusions generated by read‐through transcription of adjacent genes have been found to be important for tumor development. We screened read‐through transcription events from stomach adenocarcinoma RNA‐seq data and selected three candidates PHOSPHO2‐KLHL23, RPL17‐C18orf32, and PRR5‐ARHGAP8, to assess their biological role in gastric cancer. The expression of all three read‐through fusion transcripts was confirmed in gastric cancer cell lines and paired normal/tumor gastric cancer tissues by real‐time quantitative reverse transcription polymerase chain reaction and their expression was found to be significantly higher in the tumor (P < 0.05; n = 75). The correlation between the expression level and clinicopathological information was statistically analyzed. The level of the PHOSPHO2‐KLHL23 read‐through fusion transcript correlated with the Lauren classification and was significantly associated with the presence of perineural invasion. Overexpression of KLHL23 from PHOSPHO2‐KLHL23 read‐through transcript led to a significant increase in cell proliferation and resistance to anticancer drug treatment. Silencing of KLHL23 expression decreased cyclin D1 levels. The expression of KLHL23 from prevalent read‐through transcripts of PHOSPHO2‐KLHL23 in gastric cancer may undermine the efficacy of anticancer drug treatment.

Intergenically Spliced Chimeric RNAs in Cancer

Gene fusions and their encoded products (fusion RNAs and proteins) are viewed as one of the hallmarks of cancer. Traditionally, they were thought to be generated solely by chromosomal rearrangements. However, recent discoveries of trans-splicing and cis-splicing events between neighboring genes, suggest that there are other mechanisms to generate chimeric fusion RNAs without corresponding changes in DNA. In addition, chimeric RNAs have been detected in normal physiology, complicating the use of fusions in cancer detection and therapy. On the other hand, "intergenically spliced" fusion RNAs represent a new repertoire of biomarkers and therapeutic targets. Here, we review current knowledge on chimeric RNAs and implications for cancer detection and treatment, and discuss outstanding questions for the advancement of the field.

Open questions in the study of de novo genes: what, how and why

The study of de novo protein-coding genes is maturing from the ad hoc reporting of individual cases to the systematic analysis of extensive genomic data from several species. We identify three key challenges for this emerging field: understanding how best to identify de novo genes, how they arise and why they spread. We highlight the intellectual challenges of understanding how a de novo gene becomes integrated into pre-existing functions and becomes essential. We suggest that, as with protein sequence evolution, antagonistic co-evolution may be key to de novo gene evolution, particularly for new essential genes and new cancer-associated genes.

Expression of evolutionarily novel genes in tumors

The evolutionarily novel genes originated through different molecular mechanisms are expressed in tumors. Sometimes the expression of evolutionarily novel genes in tumors is highly specific. Moreover positive selection of many human tumor-related genes in primate lineage suggests their involvement in the origin of new functions beneficial to organisms. It is suggested to consider the expression of evolutionarily young or novel genes in tumors as a new biological phenomenon, a phenomenon of TSEEN (tumor specifically expressed, evolutionarily novel) genes.

Recurrent and pathological gene fusions in breast cancer: current advances in genomic discovery and clinical implications

Gene fusions have long been considered principally as the oncogenic events of hematologic malignancies, but have recently gained wide attention in solid tumors due to several milestone discoveries and the advancement of deep sequencing technologies. With the progress in deep sequencing studies of breast cancer transcriptomes and genomes, the discovery of recurrent and pathological gene fusions in breast cancer is on the focus. Recently, driven by new deep sequencing studies, several recurrent or pathological gene fusions have been identified in breast cancer, including ESR1-CCDC170, SEC16A-NOTCH1, SEC22B-NOTCH2, and ESR1-YAP1 etc. More important, most of these gene fusions are preferentially identified in the more aggressive breast cancers, such as luminal B, basal-like, or endocrine-resistant breast cancer, suggesting recurrent gene fusions as additional key driver events in these tumors other than the known drivers such as the estrogen receptor. In this paper, we have comprehensively summarized the newly identified recurrent or pathological gene fusion events in breast cancer, reviewed the contributions of new genomic and deep sequencing technologies to new fusion discovery and the integrative bioinformatics tools to analyze these data, highlighted the biological relevance and clinical implications of these fusion discoveries, and discussed future directions of gene fusion research in breast cancer.

Transcriptome analysis revealed chimeric RNAs, single nucleotide polymorphisms and allele-specific expression in porcine prenatal skeletal muscle

Prenatal skeletal muscle development genetically determines postnatal muscle characteristics such as growth and meat quality in pigs. However, the molecular mechanisms underlying prenatal skeletal muscle development remain unclear. Here, we performed the first genome-wide analysis of chimeric RNAs, single nuclear polymorphisms (SNPs) and allele-specific expression (ASE) in prenatal skeletal muscle in pigs. We identified 14,810 protein coding genes and 163 high-confidence chimeric RNAs expressed in prenatal skeletal muscle. More than 94.5% of the chimeric RNAs obeyed the canonical GT/AG splice rule and were trans-splicing events. Ten and two RNAs were aligned to human and mouse chimeric transcripts, respectively. We detected 106,457 high-quality SNPs (6,955 novel), which were mostly (89.09%) located within QTLs for production traits. The high proportion of non-exonic SNPs revealed the incomplete annotation status of the current swine reference genome. ASE analysis revealed that 11,300 heterozygous SNPs showed allelic imbalance, whereas 131 ASE variants were located in the chimeric RNAs. Moreover, 4 ASE variants were associated with various economically relevant traits of pigs. Taken together, our data provide a source for studies of chimeric RNAs and biomarkers for pig breeding, while illuminating the complex transcriptional events underlying prenatal skeletal muscle development in mammals.

Recurrent cis-SAGe chimeric RNA, D2HGDH-GAL3ST2, in prostate cancer

Neighboring genes transcribing in the same direction can form chimeric RNAs via cis-splicing (cis-SAGe). Previously, we reported 16 novel cis-SAGe chimeras in prostate cancer cell lines, and performed in silico validation on 14 pairs of normal and tumor samples from Chinese patients. However, whether these fusions exist in different populations, as well as their clinical implications, remains unclear. To investigate, we developed a bioinformatics pipeline using modified Spliced Transcripts Alignment to a Reference (STAR) to quantify these fusion RNAs simultaneously in silico. From RNA-Seq data of 100 paired normal and prostate cancer samples from TCGA, we find that most fusions are not specific to cancer. However, D2HGDH-GAL3ST2 is more frequently seen in cancer samples, and seems to be enriched in the African American group. Further validation with our own collection as well as from commercial sources did not detect this fusion RNA in 29 normal prostate samples, but in 19 of 93 prostate cancer samples. It is more frequently detected in late stage cancer, suggesting a role in cancer progression. Consistently, silencing this fusion resulted in dramatic reduction of cell proliferation rate and cell motility.