Secondary structure is required for 3' splice site recognition in yeast

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

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.

Impact on splicing in Saccharomyces cerevisiae of random 50-base sequences inserted into an intron

Intron splicing is a key regulatory step in gene expression in eukaryotes. Three sequence elements required for splicing-5' and 3' splice sites and a branchpoint-are especially well-characterized in Saccharomyces cerevisiae, but our understanding of additional intron features that impact splicing in this organism is incomplete, due largely to its small number of introns. To overcome this limitation, we constructed a library in S. cerevisiae of random 50-nt (N50) elements individually inserted into the intron of a reporter gene and quantified canonical splicing and the use of cryptic splice sites by sequencing analysis. More than 70% of approximately 140,000 N50 elements reduced splicing by at least 20%. N50 features, including higher GC content, presence of GU repeats, and stronger predicted secondary structure of its pre-mRNA, correlated with reduced splicing efficiency. A likely basis for the reduced splicing of such a large proportion of variants is the formation of RNA structures that pair N50 bases-such as the GU repeats-with other bases specifically within the reporter pre-mRNA analyzed. However, multiple models were unable to explain more than a small fraction of the variance in splicing efficiency across the library, suggesting that complex nonlinear interactions in RNA structures are not accurately captured by RNA structure prediction methods. Our results imply that the specific context of a pre-mRNA may determine the bases allowable in an intron to prevent secondary structures that reduce splicing. This large data set can serve as a resource for further exploration of splicing mechanisms.

Splicing factor Prp18p promotes genome-wide fidelity of consensus 3'-splice sites

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is challenging considering the low complexity of the 3'SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3'SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3'SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3'SS. The role of Prp18p in 3'SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3'SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3'SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3'SS sequences.

α-tubulin regulation by 5' introns in Saccharomyces cerevisiae

Across eukaryotic genomes, multiple α- and β-tubulin genes require regulation to ensure sufficient production of tubulin heterodimers. Features within these gene families that regulate expression remain underexplored. Here, we investigate the role of the 5' intron in regulating α-tubulin expression in Saccharomyces cerevisiae. We find that the intron in the α-tubulin, TUB1, promotes α-tubulin expression and cell fitness during microtubule stress. The role of the TUB1 intron depends on proximity to the TUB1 promoter and sequence features that are distinct from the intron in the alternative α-tubulin isotype, TUB3. These results lead us to perform a screen to identify genes that act with the TUB1 intron. We identified several genes involved in chromatin remodeling, α/β-tubulin heterodimer assembly, and the spindle assembly checkpoint. We propose a model where the TUB1 intron promotes expression from the chromosomal locus and that this may represent a conserved mechanism for tubulin regulation under conditions that require high levels of tubulin production.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

Spliceosomal mutations decouple 3' splice site fidelity from cellular fitness

The fidelity of splice site selection is thought to be critical for proper gene expression and cellular fitness. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is a daunting task considering the low complexity of the 3'SS consensus sequence YAG. Here we show that inactivating the near-essential splicing factor Prp18p results in a global activation of alternative 3'SS, many of which harbor sequences that highly diverge from the YAG consensus, including some highly unusual non-AG 3'SS. We show that the role of Prp18p in 3'SS fidelity is promoted by physical interactions with the essential splicing factors Slu7p and Prp8p and synergized by the proofreading activity of the Prp22p helicase. Strikingly, structure-guided point mutations that disrupt Prp18p-Slu7p and Prp18p-Prp8p interactions mimic the loss of 3'SS fidelity without any impact on cellular growth, suggesting that accumulation of incorrectly spliced transcripts does not have a major deleterious effect on cellular viability. These results show that spliceosomes exhibit remarkably relaxed fidelity in the absence of Prp18p, and that new 3'SS sampling can be achieved genome-wide without a major negative impact on cellular fitness, a feature that could be used during evolution to explore new productive alternative splice sites.

Splicing of branchpoint-distant exons is promoted by Cactin, Tls1 and the ubiquitin-fold-activated Sde2

