Creation and disruption of protein features by alternative splicing -- a novel mechanism to modulate function

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high-throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing Translation > Regulation RNA Evolution and Genomics > Computational Analyses of RNA.

The WRKY transcription factor family in cowpea: Genomic characterization and transcriptomic profiling under root dehydration

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most tolerant legume crops to drought and salt stresses. WRKY transcription factor (TF) family members stand out among plant transcriptional regulators related to abiotic stress tolerance. However, little information is currently available on the expression of the cowpea WRKY gene family (VuWRKY) in response to water deficit. Thus, we analyzed genomic and transcriptomic data from cowpea to identify VuWRKY members and characterize their structure and transcriptional response under root dehydration stress. Ninety-two complete VuWRKY genes were found in the cowpea genome based on their domain characteristics. They were clustered into three groups: I (15 members), II (58), and III (16), while three genes were unclassified. Domain analysis of the encoded proteins identified four major variants of the conserved heptapeptide motif WRKYGQK. In silico analysis of VuWRKY gene promoters identified eight candidate binding motifs of cis-regulatory elements, regulated mainly by six TF families associated with abiotic stress responses. Ninety-seven VuWRKY modulated splicing variants associated with 55 VuWRKY genes were identified via RNA-Seq analysis available at the Cowpea Genomics Consortium (CpGC) database. qPCR analyses showed that 22 genes are induced under root dehydration, with VuWRKY18, 21, and 75 exhibiting the most significant induction levels. Given their central role in activating signal transduction cascades in abiotic stress response, the data provide a foundation for the targeted modification of specific VuWRKY family members to improve drought tolerance in this important climate-resilient legume in the developing world and beyond.

Transcriptome analysis reveals the molecular mechanisms of heterosis on thermal resistance in hybrid abalone

BACKGROUND

Heterosis has been exploited for decades in different animals and crops due to it resulting in dramatic increases in yield and adaptability. Hybridization is a classical breeding method that can effectively improve the genetic characteristics of organisms through heterosis. Abalone has become an increasingly economically important aquaculture resource with high commercial value. However, due to changing climate, abalone is now facing serious threats of high temperature in summer. Interspecific hybrid abalone (Haliotis gigantea ♀ × H. discus hannai ♂, SD) has been cultured at large scale in southern China and has been shown high survival rates under heat stress in summer. Therefore, SD has become a good model material for heterosis research, but the molecular basis of heterosis remains elusive.

RESULTS

Heterosis in thermal tolerance of SD was verified through Arrhenius break temperatures (ABT) of cardiac performance in this study. Then RNA-Sequencing was conducted to obtain gene expression patterns and alternative splicing events at control temperature (20 °C) and heat stress temperature (30 °C). A total of 356 (317 genes), 476 (435genes), and 876 (726 genes) significantly diverged alternative splicing events were identified in H. discus hannai (DD), H. gigantea (SS), and SD in response to heat stress, respectively. In the heat stress groups, 93.37% (20,512 of 21,969) of the expressed genes showed non-additive expression patterns, and over-dominance expression patterns of genes account for the highest proportion (40.15%). KEGG pathway enrichment analysis showed that the overlapping genes among common DEGs and NAGs were significantly enriched in protein processing in the endoplasmic reticulum, mitophagy, and NF-κB signaling pathway. In addition, we found that among these overlap genes, 39 genes had undergone alternative splicing events in SD. These pathways and genes may play an important role in the thermal resistance of hybrid abalone.

CONCLUSION

More alternative splicing events and non-additive expressed genes were detected in hybrid under heat stress and this may contribute to its thermal heterosis. These results might provide clues as to how hybrid abalone has a better physiological regulation ability than its parents under heat stress, to increase our understanding of heterosis in abalone.

An isoform of the plastid RNA polymerase-associated protein FSD3 negatively regulates chloroplast development

BACKGROUND

