Branch nucleophile selection in pre-mRNA splicing: evidence for the bulged duplex model

Splicing efficiency of minor introns in a mouse model of SMA predominantly depends on their branchpoint sequence and can involve the contribution of major spliceosome components

Spinal muscular atrophy (SMA) is a devastating neurodegenerative disease caused by reduced amounts of the ubiquitously expressed Survival of Motor Neuron (SMN) protein. In agreement with its crucial role in the biogenesis of spliceosomal snRNPs, SMN-deficiency is correlated to numerous splicing alterations in patient cells and various tissues of SMA mouse models. Among the snRNPs whose assembly is impacted by SMN-deficiency, those involved in the minor spliceosome are particularly affected. Importantly, splicing of several, but not all U12-dependent introns has been shown to be affected in different SMA models. Here, we have investigated the molecular determinants of this differential splicing in spinal cords from SMA mice. We show that the branchpoint sequence (BPS) is a key element controlling splicing efficiency of minor introns. Unexpectedly, splicing of several minor introns with suboptimal BPS is not affected in SMA mice. Using in vitro splicing experiments and oligonucleotides targeting minor or major snRNAs, we show for the first time that splicing of these introns involves both the minor and major machineries. Our results strongly suggest that splicing of a subset of minor introns is not affected in SMA mice because components of the major spliceosome compensate for the loss of minor splicing activity.

Molecular mechanisms in governing genomic stability and tumor suppression by the SETD2 H3K36 methyltransferase

Epigenetic dysregulation is an important contributor to carcinogenesis. This is not surprising, as chromatin-genomic DNA organized around structural histone scaffolding-serves as the template on which occurs essential nuclear processes, such as transcription, DNA replication and DNA repair. Histone H3 lysine 36 (H3K36) methyltransferases, such as the SET-domain 2 protein (SETD2), have emerged as critical tumor suppressors. Previous work on mammalian SETD2 and its counterpart in model organisms, Set2, has highlighted the role of this protein in governing genomic stability through transcriptional elongation and splicing, as well as in DNA damage response processes and cell cycle progression. A compendium of SETD2 mutations have been documented, garnered from sequenced cancer patient genome data, and these findings underscore the cancer-driving properties of SETD2 loss-of-function. In this review, we consolidate the molecular mechanisms regulated by SETD2/Set2 and discuss evidence of its dysregulation in tumorigenesis. Insight into the genetic interactions that exist between SETD2 and various canonical intracellular signaling pathways has not only empowered pharmacological intervention by taking advantage of synthetic lethality but underscores SETD2 as a druggable target for precision cancer therapy.

The "guiding" principles of noncoding RNA function

The human genome is pervasively transcribed and yet only a small fraction of these RNAs (less than 2%) are known to code for proteins. The vast majority of the RNAs are classified as noncoding RNAs (ncRNAs) and are further subgrouped as small (shorter than 200 bases) and long noncoding RNAs. The ncRNAs have been identified in all three domains of life and regulate diverse cellular processes through transcriptional and posttranscriptional gene regulation. Most of these RNAs work in conjunction with proteins forming a wide array of base pairing interactions. The determinants of these base pairing interactions are now becoming more evident and show striking similarities among the diverse group of ncRNAs. Here we present a mechanistic overview of pairing between RNA-RNA or RNA-DNA that dictates the function of ncRNAs; we provide examples to illustrate that ncRNAs work through shared evolutionary mechanisms that encompasses a guide-target interaction, involving not only classical Watson-Crick but also noncanonical Wobble and Hoogsteen base pairing. We also highlight the similarities in target selection, proofreading, and the ruler mechanism of ncRNA-protein complexes that confers target specificity and target site selection. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA-Based Catalysis > RNA-Mediated Cleavage RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

The Chloroplast Trans-Splicing RNA-Protein Supercomplex from the Green Alga Chlamydomonas reinhardtii

In eukaryotes, RNA trans-splicing is a significant RNA modification process for the end-to-end ligation of exons from separately transcribed primary transcripts to generate mature mRNA. So far, three different categories of RNA trans-splicing have been found in organisms within a diverse range. Here, we review trans-splicing of discontinuous group II introns, which occurs in chloroplasts and mitochondria of lower eukaryotes and plants. We discuss the origin of intronic sequences and the evolutionary relationship between chloroplast ribonucleoprotein complexes and the nuclear spliceosome. Finally, we focus on the ribonucleoprotein supercomplex involved in trans-splicing of chloroplast group II introns from the green alga Chlamydomonas reinhardtii. This complex has been well characterized genetically and biochemically, resulting in a detailed picture of the chloroplast ribonucleoprotein supercomplex. This information contributes substantially to our understanding of the function of RNA-processing machineries and might provide a blueprint for other splicing complexes involved in trans- as well as cis-splicing of organellar intron RNAs.

Structural and functional modularity of the U2 snRNP in pre-mRNA splicing

The U2 small nuclear ribonucleoprotein (snRNP) is an essential component of the spliceosome, the cellular machine responsible for removing introns from precursor mRNAs (pre-mRNAs) in all eukaryotes. U2 is an extraordinarily dynamic splicing factor and the most frequently mutated in cancers. Cryo-electron microscopy (cryo-EM) has transformed our structural and functional understanding of the role of U2 in splicing. In this review, we synthesize these and other data with respect to a view of U2 as an assembly of interconnected functional modules. These modules are organized by the U2 small nuclear RNA (snRNA) for roles in spliceosome assembly, intron substrate recognition, and protein scaffolding. We describe new discoveries regarding the structure of U2 components and how the snRNP undergoes numerous conformational and compositional changes during splicing. We specifically highlight large scale movements of U2 modules as the spliceosome creates and rearranges its active site. U2 serves as a compelling example for how cellular machines can exploit the modular organization and structural plasticity of an RNP.

The regulation properties of RNA secondary structure in alternative splicing

