Lariat sequencing in a unicellular yeast identifies regulated alternative splicing of exons that are evolutionarily conserved with humans

Emerging roles of noncoding RNAs in cardiovascular pathophysiology

This review comprehensively examines the diverse roles of noncoding RNAs (ncRNAs) in the pathogenesis and treatment of cardiovascular disease (CVD), focusing on microRNA (miRNA), long noncoding RNA (lncRNA), piwi-interacting RNA (piRNA), small-interfering RNA (siRNA), circular RNA (circRNA), and vesicle-associated RNAs. These ncRNAs are integral regulators of key cellular processes, including gene expression, inflammation, and fibrosis, and they hold great potential as both diagnostic biomarkers and therapeutic targets. The review highlights novel insights into how these RNA species, particularly miRNAs, lncRNAs, and piRNAs, contribute to various CVDs such as hypertension, atherosclerosis, and myocardial infarction. In addition, it explores the emerging role of extracellular vesicles (EVs) in intercellular communication and their therapeutic potential in cardiovascular health. The review underscores the need for continued research into ncRNAs and RNA-based therapies, with a focus on advancing delivery systems and expanding personalized medicine approaches to improve cardiovascular outcomes.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

The fission yeast methyl phosphate capping enzyme Bmc1 guides 2'-O-methylation of the U6 snRNA

Splicing requires the tight coordination of dynamic spliceosomal RNAs and proteins. U6 is the only spliceosomal RNA transcribed by RNA Polymerase III and undergoes an extensive maturation process. In humans and fission yeast, this includes addition of a 5' γ-monomethyl phosphate cap by members of the Bin3/MePCE family as well as snoRNA guided 2'-O-methylation. Previously, we have shown that the Bin3/MePCE homolog Bmc1 is recruited to the S. pombe telomerase holoenzyme by the LARP7 family protein Pof8, where it acts in a catalytic-independent manner to protect the telomerase RNA and facilitate holoenzyme assembly. Here, we show that Bmc1 and Pof8 are required for the formation of a distinct U6 snRNP that promotes 2'-O-methylation of U6, and identify a non-canonical snoRNA that guides this methylation. We also show that the 5' γ-monomethyl phosphate capping activity of Bmc1 is not required for its role in promoting snoRNA guided 2'-O-methylation, and that this role relies on different regions of Pof8 from those required for Pof8 function in telomerase. Our results are consistent with a novel role for Bmc1/MePCE family members in stimulating 2'-O-methylation and a more general role for Bmc1 and Pof8 in guiding noncoding RNP assembly beyond the telomerase RNP.

The Role of Circular RNAs in Ischemic Stroke

Ischemic stroke (IS), a devastating condition characterized by intracranial artery stenosis and middle cerebral artery occlusion leading to insufficient oxygen supply to the brain, is a major cause of death and physical disability worldwide. Recent research has demonstrated the critical role of circular RNAs (circRNAs), a class of covalently enclosed noncoding RNAs that are widespread in eukaryotic cells, in regulating various physiological and pathophysiological cellular processes, including cell apoptosis, autophagy, synaptic plasticity, and neuroinflammation. In the past few years, circRNAs have attracted extensive attention in the field of IS research. This review summarizes the current understanding of the mechanisms underlying the involvement of circRNAs in IS development. A better understanding of circRNA-mediated pathogenic mechanisms in IS may pave the way for translating circRNA research into clinical practice, ultimately improving the clinical outcomes of IS patients.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Native RNA sequencing in fission yeast reveals frequent alternative splicing isoforms

The unicellular yeast Schizosaccharomyces pombe (fission yeast) retains many of the splicing features observed in humans and is thus an excellent model to study the basic mechanisms of splicing. Nearly half the genes contain introns, but the impact of alternative splicing in gene regulation and proteome diversification remains largely unexplored. Here we leverage Oxford Nanopore Technologies native RNA sequencing (dRNA), as well as ribosome profiling data, to uncover the full range of polyadenylated transcripts and translated open reading frames. We identify 332 alternative isoforms affecting the coding sequences of 262 different genes, 97 of which occur at frequencies >20%, indicating that functional alternative splicing in S. pombe is more prevalent than previously suspected. Intron retention events make ∼80% of the cases; these events may be involved in the regulation of gene expression and, in some cases, generate novel protein isoforms, as supported by ribosome profiling data in 18 of the intron retention isoforms. One example is the rpl22 gene, in which intron retention is associated with the translation of a protein of only 13 amino acids. We also find that lowly expressed transcripts tend to have longer poly(A) tails than highly expressed transcripts, highlighting an interdependence between poly(A) tail length and transcript expression level. Finally, we discover 214 novel transcripts that are not annotated, including 158 antisense transcripts, some of which also show translation evidence. The methodologies described in this work open new opportunities to study the regulation of splicing in a simple eukaryotic model.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Circulating expression levels of CircHIPK3 and CDR1as circular-RNAs in type 2 diabetes patients

BACKGROUND

Recent investigations suggested that deregulated levels of Circular RNAs (circRNAs) could be associated with type 2 diabetes mellitus (T2DM) pathogenesis. Accordingly, this study aimed to determine the expression levels of circulating CircHIPK3, CDR1as and their correlation with biochemical parameters in patients with T2DM, pre-diabetes and control subjects.

METHODS AND RESULTS

The expression of circRNAs in peripheral blood was determined using QRT-PCR in 70 patients with T2DM, 60 pre-diabetes and in 69 age and sex matched healthy controls. Moreover, bioinformatics tools were applied to explore and predict the potential interactions between circRNAs and other non-coding RNAs (ncRNAs). Our analysis revealed that the expression level of CircHIPK3 was significantly elevated in T2DM patients compared to healthy participants (P < 0.001) and pre-diabetes subjects (P = 0.018). In addition, ROC analysis suggested that at the cutoff value of 0.24 and the sensitivity and specificity of 50% and 88.4%, respectively, CircHIPK3 could distinguish between T2DM patients and control subjects. Furthermore, it was observed that the expression level of CDR1as is higher in pre-diabetic individuals than healthy individuals (P = 0.004). Finally, Spearman correlation analysis showed that there was a significant correlation between CircHIPK3 and CDR1as expression levels and clinical and anthropometrical parameters such as BMI, systolic and diastolic blood pressure, HbA1c and fasting blood glucose (P < 0.005).

CONCLUSIONS

