基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序（RNA-seq）数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件（其中新识别的可变剪接事件3922个）,包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.     

Wilms’ tumor 1-associating protein complex regulates alternative splicing and polyadenylation at potential G-quadruplex-forming splice site sequences

Wilms’ tumor 1-associating protein (WTAP) is a core component of the N6-methyladenosine (m6A)-methyltransferase complex, along with VIRMA, CBLL1, ZC3H13 (KIAA0853), RBM15/15B, and METTL3/14, which generate m6A, a key RNA modification that affects various processes of RNA metabolism. WTAP also interacts with splicing factors; however, despite strong evidence suggesting a role of Drosophila WTAP homolog fl(2)d in alternative splicing (AS), its role in splicing regulation in mammalian cells remains elusive. Here we demonstrate using RNAi coupled with RNA-seq that WTAP, VIRMA, CBLL1, and ZC3H13 modulate AS, promoting exon skipping and intron retention in AS events that involve short introns/exons with higher GC content and introns with weaker polypyrimidine-tract and branch points. Further analysis of GC-rich sequences involved in AS events regulated by WTAP, together with minigene assay analysis, revealed potential G-quadruplex formation at splice sites where WTAP has an inhibitory effect. We also found that several AS events occur in the last exon of one isoform of MSL1 and WTAP, leading to competition for polyadenylation. Proteomic analysis also suggested that WTAP/CBLL1 interaction promotes recruitment of the 3′-end processing complex. Taken together, our results indicate that the WTAP complex regulates AS and alternative polyadenylation via inhibitory mechanisms in GC-rich sequences.     

Defective splicing induced by 4NQO in the hamster hprt gene

The molecular analysis of mutations affecting mRNA processing may contribute to a better understanding of the splicing mechanism through the identification of genomic sequences necessary for the recognition of splice sites. In this paper we report the sequence analysis of 14 splice mutants induced by 4-nitroquinoline 1-oxide (4NQO) at the hamster hypoxanthine-guanine-phosphoribosyltransferase (hprt) locus. We show that mutations at the 3′ acceptor splice site or at the first or fifth base of the 5′ donor splice site are responsible for exon skipping. In addition, mutations in exon sequences also determine the skipping of one or more exons. Our data indicate that point mutations in intron regions at either side of an internal exon may induce the skipping of the same exon, supporting a model where the exon is the unit of early spliceosome assembly. Furthermore, they suggest that the splicing of hprt mRNA precursors may proceed through a clustering of exons 2, 3 and 4 which are then spliced in a concerted way.     

Autoregulation of Fox protein expression to produce dominant negative splicing factors

The Fox proteins are a family of regulators that control the alternative splicing of many exons in neurons, muscle, and other tissues. Each of the three mammalian paralogs, Fox-1 (A2BP1), Fox-2 (RBM9), and Fox-3 (HRNBP3), produces proteins with a single RNA-binding domain (RRM) flanked by N- and C-terminal domains that are highly diversified through the use of alternative promoters and alternative splicing patterns. These genes also express protein isoforms lacking the second half of the RRM (FoxΔRRM), due to the skipping of a highly conserved 93-nt exon. Fox binding elements overlap the splice sites of these exons in Fox-1 and Fox-2, and the Fox proteins themselves inhibit exon inclusion. Unlike other cases of splicing autoregulation by RNA-binding proteins, skipping the RRM exon creates an in-frame deletion in the mRNA to produce a stable protein. These FoxΔRRM isoforms expressed from cDNA exhibit highly reduced binding to RNA in vivo. However, we show that they can act as repressors of Fox-dependent splicing, presumably by competing with full-length Fox isoforms for interaction with other splicing factors. Interestingly, the Drosophila Fox homolog contains a nearly identical exon in its RRM domain that also has flanking Fox-binding sites. Thus, rather than autoregulation of splicing controlling the abundance of the regulator, the Fox proteins use a highly conserved mechanism of splicing autoregulation to control production of a dominant negative isoform.     

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Conserved sequence elements associated with exon skipping