The RNA secondary structure is important for many functional processes in the cell. The secondary and tertiary structures of cellular RNAs are essential for the activity of these molecules in processes such as transcription, splicing, translation, and localization. New high-throughput analytical methods, including next generation sequencing, have allowed for the in-depth characterization of the 'RNA structurome': a new term describing how the RNA structure controls the activity of RNA by itself and how it regulates the expression of genes. In this review, we present many examples of the influence of structural motifs of RNA, long range interactions and global RNA structure on the alternative splicing processes. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Global donor and acceptor splicing site kinetics in human cells

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

Dynamic stacking of an expected branch point adenosine in duplexes containing pseudouridine-modified or unmodified U2 snRNA sites

The pre-mRNA branch point sequence (BPS) anneals with a pseudouridine-modified region of the U2 small nuclear (sn)RNA, and offers a 2' hydroxyl group of a bulged adenosine as the nucleophile for the first catalytic step of pre-mRNA splicing. To increase our structural understanding of branch site selection, we characterized a duplex containing a BPS sequence that is common among multicellular eukaryotes (5'-UACUGAC-3') and the complementary U2 snRNA site using NMR. A major conformation of the expected branch site adenosine stacked within the duplex and paired with the conserved pseudouridine of the U2 snRNA strand. In contrast, the guanosine preceding the branch site appeared flexible and had weak contacts with the surrounding nucleotides. Pseudouridine-modified and unmodified U2 snRNA-BPS-containing duplexes remained structurally similar. These results highlight the importance of auxiliary factors to achieve the active bulged conformation of the branch site nucleophile for the first step of pre-mRNA splicing.

Chemical Inhibition of Pre-mRNA Splicing in Living Saccharomyces cerevisiae

The spliceosome mediates precursor mRNA splicing in eukaryotes, including the model organism Saccharomyces cerevisiae (yeast). Despite decades of study, no chemical inhibitors of yeast splicing in vivo are available. We have developed a system to efficiently inhibit splicing and block proliferation in living yeast cells using compounds that target the human spliceosome protein SF3B1. Potent inhibition is observed in yeast expressing a chimeric protein containing portions of human SF3B1. However, only a single point mutation in the yeast homolog of SF3B1 is needed for selective inhibition of splicing by pladienolide B, herboxidiene, or meayamycin in liquid culture. Mutations that enable inhibition also improve splicing of branch sites containing mismatches between the intron and small nuclear RNA-suggesting a link between inhibitor sensitivity and usage of weak branch sites in humans. This approach provides powerful new tools for manipulating splicing in live yeast and studies of spliceosome inhibitors.

Functional analysis of Hsh155/SF3b1 interactions with the U2 snRNA/branch site duplex

SF3b1 is an essential component of the U2 snRNP implicated in branch site (BS) recognition and found to be frequently mutated in several human cancers. While recent structures of yeast and human SF3b1 have revealed its molecular architecture, the importance of specific RNA:protein contacts and conformational changes remains largely uncharacterized. Here, we performed mutational analysis of yeast SF3b1, guided by recent structures of the spliceosome. We find that conserved amino acids contacting the U2 snRNA backbone of the U2/BS duplex are nonessential, and that yeast can tolerate truncation of the HEAT repeats containing these amino acids. The pocket housing the branchpoint adenosine (BP-A) is also amenable to mutation despite strong conservation. However, mutations that support viability can still lead to defects in splicing pre-mRNAs with nonconsensus BS substitutions found at -3, -2, -1, and +1 positions relative to the BP-A or at the branchpoint position. Through the generation of yeast and human chimeric proteins, we further defined the functionally conserved regions of Hsh155 as well as identify changes in BS usage resulting from inclusion of human SF3b1 HEAT repeats. Moreover, these chimeric proteins confer a sensitivity to small molecule inhibition by pladienolide B to yeast splicing. Together, these data reveal the importance of individual contacts of Hsh155/SF3b1 to the U2/BS duplex and define their contribution to BS usage by the spliceosome.

BPP: a sequence-based algorithm for branch point prediction

Motivation

Although high-throughput sequencing methods have been proposed to identify splicing branch points in the human genome, these methods can only detect a small fraction of the branch points subject to the sequencing depth, experimental cost and the expression level of the mRNA. An accurate computational model for branch point prediction is therefore an ongoing objective in human genome research.

Results

We here propose a novel branch point prediction algorithm that utilizes information on the branch point sequence and the polypyrimidine tract. Using experimentally validated data, we demonstrate that our proposed method outperforms existing methods. Availability and implementation: https://github.com/zhqingit/BPP.

Contact

djguo@cuhk.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

SF3b1 mutations associated with myelodysplastic syndromes alter the fidelity of branchsite selection in yeast

RNA and protein components of the spliceosome work together to identify the 5΄ splice site, the 3΄ splice site, and the branchsite (BS) of nascent pre-mRNA. SF3b1 plays a key role in recruiting the U2 snRNP to the BS. Mutations in human SF3b1 have been linked to many diseases such as myelodysplasia (MDS) and cancer. We have used SF3b1 mutations associated with MDS to interrogate the role of the yeast ortholog, Hsh155, in BS selection and splicing. These alleles change how the spliceosome recognizes the BS and alter splicing when nonconsensus nucleotides are present at the −2, −1 and +1 positions relative to the branchpoint adenosine. This indicates that changes in BS usage observed in humans with SF3b1 mutations may result from perturbation of a conserved mechanism of BS recognition. Notably, different HSH155 alleles elicit disparate effects on splicing: some increase the fidelity of BS selection while others decrease fidelity. Our data support a model wherein conformational changes in SF3b1 promote U2 association with the BS independently of the action of the DEAD-box ATPase Prp5. We propose that SF3b1 functions to stabilize weak U2/BS duplexes to drive spliceosome assembly and splicing.

