Predicting functional alternative splicing by measuring RNA selection pressure from multigenome alignments

Src dependency of the regulation of LTP by alternative splicing of GRIN1 exon 5

Alternative splicing of Grin1 exon 5 regulates induction of long-term potentiation (LTP) at Schaffer collateral-CA1 synapses: LTP in mice lacking the GluN1 exon 5-encoded N1 cassette (GluN1a mice) is significantly increased compared with that in mice compulsorily expressing this exon (GluN1b mice). The mechanism underlying this difference is unknown. Here, we report that blocking the non-receptor tyrosine kinase Src prevents induction of LTP in GluN1a mice but not in GluN1b. We find that activating Src enhances pharmacologically isolated synaptic N-methyl-d-aspartate receptor (NMDAR) currents in GluN1a mice but not in GluN1b. Moreover, we observe that Src activation increases the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor component of Schaffer collateral-evoked excitatory post-synaptic potentials in GluN1a mice, but this increase is prevented by blocking NMDARs. We conclude that at these synapses, NMDARs in GluN1a mice are subject to upregulation by Src that mediates induction of LTP, whereas NMDARs in GluN1b mice are not regulated by Src, leading to Src-resistance of LTP. Thus, we have uncovered that a key regulatory mechanism for synaptic potentiation is gated by differential splicing of exon 5 of Grin1. This article is part of a discussion meeting issue 'Long-term potentiation: 50 years on'.

The Impact of Pro-Inflammatory Cytokines on Alternative Splicing Patterns in Human Islets

Alternative splicing (AS) within the β cell has been proposed as one potential pathway that may exacerbate autoimmunity and unveil novel immunogenic epitopes in type 1 diabetes (T1D). We employed a computational strategy to prioritize pathogenic splicing events in human islets treated with IL-1β + IFN-γ as an ex vivo model of T1D and coupled this analysis with a k-mer based approach to predict RNA binding proteins involved in AS. In total, 969 AS events were identified in cytokine-treated islets, with the majority (44.8%) involving a skipped exon. ExonImpact identified 129 events predicted to impact protein structure. AS occurred with high frequency in MHC Class II-related mRNAs, and targeted qPCR validated reduced inclusion of Exon5 in the MHC Class II gene HLA-DMB. Single molecule RNA FISH confirmed increased HLA-DMB splicing in pancreatic sections from human donors with established T1D and autoantibody positivity. Serine and Arginine Rich Splicing Factor 2 was implicated in 37.2% of potentially pathogenic events, including Exon5 exclusion in HLA-DMB. Together, these data suggest that dynamic control of AS plays a role in the β cell response to inflammatory signals during T1D evolution.

Control of Long-Term Synaptic Potentiation and Learning by Alternative Splicing of the NMDA Receptor Subunit GluN1

NMDA receptors (NMDARs) are critical for physiological synaptic plasticity, learning, and memory and for pathological plasticity and neuronal death. The GluN1 subunit is encoded by a single gene, GRIN1, with 8 splice variants, but whether the diversity generated by this splicing has physiological consequences remains enigmatic. Here, we generate mice lacking from the GluN1 exon 5-encoded N1 cassette (GluN1a mice) or compulsorily expressing this exon (GluN1b mice). Despite no differences in basal synaptic transmission, long-term potentiation in the hippocampus is significantly enhanced in GluN1a mice compared with that in GluN1b mice. Furthermore, GluN1a mice learn more quickly and have significantly better spatial memory performance than do GluN1b mice. In addition, in human iPSC-derived neurons in autism spectrum disorder NMDARs show characteristics of N1-lacking GluN1. Our findings indicate that alternative splicing of GluN1 is a mechanism for controlling physiological long-lasting synaptic potentiation, learning, and memory.

ExonImpact: Prioritizing Pathogenic Alternative Splicing Events