The data of this study provided evidence that the expression levels of CircHIPK3, CDR1as increased in T2DM and pre-diabetes subjects, respectively.

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

Ascorbic acid supplementation suppresses cadmium-derived alterations in the fission yeast Schizosaccharomyces pombe

Cadmium (Cd) a highly toxic environmental pollutant, that does not have any physiological function in the organism, represents a great concern for human health as it can be easily transported from its environmental sources to the food chain. Food, water, and air are the major sources of Cd exposure to the population. Cd-mediated impairments of the basic cellular properties largely depend on its ability to enhance the formation of reactive oxygen species (ROS) and thus triggers oxidative stress to the cell. With the use of fission yeast Schizosaccharomyces pombe (S. pombe) as the model organism, we have analyzed the impact of Cd on the cell growth intensity, as it represents the fundamental feature of all living organisms. Cells were incubated with different Cd concentrations for 3, 6, and 9 hours to investigate the effect of Cd on cell growth in a time and dose-dependent manner. Further possible Cd-derived alterations, as the peroxidation of membrane lipids or the functional impairment of the enzymatic antioxidant protection mechanisms, were investigated by determination of the MDA content and via catalase (CAT) activity detection. Moreover, ascorbic acid (AsA) pre-treatment was subjected to investigate the assumed positive effect of AsA against Cd toxicity. We show here on one hand that cells suffer under the influence of Cd, but on the other hand, they substantially profit from AsA supplementation. Because S. pombe is known to shares many molecular, and biochemical similarities with higher organisms, the effect of AsA in cadmium toxicity elimination might be expected to a similar extent also in other cell types.

Current research on circular RNAs and their potential clinical implications in breast cancer

Breast cancer (BC) is one of the most common cancers and the leading causes of death among women worldwide, and its morbidity rate is growing. Discovery of novel biomarkers is necessary for early BC detection, treatment, and prognostication. Circular RNAs (circRNAs), a novel type of endogenous non-coding RNAs with covalently closed continuous loops, have been found to have a crucial role in tumorigenesis. Studies have demonstrated that circRNAs are aberrantly expressed in the tumor tissues and plasma of patients with BC, and they modulate gene expression affecting the proliferation, metastasis, and chemoresistance of BC by specifically binding and regulating the expression of microRNAs (miRNAs). Therefore, circRNAs can be used as novel potential diagnostic and prognostic markers, and therapeutic targets for BC. This article summarizes the properties, functions, and regulatory mechanisms of circRNAs, particularly current research on their association with BC proliferation, metastasis, and chemoresistance.

Saccharomyces cerevisiae RNA lariat debranching enzyme, Dbr1p, is required for completion of reverse transcription by the retrovirus-like element Ty1 and cleaves branched Ty1 RNAs

RNA debranching enzymes are 2'-5' phosphodiesterases found in all eukaryotes. Their main role is cleavage of intron RNA lariat branch points, promoting RNA turnover via exonucleases. Consistent with this role, cells with reduced RNA debranching enzyme activity accumulate intron RNA lariats. The Saccharomyces cerevisiae RNA debranching enzyme Dbr1p is also a host factor for the yeast long terminal repeat (LTR) retrotransposon Ty1, a model for many aspects of retroviral replication. Fittingly, the human RNA debranching enzyme Dbr1 is a host factor for the human immunodeficiency virus, HIV-1. The yeast and human RNA debranching enzymes act at the reverse transcription stages for Ty1 and HIV-1, respectively. Although efficient production of full-length Ty1 cDNA requires Dbr1p, the findings reported here indicate that production of the earliest distinct cDNA product, minus strand strong stop DNA (-sssDNA), is equivalent in wild type and dbr1∆ mutant cells. Several branched Ty1 RNAs are shown to accumulate in dbr1∆ cells during retrotransposition. These data are consistent with creation of Ty1 RNA branches prior to Ty1 reverse transcription and their removal by Dbr1p to allow efficient extension of early cDNA products. The data support the possibility that RNA branch formation and cleavage play broadly shared, but unknown roles in retroviral and LTR retrotransposon reverse transcription.

RNA Lariat Debranching Enzyme as a Retroviral and Long-Terminal-Repeat Retrotransposon Host Factor

Host cell factors are integral to viral replication. Human immunodeficiency virus 1 (HIV-1), the retroviral agent of acquired immune deficiency syndrome, requires several host factors for reverse transcription of the viral genomic RNA (gRNA) into DNA shortly after viral entry. One of these host factors is the RNA lariat debranching enzyme (Dbr1), which cleaves the 2'-5' bond of branched and lariat RNAs. A recent study has revealed that Dbr1 cleaves HIV-1 gRNA lariats that form early after viral entry. Without Dbr1 activity, HIV-1 reverse transcription stalls, consistent with blockage of viral reverse transcriptase at gRNA branch points. These findings echo an earlier study with the long-terminal-repeat retrotransposon of Saccharomyces cerevisiae, Ty1, which is a retrovirus model. Currently, branching and debranching of viral gRNA are not widely recognized as features of HIV-1 replication, and the role of a gRNA lariat is not known. Future studies will determine whether these gRNA dynamics represent fundamental features of retroviral biology and whether they occur for other positive-sense RNA viruses.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.

Functional Analysis of Circular RNAs

Circular RNAs characterize a class of widespread and diverse endogenous RNAs which are non-coding RNAs that are made by back-splicing events and have covalently closed loops with no polyadenylated tails. Various indications specify that circular RNAs (circRNAs) are plentiful in the human transcriptome. However, their participation in biological processes remains mostly undescribed. To date thousands of circRNAs have been revealed in organisms ranging from Drosophila melanogaster to Homo sapiens. Functional studies specify that these transcripts control expression of protein-coding linear transcripts and thus encompass a key component of gene expression regulation. This chapter provide a comprehensive overview on functional validation of circRNAs. Furthermore, we discuss the recent modern methodologies for the functional validation of circRNAs such as RNA interference (RNAi) gene silencing assay, luciferase reporter assays, circRNA gain-of-function investigation via overexpression of circular transcript assay, RT-q-PCR quantification, and other latest applicable assays. The methods described in this chapter are demonstrated on the cellular model.

Bioinformatics and Translation Elongation

Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.

Conservation of mRNA quality control factor Ski7 and its diversification through changes in alternative splicing and gene duplication