CryoEM structures of spliceosomal complexes reveal the molecular mechanism of pre-mRNA splicing

The spliceosome is an intricate molecular machine which catalyses the removal of introns from eukaryotic mRNA precursors by two trans-esterification reactions (branching and exon ligation) to produce mature mRNA with uninterrupted protein coding sequences. The structures of the spliceosome in several key states determined by electron cryo-microscopy have greatly advanced our understanding of its molecular mechanism. The catalytic RNA core is formed during the activation of the fully assembled B to Bact complex and remains largely unchanged throughout the splicing cycle. RNA helicases and step specific factors regulate docking and undocking of the substrates (branch site and 3' splice site) to the single RNA-based active site to catalyse the two trans-esterification reactions.

Radical probing of spliceosome assembly

Here we describe the synthesis and use of a directed hydroxyl radical probe, tethered to a pre-mRNA substrate, to map the structure of this substrate during the spliceosome assembly process. These studies indicate an early organization and proximation of conserved pre-mRNA sequences during spliceosome assembly. This methodology may be adapted to the synthesis of a wide variety of modified RNAs for use as probes of RNA structure and RNA-protein interaction.

Reversibly constraining spliceosome-substrate complexes by engineering disulfide crosslinks

The spliceosome is a highly dynamic mega-Dalton enzyme, formed in part by assembly of U snRNPs onto its pre-mRNA substrate transcripts. Early steps in spliceosome assembly are challenging to study biochemically and structurally due to compositional and conformational dynamics. We detail an approach to covalently and reversibly constrain or trap non-covalent pre-mRNA/protein spliceosome complexes. This approach involves engineering a single disulfide bond between a thiol-bearing cysteine sidechain and a proximal backbone phosphate of the pre-mRNA, site-specifically modified with an N-thioalkyl moiety. When distance and angle between reactants is optimal, the sidechain will react with the single N-thioalkyl to form a crosslink upon oxidation. We provide protocols detailing how this has been applied successfully to trap an 11-subunit RNA-protein assembly, the human U1 snRNP, in complex with a pre-mRNA.

RIsearch2: suffix array-based large-scale prediction of RNA-RNA interactions and siRNA off-targets

Intermolecular interactions of ncRNAs are at the core of gene regulation events, and identifying the full map of these interactions bears crucial importance for ncRNA functional studies. It is known that RNA-RNA interactions are built up by complementary base pairings between interacting RNAs and high level of complementarity between two RNA sequences is a powerful predictor of such interactions. Here, we present RIsearch2, a large-scale RNA-RNA interaction prediction tool that enables quick localization of potential near-complementary RNA-RNA interactions between given query and target sequences. In contrast to previous heuristics which either search for exact matches while including G-U wobble pairs or employ simplified energy models, we present a novel approach using a single integrated seed-and-extend framework based on suffix arrays. RIsearch2 enables fast discovery of candidate RNA-RNA interactions on genome/transcriptome-wide scale. We furthermore present an siRNA off-target discovery pipeline that not only predicts the off-target transcripts but also computes the off-targeting potential of a given siRNA. This is achieved by combining genome-wide RIsearch2 predictions with target site accessibilities and transcript abundance estimates. We show that this pipeline accurately predicts siRNA off-target interactions and enables off-targeting potential comparisons between different siRNA designs. RIsearch2 and the siRNA off-target discovery pipeline are available as stand-alone software packages from http://rth.dk/resources/risearch.

Genetics and biochemistry remain essential in the structural era of the spliceosome

The spliceosome is not a single macromolecular machine. Rather it is a collection of dynamic heterogeneous subcomplexes that rapidly interconvert throughout the course of a typical splicing cycle. Because of this, for many years the only high resolution structures of the spliceosome available were of smaller, isolated protein or RNA components. Consequently much of our current understanding of the spliceosome derives from biochemical and genetic techniques. Now with the publication of multiple, high resolution structures of the spliceosome, some question the relevance of traditional biochemical and genetic techniques to the splicing field. We argue such techniques are not only relevant, but vital for an in depth mechanistic understanding of pre-mRNA splicing.

Spliceosomal DEAH-Box ATPases Remodel Pre-mRNA to Activate Alternative Splice Sites

During pre-mRNA splicing, a central step in the expression and regulation of eukaryotic genes, the spliceosome selects splice sites for intron excision and exon ligation. In doing so, the spliceosome must distinguish optimal from suboptimal splice sites. At the catalytic stage of splicing, suboptimal splice sites are repressed by the DEAH-box ATPases Prp16 and Prp22. Here, using budding yeast, we show that these ATPases function further by enabling the spliceosome to search for and utilize alternative branch sites and 3' splice sites. The ATPases facilitate this search by remodeling the splicing substrate to disengage candidate splice sites. Our data support a mechanism involving 3' to 5' translocation of the ATPases along substrate RNA and toward a candidate site, but, surprisingly, not across the site. Thus, our data implicate DEAH-box ATPases in acting at a distance by pulling substrate RNA from the catalytic core of the spliceosome.

Towards understanding pre-mRNA splicing mechanisms and the role of SR proteins

Alternative pre-mRNA splicing provides a source of vast protein diversity by removing non-coding sequences (introns) and accurately linking different exonic regions in the correct reading frame. The regulation of alternative splicing is essential for various cellular functions in both pathological and physiological conditions. In eukaryotic cells, this process is commonly used to increase proteomic diversity and to control gene expression either co- or post-transcriptionally. Alternative splicing occurs within a megadalton-sized, multi-component machine consisting of RNA and proteins; during the splicing process, this complex undergoes dynamic changes via RNA-RNA, protein-protein and RNA-protein interactions. Co-transcriptional splicing functionally integrates the transcriptional machinery, thereby enabling the two processes to influence one another, whereas post-transcriptional splicing facilitates the coupling of RNA splicing with post-splicing events. This review addresses the structural aspects of spliceosomes and the mechanistic implications of their stepwise assembly on the regulation of pre-mRNA splicing. Moreover, the role of phosphorylation-based, signal-induced changes in the regulation of the splicing process is demonstrated.

Prp4 Kinase Grants the License to Splice: Control of Weak Splice Sites during Spliceosome Activation

The genome of the fission yeast Schizosaccharomyces pombe encodes 17 kinases that are essential for cell growth. These include the cell-cycle regulator Cdc2, as well as several kinases that coordinate cell growth, polarity, and morphogenesis during the cell cycle. In this study, we further characterized another of these essential kinases, Prp4, and showed that the splicing of many introns is dependent on Prp4 kinase activity. For detailed characterization, we chose the genes res1 and ppk8, each of which contains one intron of typical size and position. Splicing of the res1 intron was dependent on Prp4 kinase activity, whereas splicing of the ppk8 intron was not. Extensive mutational analyses of the 5’ splice site of both genes revealed that proper transient interaction with the 5’ end of snRNA U1 governs the dependence of splicing on Prp4 kinase activity. Proper transient interaction between the branch sequence and snRNA U2 was also important. Therefore, the Prp4 kinase is required for recognition and efficient splicing of introns displaying weak exon1/5’ splice sites and weak branch sequences.

Prp4 is an essential protein kinase that is involved in the splicing of some introns. Using a conditional mutant of Prp4, we showed that a subset of genes, including several cell cycle–regulatory genes, are dependent on Prp4 for splicing. Furthermore, we could convert genes between Prp4-dependent and -independent states by introducing single-nucleotide mutations in the exon1/5’ splice sites and branch sequence of introns. This work shows that Prp4 activity is required for splicing surveillance in a subset of mRNAs.

Mammalian splicing factor SF1 interacts with SURP domains of U2 snRNP-associated proteins

Splicing factor 1 (SF1) recognizes the branch point sequence (BPS) at the 3′ splice site during the formation of early complex E, thereby pre-bulging the BPS adenosine, thought to facilitate subsequent base-pairing of the U2 snRNA with the BPS. The 65-kDa subunit of U2 snRNP auxiliary factor (U2AF65) interacts with SF1 and was shown to recruit the U2 snRNP to the spliceosome. Co-immunoprecipitation experiments of SF1-interacting proteins from HeLa cell extracts shown here are consistent with the presence of SF1 in early splicing complexes. Surprisingly almost all U2 snRNP proteins were found associated with SF1. Yeast two-hybrid screens identified two SURP domain-containing U2 snRNP proteins as partners of SF1. A short, evolutionarily conserved region of SF1 interacts with the SURP domains, stressing their role in protein–protein interactions. A reduction of A complex formation in SF1-depleted extracts could be rescued with recombinant SF1 containing the SURP-interaction domain, but only partial rescue was observed with SF1 lacking this sequence. Thus, SF1 can initially recruit the U2 snRNP to the spliceosome during E complex formation, whereas U2AF65 may stabilize the association of the U2 snRNP with the spliceosome at later times. In addition, these findings may have implications for alternative splicing decisions.

Dynamic Contacts of U2, RES, Cwc25, Prp8 and Prp45 Proteins with the Pre-mRNA Branch-Site and 3' Splice Site during Catalytic Activation and Step 1 Catalysis in Yeast Spliceosomes

Little is known about contacts in the spliceosome between proteins and intron nucleotides surrounding the pre-mRNA branch-site and their dynamics during splicing. We investigated protein-pre-mRNA interactions by UV-induced crosslinking of purified yeast B(act) spliceosomes formed on site-specifically labeled pre-mRNA, and analyzed their changes after conversion to catalytically-activated B* and step 1 C complexes, using a purified splicing system. Contacts between nucleotides upstream and downstream of the branch-site and the U2 SF3a/b proteins Prp9, Prp11, Hsh49, Cus1 and Hsh155 were detected, demonstrating that these interactions are evolutionarily conserved. The RES proteins Pml1 and Bud13 were shown to contact the intron downstream of the branch-site. A comparison of the B(act) crosslinking pattern versus that of B* and C complexes revealed that U2 and RES protein interactions with the intron are dynamic. Upon step 1 catalysis, Cwc25 contacts with the branch-site region, and enhanced crosslinks of Prp8 and Prp45 with nucleotides surrounding the branch-site were observed. Cwc25's step 1 promoting activity was not dependent on its interaction with pre-mRNA, indicating it acts via protein-protein interactions. These studies provide important insights into the spliceosome's protein-pre-mRNA network and reveal novel RNP remodeling events during the catalytic activation of the spliceosome and step 1 of splicing.

Intronic deletions of tva receptor gene decrease the susceptibility to infection by avian sarcoma and leukosis virus subgroup A

The group of avian sarcoma and leukosis virus (ASLV) in chickens contains six highly related subgroups, A to E and J. Four genetic loci, tva, tvb, tvc and tvj, encode for corresponding receptors that determine the susceptibility to the ASLV subgroups. The prevalence of ASLV in hosts may have imposed strong selection pressure toward resistance to ASLV infection, and the resistant alleles in all four receptor genes have been identified. In this study, two new alleles of the tva receptor gene, tva^r5 and tva^r6, with similar intronic deletions were identified in Chinese commercial broilers. These natural mutations delete the deduced branch point signal within the first intron, disrupting mRNA splicing of the tva receptor gene and leading to the retention of intron 1 and introduction of premature TGA stop codons in both the longer and shorter tva isoforms. As a result, decreased susceptibility to subgroup A ASLV in vitro and in vivo was observed in the subsequent analysis. In addition, we identified two groups of heterozygous allele pairs which exhibited quantitative differences in host susceptibility to ASLV-A. This study demonstrated that defective splicing of the tva receptor gene can confer genetic resistance to ASLV subgroup A in the host.

Crystal structure of a eukaryotic group II intron lariat

The formation of branched lariat RNA is an evolutionarily conserved feature of splicing reactions for both group II and spliceosomal introns. The lariat is important for the fidelity of 5' splice-site selection and consists of a 2'-5' phosphodiester bond between a bulged adenosine and the 5' end of the intron. To gain insight into this ubiquitous intramolecular linkage, we determined the crystal structure of a eukaryotic group IIB intron in the lariat form at 3.7 Å. This revealed that two tandem tetraloop-receptor interactions, η-η' and π-π', place domain VI in the core to position the lariat bond in the post-catalytic state. On the basis of structural and biochemical data, we propose that π-π' is a dynamic interaction that mediates the transition between the two steps of splicing, with η-η' serving an ancillary role. The structure also reveals a four-magnesium-ion cluster involved in both catalysis and positioning of the 5' end. Given the evolutionary relationship between group II and nuclear introns, it is likely that this active site configuration exists in the spliceosome as well.

Splicing factor 3B1 hypomethylation is associated with altered SF3B1 transcript expression in older humans

Ageing in man is associated with changes to the splicing factor pool. A proportion of splicing factors are regulated during ageing by mechanisms involving the Ataxia Telangiectasia Mutated (ATM) gene, but the factors that determine the remaining proportion have yet to be identified. DNA methylation is known to be an important regulatory mechanism of gene expression. We assessed age-associated methylation and expression levels for 27 splicing factor genes, in peripheral blood samples from the InCHIANTI study. Examination of splicing patterns at specific loci was examined in a second cohort, the Exeter 10000 study. 27/502 methylation probes in 17 different genes were associated with age. Most changes were not associated with transcript expression levels or splicing patterns, but hypomethylation of the SF3B1 promoter region was found to mediate 53% of the relationship between age and transcript expression at this locus (p=0.02). DNA methylation does not appear to play a major role in regulation of the splicing factors, but changes in SF3B1 expression may be attributable to promoter hypomethylation at this locus. SF3B1 encodes a critical component of the U2 snRNP; altered expression of this gene may therefore contribute to the loss of regulated mRNA splicing that occurs with age.

The thermodynamic patterns of eukaryotic genes suggest a mechanism for intron-exon recognition

The essential cis- and trans-acting elements required for RNA splicing have been defined, however, the detailed molecular mechanisms underlying intron-exon recognition are still unclear. Here we demonstrate that the ratio between stability of mRNA/DNA and DNA/DNA duplexes near 3'-spice sites is a characteristic feature that can contribute to intron-exon differentiation. Remarkably, throughout all transcripts, the most unstable mRNA/DNA duplexes, compared with the corresponding DNA/DNA duplexes, are situated upstream of the 3'-splice sites and include the polypyrimidine tracts. This characteristic instability is less pronounced in weak alternative splice sites and disease-associated cryptic 3'-splice sites. Our results suggest that this thermodynamic pattern can prevent the re-annealing of mRNA to the DNA template behind the RNA polymerase to ensure access of the splicing machinery to the polypyrimidine tract and the branch point. In support of this mechanism, we demonstrate that RNA/DNA duplex formation at this region prevents pre-spliceosome A complex assembly.

Now on display: a gallery of group II intron structures at different stages of catalysis

Group II introns are mobile genetic elements that self-splice and retrotranspose into DNA and RNA. They are considered evolutionary ancestors of the spliceosome, the ribonucleoprotein complex essential for pre-mRNA processing in higher eukaryotes. Over a 20-year period, group II introns have been characterized first genetically, then biochemically, and finally by means of X-ray crystallography. To date, 17 crystal structures of a group II intron are available, representing five different stages of the splicing cycle. This review provides a framework for classifying and understanding these new structures in the context of the splicing cycle. Structural and functional implications for the spliceosome are also discussed.

Impact of base pair identity 5' to the spliceosomal branch site adenosine on branch site conformation

The branch site helix from Saccharomyces cerevisiae with pseudouridine (ψ) incorporated in a phylogenetically conserved position of U2 snRNA features an extrahelical branch site adenosine (A) that forms a base triple interaction with the minor groove edge of a widely conserved purine(U2 strand)-pyrimidine(intron strand) (R(U2)-Y(intron)) base pair two positions upstream. In these studies, NMR spectra of a duplex in which 2-aminopurine (2ap), a fluorescent analog of adenine lacking the proposed hydrogen bond donor, was substituted for the branch site A, indicated that the substitution does not alter the extrahelical position of the branch site residue; thus, it appears that a hydrogen bond between the adenine amino group and the R-Y pair is not obligatory for stabilization of the extrahelical conformation. In contrast, reversal of the orientation of A(U2)-U(intron) to U(U2)-A(intron) resulted in an intrahelical position for the branch site A or 2ap. Fluorescence intensity of 2ap substituted for the branch site A with the original R(U2)-Y(intron) orientation (AU or GC) was high, consistent with an extrahelical position, whereas fluorescence in helices with the reversed R-Y orientation, or with a mismatched pair (A-U → G•A or U•C), was markedly quenched, implying that the residue was stacked in the helix. The A 5' to the branch site residue was not extrahelical in any of the duplexes. These findings suggest that the R(U2)-Y(intron) base pair orientation in the ψ-dependent branch site helix plays an important role in positioning the branch site A for recognition and/or function.

Intron definition and a branch site adenosine at nt 385 control RNA splicing of HPV16 E6*I and E7 expression

HPV16 E6 and E7, two viral oncogenes, are expressed from a single bicistronic pre-mRNA. In this report, we provide the evidence that the bicistronic pre-mRNA intron 1 contains three 5' splice sites (5' ss) and three 3' splice sites (3' ss) normally used in HPV16(+) cervical cancer and its derived cell lines. The choice of two novel alternative 5' ss (nt 221 5' ss and nt 191 5' ss) produces two novel isoforms of E6E7 mRNAs (E6*V and E6*VI). The nt 226 5' ss and nt 409 3' ss is preferentially selected over the other splice sites crossing over the intron to excise a minimal length of the intron in RNA splicing. We identified AACAAAC as the preferred branch point sequence (BPS) and an adenosine at nt 385 (underlined) in the BPS as a branch site to dictate the selection of the nt 409 3' ss for E6*I splicing and E7 expression. Introduction of point mutations into the mapped BPS led to reduced U2 binding to the BPS and thereby inhibition of the second step of E6E7 splicing at the nt 409 3' ss. Importantly, the E6E7 bicistronic RNA with a mutant BPS and inefficient splicing makes little or no E7 and the resulted E6 with mutations of (91)QYNK(94) to (91)PSFW(94) displays attenuate activity on p53 degradation. Together, our data provide structural basis of the E6E7 intron 1 for better understanding of how viral E6 and E7 expression is regulated by alternative RNA splicing. This study elucidates for the first time a mapped branch point in HPV16 genome involved in viral oncogene expression.