Intron diversity facilitates regulated gene expression and alternative splicing. Spliceosomes excise introns after recognizing their splicing signals: the 5'-splice site (5'ss), branchpoint (BP) and 3'-splice site (3'ss). The latter two signals are recognized by U2 small nuclear ribonucleoprotein (snRNP) and its accessory factors (U2AFs), but longer spacings between them result in weaker splicing. Here, we show that excision of introns with a BP-distant 3'ss (e.g. rap1 intron 2) requires the ubiquitin-fold-activated splicing regulator Sde2 in Schizosaccharomyces pombe. By monitoring splicing-specific ura4 reporters in a collection of S. pombe mutants, Cay1 and Tls1 were identified as additional regulators of this process. The role of Sde2, Cay1 and Tls1 was further confirmed by increasing BP-3'ss spacings in a canonical tho5 intron. We also examined BP-distant exons spliced independently of these factors and observed that RNA secondary structures possibly bridged the gap between the two signals. These proteins may guide the 3'ss towards the spliceosome's catalytic centre by folding the RNA between the BP and 3'ss. Orthologues of Sde2, Cay1 and Tls1, although missing in the intron-poor Saccharomyces cerevisiae, are present in intron-rich eukaryotes, including humans. This type of intron-specific pre-mRNA splicing appears to have evolved for regulated gene expression and alternative splicing of key heterochromatin factors.

Improved recognition of splice sites in A. thaliana by incorporating secondary structure information into sequence-derived features: a computational study

Identification of splice sites is an important aspect with regard to the prediction of gene structure. In most of the existing splice site prediction studies, machine learning algorithms coupled with sequence-derived features have been successfully employed for splice site recognition. However, the splice site identification by incorporating the secondary structure information is lacking, particularly in plant species. Thus, we made an attempt in this study to evaluate the performance of structural features on the splice site prediction accuracy in Arabidopsis thaliana. Prediction accuracies were evaluated with the sequence-derived features alone as well as by incorporating the structural features into the sequence-derived features, where support vector machine (SVM) was employed as prediction algorithm. Both short (40 base pairs) and long (105 base pairs) sequence datasets were considered for evaluation. After incorporating the secondary structure features, improvements in accuracies were observed only for the longer sequence dataset and the improvement was found to be higher with the sequence-derived features that accounted nucleotide dependencies. On the other hand, either a little or no improvement in accuracies was found for the short sequence dataset. The performance of SVM was further compared with that of LogitBoost, Random Forest (RF), AdaBoost and XGBoost machine learning methods. The prediction accuracies of SVM, AdaBoost and XGBoost were observed to be at par and higher than that of RF and LogitBoost algorithms. While prediction was performed by taking all the sequence-derived features along with the structural features, a little improvement in accuracies was found as compared to the combination of individual sequence-based features and structural features. To the best of our knowledge, this is the first attempt concerning the computational prediction of splice sites using machine learning methods by incorporating the secondary structure information into the sequence-derived features. All the source codes are available at https://github.com/meher861982/SSFeature.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03036-8.

Reciprocal hemizygosity analysis reveals that the Saccharomyces cerevisiae CGI121 gene affects lag time duration in synthetic grape must

It is standard practice to ferment white wines at low temperatures (10–18°C). However, low temperatures increase fermentation duration and risk of problem ferments, leading to significant costs. The lag duration at fermentation initiation is heavily impacted by temperature; therefore, identification of Saccharomyces cerevisiae genes influencing fermentation kinetics is of interest for winemaking. We selected 28 S. cerevisiae BY4743 single deletants, from a prior list of open reading frames (ORFs) mapped to quantitative trait loci (QTLs) on Chr. VII and XIII, influencing the duration of fermentative lag time. Five BY4743 deletants, Δapt1, Δcgi121, Δclb6, Δrps17a, and Δvma21, differed significantly in their fermentative lag duration compared to BY4743 in synthetic grape must (SGM) at 15 °C, over 72 h. Fermentation at 12.5°C for 528 h confirmed the longer lag times of BY4743 Δcgi121, Δrps17a, and Δvma21. These three candidates ORFs were deleted in S. cerevisiae RM11-1a and S288C to perform single reciprocal hemizygosity analysis (RHA). RHA hybrids and single deletants of RM11-1a and S288C were fermented at 12.5°C in SGM and lag time measurements confirmed that the S288C allele of CGI121 on Chr. XIII, encoding a component of the EKC/KEOPS complex, increased fermentative lag phase duration. Nucleotide sequences of RM11-1a and S288C CGI121 alleles differed by only one synonymous nucleotide, suggesting that intron splicing, codon bias, or positional effects might be responsible for the impact on lag phase duration. This research demonstrates a new role of CGI121 and highlights the applicability of QTL analysis for investigating complex phenotypic traits in yeast.