Eukaryotes maintain fidelity of gene expression by preferential degradation of aberrant mRNAs that arise by errors in RNA processing reactions. In Saccharomyces cerevisiae, Ski7 plays an important role in this mRNA quality control by mediating mRNA degradation by the RNA exosome. Ski7 was initially thought to be restricted to Saccharomyces cerevisiae and close relatives because the SKI7 gene and its paralog HBS1 arose by whole genome duplication (WGD) in a recent ancestor. We have recently shown that the preduplication gene was alternatively spliced and that Ski7 function predates WGD. Here, we use transcriptome analysis of diverse eukaryotes to show that diverse eukaryotes use alternative splicing of SKI7/HBS1 to encode two proteins. Although alternative splicing affects the same intrinsically disordered region of the protein, the pattern of splice site usage varies. This alternative splicing event arose in an early eukaryote that is a common ancestor of plants, animals, and fungi. Remarkably, through changes in alternative splicing and gene duplication, the Ski7 protein has diversified such that different species express one of four distinct Ski7-like proteins. We also show experimentally that the Saccharomyces cerevisiae SKI7 gene has undergone multiple changes that are incompatible with the Hbs1 function and may also have undergone additional changes to optimize mRNA quality control. The combination of transcriptome analysis in diverse eukaryotes and genetic analysis in yeast clarifies the mechanism by which a Ski7-like protein is expressed across eukaryotes and provides a unique view of changes in alternative splicing patterns of one gene over long evolutionary time.

Transcriptome-wide Interrogation of the Functional Intronome by Spliceosome Profiling

Full understanding of eukaryotic transcriptomes and how they respond to different conditions requires deep knowledge of all sites of intron excision. Although RNA sequencing (RNA-seq) provides much of this information, the low abundance of many spliced transcripts (often due to their rapid cytoplasmic decay) limits the ability of RNA-seq alone to reveal the full repertoire of spliced species. Here, we present "spliceosome profiling," a strategy based on deep sequencing of RNAs co-purifying with late-stage spliceosomes. Spliceosome profiling allows for unambiguous mapping of intron ends to single-nucleotide resolution and branchpoint identification at unprecedented depths. Our data reveal hundreds of new introns in S. pombe and numerous others that were previously misannotated. By providing a means to directly interrogate sites of spliceosome assembly and catalysis genome-wide, spliceosome profiling promises to transform our understanding of RNA processing in the nucleus, much as ribosome profiling has transformed our understanding mRNA translation in the cytoplasm.

Introduction to Fission Yeast as a Model System

Here, we briefly outline the history of fission yeast, its life cycle, and aspects of its biology that make it a useful model organism for studying problems of eukaryotic molecular and cell biology.

Functions for fission yeast splicing factors SpSlu7 and SpPrp18 in alternative splice-site choice and stress-specific regulated splicing

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5'ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5'ss in dtd1+ intron 1 and of an upstream alternative 3'ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5'ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3'ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

ARSDA: A New Approach for Storing, Transmitting and Analyzing Transcriptomic Data

Two major stumbling blocks exist in high-throughput sequencing (HTS) data analysis. The first is the sheer file size, typically in gigabytes when uncompressed, causing problems in storage, transmission, and analysis. However, these files do not need to be so large, and can be reduced without loss of information. Each HTS file, either in compressed .SRA or plain text .fastq format, contains numerous identical reads stored as separate entries. For example, among 44,603,541 forward reads in the SRR4011234.sra file (from a Bacillus subtilis transcriptomic study) deposited at NCBI’s SRA database, one read has 497,027 identical copies. Instead of storing them as separate entries, one can and should store them as a single entry with the SeqID_NumCopy format (which I dub as FASTA+ format). The second is the proper allocation of reads that map equally well to paralogous genes. I illustrate in detail a new method for such allocation. I have developed ARSDA software that implement these new approaches. A number of HTS files for model species are in the process of being processed and deposited at http://coevol.rdc.uottawa.ca to demonstrate that this approach not only saves a huge amount of storage space and transmission bandwidth, but also dramatically reduces time in downstream data analysis. Instead of matching the 497,027 identical reads separately against the B. subtilis genome, one only needs to match it once. ARSDA includes functions to take advantage of HTS data in the new sequence format for downstream data analysis such as gene expression characterization. I contrasted gene expression results between ARSDA and Cufflinks so readers can better appreciate the strength of ARSDA. ARSDA is freely available for Windows, Linux. and Macintosh computers at http://dambe.bio.uottawa.ca/ARSDA/ARSDA.aspx.

Deep learning of the splicing (epi)genetic code reveals a novel candidate mechanism linking histone modifications to ESC fate decision

Alternative splicing (AS) is a genetically and epigenetically regulated pre-mRNA processing to increase transcriptome and proteome diversity. Comprehensively decoding these regulatory mechanisms holds promise in getting deeper insights into a variety of biological contexts involving in AS, such as development and diseases. We assembled splicing (epi)genetic code, DeepCode, for human embryonic stem cell (hESC) differentiation by integrating heterogeneous features of genomic sequences, 16 histone modifications with a multi-label deep neural network. With the advantages of epigenetic features, DeepCode significantly improves the performance in predicting the splicing patterns and their changes during hESC differentiation. Meanwhile, DeepCode reveals the superiority of epigenomic features and their dominant roles in decoding AS patterns, highlighting the necessity of including the epigenetic properties when assembling a more comprehensive splicing code. Moreover, DeepCode allows the robust predictions across cell lineages and datasets. Especially, we identified a putative H3K36me3-regulated AS event leading to a nonsense-mediated mRNA decay of BARD1. Reduced BARD1 expression results in the attenuation of ATM/ATR signalling activities and further the hESC differentiation. These results suggest a novel candidate mechanism linking histone modifications to hESC fate decision. In addition, when trained in different contexts, DeepCode can be expanded to a variety of biological and biomedical fields.

The histone variant H2A.Z promotes splicing of weak introns

In this study, Nissen et al. investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast S. pombe. The findings suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases.