Reduced fidelity of branch point recognition and alternative splicing induced by the anti-tumor drug spliceostatin A

Spliceostatin A (SSA) is a stabilized derivative of a Pseudomonas bacterial fermentation product that displays potent anti-proliferative and anti-tumor activities in cancer cells and animal models. The drug inhibits pre-mRNA splicing in vitro and in vivo and binds SF3b, a protein subcomplex of U2 small nuclear ribonucleoprotein (snRNP), which is essential for recognition of the pre-mRNA branch point. We report that SSA prevents interaction of an SF3b 155-kDa subunit with the pre-mRNA, concomitant with nonproductive recruitment of U2 snRNP to sequences 5' of the branch point. Differences in base-pairing potential with U2 snRNA in this region lead to different sensitivity of 3' splice sites to SSA, and to SSA-induced changes in alternative splicing. Indeed, rather than general splicing inhibition, splicing-sensitive microarray analyses reveal specific alternative splicing changes induced by the drug that significantly overlap with those induced by knockdown of SF3b 155. These changes lead to down-regulation of genes important for cell division, including cyclin A2 and Aurora A kinase, thus providing an explanation for the anti-proliferative effects of SSA. Our results reveal a mechanism that prevents nonproductive base-pairing interactions in the spliceosome, and highlight the regulatory and cancer therapeutic potential of perturbing the fidelity of splice site recognition.