One of the major forms of alternative splicing, which generates multiple mRNA isoforms differing in the precise combinations of their exon sequences, is exon skipping. While in constitutive splicing all exons are included, in the skipped pattern(s) one or more exons are skipped. The regulation of this process is still not well understood; so far, cis- regulatory elements (such as exonic splicing enhancers) were identified in individual cases. We therefore set to investigate the possibility that exon skipping is controlled by sequences in the adjacent introns. We employed a computer analysis on 54 sequences documented as undergoing exon skipping, and identified two motifs both in the upstream and downstream introns of the skipped exons. One motif is highly enriched in pyrimidines (mostly C residues), and the other motif is highly enriched in purines (mostly G residues). The two motifs differ from the known cis-elements present at the 5′ and 3′ splice site. Interestingly, the two motifs are complementary, and their relative positional order is conserved in the flanking introns. These suggest that base pairing interactions can underlie a mechanism that involves secondary structure to regulate exon skipping. Remarkably, the two motifs are conserved in mouse orthologous genes that undergo exon skipping.     

A Multi-Exon-Skipping Detection Assay Reveals Surprising Diversity of Splice Isoforms of Spinal Muscular Atrophy Genes

Humans have two near identical copies of Survival Motor Neuron gene: SMN1 and SMN2. Loss of SMN1 coupled with the predominant skipping of SMN2 exon 7 causes spinal muscular atrophy (SMA), a neurodegenerative disease. SMA patient cells devoid of SMN1 provide a powerful system to examine splicing pattern of various SMN2 exons. Until now, similar system to examine splicing of SMN1 exons was unavailable. We have recently screened several patient cell lines derived from various diseases, including SMA, Alzheimer’s disease, Parkinson’s disease and Batten disease. Here we report a Batten disease cell line that lacks functional SMN2, as an ideal system to examine pre-mRNA splicing of SMN1. We employ a multiple-exon-skipping detection assay (MESDA) to capture simultaneously skipping of multiple exons. Our results show surprising diversity of splice isoforms and reveal novel splicing events that include skipping of exon 4 and co-skipping of three adjacent exons of SMN. Contrary to the general belief, MESDA captured oxidative-stress induced skipping of SMN1 exon 5 in several cell types, including non-neuronal cells. We further demonstrate that the predominant SMN2 exon 7 skipping induced by oxidative stress is modulated by a combinatorial control that includes promoter sequence, endogenous context, and the weak splice sites. We also show that an 8-mer antisense oligonucleotide blocking a recently described GC-rich sequence prevents SMN2 exon 7 skipping under the conditions of oxidative stress. Our findings bring new insight into splicing regulation of an essential housekeeping gene linked to neurodegeneration and infant mortality.     

Exon definition may facilitate splice site selection in RNAs with multiple exons.

Interactions at the 3' end of the intron initiate spliceosome assembly and splice site selection in vertebrate pre-mRNAs. Multiple factors, including U1 small nuclear ribonucleoproteins (snRNPs), are involved in initial recognition at the 3' end of the intron. Experiments were designed to test the possibility that U1 snRNP interaction at the 3' end of the intron during early assembly functions to recognize and define the downstream exon and its resident 5' splice site. Splicing precursor RNAs constructed to have elongated second exons lacking 5' splice sites were deficient in spliceosome assembly and splicing activity in vitro. Similar substrates including a 5' splice site at the end of exon 2 assembled and spliced normally as long as the second exon was less than 300 nucleotides long. U2 snRNPs were required for protection of the 5' splice site terminating exon 2, suggesting direct communication during early assembly between factors binding the 3' and 5' splice sites bordering an exon. We suggest that exons are recognized and defined as units during early assembly by binding of factors to the 3' end of the intron, followed by a search for a downstream 5' splice site. In this view, only the presence of both a 3' and a 5' splice site in the correct orientation and within 300 nucleotides of one another will stable exon complexes be formed. Concerted recognition of exons may help explain the 300-nucleotide-length maximum of vertebrate internal exons, the mechanism whereby the splicing machinery ignores cryptic sites within introns, the mechanism whereby exon skipping is normally avoided, and the phenotypes of 5' splice site mutations that inhibit splicing of neighboring introns.     