Mating-Type-Specific Ribosomal Proteins Control Aspects of Sexual Reproduction in Cryptococcus neoformans

The MAT locus of Cryptococcus neoformans has a bipolar organization characterized by an unusually large structure, spanning over 100 kb. MAT genes have been characterized by functional genetics as being involved in sexual reproduction and virulence. However, classical gene replacement failed to achieve mutants for five MAT genes (RPL22, RPO41, MYO2, PRT1, and RPL39), indicating that they are likely essential. In the present study, targeted gene replacement was performed in a diploid strain for both the α and a alleles of the ribosomal genes RPL22 and RPL39 Mendelian analysis of the progeny confirmed that both RPL22 and RPL39 are essential for viability. Ectopic integration of the RPL22 allele of opposite MAT identity in the heterozygous RPL22a/rpl22αΔ or RPL22α/rpl22aΔ mutant strains failed to complement their essential phenotype. Evidence suggests that this is due to differential expression of the RPL22 genes, and an RNAi-dependent mechanism that contributes to control RPL22a expression. Furthermore, via CRISPR/Cas9 technology, the RPL22 alleles were exchanged in haploid MATα and MATa strains of C. neoformans These RPL22 exchange strains displayed morphological and genetic defects during bilateral mating. These results contribute to elucidating functions of C. neoformans essential mating type genes that may constitute a type of imprinting system to promote inheritance of nuclei of both mating types.

Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns’ splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron’s splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein’s domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

An experimental evolution study involving an inefficiently spliced intron reveals that the splicing machinery, introns, and RNA quality control factors evolve in Cis and in Trans when cells optimize their transcriptome to new challenges.

RNAs in the spliceosome: Insight from cryoEM structures

Pre-mRNA splicing is catalyzed by the spliceosome, a multimegadalton RNA-protein complex. The spliceosome undergoes dramatic compositional and conformational changes through the splicing cycle, forming at least 10 distinct complexes. Recent high-resolution cryoEM structures of various spliceosomal complexes revealed unprecedented details of this large molecular machine. This review highlights insight into the structure and function of the spliceosomal RNA components obtained from these new structures, with a focus on the yeast spliceosome. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Introns provide a platform for intergenic regulatory feedback of RPL22 paralogs in yeast

Ribosomal protein genes (RPGs) in Saccharomyces cerevisiae are a remarkable regulatory group that may serve as a model for understanding genetic redundancy in evolutionary adaptations. Most RPGs exist as pairs of highly conserved functional paralogs with divergent untranslated regions and introns. We examined the roles of introns in strains with various combinations of intron and gene deletions in RPL22, RPL2, RPL16, RPL37, RPL17, RPS0, and RPS18 paralog pairs. We found that introns inhibited the expression of their genes in the RPL22 pair, with the RPL22B intron conferring a much stronger effect. While the WT RPL22A/RPL22B mRNA ratio was 93/7, the rpl22aΔi/RPL22B and RPL22A/rpl22bΔi ratios were >99/<1 and 60/40, respectively. The intron in RPL2A stimulated the expression of its own gene, but the removal of the other introns had little effect on expression of the corresponding gene pair. Rpl22 protein abundances corresponded to changes in mRNAs.

Using splicing reporters containing endogenous intron sequences, we demonstrated that these effects were due to the inhibition of splicing by Rpl22 proteins but not by their RNA-binding mutant versions. Indeed, only WT Rpl22A/Rpl22B proteins (but not the mutants) interacted in a yeast three-hybrid system with an RPL22B intronic region between bp 165 and 236. Transcriptome analysis showed that both the total level of Rpl22 and the A/B ratio were important for maintaining the WT phenotype. The data presented here support the contention that the Rpl22B protein has a paralog-specific role.

