Alternative splicing regulation at tandem 3' splice sites

Alternative splicing (AS) constitutes a major mechanism creating protein diversity in humans. Previous bioinformatics studies based on expressed sequence tag and mRNA data have identified many AS events that are conserved between humans and mice. Of these events, ~25% are related to alternative choices of 3′ and 5′ splice sites. Surprisingly, half of all these events involve 3′ splice sites that are exactly 3 nt apart. These tandem 3′ splice sites result from the presence of the NAGNAG motif at the acceptor splice site, recently reported to be widely spread in the human genome. Although the NAGNAG motif is common in human genes, only a small subset of sites with this motif is confirmed to be involved in AS. We examined the NAGNAG motifs and observed specific features such as high sequence conservation of the motif, high conservation of ~30 bp at the intronic regions flanking the 3′ splice site and overabundance of cis-regulatory elements, which are characteristic of alternatively spliced tandem acceptor sites and can distinguish them from the constitutive sites in which the proximal NAG splice site is selected. Our findings imply that AS at tandem splice sites and constitutive splicing of the distal NAG are highly regulated.     

Single-nucleotide polymorphisms in NAGNAG acceptors are highly predictive for variations of alternative splicing

Aberrant or modified splicing patterns of genes are causative for many human diseases. Therefore, the identification of genetic variations that cause changes in the splicing pattern of a gene is important. Elsewhere, we described the widespread occurrence of alternative splicing at NAGNAG acceptors. Here, we report a genomewide screen for single-nucleotide polymorphisms (SNPs) that affect such tandem acceptors. From 121 SNPs identified, we extracted 64 SNPs that most likely affect alternative NAGNAG splicing. We demonstrate that the NAGNAG motif is necessary and sufficient for this type of alternative splicing. The evolutionarily young NAGNAG alleles, as determined by the comparison with the chimpanzee genome, exhibit the same biases toward intron phase 1 and single-amino acid insertion/deletions that were already observed for all human NAGNAG acceptors. Since 28% of the NAGNAG SNPs occur in known disease genes, they represent preferable candidates for a more-detailed functional analysis, especially since the splice relevance for some of the coding SNPs is overlooked. Against the background of a general lack of methods for identifying splice-relevant SNPs, the presented approach is highly effective in the prediction of polymorphisms that are causal for variations in alternative splicing.     

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.     

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.     

Genome-wide study of NAGNAG alternative splicing in Arabidopsis

NAGNAG alternative splicing is one type of alternative splicing in mammals and plants. There are two opposite arguments regarding the mechanism of this NAGNAG event, i.e. whether splice variation is controllable by the cell or is just biological noise. In this paper, we systematically investigated NAGNAG acceptors in Arabidopsis thaliana using both cDNA/EST and RNA-Seq data. We identified 9,473 NAGNAG motifs, including 529 cDNA/EST-confirmed NAGNAG acceptors. A nomenclature tree for this type of alternative splicing was defined based on the cDNA/EST validation, location in the exon, sequence and expression level. Low expression of some NAGNAG motifs was observed in various tissues or pathogen-infected samples, indicating the existence of background splicing. Tissue-specific or treatment-specific differences in the dynamic profiles suggest that some NAGNAG acceptors are highly regulated.     

Accurate prediction of NAGNAG alternative splicing

Alternative splicing (AS) involving NAGNAG tandem acceptors is an evolutionarily widespread class of AS. Recent predictions of alternative acceptor usage reported better results for acceptors separated by larger distances, than for NAGNAGs. To improve the latter, we aimed at the use of Bayesian networks (BN), and extensive experimental validation of the predictions. Using carefully constructed training and test datasets, a balanced sensitivity and specificity of ≥92% was achieved. A BN trained on the combined dataset was then used to make predictions, and 81% (38/47) of the experimentally tested predictions were verified. Using a BN learned on human data on six other genomes, we show that while the performance for the vertebrate genomes matches that achieved on human data, there is a slight drop for Drosophila and worm. Lastly, using the prediction accuracy according to experimental validation, we estimate the number of yet undiscovered alternative NAGNAGs. State of the art classifiers can produce highly accurate prediction of AS at NAGNAGs, indicating that we have identified the major features of the ‘NAGNAG-splicing code’ within the splice site and its immediate neighborhood. Our results suggest that the mechanism behind NAGNAG AS is simple, stochastic, and conserved among vertebrates and beyond.     

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.     

