Evolution of alternative splicing regulation: changes in predicted exonic splicing regulators are not associated with changes in alternative splicing levels in primates

Towards understanding pre-mRNA splicing mechanisms and the role of SR proteins

Alternative pre-mRNA splicing provides a source of vast protein diversity by removing non-coding sequences (introns) and accurately linking different exonic regions in the correct reading frame. The regulation of alternative splicing is essential for various cellular functions in both pathological and physiological conditions. In eukaryotic cells, this process is commonly used to increase proteomic diversity and to control gene expression either co- or post-transcriptionally. Alternative splicing occurs within a megadalton-sized, multi-component machine consisting of RNA and proteins; during the splicing process, this complex undergoes dynamic changes via RNA-RNA, protein-protein and RNA-protein interactions. Co-transcriptional splicing functionally integrates the transcriptional machinery, thereby enabling the two processes to influence one another, whereas post-transcriptional splicing facilitates the coupling of RNA splicing with post-splicing events. This review addresses the structural aspects of spliceosomes and the mechanistic implications of their stepwise assembly on the regulation of pre-mRNA splicing. Moreover, the role of phosphorylation-based, signal-induced changes in the regulation of the splicing process is demonstrated.

Cold-dependent alternative splicing of a Jumonji C domain-containing gene MtJMJC5 in Medicago truncatula

Histone methylation is an epigenetic modification mechanism that regulates gene expression in eukaryotic cells. Jumonji C domain-containing demethylases are involved in removal of methyl groups at lysine or arginine residues. The JmjC domain-only member, JMJ30/JMJD5 of Arabidopsis, is a component of the plant circadian clock. Although some plant circadian clock genes undergo alternative splicing in response to external cues, there is no evidence that JMJ30/JMJD5 is regulated by alternative splicing. In this study, the expression of an Arabidopsis JMJ30/JMJD5 ortholog in Medicago truncatula, MtJMJC5, in response to circadian clock and abiotic stresses were characterized. The results showed that MtJMJC5 oscillates with a circadian rhythm, and undergoes cold specifically induced alternative splicing. The cold-induced alternative splicing could be reversed after ambient temperature returning to the normal. Sequencing results revealed four alternative splicing RNA isoforms including a full-length authentic protein encoding variant, and three premature termination condon-containing variants due to alternative 3' splice sites at the first and second intron. Under cold treatment, the variants that share a common 3' alternative splicing site at the second intron were intensively up-regulated while the authentic protein encoding variant and the premature termination condon-containing variant only undergoing a 3' alternative splicing at the first intron were down regulated. Although all the premature termination condon-harboring alternative splicing variants were sensitive to nonsense-mediated decay, the premature termination codon-harboring alternative splicing variants sharing the 3' alternative splicing site at the second intron showed less sensitivity than the one only containing the 3' alternative slicing site at the first intron under cold treatment. These results suggest that the cold-dependent alternative splicing of MtJMJC5 is likely a species or genus-specific mechanism of gene expression regulation on RNA levels, and might play a role in epigenetic regulation of the link between the circadian clock and ambient temperature fluctuation in Medicago.

Evolution of splicing regulatory networks in Drosophila

The proteome expanding effects of alternative pre-mRNA splicing have had a profound impact on eukaryotic evolution. The events that create this diversity can be placed into four major classes: exon skipping, intron retention, alternative 5′ splice sites, and alternative 3′ splice sites. Although the regulatory mechanisms and evolutionary pressures among alternative splicing classes clearly differ, how these differences affect the evolution of splicing regulation remains poorly characterized. We used RNA-seq to investigate splicing differences in D. simulans, D. sechellia, and three strains of D. melanogaster. Regulation of exon skipping and tandem alternative 3′ splice sites (NAGNAGs) were more divergent than other splicing classes. Splicing regulation was most divergent in frame-preserving events and events in noncoding regions. We further determined the contributions of cis- and trans-acting changes in splicing regulatory networks by comparing allele-specific splicing in F1 interspecific hybrids, because differences in allele-specific splicing reflect changes in cis-regulatory element activity. We find that species-specific differences in intron retention and alternative splice site usage are primarily attributable to changes in cis-regulatory elements (median ∼80% cis), whereas species-specific exon skipping differences are driven by both cis- and trans-regulatory divergence (median ∼50% cis). These results help define the mechanisms and constraints that influence splicing regulatory evolution and show that networks regulating the four major classes of alternative splicing diverge through different genetic mechanisms. We propose a model in which differences in regulatory network architecture among classes of alternative splicing affect the evolution of splicing regulation.

