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.     

Discovery of cSNPs in pig using full-length enriched cDNA libraries of the Korean native pig as a source of genetic diversity

Clones from full-length enriched cDNA libraries serve as valuable resources for functional genomic studies. We analyzed 3.210 chromatograms obtained from sequencing the 5′-ends of brainstem, liver, neocortex, and spleen clones derived from full-length enriched cDNA libraries of Korean native pigs. In addition, 50,000 pig EST sequence trace files were obtained from Genbank and combined with our sequencing information for SNP identificationin silico. For the SNP analysis, neocortex, and liver libraries were newly constructed, whereas the sequencing results from brainstem and spleen libraries were from previously constructed libraries. The putative SNPs from thein silico analysis were confirmed by genomic PCR from a group of 20 pigs of four different breeds. Using this approach, 86% of cSNPs identifiedin silico were confirmed and the SNP detection frequency was 1 SNP per 338 bp. Interestingly, we found a valine deletion at amino acid position 126 of the neuronal and endocrine protein gene in the Korean native pig. We confirmed that this deletion was caused by alternative splicing at the NAGNAG acceptors. Our study shows that large-scale EST sequencing of Korean native pigs can be effectively employed for natural polymorphism-based pig genome analysis.     

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.     

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.     

Alternative Splicing at a NAGNAG Acceptor Site as a Novel Phenotype Modifier

Approximately 30% of alleles causing genetic disorders generate premature termination codons (PTCs), which are usually associated with severe phenotypes. However, bypassing the deleterious stop codon can lead to a mild disease outcome. Splicing at NAGNAG tandem splice sites has been reported to result in insertion or deletion (indel) of three nucleotides. We identified such a mechanism as the origin of the mild to asymptomatic phenotype observed in cystic fibrosis patients homozygous for the E831X mutation (2623G>T) in the CFTR gene. Analyses performed on nasal epithelial cell mRNA detected three distinct isoforms, a considerably more complex situation than expected for a single nucleotide substitution. Structure-function studies and in silico analyses provided the first experimental evidence of an indel of a stop codon by alternative splicing at a NAGNAG acceptor site. In addition to contributing to proteome plasticity, alternative splicing at a NAGNAG tandem site can thus remove a disease-causing UAG stop codon. This molecular study reveals a naturally occurring mechanism where the effect of either modifier genes or epigenetic factors could be suspected. This finding is of importance for genetic counseling as well as for deciding appropriate therapeutic strategies.     

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.     

A network of conserved co-occurring motifs for the regulation of alternative splicing

Cis-acting short sequence motifs play important roles in alternative splicing. It is now possible to identify such sequence motifs as conserved sequence patterns in genome sequence alignments. Here, we report the systematic search for motifs in the neighboring introns of alternatively spliced exons by using comparative analysis of mammalian genome alignments. We identified 11 conserved sequence motifs that might be involved in the regulation of alternative splicing. These motifs are not only significantly overrepresented near alternatively spliced exons, but they also co-occur with each other, thus, forming a network of cis-elements, likely to be the basis for context-dependent regulation. Based on this finding, we applied the motif co-occurrence to predict alternatively skipped exons. We verified exon skipping in 29 cases out of 118 predictions (25%) by EST and mRNA sequences in the databases. For the predictions not verified by the database sequences, we confirmed exon skipping in 10 additional cases by using both RT–PCR experiments and the publicly available RNA-Seq data. These results indicate that even more alternative splicing events will be found with the progress of large-scale and high-throughput analyses for various tissue samples and developmental stages.     

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.     

Alternative splicing and expression profile analysis of expressed sequence tags in domestic pig

Domestic pig （Sus scrofa domestica） is one of the most important mammals to humans. Alternative splicing is a cellular mechanism in eukaryotes that greatly increases the diversity of gene products. Expression sequence tags （ESTs） have been widely used for gene discovery, expression profile analysis, and alternative splicing detection. In this study, a total of 712,905 ESTs extracted from 101 different nonnormalized EST libraries of the domestic pig were analyzed. These EST libraries cover the nervous system, digestive system, immune system, and meat production related tissues from embryo, newborn, and adult pigs, making contributions to the analysis of alternative splicing variants as well as expression profiles in various stages of tissues. A modified approach was designed to cluster and assemble large EST datasets, aiming to detect alternative splicing together with EST abundance of each splicing variant. Much efforts were made to classify alternative splicing into different types and apply different filters to each type to get more reliable results. Finally, a total of 1,223 genes with average 2.8 splicing variants were detected among 16,540 unique genes. The overview of expression profiles would change when we take alternative splicing into account.     