Alternative splicing (AS) is a closely regulated process that allows a single gene to encode multiple protein isoforms, thereby contributing to the diversity of the proteome. Dysregulation of the splicing process has been found to be associated with many inherited diseases. However, among the pathogenic AS events, there are numerous "passenger" events whose inclusion or exclusion does not lead to significant changes with respect to protein function. In this study, we evaluate the secondary and tertiary structural features of proteins associated with disease-causing and neutral AS events, and show that several structural features are strongly associated with the pathological impact of exon inclusion. We further develop a machine-learning-based computational model, ExonImpact, for prioritizing and evaluating the functional consequences of hitherto uncharacterized AS events. We evaluated our model using several strategies including cross-validation, and data from the Gene-Tissue Expression (GTEx) and ClinVar databases. ExonImpact is freely available at http://watson.compbio.iupui.edu/ExonImpact.

A proteogenomic approach to understand splice isoform functions through sequence and expression-based computational modeling

The products of multi-exon genes are a mixture of alternatively spliced isoforms, from which the translated proteins can have similar, different or even opposing functions. It is therefore essential to differentiate and annotate functions for individual isoforms. Computational approaches provide an efficient complement to expensive and time-consuming experimental studies. The input data of these methods range from DNA sequence, to RNA selection pressure, to expressed sequence tags, to full-length complementary DNA, to exon array, to RNA-seq expression, to proteomic data. Notably, RNA-seq technology generates quantitative profiling of transcript expression at the genome scale, with an unprecedented amount of expression data available for developing isoform function prediction methods. Integrative analysis of these data at different molecular levels enables a proteogenomic approach to systematically interrogate isoform functions. Here, we briefly review the state-of-the-art methods according to their input data sources, discuss their advantages and limitations and point out potential ways to improve prediction accuracies.

PDZD7 and hearing loss: More than just a modifier

Deafness is the most frequent sensory disorder. With over 90 genes and 110 loci causally implicated in non-syndromic hearing loss, it is phenotypically and genetically heterogeneous. Here, we investigate the genetic etiology of deafness in four families of Iranian origin segregating autosomal recessive non-syndromic hearing loss (ARNSHL). We used a combination of linkage analysis, homozygosity mapping, and a targeted genomic enrichment platform to simultaneously screen 90 known deafness-causing genes for pathogenic variants. Variant segregation was confirmed by Sanger sequencing. Linkage analysis and homozygosity mapping showed segregation with the DFNB57 locus on chromosome 10 in two families. Targeted genomic enrichment with massively parallel sequencing identified causal variants in PDZD7: a homozygous missense variant (p.Gly103Arg) in one family and compound heterozygosity for missense (p.Met285Arg) and nonsense (p.Tyr500Ter) variants in the second family. Screening of two additional families identified two more variants: (p.Gly228Arg) and (p.Gln526Ter). Variant segregation with the hearing loss phenotype was confirmed in all families by Sanger sequencing. The missense variants are predicted to be deleterious, and the two nonsense mutations produce null alleles. This report is the first to show that mutations in PDZD7 cause ARNSHL, a finding that offers addition insight into the USH2 interactome. We also describe a novel likely disease-causing mutation in CIB2 and illustrate the complexity associated with gene identification in diseases that exhibit large genetic and phenotypic heterogeneity.

Comparative genomics of grass EST libraries reveals previously uncharacterized splicing events in crop plants

BACKGROUND

Crop plants such as rice, maize and sorghum play economically-important roles as main sources of food, fuel, and animal feed. However, current genome annotations of crop plants still suffer false-positive predictions; a more comprehensive registry of alternative splicing (AS) events is also in demand. Comparative genomics of crop plants is largely unexplored.

RESULTS

