Expression of conjoined genes: another mechanism for gene regulation in eukaryotes

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Longitudinal APOE4- and amyloid-dependent changes in the blood transcriptome in cognitively intact older adults

BACKGROUND

Gene expression is dysregulated in Alzheimer's disease (AD) patients, both in peripheral blood and post mortem brain. We investigated peripheral whole-blood gene (co)expression to determine molecular changes prior to symptom onset.

METHODS

RNA was extracted and sequenced for 65 cognitively healthy F-PACK participants (65 (56-80) years, 34 APOE4 non-carriers, 31 APOE4 carriers), at baseline and follow-up (interval: 5.0 (3.4-8.6) years). Participants received amyloid PET at both time points and amyloid rate of change derived. Accumulators were defined with rate of change ≥ 2.19 Centiloids. We performed differential gene expression and weighted gene co-expression network analysis to identify differentially expressed genes and networks of co-expressed genes, respectively, with respect to traits of interest (APOE4 status, amyloid accumulation (binary/continuous)), and amyloid positivity status, followed by Gene Ontology annotation.

RESULTS

There were 166 significant differentially expressed genes at follow-up compared to baseline in APOE4 carriers only, whereas 12 significant differentially expressed genes were found only in APOE4 non-carriers, over time. Among the significant genes in APOE4 carriers, several had strong evidence for a pathogenic role in AD based on direct association scores generated from the DISQOVER platform: NGRN, IGF2, GMPR, CLDN5, SMIM24. Top enrichment terms showed upregulated mitochondrial and metabolic pathways, and an exacerbated upregulation of ribosomal pathways in APOE4 carriers compared to non-carriers. Similarly, there were 33 unique significant differentially expressed genes at follow-up compared to baseline in individuals classified as amyloid negative at baseline and positive at follow-up or amyloid positive at both time points and 32 unique significant differentially expressed genes over time in individuals amyloid negative at both time points. Among the significant genes in the first group, the top five with the highest direct association scores were as follows: RPL17-C18orf32, HSP90AA1, MBP, SIRPB1, and GRINA. Top enrichment terms included upregulated metabolism and focal adhesion pathways. Baseline and follow-up gene co-expression networks were separately built. Seventeen baseline co-expression modules were derived, with one significantly negatively associated with amyloid accumulator status (r² =  - 0.25, p = 0.046). This was enriched for proteasomal protein catabolic process and myeloid cell development. Thirty-two follow-up modules were derived, with two significantly associated with APOE4 status: one downregulated (r² =  - 0.27, p = 0.035) and one upregulated (r² = 0.26, p = 0.039) module. Top enrichment processes for the downregulated module included proteasomal protein catabolic process and myeloid cell homeostasis. Top enrichment processes for the upregulated module included cytoplasmic translation and rRNA processing.

CONCLUSIONS

We show that there are longitudinal gene expression changes that implicate a disrupted immune system, protein removal, and metabolism in cognitively intact individuals who carry APOE4 or who accumulate in cortical amyloid. This provides insight into the pathophysiology of AD, whilst providing novel targets for drug and therapeutic development.

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.

Definition of the transcriptional units of inherited retinal disease genes by meta-analysis of human retinal transcriptome data

BACKGROUND

Inherited retinal diseases (IRD) are genetically heterogeneous disorders that cause the dysfunction or loss of photoreceptor cells and ultimately lead to blindness. To date, next-generation sequencing procedures fail to detect pathogenic sequence variants in coding regions of known IRD disease genes in about 30-40% of patients. One of the possible explanations for this missing heritability is the presence of yet unidentified transcripts of known IRD genes. Here, we aimed to define the transcript composition of IRD genes in the human retina by a meta-analysis of publicly available RNA-seq datasets using an ad-hoc designed pipeline.

RESULTS

We analysed 218 IRD genes and identified 5,054 transcripts, 3,367 of which were not previously reported. We assessed their putative expression levels and focused our attention on 435 transcripts predicted to account for at least 5% of the expression of the corresponding gene. We looked at the possible impact of the newly identified transcripts at the protein level and experimentally validated a subset of them.

CONCLUSIONS

This study provides an unprecedented, detailed overview of the complexity of the human retinal transcriptome that can be instrumental in contributing to the resolution of some cases of missing heritability in IRD patients.

Transcriptomic response of bioengineered human cartilage to parabolic flight microgravity is sex-dependent

Spaceflight and simulated spaceflight microgravity induced osteoarthritic-like alterations at the transcriptomic and proteomic levels in the articular and meniscal cartilages of rodents. But little is known about the effect of spaceflight or simulated spaceflight microgravity on the transcriptome of tissue-engineered cartilage developed from human cells. In this study, we investigate the effect of simulated spaceflight microgravity facilitated by parabolic flights on tissue-engineered cartilage developed from in vitro chondrogenesis of human bone marrow mesenchymal stem cells obtained from age-matched female and male donors. The successful induction of cartilage-like tissue was confirmed by the expression of well-demonstrated chondrogenic markers. Our bulk transcriptome data via RNA sequencing demonstrated that parabolic flight altered mostly fundamental biological processes, and the modulation of the transcriptome profile showed sex-dependent differences. The secretome profile analysis revealed that two genes (WNT7B and WNT9A) from the Wnt-signaling pathway, which is implicated in osteoarthritis development, were only up-regulated for female donors. The results of this study showed that the engineered cartilage tissues responded to microgravity in a sex-dependent manner, and the reported data offers a strong foundation to further explore the underlying mechanisms.

