Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

The first intron of EIJ1 confers a specific response to wounding and herbivore stresses

Plants are constantly exposed to different kinds of biotic stress, such as herbivore attack and wounding. To deal with these stresses, plants have evolved sophisticated defence mechanisms to protect themselves. Previously, we found that EIJ1 (EDS1-interacting J protein 1) plays a negative regulatory role in plant disease resistance in Arabidopsis thaliana. Follow-up studies revealed that EIJ1 specifically responds to wounding and herbivore stresses. The expression of EIJ1 was specifically induced by wounding or herbivore stress, as demonstrated by similar results in EIJ1 protein assay. Interestingly, GUS staining found that the promoter of EIJ1 is not involved in the induction of expression under wounding stress. Instead, we identified the first intron of EIJ1 as a key factor in response to wounding stress. Deleting the first intron of EIJ1 resulted in a loss of response to wounding stress in plants. Our results broaden the role of EIJ1 in plant resistance to biotic stress and provide new insights into plant responses to biotic stress.

Broad variation in response of individual introns to splicing inhibitors in a humanized yeast strain

Intron branchpoint (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect the binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeast Saccharomyces cerevisiae by changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after the addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during cotranscriptional splicing in Plad-B using single-molecule intron tracking to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between the binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten the characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

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.

Broad variation in response of individual introns to splicing inhibitors in a humanized yeast strain

Intron branch point (BP) recognition by the U2 snRNP is a critical step of splicing, vulnerable to recurrent cancer mutations and bacterial natural product inhibitors. The BP binds a conserved pocket in the SF3B1 (human) or Hsh155 (yeast) U2 snRNP protein. Amino acids that line this pocket affect binding of splicing inhibitors like Pladienolide-B (Plad-B), such that organisms differ in their sensitivity. To study the mechanism of splicing inhibitor action in a simplified system, we modified the naturally Plad-B resistant yeast Saccharomyces cerevisiae by changing 14 amino acids in the Hsh155 BP pocket to those from human. This humanized yeast grows normally, and splicing is largely unaffected by the mutation. Splicing is inhibited within minutes after addition of Plad-B, and different introns appear inhibited to different extents. Intron-specific inhibition differences are also observed during co-transcriptional splicing in Plad-B using single-molecule intron tracking (SMIT) to minimize gene-specific transcription and decay rates that cloud estimates of inhibition by standard RNA-seq. Comparison of Plad-B intron sensitivities to those of the structurally distinct inhibitor Thailanstatin-A reveals intron-specific differences in sensitivity to different compounds. This work exposes a complex relationship between binding of different members of this class of inhibitors to the spliceosome and intron-specific rates of BP recognition and catalysis. Introns with variant BP sequences seem particularly sensitive, echoing observations from mammalian cells, where monitoring individual introns is complicated by multi-intron gene architecture and alternative splicing. The compact yeast system may hasten characterization of splicing inhibitors, accelerating improvements in selectivity and therapeutic efficacy.

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation

RNA-binding proteins (RBPs) bind at different positions of the pre-mRNA molecules to promote or reduce the usage of a particular exon. Seeking to understand the working principle of these positional effects, we develop a capture RIC-seq (CRIC-seq) method to enrich specific RBP-associated in situ proximal RNA-RNA fragments for deep sequencing. We determine hnRNPA1-, SRSF1-, and PTBP1-associated proximal RNA-RNA contacts and regulatory mechanisms in HeLa cells. Unexpectedly, the 3D RNA map analysis shows that PTBP1-associated loops in individual introns preferentially promote cassette exon splicing by accelerating asymmetric intron removal, whereas the loops spanning across cassette exon primarily repress splicing. These "positional rules" can faithfully predict PTBP1-regulated splicing outcomes. We further demonstrate that cancer-related splicing quantitative trait loci can disrupt RNA loops by reducing PTBP1 binding on pre-mRNAs to cause aberrant splicing in tumors. Our study presents a powerful method for exploring the functions of RBP-associated RNA-RNA proximal contacts in gene regulation and disease.

Secondary structures in RNA synthesis, splicing and translation

Even though the functional role of mRNA molecules is primarily decided by the nucleotide sequence, several properties are determined by secondary structure conformations. Examples of secondary structures include long range interactions, hairpins, R-loops and G-quadruplexes and they are formed through interactions of non-adjacent nucleotides. Here, we discuss advances in our understanding of how secondary structures can impact RNA synthesis, splicing, translation and mRNA half-life. During RNA synthesis, secondary structures determine RNA polymerase II (RNAPII) speed, thereby influencing splicing. Splicing is also determined by RNA binding proteins and their binding rates are modulated by secondary structures. For the initiation of translation, secondary structures can control the choice of translation start site. Here, we highlight the mechanisms by which secondary structures modulate these processes, discuss advances in technologies to detect and study them systematically, and consider the roles of RNA secondary structures in disease.

Emerging roles of RNA-RNA interactions in transcriptional regulation

Pervasive transcription of the human genome generates a massive amount of noncoding RNAs (ncRNAs) that lack protein-coding potential but play crucial roles in development, differentiation, and tumorigenesis. To achieve these biological functions, ncRNAs must first fold into intricate structures via intramolecular RNA-RNA interactions (RRIs) and then interact with different RNA substrates via intermolecular RRIs. RRIs are usually facilitated, stabilized, or mediated by RNA-binding proteins. With this guiding principle, several protein-based high-throughput methods have been developed for unbiased mapping of defined or all RNA-binding protein-mediated RRIs in various species and cell lines. In addition, some chemical-based approaches are also powerful to detect RRIs globally based on the fact that RNA duplex can be cross-linked by psoralen or its derivative 4'-aminomethyltrioxsalen. These efforts have significantly expanded our understanding of RRIs in determining the specificity and variability of gene regulation. Here, we review the current knowledge of the regulatory roles of RRI, focusing on their emerging roles in transcriptional regulation and nuclear body formation. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Intronic RNA: Ad'junk' mediator of post-transcriptional gene regulation

RNA splicing, the process through which intervening segments of noncoding RNA (introns) are excised from pre-mRNAs to allow for the formation of a mature mRNA product, has long been appreciated for its capacity to add complexity to eukaryotic proteomes. However, evidence suggests that the utility of this process extends beyond protein output and provides cells with a dynamic tool for gene regulation. In this review, we aim to highlight the role that intronic RNA plays in mediating specific splicing outcomes in pre-mRNA processing, as well as explore an emerging class of stable intronic sequences that have been observed to act in gene expression control. Building from underlying flexibility in both sequence and structure, intronic RNA provides mechanisms for post-transcriptional gene regulation that are amenable to the tissue and condition specific needs of eukaryotic cells. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

ARID5B Genetic Polymorphisms Contribute to the Susceptibility and Prognosis of Male Acute Promyelocytic Leukemia

This study was conducted using TagSNPs to systematically explore the relationship between ARID5B polymorphisms and the occurrence, clinical characterization, and prognosis of acute myeloid leukemia (AML). A total of 569 unrelated AML patients and 410 healthy individuals from West China were recruited, and ARID5B TagSNPs were genotyped using iMLDR^® (improved multiplex ligation detection reaction). It was found that the association of ARID5B polymorphisms with AML was most significant in acute promyelocytic leukemia (APL), and exclusively in males, the mutant alleles of rs6415872, rs2393726, rs7073837, rs10821936, and rs7089424 were found to increase the risk of developing APL in men, the odds ratio (OR) were 1.36, 1.74, 1.45, 1.53, and 1.56 (all p < 0.05), respectively. Haplotype analysis revealed that haplotype [AACCG] increased the risk of male APL with an OR of 1.53 (95% confidence interval: 1.10-2.14, p = 0.012). Besides, there was a strong positive additive interaction between rs6415872 and rs10821936, rs7089424, respectively, and cases attributed to the interaction of rs6415872, rs10821936, and rs7089424 accounted for 100%. Furthermore, ARID5B single nucleotide polymorphisms were found associated with clinical features of AML, and rs6415872 was shown to be an independent prognosis factor in APL patients. Besides, dual luciferase report assay showed that rs6415872 may affect the binding activity of PPARG with ARID5B. ARID5B polymorphisms contribute to male APL risk, clinical feature, and prognosis, suggesting the importance of ARDI5B in AML pathogenesis and development, and the gender and subtype preference may prompt some specific mechanisms of ARID5B. Besides, ARID5B polymorphisms might be a potential prognosis biomarker.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

Role and convergent evolution of competing RNA secondary structures in mutually exclusive splicing

Exon or cassette duplication is an important means of expanding protein and functional diversity through mutually exclusive splicing. However, the mechanistic basis of this process in non-arthropod species remains poorly understood. Here, we demonstrate that MRP1 genes underwent tandem exon duplication in Nematoda, Platyhelminthes, Annelida, Mollusca, Arthropoda, Echinodermata, and early-diverging Chordata but not in late-diverging vertebrates. Interestingly, these events were of independent origin in different phyla, suggesting convergent evolution of alternative splicing. Furthermore, we showed that multiple sets of clade-conserved RNA pairings evolved to guide species-specific mutually exclusive splicing in Arthropoda. Importantly, we also identified a similar structural code in MRP exon clusters of the annelid, Capitella teleta, and chordate, Branchiostoma belcheri, suggesting an evolutionarily conserved competing pairing-guided mechanism in bilaterians. Taken together, these data reveal the molecular determinants and RNA pairing-guided evolution of species-specific mutually exclusive splicing spanning more than 600 million years of bilaterian evolution. These findings have a significant impact on our understanding of the evolution of and mechanism underpinning isoform diversity and complex gene structure.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

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.

Introns regulate the production of ribosomal proteins by modulating splicing of duplicated ribosomal protein genes

Most budding yeast introns exist in the many duplicated ribosomal protein genes (RPGs) and it has been posited that they remain there to modulate the expression of RPGs and cell growth in response to stress. However, the mechanism by which introns regulate the expression of RPGs and their impact on the synthesis of ribosomal proteins remain unclear. In this study, we show that introns determine the ratio of ribosomal protein isoforms through asymmetric paralog-specific regulation of splicing. Exchanging the introns and 3′ untranslated regions of the duplicated RPS9 genes altered the splicing efficiency and changed the ratio of the ribosomal protein isoforms. Mutational analysis of the RPS9 genes indicated that splicing is regulated by variations in the intron structure and the 3′ untranslated region. Together these data suggest that preferential splicing of duplicated RPGs provides a means for adjusting the ratio of different ribosomal protein isoforms, while maintaining the overall expression level of each ribosomal protein.

RNA structure replaces the need for U2AF2 in splicing

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

Long-range RNA pairings contribute to mutually exclusive splicing

Mutually exclusive splicing is an important means of increasing the protein repertoire, by which the Down's syndrome cell adhesion molecule (Dscam) gene potentially generates 38,016 different isoforms in Drosophila melanogaster. However, the regulatory mechanisms remain obscure due to the complexity of the Dscam exon cluster. Here, we reveal a molecular model for the regulation of the mutually exclusive splicing of the serpent pre-mRNA based on competition between upstream and downstream RNA pairings. Such dual RNA pairings confer fine tuning of the inclusion of alternative exons. Moreover, we demonstrate that the splicing outcome of alternative exons is mediated in relative pairing strength-correlated mode. Combined comparative genomics analysis and experimental evidence revealed similar bidirectional structural architectures in exon clusters 4 and 9 of the Dscam gene. Our findings provide a novel mechanistic framework for the regulation of mutually exclusive splicing and may offer potentially applicable insights into long-range RNA-RNA interactions in gene regulatory networks.

The association of ADORA2A and ADORA2B polymorphisms with the risk and severity of chronic heart failure: a case-control study of a northern Chinese population

The causes of chronic heart failure (CHF) and its progression are likely to be due to complex genetic factors. Adenosine receptors A2A and A2B (ADORA2A and ADORA2B, respectively) play an important role in cardio-protection. Therefore, polymorphisms in the genes encoding those receptors may affect the risk and severity of CHF. This study was a case-control comparative investigation of 300 northern Chinese Han CHF patients and 400 ethnicity-matched healthy controls. Four common single-nucleotide polymorphisms (SNPs) of ADORA2A (rs2236625, rs2236624, rs4822489, and rs5751876) and one SNP of ADORA2B (rs7208480) were genotyped and an association between SNPs and clinical outcomes was evaluated. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the association. The rs4822489 was significantly associated with the severity of CHF after adjustment for traditional cardiovascular risk factors (p = 0.040, OR = 1.912, 95% CI = 1.029–3.550). However, the five SNPs as well as the haplotypes were not found to be associated with CHF susceptibility. The findings of this study suggest that rs4822489 may contribute to the severity of CHF in the northern Chinese. However, further studies performed in larger populations and aimed at better defining the role of this gene are required.

Effects of CYP3A4 polymorphisms on the plasma concentration of voriconazole

Voriconazole is frequently utilized for the prevention and treatment of invasive fungal infections (IFIs), and is extensively metabolized by the cytochrome P450 (CYP) system. The impact of activity of the genes encoding CYP3A4, CYP3A5, and CYP2C9 on the pharmacokinetics of voriconazole cannot be ignored because, second to CYP2C19, they are the most important enzymes involved in voriconazole metabolism. The influence of genetic polymorphisms in CYP3A4, CYP3A5, and CYP2C9 on the plasma concentrations of voriconazole was evaluated in the present study. The study cohort comprised 158 patients with IFIs in whom 22 single-nucleotide polymorphisms (SNPs) in CYP3A4, CYP3A5, and CYP2C9 were genotyped using the Sequenom MassARRAY RS1000 system, and voriconazole plasma concentrations were measured by high-performance liquid chromatography (HPLC). 40, 91, and 27 patients presented with low (<1 mg/L), normal (1-4 mg/L), and high (>4 mg/L) plasma voriconazole concentrations, respectively. Correlation analysis between polymorphisms and the plasma voriconazole concentration revealed an association between the presence of the rs4646437 T allele and a higher plasma voriconazole concentration [p = 0.033, odds ratio (OR) = 2.832, 95% confidence interval (CI) = 1.086-7.384]. This study has identified a new SNP related to the metabolism of voriconazole, potentially providing novel insight into the influence of CYP3A4 on the pharmacokinetics of this antifungal agent.

detectIR: a novel program for detecting perfect and imperfect inverted repeats using complex numbers and vector calculation

Inverted repeats are present in abundance in both prokaryotic and eukaryotic genomes and can form DNA secondary structures--hairpins and cruciforms that are involved in many important biological processes. Bioinformatics tools for efficient and accurate detection of inverted repeats are desirable, because existing tools are often less accurate and time consuming, sometimes incapable of dealing with genome-scale input data. Here, we present a MATLAB-based program called detectIR for the perfect and imperfect inverted repeat detection that utilizes complex numbers and vector calculation and allows genome-scale data inputs. A novel algorithm is adopted in detectIR to convert the conventional sequence string comparison in inverted repeat detection into vector calculation of complex numbers, allowing non-complementary pairs (mismatches) in the pairing stem and a non-palindromic spacer (loop or gaps) in the middle of inverted repeats. Compared with existing popular tools, our program performs with significantly higher accuracy and efficiency. Using genome sequence data from HIV-1, Arabidopsis thaliana, Homo sapiens and Zea mays for comparison, detectIR can find lots of inverted repeats missed by existing tools whose outputs often contain many invalid cases. detectIR is open source and its source code is freely available at: https://sourceforge.net/projects/detectir.

The adenosine deaminase gene polymorphism is associated with chronic heart failure risk in Chinese

Adenosine (Ado) is an important cardioprotective agent. Since endogenous Ado levels are affected by the enzyme Ado deaminase (ADA), polymorphisms within the ADA gene may exert some effect on chronic heart failure (CHF). This study applied a case-control investigation to 300 northern Chinese Han CHF patients and 400 ethnicity-matched healthy controls in which nine single-nucleotide polymorphisms (SNPs) of ADA were genotyped and association analyses were performed. Odds ratios (ORs) with 95% confidence intervals (CI) were used to assess the association. Overall, rs452159 polymorphism in ADA gene was significantly associated with susceptibility to CHF under the dominant model (p = 0.013, OR = 1.537, 95% CI = 1.10–2.16), after adjustment for age, sex, and traditional cardiovascular risk factors. No difference in genotype distribution and allele frequency for the rs452159 according to the functional New York Heart Association class was found. Furthermore, the values of left ventricular ejection fraction, left-ventricle end-diastolic diameter or left-ventricle end-systolic diameter did not differ significantly among the different rs452159 genotype CHF patients. Although further studies with larger cohorts and other ethnicities are required to validate the conclusions, the findings of this study potentially provide novel insight into the pathogenesis of CHF.

Competition between pre-mRNAs for the splicing machinery drives global regulation of splicing

During meiosis in yeast, global splicing efficiency increases and then decreases. Here we provide evidence that splicing improves due to reduced competition for the splicing machinery. The timing of this regulation corresponds to repression and reactivation of ribosomal protein genes (RPGs) during meiosis. In vegetative cells, RPG repression by rapamycin treatment also increases splicing efficiency. Downregulation of the RPG-dedicated transcription factor gene IFH1 genetically suppresses two spliceosome mutations, prp11-1 and prp4-1, and globally restores splicing efficiency in prp4-1 cells. We conclude that the splicing apparatus is limiting and that pre-messenger RNAs compete. Splicing efficiency of a pre-mRNA therefore depends not just on its own concentration and affinity for limiting splicing factor(s), but also on those of competing pre-mRNAs. Competition between RNAs for limiting processing factors appears to be a general condition in eukaryotes for a variety of posttranscriptional control mechanisms including microRNA (miRNA) repression, polyadenylation, and splicing.

Eukaryotic TPP riboswitch regulation of alternative splicing involving long-distance base pairing

Thiamin pyrophosphate (TPP) riboswitches are found in organisms from all three domains of life. Examples in bacteria commonly repress gene expression by terminating transcription or by blocking ribosome binding, whereas most eukaryotic TPP riboswitches are predicted to regulate gene expression by modulating RNA splicing. Given the widespread distribution of eukaryotic TPP riboswitches and the diversity of their locations in precursor messenger RNAs (pre-mRNAs), we sought to examine the mechanism of alternative splicing regulation by a fungal TPP riboswitch from Neurospora crassa, which is mostly located in a large intron separating protein-coding exons. Our data reveal that this riboswitch uses a long-distance (∼530-nt separation) base-pairing interaction to regulate alternative splicing. Specifically, a portion of the TPP-binding aptamer can form a base-paired structure with a conserved sequence element (α) located near a 5′ splice site, which greatly increases use of this 5′ splice site and promotes gene expression. Comparative sequence analyses indicate that many fungal species carry a TPP riboswitch with similar intron architecture, and therefore the homologous genes in these fungi are likely to use the same mechanism. Our findings expand the scope of genetic control mechanisms relying on long-range RNA interactions to include riboswitches.

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.

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

Higher order RNA structures can mask splicing signals, loop out exons, or constitute riboswitches all of which contributes to the complexity of splicing regulation. We identified a G to A substitution between branch point (BP) and 3′ splice site (3′ss) of Saccharomyces cerevisiae COF1 intron, which dramatically impaired its splicing. RNA structure prediction and in-line probing showed that this mutation disrupted a stem in the BP-3′ss region. Analyses of various COF1 intron modifications revealed that the secondary structure brought about the reduction of BP to 3′ss distance and masked potential 3′ss. We demonstrated the same structural requisite for the splicing of UBC13 intron. Moreover, RNAfold predicted stable structures for almost all distant BP introns in S. cerevisiae and for selected examples in several other Saccharomycotina species. The employment of intramolecular structure to localize 3′ss for the second splicing step suggests the existence of pre-mRNA structure-based mechanism of 3′ss recognition.

RNA secondary structure in mutually exclusive splicing

Mutually exclusive splicing is a regulated means to generate protein diversity, but the underlying mechanisms are poorly understood. Here comparative genome analysis revealed the built-in intronic elements for controlling mutually exclusive splicing of the 14-3-3ξ pre-mRNA. These elements are clade specific but are evolutionarily conserved at the secondary structure level. Combined evidence revealed the triple functions of these inter-intronic RNA pairings in synergistically ensuring the selection of only one of multiple exons, through activation of the proximal variable exon outside the loop by the approximation of cis elements, and simultaneous repression of the exon within the loop, in combination with the physical competition of RNA pairing. Additionally, under this model, we also deciphered a similar structural code in exon clusters 4 and 9 of Dscam (38,016 isoforms) and Mhc (480 isoforms). Our findings suggest a broadly applicable mechanism to ensure mutually exclusive splicing.

Selective constraints in conserved folded RNAs of drosophilid and hominid genomes

Small noncoding RNAs as well as folded RNA structures in genic regions are crucial for many cellular processes. They are involved in posttranscriptional gene regulation (microRNAs), RNA modification (small nucleolar RNAs), regulation of splicing, correct localization of proteins, and many other processes. In most cases, a distinct secondary structure of the molecule is necessary for its correct function. Hence, selection should act to retain the structure of the molecule, although the underlying sequence is allowed to vary. Here, we present the first genome-wide estimates of selective constraints in folded RNA molecules in the nuclear genomes of drosophilids and hominids. In comparison to putatively neutrally evolving sites, we observe substantially reduced rates of substitutions at paired and unpaired sites of folded molecules. We estimated evolutionary constraints to be in the ranges of (0.974,0.991) and (0.895,1.000) for paired nucleotides in drosophilids and hominids, respectively. These values are significantly higher than for constraints at nonsynonymous sites of protein-coding genes in both genera. Nonetheless, valleys of only moderately reduced fitness (s ≈ 10(-4)) are sufficient to generate the observed fraction of nucleotide changes that are removed by purifying selection. In addition, a comparison of selective coefficients between drosophilids and hominids revealed significantly higher constraints in drosophilids, which can be attributed to the difference in long-term effective population size between these two groups of species. This difference is particularly apparent at the independently evolving (unpaired) sites.

Role of RNA structure in regulating pre-mRNA splicing

Pre-mRNA splicing involves removing non-coding introns from RNA transcripts. It is carried out by the spliceosome, along with other auxiliary factors. In general, research in splicing has focused on the sequences within the pre-mRNA, without considering the structures that these sequences might form. We propose that the role of RNA structure deserves more consideration when thinking about splicing mechanisms. RNA structures can inhibit or aid binding of spliceosomal components to the pre-mRNA, or can increase splicing efficiency by bringing important sequences into close proximity. Recent reports have identified proteins and small molecules that can regulate splicing by modulating RNA structures, thereby expanding our knowledge of the mechanisms used to regulate splicing.

Epistatic interactions modulate the evolution of mammalian mitochondrial respiratory complex components

BACKGROUND

The deleterious effect of a mutation can be reverted by a second-site interacting residue. This is an epistatic compensatory process explaining why mutations that are deleterious in some species are tolerated in phylogenetically related lineages, rendering evident that those mutations are, by all means, only deleterious in the species-specific context. Although an extensive and refined theoretical framework on compensatory evolution does exist, the supporting evidence remains limited, especially for protein models. In this current study, we focused on the molecular mechanism underlying the epistatic compensatory process in mammalian mitochondrial OXPHOS proteins using a combination of in-depth structural and sequence analyses.

RESULTS

Modeled human structures were used in this study to predict the structural impairment and recovery of deleterious mutations alone and combined with an interacting compensatory partner, respectively. In two cases, COI and COIII, intramolecular interactions between spatially linked residues restore the folding pattern impaired by the deleterious mutation. In a third case, intermolecular contact between mitochondrial CYB and nuclear CYT1 encoded components of the cytochrome bc1 complex are likely to restore protein binding. Moreover, we observed different modes of compensatory evolution that have resulted in either a quasi-simultaneous occurrence of a mutation and corresponding compensatory partner, or in independent occurrences of mutations in distinct lineages that were always preceded by the compensatory site.

CONCLUSION

Epistatic interactions between individual replacements involving deleterious mutations seems to follow a parsimonious model of evolution in which genomes hold pre-compensating states that subsequently tolerate deleterious mutations. This phenomenon is likely to have been constraining the variability at coevolving sites and shaping the interaction between the mitochondrial and the nuclear genome.

Intronic Alus influence alternative splicing

Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called Alu. A large fraction of Alus are found in intronic sequences, implying extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these intronic Alus on splicing of the flanking exons is largely unknown. Here, we show that more Alus flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that Alu insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two Alu elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic Alus in influencing the splicing of flanking exons, further emphasizing the role of Alus in shaping of the human transcriptome.

Conserved secondary structures in Aspergillus

Background

Recent evidence suggests that the number and variety of functional RNAs (ncRNAs as well as cis-acting RNA elements within mRNAs ) is much higher than previously thought; thus, the ability to computationally predict and analyze RNAs has taken on new importance. We have computationally studied the secondary structures in an alignment of six Aspergillus genomes. Little is known about the RNAs present in this set of fungi, and this diverse set of genomes has an optimal level of sequence conservation for observing the correlated evolution of base-pairs seen in RNAs.

Methodology/Principal Findings

We report the results of a whole-genome search for evolutionarily conserved secondary structures, as well as the results of clustering these predicted secondary structures by structural similarity. We find a total of 7450 predicted secondary structures, including a new predicted ∼60 bp long hairpin motif found primarily inside introns. We find no evidence for microRNAs. Different types of genomic regions are over-represented in different classes of predicted secondary structures. Exons contain the longest motifs (primarily long, branched hairpins), 5′ UTRs primarily contain groupings of short hairpins located near the start codon, and 3′ UTRs contain very little secondary structure compared to other regions. There is a large concentration of short hairpins just inside the boundaries of exons. The density of predicted intronic RNAs increases with the length of introns, and the density of predicted secondary structures within mRNA coding regions increases with the number of introns in a gene.

Conclusions/Sigificance

There are many conserved, high-confidence RNAs of unknown function in these Aspergillus genomes, as well as interesting spatial distributions of predicted secondary structures. This study increases our knowledge of secondary structure in these aspergillus organisms.

Correlation between the secondary structure of pre-mRNA introns and the efficiency of splicing in Saccharomyces cerevisiae

BACKGROUND

Secondary structure interactions within introns have been shown to be essential for efficient splicing of several yeast genes. The nature of these base-pairing interactions and their effect on splicing efficiency were most extensively studied in ribosomal protein gene RPS17B (previously known as RP51B). It was determined that complementary pairing between two sequence segments located downstream of the 5' splice site and upstream of the branchpoint sequence promotes efficient splicing of the RPS17B pre-mRNA, presumably by shortening the branchpoint distance. However, no attempts were made to compute a shortened, 'structural' branchpoint distance and thus the functional relationship between this distance and the splicing efficiency remains unknown.

RESULTS

In this paper we use computational RNA secondary structure prediction to analyze the secondary structure of the RPS17B intron. We show that it is necessary to consider suboptimal structure predictions and to compute the structural branchpoint distances in order to explain previously published splicing efficiency results. Our study reveals that there is a tight correlation between this distance and splicing efficiency levels of intron mutants described in the literature. We experimentally test this correlation on additional RPS17B mutants and intron mutants within two other yeast genes.

CONCLUSION

The proposed model of secondary structure requirements for efficient splicing is the first attempt to specify the functional relationship between pre-mRNA secondary structure and splicing. Our findings provide further insights into the role of pre-mRNA secondary structure in gene splicing in yeast and also offer basis for improvement of computational methods for splice site identification and gene-finding.

Cross-kingdom patterns of alternative splicing and splice recognition

BACKGROUND

Variations in transcript splicing can reveal how eukaryotes recognize intronic splice sites. Retained introns (RIs) commonly appear when the intron definition (ID) mechanism of splice site recognition inconsistently identifies intron-exon boundaries, and cassette exons (CEs) are often caused by variable recognition of splice junctions by the exon definition (ED) mechanism. We have performed a comprehensive survey of alternative splicing across 42 eukaryotes to gain insight into how spliceosomal introns are recognized.

RESULTS

All eukaryotes we studied exhibit RIs, which appear more frequently than previously thought. CEs are also present in all kingdoms and most of the organisms in our analysis. We observe that the ratio of CEs to RIs varies substantially among kingdoms, while the ratio of competing 3' acceptor and competing 5' donor sites remains nearly constant. In addition, we find the ratio of CEs to RIs in each organism correlates with the length of its introns. In all 14 fungi we examined, as well as in most of the 9 protists, RIs far outnumber CEs. This differs from the trend seen in 13 multicellular animals, where CEs occur much more frequently than RIs. The six plants we analyzed exhibit intermediate proportions of CEs and RIs.

CONCLUSION

Our results suggest that most extant eukaryotes are capable of recognizing splice sites via both ID and ED, although ED is most common in multicellular animals and ID predominates in fungi and most protists.

Ligand-induced sequestering of branchpoint sequence allows conditional control of splicing

Background

Despite tremendous progress in understanding the mechanisms of constitutive and alternative splicing, an important and widespread step along the gene expression pathway, our ability to deliberately regulate gene expression at this step remains rudimentary. The present study was performed to investigate whether a theophylline-dependent "splice switch" that sequesters the branchpoint sequence (BPS) within RNA-theophylline complex can regulate alternative splicing.

Results

We constructed a series of pre-mRNAs in which the BPS was inserted within theophylline aptamer. We show that theophylline-induced sequestering of BPS inhibits pre-mRNA splicing both in vitro and in vivo in a dose-dependent manner. Several lines of evidence suggest that theophylline-dependent inhibition of splicing is highly specific, and thermodynamic stability of RNA-theophylline complex as well as the location of BPS within this complex affects the efficiency of splicing inhibition. Finally, we have constructed an alternative splicing model pre-mRNA substrate in which theophylline caused exon skipping both in vitro and in vivo, suggesting that a small molecule-RNA interaction can modulate alternative splicing.

Conclusion

These findings provide the ability to control splicing pattern at will and should have important implications for basic, biotechnological, and biomedical research.

RNA structure is a key regulatory element in pathological ATM and CFTR pseudoexon inclusion events

Genomic variations deep in the intronic regions of pre-mRNA molecules are increasingly reported to affect splicing events. However, there is no general explanation why apparently similar variations may have either no effect on splicing or cause significant splicing alterations. In this work we have examined the structural architecture of pseudoexons previously described in ATM and CFTR patients. The ATM case derives from the deletion of a repressor element and is characterized by an aberrant 5′ss selection despite the presence of better alternatives. The CFTR pseudoexon instead derives from the creation of a new 5′ss that is used while a nearby pre-existing donor-like sequence is never selected. Our results indicate that RNA structure is a major splicing regulatory factor in both cases. Furthermore, manipulation of the original RNA structures can lead to pseudoexon inclusion following the exposure of unused 5′ss already present in their wild-type intronic sequences and prevented to be recognized because of their location in RNA stem structures. Our data show that intrinsic structural features of introns must be taken into account to understand the mechanism of pseudoexon activation in genetic diseases. Our observations may help to improve diagnostics prediction programmes and eventual therapeutic targeting.

The position of yeast snoRNA-coding regions within host introns is essential for their biosynthesis and for efficient splicing of the host pre-mRNA

Genomic location of sequences encoding small nucleolar RNAs (snoRNAs) is peculiar in all eukaryotes from yeast to mammals: most of them are encoded within the introns of host genes. In Saccharomyces cerevisiae, seven snoRNAs show this location. In this work we demonstrate that the position of snoRNA-coding regions with respect to splicing consensus sequences is critical: yeast strains expressing mutant constructs containing shorter or longer spacers (the regions between snoRNA ends and intron splice sites) show a drop in accumulation of U24 and U18 snoRNAs. Further mutational analysis demonstrates that altering the distance between the 3' end of the snoRNA and the branch point is the most important constraint for snoRNA biosynthesis, and that stable external stems, which are sometimes present in introns containing snoRNAs, can overcome the positional effect. Surprisingly enough, splicing of the host introns is clearly affected in most of these constructs indicating that, at least in S. cerevisiae, an incorrect location of snoRNA-coding sequences within the host intron is detrimental to the splicing process. This is different with respect to what was demonstrated in mammals, where the activity of the splicing machinery seems to be dominant with respect to the assembly of snoRNPs, and it is not affected by the location of snoRNA sequences. We also show that intronic box C/D snoRNA recognition and assembly of snoRNPs occur during transcription when splicing sequences are recognized.

An extended RNA binding site for the yeast branch point-binding protein and the role of its zinc knuckle domains in RNA binding

The highly conserved branch point sequence (BPS) of UACUAAC in Saccharomyces cerevisiae is initially recognized by the branch point-binding protein (BBP). Using systematic evolution of ligands by exponential enrichment we have determined that yeast BBP binds the branch point sequence UACUAAC with highest affinity and prefers an additional adenosine downstream of the BPS. Furthermore, we also found that a stem-loop upstream of the BPS enhances binding both to an artificially designed RNA (30-fold effect) and to an RNA from a yeast intron (3-fold effect). The zinc knuckles of BBP are partially responsible for the enhanced binding to the stem-loop but do not appear to have a significant role in the binding of BBP to single-strand RNA substrates. C-terminal deletions of BBP reveal that the linker regions between the two zinc knuckles and between the N-terminal RNA binding domains (KH and QUA2 domains) and the first zinc knuckle are important for binding to RNA. The lack of involvement of the second highly conserved zinc knuckle in RNA binding suggests that this zinc knuckle plays a different role in RNA processing than enhancing the binding of BBP to the BPS.

Intronic binding sites for hnRNP A/B and hnRNP F/H proteins stimulate pre-mRNA splicing

hnRNP A/B proteins modulate the alternative splicing of several mammalian and viral pre-mRNAs, and are typically viewed as proteins that enforce the activity of splicing silencers. Here we show that intronic hnRNP A/B–binding sites (ABS) can stimulate the in vitro splicing of pre-mRNAs containing artificially enlarged introns. Stimulation of in vitro splicing could also be obtained by providing intronic ABS in trans through the use of antisense oligonucleotides containing a non-hybridizing ABS-carrying tail. ABS-tailed oligonucleotides also improved the in vivo inclusion of an alternative exon flanked by an enlarged intron. Notably, binding sites for hnRNP F/H proteins (FBS) replicate the activity of ABS by improving the splicing of an enlarged intron and by modulating 5′ splice-site selection. One hypothesis formulated to explain these effects is that bound hnRNP proteins self-interact to bring in closer proximity the external pair of splice sites. Consistent with this model, positioning FBS or ABS at both ends of an intron was required to stimulate splicing of some pre-mRNAs. In addition, a computational analysis of the configuration of putative FBS and ABS located at the ends of introns supports the view that these motifs have evolved to support cooperative interactions. Our results document a positive role for the hnRNP A/B and hnRNP F/H proteins in generic splicing, and suggest that these proteins may modulate the conformation of mammalian pre-mRNAs.

Typically viewed as enforcing splicing silencers, hnRNP A/B proteins may facilitate splicing by modulating the conformation of mammalian pre-mRNAs.

The iStem, a long-range RNA secondary structure element required for efficient exon inclusion in the Drosophila Dscam pre-mRNA

The Drosophila Dscam gene encodes 38,016 different proteins, due to alternative splicing of 95 of its 115 exons, that function in axon guidance and innate immunity. The alternative exons are organized into four clusters, and the exons within each cluster are spliced in a mutually exclusive manner. Here we describe an evolutionarily conserved RNA secondary structure we call the Inclusion Stem (iStem) that is required for efficient inclusion of all 12 variable exons in the exon 4 cluster. Although the iStem governs inclusion or exclusion of the entire exon 4 cluster, it does not play a significant role in determining which variable exon is selected. Thus, the iStem is a novel type of regulatory element that simultaneously controls the splicing of multiple alternative exons.

Identification of splicing silencers and enhancers in sense Alus: a role for pseudoacceptors in splice site repression

Auxiliary splicing signals in introns play an important role in splice site selection, but these elements are poorly understood. We show that a subset of serine/arginine (SR)-rich proteins activate a cryptic 3' splice site in a sense Alu repeat located in intron 4 of the human LST1 gene. Utilization of this cryptic splice site is controlled by juxtaposed Alu-derived splicing silencers and enhancers between closely linked short tandem repeats TNFd and TNFe. Systematic mutagenesis of these elements showed that AG dinucleotides that were not preceded by purine residues were critical for repressing exon inclusion of a chimeric splicing reporter. Since the splice acceptor-like sequences are present in excess in exonic splicing silencers, these signals may contribute to inhibition of a large number of pseudosites in primate genomes.

Exonization of AluYa5 in the human ACE gene requires mutations in both 3' and 5' splice sites and is facilitated by a conserved splicing enhancer

Ancient Alu elements have been shown to be included in mature transcripts by point mutations that improve their 5′ or 3′ splice sites. We have examined requirements for exonization of a younger, disease-associated AluYa5 in intron 16 of the human ACE gene. A single G>C transversion in position −3 of the new Alu exon was insufficient for Alu exonization and a significant inclusion in mRNA was only observed when improving several potential splice donor sites in the presence of 3′ CAG. Since complete Alu exonization was not achieved by optimizing traditional splicing signals, including the branch site, we tested whether auxiliary elements in AluYa5 were required for constitutive inclusion. Exonization was promoted by a SELEX-predicted heptamer in Alu consensus sequence 222–228 and point mutations in highly conserved nucleotides of this heptamer decreased Alu inclusion. In addition, we show that Alu exonization was facilitated by a subset of serine/arginine-rich (SR) proteins through activation of the optimized 3′ splice site. Finally, the haplotype- and allele-specific ACE minigenes generated similar splicing patterns in both ACE-expressing and non-expressing cells, suggesting that previously reported allelic association with plasma ACE activity and cardiovascular disease is not attributable to differential splicing of introns 16 and 17.

Ultraconserved elements in insect genomes: a highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Rare autosomal recessive cardiac valvular form of Ehlers-Danlos syndrome results from mutations in the COL1A2 gene that activate the nonsense-mediated RNA decay pathway

Splice site mutations in the COL1A2 gene of type I collagen can give rise to forms of Ehlers-Danlos syndrome (EDS) because of partial or complete skipping of exon 6, as well as to mild, moderate, or lethal forms of osteogenesis imperfecta as a consequence of skipping of other exons. We identified three unrelated individuals with a rare recessively inherited form of EDS (characterized by joint hypermobility, skin hyperextensibility, and cardiac valvular defects); in two of them, COL1A2 messenger RNA (mRNA) instability results from compound heterozygosity for splice site mutations in the COL1A2 gene, and, in the third, it results from homozygosity for a nonsense codon. The splice site mutations led to use of cryptic splice donor sites, creation of a downstream premature termination codon, and extremely unstable mRNA. In the wild-type allele, the two introns (IVS11 and IVS24) in which these mutations occurred were usually spliced slowly in relation to their respective immediate upstream introns. In the mutant alleles, the upstream intron was removed, so that exon skipping could not occur. In the context of the mutation in IVS24, computer-generated folding of a short stretch of mRNA surrounding the mutation site demonstrated realignment of the relationships between the donor and acceptor sites that could facilitate use of a cryptic donor site. These findings suggest that the order of intron removal is an important variable in prediction of mutation outcome at splice sites and that folding of the nascent mRNA could be one element that contributes to determination of order of splicing. The complete absence of pro alpha 2(I) chains has the surprising effect of producing cardiac valvular disease without bone involvement.

A stem structure in fibroblast growth factor receptor 2 transcripts mediates cell-type-specific splicing by approximating intronic control elements

Alternative splicing of fibroblast growth factor receptor 2 (FGFR2) occurs in a cell-type-specific manner with the mutually exclusive use of exon IIIb or exon IIIc. Specific inclusion of exon IIIb is observed in epithelial cells, whereas exon IIIc inclusion is seen in mesenchymal cells. Epithelium-specific activation of exon IIIb and repression of exon IIIc are coordinately regulated by intronic activating sequence 2 (IAS2) and intronic splicing activator and repressor (ISAR) elements in FGFR2 pre-mRNA. Previously, it has been suggested that IAS2 and a 20-nucleotide core sequence of ISAR form a stem structure that allows for the proper regulation of FGFR2 alternative splicing. Replacement of IAS2 and the ISAR core with random sequences capable of stem formation resulted in the proper activation of exon IIIb and repression of exon IIIc in epithelial cells. Given the high degree of phylogenetic conservation of the IAS2-ISAR core structure and the fact that unrelated stem-forming sequences could functionally substitute for IAS2 and ISAR elements, we postulated that the stem structure facilitated the approximation of intronic control elements. Indeed, deletion of the entire stem-loop region and juxtaposition of sequences immediately upstream of IAS2 with sequences immediately downstream of the ISAR core maintained proper cell-type-specific inclusion of exon IIIb. These data demonstrate that IAS2 and the ISAR core are dispensable for the cell-type-specific activation of exon IIIb; thus, the major, if not the sole, role of the IAS2-ISAR stem in exon IIIb activation is to approximate sequences upstream of IAS2 with sequences downstream of the ISAR core. The downstream sequence is very likely a highly conserved GCAUG element, which we show was required for efficient exon IIIb activation.

Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene

Evidence for the evolutionary maintenance of a hairpin structure possibly involved in intron processing had been found in intron 1 of the alcohol dehydrogenase gene (Adh) in diverse Drosophila species. In this study, the putative hairpin structure was evaluated systematically in Drosophila melanogaster by elimination of either side of the stem using site-directed mutagenesis. The effects of these mutations and the compensatory double mutant on intron splicing efficiency and ADH protein production were assayed in Drosophila melanogaster Schneider L2 cells and germ-line transformed adult flies. Mutations that disrupt the putative hairpin structure right upstream of the intron branch point were found to cause a significant reduction in both splicing efficiency and ADH protein production. In contrast, the compensatory double mutant that restores the putative hairpin structure was indistinguishable from the WT in both splicing efficiency and ADH level. It was also observed by mutational analysis that a more stable secondary structure (with a longer stem) in this intron decreases both splicing efficiency and ADH protein production. Implications for RNA secondary structure and intron evolution are discussed.