We performed a large-scale comparative analysis (ExonFinder) of the expressed sequence tag (EST) library from nine grass plants against three crop genomes (rice, maize, and sorghum) and identified 2,879 previously-unannotated exons (i.e., novel exons) in the three crops. We validated 81% of the tested exons by RT-PCR-sequencing, supporting the effectiveness of our in silico strategy. Evolutionary analysis reveals that the novel exons, comparing with their flanking annotated ones, are generally under weaker selection pressure at the protein level, but under stronger pressure at the RNA level, suggesting that most of the novel exons also represent novel alternatively spliced variants (ASVs). However, we also observed the consistency of evolutionary rates between certain novel exons and their flanking exons, which provided further evidence of their co-occurrence in the transcripts, suggesting that previously-annotated isoforms might be subject to erroneous predictions. Our validation showed that 54% of the tested genes expressed the newly-identified isoforms that contained the novel exons, rather than the previously-annotated isoforms that excluded them. The consistent results were steadily observed across cultivated (Oryza sativa and O. glaberrima) and wild (O. rufipogon and O. nivara) rice species, asserting the necessity of our curation of the crop genome annotations. Our comparative analyses also inferred the common ancestral transcriptome of grass plants and gain- and loss-of-ASV events.

CONCLUSIONS

We have reannotated the rice, maize, and sorghum genomes, and showed that evolutionary rates might serve as an indicator for determining whether the identified exons were alternatively spliced. This study not only presents an effective in silico strategy for the improvement of plant annotations, but also provides further insights into the role of AS events in the evolution and domestication of crop plants. ExonFinder and the novel exons/ASVs identified are publicly accessible at http://exonfinder.sourceforge.net/ .

Unraveling patterns of site-to-site synonymous rates variation and associated gene properties of protein domains and families

In protein-coding genes, synonymous mutations are often thought not to affect fitness and therefore are not subject to natural selection. Yet increasingly, cases of non-neutral evolution at certain synonymous sites were reported over the last decade. To evaluate the extent and the nature of site-specific selection on synonymous codons, we computed the site-to-site synonymous rate variation (SRV) and identified gene properties that make SRV more likely in a large database of protein-coding gene families and protein domains. To our knowledge, this is the first study that explores the determinants and patterns of the SRV in real data. We show that the SRV is widespread in the evolution of protein-coding sequences, putting in doubt the validity of the synonymous rate as a standard neutral proxy. While protein domains rarely undergo adaptive evolution, the SRV appears to play important role in optimizing the domain function at the level of DNA. In contrast, protein families are more likely to evolve by positive selection, but are less likely to exhibit SRV. Stronger SRV was detected in genes with stronger codon bias and tRNA reusage, those coding for proteins with larger number of interactions or forming larger number of structures, located in intracellular components and those involved in typically conserved complex processes and functions. Genes with extreme SRV show higher expression levels in nearly all tissues. This indicates that codon bias in a gene, which often correlates with gene expression, may often be a site-specific phenomenon regulating the speed of translation along the sequence, consistent with the co-translational folding hypothesis. Strikingly, genes with SRV were strongly overrepresented for metabolic pathways and those associated with several genetic diseases, particularly cancers and diabetes.

Evolution at protein ends: major contribution of alternative transcription initiation and termination to the transcriptome and proteome diversity in mammals

Alternative splicing (AS), alternative transcription initiation (ATI) and alternative transcription termination (ATT) create the extraordinary complexity of transcriptomes and make key contributions to the structural and functional diversity of mammalian proteomes. Analysis of mammalian genomic and transcriptomic data shows that contrary to the traditional view, the joint contribution of ATI and ATT to the transcriptome and proteome diversity is quantitatively greater than the contribution of AS. Although the mean numbers of protein-coding constitutive and alternative nucleotides in gene loci are nearly identical, their distribution along the transcripts is highly non-uniform. On average, coding exons in the variable 5' and 3' transcript ends that are created by ATI and ATT contain approximately four times more alternative nucleotides than core protein-coding regions that diversify exclusively via AS. Short upstream exons that encompass alternative 5'-untranslated regions and N-termini of proteins evolve under strong nucleotide-level selection whereas in 3'-terminal exons that encode protein C-termini, protein-level selection is significantly stronger. The groups of genes that are subject to ATI and ATT show major differences in biological roles, expression and selection patterns.