Skin Phototype and Disease: A Comprehensive Genetic Approach to Pigmentary Traits Pleiotropy Using PRS in the GCAT Cohort

Human pigmentation has largely been associated with different disease prevalence among populations, but most of these studies are observational and inconclusive. Known to be genetically determined, pigmentary traits have largely been studied by Genome-Wide Association Study (GWAS), mostly in Caucasian ancestry cohorts from North Europe, identifying robustly, several loci involved in many of the pigmentary traits. Here, we conduct a detailed analysis by GWAS and Polygenic Risk Score (PRS) of 13 pigmentary-related traits in a South European cohort of Caucasian ancestry (n = 20,000). We observed fair phototype strongly associated with non-melanoma skin cancer and other dermatoses and confirmed by PRS-approach the shared genetic basis with skin and eye diseases, such as melanoma (OR = 0.95), non-melanoma skin cancer (OR = 0.93), basal cell carcinoma (OR = 0.97) and darker phototype with vitiligo (OR = 1.02), cataracts (OR = 1.04). Detailed genetic analyses revealed 37 risk loci associated with 10 out of 13 analyzed traits, and 16 genes significantly associated with at least two pigmentary traits. Some of them have been widely reported, such as MC1R, HERC2, OCA2, TYR, TYRP1, SLC45A2, and some novel candidate genes C1QTNF3, LINC02876, and C1QTNF3-AMACR have not been reported in the GWAS Catalog, with regulatory potential. These results highlight the importance of the assess phototype as a genetic proxy of skin functionality and disease when evaluating open mixed populations.

Long-Read Transcriptome of Equine Bronchoalveolar Cells

We used Pacific Biosciences long-read isoform sequencing to generate full-length transcript sequences in equine bronchoalveolar lavage fluid (BALF) cells. Our dataset consisted of 313,563 HiFi reads comprising 805 Mb of polished sequence information. The resulting equine BALF transcriptome consisted of 14,234 full-length transcript isoforms originating from 7017 unique genes. These genes consisted of 6880 previously annotated genes and 137 novel genes. We identified 3428 novel transcripts in addition to 10,806 previously known transcripts. These included transcripts absent from existing genome annotations, transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. We provide transcript-level data for equine BALF cells as a resource to the scientific community.

Identification of novel prognostic risk signature of breast cancer based on ferroptosis-related genes

Ferroptosis is a type of cell regulated necrosis triggered by intracellular phospholipid peroxidation, which is more immunogenic than apoptosis. Therefore, genes controlling ferroptosis may be promising candidate biomarkers for tumor therapy. In this study, we investigate the function of genes associated with ferroptosis in breast cancer (BC) and systematically evaluate the relationship between ferroptosis-related gene expression and prognosis of BC patients from the Cancer Genome Atlas database. By using the consensus clustering method, 1203 breast cancer samples were clustered into two clearly divided subgroups based on the expression of 237 ferroptosis-related genes. Then differentially expressed analysis and least absolute shrinkage and selection operator were used to identify the prognosis-related genes. Furthermore, the genetic risk signature was constructed using the expression of prognosis-related genes. Our results showed that the genetic risk signature can identify patient subgroups with distinct prognosis in either training cohort or validation, and the genetic risk signature was associated with the tumor immune microenvironment. Finally, the Cox regression analysis indicated that our risk signature was an independent prognostic factor for BC patients and this signature was verified by the polymerase chain reaction and western blot. Within this study, we identified a novel prognostic classifier based on five ferroptosis-related genes which may provide a new reference for the treatment of BRCA patients.

Conjoined Genes as Common Events in Childhood Acute Lymphoblastic Leukemia

Simple Summary

Acute lymphoblastic leukemia (ALL) is the most frequent childhood cancer. In recent years, broad application of NGS technologies enabled the discovery of novel genomically defined ALL. In this study, as a proof-of-principle, we applied RNA-seq technology to comprehensively profile the transcriptional landscape of a collection of 10 childhood BCP-ALL cases, and performed a deep bioinformatics analysis including several publicly available datasets, in order to characterize their full spectrum of transcriptional events. The paired-end RNA sequencing of our BCP-ALL pediatric cohort revealed a total of 9001 raw fusion events, which, after filtering, resulted in 245 candidate fusions. Overall, 235 out of 245 events were intra-chromosomal fusions, among which 229 involved two contiguous or overlapping genes, also known as conjoined genes (CGs). Among them, we identified a subset of 14 CGs (6.1%) exclusively expressed in leukemic cases but neither in solid cancers nor in normal samples. These events could be suggestive of a novel mechanism of transcriptional regulation in childhood leukemia and may represent novel potential leukemia-specific biomarkers.

Abstract