Characterization of exon skipping mutants of the COP1 gene from Arabidopsis

The removal of introns from pre-mRNA requires accurate recognition and selection of the intron splice sites. Mutations which alter splice site selection and which lead to skipping of specific exons are indicative of intron/exon recognition mechanisms involving an exon definition process. In this paper, three independent mutants to the COP1 gene in Arabidopsis which show exon skipping were identified and the mutations which alter the normal splicing pattern were characterized. The mutation in cop1–1 was a G→A change 4 nt upstream from the 3′ splice site of intron 5, while the mutation in cop1–2 was a G→A at the first nucleotide of intron 6, abolishing the conserved G within the 5′ splice site consensus. The effect of these mutations was skipping of exon 6. The mutation in cop1–8 was G→A in the final nucleotide of intron 10 abolishing the conserved G within the 3′ splice site consensus and leading to skipping of exon 11. The splicing patterns surrounding exons 6 and 11 of COP1 in these three mutant lines of Arabidopsis provide evidence for exon definition mechanisms operating in plant splicing.     

During in vivo maturation of eukaryotic nuclear mRNA, splicing yields excised exon circles.

Circular splicing has already been described on nuclear pre-mRNA for certain splice sites far apart in the multi exonic ETS-1 gene and in the single 1.2 kb exon of the Sry locus. To date, it is unclear how splice site juxtaposition occurs in normal and circular splicing. The splice site selection of an internal exon is likely to involve pairing between splice sites across that exon. Based on this, we predict that, albeit at low frequency, internal exons yield circular RNA by splicing as an error-prone mechanism of exon juxtaposition or, perhaps more interestingly, as a regulated mechanism on alternative exons. To address this question, the circular exon formation was analyzed at three ETS-1 internal exons (one alternative spliced exon and two constitutive), in human cell line and blood cell samples. Here, we show by RT-PCR and sequencing that exon circular splicing occurs at the three individual exons that we examined. RNase protection experiments suggest that there is no correlation between exon circle expression and exon skipping.     

Atypical 5' splice sites cause CFTR exon 9 to be vulnerable to skipping

The molecular basis of the skipping of constitutive exons in many messenger RNAs is not fully understood. A well-studied example is exon 9 of the human cystic fibrosis transmembrane conductance regulator gene (CFTR), in which an abbreviated polypyrimidine tract between the branch point A and the 3' splice site is associated with increased exon skipping and disease. However, many exons, both in CFTR and in other genes and have short polypyrimidine tracts in their 3' splice sites, yet they are not skipped. Inspection of the 5' splice sites immediately up- and downstream of exon 9 revealed deviations from consensus sequence, so we hypothesized that this exon may be inherently vulnerable to skipping. To test this idea, we constructed a CFTR minigene and replicated exon 9 skipping associated with the length of the polypyrimidine tract upstream of exon 9. We then mutated the flanking 5' splice sites and determined the effect on exon skipping. Conversion of the upstream 5' splice site to consensus by replacing a pyrimidine at position +3 with a purine resulted in increased exon skipping. In contrast, conversion of the downstream 5' splice site to consensus by insertion of an adenine at position +4 resulted in a substantial reduction in exon 9 skipping, regardless of whether the upstream 5' splice site was consensus or not. These results suggested that the native downstream 5' splice site plays an important role in CFTR exon 9 skipping, a hypothesis that was supported by data from sheep and mouse genomes. Although CFTR exon 9 in sheep is preceded by a long polypyrimidine tract (Y(14)), it skips exon 9 in vivo and has a nonconsensus downstream 5' splice site identical to that in humans. On the other hand, CFTR exon 9 in mice is preceded by a short polypyrimidine tract (Y(5)) but is not skipped in vivo. Its downstream 5' splice site differs from that in humans by a 2-nt insertion, which, when introduced into the human CFTR minigene, abolished exon 9 skipping. Taken together, these observations place renewed emphasis on deviations at 5' splice sites in nucleotides other than the invariant GT, particularly when such changes are found in conjunction with other altered splicing sequences, such as a shortened polypyrimidine tract. Thus, careful inspection of entire 5' splice sites may identify constitutive exons that are vulnerable to skipping.     