Multiple lines of evidence implicate chromatin in the regulation of premessenger RNA (pre-mRNA) splicing. However, the influence of chromatin factors on cotranscriptional splice site usage remains unclear. Here we investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast Schizosaccharomyces pombe. Using epistatic miniarray profiles (EMAPs) to survey the genetic interaction landscape of the Swr1 nucleosome remodeling complex, which deposits H2A.Z, we uncovered evidence for functional interactions with components of the spliceosome. In support of these genetic connections, splicing-specific microarrays show that H2A.Z and the Swr1 ATPase are required during temperature stress for the efficient splicing of a subset of introns. Notably, affected introns are enriched for H2A.Z occupancy and more likely to contain nonconsensus splice sites. To test the significance of the latter correlation, we mutated the splice sites in an affected intron to consensus and found that this suppressed the requirement for H2A.Z in splicing of that intron. These data suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases. Consistent with this model, we show that overexpression of splicing ATPase Prp16 suppresses both the growth and splicing defects seen in the absence of H2A.Z.

Circular RNAs: analysis, expression and potential functions

Just a few years ago, it had been assumed that the dominant RNA isoforms produced from eukaryotic genes were variants of messenger RNA, functioning as intermediates in gene expression. In early 2012, however, a surprising discovery was made: circular RNA (circRNA) was shown to be a transcriptional product in thousands of human and mouse genes and in hundreds of cases constituted the dominant RNA isoform. Subsequent studies revealed that the expression of circRNAs is developmentally regulated, tissue and cell-type specific, and shared across the eukaryotic tree of life. These features suggest important functions for these molecules. Here, we describe major advances in the field of circRNA biology, focusing on the regulation of and functional roles played by these molecules.

Analysis of RNA Metabolism in Fission Yeast

Here we focus on the biogenesis and function of messenger RNA (mRNA) in fission yeast cells. Following a general introduction that also briefly touches on other classes of RNA, we provide an overview of methods used to analyze mRNAs throughout their life cycles.

Structural toggle in the RNaseH domain of Prp8 helps balance splicing fidelity and catalytic efficiency

Pre-mRNA splicing is an essential step of eukaryotic gene expression that requires both high efficiency and high fidelity. Prp8 has long been considered the "master regulator" of the spliceosome, the molecular machine that executes pre-mRNA splicing. Cross-linking and structural studies place the RNaseH domain (RH) of Prp8 near the spliceosome's catalytic core and demonstrate that prp8 alleles that map to a 17-aa extension in RH stabilize it in one of two mutually exclusive structures, the biological relevance of which are unknown. We performed an extensive characterization of prp8 alleles that map to this extension and, using in vitro and in vivo reporter assays, show they fall into two functional classes associated with the two structures: those that promote error-prone/efficient splicing and those that promote hyperaccurate/inefficient splicing. Identification of global locations of endogenous splice-site activation by lariat sequencing confirms the fidelity effects seen in our reporter assays. Furthermore, we show that error-prone/efficient RH alleles suppress a prp2 mutant deficient at promoting the first catalytic step of splicing, whereas hyperaccurate/inefficient RH alleles exhibit synthetic sickness. Together our data indicate that prp8 RH alleles link splicing fidelity with catalytic efficiency by biasing the relative stabilities of distinct spliceosome conformations. We hypothesize that the spliceosome "toggles" between such error-prone/efficient and hyperaccurate/inefficient conformations during the splicing cycle to regulate splicing fidelity.

Large-scale analysis of branchpoint usage across species and cell lines

The coding sequence of each human pre-mRNA is interrupted, on average, by 11 introns that must be spliced out for proper gene expression. Each intron contains three obligate signals: a 5′ splice site, a branch site, and a 3′ splice site. Splice site usage has been mapped exhaustively across different species, cell types, and cellular states. In contrast, only a small fraction of branch sites have been identified even once. The few reported annotations of branch site are imprecise as reverse transcriptase skips several nucleotides while traversing a 2–5 linkage. Here, we report large-scale mapping of the branchpoints from deep sequencing data in three different species and in the SF3B1 K700E oncogenic mutant background. We have developed a novel method whereby raw lariat reads are refined by U2snRNP/pre-mRNA base-pairing models to return the largest current data set of branchpoint sequences with quality metrics. This analysis discovers novel modes of U2snRNA:pre-mRNA base-pairing conserved in yeast and provides insight into the biogenesis of intron circles. Finally, matching branch site usage with isoform selection across the extensive panel of ENCODE RNA-seq data sets offers insight into the mechanisms by which branchpoint usage drives alternative splicing.

The dynamic landscape of fission yeast meiosis alternative-splice isoforms

Alternative splicing increases the diversity of transcriptomes and proteomes in metazoans. The extent to which alternative splicing is active and functional in unicellular organisms is less understood. Here, we exploit a single-molecule long-read sequencing technique and develop an open-source software program called SpliceHunter to characterize the transcriptome in the meiosis of fission yeast. We reveal 14,353 alternative splicing events in 17,669 novel isoforms at different stages of meiosis, including antisense and read-through transcripts. Intron retention is the major type of alternative splicing, followed by alternate “intron in exon.” Seven hundred seventy novel transcription units are detected; 53 of the predicted proteins show homology in other species and form theoretical stable structures. We report the complexity of alternative splicing along isoforms, including 683 intra-molecularly co-associated intron pairs. We compare the dynamics of novel isoforms based on the number of supporting full-length reads with those of annotated isoforms and explore the translational capacity and quality of novel isoforms. The evaluation of these factors indicates that the majority of novel isoforms are unlikely to be both condition-specific and translatable but consistent with the possibility of biologically functional novel isoforms. Moreover, the co-option of these unusual transcripts into newly born genes seems likely. Together, the results of this study highlight the diversity and dynamics at the isoform level in the sexual development of fission yeast.

Identification of new branch points and unconventional introns in Saccharomyces cerevisiae

Spliced messages constitute one-fourth of expressed mRNAs in the yeast Saccharomyces cerevisiae, and most mRNAs in metazoans. Splicing requires 5' splice site (5'SS), branch point (BP), and 3' splice site (3'SS) elements, but the role of the BP in splicing control is poorly understood because BP identification remains difficult. We developed a high-throughput method, Branch-seq, to map BPs and 5'SSs of isolated RNA lariats. Applied to S. cerevisiae, Branch-seq detected 76% of expressed, annotated BPs and identified a comparable number of novel BPs. We performed RNA-seq to confirm associated 3'SS locations, identifying some 200 novel splice junctions, including an AT-AC intron. We show that several yeast introns use two or even three different BPs, with effects on 3'SS choice, protein coding potential, or RNA stability, and identify novel introns whose splicing changes during meiosis or in response to stress. Together, these findings show unanticipated complexity of splicing in yeast.

The power of fission: yeast as a tool for understanding complex splicing