The RPL22 singleton of Kluyveromyces lactis, which did not undergo whole genome duplication, also responded to Rpl22-mediated inhibition in K. lactis cells. Vice versa, the overproduction of the K. lactis protein reduced the expression of RPL22A/B in S. cerevisiae. The extraribosomal function of of the K. lactis Rpl22 suggests that the loop regulating RPL22 paralogs of S. cerevisiae evolved from autoregulation.

Structure of the yeast spliceosomal postcatalytic P complex

The spliceosome undergoes dramatic changes in a splicing cycle. Structures of B, Bact, C, C*, and intron lariat spliceosome complexes revealed mechanisms of 5'-splice site (ss) recognition, branching, and intron release, but lacked information on 3'-ss recognition, exon ligation, and exon release. Here we report a cryo-electron microscopy structure of the postcatalytic P complex at 3.3-angstrom resolution, revealing that the 3' ss is mainly recognized through non-Watson-Crick base pairing with the 5' ss and branch point. Furthermore, one or more unidentified proteins become stably associated with the P complex, securing the 3' exon and potentially regulating activity of the helicase Prp22. Prp22 binds nucleotides 15 to 21 in the 3' exon, enabling it to pull the intron-exon or ligated exons in a 3' to 5' direction to achieve 3'-ss proofreading or exon release, respectively.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

The effects of structure on pre-mRNA processing and stability

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

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.

Links between mRNA splicing, mRNA quality control, and intellectual disability

In recent years, the impairment of RNA binding proteins that play key roles in the post-transcriptional regulation of gene expression has been linked to numerous neurological diseases. These RNA binding proteins perform critical mRNA processing steps in the nucleus, including splicing, polyadenylation, and export. In many cases, these RNA binding proteins are ubiquitously expressed raising key questions about why only brain function is impaired. Recently, mutations in the ZC3H14 gene, encoding an evolutionarily conserved, polyadenosine RNA binding protein, have been linked to a nonsyndromic form of autosomal recessive intellectual disability. Thus far, research on ZC3H14 and its Nab2 orthologs in budding yeast and Drosophila reveals that ZC3H14/Nab2 is important for mRNA processing and neuronal patterning. Two recent studies now provide evidence that ZC3H14/Nab2 may function in the quality control of mRNA splicing and export and could help to explain the molecular defects that cause neuronal dysfunction and lead to an inherited form of intellectual disability. These studies on ZC3H14/Nab2 reveal new clues to the puzzle of why loss of the ubiquitously expressed ZC3H14 protein specifically affects neurons.

RNA structure in splicing: An evolutionary perspective

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic β sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

Parallel computation of genome-scale RNA secondary structure to detect structural constraints on human genome

Background

RNA secondary structure around splice sites is known to assist normal splicing by promoting spliceosome recognition. However, analyzing the structural properties of entire intronic regions or pre-mRNA sequences has been difficult hitherto, owing to serious experimental and computational limitations, such as low read coverage and numerical problems.

Results

Our novel software, “ParasoR”, is designed to run on a computer cluster and enables the exact computation of various structural features of long RNA sequences under the constraint of maximal base-pairing distance. ParasoR divides dynamic programming (DP) matrices into smaller pieces, such that each piece can be computed by a separate computer node without losing the connectivity information between the pieces. ParasoR directly computes the ratios of DP variables to avoid the reduction of numerical precision caused by the cancellation of a large number of Boltzmann factors. The structural preferences of mRNAs computed by ParasoR shows a high concordance with those determined by high-throughput sequencing analyses.

Using ParasoR, we investigated the global structural preferences of transcribed regions in the human genome. A genome-wide folding simulation indicated that transcribed regions are significantly more structural than intergenic regions after removing repeat sequences and k-mer frequency bias. In particular, we observed a highly significant preference for base pairing over entire intronic regions as compared to their antisense sequences, as well as to intergenic regions. A comparison between pre-mRNAs and mRNAs showed that coding regions become more accessible after splicing, indicating constraints for translational efficiency. Such changes are correlated with gene expression levels, as well as GC content, and are enriched among genes associated with cytoskeleton and kinase functions.

