Assessing the fraction of short-distance tandem splice sites under purifying selection

Nucleotides in both donor and acceptor splice sites are responsible for choice in NAGNAG tandem splice sites

Among alternative splicing events in the human transcriptome, tandem NAGNAG acceptor splice sites represent an appreciable proportion. Both proximal and distal NAG can be used to produce two splicing isoforms differing by three nucleotides. In some cases, the upstream exon can be alternatively spliced as well, which further increases the number of possible transcripts. In this study, we showed that NAG choice in tandem splice site depends considerably not only on the concerned acceptor, but also on the upstream donor splice site sequence. Using an extensive set of experiments with systematically modified two-exonic minigene systems of AFAP1L2 or CSTD gene, we recognized the third and fifth intronic upstream donor splice site position and the tandem acceptor splice site region spanning from -10 to +2, including NAGNAG itself, as the main drivers. In addition, competition between different branch points and their composition were also shown to play a significant role in NAG choice. All these nucleotide effects appeared almost additive, which explained the high variability in proximal versus distal NAG usage.

Selection for reduced translation costs at the intronic 5' end in fungi

It is generally believed that introns are not translated; therefore, the potential intronic features that may be related to the translation step (occurring after splicing) have yet to be thoroughly studied. Here, focusing on four fungi, we performed for the first time a comprehensive study aimed at characterizing how translation efficiency is encoded in introns and affects their evolution. By analysing their intronome we provide evidence of selection for STOP codons close to the intronic 5′ end, and show that the beginning of introns are selected for significantly high translation, presumably to reduce translation and metabolic costs in cases of non-spliced introns. Ribosomal profiling data analysis in Saccharomyces cerevisiae supports the conjecture that in this organism intron retention frequently occurs, introns are partially translated, and their translation efficiency affects organismal fitness. We show that the reported results are more significant in highly translated and highly spliced genes, but are not associated only with genes with a specific function. We also discuss the potential relation of the reported signals to efficient nonsense-mediated decay due to splicing errors. These new discoveries are supported by population-genetics considerations. In addition, they are contributory steps towards a broader understanding of intron evolution and the effect of silent mutations on gene expression and organismal fitness.

It's a bit over, is that ok? The subtle surplus from tandem alternative splicing

Tandem alternative splice sites (TASS) form a defined class of alternative splicing and give rise to mRNA insertion/deletion variants with only small size differences. Previous work has confirmed evolutionary conservation of TASS elements while many cases show only low tissue specificity of isoform ratios. We pinpoint stochasticity and noise as important methodological issues for the dissection of TASS isoform patterns. Resolving such uncertainties, a recent report showed regulation in a cell culture system, with shifts of alternative splicing isoform ratios dependent on cell density. This novel type of regulation affects not only multiple TASS isoforms, but also other alternative splicing classes, in a concerted manner. Here, we discuss how specific regulatory network architectures may be realized through the novel regulation type and highlight the role of differential isoform functions as a key step in order to better understand the functional role of TASS.

Alternative splicing at GYNNGY 5' splice sites: more noise, less regulation

Numerous eukaryotic genes are alternatively spliced. Recently, deep transcriptome sequencing has skyrocketed proportion of alternatively spliced genes; over 95% human multi-exon genes are alternatively spliced. One fundamental question is: are all these alternative splicing (AS) events functional? To look into this issue, we studied the most common form of alternative 5' splice sites-GYNNGYs (Y = C/T), where both GYs can function as splice sites. Global analyses suggest that splicing noise (due to stochasticity of splicing process) can cause AS at GYNNGYs, evidenced by higher AS frequency in non-coding than in coding regions, in non-conserved than in conserved genes and in lowly expressed than in highly expressed genes. However, ∼20% AS GYNNGYs in humans and ∼3% in mice exhibit tissue-dependent regulation. Consistent with being functional, regulated GYNNGYs are more conserved than unregulated ones. And regulated GYNNGYs have distinctive sequence features which may confer regulation. Particularly, each regulated GYNNGY comprises two splice sites more resembling each other than unregulated GYNNGYs, and has more conserved downstream flanking intron. Intriguingly, most regulated GYNNGYs may tune gene expression through coupling with nonsense-mediated mRNA decay, rather than encode different proteins. In summary, AS at GYNNGY 5' splice sites is primarily splicing noise, and secondarily a way of regulation.

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.

A subtle alternative splicing event gives rise to a widely expressed human RNase k isoform

