Branch Point Sequence Research Articles

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene

BackgroundThe BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates.ResultsThe results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3′ splice site (AG) and 5′ splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination.ConclusionsWe show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

Disclosing the Impact of Carcinogenic SF3b Mutations on Pre-mRNA Recognition Via All-Atom Simulations.

The spliceosome accurately promotes precursor messenger-RNA splicing by recognizing specific noncoding intronic tracts including the branch point sequence (BPS) and the 3’-splice-site (3’SS). Mutations of Hsh155 (yeast)/SF3B1 (human), which is a protein of the SF3b factor involved in BPS recognition and induces altered BPS binding and 3’SS selection, lead to mis-spliced mRNA transcripts. Although these mutations recur in hematologic malignancies, the mechanism by which they change gene expression remains unclear. In this study, multi-microsecond-long molecular-dynamics simulations of eighth distinct ∼700,000 atom models of the spliceosome Bact complex, and gene sequencing of SF3B1, disclose that these carcinogenic isoforms destabilize intron binding and/or affect the functional dynamics of Hsh155/SF3B1 only when binding non-consensus BPSs, as opposed to the non-pathogenic variants newly annotated here. This pinpoints a cross-talk between the distal Hsh155 mutation and BPS recognition sites. Our outcomes unprecedentedly contribute to elucidating the principles of pre-mRNA recognition, which provides critical insights on the mechanism underlying constitutive/alternative/aberrant splicing.

Activation of Cryptic 3' Splice-Sites by SRSF2 Contributes to Cassette Exon Skipping.

Here we show that the serine/arginine rich splicing factor 2 (SRSF2) promotes cryptic 3′ splice-site (3′AG′) usage during cassette exon exclusion in survival of motor neuron (SMN2) minigenes. Deletion of the 3′AG′ (3′AG′1), its associated branch point (BP′) and polypyrimidine tract (PPT′) sequences directs SRSF2 to promote a second 3′AG′ (3′AG′2) with less conserved associated region for intron splicing. Furthermore, deletion of both 3′AG′1 and 3′AG′2 and their associated sequences triggered usage of a third 3′AG′3 that has very weak associated sequences. Interestingly, when intron splicing was directed to the 3′AG′ cryptic splice-sites, intron splicing from the canonical 3′AG splice-site was reduced along with a decrease in cassette exon inclusion. Moreover, multiple SRSF2 binding sites within the intron are responsible for 3′AG′ activation. We conclude that SRSF2 facilitates exon exclusion by activating a cryptic 3′AG′ and inhibiting downstream intron splicing.

HPV18 Utilizes Two Alternative Branch Sites for E6*I Splicing to Produce E7 Protein

Human papillomavirus 18 (HPV18) E6 and E7 oncogenes are transcribed as a single bicistronic E6E7 pre-mRNA. The E6 ORF region in the bicistronic E6E7 pre-mRNA contains an intron. Splicing of this intron disrupts the E6 ORF integrity and produces a spliced E6*I RNA for efficient E7 translation. Here we report that the E6 intron has two overlapped branch point sequences (BPS) upstream of its 3′ splice site, with an identical heptamer AACUAAC, for E6*I splicing. One heptamer has a branch site adenosine (underlined) at nt 384 and the other at nt 388. E6*I splicing efficiency correlates to the expression level of E6 and E7 proteins and depends on the selection of which branch site. In general, E6*I splicing prefers the 3′ss-proximal branch site at nt 388 over the distal branch site at nt 384. Inactivation of the nt 388 branch site was found to activate a cryptic acceptor site at nt 636 for aberrant RNA splicing. Together, these data suggest that HPV18 modulates its production ratio of E6 and E7 proteins by alternative selection of the two mapped branch sites for the E6*I splicing, which could be beneficial in its productive or oncogenic infection according to the host cell environment.

Identification of Genetic Interactions Associated with Splicing Factor Prp5

