The adaptive significance of unproductive alternative splicing in primates

Isoforms of Base Excision Repair Enzymes Produced by Alternative Splicing

Transcripts of many enzymes involved in base excision repair (BER) undergo extensive alternative splicing, but functions of the corresponding alternative splice variants remain largely unexplored. In this review, we cover the studies describing the common alternatively spliced isoforms and disease-associated variants of DNA glycosylases, AP-endonuclease 1, and DNA polymerase beta. We also discuss the roles of alternative splicing in the regulation of their expression, catalytic activities, and intracellular transport.

Defining the genetic and evolutionary architecture of alternative splicing in response to infection

Host and environmental factors contribute to variation in human immune responses, yet the genetic and evolutionary drivers of alternative splicing in response to infection remain largely uncharacterised. Leveraging 970 RNA-sequencing profiles of resting and stimulated monocytes from 200 individuals of African- and European-descent, we show that immune activation elicits a marked remodelling of the isoform repertoire, while increasing the levels of erroneous splicing. We identify 1,464 loci associated with variation in isoform usage (sQTLs), 9% of them being stimulation-specific, which are enriched in disease-related loci. Furthermore, we detect a longstanding increased plasticity of immune gene splicing, and show that positive selection and Neanderthal introgression have both contributed to diversify the splicing landscape of human populations. Together, these findings suggest that differential isoform usage has been an important substrate of innovation in the long-term evolution of immune responses and a more recent vehicle of population local adaptation.

Genetic ancestry might influence immunological response to infection at different regulatory levels. Here, the authors use RNA-Seq to investigate the variability of alternative splicing patterns in resting and stimulated monocytes of African- and European-descent.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

α3-Deletion Isoform of HLA-A11 Modulates Cytotoxicity of NK Cells: Correlations with HIV-1 Infection of Cells

Alternative splicing occurs frequently in many genes, especially those involved in immunity. Unfortunately, the functions of many alternatively spliced molecules from immunologically relevant genes remain unknown. Classical HLA-I molecules are expressed on almost all nucleated cells and play a pivotal role in both innate and adaptive immunity. Although splice variants of HLA-I genes have been reported, the details of their functions have not been reported. In the current study, we determined the characteristics, expression, and function of a novel splice variant of HLA-A11 named HLA-A11svE4 HLA-A11svE4 is located on the cell surface without β2-microglobulin (β2m). Additionally, HLA-A11svE4 forms homodimers as well as heterodimers with HLA-A open conformers, instead of combining with β2m. Moreover, HLA-A11svE4 inhibits the activation of NK cells to protect target cells. Compared with β2m and HLA-A11, the heterodimer of HLA-A11svE4 and HLA-A11 protected target cells from lysis by NK cells more effectively. Furthermore, HLA-AsvE4 expression was upregulated by HIV-1 in vivo and by HSV, CMV, and hepatitis B virus in vitro. In addition, our findings indicated that HLA-A11svE4 molecules were functional in activating CD8⁺ T cells through Ag presentation. Taken together, these results suggested that HLA-A11svE4 can homodimerize and form a novel heterodimeric complex with HLA-A11 open conformers. Furthermore, the data are consistent with HLA-A11svE4 playing a role in the immune escape of HIV-1.

Estimation of the minimum mRNA splicing error rate in vertebrates

The majority of protein coding genes in vertebrates contain several introns that are removed by the mRNA splicing machinery. Errors during splicing can generate aberrant transcripts and degrade the transmission of genetic information thus contributing to genomic instability and disease. However, estimating the error rate of constitutive splicing is complicated by the process of alternative splicing which can generate multiple alternative transcripts per locus and is particularly active in humans. In order to estimate the error frequency of constitutive mRNA splicing and avoid bias by alternative splicing we have characterized the frequency of splice variants at three loci, HPRT, POLB, and TRPV1 in multiple tissues of six vertebrate species. Our analysis revealed that the frequency of splice variants varied widely among loci, tissues, and species. However, the lowest observed frequency is quite constant among loci and approximately 0.1% aberrant transcripts per intron. Arguably this reflects the "irreducible" error rate of splicing, which consists primarily of the combination of replication errors by RNA polymerase II in splice consensus sequences and spliceosome errors in correctly pairing exons.