Does distance matter? Variations in alternative 3′ splicing regulation

Alternative splicing constitutes a major mechanism creating protein diversity in humans. This diversity can result from the alternative skipping of entire exons or by alternative selection of the 5′ or 3′ splice sites that define the exon boundaries. In this study, we analyze the sequence and evolutionary characteristics of alternative 3′ splice sites conserved between human and mouse genomes for distances ranging from 3 to 100 nucleotides. We show that alternative splicing events can be distinguished from constitutive splicing by a combination of properties which vary depending on the distance between the splice sites. Among the unique features of alternative 3′ splice sites, we observed an unexpectedly high occurrence of events in which a polypyrimidine tract was found to overlap the upstream splice site. By applying a machine-learning approach, we show that we can successfully discriminate true alternative 3′ splice sites from constitutive 3′ splice sites. Finally, we propose that the unique features of the intron flanking alternative splice sites are indicative of a regulatory mechanism that is involved in splice site selection. We postulate that the process of splice site selection is influenced by the distance between the competitive splice sites.     

Modulation of exon skipping and inclusion by heterogeneous nuclear ribonucleoprotein A1 and pre-mRNA splicing factor SF2/ASF.

The essential splicing factor SF2/ASF and the heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) modulate alternative splicing in vitro of pre-mRNAs that contain 5' splice sites of comparable strengths competing for a common 3' splice site. Using natural and model pre-mRNAs, we have examined whether the ratio of SF2/ASF to hnRNP A1 also regulates other modes of alternative splicing in vitro. We found that an excess of SF2/ASF effectively prevents inappropriate exon skipping and also influences the selection of mutually exclusive tissue-specific exons in natural beta-tropomyosin pre-mRNA. In contrast, an excess of hnRNP A1 does not cause inappropriate exon skipping in natural constitutively or alternatively spliced pre-mRNAs. Although hnRNP A1 can promote alternative exon skipping, this effect is not universal and is dependent, e.g., on the size of the internal alternative exon and on the strength of the polypyrimidine tract in the preceding intron. With appropriate alternative exons, an excess of SF2/ASF promotes exon inclusion, whereas an excess of hnRNP A1 causes exon skipping. We propose that in some cases the ratio of SF2/ASF to hnRNP A1 may play a role in regulating alternative splicing by exon inclusion or skipping through the antagonistic effects of these proteins on alternative splice site selection.     

Intron Retention in the Alternatively Spliced Region of RON Results from Weak 3’ Splice Site Recognition

The RON gene encodes a tyrosine kinase receptor for macrophage-stimulating protein. A constitutively active isoform that arises by skipping of exon 11 is expressed in carcinomas and contributes to an invasive phenotype. However, a high proportion of the mRNA expressed from the endogenous gene, or from transfected minigenes, appears to retain introns 10 and 11. It is not known whether this represents specific repression or the presence of weak splicing signals. We have used chimeric pre-mRNAs spliced in vitro to investigate the reason for intron retention. A systematic test showed that, surprisingly, the exon sequences known to modulate exon 11 skipping were not limiting, but the 3’ splice site regions adjacent to exons 11 and 12 were too weak to support splicing when inserted into a globin intron. UV-crosslinking experiments showed binding of hnRNP F/H just 5’ of these regions, but the hnRNP F/H target sequences did not mediate inhibition. Instead, the failure of splicing is linked to weak binding of U2AF65, and spliceosome assembly stalls prior to formation of any of the ATP-dependent complexes. We discuss mechanisms by which U2AF65 binding is facilitated in vivo.     