Conclusions

We have shown that ParasoR is very useful for analyzing the structural properties of long RNA sequences such as mRNAs, pre-mRNAs, and long non-coding RNAs whose lengths can be more than a million bases in the human genome. In our analyses, transcribed regions including introns are indicated to be subject to various types of structural constraints that cannot be explained from simple sequence composition biases. ParasoR is freely available at https://github.com/carushi/ParasoR.

Electronic supplementary material

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

Splicing of Nascent RNA Coincides with Intron Exit from RNA Polymerase II

Protein-coding genes in eukaryotes are transcribed by RNA polymerase II (Pol II) and introns are removed from pre-mRNA by the spliceosome. Understanding the time lag between Pol II progression and splicing could provide mechanistic insights into the regulation of gene expression. Here, we present two single-molecule nascent RNA sequencing methods that directly determine the progress of splicing catalysis as a function of Pol II position. Endogenous genes were analyzed on a global scale in budding yeast. We show that splicing is 50% complete when Pol II is only 45 nt downstream of introns, with the first spliced products observed as introns emerge from Pol II. Perturbations that slow the rate of spliceosome assembly or speed up the rate of transcription caused splicing delays, showing that regulation of both processes determines in vivo splicing profiles. We propose that matched rates streamline the gene expression pathway, while allowing regulation through kinetic competition.

IntSplice: prediction of the splicing consequences of intronic single-nucleotide variations in the human genome