Subtle alternative splicing leads to the formation of RNA variants lacking or including a small number of nucleotides. To date, the impact of subtle alternative splicing phenomena on protein biosynthesis has been studied in frame-preserving incidents. On the contrary, mRNA isoforms derived from frame-shifting events were poorly studied and generally characterized as non-coding. This work provides evidence for a frame-shifting subtle alternative splicing event which results in the production of a novel protein isoform. We applied a combined molecular approach for the cloning and expression analysis of a human RNase κ transcript (RNase κ-02) which lacks four consecutive bases compared to the previously isolated RNase κ isoform. RNase κ-02 mRNA is expressed in all human cell lines tested end encodes the synthesis of a 134-amino-acid protein by utilizing an alternative initiation codon. The expression of RNase κ-02 in the cytoplasm of human cells was verified by Western blot and immunofluorescence analysis using a specific polyclonal antibody developed on the basis of the amino-acid sequence difference between the two protein isoforms. The results presented here show that subtle changes during mRNA splicing can lead to the expression of significantly altered protein isoforms.

Cloning and transcriptional expression of mouse mannosyltransferase IV/V cDNA, which is involved in the synthesis of lipid-linked oligosaccharides

We cloned the mouse mannosyltransferase IV/V gene (mALG11) from FM3A cells by a bioinformatic approach. The ORF contained 1476 bp encoding 492 amino acids. The cloned mALG11 complemented the growth defect of the Saccharomyces cerevisiae ALG11Δ mutant. In addition, we detected a variant cDNA by alternate splicing that had an additional four-nucleotide ATGC insertion at base 276 of the ORF. Consequently the variant cDNA encoded a truncated protein with 92 amino acids, lacking the glycosyltransferase group-1 domain. The variant cDNA occurs in many mouse strains according to EST database searches. Moreover, we detected it in FM3A cDNA, but we did not detect any such variants in the human EST database or in HeLa cDNA, although human ALG11 (hALG11) genomic DNA has the same sequence around the intron-exon boundaries as those of mALG11 genomic DNA. Hence, we concluded that there is different transcriptional control mechanism between mALG11 and hALG11.

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ∼15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.

Sequence features involved in the mechanism of 3' splice junction wobbling

BACKGROUND

Alternative splicing is an important mechanism mediating the diversified functions of genes in multicellular organisms, and such event occurs in around 40-60% of human genes. Recently, a new splice-junction wobbling mechanism was proposed that subtle modifications exist in mRNA maturation by alternatively choosing at 5'- GTNGT and 3'- NAGNAG, which created single amino acid insertion and deletion isoforms.

RESULTS

By browsing the Alternative Splicing Database information, we observed that most 3' alternative splice site choices occur within six nucleotides of the dominant splice site and the incidence significantly decreases further away from the dominant acceptor site. Although a lower frequency of alternative splicing occurs within the intronic region (alternative splicing at the proximal AG) than in the exonic region (alternative splicing at the distal AG), alternative AG sites located within the intronic region show stronger potential as the acceptor. These observations revealed that the choice of 3' splice sites during 3' splicing junction wobbling could depend on the distance between the duplicated AG and the branch point site (BPS). Further mutagenesis experiments demonstrated that the distance of AG-to-AG and BPS-to-AG can greatly influence 3' splice site selection. Knocking down a known alternative splicing regulator, hSlu7, failed to affect wobble splicing choices.

CONCLUSION

Our results implied that nucleotide distance between proximal and distal AG sites has an important regulatory function. In this study, we showed that occurrence of 3' wobble splicing occurs in a distance-dependent manner and that most of this wobble splicing is probably caused by steric hindrance from a factor bound at the neighboring tandem motif sequence.

TassDB2 - A comprehensive database of subtle alternative splicing events

BACKGROUND

Subtle alternative splicing events involving tandem splice sites separated by a short (2-12 nucleotides) distance are frequent and evolutionarily widespread in eukaryotes, and a major contributor to the complexity of transcriptomes and proteomes. However, these events have been either omitted altogether in databases on alternative splicing, or only the cases of experimentally confirmed alternative splicing have been reported. Thus, a database which covers all confirmed cases of subtle alternative splicing as well as the numerous putative tandem splice sites (which might be confirmed once more transcript data becomes available), and allows to search for tandem splice sites with specific features and download the results, is a valuable resource for targeted experimental studies and large-scale bioinformatics analyses of tandem splice sites. Towards this goal we recently set up TassDB (Tandem Splice Site DataBase, version 1), which stores data about alternative splicing events at tandem splice sites separated by 3 nt in eight species.