Genome-wide detection of condition-sensitive alternative splicing in Arabidopsis roots

Iron (Fe) deficiency is a world-wide nutritional disorder in both plants and humans, resulting from its restricted bioavailability for plants and, subsequently, low Fe concentration in edible plant parts. Plants have evolved sophisticated mechanisms to alleviate Fe deficiency, with the aim of recalibrating metabolic fluxes and maintaining cellular Fe homeostasis. To analyze condition-sensitive changes in precursor mRNA (pre-mRNA) splicing pattern, we mapped the transcriptome of Fe-deficient and Fe-sufficient Arabidopsis (Arabidopsis thaliana) roots using the RNA sequencing technology and a newly developed software toolbox, the Read Analysis & Comparison Kit in Java (RACKJ). In alternatively spliced genes, stress-related Gene Ontology categories were overrepresented, while housekeeping cellular functions were mainly transcriptionally controlled. Fe deficiency increased the complexity of the splicing pattern and triggered the differential alternative splicing of 313 genes, the majority of which had differentially retained introns. Several genes with important functions in Fe acquisition and homeostasis were both differentially expressed and differentially alternatively spliced upon Fe deficiency, indicating a complex regulation of gene activity in Fe-deficient conditions. A comparison with a data set for phosphate-deficient plants suggests that changes in splicing patterns are nutrient specific and not or not chiefly caused by stochastic fluctuations. In sum, our analysis identified extensive posttranscriptional control, biasing the abundance and activity of proteins in a condition-dependent manner. The production of a mixture of functional and nonfunctional transcripts may provide a means to fine-tune the abundance of transcripts with critical importance in cellular Fe homeostasis. It is assumed that differential gene expression and nutrient deficiency-induced changes in pre-mRNA splicing represent parallel, but potentially interacting, regulatory mechanisms.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.

Histone deacetylase activity modulates alternative splicing

There is increasing evidence to suggest that splicing decisions are largely made when the nascent RNA is still associated with chromatin. Here we demonstrate that activity of histone deacetylases (HDACs) influences splice site selection. Using splicing-sensitive microarrays, we identified ∼700 genes whose splicing was altered after HDAC inhibition. We provided evidence that HDAC inhibition induced histone H4 acetylation and increased RNA Polymerase II (Pol II) processivity along an alternatively spliced element. In addition, HDAC inhibition reduced co-transcriptional association of the splicing regulator SRp40 with the target fibronectin exon. We further showed that the depletion of HDAC1 had similar effect on fibronectin alternative splicing as global HDAC inhibition. Importantly, this effect was reversed upon expression of mouse HDAC1 but not a catalytically inactive mutant. These results provide a molecular insight into a complex modulation of splicing by HDACs and chromatin modifications.

TIA1 prevents skipping of a critical exon associated with spinal muscular atrophy

Prevention of skipping of exon 7 during pre-mRNA splicing of Survival Motor Neuron 2 (SMN2) holds the promise for cure of spinal muscular atrophy (SMA), a leading genetic cause of infant mortality. Here, we report T-cell-restricted intracellular antigen 1 (TIA1) and TIA1-related (TIAR) proteins as intron-associated positive regulators of SMN2 exon 7 splicing. We show that TIA1/TIAR stimulate exon recognition in an entirely novel context in which intronic U-rich motifs are separated from the 5' splice site by overlapping inhibitory elements. TIA1 and TIAR are modular proteins with three N-terminal RNA recognition motifs (RRMs) and a C-terminal glutamine-rich (Q-rich) domain. Our results reveal that any one RRM in combination with a Q domain is necessary and sufficient for TIA1-associated regulation of SMN2 exon 7 splicing in vivo. We also show that increased expression of TIA1 counteracts the inhibitory effect of polypyrimidine tract binding protein, a ubiquitously expressed factor recently implicated in regulation of SMN exon 7 splicing. Our findings expand the scope of TIA1/TIAR in genome-wide regulation of alternative splicing under normal and pathological conditions.