Acute lymphoblastic leukemia (ALL) is the most frequent childhood cancer. For the last three decades, conventional cytogenetic and molecular approaches allowed the identification of genetic abnormalities having prognostic and therapeutic relevance. Although the current cure rate in pediatric B cell acute leukemia is approximately 90%, it remains one of the leading causes of mortality in childhood. Furthermore, in the contemporary protocols, chemotherapy intensity was raised to the maximal levels of tolerability, and further improvements in the outcome will depend on the characterization and reclassification of the disease, as well as on the development of new targeted drugs. The recent technological advances in genome-wide profiling techniques have allowed the exploration of the molecular heterogeneity of this disease, even though some potentially interesting biomarkers such as conjoined genes have not been deeply investigated yet. In the present study, we performed the transcriptome sequencing (RNA-seq) of 10 pediatric B cell precursor (BCP)-ALL cases with different risk (four standard- and six high-risk patients) enrolled in the Italian AIEOP-BFM ALL2000 protocol, in order to characterize the full spectrum of transcriptional events and to identify novel potential genetic mechanisms sustaining their different early response to therapy. Total RNA was extracted from primary leukemic blasts and RNA-seq was performed by Illumina technology. Bioinformatics analysis focused on fusion transcripts, originated from either inter- or intra-chromosomal structural rearrangements. Starting from a raw list of 9001 candidate events, by employing a custom-made bioinformatics pipeline, we obtained a short list of 245 candidate fusions. Among them, 10 events were compatible with chromosomal translocations. Strikingly, 235/245 events were intra-chromosomal fusions, 229 of which involved two contiguous or overlapping genes, resulting in the so-called conjoined genes (CGs). To explore the specificity of these events in leukemia, we performed an extensive bioinformatics meta-analysis and evaluated the presence of the fusions identified in our 10 BCP-ALL cohort in several other publicly available RNA-seq datasets, including leukemic, solid tumor and normal sample collections. Overall, 14/229 (6.1%) CGs were found to be exclusively expressed in leukemic cases, suggesting an association between CGs and leukemia. Moreover, CGs were found to be common events both in standard- and high-risk BCP-ALL patients and it might be suggestive of a novel potential transcriptional regulation mechanism active in leukemic cells.

GraphChrom: A Novel Graph-Based Framework for Cancer Classification Using Chromosomal Rearrangement Endpoints

Chromosomal rearrangements are generally a consequence of improperly repaired double-strand breaks in DNA. These genomic aberrations can be a driver of cancers. Here, we investigated the use of chromosomal rearrangements for classification of cancer tumors and the effect of inter- and intrachromosomal rearrangements in cancer classification. We used data from the Catalogue of Somatic Mutations in Cancer (COSMIC) for breast, pancreatic, and prostate cancers, for which the COSMIC dataset reports the highest number of chromosomal aberrations. We developed a framework known as GraphChrom for cancer classification. GraphChrom was developed using a graph neural network which models the complex structure of chromosomal aberrations (CA) and provides local connectivity between the aberrations. The proposed framework illustrates three important contributions to the field of cancers. Firstly, it successfully classifies cancer types and subtypes. Secondly, it evolved into a novel data extraction technique which can be used to extract more informative graphs (informative aberrations associated with a sample); and thirdly, it predicts that interCAs (rearrangements between two or more chromosomes) are more effective in cancer prediction than intraCAs (rearrangements within the same chromosome), although intraCAs are three times more likely to occur than intraCAs.

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.

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

The Diverse Roles of TNNI3K in Cardiac Disease and Potential for Treatment

In the two decades since the discovery of TNNI3K it has been implicated in multiple cardiac phenotypes and physiological processes. TNNI3K is an understudied kinase, which is mainly expressed in the heart. Human genetic variants in TNNI3K are associated with supraventricular arrhythmias, conduction disease, and cardiomyopathy. Furthermore, studies in mice implicate the gene in cardiac hypertrophy, cardiac regeneration, and recovery after ischemia/reperfusion injury. Several new papers on TNNI3K have been published since the last overview, broadening the clinical perspective of TNNI3K variants and our understanding of the underlying molecular biology. We here provide an overview of the role of TNNI3K in cardiomyopathy and arrhythmia covering both a clinical perspective and basic science advancements. In addition, we review the potential of TNNI3K as a target for clinical treatments in different cardiac diseases.

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.

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.

annoFuse: an R Package to annotate, prioritize, and interactively explore putative oncogenic RNA fusions

BACKGROUND

Gene fusion events are significant sources of somatic variation across adult and pediatric cancers and are some of the most clinically-effective therapeutic targets, yet low consensus of RNA-Seq fusion prediction algorithms makes therapeutic prioritization difficult. In addition, events such as polymerase read-throughs, mis-mapping due to gene homology, and fusions occurring in healthy normal tissue require informed filtering, making it difficult for researchers and clinicians to rapidly discern gene fusions that might be true underlying oncogenic drivers of a tumor and in some cases, appropriate targets for therapy.

RESULTS