Plastid-encoded RNA polymerase (PEP) plays an essential role in chloroplast development by governing the expression of genes involved in photosynthesis. At least 12 PEP-associated proteins (PAPs), including FSD3/PAP4, regulate PEP activity and chloroplast development by modulating formation of the PEP complex.

RESULTS

In this study, we identified FSD3S, a splicing variant of FSD3; the FSD3 and FSD3S transcripts encode proteins with identical N-termini, but different C-termini. Characterization of FSD3 and FSD3S proteins showed that the C-terminal region of FSD3S contains a transmembrane domain, which promotes FSD3S localization to the chloroplast membrane but not to nucleoids, in contrast to FSD3, which localizes to the chloroplast nucleoid. We also found that overexpression of FSD3S negatively affects photosynthetic activity and chloroplast development by reducing expression of genes involved in photosynthesis. In addition, FSD3S failed to complement the chloroplast developmental defects in the fsd3 mutant.

CONCLUSION

These results suggest FSD3 and FSD3S, with their distinct localization patterns, have different functions in chloroplast development, and FSD3S negatively regulates expression of PEP-dependent chloroplast genes, and development of chloroplasts.

Search for potential reading frameshifts in cds from Arabidopsis thaliana and other genomes

A new mathematical method for potential reading frameshift detection in protein-coding sequences (cds) was developed. The algorithm is adjusted to the triplet periodicity of each analysed sequence using dynamic programming and a genetic algorithm. This does not require any preliminary training. Using the developed method, cds from the Arabidopsis thaliana genome were analysed. In total, the algorithm found 9,930 sequences containing one or more potential reading frameshift(s). This is ∼21% of all analysed sequences of the genome. The Type I and Type II error rates were estimated as 11% and 30%, respectively. Similar results were obtained for the genomes of Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Rattus norvegicus and Xenopus tropicalis. Also, the developed algorithm was tested on 17 bacterial genomes. We compared our results with the previously obtained data on the search for potential reading frameshifts in these genomes. This study discussed the possibility that the reading frameshift seems like a relatively frequently encountered mutation; and this mutation could participate in the creation of new genes and proteins.

Fluorescent Protein-Based Quantification of Alternative Splicing of a Target Cassette Exon in Mammalian Cells

Alternative splicing is an important mechanism of regulation of gene expression and expansion of proteome complexity. Recently we developed a new fluorescence reporter for quantitative analysis of alternative splicing of a target cassette exon in live cells (Gurskaya et al., 2012). It consists of a specially designed minigene encoding red and green fluorescent proteins (Katushka and TagGFP2) and a fragment of the target gene between them. Skipping or inclusion of the alternative exon induces a frameshift; ie, alternative exon length must not be a multiple of 3. Finally, red and green fluorescence intensities of cells expressing this reporter are used to estimate the percentage of alternative (exon-skipped) and normal (exon-retained) transcripts. Here, we provide a detailed description of design and application of the fluorescence reporter of a target alternative exon splicing in mammalian cell lines.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1914-5) contains supplementary material, which is available to authorized users.

Evolution of SH2 domains and phosphotyrosine signalling networks

Src homology 2 (SH2) domains mediate selective protein-protein interactions with tyrosine phosphorylated proteins, and in doing so define specificity of phosphotyrosine (pTyr) signalling networks. SH2 domains and protein-tyrosine phosphatases expand alongside protein-tyrosine kinases (PTKs) to coordinate cellular and organismal complexity in the evolution of the unikont branch of the eukaryotes. Examination of conserved families of PTKs and SH2 domain proteins provides fiduciary marks that trace the evolutionary landscape for the development of complex cellular systems in the proto-metazoan and metazoan lineages. The evolutionary provenance of conserved SH2 and PTK families reveals the mechanisms by which diversity is achieved through adaptations in tissue-specific gene transcription, altered ligand binding, insertions of linear motifs and the gain or loss of domains following gene duplication. We discuss mechanisms by which pTyr-mediated signalling networks evolve through the development of novel and expanded families of SH2 domain proteins and the elaboration of connections between pTyr-signalling proteins. These changes underlie the variety of general and specific signalling networks that give rise to tissue-specific functions and increasingly complex developmental programmes. Examination of SH2 domains from an evolutionary perspective provides insight into the process by which evolutionary expansion and modification of molecular protein interaction domain proteins permits the development of novel protein-interaction networks and accommodates adaptation of signalling networks.

