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 ⁹¹QYNK⁹⁴ to ⁹¹PSFW⁹⁴ 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.     

Functional studies on the ATM intronic splicing processing element

In disease-associated genes, the understanding of the functional significance of deep intronic nucleotide variants may represent a difficult challenge. We have previously reported a new disease-causing mechanism that involves an intronic splicing processing element (ISPE) in ATM, composed of adjacent consensus 5′ and 3′ splice sites. A GTAA deletion within ISPE maintains potential adjacent splice sites, disrupts a non-canonical U1 snRNP interaction and activates an aberrant exon. In this paper, we demonstrate that binding of U1 snRNA through complementarity within a ~40 nt window downstream of the ISPE prevents aberrant splicing. By selective mutagenesis at the adjacent consensus ISPE splice sites, we show that this effect is not due to a resplicing process occurring at the ISPE. Functional comparison of the ATM mouse counterpart and evaluation of the pre-mRNA splicing intermediates derived from affected cell lines and hybrid minigene assays indicate that U1 snRNP binding at the ISPE interferes with the cryptic acceptor site. Activation of this site results in a stringent 5′–3′ order of intron sequence removal around the cryptic exon. Artificial U1 snRNA loading by complementarity to heterologous exonic sequences represents a potential therapeutic method to prevent the usage of an aberrant CFTR cryptic exon. Our results suggest that ISPE-like intronic elements binding U1 snRNPs may regulate correct intron processing.     

In vitro and in silico analysis reveals an efficient algorithm to predict the splicing consequences of mutations at the 5' splice sites

We have found that two previously reported exonic mutations in the PINK1 and PARK7 genes affect pre-mRNA splicing. To develop an algorithm to predict underestimated splicing consequences of exonic mutations at the 5′ splice site, we constructed and analyzed 31 minigenes carrying exonic splicing mutations and their derivatives. We also examined 189 249 U2-dependent 5′ splice sites of the entire human genome and found that a new variable, the SD-Score, which represents a common logarithm of the frequency of a specific 5′ splice site, efficiently predicts the splicing consequences of these minigenes. We also employed the information contents (R_i) to improve the prediction accuracy. We validated our algorithm by analyzing 32 additional minigenes as well as 179 previously reported splicing mutations. The SD-Score algorithm predicted aberrant splicings in 198 of 204 sites (sensitivity = 97.1%) and normal splicings in 36 of 38 sites (specificity = 94.7%). Simulation of all possible exonic mutations at positions −3, −2 and −1 of the 189 249 sites predicts that 37.8, 88.8 and 96.8% of these mutations would affect pre-mRNA splicing, respectively. We propose that the SD-Score algorithm is a practical tool to predict splicing consequences of mutations affecting the 5′ splice site.     

An exonic splicing silencer in the testes-specific DNA ligase III beta exon

Alternative pre-mRNA splicing of two terminal exons (α and β) regulates the expression of the human DNA ligase III gene. In most tissues, the α exon is expressed. In testes and during spermatogenesis, the β exon is used instead. The α exon encodes the interaction domain with a scaffold DNA repair protein, XRCC1, while the β exon-encoded C-terminal does not. Sequence elements regulating the alternative splicing pattern were mapped by in vitro splicing assays in HeLa nuclear extracts. Deletion of a region beginning in the β exon and extending into the downstream intron derepressed splicing to the β exon. Two silencing elements were found within this 101 nt region: a 16 nt exonic splicing silencer immediately upstream of the β exon polyadenylation signal and a 45 nt intronic splicing silencer. The exonic splicing silencer inhibited splicing, even when the polyadenylation signal was deleted or replaced by a 5′ splice site. This element also enhanced polyadenylation under conditions unfavourable to splicing. The splicing silencer partially inhibited assembly of spliceosomal complexes and functioned in an adenoviral pre-mRNA context. Silencing of splicing by the element was associated with cross-linking of a 37 kDa protein to the RNA substrate. The element exerts opposite functions in splicing and polyadenylation.     

Evidence that RNA editing modulates splice site selection in the 5-HT2C receptor gene