RNA labeling, conjugation and ligation

Advances in RNA nanotechnology will depend on the ability to manipulate, probe the structure and engineer the function of RNA with high precision. This article reviews current abilities to incorporate site-specific labels or to conjugate other useful molecules to RNA either directly or indirectly through post-synthetic labeling methodologies that have enabled a broader understanding of RNA structure and function. Readily applicable modifications to RNA can range from isotopic labels and fluorescent or other molecular probes to protein, lipid, glycoside or nucleic acid conjugates that can be introduced using combinations of synthetic chemistry, enzymatic incorporation and various conjugation chemistries. These labels, conjugations and ligations to RNA are quintessential for further investigation and applications of RNA as they enable the visualization, structural elucidation, localization, and biodistribution of modified RNA.

Structural model of the p14/SF3b155 · branch duplex complex

Human p14 (SF3b14), a component of the spliceosomal U2 snRNP, interacts directly with the pre-mRNA branch adenosine within the context of the bulged duplex formed between the pre-mRNA branch region and U2 snRNA. This association occurs early in spliceosome assembly and persists within the fully assembled spliceosome. Analysis of the crystal structure of a complex containing p14 and a peptide derived from p14-associated SF3b155 combined with the results of cross-linking studies has suggested that the branch nucleotide interacts with a pocket on a non-canonical RNA binding surface formed by the complex. Here we report a structural model of the p14 · bulged duplex interaction based on a combination of X-ray crystallography of an adenine p14/SF3b155 peptide complex, biochemical comparison of a panel of disulfide cross-linked protein-RNA complexes, and small-angle X-ray scattering (SAXS). These studies reveal specific recognition of the branch adenosine within the p14 pocket and establish the orientation of the bulged duplex RNA bound on the protein surface. The intimate association of one surface of the bulged duplex with the p14/SF3b155 peptide complex described by this model buries the branch nucleotide at the interface and suggests that p14 · duplex interaction must be disrupted before the first step of splicing.