We developed annoFuse, an R package, and shinyFuse, a companion web application, to annotate, prioritize, and explore biologically-relevant expressed gene fusions, downstream of fusion calling. We validated annoFuse using a random cohort of TCGA RNA-Seq samples (N = 160) and achieved a 96% sensitivity for retention of high-confidence fusions (N = 603). annoFuse uses FusionAnnotator annotations to filter non-oncogenic and/or artifactual fusions. Then, fusions are prioritized if previously reported in TCGA and/or fusions containing gene partners that are known oncogenes, tumor suppressor genes, COSMIC genes, and/or transcription factors. We applied annoFuse to fusion calls from pediatric brain tumor RNA-Seq samples (N = 1028) provided as part of the Open Pediatric Brain Tumor Atlas (OpenPBTA) Project to determine recurrent fusions and recurrently-fused genes within different brain tumor histologies. annoFuse annotates protein domains using the PFAM database, assesses reciprocality, and annotates gene partners for kinase domain retention. As a standard function, reportFuse enables generation of a reproducible R Markdown report to summarize filtered fusions, visualize breakpoints and protein domains by transcript, and plot recurrent fusions within cohorts. Finally, we created shinyFuse for algorithm-agnostic interactive exploration and plotting of gene fusions.

CONCLUSIONS

annoFuse provides standardized filtering and annotation for gene fusion calls from STAR-Fusion and Arriba by merging, filtering, and prioritizing putative oncogenic fusions across large cancer datasets, as demonstrated here with data from the OpenPBTA project. We are expanding the package to be widely-applicable to other fusion algorithms and expect annoFuse to provide researchers a method for rapidly evaluating, prioritizing, and translating fusion findings in patient tumors.

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.

Emergent Coordination of the CHKB and CPT1B Genes in Eutherian Mammals: Implications for the Origin of Brown Adipose Tissue

Mitochondrial fatty acid oxidation (FAO) contributes to the proton motive force that drives ATP synthesis in many mammalian tissues. In eutherian (placental) mammals, brown adipose tissue (BAT) can also dissipate this proton gradient through uncoupling protein 1 (UCP1) to generate heat, but the evolutionary events underlying the emergence of BAT are unknown. An essential step in FAO is the transport of cytoplasmic long chain acyl-coenzyme A (acyl-CoA) into the mitochondrial matrix, which requires the action of carnitine palmitoyltransferase 1B (CPT1B) in striated muscle and BAT. In eutherians, the CPT1B gene is closely linked to the choline kinase beta (CHKB) gene, which is transcribed from the same DNA strand and terminates just upstream of CPT1B. CHKB is a rate-limiting enzyme in the synthesis of phosphatidylcholine (PC), a predominant mitochondrial membrane phospholipid, suggesting that the coordinated expression of CHKB and CPT1B may cooperatively enhance mitochondrial FAO. The present findings show that transcription of the eutherian CHKB and CPT1B genes is linked within a unitary epigenetic domain targeted to the CHKB gene, and that that this regulatory linkage appears to have resulted from an intergenic deletion in eutherians that significantly altered the distribution of CHKB and CPT1B expression. Informed by the timing of this event relative to the emergence of BAT, the phylogeny of CHKB-CPT1B synteny, and the insufficiency of UCP1 to account for eutherian BAT, these data support a mechanism for the emergence of BAT based on the acquisition of a novel capacity for adipocyte FAO in a background of extant UCP1.

Progress, Challenges, and Surprises in Annotating the Human Genome

Our understanding of the human genome has continuously expanded since its draft publication in 2001. Over the years, novel assays have allowed us to progressively overlay layers of knowledge above the raw sequence of A's, T's, G's, and C's. The reference human genome sequence is now a complex knowledge base maintained under the shared stewardship of multiple specialist communities. Its complexity stems from the fact that it is simultaneously a template for transcription, a record of evolution, a vehicle for genetics, and a functional molecule. In short, the human genome serves as a frame of reference at the intersection of a diversity of scientific fields. In recent years, the progressive fall in sequencing costs has given increasing importance to the quality of the human reference genome, as hundreds of thousands of individuals are being sequenced yearly, often for clinical applications. Also, novel sequencing-based assays shed light on novel functions of the genome, especially with respect to gene expression regulation. Keeping the human genome annotation up to date and accurate is therefore an ongoing partnership between reference annotation projects and the greater community worldwide.

Arsenic Metabolism in Mice Carrying a BORCS7/AS3MT Locus Humanized by Syntenic Replacement

BACKGROUND

Chronic exposure to inorganic arsenic (iAs) is a significant public health problem. Methylation of iAs by arsenic methyltransferase (AS3MT) controls iAs detoxification and modifies risks of iAs-induced diseases. Mechanisms underlying these diseases have been extensively studied using animal models. However, substantive differences between humans and laboratory animals in efficiency of iAs methylation have hindered the translational potential of the laboratory studies.

OBJECTIVES

The goal of this study was to determine whether humanization of the As3mt gene confers a human-like pattern of iAs metabolism in mice.

METHODS

We generated a mouse strain in which the As3mt gene along with the adjacent Borcs7 gene was humanized by syntenic replacement. We compared expression of the mouse As3mt and the human AS3MT and the rate and pattern of iAs metabolism in the wild-type and humanized mice.

RESULTS

AS3MT expression in mouse tissues closely modeled that of human and differed substantially from expression of As3mt. Detoxification of iAs was much less efficient in the humanized mice than in wild-type mice. Profiles for iAs and its methylated metabolites in tissues and excreta of the humanized mice were consistent with those reported in humans. Notably, the humanized mice expressed both the full-length AS3MT that catalyzes iAs methylation and the human-specific AS3MTd2d3 splicing variant that has been linked to schizophrenia.

CONCLUSIONS