Adenosine to inosine editing of mRNA from the human 5-HT2C receptor gene (HTR2C) occurs at five exonic positions (A–E) in a stable stem–loop that includes the normal 5′ splice site of intron 5 and is flanked by two alternative splice sites. Using in vitro editing, we identified a novel editing site (F) located in the intronic part of the stem–loop and demonstrated editing at this site in human brain. We have shown that in cell culture, base substitutions to mimic editing at different combinations of the six sites profoundly affect relative splicing at the normal and the upstream alternative splice site, but splicing at the downstream alternative splice site was consistently rare. Editing combinations in different splice variants from human brain were determined and are consistent with the effects of editing on splicing observed in cell culture. As RNA editing usually occurs close to exon/intron boundaries, this is likely to be a general phenomenon and suggests an important novel role for RNA editing.     

Function of a Bovine Papillomavirus Type 1 Exonic Splicing Suppressor Requires a Suboptimal Upstream 3′ Splice Site

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. Previous studies in our laboratory have identified two purine-rich exonic splicing enhancers (ESEs), SE1 and SE2, located between two alternative 3′ splice sites at nucleotide (nt) 3225 and nt 3605. Further analysis of BPV-1 late-pre-mRNA splicing in vitro revealed a 48-nt pyrimidine-rich region immediately downstream of SE1 that inhibits utilization of the nt 3225 3′ splice site. This inhibitory element, which we named an exonic splicing suppressor (ESS), has a U-rich 5′ end, a C-rich central part, and an AG-rich 3′ end (Z. M. Zheng, P. He, and C. C. Baker, J. Virol. 70:4691–4699, 1996). The present study utilized in vitro splicing of both homologous and heterologous pre-mRNAs to further characterize the ESS. The BPV-1 ESS was inserted downstream of the 3′ splice site in the BPV-1 late pre-mRNA, Rous sarcoma virus src pre-mRNA, human immunodeficiency virus tat-rev pre-mRNA, and Drosophila dsx pre-mRNA, all containing a suboptimal 3′ splice site, and in the human β-globin pre-mRNA, which contains a constitutive 3′ splice site. These studies demonstrated that suppression of splicing by the BPV-1 ESS requires an upstream suboptimal 3′ splice site but not an upstream ESE. Furthermore, the ESS functions when located either upstream or downstream of BPV-1 SE1. Mutational analyses demonstrated that the function of the ESS is sequence dependent and that only the C-rich region of the ESS is essential for suppression of splicing in all the pre-mRNAs tested.     

Functionally antagonistic sequences are required for normal autoregulation of Drosophila tra-2 pre-mRNA splicing

Expression of functional TRA-2 protein in the male germline of Drosophila is regulated through a negative feedback mechanism in which a specific TRA-2 isoform represses splicing of the M1 intron in the TRA-2 pre-mRNA. We have previously shown that the mechanism of M1 splicing repression is conserved between distantly related Drosophila species. Using transgenic fly strains, we have examined the effects on regulation of mutations in two conserved features of the M1 intron. Our results show that TRA-2-dependent repression of M1 splicing depends on the presence of a suboptimal non-consensus 3′ splice site. Substitution of this 3′ splice site with a strong splice site resulted in TRA-2 independent splicing, while substitution with an unrelated weak 3′ splice site was compatible with repression, implying that reduced basal splicing efficiency is important for regulation. A second conserved element internal to the intron was found to be essential for efficient M1 splicing in the soma where the intron is not normally retained. We show that the role of this element is to enhance splicing and overcome the reduction in efficiency caused by the intron’s suboptimal 3′ splice site. Our results indicate that antagonistic elements in the M1 intron act together to establish a context that is permissive for repression of splicing by TRA-2 while allowing efficient splicing in the absence of a repressor.     

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.     

The Pre-mRNA 5′ Cap Determines Whether U6 Small Nuclear RNA Succeeds U1 Small Nuclear Ribonucleoprotein Particle at 5′ Splice Sites

Efficient splicing of the 5′-most intron of pre-mRNA requires a 5′ m⁷G(5′)ppp(5′)N cap, which has been implicated in U1 snRNP binding to 5′ splice sites. We demonstrate that the cap alters the kinetic profile of U1 snRNP binding, but its major effect is on U6 snRNA binding. With two alternative wild-type splice sites in an adenovirus pre-mRNA, the cap selectively alters U1 snRNA binding at the site to which cap-independent U1 snRNP binding is stronger and that is used predominantly in splicing; with two consensus sites, the cap acts on both, even though one is substantially preferred for splicing. However, the most striking quantitative effect of the 5′ cap is neither on U1 snRNP binding nor on the assembly of large complexes but on the replacement of U1 snRNP by U6 snRNA at the 5′ splice site. Inhibition of splicing by a cap analogue is correlated with the loss of U6 interactions at the 5′ splice site and not with any loss of U1 snRNP binding.     