Precise spatiotemporal regulation of splicing is mediated by splicing cis-elements on pre-mRNA. Single-nucleotide variations (SNVs) affecting intronic cis-elements possibly compromise splicing, but no efficient tool has been available to identify them. Following an effect-size analysis of each intronic nucleotide on annotated alternative splicing, we extracted 105 parameters that could affect the strength of the splicing signals. However, we could not generate reliable support vector regression models to predict the percent-splice-in (PSI) scores for normal human tissues. Next, we generated support vector machine (SVM) models using 110 parameters to directly differentiate pathogenic SNVs in the Human Gene Mutation Database and normal SNVs in the dbSNP database, and we obtained models with a sensitivity of 0.800±0.041 (mean and s.d.) and a specificity of 0.849±0.021. Our IntSplice models were more discriminating than SVM models that we generated with Shapiro-Senapathy score and MaxEntScan::score3ss. We applied IntSplice to a naturally occurring and nine artificial intronic mutations in RAPSN causing congenital myasthenic syndrome. IntSplice correctly predicted the splicing consequences for nine of the ten mutants. We created a web service program, IntSplice (http://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice) to predict splicing-affecting SNVs at intronic positions from -50 to -3.

An intronic RNA structure modulates expression of the mRNA biogenesis factor Sus1

Sus1 is a conserved protein involved in chromatin remodeling and mRNA biogenesis. Unlike most yeast genes, the SUS1 pre-mRNA of Saccharomyces cerevisiae contains two introns and is alternatively spliced, retaining one or both introns in response to changes in environmental conditions. SUS1 splicing may allow the cell to control Sus1 expression, but the mechanisms that regulate this process remain unknown. Using in silico analyses together with NMR spectroscopy, gel electrophoresis, and UV thermal denaturation experiments, we show that the downstream intron (I2) of SUS1 forms a weakly stable, 37-nucleotide stem-loop structure containing the branch site near its apical loop and the 3' splice site after the stem terminus. A cellular assay revealed that two of four mutants containing altered I2 structures had significantly impaired SUS1 expression. Semiquantitative RT-PCR experiments indicated that all mutants accumulated unspliced SUS1 pre-mRNA and/or induced distorted levels of fully spliced mRNA relative to wild type. Concomitantly, Sus1 cellular functions in histone H2B deubiquitination and mRNA export were affected in I2 hairpin mutants that inhibited splicing. This work demonstrates that I2 structure is relevant for SUS1 expression, and that this effect is likely exerted through modulation of splicing.

Diverse regulation of 3' splice site usage

The regulation of splice site (SS) usage is important for alternative pre-mRNA splicing and thus proper expression of protein isoforms in cells; its disruption causes diseases. In recent years, an increasing number of novel regulatory elements have been found within or nearby the 3'SS in mammalian genes. The diverse elements recruit a repertoire of trans-acting factors or form secondary structures to regulate 3'SS usage, mostly at the early steps of spliceosome assembly. Their mechanisms of action mainly include: (1) competition between the factors for RNA elements, (2) steric hindrance between the factors, (3) direct interaction between the factors, (4) competition between two splice sites, or (5) local RNA secondary structures or longer range loops, according to the mode of protein/RNA interactions. Beyond the 3'SS, chromatin remodeling/transcription, posttranslational modifications of trans-acting factors and upstream signaling provide further layers of regulation. Evolutionarily, some of the 3'SS elements seem to have emerged in mammalian ancestors. Moreover, other possibilities of regulation such as that by non-coding RNA remain to be explored. It is thus likely that there are more diverse elements/factors and mechanisms that influence the choice of an intron end. The diverse regulation likely contributes to a more complex but refined transcriptome and proteome in mammals.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Nucleotide sequence composition adjacent to intronic splice sites improves splicing efficiency via its effect on pre-mRNA local folding in fungi

RNA splicing is the central process of intron removal in eukaryotes known to regulate various cellular functions such as growth, development, and response to external signals. The canonical sequences indicating the splicing sites needed for intronic boundary recognition are well known. However, the roles and evolution of the local folding of intronic and exonic sequence features adjacent to splice sites has yet to be thoroughly studied. Here, focusing on four fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, and Candida albicans), we performed for the first time a comprehensive high-resolution study aimed at characterizing the encoding of intronic splicing efficiency in pre-mRNA transcripts and its effect on intron evolution. Our analysis supports the conjecture that pre-mRNA local folding strength at intronic boundaries is under selective pressure, as it significantly affects splicing efficiency. Specifically, we show that in the immediate region of 12-30 nucleotides (nt) surrounding the intronic donor site there is a preference for weak pre-mRNA folding; similarly, in the region of 15-33 nt surrounding the acceptor and branch sites there is a preference for weak pre-mRNA folding. We also show that in most cases there is a preference for strong pre-mRNA folding further away from intronic splice sites. In addition, we demonstrate that these signals are not associated with gene-specific functions, and they correlate with splicing efficiency measurements (r = 0.77, P = 2.98 × 10(-21)) and with expression levels of the corresponding genes (P = 1.24 × 10(-19)). We suggest that pre-mRNA folding strength in the above-mentioned regions has a direct effect on splicing efficiency by improving the recognition of intronic boundaries. These new discoveries are contributory steps toward a broader understanding of splicing regulation and intronic/transcript evolution.

Cwc21p promotes the second step conformation of the spliceosome and modulates 3' splice site selection

Pre-mRNA splicing involves two transesterification steps catalyzed by the spliceosome. How RNA substrates are positioned in each step and the molecular rearrangements involved, remain obscure. Here, we show that mutations in PRP16, PRP8, SNU114 and the U5 snRNA that affect this process interact genetically with CWC21, that encodes the yeast orthologue of the human SR protein, SRm300/SRRM2. Our microarray analysis shows changes in 3′ splice site selection at elevated temperature in a subset of introns in cwc21Δ cells. Considering all the available data, we propose a role for Cwc21p positioning the 3′ splice site at the transition to the second step conformation of the spliceosome, mediated through its interactions with the U5 snRNP. This suggests a mechanism whereby SRm300/SRRM2, might influence splice site selection in human cells.

Approaches to link RNA secondary structures with splicing regulation

In higher eukaryotes, alternative splicing is usually regulated by protein factors, which bind to the pre-mRNA and affect the recognition of splicing signals. There is recent evidence that the secondary structure of the pre-mRNA may also play an important role in this process, either by facilitating or hindering the interaction with factors and small nuclear ribonucleoproteins (snRNPs) that regulate splicing. Moreover, the secondary structure could play a fundamental role in the splicing of yeast species, which lack many of the regulatory splicing factors present in metazoans. This chapter describes the steps in the analysis of the secondary structure of the pre-mRNA and its possible relation to splicing. As a working example, we use the case of yeast and the problem of the recognition of the 3' splice site (3'ss).