Pre-mRNA splicing is an essential component of eukaryotic gene expression. Many metazoans, including humans, regulate alternative splicing patterns to generate expansions of their proteome from a limited number of genes. Importantly, a considerable fraction of human disease causing mutations manifest themselves through altering the sequences that shape the splicing patterns of genes. Thus, understanding the mechanistic bases of this complex pathway will be an essential component of combating these diseases. Dating almost to the initial discovery of splicing, researchers have taken advantage of the genetic tractability of budding yeast to identify the components and decipher the mechanisms of splicing. However, budding yeast lacks the complex splicing machinery and alternative splicing patterns most relevant to humans. More recently, many researchers have turned their efforts to study the fission yeast, Schizosaccharomyces pombe, which has retained many features of complex splicing, including degenerate splice site sequences, the usage of exonic splicing enhancers, and SR proteins. Here, we review recent work using fission yeast genetics to examine pre-mRNA splicing, highlighting its promise for modeling the complex splicing seen in higher eukaryotes.

The fission yeast MTREC and EJC orthologs ensure the maturation of meiotic transcripts during meiosis

Meiosis is a highly regulated process by which genetic information is transmitted through sexual reproduction. It encompasses unique mechanisms that do not occur in vegetative cells, producing a distinct, well-regulated meiotic transcriptome. During vegetative growth, many meiotic genes are constitutively transcribed, but most of the resulting mRNAs are rapidly eliminated by the Mmi1-MTREC (Mtl1-Red1 core) complex. While Mmi1-MTREC targets premature meiotic RNAs for degradation by the nuclear 3'-5' exoribonuclease exosome during mitotic growth, its role in meiotic gene expression during meiosis is not known. Here, we report that Red5, an essential MTREC component, interacts with pFal1, an ortholog of eukaryotic translation initiation factor eIF4aIII in the fission yeast Schizosaccharomyces pombe In mammals, together with MAGO (Mnh1), Rnps1, and Y14, elF4AIII (pFal1) forms the core of the exon junction complex (EJC), which is essential for transcriptional surveillance and localization of mature mRNAs. In fission yeast, two EJC orthologs, pFal1 and Mnh1, are functionally connected with MTREC, specifically in the process of meiotic gene expression during meiosis. Although pFal1 interacts with Mnh1, Y14, and Rnps1, its association with Mnh1 is not disrupted upon loss of Y14 or Rnps1. Mutations of Red1, Red5, pFal1, or Mnh1 produce severe meiotic defects; the abundance of meiotic transcripts during meiosis decreases; and mRNA maturation processes such as splicing are impaired. Since studying meiosis in mammalian germline cells is difficult, our findings in fission yeast may help to define the general mechanisms involved in accurate meiotic gene expression in higher eukaryotes.

Circular RNAs: Unexpected outputs of many protein-coding genes

Pre-mRNAs from thousands of eukaryotic genes can be non-canonically spliced to generate circular RNAs, some of which accumulate to higher levels than their associated linear mRNA. Recent work has revealed widespread mechanisms that dictate whether the spliceosome generates a linear or circular RNA. For most genes, circular RNA biogenesis via backsplicing is far less efficient than canonical splicing, but circular RNAs can accumulate due to their long half-lives. Backsplicing is often initiated when complementary sequences from different introns base pair and bring the intervening splice sites close together. This process is further regulated by the combinatorial action of RNA binding proteins, which allow circular RNAs to be expressed in unique patterns. Some genes do not require complementary sequences to generate RNA circles and instead take advantage of exon skipping events. It is still unclear what most mature circular RNAs do, but future investigations into their functions will be facilitated by recently described methods to modulate circular RNA levels.

Interconnections Between RNA-Processing Pathways Revealed by a Sequencing-Based Genetic Screen for Pre-mRNA Splicing Mutants in Fission Yeast

Pre-mRNA splicing is an essential component of eukaryotic gene expression and is highly conserved from unicellular yeasts to humans. Here, we present the development and implementation of a sequencing-based reverse genetic screen designed to identify nonessential genes that impact pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, an organism that shares many of the complex features of splicing in higher eukaryotes. Using a custom-designed barcoding scheme, we simultaneously queried ∼3000 mutant strains for their impact on the splicing efficiency of two endogenous pre-mRNAs. A total of 61 nonessential genes were identified whose deletions resulted in defects in pre-mRNA splicing; enriched among these were factors encoding known or predicted components of the spliceosome. Included among the candidates identified here are genes with well-characterized roles in other RNA-processing pathways, including heterochromatic silencing and 3ʹ end processing. Splicing-sensitive microarrays confirm broad splicing defects for many of these factors, revealing novel functional connections between these pathways.

Sequencing of lariat termini in S. cerevisiae reveals 5' splice sites, branch points, and novel splicing events

Pre-mRNA splicing is a central step in the shaping of the eukaryotic transcriptome and in the regulation of gene expression. Yet, due to a focus on fully processed mRNA, common approaches for defining pre-mRNA splicing genome-wide are suboptimal-especially with respect to defining the branch point sequence, a key cis-element that initiates the chemistry of splicing. Here, we report a complementary intron-centered approach designed to more efficiently, simply, and directly define splicing events genome-wide. Specifically, we developed a method distinguished by deep sequencing of lariat intron termini (LIT-seq). In a test of LIT-seq using the budding yeast Saccharomyces cerevisiae, we not only successfully captured the majority of annotated, expressed splicing events but also uncovered 45 novel splicing events, establishing the sensitivity of LIT-seq. Moreover, our libraries were highly enriched with reads that reported on splice sites; by a simple and direct inspection of sequencing reads, we empirically defined both 5' splice sites and branch sites, as well as their consensus sequences, with nucleotide resolution. Additionally, our study revealed that the 3' termini of lariat introns are subject to nontemplated addition of adenosines, characteristic of signals sensed by 3' to 5' RNA turnover machinery. Collectively, this work defines a novel, genome-wide approach for analyzing splicing with unprecedented depth, specificity, and resolution.

Transcriptome-wide RNA processing kinetics revealed using extremely short 4tU labeling

BACKGROUND

RNA levels detected at steady state are the consequence of multiple dynamic processes within the cell. In addition to synthesis and decay, transcripts undergo processing. Metabolic tagging with a nucleotide analog is one way of determining the relative contributions of synthesis, decay and conversion processes globally.

RESULTS