These results suggest that AS3MT is the primary genetic locus responsible for the unique pattern of iAs metabolism in humans. Thus, the humanized mouse strain can be used to study the role of iAs methylation in the pathogenesis of iAs-induced diseases, as well as to evaluate the role of AS3MTd2d3 in schizophrenia. https://doi.org/10.1289/EHP6943.

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.

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 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.

A Simple RNA Target Capture NGS Strategy for Fusion Genes Assessment in the Diagnostics of Pediatric B-cell Acute Lymphoblastic Leukemia

Supplemental Digital Content is available in the text

Acute lymphoblastic leukemia (ALL) is the most frequent pediatric cancer. Fusion genes are hallmarks of ALL, and they are used as biomarkers for risk stratification as well as targets for precision medicine. Hence, clinical diagnostics pursues broad and comprehensive strategies for accurate discovery of fusion genes. Currently, the gold standard methodologies for fusion gene detection are fluorescence in situ hybridization and polymerase chain reaction; these, however, lack sensitivity for the identification of new fusion genes and breakpoints. In this study, we implemented a simple operating procedure (OP) for detecting fusion genes. The OP employs RNA CaptureSeq, a versatile and effortless next-generation sequencing assay, and an in-house as well as a purpose-built bioinformatics pipeline for the subsequent data analysis. The OP was evaluated on a cohort of 89 B-cell precursor ALL (BCP-ALL) pediatric samples annotated as negative for fusion genes by the standard techniques. The OP confirmed 51 samples as negative for fusion genes, and, more importantly, it identified known (KMT2A rearrangements) as well as new fusion events (JAK2 rearrangements) in the remaining 38 investigated samples, of which 16 fusion genes had prognostic significance. Herein, we describe the OP and its deployment into routine ALL diagnostics, which will allow substantial improvements in both patient risk stratification and precision medicine.

Detection of novel fusion-transcripts by RNA-Seq in T-cell lymphoblastic lymphoma

Fusions transcripts have been proven to be strong drivers for neoplasia-associated mutations, although their incidence in T-cell lymphoblastic lymphoma needs to be determined yet. Using RNA-Seq we have selected 55 fusion transcripts identified by at least two of three detection methods in the same tumour. We confirmed the existence of 24 predicted novel fusions that had not been described in cancer or normal tissues yet, indicating the accuracy of the prediction. Of note, one of them involves the proto oncogene TAL1. Other confirmed fusions could explain the overexpression of driver genes such as COMMD3-BMI1, LMO1 or JAK3. Five fusions found exclusively in tumour samples could be considered pathogenic (NFYG-TAL1, RIC3-TCRBC2, SLC35A3-HIAT1, PICALM MLLT10 and MLLT10-PICALM). However, other fusions detected simultaneously in normal and tumour samples (JAK3-INSL3, KANSL1-ARL17A/B and TFG-ADGRG7) could be germ-line fusions genes involved in tumour-maintaining tasks. Notably, some fusions were confirmed in more tumour samples than predicted, indicating that the detection methods underestimated the real number of existing fusions. Our results highlight the potential of RNA-Seq to identify new cryptic fusions, which could be drivers or tumour-maintaining passenger genes. Such novel findings shed light on the searching for new T-LBL biomarkers in these haematological disorders.

Decoding Transcriptome Dynamics of Genome-Encoded Polyadenylation and Autoregulation with Small-Molecule Modulators of Alternative Polyadenylation

Alternative polyadenylation (APA) plays a critical role in regulating gene expression. However, the balance between genome-encoded APA processing and autoregulation by APA modulating RNA binding protein (RBP) factors is not well understood. We discovered two potent small-molecule modulators of APA (T4 and T5) that promote distal-to-proximal (DtoP) APA usage in multiple transcripts. Monotonically responsive APA events, induced by short exposure to T4 or T5, were defined in the transcriptome, allowing clear isolation of the genomic sequence features and RBP motifs associated with DtoP regulation. We found that longer vulnerable introns, enriched with distinctive A-rich motifs, were preferentially affected by DtoP APA, thus defining a core set of genes with genomically encoded DtoP regulation. Through APA response pattern and compound-small interfering RNA epistasis analysis of APA-associated RBP factors, we further demonstrated that DtoP APA usage is partly modulated by altered autoregulation of polyadenylate binding nuclear protein-1 signaling.

A critical role of the thioredoxin domain containing protein 5 (TXNDC5) in redox homeostasis and cancer development

Correct folding of nascent peptides occurs in the endoplasmic reticulum (ER). It is a complicate process primarily accomplished by the coordination of multiple redox proteins including members of the protein disulfide isomerase (PDI) family. As a critical member of the PDI family, thioredoxin domain containing protein 5 (TXNDC5) assists the folding of newly synthesized peptides to their mature form through series of disulfide bond exchange reactions. Interestingly, TXNDC5 is frequently found overexpressed in specimens of many human diseases including various types of cancer. In this review, we summarized the biochemical function of TXNDC5 in mammalian cells and the recent progress on the understanding of its role and molecular mechanisms in cancer development. Findings of TXNDC5 in the activation of intracellular signaling pathways, stimulation of cell growth & proliferation, facilitation of cell survival and modulation of extracellular matrix to affect cancer cell invasion and metastasis are reviewed. These published studies suggest that strategies of targeting TXNDC5 can be developed as potentially valuable methods for the treatment of certain types of cancer in patients.