Diversity in mRNA expression of the serine-type carboxypeptidase ocpG in Aspergillus oryzae through intron retention

Alternative splicing is thought to be a means for diversification of products by mRNA modification. Although some intron retentions are predicted by transcriptome analysis in Aspergillus oryzae, its physiological significance remains unknown. We found that intron retention occurred occasionally in the serine-type carboxypeptidase gene, ocpG. Analysis under various culture conditions revealed that extracellular nitrogen conditions influence splicing patterns; this suggested that there might be a correlation between splicing efficiency and the necessity of OcpG activity for obtaining a nitrogen source. Since further analysis showed that splicing occurred independently in each intron, we constructed ocpG intron-exchanging strain by interchanging the positions of intron-1 and intron-2. The splicing pattern indicated the probability that ocpG intron retention was affected by the secondary structures of intronic mRNA.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

Accurate, model-based tuning of synthetic gene expression using introns in S. cerevisiae

Introns are key regulators of eukaryotic gene expression and present a potentially powerful tool for the design of synthetic eukaryotic gene expression systems. However, intronic control over gene expression is governed by a multitude of complex, incompletely understood, regulatory mechanisms. Despite this lack of detailed mechanistic understanding, here we show how a relatively simple model enables accurate and predictable tuning of synthetic gene expression system in yeast using several predictive intron features such as transcript folding and sequence motifs. Using only natural Saccharomyces cerevisiae introns as regulators, we demonstrate fine and accurate control over gene expression spanning a 100 fold expression range. These results broaden the engineering toolbox of synthetic gene expression systems and provide a framework in which precise and robust tuning of gene expression is accomplished.

Synthetic biology is gradually expanding our capability to engineer biology through rational genetic engineering of synthetic gene expression systems. These developments are already paving the way for the accelerated study of biology and applying engineered biological systems to major environmental and health problems. However, our capacity to intelligently modify and control gene expression depends on our ability to apply a broad range of genetic regulators in the engineering process. Here we show that Introns, pivotal regulators of Eukaryotic gene expression, can be rationally engineered to control a synthetic gene expression system of a Eukaryote. We developed a unique reporter-based system to evaluate the effects of engineering splicing in synthetic biology and show that the entire intron repertoire of S. cerevisiae can be accurately used to rationally engineer gene expression. Our results provide both a proof-of-concept for the integration of splicing into synthetic biology designs and a model that can be used by the scientific community for integrating splicing into their own designs. Following the extensive use of transcriptional (promoter) and translational (UTR) elements in synthetic constructs, our results introduce a new major regulatory system, splicing, that can be used to rationally engineer genetic systems.

Quantifying variances in comparative RNA secondary structure prediction

Background

With the advancement of next-generation sequencing and transcriptomics technologies, regulatory effects involving RNA, in particular RNA structural changes are being detected. These results often rely on RNA secondary structure predictions. However, current approaches to RNA secondary structure modelling produce predictions with a high variance in predictive accuracy, and we have little quantifiable knowledge about the reasons for these variances.

Results

In this paper we explore a number of factors which can contribute to poor RNA secondary structure prediction quality. We establish a quantified relationship between alignment quality and loss of accuracy. Furthermore, we define two new measures to quantify uncertainty in alignment-based structure predictions. One of the measures improves on the “reliability score” reported by PPfold, and considers alignment uncertainty as well as base-pair probabilities. The other measure considers the information entropy for SCFGs over a space of input alignments.

Conclusions

Our predictive accuracy improves on the PPfold reliability score. We can successfully characterize many of the underlying reasons for and variances in poor prediction. However, there is still variability unaccounted for, which we therefore suggest comes from the RNA secondary structure predictive model itself.

Brr2p-mediated conformational rearrangements in the spliceosome during activation and substrate repositioning