By improving 4-thiouracil labeling of RNA in Saccharomyces cerevisiae we were able to isolate RNA produced during as little as 1 minute, allowing the detection of nascent pervasive transcription. Nascent RNA labeled for 1.5, 2.5 or 5 minutes was isolated and analyzed by reverse transcriptase-quantitative polymerase chain reaction and RNA sequencing. High kinetic resolution enabled detection and analysis of short-lived non-coding RNAs as well as intron-containing pre-mRNAs in wild-type yeast. From these data we measured the relative stability of pre-mRNA species with different high turnover rates and investigated potential correlations with sequence features.

CONCLUSIONS

Our analysis of non-coding RNAs reveals a highly significant association between non-coding RNA stability, transcript length and predicted secondary structure. Our quantitative analysis of the kinetics of pre-mRNA splicing in yeast reveals that ribosomal protein transcripts are more efficiently spliced if they contain intron secondary structures that are predicted to be less stable. These data, in combination with previous results, indicate that there is an optimal range of stability of intron secondary structures that allows for rapid splicing.

Alternative splicing acting as a bridge in evolution

BACKGROUND

Alternative splicing (AS) regulates diverse cellular and developmental functions through alternative protein structures of different isoforms. Alternative exons dominate AS in vertebrates; however, very little is known about the extent and function of AS in lower eukaryotes. To understand the role of introns in gene evolution, we examined AS from a green algal and five fungal genomes using a novel EST-based gene-modeling algorithm (COMBEST).

METHODS

AS from each genome was classified with COMBEST that maps EST sequences to genomes to build gene models. Various aspects of AS were analyzed through statistical methods. The interplay of intron 3n length, phase, coding property, and intron retention (RI) were examined with Chi-square testing.

RESULTS

With 3 to 834 times EST coverage, we identified up to 73% of AS in intron-containing genes and found preponderance of RI among 11 types of AS. The number of exons, expression level, and maximum intron length correlated with number of AS per gene (NAG), and intron-rich genes suppressed AS. Genes with AS were more ancient, and AS was conserved among fungal genomes. Among stopless introns, non-retained introns (NRI) avoided, but major RI preferred 3n length. In contrast, stop-containing introns showed uniform distribution among 3n, 3n+1, and 3n+2 lengths. We found a clue to the intron phase enigma: it was the coding function of introns involved in AS that dictates the intron phase bias.

CONCLUSIONS

Majority of AS is non-functional, and the extent of AS is suppressed for intron-rich genes. RI through 3n length, stop codon, and phase bias bridges the transition from functionless to functional alternative isoforms.

Cwf16p Associating with the Nineteen Complex Ensures Ordered Exon Joining in Constitutive Pre-mRNA Splicing in Fission Yeast

Exons are ligated in an ordered manner without the skipping of exons in the constitutive splicing of pre-mRNAs with multiple introns. To identify factors ensuring ordered exon joining in constitutive pre-mRNA splicing, we previously screened for exon skipping mutants in Schizosaccharomyces pombe using a reporter plasmid, and characterized three exon skipping mutants named ods1 (ordered splicing 1), ods2, and ods3, the responsible genes of which encode Prp2/U2AF59, U2AF23, and SF1, respectively. They form an SF1-U2AF59-U2AF23 complex involved in recognition of the branch and 3' splice sites in pre-mRNA. In the present study, we identified a fourth ods mutant, ods4, which was isolated in an exon-skipping screen. The ods4+ gene encodes Cwf16p, which interacts with the NineTeen Complex (NTC), a complex thought to be involved in the first catalytic step of the splicing reaction. We isolated two multi-copy suppressors for the ods4-1 mutation, Srp2p, an SR protein essential for pre-mRNA splicing, and Tif213p, a translation initiation factor, in S. pombe. The overexpression of Srp2p suppressed the exon-skipping phenotype of all ods mutants, whereas Tif213p suppressed only ods4-1, which has a mutation in the translational start codon of the cwf16 gene. We also showed that the decrease in the transcriptional elongation rate induced by drug treatment suppressed exon skipping in ods4-1. We propose that Cwf16p/NTC participates in the early recognition of the branch and 3' splice sites and cooperates with the SF1-U2AF59-U2AF23 complex to maintain ordered exon joining.

Widespread alternative and aberrant splicing revealed by lariat sequencing

Alternative splicing is an important and ancient feature of eukaryotic gene structure, the existence of which has likely facilitated eukaryotic proteome expansions. Here, we have used intron lariat sequencing to generate a comprehensive profile of splicing events in Schizosaccharomyces pombe, amongst the simplest organisms that possess mammalian-like splice site degeneracy. We reveal an unprecedented level of alternative splicing, including alternative splice site selection for over half of all annotated introns, hundreds of novel exon-skipping events, and thousands of novel introns. Moreover, the frequency of these events is far higher than previous estimates, with alternative splice sites on average activated at ∼3% the rate of canonical sites. Although a subset of alternative sites are conserved in related species, implying functional potential, the majority are not detectably conserved. Interestingly, the rate of aberrant splicing is inversely related to expression level, with lowly expressed genes more prone to erroneous splicing. Although we validate many events with RNAseq, the proportion of alternative splicing discovered with lariat sequencing is far greater, a difference we attribute to preferential decay of aberrantly spliced transcripts. Together, these data suggest the spliceosome possesses far lower fidelity than previously appreciated, highlighting the potential contributions of alternative splicing in generating novel gene structures.

From Glacier to Sauna: RNA-Seq of the Human Pathogen Black Fungus Exophiala dermatitidis under Varying Temperature Conditions Exhibits Common and Novel Fungal Response

Exophiala dermatitidis (Wangiella dermatitidis) belongs to the group of the so-called black yeasts. Thanks in part to its thick and strongly melanized cell walls, E. dermatitidis is extremely tolerant to various kinds of stress, including extreme pH, temperature and desiccation. E. dermatitidis is also the agent responsible for various severe illnesses in humans, such as pneumonia and keratitis, and might lead to fatal brain infections. Due to its association with the human environment, its poly-extremophilic lifestyle and its pathogenicity in humans, E. dermatitidis has become an important model organism. In this study we present the functional analysis of the transcriptional response of the fungus at 1°C and 45°C, in comparison with that at 37°C, for two different exposition times, i.e. 1 hour and 1 week. At 1°C, E. dermatitidis uses a large repertoire of tools to acclimatize, such as lipid membrane fluidization, trehalose production or cytoskeleton rearrangement, which allows the fungus to remain metabolically active. At 45°C, the fungus drifts into a replicative state and increases the activity of the Golgi apparatus. As a novel finding, our study provides evidence that, apart from the protein coding genes, non-coding RNAs, circular RNAs as well as fusion-transcripts are differentially regulated and that the function of the fusion-transcripts can be related to the corresponding temperature condition. This work establishes that E. dermatitidis adapts to its environment by modulating coding and non-coding gene transcription levels and through the regulation of chimeric and circular RNAs.