Genome-wide determination of poly(A) sites in Medicago truncatula: evolutionary conservation of alternative poly(A) site choice

Background

Alternative polyadenylation (APA) plays an important role in the post-transcriptional regulation of gene expression. Little is known about how APA sites may evolve in homologous genes in different plant species. To this end, comparative studies of APA sites in different organisms are needed. In this study, a collection of poly(A) sites in Medicago truncatula, a model system for legume plants, has been generated and compared with APA sites in Arabidopsis thaliana.

Results

The poly(A) tags from a deep-sequencing protocol were mapped to the annotated M. truncatula genome, and the identified poly(A) sites used to update the annotations of 14,203 genes. The results show that 64% of M. truncatula genes possess more than one poly(A) site, comparable to the percentages reported for Arabidopsis and rice. In addition, the poly(A) signals associated with M. truncatula genes were similar to those seen in Arabidopsis and other plants. The 3′-UTR lengths are correlated in pairs of orthologous genes between M. truncatula and Arabidopsis. Very little conservation of intronic poly(A) sites was found between Arabidopsis and M. truncatula, which suggests that such sites are likely to be species-specific in plants. In contrast, there is a greater conservation of CDS-localized poly(A) sites in these two species. A sizeable number of M. truncatula antisense poly(A) sites were found. A high percentage of the associated target genes possess Arabidopsis orthologs that are also associated with antisense sites. This is suggestive of important roles for antisense regulation of these target genes.

Conclusions

Our results reveal some distinct patterns of sense and antisense poly(A) sites in Arabidopsis and M. truncatula. In so doing, this study lends insight into general evolutionary trends of alternative polyadenylation in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-615) contains supplementary material, which is available to authorized users.

Identification of novel point mutations in splicing sites integrating whole-exome and RNA-seq data in myeloproliferative diseases

Point mutations in intronic regions near mRNA splice junctions can affect the splicing process. To identify novel splicing variants from exome sequencing data, we developed a bioinformatics splice-site prediction procedure to analyze next-generation sequencing (NGS) data (SpliceFinder). SpliceFinder integrates two functional annotation tools for NGS, ANNOVAR and MutationTaster and two canonical splice site prediction programs for single mutation analysis, SSPNN and NetGene2. By SpliceFinder, we identified somatic mutations affecting RNA splicing in a colon cancer sample, in eight atypical chronic myeloid leukemia (aCML), and eight CML patients. A novel homozygous splicing mutation was found in APC (NM_000038.4:c.1312+5G>A) and six heterozygous in GNAQ (NM_002072.2:c.735+1C>T), ABCC3 (NM_003786.3:c.1783-1G>A), KLHDC1 (NM_172193.1:c.568-2A>G), HOOK1 (NM_015888.4:c.1662-1G>A), SMAD9 (NM_001127217.2:c.1004-1C>T), and DNAH9 (NM_001372.3:c.10242+5G>A). Integrating whole-exome and RNA sequencing in aCML and CML, we assessed the phenotypic effect of mutations on mRNA splicing for GNAQ, ABCC3, HOOK1. In ABCC3 and HOOK1, RNA-Seq showed the presence of aberrant transcripts with activation of a cryptic splice site or intron retention, validated by the reverse transcription-polymerase chain reaction (RT-PCR) in the case of HOOK1. In GNAQ, RNA-Seq showed 22% of wild-type transcript and 78% of mRNA skipping exon 5, resulting in a 4–6 frameshift fusion confirmed by RT-PCR. The pipeline can be useful to identify intronic variants affecting RNA sequence by complementing conventional exome analysis.

Transcriptome analysis of human tissues and cell lines reveals one dominant transcript per gene