Rapid storage and retrieval of genomic intervals from a relational database system using nested containment lists

Efficient storage and retrieval of genomic annotations based on range intervals is necessary, given the amount of data produced by next-generation sequencing studies. The indexing strategies of relational database systems (such as MySQL) greatly inhibit their use in genomic annotation tasks. This has led to the development of stand-alone applications that are dependent on flat-file libraries. In this work, we introduce MyNCList, an implementation of the NCList data structure within a MySQL database. MyNCList enables the storage, update and rapid retrieval of genomic annotations from the convenience of a relational database system. Range-based annotations of 1 million variants are retrieved in under a minute, making this approach feasible for whole-genome annotation tasks.

Database URL:https://github.com/bushlab/mynclist

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Selective constraint on human pre-mRNA splicing by protein structural properties

Alternative splicing (AS) is a major mechanism of increasing proteome diversity in complex organisms. Different AS transcript isoforms may be translated into peptide sequences of significantly different lengths and amino acid compositions. One important question, then, is how AS is constrained by protein structural requirements while peptide sequences may be significantly changed in AS events. Here, we address this issue by examining whether the intactness of three-dimensional protein structural units (compact units in protein structures, namely protein units [PUs]) tends to be preserved in AS events in human. We show that PUs tend to occur in constitutively spliced exons and to overlap constitutive exon boundaries. Furthermore, when PUs are located at the boundaries between two alternatively spliced exons (ASEs), these neighboring ASEs tend to co-occur in different transcript isoforms. In addition, such PU-spanned ASE pairs tend to have a higher frequency of being included in transcript isoforms. ASE regions that overlap with PUs also have lower nonsynonymous-to-synonymous substitution rate ratios than those that do not overlap with PUs, indicating stronger negative selection pressure in PU-overlapped ASE regions. Of note, we show that PUs have protein domain- and structural orderness-independent effects on messenger RNA (mRNA) splicing. Overall, our results suggest that fine-scale protein structural requirements have significant influences on the splicing patterns of human mRNAs.

The function of introns

The intron–exon architecture of many eukaryotic genes raises the intriguing question of whether this unique organization serves any function, or is it simply a result of the spread of functionless introns in eukaryotic genomes. In this review, we show that introns in contemporary species fulfill a broad spectrum of functions, and are involved in virtually every step of mRNA processing. We propose that this great diversity of intronic functions supports the notion that introns were indeed selfish elements in early eukaryotes, but then independently gained numerous functions in different eukaryotic lineages. We suggest a novel criterion of evolutionary conservation, dubbed intron positional conservation, which can identify functional introns.

Identification of a new susceptibility variant for multiple sclerosis in OAS1 by population genetics analysis

Contrasting results have been reported concerning the association of a splice-site polymorphism (rs10774671) in OAS1 with multiple sclerosis (MS). We analysed two OAS1 regions encompassing alternatively spliced exons. While the region carrying the splice-site variant is neutrally evolving, a signature of long-standing balancing selection was observed across an alternative exon 7. Analysis of variants in this exon identified an insertion/deletion polymorphism (rs11352835, A/−) that originates predicted products with distinct C termini. This variant is located along the major branch of the haplotype genealogy, suggesting that it may represent the selection target. A case/control study for MS indicated that rs11352835 is associated with disease susceptibility (for an allelic model with the deleted allele predisposing to MS, OR 1.27, 95% CI 1.072–1.513, p = 0.010). No association was found between rs10774671 and MS. As the two SNPs are in linkage disequilibrium in Europeans, the previously reported association between rs10774671 and MS susceptibility might be driven by rs11352835, possibly explaining the contrasting results previously observed for the splice-site polymorphism. Thus, we describe a novel susceptibility variant for MS in OAS1 and show that population genetic analyses can be instrumental to the identification of selection targets and, consequently, of functional polymorphisms with an effect on phenotypic traits.