Bioinformatics analysis of alternative splicing

Over the past few years, the analysis of alternative splicing using bioinformatics has emerged as an important new field, and has significantly changed our view of genome function. One exciting front has been the analysis of microarray data to measure alternative splicing genome-wide. Pioneering studies of both human and mouse data have produced algorithms for discerning evidence of alternative splicing and clustering genes and samples by their alternative splicing patterns. Moreover, these data indicate the presence of alternative splice forms in up to 80 per cent of human genes. Comparative genomics studies in both mammals and insects have demonstrated that alternative splicing can in some cases be predicted directly from comparisons of genome sequences, based on heightened sequence conservation and exon length. Such studies have also provided new insights into the connection between alternative splicing and a variety of evolutionary processes such as Alu-based exonisation, exon creation and loss. A number of groups have used a combination of bioinformatics, comparative genomics and experimental validation to identify new motifs for splice regulatory factors, analyse the balance of factors that regulate alternative splicing, and propose a new mechanism for regulation based on the interaction of alternative splicing and nonsense-mediated decay. Bioinformatics studies of the functional impact of alternative splicing have revealed a wide range of regulatory mechanisms, from NAGNAG sites that add a single amino acid; to short peptide segments that can play surprisingly complex roles in switching protein conformation and function (as in the Piccolo C2A domain); to events that entirely remove a specific protein interaction domain or membrane anchoring domain. Common to many bioinformatics studies is a new emphasis on graph representations of alternative splicing structures, which have many advantages for analysis.     

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.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

Identification and characterization of NAGNAG alternative splicing in the moss <Emphasis Type="Italic">Physcomitrella patens</Emphasis>

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.     

Alternate Promoters in the Rat Aromatic l-Amino Acid Decarboxylase Gene for Neuronal and Nonneuronal Expression: An In Situ Hybridization Study

Abstract: Aromatic l -amino acid decarboxylase (AADC) is found in both neuronal cells and nonneuronal cells, and a single gene encodes rat AADC in both neuronal and nonneuronal tissues. However, two cDNAs for this enzyme have been identified: one from the liver and the other from pheochromocytoma. Exons 1a and 1b are found in the liver cDNA and the pheochromocytoma cDNA, respectively. In the third exon (exon 2), there are two alternatively utilized splicing acceptors specific to these exons, 1a and 1b. Structural analysis of the rat AADC gene showed that both alternative promoter usage and alternative splicing are operative for the differential expression of this gene. To demonstrate whether alternative promoter usage and splicing are tissue specific and whether the exons 1a and 1b are differentially and specifically transcribed in nonneuronal and neuronal cells, respectively, in situ hybridization histochemistry for the rat brain, adrenal gland, liver, and kidney was carried out using these two exon probes. The exon 1a probe specifically identified AADC mRNA only in nonneuronal cells, including the liver and kidney, and the exon 1b probe localized AADC mRNA to monoaminergic neurons in the CNS and the adrenal medulla. Thus, both alternative promoter usage and differential splicing are in fact operative for the tissue-specific expression of the rat AADC gene.     

Pressure-Overload Cardiac Hypertrophy Is Associated with Distinct Alternative Splicing Due to Altered Expression of Splicing Factors

Chronic pressure-overload cardiac hypertrophy is associated with an increased risk of morbidity/mortality, largely due to maladaptive remodeling and dilatation that progresses to dilated cardiomyopathy. Alternative splicing is an important biological mechanism that generates proteomic complexity and diversity. The recent development of next-generation RNA sequencing has improved our understanding of the qualitative signatures associated with alternative splicing in various biological conditions. However, the role of alternative splicing in cardiac hypertrophy is yet unknown. The present study employed RNA-Seq and a bioinformatic approach to detect the RNA splicing regulatory elements involved in alternative splicing during pressure-overload cardiac hypertrophy. We found GC-rich exonic motifs that regulate intron retention in 5′ UTRs and AT-rich exonic motifs that are involved in exclusion of the AT-rich elements that cause mRNA instability in 3′ UTRs. We also identified motifs in the intronic regions involved in exon exclusion and inclusion, which predicted splicing factors that bind to these motifs. We found, through Western blotting, that the expression levels of three splicing factors, ESRP1, PTB and SF2/ASF, were significantly altered during cardiac hypertrophy. Collectively, the present results suggest that chronic pressure-overload hypertrophy is closely associated with distinct alternative splicing due to altered expression of splicing factors.