Regulation of gene expression occurs at multiple levels, including the removal of non‐coding sequences from protein‐coding pre‐messenger RNA transcripts, commonly known as pre‐mRNA splicing. Correct gene expression requires high fidelity of pre‐mRNA splicing, as up to 15% of mutations leading to disease may be a result of aberrant pre‐mRNA splicing, including cystic fibrosis, spinal muscular atrophy (SMA) and some leukemias. Pre‐mRNA splicing is a multistep process that requires several proteins and RNAs for ATP hydrolysis, including multiple essential DEAD‐box proteins such as Prp5. Here we investigate the role of the DEAD‐box protein Prp5 in splicing using the yeast Saccharomyces cerevisiae as a model. Prp5 has previously been shown to be required for ATP‐dependent RNA remodeling in splicing, facilitating binding of the U2 snRNA to the branchpoint sequence in the pre‐mRNA intron. In addition, unpublished research has revealed genetic interactions between Prp5 and non‐splicing genes such as Gcn5 and Spt8. In this study, additional potential genetic interactions involving Prp5 were investigated. To identify these interactions, a transposon mutagenesis screen was employed. Approximately 20,000 mutations were screened, 61 of which showed putative genetic suppression of the temperature sensitive prp5‐1 mutation. Subsequent growth assays confirmed that many of these mutations exhibited genetic suppression, indicating important genetic interactions between these genes and Prp5. hiTAIL PCR was utilized to identify a number of these genetic interactions, including an interaction between Prp5 and the MMO1 gene, a verified protein in S. cerevisiae with unknown function. Additional experiments to characterize this interaction with MMO1 are currently being performed. Furthermore, we have used high‐throughput sequencing to identify an additional subset of the interactions from the transposon screen. This study will help shed light upon the role of the DEAD‐box protein Prp5 in pre‐mRNA splicing.This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Structures of the Catalytically Activated Yeast Spliceosome Reveal the Mechanism of Branching

Pre-mRNA splicing is executed by the spliceosome. Structural characterization of the catalytically activated complex (B∗) is pivotal for understanding thebranching reaction. In this study, we assembled theB∗ complexes on two different pre-mRNAs from Saccharomyces cerevisiae and determined the cryo-EM structures of four distinct B∗ complexes at overall resolutions of 2.9-3.8Å. The duplex between U2 small nuclear RNA (snRNA) and the branch pointsequence (BPS) is discretely away from the 5'-splice site (5'SS) in the three B∗ complexes that are devoid of the step I splicing factors Yju2 and Cwc25. Recruitment of Yju2 into the active site brings the U2/BPS duplex into the vicinity of 5'SS, with theBPS nucleophile positioned 4Å away from the catalytic metal M2. This analysis reveals the functional mechanism of Yju2 and Cwc25 in branching. These structures on different pre-mRNAs reveal substrate-specific conformations of the spliceosome in a major functional state.

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.

Structures of the human spliceosomes before and after release of the ligated exon

Pre-mRNA splicing is executed by the spliceosome, which has eight major functional states each with distinct composition. Five of these eight human spliceosomal complexes, all preceding exon ligation, have been structurally characterized. In this study, we report the cryo-electron microscopy structures of the human post-catalytic spliceosome (P complex) and intron lariat spliceosome (ILS) at average resolutions of 3.0 and 2.9 Å, respectively. In the P complex, the ligated exon remains anchored to loop I of U5 small nuclear RNA, and the 3′-splice site is recognized by the junction between the 5′-splice site and the branch point sequence. The ATPase/helicase Prp22, along with the ligated exon and eight other proteins, are dissociated in the P-to-ILS transition. Intriguingly, the ILS complex exists in two distinct conformations, one with the ATPase/helicase Prp43 and one without. Comparison of these three late-stage human spliceosomes reveals mechanistic insights into exon release and spliceosome disassembly.