On programmed ribosomal frameshifting: the alternative proteomes

Frameshifting results from two main mechanisms: genomic insertions or deletions (indels) or programmed ribosomal frameshifting. Whereas indels can disrupt normal protein function, programmed ribosomal frameshifting can result in dual-coding genes, each of which can produce multiple functional products. Here, I summarize technical advances that have made it possible to identify programmed ribosomal frameshifting events in a systematic way. The results of these studies suggest that such frameshifting occurs in all genomes, and I will discuss methods that could help characterize the resulting alternative proteomes.

Alternatively spliced domains interact to regulate BK potassium channel gating

Most human genes contain multiple alternative splice sites believed to extend the complexity and diversity of the proteome. However, little is known about how interactions among alternative exons regulate protein function. We used the Caenorhabditis elegans slo-1 large-conductance calcium and voltage-activated potassium (BK) channel gene, which contains three alternative splice sites (A, B, and C) and encodes at least 12 splice variants, to investigate the functional consequences of alternative splicing. These splice sites enable the insertion of exons encoding part of the regulator of K(+) conductance (RCK)1 Ca(2+) coordination domain (exons A1 and A2) and portions of the RCK1-RCK2 linker (exons B0, B1, B2, C0, and C1). Exons A1 and A2 are used in a mutually exclusive manner and are 67% identical. The other exons can extend the RCK1-RCK2 linker by up to 41 residues. Electrophysiological recordings of all isoforms show that the A1 and A2 exons regulate activation kinetics and Ca(2+) sensitivity, but only if alternate exons are inserted at site B or C. Thus, RCK1 interacts with the RCK1-RCK2 linker, and the effect of exon variation on gating depends on the combination of alternate exons present in each isoform.

Alternative splicing of transcription factors' genes: beyond the increase of proteome diversity

Functional modification of transcription regulators may lead to developmental changes and phenotypical differences between species. In this work, we study the influence of alternative splicing on transcription factors in human and mouse. Our results show that the impact of alternative splicing on transcription factors is similar in both species, meaning that the ways to increase variability should also be similar. However, when looking at the expression patterns of transcription factors, we observe that they tend to diverge regardless of the role of alternative splicing. Finally, we hypothesise that transcription regulation of alternatively spliced transcription factors could play an important role in the phenotypical differences between species, without discarding other phenomena or functional families.

Genome and proteome annotation: organization, interpretation and integration

Recent years have seen a huge increase in the generation of genomic and proteomic data. This has been due to improvements in current biological methodologies, the development of new experimental techniques and the use of computers as support tools. All these raw data are useless if they cannot be properly analysed, annotated, stored and displayed. Consequently, a vast number of resources have been created to present the data to the wider community. Annotation tools and databases provide the means to disseminate these data and to comprehend their biological importance. This review examines the various aspects of annotation: type, methodology and availability. Moreover, it puts a special interest on novel annotation fields, such as that of phenotypes, and highlights the recent efforts focused on the integrating annotations.

Phylogenetically widespread alternative splicing at unusual GYNGYN donors

BACKGROUND

Splice donor sites have a highly conserved GT or GC dinucleotide and an extended intronic consensus sequence GTRAGT that reflects the sequence complementarity to the U1 snRNA. Here, we focus on unusual donor sites with the motif GYNGYN (Y stands for C or T; N stands for A, C, G, or T).

RESULTS