New isoforms of human mitochondrial transcription factor A detected in normal and tumoral cells

Novel alternatively spliced variants of the human mitochondrial transcription factor A predicted by the computational tool ASPic were experimentally validated in different normal and tumoral human tissues by RT-PCR and DNA sequencing. The comparison between the 5'UTR length and the distribution of the different transcripts showed that the transcripts with the shortest 5'UTR are present in all the investigated tissues, while the longest 5'UTR seems to be related to tissue-specificity. Studies about the localization and function of the most widely diffuse alternative isoform Tr6 were carried out.

Human brain evolution: harnessing the genomics (r)evolution to link genes, cognition, and behavior

The evolution of the human brain has resulted in numerous specialized features including higher cognitive processes such as language. Knowledge of whole-genome sequence and structural variation via high-throughput sequencing technology provides an unprecedented opportunity to view human evolution at high resolution. However, phenotype discovery is a critical component of these endeavors and the use of nontraditional model organisms will also be critical for piecing together a complete picture. Ultimately, the union of developmental studies of the brain with studies of unique phenotypes in a myriad of species will result in a more thorough model of the groundwork the human brain was built upon. Furthermore, these integrative approaches should provide important insights into human diseases.

Mutation analysis and characterization of alternative splice variants of the Wilson disease gene ATP7B

Wilson disease is a copper metabolism disorder caused by mutations in ATP7B, a copper-transporting adenosine triphosphatase. A molecular diagnosis was performed on 135 patients with Wilson disease in Taiwan. We identified 36 different mutations, eight of which were novel: five missense mutations (Ser986Phe, Ile1348Asn, Gly1355Asp, Met1392Lys, and Ala1445Pro), one deletion (2810delT) in the coding region, and two nucleotide substitutions (-133A→C and -215A→T) in the promoter region. These mutations were not observed in 100 control subjects and reduced the activity of the mutated protein by at least 50% when compared with wild-type ATP7B. In addition to exon 8, our data indicate another mutation hotspot in exon 12 where 9.62% of all mutations occurred. An alternative splice variant of ATP7B lacking exon 12 was observed in one patient who had a homozygous 2810delT mutation and very mild clinical symptoms. Clinical examination and functional characterization of alternative splice variants of ATP7B lacking exon 12 showed that they retained 80% of their biological activity. The 2810delT mutation increased the expression of these variants, which may have explained the mild symptoms in the patient with the 2810delT mutation. We also discovered that treating liver cancer cells with a Na(+)/H(+) exchanger inhibitor, 5-(N-ethyl-N-isopropyl)-amiloride, significantly enhanced the expression of the alternative splice variant of ATP7B lacking exon 12.

CONCLUSION

This study suggests a novel therapeutic strategy for patients with mutations in exon 12.

The adaptive significance of unproductive alternative splicing in primates

Alternative gene splicing is pervasive in metazoa, particularly in humans, where the majority of genes generate splice variant transcripts. Characterizing the biological significance of alternative transcripts is methodologically difficult since it is impractical to assess thousands of splice variants as to whether they actually encode proteins, whether these proteins are functional, or whether transcripts have a function independent of protein synthesis. Consequently, to elucidate the functional significance of splice variants and to investigate mechanisms underlying the fidelity of mRNA splicing, we used an indirect approach based on analyzing the evolutionary conservation of splice variants among species. Using DNA polymerase β as an indicator locus, we cloned and characterized the types and frequencies of transcripts generated in primary cell lines of five primate species. Overall, we found that in addition to the canonical DNA polymerase β transcript, there were 25 alternative transcripts generated, most containing premature terminating codons. We used a statistical method borrowed from community ecology to show that there is significant diversity and little conservation in alternative splicing patterns among species, despite high sequence similarity in the underlying genomic (exonic) sequences. However, the frequency of alternative splicing at this locus correlates well with life history parameters such as the maximal longevity of each species, indicating that the alternative splicing of unproductive splice variants may have adaptive significance, even if the specific RNA transcripts themselves have no function. These results demonstrate the validity of the phylogenetic conservation approach in elucidating the biological significance of alternative splicing.

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