BACKGROUND

RNA sequencing has opened new avenues for the study of transcriptome composition. Significant evidence has accumulated showing that the human transcriptome contains in excess of a hundred thousand different transcripts. However, it is still not clear to what extent this diversity prevails when considering the relative abundances of different transcripts from the same gene.

RESULTS

Here we show that, in a given condition, most protein coding genes have one major transcript expressed at significantly higher level than others, that in human tissues the major transcripts contribute almost 85 percent to the total mRNA from protein coding loci, and that often the same major transcript is expressed in many tissues. We detect a high degree of overlap between the set of major transcripts and a recently published set of alternatively spliced transcripts that are predicted to be translated utilizing proteomic data. Thus, we hypothesize that although some minor transcripts may play a functional role, the major ones are likely to be the main contributors to the proteome. However, we still detect a non-negligible fraction of protein coding genes for which the major transcript does not code a protein.

CONCLUSIONS

Overall, our findings suggest that the transcriptome from protein coding loci is dominated by one transcript per gene and that not all the transcripts that contribute to transcriptome diversity are equally likely to contribute to protein diversity. This observation can help to prioritize candidate targets in proteomics research and to predict the functional impact of the detected changes in variation studies.

Alternative splicing in the regulation of cholesterol homeostasis

PURPOSE OF REVIEW

With the advent of whole-transcriptome sequencing, or RNA-seq, we now know that alternative splicing is a generalized phenomenon, with nearly all multiexonic genes subject to alternative splicing. In this review, we highlight recent studies examining alternative splicing as a modulator of cellular cholesterol homeostasis and as an underlying mechanism of dyslipidemia.

RECENT FINDINGS

A number of key genes involved in cholesterol metabolism are known to undergo functionally relevant alternative splicing. Recently, we have identified coordinated changes in alternative splicing in multiple genes in response to alterations in cellular sterol content. We and others have implicated several splicing factors as regulators of lipid metabolism. Furthermore, a number of cis-acting human gene variants that modulate alternative splicing have been implicated in a variety of human metabolic diseases.

SUMMARY

Alternative splicing is of importance in various types of genetically influenced dyslipidemias and in the regulation of cellular cholesterol metabolism.

Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3′UTR in AmoebaDB, providing valuable resources to the community.

Functional implications of novel human acid sphingomyelinase splice variants

BACKGROUND

Acid sphingomyelinase (ASM) hydrolyses sphingomyelin and generates the lipid messenger ceramide, which mediates a variety of stress-related cellular processes. The pathological effects of dysregulated ASM activity are evident in several human diseases and indicate an important functional role for ASM regulation. We investigated alternative splicing as a possible mechanism for regulating cellular ASM activity.

METHODOLOGY/PRINCIPAL FINDINGS

We identified three novel ASM splice variants in human cells, termed ASM-5, -6 and -7, which lack portions of the catalytic- and/or carboxy-terminal domains in comparison to full-length ASM-1. Differential expression patterns in primary blood cells indicated that ASM splicing might be subject to regulatory processes. The newly identified ASM splice variants were catalytically inactive in biochemical in vitro assays, but they decreased the relative cellular ceramide content in overexpression studies and exerted a dominant-negative effect on ASM activity in physiological cell models.

CONCLUSIONS/SIGNIFICANCE

These findings indicate that alternative splicing of ASM is of functional significance for the cellular stress response, possibly representing a mechanism for maintaining constant levels of cellular ASM enzyme activity.

The importance of identifying alternative splicing in vertebrate genome annotation