DESCRIPTION

We have substantially revised and extended TassDB. The currently available version 2 contains extensive information about tandem splice sites separated by 2-12 nt for the human and mouse transcriptomes including data on the conservation of the tandem motifs in five vertebrates. TassDB2 offers a user-friendly interface to search for specific genes or for genes containing tandem splice sites with specific features as well as the possibility to download result datasets. For example, users can search for cases of alternative splicing where the proportion of EST/mRNA evidence supporting the minor isoform exceeds a specific threshold, or where the difference in splice site scores is specified by the user. The predicted impact of each event on the protein is also reported, along with information about being a putative target for the nonsense-mediated decay (NMD) pathway. Links are provided to the UCSC genome browser and other external resources.

CONCLUSION

TassDB2, available via http://www.tassdb.info, provides comprehensive resources for researchers interested in both targeted experimental studies and large-scale bioinformatics analyses of short distance tandem splice sites.

Identification and characterization of NAGNAG alternative splicing in the moss Physcomitrella patens

BACKGROUND

Alternative splicing (AS) involving tandem acceptors that are separated by three nucleotides (NAGNAG) is an evolutionarily widespread class of AS, which is well studied in Homo sapiens (human) and Mus musculus (mouse). It has also been shown to be common in the model seed plants Arabidopsis thaliana and Oryza sativa (rice). In one of the first studies involving sequence-based prediction of AS in plants, we performed a genome-wide identification and characterization of NAGNAG AS in the model plant Physcomitrella patens, a moss.

RESULTS

Using Sanger data, we found 295 alternatively used NAGNAG acceptors in P. patens. Using 31 features and training and test datasets of constitutive and alternative NAGNAGs, we trained a classifier to predict the splicing outcome at NAGNAG tandem splice sites (alternative splicing, constitutive at the first acceptor, or constitutive at the second acceptor). Our classifier achieved a balanced specificity and sensitivity of >or= 89%. Subsequently, a classifier trained exclusively on data well supported by transcript evidence was used to make genome-wide predictions of NAGNAG splicing outcomes. By generation of more transcript evidence from a next-generation sequencing platform (Roche 454), we found additional evidence for NAGNAG AS, with altogether 664 alternative NAGNAGs being detected in P. patens using all currently available transcript evidence. The 454 data also enabled us to validate the predictions of the classifier, with 64% (80/125) of the well-supported cases of AS being predicted correctly.

CONCLUSION

NAGNAG AS is just as common in the moss P. patens as it is in the seed plants A. thaliana and O. sativa (but not conserved on the level of orthologous introns), and can be predicted with high accuracy. The most informative features are the nucleotides in the NAGNAG and in its immediate vicinity, along with the splice sites scores, as found earlier for NAGNAG AS in animals. Our results suggest that the mechanism behind NAGNAG AS in plants is similar to that in animals and is largely dependent on the splice site and its immediate neighborhood.

Decrypting the genome's alternative messages

Alternative splicing of messenger RNA (mRNA) precursors affects the majority of human genes, has a considerable impact on eukaryotic gene function and offers distinct opportunities for regulation. Alterations in alternative splicing can cause or modify the progression of a significant number of pathologies. Recent high-throughput technologies have uncovered a wealth of transcript diversity generated by alternative splicing, as well as examples for how this diversity can be established and become misregulated. A variety of mechanisms modulate splice site choice coordinately with other cellular processes, from transcription and mRNA editing or decay to miRNA-based regulation and telomerase function. Alternative splicing studies can contribute to our understanding of multiple biological processes, including genetic diversity, speciation, cell/stem cell differentiation, nervous system function, neuromuscular disorders and tumour progression.

Molecular determinants and evolutionary dynamics of wobble splicing