Degenerate minigene library analysis enables identification of altered branch point utilization by mutant splicing factor 3B1 (SF3B1).

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3′ splice sites (3′SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3′SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3′ intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1K700E and SF3B1WT. Using this model, we show that SF3B1K700E spliceosomes utilize non-canonical sequence variants (at position −1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

ShapeShifter: a novel approach for identifying and quantifying stable lariat intronic species in RNAseq data.

Most intronic lariats are rapidly turned over after splicing. However, new research suggests that some introns may have additional post-splicing functions. Current bioinformatics methods used to identify lariats require a sequencing read that traverses the lariat branchpoint. This method provides precise branchpoint sequence and position information, but is limited in its ability to quantify abundance of stabilized lariat species in a given RNAseq sample. Bioinformatic tools are needed to better address these emerging biological questions. We used an unsupervised machine learning approach on sequencing reads from publicly available ENCODE data to learn to identify and quantify lariats based on RNAseq read coverage shape. We developed ShapeShifter, a novel approach for identifying and quantifying stable lariat species in RNAseq datasets. We learned a characteristic "lariat" curve from ENCODE RNAseq data and were able to estimate abundances for introns based on read coverage. Using this method we discovered new stable introns in these samples that were not represented using the older, branchpoint-traversing read method. ShapeShifter provides a robust approach towards detecting and quantifying stable lariat species.

Prespliceosome structure provides insights into spliceosome assembly and regulation

The spliceosome catalyzes the excision of introns from pre-mRNA in two steps, branching and exon ligation, and is assembled from five small nuclear ribonucleoprotein particles (snRNPs; U1, U2, U4, U5, U6) and numerous non-snRNP factors1. For branching, the intron 5'-splice site (5'SS) and the branch point (BP) sequence are selected and brought into the prespliceosome by the U1 and U2 snRNPs1, which is a focal point for the regulation by alternative splicing factors2. The U4/U6.U5 tri-snRNP subsequently joins the prespliceosome to form the complete pre-catalytic spliceosome. Recent studies have revealed the structural basis of the branching and exon-ligation reactions3. However, the structural basis of early spliceosome assembly events remains poorly understood4. Here we report the cryo-electron microscopy structure of the yeast Saccharomyces cerevisiae prespliceosome at near-atomic resolution. The structure reveals an induced stabilization of the 5'SS in the U1 snRNP, and provides structural insights into the functions of the human alternative splicing factors LUC7-like (yeast Luc7) and TIA-1 (yeast Nam8) that are linked to human disease5,6. In the prespliceosome, the U1 snRNP associates with the U2 snRNP through a stable contact with the U2 3' domain and a transient yeast-specific contact with the U2 SF3b-containing 5' region, leaving its tri-snRNP-binding interface fully exposed. The results suggest mechanisms for 5'SS transfer to the U6 ACAGAGA region within the assembled spliceosome and for its subsequent conversion to the activation-competent B complex spliceosome7,8. Taken together, the data provide a working model to investigate the early steps of spliceosome assembly.

Structure of the Post-catalytic Spliceosome from Saccharomyces cerevisiae

Removal of an intron from a pre-mRNA by the spliceosome results in the ligation of two exons in the post-catalytic spliceosome (known as the P complex). Here, we present a cryo-EM structure of the Pcomplex from Saccharomyces cerevisiae at an average resolution of 3.6Å. The ligated exon is held in the active site through RNA-RNA contacts. Three bases at the 3' end of the 5' exon remain anchored to loop I of U5 small nuclear RNA, and the conserved AG nucleotides of the 3'-splice site (3'SS) are specifically recognized by the invariant adenine of the branch point sequence, the guanine base at the 5'end of the 5'SS, and an adenine base of U6 snRNA. The 3'SS is stabilized through an interaction with the1585-loop of Prp8. The P complex structure provides a view on splice junction formation critical for understanding the complete splicing cycle.

Rules and tools to predict the splicing effects of exonic and intronic mutations.

Development of next generation sequencing technologies has enabled detection of extensive arrays of germline and somatic single nucleotide variations (SNVs) in human diseases. SNVs affecting intronic GT-AG dinucleotides invariably compromise pre-mRNA splicing. Most exonic SNVs introduce missense/nonsense codons, but some affect auxiliary splicing cis-elements or generate cryptic GT-AG dinucleotides. Similarly, most intronic SNVs are silent, but some affect canonical and auxiliary splicing cis-elements or generate cryptic GT-AG dinucleotides. However, prediction of the splicing effects of SNVs is challenging. The splicing effects of SNVs generating cryptic AG or disrupting canonical AG can be inferred from the AG-scanning model. Similarly, the splicing effects of SNVs affecting the first nucleotide G of an exon can be inferred from AG-dependence of the 3' splice site (ss). A variety of tools have been developed for predicting the splicing effects of SNVs affecting the 5' ss, as well as exonic and intronic splicing enhancers/silencers. In contrast, only two tools, the Human Splicing Finder and the SVM-BP finder, are available for predicting the position of the branch point sequence. Similarly, IntSplice and Splicing based Analysis of Variants (SPANR) are the only tools to predict the splicing effects of intronic SNVs. The rules and tools introduced in this review are mostly based on observations of a limited number of genes, and no rule or tool can ensure 100% accuracy. Experimental validation is always required before any clinically relevant conclusions are drawn. Development of efficient tools to predict aberrant splicing, however, will facilitate our understanding of splicing pathomechanisms in human diseases. WIREs RNA 2018, 9:e1451. doi: 10.1002/wrna.1451 This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease RNA Methods > RNA Analyses In Vitro and In Silico.

Error-Prone Splicing Controlled by the Ubiquitin Relative Hub1

Accurate pre-mRNA splicing is needed for correct gene expression and relies on faithful splice site recognition. Here, we show that the ubiquitin-like protein Hub1 binds to the DEAD-box helicase Prp5, a key regulator of early spliceosome assembly, and stimulates its ATPase activity thereby enhancing splicing and relaxing fidelity. High Hub1 levels enhance splicing efficiency but also cause missplicing by tolerating suboptimal splice sites and branchpoint sequences. Notably, Prp5 itself is regulated by a Hub1-dependent negative feedback loop. Since Hub1-mediated splicing activation induces cryptic splicing of Prp5, it also represses Prp5 protein levels and thus curbs excessive missplicing. Our findings indicate that Hub1 mediates enhanced, but error-prone splicing, a mechanism that is tightly controlled by a feedback loop of PRP5 cryptic splicing activation.

Rare splicing defects of FAS underly severe recessive autoimmune lymphoproliferative syndrome