Genetic variants in nucleotide excision repair pathway predict survival of esophageal squamous cell cancer patients receiving platinum-based chemotherapy

The benefits of platinum-based chemotherapy (PBC) on survival of esophageal squamous cell carcinoma (ESCC) patients are inexplicit due to the varied therapeutic effects. Nucleotide excision repair (NER) pathway plays a vital role in removing platinum-DNA adducts in tumor cells and hence may modulate the therapeutic effect and survival outcome. The present study assessed the associations of 26 potentially functional regulatory single nucleotide polymorphisms (rSNPs) in nine core NER genes with disease-free survival (DFS) and overall survival (OS) in 339 ESCC patients. We found that ERCC2 rs2097215 T and rs3916788 A, ERCC5 rs3759497 A and XPC rs3731054 C alleles were associated with unfavorable DFS. Patients carrying high-risk allele group (HRG, 5-8 risk alleles) had a significantly shorter DFS, compared with those carrying low-risk alleles (LRG, 0-4 risk alleles) [adjusted hazards ratio (HR_adj ) = 1.64, 95%CI = 1.23-2.19, P_adj < 0.001]. Three of these SNPs (ie, ERCC2 rs2097215 T and rs3916788 A and ERCC5 rs3759497 A) were also significantly associated with a poorer OS (HRG vs LRG: HR_adj  = 1.75, 95%CI = 1.23-2.47, P_adj  = 0.002). The expression quantitative trait loci (eQTL) analysis revealed significant genotype-expression correlations for ERCC5 rs3759497 and ERCC2 2097215 and rs3916788, which suggest regulatory roles of these SNPs. It appears that these NER variants may independently or jointly exert an impact on survival outcome of Chinese ESCC patients undergoing adjuvant platinum-based therapy. Large studies are warranted to validate these findings.

FusionHub: A unified web platform for annotation and visualization of gene fusion events in human cancer

Gene fusion is a chromosomal rearrangement event which plays a significant role in cancer due to the oncogenic potential of the chimeric protein generated through fusions. At present many databases are available in public domain which provides detailed information about known gene fusion events and their functional role. Existing gene fusion detection tools, based on analysis of transcriptomics data usually report a large number of fusion genes as potential candidates, which could be either known or novel or false positives. Manual annotation of these putative genes is indeed time-consuming. We have developed a web platform FusionHub, which acts as integrated search engine interfacing various fusion gene databases and simplifies large scale annotation of fusion genes in a seamless way. In addition, FusionHub provides three ways of visualizing fusion events: circular view, domain architecture view and network view. Design of potential siRNA molecules through ensemble method is another utility integrated in FusionHub that could aid in siRNA-based targeted therapy. FusionHub is freely available at https://fusionhub.persistent.co.in.

Perspective Insight into Future Potential Fusion Gene Transcript Biomarker Candidates in Breast Cancer

Next generation sequencing has accelerated the discovery of a variety of new fusion gene types in clinical breast cancer samples by analyzing cancer genomes and transcriptomes. Although previous studies have focused on a few clinically validated oncogenic fusion genes as diagnostic and therapeutic targets in breast cancer, a perspective consideration has not been given thus far for a plethora of breast cancer fusion genes, which are being newly identified at an overwhelmingly increasing pace. In this perspective review, we discuss diverse fusion gene types recently identified in a variety of breast cancer subtypes, including breast clinical cancer samples in TCGA (The Cancer Genome Atlas) database. This perspective review will confer fresh and promising guidance onto breast cancer surgeons, clinical oncologists, and tumor biologists in determining research directions for seeking and developing novel fusion gene biomarkers for breast cancer diagnostics and therapeutic treatment in upcoming years.

Absence of Correlation between Chimeric RNA and Aging

Chimeric RNAs have been recognized as a phenomenon not unique to cancer cells. They also exist in normal physiology. Aging is often characterized by deregulation of molecular and cellular mechanisms, including loss of heterochromatin, increased transcriptional noise, less tight control on alternative splicing, and more stress-induced changes. It is thus assumed that chimeric RNAs are more abundant in older people. In this study, we conducted a preliminary investigation to identify any chimeric RNAs with age-based trends in their expression levels in blood samples. A chimeric RNA candidate list generated by bioinformatic analysis indicated the possibility of both negative and positive trends in the expression of chimeric RNAs. Out of this candidate list, five novel chimeric RNAs were successfully amplified in multiple blood samples and then sequenced. Although primary smaller sample sizes displayed some weak trends with respect to age, analysis of quantitative PCR data from larger sample sizes showed essentially no relationship between expression levels and age. Altogether, these results indicate that, contradictory to the common assumption, chimeric RNAs as a group are not all higher in older individuals and that placing chimeric RNAs in the context of aging will be a much more complex task than initially anticipated.

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.

Conserved Gene Microsynteny Unveils Functional Interaction Between Protein Disulfide Isomerase and Rho Guanine-Dissociation Inhibitor Families

