Regulation of the kinase activity of the MIK GCK-like MAP4K by alternative splicing

Alternative splicing of introns is essential to ensure the complexity of mammalian genome functions. In particular, the generation of a high number of different isoforms by alternative splicing is an important characteristic of genes coding for signalling proteins such as mitogen activated protein kinases (MAPKs). This is thought to allow these proteins to transduce multiple stimuli in a highly regulated manner. Plant genes are also subjected to alternative splicing. Nevertheless, clear examples of the functional consequences of this phenomenon are still scarce in plants. MIK is a maize gene coding for a GCK-like MAP4K that can be activated by interaction with maize atypical receptor kinase (MARK), an atypical receptor kinase. Here we show that MIK is subjected to alternative splicing. Expression of MIK leads to, at least, 4 different mature mRNAs that accumulate with particular expression profiles during maize development. Our results show that the polypeptides encoded by the different MIK mRNAs display different kinase activity and are differentially activated by interaction with the MARK receptor. Two MIK isoforms display constitutive kinase activity, one isoform is inactive but can be activated by MARK, and the fourth MIK isoform is inactive and cannot be activated by MARK. Our results constitute a clear example of the biochemical consequences of alternative splicing in plants. The selective conservation during evolution of the intron–exon structure of the region coding for the regulator domain of MIK, as well as the maintenance in maize, rice and Arabidopsis of the alternative splicing of some of these introns, are strong indications of its functional importance.     

Alternative splicing regulates targeting of malate dehydrogenase in Yarrowia lipolytica

Alternative pre-mRNA splicing is a major mechanism contributing to the proteome complexity of most eukaryotes, especially mammals. In less complex organisms, such as yeasts, the numbers of genes that contain introns are low and cases of alternative splicing (AS) with functional implications are rare. We report the first case of AS with functional consequences in the yeast Yarrowia lipolytica. The splicing pattern was found to govern the cellular localization of malate dehydrogenase, an enzyme of the central carbon metabolism. This ubiquitous enzyme is involved in the tricarboxylic acid cycle in mitochondria and in the glyoxylate cycle, which takes place in peroxisomes and the cytosol. In Saccharomyces cerevisiae, three genes encode three compartment-specific enzymes. In contrast, only two genes exist in Y. lipolytica. One gene (YlMDH1, YALI0D16753g) encodes a predicted mitochondrial protein, whereas the second gene (YlMDH2, YALI0E14190g) generates the cytosolic and peroxisomal forms through the alternative use of two 3'-splice sites in the second intron. Both splicing variants were detected in cDNA libraries obtained from cells grown under different conditions. Mutants expressing the individual YlMdh2p isoforms tagged with fluorescent proteins confirmed that they localized to either the cytosolic or the peroxisomal compartment.     

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.     

Evolution of alternative splicing: deletions,insertions and origin of functional parts of proteins from intron sequences

Alternative splicing is thought to be a major source of functional diversity in animal proteins. We analyzed the evolutionary conservation of proteins encoded by alternatively spliced genes and predicted the ancestral state for 73 cases of alternative splicing (25 insertions and 48 deletions). The amino acid sequences of most of the inserts in proteins produced by alternative splicing are as conserved as the surrounding sequences. Thus, alternative splicing often creates novel isoforms by the insertion of new, functional protein sequences that probably originated from noncoding sequences of introns.     

Integration of RNA processing and expression level control modulates the function of the Drosophila Hox gene Ultrabithorax during adult development

Although most metazoan genes undergo alternative splicing, the functional relevance of the majority of alternative splicing products is still unknown. Here we explore this problem in the Drosophila Hox gene Ultrabithorax (Ubx). Ubx produces a family of six protein isoforms through alternative splicing. To investigate the functional specificity of the Ubx isoforms, we studied their role during the formation of the Drosophila halteres, small dorsal appendages that are essential for normal flight. Our work shows that isoform Ia, which is encoded by all Ubx exons, is more efficient than isoform IVa, which lacks the amino acids coded by two small exons, in controlling haltere development and regulating Ubx downstream targets. However, our experiments also demonstrate that the functional differences among the Ubx isoforms can be compensated for by increasing the expression levels of the less efficient form. The analysis of the DNA-binding profiles of Ubx isoforms to a natural Ubx target, spalt, shows no major differences in isoform DNA-binding activities, suggesting that alternative splicing might primarily affect the regulatory capacity of the isoforms rather than their DNA-binding patterns. Our results suggest that to obtain distinct functional outputs during normal development genes must integrate the generation of qualitative differences by alternative splicing to quantitative processes affecting isoform protein expression levels.     

Analysis and design of RNA sequencing experiments for identifying isoform regulation

Through alternative splicing, most human genes express multiple isoforms that often differ in function. To infer isoform regulation from high-throughput sequencing of cDNA fragments (RNA-seq), we developed the mixture-of-isoforms (MISO) model, a statistical model that estimates expression of alternatively spliced exons and isoforms and assesses confidence in these estimates. Incorporation of mRNA fragment length distribution in paired-end RNA-seq greatly improved estimation of alternative-splicing levels. MISO also detects differentially regulated exons or isoforms. Application of MISO implicated the RNA splicing factor hnRNP H1 in the regulation of alternative cleavage and polyadenylation, a role that was supported by UV cross-linking-immunoprecipitation sequencing (CLIP-seq) analysis in human cells. Our results provide a probabilistic framework for RNA-seq analysis, give functional insights into pre-mRNA processing and yield guidelines for the optimal design of RNA-seq experiments for studies of gene and isoform expression.     

Identification of alternatively spliced dab1 isoforms in zebrafish

We have investigated the genomic organization, the occurrence of alternative splicing and the differential expression of the zebrafish disabled1 (dab1) gene. Dab1 is a key effector of the Reelin pathway, which regulates neuronal migration during brain development in vertebrates. The coding region of the zebrafish dab1 gene spans over 600 kb of genomic DNA and is composed of 15 exons. Alternative splicing in a region enriched for tyrosine residues generates at least three different isoforms. These isoforms are developmentally regulated and show differential tissue expression. Comparison with mouse and human data shows an overall conservation of the genomic organization with different alternative splicing events generating species-specific isoforms. Because these alternative splicing events give rise to isoforms with different numbers of phosphorylateable tyrosines, we speculate that alternative splicing of the dab1 gene in zebrafish and in other vertebrates regulates the nature of the cellular response to the Reelin signal.Electronic supplementary material Supplementary material is available for this article at and accessible for authorised users.     

Recognition of unknown conserved alternatively spliced exons

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.