Autoimmune lymphoproliferative syndrome (ALPS) is a prototypic disorder of impaired apoptosis characterized by autoimmune features and lymphoproliferation. Heterozygous germline or somatic FAS mutations associated with preserved protein expression have been described. Very rare cases of homozygous germline FAS mutations causing severe autosomal recessive form of ALPS with a complete defect of Fas expression have been reported.We report two unrelated patients from highly inbred North African population showing a severe ALPS phenotype and an undetectable Fas surface expression. Two novel homozygous mutations have been identified underlying rare splicing defects mechanisms. The first mutation breaks a branch point sequence and the second alters a regulatory exonic splicing site. These splicing defects induce the skipping of exon 6 encoding the transmembrane domain of CD95. Our findings highlight the requirement of tight regulation of FAS exon 6 splicing for balanced alternative splicing and illustrate the importance of such studies in highly consanguineous populations.

BPP: a sequence-based algorithm for branch point prediction.

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. 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. djguo@cuhk.edu.hk. Supplementary data are available at Bioinformatics online.

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

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

A homozygous intronic branch-point deletion in the ALPL gene causes infantile hypophosphatasia

Hypophosphatasia (HPP) is a multi-systemic inborn disease with an extraordinary spectrum of severity, ranging from the absence of mineralization to high lethality and it involves different organs including bone, muscle, kidney, lung, gastrointestinal tract and the nervous system. The disease is characterized by low levels of serum alkaline phosphatase, caused by loss-of-function mutations within the ALPL gene that encodes the tissue-nonspecific alkaline phosphatase TNAP. Here we present the functional characterization of a gene mutation, detected in intron 7 of the ALPL gene of a boy with infantile HPP in whom routine sequencing of the coding region failed to detect any mutation. The homozygous c.793del-14_33 mutation results in the loss of the branch-point motif, relevant for correct ALPL pre-mRNA splicing. The main transcript skips exon 8 and codes for a C-terminally truncated TNAP protein of 275 amino acids, which was detected in peripheral blood mononuclear cells and serum from the patient. The functional characterization of recombinant TNAP275 revealed no enzymatic activity nor any dominant-negative effect, relevant for the heterozygous parents. Nevertheless correct pre-mRNA splicing can take place without the branch-point sequence to a limited extend, as concluded from the ALPL cDNA, obtained from patient's PBMC, and from the low serum AP activity.These data reaffirm that in clear cut clinical cases, where conventional sequencing including the coding sequence and direct exon-intron-boundaries fails to detect mutations, deeper analyses of regulatory important motifs like branch-point sequences are required to establish a genetic diagnosis.

A recurrent deletion in the SLC5A2 gene including the intron 7 branch site responsible for familial renal glucosuria.

Familial renal glycosuria (FRG) is caused by mutations in the SLC5A2 gene, which codes for Na+-glucose co-transporters 2 (SGLT2). The aim of this study was to analyze and identify the mutations in 16 patients from 8 families with FRG. All coding regions, including intron-exon boundaries, were analyzed using PCR followed by direct sequence analysis. Six mutations in SLC5A2 gene were identified, including five missense mutations (c.393G > C, p.K131N; c.1003A > G, p.S335G; c.1343A > G, p.Q448R; c.1420G > C, p.A474P; c.1739G > A, p.G580D) and a 22-bp deletion in intron 7 (c.886(-10_-31)del) removing the putative branch point sequence. By the minigene studies using the pSPL3 plasmids, we confirmed that the deletion c.886(-10_-31)del acts as a splicing mutation. Furthermore, we found that this deletion causes exclusion of exon 8 in the SCL5A2 transcript in patients. The mutation c.886(-10_-31)del was present in 5 (62.5%) of 8 families, and accounts for about 37.5% of the total alleles (6/16). In conclusion, six mutations resulting in FRG were found, and the c.886(-10_-31)del may be the high frequency mutation that can be screened in FRG patients with uniallelic or negative SLC5A2 mutations.

Reprogramming Msl5 to recognize mutant branch point sequences during spliceosome assembly

Saccharomyces cerevisiae Msl5 is an essential protein responsible for the assembly of a catalytically active spliceosome. Msl5 mediates spliceosome assembly through the recognition of the intronic branchpoint sequence 5′‐UACUAAC. Upon binding, Msl5 establishes an interaction with U1 bound to the 5′ splice site before recruiting the U2 snRNP to the intronic branchpoint sequence. We have generated and expressed mutant Msl5 proteins containing L176N and L167A/Q255A mutations, designed to bind specifically to mutant branch point sequences 5′‐UACUACC and 5′‐UACAAAC, respectively. Electrophoretic mobility shift assays will be performed to determine if reprogrammed Msl5 can specifically recognize and bind to a mutant branch point sequence. The reprogramming of Msl5 can be used to generate an orthogonal Msl5 system in vivo, and to analyze the effectiveness of intron splicing containing different branchpoint sequences.