An SF1 affinity model to identify branch point sequences in human introns

Splicing factor 1 (SF1) binds to the branch point sequence (BPS) of mammalian introns and is believed to be important for the splicing of some, but not all, introns. To help identify BPSs, particularly those that depend on SF1, we generated a BPS profile model in which SF1 binding affinity data, validated by branch point mapping, were iteratively incorporated into computational models. We searched a data set of 117 499 human introns for best matches to the SF1 Affinity Model above a threshold, and counted the number of matches at each intronic position. After subtracting a background value, we found that 87.9% of remaining high-scoring matches identified were located in a region upstream of 3′-splice sites where BPSs are typically found. Since U2AF65 recognizes the polypyrimidine tract (PPT) and forms a cooperative RNA complex with SF1, we combined the SF1 model with a PPT model computed from high affinity binding sequences for U2AF65. The combined model, together with binding site location constraints, accurately identified introns bound by SF1 that are candidates for SF1-dependent splicing.

Regulation of pre-mRNA splicing in Xenopus oocytes by targeted 2'-O-methylation

The 2'-OH group of the branch point adenosine is a key moiety to initiate pre-mRNA splicing. We use RNA-guided RNA modification to target the pre-mRNA branch point adenosine for 2'-O-methylation, with the aim of blocking pre-mRNA splicing in vertebrate cells. We show that, under certain conditions, injection of a branch point-specific artificial box C/D RNA into Xenopus oocytes effectively 2'-O-methylates adenovirus pre-mRNA at the target nucleotide. However, 2'-O-methylation at the authentic branch point activates a host of cryptic branch points, thus allowing splicing to continue. These cryptic sites are mapped, and mutated. Upon injection, pre-mRNA free of cryptic branch points fails to splice when the branch point-specific box C/D RNA is present. However, 2'-O-methylation at the branch point does not prevent pre-mRNA from being assembled into pre-catalytic spliceosome-like complexes prior to the first chemical step of splicing. Our results demonstrate that RNA-guided pre-mRNA modification can occur in the nucleoplasm of vertebrate cells, thus offering a powerful tool for molecular biology research.