A genetic screen in C. elegans reveals roles for KIN17 and PRCC in maintaining 5’ splice site identity

Pre-mRNA splicing is an essential step of eukaryotic gene expression carried out by a series of dynamic macromolecular protein/RNA complexes, known collectively and individually as the spliceosome. This series of spliceosomal complexes define, assemble on, and catalyze the removal of introns. Molecular model snapshots of intermediates in the process have been created from cryo-EM data, however, many aspects of the dynamic changes that occur in the spliceosome are not fully understood. Caenorhabditis elegans follow the GU-AG rule of splicing, with almost all introns beginning with 5’ GU and ending with 3’ AG. These splice sites are identified early in the splicing cycle, but as the cycle progresses and “custody” of the pre-mRNA splice sites is passed from factor to factor as the catalytic site is built, the mechanism by which splice site identity is maintained or re-established through these dynamic changes is unclear. We performed a genetic screen in C. elegans for factors that are capable of changing 5’ splice site choice. We report that KIN17 and PRCC are involved in splice site choice, the first functional splicing role proposed for either of these proteins. Previously identified suppressors of cryptic 5’ splicing promote distal cryptic GU splice sites, however, mutations in KIN17 and PRCC instead promote usage of an unusual proximal 5’ splice site which defines an intron beginning with UU, separated by 1nt from a GU donor. We performed high-throughput mRNA sequencing analysis and found that mutations in PRCC, and to a lesser extent KIN17, changed alternative 5’ splice site usage at native sites genome-wide, often promoting usage of nearby non-consensus sites. Our work has uncovered both fine and coarse mechanisms by which the spliceosome maintains splice site identity during the complex assembly process.     

Regulation of glutamate receptor B pre-mRNA splicing by RNA editing

RNA-editing enzymes of the ADAR family convert adenosines to inosines in double-stranded RNA substrates. Frequently, editing sites are defined by base-pairing of the editing site with a complementary intronic region. The glutamate receptor subunit B (GluR-B) pre-mRNA harbors two such exonic editing sites termed Q/R and R/G. Data from ADAR knockout mice and in vitro editing assays suggest an intimate connection between editing and splicing of GluR-B pre-mRNA.

By comparing the events at the Q/R and R/G sites, we can show that editing can both stimulate and repress splicing efficiency. The edited nucleotide, but not ADAR binding itself, is sufficient to exert this effect. The presence of an edited nucleotide at the R/G site reduces splicing efficiency of the adjacent intron facilitating alternative splicing events occurring downstream of the R/G site.

Lack of editing inhibits splicing at the Q/R site. Editing of both the Q/R nucleotide and an intronic editing hotspot are required to allow efficient splicing. Inefficient intron removal may ensure that only properly edited mRNAs become spliced and exported to the cytoplasm.

     

Human ribosomal protein S13 regulates expression of its own gene at the splicing step by a feedback mechanism

The expression of ribosomal protein (rp) genes is regulated at multiple levels. In yeast, two genes are autoregulated by feedback effects of the protein on pre-mRNA splicing. Here, we have investigated whether similar mechanisms occur in eukaryotes with more complicated and highly regulated splicing patterns. Comparisons of the sequences of ribosomal protein S13 gene (RPS13) among mammals and birds revealed that intron 1 is more conserved than the other introns. Transfection of HEK 293 cells with a minigene-expressing ribosomal protein S13 showed that the presence of intron 1 reduced expression by a factor of four. Ribosomal protein S13 was found to inhibit excision of intron 1 from rpS13 pre-mRNA fragment in vitro. This protein was shown to be able to specifically bind the fragment and to confer protection against ribonuclease cleavage at sequences near the 5′ and 3′ splice sites. The results suggest that overproduction of rpS13 in mammalian cells interferes with splicing of its own pre-mRNA by a feedback mechanism.     

Temperature-dependent splicing of beta-globin pre-mRNA