Characterization of an alternative splicing by a NAGNAG splice acceptor site in the porcine KIT gene

The KIT gene has been shown to have multiple functions in hematopoiesis, melanogenesis, and gametogenesis. In addition, mutations of this gene cause pigmentation disorders in humans and mice and are responsible for coat color differences in pigs. While characterizing polymorphisms in the porcine KIT gene, we detected alternative splicing (AS) of the NAGNAG splice acceptor site at the boundary of intron 4 and exon 5. This AS event generated the E and I isoforms, characterized by insertion or deletion, respectively, of CAG at the borders of coding sequence. AS patterns measured in tissue samples from two randomly selected animals did not identified any tissue-specific outcomes. Analysis of AS patterns using three breeds demonstrated that Landrace and Large White pigs expressed both the E and I isoforms. In contrast, a subset of specimens from Korean Native Pigs (KNP) yielded a single I isoform. Alignment of the sequence from several species revealed that the region between the branch point sequence (BPS) and 3′ acceptor site is conserved. However, it is appeared that the selection of either the proximal or distal splice site varied between species. To test the breed specificity the NAGNAG splice acceptor site, we constructed two lineages of minigenes from KNP and Landrace pigs harboring breed-specific mutations. The minigene splicing assay demonstrated that both types of minigenes expressed both the E and I isoforms in two host cell lines, and no differences were detected in the AS pattern between the two breeds. We conclude that the AS at the NAGNAG splice acceptor site on intron 4/exon 5 in the porcine KIT gene is the result of noise selection at the splice site by the splicing machinery. Therefore, this AS event in the porcine KIT gene is unlikely to have any relationship with the coat color variations of Landrace and KNP breeds.     

Physiological state co-regulates thousands of mammalian mRNA splicing events at tandem splice sites and alternative exons

Thousands of tandem alternative splice sites (TASS) give rise to mRNA insertion/deletion variants with small size differences. Recent work has concentrated on the question of biological relevance in general, and the physiological regulation of TASS in particular. We have quantitatively studied 11 representative TASS cases in comparison to one mutually exclusive exon case and two cassette exons (CEs) using a panel of human and mouse tissues, as well as cultured cell lines. Tissues show small but significant differences in TASS isoform ratios, with a variance 4- to 20-fold lower than seen for CEs. Remarkably, in cultured cells, all studied alternative splicing (AS) cases showed a cell-density-dependent shift of isoform ratios with similar time series profiles. A respective genome-wide co-regulation of TASS splicing was shown by next-generation mRNA sequencing data. Moreover, data from human and mouse organs indicate that this co-regulation of TASS occurs in vivo, with brain showing the strongest difference to other organs. Together, the results indicate a physiological AS regulation mechanism that functions almost independently from the splice site context and sequence.     

Switch in 3' splice site recognition between exon definition and splicing catalysis is important for sex-lethal autoregulation

Maintenance of female sexual identity in Drosophila melanogaster involves an autoregulatory loop in which the protein Sex-lethal (SXL) promotes skipping of exon 3 from its own pre-mRNA. We have used transient transfection of Drosophila Schneider cells to analyze the role of exon 3 splice sites in regulation. Our results indicate that exon 3 repression requires competition between the 5' splice sites of exons 2 and 3 but is independent of their relative strength. Two 3' splice site AG's precede exon 3. We report here that, while the distal site plays a critical role in defining the exon, the proximal site is preferentially used for the actual splicing reaction, arguing for a switch in 3' splice site recognition between exon definition and splicing catalysis. Remarkably, the presence of the two 3' splice sites is important for the efficient regulation by SXL, suggesting that SXL interferes with molecular events occurring between initial splice site communication across the exon and the splice site pairing that leads to intron removal.