Isotope labeling strategies for NMR studies of RNA

The known biological functions of RNA have expanded in recent years and now include gene regulation, maintenance of sub-cellular structure, and catalysis, in addition to propagation of genetic information. As for proteins, RNA function is tightly correlated with structure. Unlike proteins, structural information for larger, biologically functional RNAs is relatively limited. NMR signal degeneracy, relaxation problems, and a paucity of long-range (1)H-(1)H dipolar contacts have limited the utility of traditional NMR approaches. Selective isotope labeling, including nucleotide-specific and segmental labeling strategies, may provide the best opportunities for obtaining structural information by NMR. Here we review methods that have been developed for preparing and purifying isotopically labeled RNAs, as well as NMR strategies that have been employed for signal assignment and structure determination.

Spliceosome structure: piece by piece

Processing of pre-mRNAs by RNA splicing is an essential step in the maturation of protein coding RNAs in eukaryotes. Structural studies of the cellular splicing machinery, the spliceosome, are a major challenge in structural biology due to the size and complexity of the splicing ensemble. Specifically, the structural details of splice site recognition and the architecture of the spliceosome active site are poorly understood. X-ray and NMR techniques have been successfully used to address these questions defining the structure of individual domains, isolated splicing proteins, spliceosomal RNA fragments and recently the U1 snRNP multiprotein.RNA complex. These results combined with extant biochemical and genetic data have yielded important insights as well as posing fresh questions with respect to the regulation and mechanism of this critical gene regulatory process.

Insights into branch nucleophile positioning and activation from an orthogonal pre-mRNA splicing system in yeast

The duplex formed between the branch site (BS) of a spliceosomal intron and its cognate sequence in U2 snRNA is important for spliceosome assembly and the first catalytic step of splicing. We describe the development of an orthogonal BS-U2 system in S. cerevisiae in which spliceosomes containing a grossly substituted second-copy U2 snRNA mediate the in vivo splicing of a single reporter transcript carrying a cognate substitution. Systematic use of this approach to investigate requirements for branching catalysis reveals considerable flexibility in the sequence of the BS-U2 duplex and its positioning relative to the catalytic center. Branching efficiency depends on the identity of the branch nucleotide, its position within the BS-U2 duplex, and its distance from U2/U6 helix Ia. These results provide insights into substrate selection during spliceosomal branching catalysis; additionally, this system provides a foundation and tool for future mechanistic splicing research.

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5' splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5' and 3' splice sites, and the activation of cryptic U12-dependent branch/3' splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3' splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.

Dimerization and protein binding specificity of the U2AF homology motif of the splicing factor Puf60

PUF60 is an essential splicing factor functionally related and homologous to U2AF(65). Its C-terminal domain belongs to the family of U2AF (U2 auxiliary factor) homology motifs (UHM), a subgroup of RNA recognition motifs that bind to tryptophan-containing linear peptide motifs (UHM ligand motifs, ULMs) in several nuclear proteins. Here, we show that the Puf60 UHM is mainly monomeric in physiological buffer, whereas its dimerization is induced upon the addition of SDS. The crystal structure of PUF60-UHM at 2.2 angstroms resolution, NMR data, and mutational analysis reveal that the dimer interface is mediated by electrostatic interactions involving a flexible loop. Using glutathione S-transferase pulldown experiments, isothermal titration calorimetry, and NMR titrations, we find that Puf60-UHM binds to ULM sequences in the splicing factors SF1, U2AF65, and SF3b155. Compared with U2AF65-UHM, Puf60-UHM has distinct binding preferences to ULMs in the N terminus of SF3b155. Our data suggest that the functional cooperativity between U2AF65 and Puf60 may involve simultaneous interactions of the two proteins with SF3b155.

Splice site requirements and switches in plants

Intron sequences in nuclear pre-mRNAs are excised with either the major U2 snRNA-dependent spliceosomal pathway or the minor U12 snRNA-dependent spliceosomal pathway that exist in most eukaryotic organisms. While the predominant dinucleotides bordering each of these types of introns and the catalytic mechanism used in their excision are conserved in plants and animals, several features aiding in the recognition of plant introns are distinct from those in animals and yeast. Along with their short length, high AU content and high variation in their 5' and 3' splice sites and branchpoint consensus sequences are the most prominent characteristics of plant introns. Detailed surveys of site-directed mutant introns tested in vivo and chemically induced and naturally mutant introns analyzed in planta emphasize the effects of changing individual nucleotides in these splice site consensus sequences and highlight a number of noncanonical dinucleotides that are functional in plant systems.

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.

Spliceosomal immunophilins

The spliceosome is a dynamic, macromolecular complex, which removes non-protein-coding introns from pre-mRNA to form mature mRNA in a process known as splicing. This ribonucleoprotein assembly is comprised of five uridine-rich small nuclear RNAs (snRNAs) as well as over 300 proteins. In humans, several of the known proteinaceous splicing factors are members of the immunophilin superfamily. Immunophilins are peptidyl-prolyl cis-trans isomerases that catalyze the conversion of proteins from cis to trans at Xaa-Pro bonds. Our review of the data indicates that some members of this protein family are activators of spliceosomal proteins by way of folding and transport.