While alternative splicing (AS) can potentially expand the functional repertoire of vertebrate genomes, relatively few AS transcripts have been experimentally characterized. We describe our detailed manual annotation of vertebrate genomes, which is generating a publicly available geneset rich in AS. In order to achieve this we have adopted a highly sensitive approach to annotating gene models supported by correctly mapped, canonically spliced transcriptional evidence combined with a highly cautious approach to adding unsupported extensions to models and making decisions on their functional potential. We use information about the predicted functional potential and structural properties of every AS transcript annotated at a protein-coding or non-coding locus to place them into one of eleven subclasses. We describe the incorporation of new sequencing and proteomics technologies into our annotation pipelines, which are used to identify and validate AS. Combining all data sources has led to the production of a rich geneset containing an average of 6.3 AS transcripts for every human multi-exon protein-coding gene. The datasets produced have proved very useful in providing context to studies investigating the functional potential of genes and the effect of variation may have on gene structure and function.

Database URL:http://www.ensembl.org/index.html, http://vega.sanger.ac.uk/index.html

The origins, evolution, and functional potential of alternative splicing in vertebrates

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

Alternative splicing in self-renewal of embryonic stem cells

Much of embryonic stem cell biology has focused on transcriptional expression and regulation of genes that could mediate its unique potential in self-renewal or pluripotency. In alignment with our present understanding on the genetic, protein, and epigenetic factors that may direct cell fate, we present a short overview of the often overlooked contribution of alternative splice variants to regulatory diversity. Progressing beyond the limitations of a fixed genomic sequence, alternative splicing offers an additional layer of complexity to produce protein variants that may differ in function and localization that can direct embryonic stem cells to specific differentiation pathways. In light of the number of variants that can be produced at key ES cell genes alone, it is challenging to consider how much more multifaceted transcriptional regulation truly is, and if this can be captured more fully in future works.

Constant splice-isoform ratios in human lymphoblastoid cells support the concept of a splico-stat

Splicing generates mature transcripts from genes in pieces in eukaryotic cells. Overwhelming evidence has accumulated that alternative routes in splicing are possible for most human and mammalian genes, thereby allowing formation of different transcripts from one gene. No function has been assigned to the majority of identified alternative splice forms, and it has been assumed that they compose inert or tolerated waste from aberrant or noisy splicing. Here we demonstrate that five human transcription units (WT1, NOD2, GNAS, RABL2A, RABL2B) have constant splice-isoform ratios in genetically diverse lymphoblastoid cell lines independent of the type of alternative splicing (exon skipping, alternative donor/acceptor, tandem splice sites) and gene expression level. Even splice events that create premature stop codons and potentially trigger nonsense-mediated mRNA decay are found at constant fractions. The analyzed alternative splicing events were qualitatively but not quantitatively conserved in corresponding chimpanzee cell lines. Additionally, subtle splicing at tandem acceptor splice sites (GNAS, RABL2A/B) was highly constrained and strongly depends on the upstream donor sequence content. These results also demonstrate that unusual and unproductive splice variants are produced in a regulated manner.

Expression of insulin-like 3 (INSL3) and differential splicing of its receptor in the ovary of rhesus macaques

BACKGROUND

Although insulin-like 3 (INSL3) has been identified in the gonad of both sexes in many species, there are only limited reports on the distribution of INSL3 and its receptor, relaxin/insulin-like family peptide receptor 2 (RXFP2), in the primate ovary. Since the hormone-receptor pair is believed to play a role in female reproduction, investigating the transcription of INSL3/RXFP2 genes and the spatiotemporal expression of INSL3 in the nonhuman primate may shed light on the functional aspects of the system in humans.

METHODS

Database mining, molecular and immunological methods were applied.

RESULTS

One single INSL3 transcript and three novel splice variant transcripts of RXFP2 were identified in the ovary of rhesus macaques. While the full-length RXFP2 transcript is barely detectable in granulosa cells during the periovulatory period, INSL3 transcript and protein are highly abundant in theca cells surrounding antral follicles. Moreover, the INSL3 level in follicular fluid is 3-4 times higher than that in female serum which remains low throughout the menstrual cycle.

CONCLUSIONS

The presence of INSL3 and its receptor in the ovary implies a potential role of the ligand-receptor pair in female reproduction in nonhuman primates. However, the existence of multiple splice variants of RXFP2 indicates a very complex nature of the hormone-receptor system.