Circular RNA biogenesis can proceed through an exon-containing lariat precursor

Pervasive expression of circular RNA is a recently discovered feature of eukaryotic gene expression programs, yet its function remains largely unknown. The presumed biogenesis of these RNAs involves a non-canonical 'backsplicing' event. Recent studies in mammalian cell culture posit that backsplicing is facilitated by inverted repeats flanking the circularized exon(s). Although such sequence elements are common in mammals, they are rare in lower eukaryotes, making current models insufficient to describe circularization. Through systematic splice site mutagenesis and the identification of splicing intermediates, we show that circular RNA in Schizosaccharomyces pombe is generated through an exon-containing lariat precursor. Furthermore, we have performed high-throughput and comprehensive mutagenesis of a circle-forming exon, which enabled us to discover a systematic effect of exon length on RNA circularization. Our results uncover a mechanism for circular RNA biogenesis that may account for circularization in genes that lack noticeable flanking intronic secondary structure.

Deep developmental transcriptome sequencing uncovers numerous new genes and enhances gene annotation in the sponge Amphimedon queenslandica

Background

The demosponge Amphimedon queenslandica is amongst the few early-branching metazoans with an assembled and annotated draft genome, making it an important species in the study of the origin and early evolution of animals. Current gene models in this species are largely based on in silico predictions and low coverage expressed sequence tag (EST) evidence.

Results

Amphimedon queenslandica protein-coding gene models are improved using deep RNA-Seq data from four developmental stages and CEL-Seq data from 82 developmental samples. Over 86% of previously predicted genes are retained in the new gene models, although 24% have additional exons; there is also a marked increase in the total number of annotated 3’ and 5’ untranslated regions (UTRs). Importantly, these new developmental transcriptome data reveal numerous previously unannotated protein-coding genes in the Amphimedon genome, increasing the total gene number by 25%, from 30,060 to 40,122. In general, Amphimedon genes have introns that are markedly smaller than those in other animals and most of the alternatively spliced genes in Amphimedon undergo intron-retention; exon-skipping is the least common mode of alternative splicing. Finally, in addition to canonical polyadenylation signal sequences, Amphimedon genes are enriched in a number of unique AT-rich motifs in their 3’ UTRs.

Conclusions

The inclusion of developmental transcriptome data has substantially improved the structure and composition of protein-coding gene models in Amphimedon queenslandica, providing a more accurate and comprehensive set of genes for functional and comparative studies. These improvements reveal the Amphimedon genome is comprised of a remarkably high number of tightly packed genes. These genes have small introns and there is pervasive intron retention amongst alternatively spliced transcripts. These aspects of the sponge genome are more similar unicellular opisthokont genomes than to other animal genomes.

Electronic supplementary material

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

Widespread exon skipping triggers degradation by nuclear RNA surveillance in fission yeast

Exon skipping is considered a principal mechanism by which eukaryotic cells expand their transcriptome and proteome repertoires, creating different splice variants with distinct cellular functions. Here we analyze RNA-seq data from 116 transcriptomes in fission yeast (Schizosaccharomyces pombe), covering multiple physiological conditions as well as transcriptional and RNA processing mutants. We applied brute-force algorithms to detect all possible exon-skipping events, which were widespread but rare compared to normal splicing events. Exon-skipping events increased in cells deficient for the nuclear exosome or the 5′-3′ exonuclease Dhp1, and also at late stages of meiotic differentiation when nuclear-exosome transcripts decreased. The pervasive exon-skipping transcripts were stochastic, did not increase in specific physiological conditions, and were mostly present at less than one copy per cell, even in the absence of nuclear RNA surveillance and during late meiosis. These exon-skipping transcripts are therefore unlikely to be functional and may reflect splicing errors that are actively removed by nuclear RNA surveillance. The average splicing rate by exon skipping was ∼0.24% in wild type and ∼1.75% in nuclear exonuclease mutants. We also detected approximately 250 circular RNAs derived from single or multiple exons. These circular RNAs were rare and stochastic, although a few became stabilized during quiescence and in splicing mutants. Using an exhaustive search algorithm, we also uncovered thousands of previously unknown splice sites, indicating pervasive splicing; yet most of these splicing variants were cryptic and increased in nuclear degradation mutants. This study highlights widespread but low frequency alternative or aberrant splicing events that are targeted by nuclear RNA surveillance.

Genetic interaction mapping reveals a role for the SWI/SNF nucleosome remodeler in spliceosome activation in fission yeast

Although numerous regulatory connections between pre-mRNA splicing and chromatin have been demonstrated, the precise mechanisms by which chromatin factors influence spliceosome assembly and/or catalysis remain unclear. To probe the genetic network of pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, we constructed an epistatic mini-array profile (E-MAP) and discovered many new connections between chromatin and splicing. Notably, the nucleosome remodeler SWI/SNF had strong genetic interactions with components of the U2 snRNP SF3 complex. Overexpression of SF3 components in ΔSWI/SNF cells led to inefficient splicing of many fission yeast introns, predominantly those with non-consensus splice sites. Deletion of SWI/SNF decreased recruitment of the splicing ATPase Prp2, suggesting that SWI/SNF promotes co-transcriptional spliceosome assembly prior to first step catalysis. Importantly, defects in SWI/SNF as well as SF3 overexpression each altered nucleosome occupancy along intron-containing genes, illustrating that the chromatin landscape both affects—and is affected by—co-transcriptional splicing.