A molecular analysis of mutations at the complex dumpy locus in Drosophila melanogaster

The Drosophila dumpy gene consists of seventy eight coding exons and encodes a huge extracellular matrix protein containing large numbers of epidermal growth factor-like (EGF) modules and a novel module called dumpy (DPY). A molecular analysis of forty five mutations in the dumpy gene of Drosophila melanogaster was carried out. Mutations in this gene affect three phenotypes: wing shape, thoracic cuticular defects, and lethality. Most of the mutations were chemically induced in a single dumpy allele and were analyzed using a nuclease that cleaves single base pair mismatches in reannealed duplexes followed by dHPLC. Additionally, several spontaneous mutations were analyzed. Virtually all of the chemically induced mutations, except for several in a single exon, either generate nonsense codons or lesions that result in downstream stop codons in the reading frame. The remaining chemically induced mutations remove splice sites in the nascent dumpy message. We propose that the vast majority of nonsense mutations that affect all three basic dumpy phenotypes are in constitutive exons, whereas nonsense mutants that remove only one or two of the basic functions are in alternatively spliced exons. Evolutionary comparisons of the dumpy gene from seven Drosophila species show strong conservation of the 5′ ends of exons where mutants with partial dumpy function are found. In addition, reverse transcription PCR analyses reveal transcripts in which exons marked by nonsense mutations with partial dumpy function are absent.

An ESRP-regulated splicing programme is abrogated during the epithelial-mesenchymal transition

Alternative splicing achieves coordinated changes in post-transcriptional gene expression programmes through the activities of diverse RNA-binding proteins. Epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) are cell-type-specific regulators of transcripts that switch splicing during the epithelial-mesenchymal transition (EMT). To define a comprehensive programme of alternative splicing that is regulated during the EMT, we identified an extensive ESRP-regulated splicing network of hundreds of alternative splicing events within numerous genes with functions in cell-cell adhesion, polarity, and migration. Loss of this global ESRP-regulated epithelial splicing programme induces the phenotypic changes in cell morphology that are observed during the EMT. Components of this splicing signature provide novel molecular markers that can be used to characterize the EMT. Bioinformatics and experimental approaches revealed a high-affinity ESRP-binding motif and a predictive RNA map that governs their activity. This work establishes the ESRPs as coordinators of a complex alternative splicing network that adds an important post-transcriptional layer to the changes in gene expression that underlie epithelial-mesenchymal transitions during development and disease.

Evolution of alternative splicing in primate brain transcriptomes

Alternative splicing is a predominant form of gene regulation in higher eukaryotes. The evolution of alternative splicing provides an important mechanism for the acquisition of novel gene functions. In this work, we carried out a genome-wide phylogenetic survey of lineage-specific splicing patterns in the primate brain, via high-density exon junction array profiling of brain transcriptomes of humans, chimpanzees and rhesus macaques. We identified 509 genes showing splicing differences among these species. RT-PCR analysis of 40 exons confirmed the predicted splicing evolution of 33 exons. Of these 33 exons, outgroup analysis using rhesus macaques confirmed 13 exons with human-specific increase or decrease in transcript inclusion levels after humans diverged from chimpanzees. Some of the human-specific brain splicing patterns disrupt domains critical for protein-protein interactions, and some modulate translational efficiency of their host genes. Strikingly, for exons showing splicing differences across species, we observed a significant increase in the rate of silent substitutions within exons, coupled with accelerated sequence divergence in flanking introns. This indicates that evolution of cis-regulatory signals is a major contributor to the emergence of human-specific splicing patterns. In one gene (MAGOH), using minigene reporter assays, we demonstrated that the combination of two human-specific cis-sequence changes created its human-specific splicing pattern. Together, our data reveal widespread human-specific changes of alternative splicing in the brain and suggest an important role of splicing in the evolution of neuronal gene regulation and functions.