Alternative splicing at tandem splice sites (wobble splicing) is widespread in many species, but the mechanisms specifying the tandem sites remain poorly understood. Here, we used synaptotagmin I as a model to analyze the phylogeny of wobble splicing spanning more than 300 My of insect evolution. Phylogenetic analysis indicated that the occurrence of species-specific wobble splicing was related to synonymous variation at tandem splice sites. Further mutagenesis experiments demonstrated that wobble splicing could be lost by artificially induced synonymous point mutations due to destruction of splice acceptor sites. In contrast, wobble splicing could not be correctly restored through mimicking an ancestral tandem acceptor by artificial synonymous mutation in in vivo splicing assays, which suggests that artificial tandem splice sites might be incompatible with normal wobble splicing. Moreover, combining comparative genomics with hybrid minigene analysis revealed that alternative splicing has evolved from the 3′ tandem donor to the 5′ tandem acceptor in Culex pipiens, as a result of an evolutionary shift of cis element sequences from 3′ to 5′ splice sites. These data collectively suggest that the selection of tandem splice sites might not simply be an accident of history but rather in large part the result of coevolution between splice site and cis element sequences as a basis for wobble splicing. An evolutionary model of wobble splicing is proposed.

Alternative splicing at NAGNAG acceptors in Arabidopsis thaliana SR and SR-related protein-coding genes

Background

Several recent studies indicate that alternative splicing in Arabidopsis and other plants is a common mechanism for post-transcriptional modulation of gene expression. However, few analyses have been done so far to elucidate the functional relevance of alternative splicing in higher plants. Representing a frequent and universal subtle alternative splicing event among eukaryotes, alternative splicing at NAGNAG acceptors contributes to transcriptome diversity and therefore, proteome plasticity. Alternatively spliced NAGNAG acceptors are overrepresented in genes coding for proteins with RNA-recognition motifs (RRMs). As SR proteins, a family of RRM-containing important splicing factors, are known to be extensively alternatively spliced in Arabidopsis, we analyzed alternative splicing at NAGNAG acceptors in SR and SR-related genes.

Results

In a comprehensive analysis of the Arabidopsis thaliana genome, we identified 6,772 introns that exhibit a NAGNAG acceptor motif. Alternative splicing at these acceptors was assessed using available EST data, complemented by a sequence-based prediction method. Of the 36 identified introns within 30 SR and SR-related protein-coding genes that have a NAGNAG acceptor, we selected 15 candidates for an experimental analysis of alternative splicing under several conditions. We provide experimental evidence for 8 of these candidates being alternatively spliced. Quantifying the ratio of NAGNAG-derived splice variants under several conditions, we found organ-specific splicing ratios in adult plants and changes in seedlings of different ages. Splicing ratio changes were observed in response to heat shock and most strikingly, cold shock. Interestingly, the patterns of differential splicing ratios are similar for all analyzed genes.

Conclusion

NAGNAG acceptors frequently occur in the Arabidopsis genome and are particularly prevalent in SR and SR-related protein-coding genes. A lack of extensive EST coverage can be compensated by using the proposed sequence-based method to predict alternative splicing at these acceptors. Our findings indicate that the differential effects on NAGNAG alternative splicing in SR and SR-related genes are organ- and condition-specific rather than gene-specific.

Widespread and subtle: alternative splicing at short-distance tandem sites

Alternative splicing at donor or acceptor sites located just a few nucleotides apart is widespread in many species. It results in subtle changes in the transcripts and often in the encoded proteins. Several of these tandem splice events contribute to the repertoire of functionally different proteins, whereas many are neutral or deleterious. Remarkably, some of the functional events are differentially spliced in tissues or developmental stages, whereas others exhibit constant splicing ratios, indicating that function is not always associated with differential splicing. Stochastic splice site selection seems to play a major role in these processes. Here, we review recent progress in understanding functional and evolutionary aspects as well as the mechanism of splicing at short-distance tandem sites.

Selection against tandem splice sites affecting structured protein regions

BACKGROUND

Alternative selection of splice sites in tandem donors and acceptors is a major mode of alternative splicing. Here, we analyzed whether in-frame tandem sites leading to subtle mRNA insertions/deletions of 3, 6, or 9 nucleotides are under natural selection.

RESULTS

We found multiple lines of evidence that the human protein coding sequences are under selection against such in-frame tandem splice events, indicating that these events are often deleterious. The strength of selection is not homogeneous within the coding sequence as protein regions that fold into a fixed 3D structure (intrinsically ordered) are under stronger selection, especially against sites with a strong minor splice site. Investigating structures of functional protein domains, we found that tandem acceptors are preferentially located at the domain surface and outside structural elements such as helices and sheets. Using three-species comparisons, we estimate that more than half of all mutations that create NAGNAG acceptors in the coding region have been eliminated by selection.

CONCLUSION

We estimate that ~2,400 introns are under selection against possessing a tandem site.