While only one GY functions as a splice donor for the majority of these splice sites in human, we provide computational and experimental evidence that 110 (1.3%) allow alternative splicing at both GY donors. The resulting splice forms differ in only three nucleotides, which results mostly in the insertion/deletion of one amino acid. However, we also report the insertion of a stop codon in four cases. Investigating what distinguishes alternatively from not alternatively spliced GYNGYN donors, we found differences in the binding to U1 snRNA, a strong correlation between U1 snRNA binding strength and the preferred donor, over-represented sequence motifs in the adjacent introns, and a higher conservation of the exonic and intronic flanks between human and mouse. Extending our genome-wide analysis to seven other eukaryotic species, we found alternatively spliced GYNGYN donors in all species from mouse to Caenorhabditis elegans and even in Arabidopsis thaliana. Experimental verification of a conserved GTAGTT donor of the STAT3 gene in human and mouse reveals a remarkably similar ratio of alternatively spliced transcripts in both species.

CONCLUSION

In contrast to alternative splicing in general, GYNGYN donors in addition to NAGNAG acceptors enable subtle protein variations.

Non-EST-based prediction of novel alternatively spliced cassette exons with cell signaling function in Caenorhabditis elegans and human

To better understand the complex role that alternative splicing plays in intracellular signaling, it is important to catalog the numerous splice variants involved in signal transduction. Therefore, we developed PASE (Prediction of Alternative Signaling Exons), a computational tool to identify novel alternative cassette exons that code for kinase phosphorylation or signaling protein-binding sites. We first applied PASE to the Caenorhabditis elegans genome. In this organism, our algorithm had an overall specificity of ≥76.4%, including 33 novel cassette exons that we experimentally verified. We then used PASE to analyze the human genome and made 804 predictions, of which 308 were found as alternative exons in the transcript database. We experimentally tested 384 of the remaining unobserved predictions and discovered 26 novel human exons for a total specificity of ≥41.5% in human. By using a test set of known alternatively spliced signaling exons, we determined that the sensitivity of PASE is ∼70%. GO term analysis revealed that our exon predictions were found in the introns of known signal transduction genes more often than expected by chance, indicating PASE enriches for splice variants that function in signaling pathways. Overall, PASE was able to uncover 59 novel alternative cassette exons in C. elegans and humans through a genome-wide ab initio prediction method that enriches for exons involved in signaling.

Genetic resources, genome mapping and evolutionary genomics of the pig (Sus scrofa)

The pig, a representative of the artiodactyla clade, is one of the first animals domesticated, and has become an important agriculture animal as one of the major human nutritional sources of animal based protein. The pig is also a valuable biomedical model organism for human health. The pig's importance to human health and nutrition is reflected in the decision to sequence its genome (3X). As an animal species with its wild ancestors present in the world, the pig provides a unique opportunity for tracing mammalian evolutionary history and defining signatures of selection resulting from both domestication and natural selection. Completion of the pig genome sequencing project will have significant impacts on both agriculture and human health. Following the pig whole genome sequence drafts, along with large-scale polymorphism data, it will be possible to conduct genome sweeps using association mapping, and identify signatures of selection. Here, we provide a description of the pig genome sequencing project and perspectives on utilizing genomic technologies to exploit pig genome evolution and the molecular basis for phenotypic traits for improving pig production and health.

Relating alternative splicing to proteome complexity and genome evolution

Prior to genomics, studies of alternative splicing primarily focused on the function and mechanism of alternative splicing in individual genes and exons. This has changed dramatically since the late 1990s. High-throughput genomics technologies, such as EST sequencing and microarrays designed to detect changes in splicing, led to genome-wide discoveries and quantification of alternative splicing in a wide range of species from human to Arabidopsis. Consensus estimates of AS frequency in the human genome grew from less than 5% in mid-1990s to as high as 60-74% now. The rapid growth in sequence and microarray data for alternative splicing has made it possible to look into the global impact of alternative splicing on protein function and evolution of genomes. In this chapter, we review recent research on alternative splicing's impact on proteomic complexity and its role in genome evolution.