Genetic evidence that both dNTP-stabilized and strand slippage mechanisms may dictate DNA polymerase errors within mononucleotide microsatellites

Mononucleotide microsatellites are tandem repeats of a single base pair, abundant within coding exons and frequent sites of mutation in the human genome. Because the repeated unit is one base pair, multiple mechanisms of insertion/deletion (indel) mutagenesis are possible, including strand-slippage, dNTP-stabilized, and misincorportion-misalignment. Here, we examine the effects of polymerase identity (mammalian Pols α, β, κ, and η), template sequence, dNTP pool size, and reaction temperature on indel errors during in vitro synthesis of mononucleotide microsatellites. We utilized the ratio of insertion to deletion errors as a genetic indicator of mechanism. Strikingly, we observed a statistically significant bias toward deletion errors within mononucleotide repeats for the majority of the 28 DNA template and polymerase combinations examined, with notable exceptions based on sequence and polymerase identity. Using mutator forms of Pol β did not substantially alter the error specificity, suggesting that mispairing-misalignment mechanism is not a primary mechanism. Based on our results for mammalian DNA polymerases representing three structurally distinct families, we suggest that dNTP-stabilized mutagenesis may be an alternative mechanism for mononucleotide microsatellite indel mutation. The change from a predominantly dNTP-stabilized mechanism to a strand-slippage mechanism with increasing microsatellite length may account for the differential rates of tandem repeat mutation that are observed genome-wide.     

AUG sequences are required to sustain nonsense-codon-mediated suppression of splicing

More than 90% of human genes are rich in intronic latent 5′ splice sites whose utilization in pre-mRNA splicing would introduce in-frame stop codons into the resultant mRNAs. We have therefore hypothesized that suppression of splicing (SOS) at latent 5′ splice sites regulates alternative 5′ splice site selection in a way that prevents the production of toxic nonsense mRNAs and verified this idea by showing that the removal of such in-frame stop codons is sufficient to activate latent splicing. Splicing control by SOS requires recognition of the mRNA reading frame, presumably recognizing the start codon sequence. Here we show that AUG sequences are indeed essential for SOS. Although protein translation does not seem to be required for SOS, the first AUG is shown here to be necessary but not sufficient. We further show that latent splicing can be elicited upon treatment with pactamycin—a drug known to block translation by its ability to recognize an RNA fold—but not by treatment with other drugs that inhibit translation through other mechanisms. The effect of pactamycin on SOS is dependent neither on steady-state translation nor on the pioneer round of translation. This effect is found for both transfected and endogenous genes, indicating that SOS is a natural mechanism.     

Genetic analysis of L123 of the tRNA-mimicking eukaryote release factor eRF1, an amino acid residue critical for discrimination of stop codons

In eukaryotes, the tRNA-mimicking polypeptide-chain release factor, eRF1, decodes stop codons on the ribosome in a complex with eRF3; this complex exhibits striking structural similarity to the tRNA–eEF1A–GTP complex. Although amino acid residues or motifs of eRF1 that are critical for stop codon discrimination have been identified, the details of the molecular mechanisms involved in the function of the ribosomal decoding site remain obscure. Here, we report analyses of the position-123 amino acid of eRF1 (L123 in Saccharomyces cerevisiae eRF1), a residue that is phylogenetically conserved among species with canonical and variant genetic codes. In vivo readthrough efficiency analysis and genetic growth complementation analysis of the residue-123 systematic mutants suggested that this amino acid functions in stop codon discrimination in a manner coupled with eRF3 binding, and distinctive from previously reported adjacent residues. Furthermore, aminoglycoside antibiotic sensitivity analysis and ribosomal docking modeling of eRF1 in a quasi-A/T state suggested a functional interaction between the side chain of L123 and ribosomal residues critical for codon recognition in the decoding site, as a molecular explanation for coupling with eRF3. Our results provide insights into the molecular mechanisms underlying stop codon discrimination by a tRNA-mimicking protein on the ribosome.     