High-throughput methods such as EST sequencing, microarrays and deep sequencing have identified large numbers of alternative splicing (AS) events, but studies have shown that only a subset of these may be functional. Here we report a sensitive bioinformatics approach that identifies exons with evidence of a strong RNA selection pressure ratio (RSPR)--i.e., evolutionary selection against mutations that change only the mRNA sequence while leaving the protein sequence unchanged--measured across an entire evolutionary family, which greatly amplifies its predictive power. Using the UCSC 28 vertebrate genome alignment, this approach correctly predicted half to three-quarters of AS exons that are known binding targets of the NOVA splicing regulatory factor, and predicted 345 strongly selected alternative splicing events in human, and 262 in mouse. These predictions were strongly validated by several experimental criteria of functional AS such as independent detection of the same AS event in other species, reading frame-preservation, and experimental evidence of tissue-specific regulation: 75% (15/20) of a sample of high-RSPR exons displayed tissue specific regulation in a panel of ten tissues, vs. only 20% (4/20) among a sample of low-RSPR exons. These data suggest that RSPR can identify exons with functionally important splicing regulation, and provides biologists with a dataset of over 600 such exons. We present several case studies, including both well-studied examples (GRIN1) and novel examples (EXOC7). These data also show that RSPR strongly outperforms other approaches such as standard sequence conservation (which fails to distinguish amino acid selection pressure from RNA selection pressure), or pairwise genome comparison (which lacks adequate statistical power for predicting individual exons).

Functional properties and evolutionary splicing constraints on a composite exonic regulatory element of splicing in CFTR exon 12

In general, splicing regulatory elements are defined as Enhancers or Silencers depending on their positive or negative effect upon exon inclusion. Often, these sequences are usually present separate from each other in exonic/intronic sequences. The Composite Exonic Splicing Regulatory Elements (CERES) represent an extreme physical overlap of enhancer/silencer activity. As a result, when CERES elements are mutated the consequences on the splicing process are difficult to predict. Here, we show that the functional activity of the CERES2 sequence in CFTR exon 12 is regulated by the binding, in very close proximity to each other, of several SR and hnRNP proteins. Moreover, our results show that practically the entire exon 12 sequence context participate in its definition. The consequences of this situation can be observed at the evolutionary level by comparing changes in conservation of different splicing elements in different species. In conclusion, our study highlights how it is increasingly difficult to define many exonic sequences by simply breaking them down in isolated enhancer/silencer or even neutral elements. The real picture is close to one of continuous competition between positive and negative factors where affinity for the target sequences and other dynamic factors decide the inclusion or exclusion of the exon.

Divergence of exonic splicing elements after gene duplication and the impact on gene structures

An analysis of human exonic splicing elements in duplicated genes reveals their important role in the generation of new gene structures.

Background

The origin of new genes and their contribution to functional novelty has been the subject of considerable interest. There has been much progress in understanding the mechanisms by which new genes originate. Here we examine a novel way that new gene structures could originate, namely through the evolution of new alternative splicing isoforms after gene duplication.

Results

We studied the divergence of exonic splicing enhancers and silencers after gene duplication and the contributions of such divergence to the generation of new splicing isoforms. We found that exonic splicing enhancers and exonic splicing silencers diverge especially fast shortly after gene duplication. About 10% and 5% of paralogous exons undergo significantly asymmetric evolution of exonic splicing enhancers and silencers, respectively. When compared to pre-duplication ancestors, we found that there is a significant overall loss of exonic splicing enhancers and the magnitude increases with duplication age. Detailed examination reveals net gains and losses of exonic splicing enhancers and silencers in different copies and paralog clusters after gene duplication. Furthermore, we found that exonic splicing enhancer and silencer changes are mainly caused by synonymous mutations, though nonsynonymous changes also contribute. Finally, we found that exonic splicing enhancer and silencer divergence results in exon splicing state transitions (from constitutive to alternative or vice versa), and that the proportion of paralogous exon pairs with different splicing states also increases over time, consistent with previous predictions.

Conclusions

Our results suggest that exonic splicing enhancer and silencer changes after gene duplication have important roles in alternative splicing divergence and that these changes contribute to the generation of new gene structures.