Protein disulfide isomerases (PDIs) support endoplasmic reticulum redox protein folding and cell-surface thiol-redox control of thrombosis and vascular remodeling. The family prototype PDIA1 regulates NADPH oxidase signaling and cytoskeleton organization, however the related underlying mechanisms are unclear. Here we show that genes encoding human PDIA1 and its two paralogs PDIA8 and PDIA2 are each flanked by genes encoding Rho guanine-dissociation inhibitors (GDI), known regulators of RhoGTPases/cytoskeleton. Evolutionary histories of these three microsyntenic regions reveal their emergence by two successive duplication events of a primordial gene pair in the last common vertebrate ancestor. The arrangement, however, is substantially older, detectable in echinoderms, nematodes, and cnidarians. Thus, PDI/RhoGDI pairing in the same transcription orientation emerged early in animal evolution and has been largely maintained. PDI/RhoGDI pairs are embedded into conserved genomic regions displaying common cis-regulatory elements. Analysis of gene expression datasets supports evidence for PDI/RhoGDI coexpression in developmental/inflammatory contexts. PDIA1/RhoGDIα were co-induced in endothelial cells upon CRISP-R-promoted transcription activation of each pair component, and also in mouse arterial intima during flow-induced remodeling. We provide evidence for physical interaction between both proteins. These data support strong functional links between PDI and RhoGDI families, which likely maintained PDI/RhoGDI microsynteny along > 800-million years of evolution.

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.

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.

Alternative exon definition events control the choice between nuclear retention and cytoplasmic export of U11/U12-65K mRNA

Cellular homeostasis of the minor spliceosome is regulated by a negative feed-back loop that targets U11-48K and U11/U12-65K mRNAs encoding essential components of the U12-type intron-specific U11/U12 di-snRNP. This involves interaction of the U11 snRNP with an evolutionarily conserved splicing enhancer giving rise to unproductive mRNA isoforms. In the case of U11/U12-65K, this mechanism controls the length of the 3′ untranslated region (3′UTR). We show that this process is dynamically regulated in developing neurons and some other cell types, and involves a binary switch between translation-competent mRNAs with a short 3′UTR to non-productive isoforms with a long 3′UTR that are retained in the nucleus or/and spliced to the downstream amylase locus. Importantly, the choice between these alternatives is determined by alternative terminal exon definition events regulated by conserved U12- and U2-type 5′ splice sites as well as sequence signals used for pre-mRNA cleavage and polyadenylation. We additionally show that U11 snRNP binding to the U11/U12-65K mRNA species with a long 3′UTR is required for their nuclear retention. Together, our studies uncover an intricate molecular circuitry regulating the abundance of a key spliceosomal protein and shed new light on the mechanisms limiting the export of non-productively spliced mRNAs from the nucleus to the cytoplasm.

The cellular homeostasis of many components of the eukaryotic RNA processing machinery is regulated via negative feed-back pathways that result in the formation of both productive and non-productive mRNA species. Typically, the formation of non-productive mRNAs species results from changes in alternative splicing that disrupt the reading frame of the protein coding region and leads to destabilization of the mRNA. Here, we have investigated the homeostasis regulation of the U11/U12-65K mRNA that encodes an essential protein component of the minor (U12-dependent) spliceosome intron recognition complex. We show that homeostasis is regulated at the level of nuclear mRNA export and mRNA 3′-end formation, and that it can be further regulated during neuronal differentiation. We describe a multilayered regulatory system utilizing alternative exon definition interactions that use the input from both spliceosomes and the polyadenylation machinery to decide between productive and non-productive mRNA formation. Because the 65K protein is an essential component of the minor spliceosome, this regulatory pathway can potentially affect the expression of ~700 genes containing U12-type introns.

CLK-dependent exon recognition and conjoined gene formation revealed with a novel small molecule inhibitor

CDC-like kinase phosphorylation of serine/arginine-rich proteins is central to RNA splicing reactions. Yet, the genomic network of CDC-like kinase-dependent RNA processing events remains poorly defined. Here, we explore the connectivity of genomic CDC-like kinase splicing functions by applying graduated, short-exposure, pharmacological CDC-like kinase inhibition using a novel small molecule (T3) with very high potency, selectivity, and cell-based stability. Using RNA-Seq, we define CDC-like kinase-responsive alternative splicing events, the large majority of which monotonically increase or decrease with increasing CDC-like kinase inhibition. We show that distinct RNA-binding motifs are associated with T3 response in skipped exons. Unexpectedly, we observe dose-dependent conjoined gene transcription, which is associated with motif enrichment in the last and second exons of upstream and downstream partners, respectively. siRNA knockdown of CLK2-associated genes significantly increases conjoined gene formation. Collectively, our results reveal an unexpected role for CDC-like kinase in conjoined gene formation, via regulation of 3′-end processing and associated splicing factors.

The phosphorylation of serine/arginine-rich proteins by CDC-like kinase is a central regulatory mechanism for RNA splicing reactions. Here, the authors synthesize a novel small molecule CLK inhibitor and map CLK-responsive alternative splicing events and discover an effect on conjoined gene transcription.

A novel read-through transcript JMJD7-PLA2G4B regulates head and neck squamous cell carcinoma cell proliferation and survival