Frequent appearance of novel protein-coding sequences by frameshift translation

Genomic duplication, followed by divergence, contributes to organismal evolution. Several mechanisms, such as exon shuffling and alternative splicing, are responsible for novel gene functions, but they generate homologous domains and do not usually lead to drastic innovation. Major novelties can potentially be introduced by frameshift mutations and this idea can explain the creation of novel proteins. Here, we employ a strategy using simulated protein sequences and identify 470 human and 108 mouse frameshift events that originate new gene segments. No obvious interspecies overlap was observed, suggesting high rates of acquisition of evolutionary events. This inference is supported by a deficiency of TpA dinucleotides in the protein-coding sequences, which decreases the occurrence of translational termination, even on the complementary strand. Increased usage of the TGA codon as the termination signal in newer genes also supports our inference. This suggests that tolerated frameshift changes are a prevalent mechanism for the rapid emergence of new genes and that protein-coding sequences can be derived from existing or ancestral exons rather than from events that result in noncoding sequences becoming exons.

Identification and evolutionary analysis of novel exons and alternative splicing events using cross-species EST-to-genome comparisons in human, mouse and rat

Background

Alternative splicing (AS) is important for evolution and major biological functions in complex organisms. However, the extent of AS in mammals other than human and mouse is largely unknown, making it difficult to study AS evolution in mammals and its biomedical implications.

Results

Here we describe a cross-species EST-to-genome comparison algorithm (ENACE) that can identify novel exons for EST-scanty species and distinguish conserved and lineage-specific exons. The identified exons represent not only novel exons but also evolutionarily meaningful AS events that are not previously annotated. A genome-wide AS analysis in human, mouse and rat using ENACE reveals a total of 758 novel cassette-on exons and 167 novel retained introns that have no EST evidence from the same species. RT-PCR-sequencing experiments validated ~50 ~80% of the tested exons, indicating high presence of exons predicted by ENACE. ENACE is particularly powerful when applied to closely related species. In addition, our analysis shows that the ENACE-identified AS exons tend not to pass the nonsynonymous-to-synonymous substitution ratio test and not to contain protein domain, implying that such exons may be under positive selection or relaxed negative selection. These AS exons may contribute to considerable inter-species functional divergence. Our analysis further indicates that a large number of exons may have been gained or lost during mammalian evolution. Moreover, a functional analysis shows that inter-species divergence of AS events may be substantial in protein carriers and receptor proteins in mammals. These exons may be of interest to studies of AS evolution. The ENACE programs and sequences of the ENACE-identified AS events are available for download.

Conclusion

ENACE can identify potential novel cassette exons and retained introns between closely related species using a comparative approach. It can also provide information regarding lineage- or species-specificity in transcript isoforms, which are important for evolutionary and functional studies.

Non-EST based prediction of exon skipping and intron retention events using Pfam information

Most of the known alternative splice events have been detected by the comparison of expressed sequence tags (ESTs) and cDNAs. However, not all splice events are represented in EST databases since ESTs have several biases. Therefore, non-EST based approaches are needed to extend our view of a transcriptome. Here, we describe a novel method for the ab initio prediction of alternative splice events that is solely based on the annotation of Pfam domains. Furthermore, we applied this approach in a genome-wide manner to all human RefSeq transcripts and predicted a total of 321 exon skipping and intron retention events. We show that this method is very reliable as 78% (250 of 321) of our predictions are confirmed by ESTs or cDNAs. Subsequent analyses of splice events within Pfam domains revealed a significant preference of alternative exon junctions to be located at the protein surface and to avoid secondary structure elements. Thus, splice events within Pfams are probable to alter the structure and function of a domain which makes them highly interesting for detailed biological investigation. As Pfam domains are annotated in many other species, our strategy to predict exon skipping and intron retention events might be important for species with a lower number of ESTs.