A T→G mutation at nucleotide 705 of human β-globin intron 2 creates an aberrant 5′ splice site and activates a cryptic 3′ splice site upstream. In consequence, the pre-mRNA is spliced via aberrant splice sites, despite the presence of the still functional correct sites. Surprisingly, when IVS2-705 HeLa or K562 cells were cultured at temperatures below 30°C, aberrant splicing was inhibited and correct splicing was restored. Similar temperature effects were seen for another β-globin pre-mRNA, IVS2-745, and in a construct in which a β-globin intron was inserted into a coding sequence of EGFP. Temperature-induced alternative splicing was affected by the nature of the internal aberrant splice sites flanking the correct sites and by exonic sequences. The results indicate that in the context of thalassemic splicing mutations and possibly in other alternatively spliced pre-mRNAs, temperature is one of the parameters that affect splice site selection.     

Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

The frequency distribution of mutation-induced aberrant 3′ splice sites (3′ss) in exons and introns is more complex than for 5′ splice sites, largely owing to sequence constraints upstream of intron/exon boundaries. As a result, prediction of their localization remains a challenging task. Here, nucleotide sequences of previously reported 218 aberrant 3′ss activated by disease-causing mutations in 131 human genes were compared with their authentic counterparts using currently available splice site prediction tools. Each tested algorithm distinguished authentic 3′ss from cryptic sites more effectively than from de novo sites. The best discrimination between aberrant and authentic 3′ss was achieved by the maximum entropy model. Almost one half of aberrant 3′ss was activated by AG-creating mutations and ~95% of the newly created AGs were selected in vivo. The overall nucleotide structure upstream of aberrant 3′ss was characterized by higher purine content than for authentic sites, particularly in position −3, that may be compensated by more stringent requirements for positive and negative nucleotide signatures centred around position −11. A newly developed online database of aberrant 3′ss will facilitate identification of splicing mutations in a gene or phenotype of interest and future optimization of splice site prediction tools.     

An Alu-derived intronic splicing enhancer facilitates intronic processing and modulates aberrant splicing in ATM

We have previously reported a natural GTAA deletion within an intronic splicing processing element (ISPE) of the ataxia telangiectasia mutated (ATM) gene that disrupts a non-canonical U1 snRNP interaction and activates the excision of the upstream portion of the intron. The resulting pre-mRNA splicing intermediate is then processed to a cryptic exon, whose aberrant inclusion in the final mRNA is responsible for ataxia telangiectasia. We show here that the last 40 bases of a downstream intronic antisense Alu repeat are required for the activation of the cryptic exon by the ISPE deletion. Evaluation of the pre-mRNA splicing intermediate by a hybrid minigene assay indicates that the identified intronic splicing enhancer represents a novel class of enhancers that facilitates processing of splicing intermediates possibly by recruiting U1 snRNP to defective donor sites. In the absence of this element, the splicing intermediate accumulates and is not further processed to generate the cryptic exon. Our results indicate that Alu-derived sequences can provide intronic splicing regulatory elements that facilitate pre-mRNA processing and potentially affect the severity of disease-causing splicing mutations.     

The Peculiarities of Large Intron Splicing in Animals

In mammals a considerable 92% of genes contain introns, with hundreds and hundreds of these introns reaching the incredible size of over 50,000 nucleotides. These “large introns” must be spliced out of the pre-mRNA in a timely fashion, which involves bringing together distant 5′ and 3′ acceptor and donor splice sites. In invertebrates, especially Drosophila, it has been shown that larger introns can be spliced efficiently through a process known as recursive splicing—a consecutive splicing from the 5′-end at a series of combined donor-acceptor splice sites called RP-sites. Using a computational analysis of the genomic sequences, we show that vertebrates lack the proper enrichment of RP-sites in their large introns, and, therefore, require some other method to aid splicing. We analyzed over 15,000 non-redundant, large introns from six mammals, 1,600 from chicken and zebrafish, and 560 non-redundant large introns from five invertebrates. Our bioinformatic investigation demonstrates that, unlike the studied invertebrates, the studied vertebrate genomes contain consistently abundant amounts of direct and complementary strand interspersed repetitive elements (mainly SINEs and LINEs) that may form stems with each other in large introns. This examination showed that predicted stems are indeed abundant and stable in the large introns of mammals. We hypothesize that such stems with long loops within large introns allow intron splice sites to find each other more quickly by folding the intronic RNA upon itself at smaller intervals and, thus, reducing the distance between donor and acceptor sites.