Recent findings on the existence of oncogenic fusion genes in a wide array of solid tumors, including head and neck squamous cell carcinoma (HNSCC), suggests that fusion genes have become attractive targets for cancer diagnosis and treatment. In this study, we showed for the first time that a read-through fusion gene JMJD7-PLA2G4B is presented in HNSCC, splicing neighboring jumonji domain containing 7 (JMJD7) and phospholipase A2, group IVB (PLA2G4B) genes together. Ablation of JMJD7-PLA2G4B significantly inhibited proliferation of HNSCC cells by promoting G1 cell cycle arrest and increased starvation-induced cell death compared to JMJD7-only knockdown HNSCC cells. Mechanistically, we found that JMJD7-PLA2G4B modulates phosphorylation of Protein Kinase B (AKT) to promote HNSCC cell survival. Moreover, JMJD7-PLA2G4B also regulated an E3 ligase S-phase kinase-associated protein 2 (SKP2) to control the cell cycle progression from G1 phase to S phase by inhibiting Cyclin-dependent kinase inhibitor 1 (p21) and 1B (p27) expression. Our study provides novel insights into the oncogenic control of JMJD7-PLA2G4B in HNSCC cell proliferation and survival, and suggests that JMJD7-PLA2G4B may serve as an important therapeutic target and prognostic marker for HNSCC development and progression.

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.

Integrated genomic analyses identify frequent gene fusion events and VHL inactivation in gastrointestinal stromal tumors

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors of the gastrointestinal tract. We sequenced nine exomes and transcriptomes, and two genomes of GISTs for integrated analyses. We detected 306 somatic variants in nine GISTs and recurrent protein-altering mutations in 29 genes. Transcriptome sequencing revealed 328 gene fusions, and the most frequently involved fusion events were associated with IGF2 fused to several partner genes including CCND1, FUS, and LASP1. We additionally identified three recurrent read-through fusion transcripts: POLA2-CDC42EP2, C8orf42-FBXO25, and STX16-NPEPL1. Notably, we found intragenic deletions in one of three exons of the VHL gene and increased mRNAs of VEGF, PDGF-β, and IGF-1/2 in 56% of GISTs, suggesting a mechanistic link between VHL inactivation and overexpression of hypoxia-inducible factor target genes in the absence of hypoxia. We also identified copy number gain and increased mRNA expression of AMACR, CRIM1, SKP2, and CACNA1E. Mapping of copy number and gene expression results to the KEGG pathways revealed activation of the JAK-STAT pathway in small intestinal GISTs and the MAPK pathway in wild-type GISTs. These observations will allow us to determine the genetic basis of GISTs and will facilitate further investigation to develop new therapeutic options.

Upregulation of INS-IGF2 read-through expression and identification of a novel INS-IGF2 splice variant in insulinomas

Fusion transcripts arising from the combination of exons residing on neighboring genes on the same chromosome may give rise to chimeric or novel proteins. Such read-through transcripts have been detected in different cancers where they may be of pathogenetic interest. In this study, we describe for the first time the expression of a read-through transcript in insulinomas, a functioning neuroendocrine pancreatic neoplasm. The read-through transcript INS-IGF2, composed of exons from the two genes proinsulin precursor (INS) and insulin‑like growth factor 2 (IGF2), both mapping to chromosomal subband 11p15.5, was highly expressed in the two insulinomas analyzed. More precisely, version 2 of the INS-IGF2 transcript was expressed, indicating possible expression of the chimeric INS-IGF2 protein. We further identified a novel splice variant of the INS-IGF2 read-through transcript in one of the insulinomas, composed of exon 1 of INS3 and exons of IGF2. In the same tumor, we found high expression of INS3 and the presence of the A allele at SNP rs689. SNP rs689 has been previously described to regulate splicing of the INS transcript, indicating that this regulatory mechanism also affects splicing of INS-IGF2. The identification of the INS-IGF2 read-through transcript specifically in tumor tissue but not in normal pancreatic tissue suggests that high expression of INS-IGF2 could be neoplasia‑specific. These results may have potential clinical applications given that the read-through transcript could be used as a biomarker in insulinoma patients.

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.

Novel fusion transcripts associate with progressive prostate cancer

The mechanisms underlying the potential for aggressive behavior of prostate cancer (PCa) remain elusive. In this study, whole genome and/or transcriptome sequencing was performed on 19 specimens of PCa, matched adjacent benign prostate tissues, matched blood specimens, and organ donor prostates. A set of novel fusion transcripts was discovered in PCa. Eight of these fusion transcripts were validated through multiple approaches. The occurrence of these fusion transcripts was then analyzed in 289 prostate samples from three institutes, with clinical follow-up ranging from 1 to 15 years. The analyses indicated that most patients [69 (91%) of 76] positive for any of these fusion transcripts (TRMT11-GRIK2, SLC45A2-AMACR, MTOR-TP53BP1, LRRC59-FLJ60017, TMEM135-CCDC67, KDM4-AC011523.2, MAN2A1-FER, and CCNH-C5orf30) experienced PCa recurrence, metastases, and/or PCa-specific death after radical prostatectomy. These outcomes occurred in only 37% (58/157) of patients without carrying those fusion transcripts. Three fusion transcripts occurred exclusively in PCa samples from patients who experienced recurrence or PCaerelated death. The formation of these fusion transcripts may be the result of genome recombination. A combination of these fusion transcripts in PCa with Gleason's grading or with nomogram significantly improves the prediction rate of PCa recurrence. Our analyses suggest that formation of these fusion transcripts may underlie the aggressive behavior of PCa.