Solution structure of the second RNA recognition motif (RRM) domain of murine T cell intracellular antigen-1 (TIA-1) and its RNA recognition mode

T cell intracellular antigen-1 (TIA-1), an apoptosis promoting factor, functions as a splicing regulator for the Fas pre-mRNA. TIA-1 possesses three RNA recognition motifs (RRMs) and a glutamine-rich domain. The second RRM (RRM2) is necessary and sufficient for tight, sequence-specific binding to the uridine-rich sequences buried around the 5' splice sites. In the present study, we solved the solution structure of the murine TIA-1 RRM2 by heteronuclear-nuclear magnetic resonance spectroscopy. The TIA-1 RRM2 adopts the RRM fold (betaalphabetabetaalphabeta) and possesses an extra beta-strand between beta2 and beta3, which forms an additional beta-sheet with the C-terminal part of beta2. We refer to this structure as the beta2-beta2' beta-loop. Interestingly, this characteristic beta-loop structure is conserved among a number of RRMs, including the U2AF65 RRM2 and the Sex-lethal RRM1 and RRM2, which also bind to uridine-rich RNAs. Furthermore, we identified a new sequence motif in the beta2-beta2' beta-loop, the DxxT motif. Chemical shift perturbation analyses of both the main and side chains upon binding to the uridine pentamer RNA revealed that most of the beta-sheet surface, including the beta2-beta2' beta-loop, is involved in the RNA binding. An investigation of the chemical shift perturbation revealed similarity in the RNA recognition modes between the TIA-1 and U2AF65 RRMs.

X-ray structures of U2 snRNA-branchpoint duplexes containing conserved pseudouridines

A pseudouridine-modified region of the U2 small nuclear (sn)RNA anneals with the intronic branchpoint sequence and positions a bulged adenosine to serve as the nucleophile in the first chemical step of pre-mRNA splicing. We have determined three X-ray structures of RNA oligonucleotides containing the pseudouridylated U2 snRNA and the branchpoint consensus sequences. The expected adenosine branchpoint is extrahelical in a 1.65 A resolution structure containing the mammalian consensus sequence variant and in a 2.10 A resolution structure containing a shortened Saccharomyces cerevisiae consensus sequence. The adenosine adjacent to the expected branchpoint is extrahelical in a third structure, which contains the intact yeast consensus sequence at 1.57 A resolution. The hydration and base stacking interactions mediated by the U2 snRNA pseudouridines correlate with the identity of the unpaired adenosine. The expected adenosine bulge is associated with a well-stacked pseudouridine, which is linked via an ordered water molecule to a neighboring nucleotide. In contrast, the bulge of the adjacent adenosine shifts the base stacking and disrupts the water-mediated interactions of the pseudouridine. These structural differences may contribute to the ability of the pseudouridine modification to promote the bulged conformation of the branch site adenosine and to enhance catalysis by snRNAs. Furthermore, iodide binding sites are identified adjacent to the unconventional bulged adenosine, and the structure of the mammalian consensus sequence variant provides a high-resolution view of a hydrated magnesium ion bound in a similar manner to a divalent cation binding site of the group II intron.

The CTD role in cotranscriptional RNA processing and surveillance

In higher eukaryotes, the production of mature messenger RNA that exits the nucleus to be translated into protein requires precise and extensive processing of the nascent transcript. The processing steps include 5'-end capping, splicing, and 3'-end formation. Pre-mRNA processing is coupled to transcription by mechanisms that are not well understood but involve the carboxyl-terminal domain (CTD) of the largest subunit of RNA polymerase II. This review focuses on recent findings that provide novel insight into the role of the CTD in promoting RNA processing and surveillance.

Human branch point consensus sequence is yUnAy

Yeast carries a strictly conserved branch point sequence (BPS) of UACUAAC, whereas the human BPS is degenerative and is less well characterized. The human consensus BPS has never been extensively explored in vitro to date. Here, we sequenced 367 clones of lariat RT-PCR products arising from 52 introns of 20 human housekeeping genes. Among the 367 clones, a misincorporated nucleotide at the branch point was observed in 181 clones, for which we can precisely pinpoint the branch point. The branch points were comprised of 92.3% A, 3.3% C, 1.7% G and 2.8% U. Our analysis revealed that the human consensus BPS is simply yUnAy, where the underlined is the branch point at position zero and the lowercase pyrimidines ('y') are not as well conserved as the uppercase U and A. We found that the branch points are located 21-34 nucleotides upstream of the 3' end of an intron in 83% clones. We also found that the polypyrimidine tract spans 4-24 nucleotides downstream of the branch point. Our analysis demonstrates that the human BPSs are more degenerative than we have expected and that the human BPSs are likely to be recognized in combination with the polypyrimidine tract and/or the other splicing cis-elements.

Transposition of two amino acids changes a promiscuous RNA binding protein into a sequence-specific RNA binding protein

In yeast (Saccharomyces cerevisiae), the branchpoint binding protein (BBP) recognizes the conserved yeast branchpoint sequence (UACUAAC) with a high level of specificity and affinity, while the human branchpoint binding protein (SF1) binds the less-conserved consensus branchpoint sequence (CURAY) in human introns with a lower level of specificity and affinity. To determine which amino acids in BBP provide the additional specificity and affinity absent in SF1, a panel of chimeric SF1 proteins was tested in RNA binding assays with wild-type and mutant RNA substrates. This approach revealed that the QUA2 domain of BBP is responsible for the enhanced RNA binding affinity and specificity displayed by BBP compared with SF1. Within the QUA2 domain, a transposition of adjacent arginine and lysine residues is primarily responsible for the switch in RNA binding between BBP and SF1. Alignment of multiple branchpoint binding proteins and the related STAR/GSG proteins suggests that the identity of these two amino acids and the RNA target sequences of all of these proteins are correlated.