Convergence and constraint in eukaryotic release factor 1 (eRF1) domain 1: the evolution of stop codon specificity

Class 1 release factor in eukaryotes (eRF1) recognizes stop codons and promotes peptide release from the ribosome. The ‘molecular mimicry’ hypothesis suggests that domain 1 of eRF1 is analogous to the tRNA anticodon stem–loop. Recent studies strongly support this hypothesis and several models for specific interactions between stop codons and residues in domain 1 have been proposed. In this study we have sequenced and identified novel eRF1 sequences across a wide diversity of eukaryotes and re-evaluated the codon-binding site by bioinformatic analyses of a large eRF1 dataset. Analyses of the eRF1 structure combined with estimates of evolutionary rates at amino acid sites allow us to define the residues that are under structural (i.e. those involved in intramolecular interactions) versus non-structural selective constraints. Furthermore, we have re-assessed convergent substitutions in the ciliate variant code eRF1s using maximum likelihood-based phylogenetic approaches. Our results favor the model proposed by Bertram et al. that stop codons bind to three ‘cavities’ on the protein surface, although we suggest that the stop codon may bind in the opposite orientation to the original model. We assess the feasibility of this alternative binding orientation with a triplet stop codon and the eRF1 domain 1 structures using molecular modeling techniques.     

Regulation of BCL-X splicing reveals a role for the polypyrimidine tract binding protein (PTBP1/hnRNP I) in alternative 5′ splice site selection

Alternative splicing (AS) modulates many physiological and pathological processes. For instance, AS of the BCL-X gene balances cell survival and apoptosis in development and cancer. Herein, we identified the polypyrimidine tract binding protein (PTBP1) as a direct regulator of BCL-X AS. Overexpression of PTBP1 promotes selection of the distal 5′ splice site in BCL-X exon 2, generating the pro-apoptotic BCL-Xs splice variant. Conversely, depletion of PTBP1 enhanced splicing of the anti-apoptotic BCL-XL variant. In vivo cross-linking experiments and site-directed mutagenesis restricted the PTBP1 binding site to a polypyrimidine tract located between the two alternative 5′ splice sites. Binding of PTBP1 to this site was required for its effect on splicing. Notably, a similar function of PTBP1 in the selection of alternative 5′ splice sites was confirmed using the USP5 gene as additional model. Mechanistically, PTBP1 displaces SRSF1 binding from the proximal 5′ splice site, thus repressing its selection. Our study provides a novel mechanism of alternative 5′ splice site selection by PTBP1 and indicates that the presence of a PTBP1 binding site between two alternative 5′ splice sites promotes selection of the distal one, while repressing the proximal site by competing for binding of a positive regulator.     

Ribonucleotide reductase genes of Bacillus prophages: a refuge to introns and intein coding sequences

The ribonucleotide reductase gene tandem bnrdE/bnrdF in SPβ-related prophages of different Bacillus spp. isolates presents different configurations of intervening sequences, comprising one to three of six non-homologous splicing elements. Insertion sites of group I introns and intein DNA are clustered in three relatively short segments encoding functionally important domains of the ribonucleotide reductase. Comparison of the bnrdE homologs reveals mutual exclusion of a group I intron and an intein coding sequence flanking the codon that specifies a conserved cysteine. In vivo splicing was demonstrated for all introns. However, for two of them a part of the mRNA precursor molecules remains unspliced. Intergenic bnrdE–bnrdF regions are unexpectedly long, comprising between 238 and 541 nt. The longest encodes a putative polypeptide related to HNH homing endonucleases.