It has recently become apparent that most introns are removed from pre-mRNA while the transcript is still engaged with RNA polymerase II (RNAPII). To gain insight into possible roles for chromatin in co-transcriptional splicing, we generated a genome-wide genetic interaction map in fission yeast and uncovered numerous connections between splicing and chromatin. The SWI/SNF remodeling complex is typically thought to activate gene expression by relieving barriers to polymerase elongation imposed by nucleosomes. Here we show that this remodeler is important for an early step in splicing in which Prp2, an RNA-dependent ATPase, is recruited to the assembling spliceosome to promote catalytic activation. Interestingly, introns with sub-optimal splice sites are particularly dependent on SWI/SNF, suggesting the impact of nucleosome dynamics on the kinetics of spliceosome assembly and catalysis. By monitoring nucleosome occupancy, we show significant alterations in nucleosome density in particular splicing and chromatin mutants, which generally paralleled the levels of RNAPII. Taken together, our findings challenge the notion that nucleosomes simply act as barriers to elongation; rather, we suggest that polymerase pausing at nucleosomes can activate gene expression by allowing more time for co-transcriptional splicing.

Regulation of splicing factors by alternative splicing and NMD is conserved between kingdoms yet evolutionarily flexible

Ultraconserved elements, unusually long regions of perfect sequence identity, are found in genes encoding numerous RNA-binding proteins including arginine-serine rich (SR) splicing factors. Expression of these genes is regulated via alternative splicing of the ultraconserved regions to yield mRNAs that are degraded by nonsense-mediated mRNA decay (NMD), a process termed unproductive splicing (Lareau et al. 2007; Ni et al. 2007). As all human SR genes are affected by alternative splicing and NMD, one might expect this regulation to have originated in an early SR gene and persisted as duplications expanded the SR family. But in fact, unproductive splicing of most human SR genes arose independently (Lareau et al. 2007). This paradox led us to investigate the origin and proliferation of unproductive splicing in SR genes. We demonstrate that unproductive splicing of the splicing factor SRSF5 (SRp40) is conserved among all animals and even observed in fungi; this is a rare example of alternative splicing conserved between kingdoms, yet its effect is to trigger mRNA degradation. As the gene duplicated, the ancient unproductive splicing was lost in paralogs, and distinct unproductive splicing evolved rapidly and repeatedly to take its place. SR genes have consistently employed unproductive splicing, and while it is exceptionally conserved in some of these genes, turnover in specific events among paralogs shows flexible means to the same regulatory end.

Genome-wide discovery of human splicing branchpoints

During the splicing reaction, the 5′ intron end is joined to the branchpoint nucleotide, selecting the next exon to incorporate into the mature RNA and forming an intron lariat, which is excised. Despite a critical role in gene splicing, the locations and features of human splicing branchpoints are largely unknown. We use exoribonuclease digestion and targeted RNA-sequencing to enrich for sequences that traverse the lariat junction and, by split and inverted alignment, reveal the branchpoint. We identify 59,359 high-confidence human branchpoints in >10,000 genes, providing a first map of splicing branchpoints in the human genome. Branchpoints are predominantly adenosine, highly conserved, and closely distributed to the 3′ splice site. Analysis of human branchpoints reveals numerous novel features, including distinct features of branchpoints for alternatively spliced exons and a family of conserved sequence motifs overlapping branchpoints we term B-boxes, which exhibit maximal nucleotide diversity while maintaining interactions with the keto-rich U2 snRNA. Different B-box motifs exhibit divergent usage in vertebrate lineages and associate with other splicing elements and distinct intron–exon architectures, suggesting integration within a broader regulatory splicing code. Lastly, although branchpoints are refractory to common mutational processes and genetic variation, mutations occurring at branchpoint nucleotides are enriched for disease associations.

Widespread use of non-productive alternative splice sites in Saccharomyces cerevisiae

Saccharomyces cerevisiae has been used as a model system to investigate the mechanisms of pre-mRNA splicing but only a few examples of alternative splice site usage have been described in this organism. Using RNA-Seq analysis of nonsense-mediated mRNA decay (NMD) mutant strains, we show that many S. cerevisiae intron-containing genes exhibit usage of alternative splice sites, but many transcripts generated by splicing at these sites are non-functional because they introduce premature termination codons, leading to degradation by NMD. Analysis of splicing mutants combined with NMD inactivation revealed the role of specific splicing factors in governing the use of these alternative splice sites and identified novel functions for Prp17p in enhancing the use of branchpoint-proximal upstream 3′ splice sites and for Prp18p in suppressing the usage of a non-canonical AUG 3′-splice site in GCR1. The use of non-productive alternative splice sites can be increased in stress conditions in a promoter-dependent manner, contributing to the down-regulation of genes during stress. These results show that alternative splicing is frequent in S. cerevisiae but masked by RNA degradation and that the use of alternative splice sites in this organism is mostly aimed at controlling transcript levels rather than increasing proteome diversity.

Accurate gene expression requires the transfer of gene information from DNA to RNA. When DNA is transcribed into RNA, part of the RNA needs to be removed (spliced) to generate a proper copy of the genetic information. This process needs to be very accurate to preserve the genetic information that will be transferred into proteins. Our study shows that in baker's yeast, the splicing process does not always produce the correctly spliced products, as RNA splicing events frequently utilize incorrect splice sites. However, these deficient RNA molecules are eliminated from cells by a quality control mechanism to preserve the integrity of the genetic information. However, incorrect splicing is not useless, as it can be used to regulate the quantity of RNA that is generated.

LaSSO, a strategy for genome-wide mapping of intronic lariats and branch points using RNA-seq

Both canonical and alternative splicing of RNAs are governed by intronic sequence elements and produce transient lariat structures fastened by branch points within introns. To map precisely the location of branch points on a genomic scale, we developed LaSSO (Lariat Sequence Site Origin), a data-driven algorithm which utilizes RNA-seq data. Using fission yeast cells lacking the debranching enzyme Dbr1, LaSSO not only accurately identified canonical splicing events, but also pinpointed novel, but rare, exon-skipping events, which may reflect aberrantly spliced transcripts. Compromised intron turnover perturbed gene regulation at multiple levels, including splicing and protein translation. Notably, Dbr1 function was also critical for the expression of mitochondrial genes and for the processing of self-spliced mitochondrial introns. LaSSO showed better sensitivity and accuracy than algorithms used for computational branch-point prediction or for empirical branch-point determination. Even when applied to a human data set acquired in the presence of debranching activity, LaSSO identified both canonical and exon-skipping branch points. LaSSO thus provides an effective approach for defining high-resolution maps of branch-site sequences and intronic elements on a genomic scale. LaSSO should be useful to validate introns and uncover branch-point sequences in any eukaryote, and it could be integrated into RNA-seq pipelines.