Brr2p is one of eight RNA helicases involved in pre-mRNA splicing. Detailed understanding of the functions of Brr2p and other spliceosomal helicases has been limited by lack of knowledge of their in vivo substrates. To address this, sites of direct Brr2p-RNA interaction were identified by in vivo UV cross-linking in budding yeast. Cross-links identified in the U4 and U6 small nuclear RNAs (snRNAs) suggest U4/U6 stem I as a Brr2p substrate during spliceosome activation. Further Brr2p cross-links were identified in loop 1 of the U5 snRNA and near splice sites and 3' ends of introns, suggesting the possibility of a previously uncharacterized function for Brr2p in the catalytic center of the spliceosome. Consistent with this, mutant brr2-G858R reduced second-step splicing efficiency and enhanced cross-linking to 3' ends of introns. Furthermore, RNA sequencing indicated preferential inhibition of splicing of introns with structured 3' ends. The Brr2-G858Rp cross-linking pattern in U6 was consistent with an open conformation for the catalytic center of the spliceosome during first-to-second-step transition. We propose a previously unsuspected function for Brr2p in driving conformational rearrangements that lead to competence for the second step of splicing.

Intronic features that determine the selection of the 3' splice site

Most eukaryotic primary transcripts include segments, or introns, that will be accurately removed during RNA biogenesis. This process, known as pre-messenger RNA splicing, is catalyzed by the spliceosome, accurately selecting a set of intronic marks from others apparently equivalent. This identification is critical, as incorrectly spliced RNAs can be toxic for the organism. One of these marks, the dinucleotide AG, signals the intronic 3' end, or 3' splice site (ss). In this review we will focus on those intronic features that have an impact on 3' ss selection. These include the location and type of neighboring sequences, and their distance to the 3' end. We will see that their interplay is needed to select the right intronic end, and that this can be modulated by additional intronic elements that contribute to alternative splicing, whereby diverse RNAs can be generated from identical precursors. This complexity, still poorly understood, is fundamental for the accuracy of gene expression. In addition, a clear knowledge of 3' ss selection is needed to fully decipher the coding potential of genomes.

RNA secondary structure mediates alternative 3'ss selection in Saccharomyces cerevisiae

Alternative splicing is the mechanism by which different combinations of exons in the pre-mRNA give rise to distinct mature mRNAs. This process is mediated by splicing factors that bind the pre-mRNA and affect the recognition of its splicing signals. Saccharomyces species lack many of the regulatory factors present in metazoans. Accordingly, it is generally assumed that the amount of alternative splicing is limited. However, there is recent compelling evidence that yeast have functional alternative splicing, mainly in response to environmental conditions. We have previously shown that sequence and structure properties of the pre-mRNA could explain the selection of 3' splice sites (ss) in Saccharomyces cerevisiae. In this work, we extend our previous observations to build a computational classifier that explains most of the annotated 3'ss in the CDS and 5' UTR of this organism. Moreover, we show that the same rules can explain the selection of alternative 3'ss. Experimental validation of a number of predicted alternative 3'ss shows that their usage is low compared to annotated 3'ss. The majority of these alternative 3'ss introduce premature termination codons (PTCs), suggesting a role in expression regulation. Furthermore, a genome-wide analysis of the effect of temperature, followed by experimental validation, yields only a small number of changes, indicating that this type of regulation is not widespread. Our results are consistent with the presence of alternative 3'ss selection in yeast mediated by the pre-mRNA structure, which can be responsive to external cues, like temperature, and is possibly related to the control of gene expression.

Regulated pre-mRNA splicing: the ghostwriter of the eukaryotic genome

Intron removal is at the heart of mRNA synthesis. It is mediated by one of the cell's largest complexes, the spliceosome. Yet, the fundamental chemistry involved is simple. In this review we will address how the spliceosome acts in diverse ways to optimize gene expression in order to meet the cell's needs. This is done largely by regulating the splicing of key transcripts encoding products that control gene expression pathways. This widespread role is evident even in the yeast Saccharomyces cerevisiae, where many introns appear to have been lost; yet how this control is being achieved is known only in a few cases. Here we explore the relevant examples and posit hypotheses whereby regulated splicing fine-tunes gene expression pathways to maintain cell homeostasis. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.