Alternative splicing mechanisms for the modulation of protein function: conservation between human and other species

Molecular cloning of a hyaluronidase from Bothrops pauloensis venom gland

Background

Hyaluronate is one of the major components of extracellular matrix from vertebrates whose breakdown is catalyzed by the enzyme hyaluronidase. These enzymes are widely described in snake venoms, in which they facilitate the spreading of the main toxins in the victim’s body during the envenoming. Snake venoms also present some variants (hyaluronidases-like substances) that are probably originated by alternative splicing, even though their relevance in envenomation is still under investigation. Hyaluronidases-like proteins have not yet been purified from any snake venom, but the cDNA that encodes these toxins was already identified in snake venom glands by transcriptomic analysis. Herein, we report the cloning and in silico analysis of the first hyaluronidase-like proteins from a Brazilian snake venom.

Methods

The cDNA sequence of hyaluronidase was cloned from the transcriptome of Bothrops pauloensis venom glands. This sequence was submitted to multiple alignment with other related sequences by ClustalW. A phylogenetic analysis was performed using MEGA 4 software by the neighbor joining (NJ) method.

Results

The cDNA from Bothrops pauloensis venom gland that corresponds to hyaluronidase comprises 1175 bp and codifies a protein containing 194 amino acid residues. The sequence, denominated BpHyase, was identified as hyaluronidase-like since it shows high sequence identities (above 83%) with other described snake venom hyaluronidase-like sequences. Hyaluronidases-like proteins are thought to be products of alternative splicing implicated in deletions of central amino acids, including the catalytic residues. Structure-based sequence alignment of BpHyase to human hyaluronidase hHyal-1 demonstrates a loss of some key secondary structures. The phylogenetic analysis indicates an independent evolution of BpHyal when compared to other hyaluronidases. However, these toxins might share a common ancestor, thus suggesting a broad hyaluronidase-like distribution among venomous snakes.

Conclusions

This work is the first report of a cDNA sequence of hyaluronidase from Brazilian snake venoms. Moreover, the in silico analysis of its deduced amino acid sequence opens new perspectives about the biological function of hyaluronidases-like proteins and may direct further studies comprising their isolation and/or recombinant production, as well as their structural and functional characterization.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.

Characterization of the impact of alternative splicing on protein dynamics: the cases of glutathione S-transferase and ectodysplasin-A isoforms

Recent studies have shown how alternative splicing (AS), the process by which eukaryotic genes express more than one product, affects protein sequence and structure. However, little information is available on the impact of AS on protein dynamics, a property fundamental for protein function. In this work, we have addressed this issue using molecular dynamics simulations of the isoforms of two model proteins: glutathione S-transferase and ectodysplasin-A. We have found that AS does not have a noticeable impact on global or local structure fluctuations. We have also found that, quite interestingly, AS has a significant effect on the coupling between key structural elements such as surface cavities. Our results provide the first atom-level view of the impact of AS on protein dynamics, as far as we know. They can contribute to refine our present view of the relationship between AS and protein disorder and, more importantly, they reveal how AS may modify structural dynamic couplings in proteins.

Identification of a novel isoform of DHRS4 protein with a nuclear localization signal

The DHRS4 gene encodes an NADP(H)-dependent retinol dehydrogenase/reductase (NRDR) and plays an important role in regulating the synthesis of retinoic acid. In the present study, we identified a novel splice RNA variant, designated NRDRA2, of the human DHRS4 gene by RT-PCR, 3' RACE, and 5' RACE. NRDRA2 mRNA lacked exons 4 and 6, and had a shift in the reading frame when compared to DHRS4 mRNA, resulting in loss of the peroxisomal targeting signal of NRDR and gain of a nuclear localization signal in the predicted NRDRA2 protein. Endogenous NRDRA2 protein was identified in the human cervical carcinoma cell line HeLa by immunoprecipitation and mass spectrometric assay. A green fluorescent protein reporter assay showed that NRDRA2 protein mainly localized to the nuclei, confirming the sequence at its C-terminus as a legitimate nuclear localization signal sequence. This study identifies the alternative transcript variant NRDRA2 encoding a subcellular nuclear localized NRDRA2 protein.

The Epithelial Sodium Channel alpha subunit (alpha ENaC) alternatively spliced form "b" in Dahl rats: What's next?

Background

The amiloride-sensitive Epithelial Sodium Channel (ENaC) is critical in maintaining Na⁺ balance, extracellular fluid volume and long term blood pressure control. ENaC is composed of three main subunits α, β, & γ. While α ENaC is critical for channel functionality, β & γ ENaC maximize channel function. To date, there are four alternatively spliced forms of the α subunit of ENaC (α ENaC-a, -b, -c, & -d) that have been published in rats, in addition to the major α ENaC transcript. While α ENaC-a, -c & -d transcripts are low abundance transcripts compared to full-length α ENaC, α ENaC-b is a higher abundance and salt-sensitive transcript compared to full-length α ENaC.

Presentation of the hypothesis

α ENaC-b protein, which is preferentially produced in Dahl R rats, to a greater extent on high salt diet, exerts a dominant negative effect on full-length α ENaC subunit by physically binding to and trapping full-length α ENaC subunit in the endoplasmic reticulum, and finally accelerating full-length α ENaC proteolytic degradation in a dose-dependent manner.

Testing the hypothesis

1) To examine the mRNA and protein abundance of α ENaC-b relative to α ENaC full-length in kidney, lung, and taste tissues of Dahl rats. 2) To compare the expression (mRNA and protein) of α ENaC-b in kidneys of Dahl S and R rats on regular and high salt diet. 3) To examine the putative binding of α ENaC-b proteins to full-length α ENaC in vitro and to determine the impact of such binding on full-length α ENaC expression in vitro.

Implications of the hypothesis

Our studies will be the first to demonstrate the over-expression of salt-sensitive α ENaC-b spliced form in kidney tissues of Dahl R rats at the expense of full-length α ENaC. The current proposal will provide highly novel insights into the putative mechanisms leading to ENaC hypoactivity in high-salt-fed Dahl R rats. Finally, findings from the present proposal will uncover a new mechanism by which alternative splicing may control the regulation of ENaC expression/function.

Regulation of the epithelial sodium channel [ENaC] in kidneys of salt-sensitive Dahl rats: insights on alternative splicing

The epithelial sodium channel [ENaC] is critical for the maintenance of sodium balance, extracellular fluid volume and long term blood pressure control. Monogenic disorders causing ENaC hyperactivity have led to a severe form of hereditary hypertension in humans, known as Liddle's syndrome. Similarly, in animal models, ENaC hyperactivity has been well documented in kidneys of salt-sensitive [S] Dahl rats [a genetic model of salt-sensitive hypertension] versus their normotensive control [Dahl salt-resistant [R] rats]. The purpose of the present review is to highlight the differential regulation of ENaC in kidneys of Dahl S versus R rats. A systematic overview of the putative role of alternative splicing of the main alpha subunit of ENaC [alpha ENaC] in modulating ENaC expression in kidneys of Dahl rats will be discussed. Finally, a better understanding of the meaningful contribution of ENaC in the pathogenesis of salt-sensitive hypertension will be achieved upon completion of this review.

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.

A Novel Mechanism in Regulating the Alpha-Subunit of the Epithelial Sodium Channel (α ENaC) by the Alternatively Spliced Form α ENaC-b

Introduction

In Dahl rats’ kidney cortex, the alternatively spliced form of the epithelial sodium channel α subunit (α ENaC-b) is the most abundant mRNA transcript (32+/-3 fold > α ENaC-wt) as was investigated by quantitative RT-PCR analysis. α ENaC-b mRNA levels were significantly higher in Dahl R versus S rats, and were further augmented by high salt diet.

Objectives

In the present study, we described the molecular cloning and searched for a possible role of α ENaC-b by testing its potential expression in COS7 cells as well as its impact on α ENaC-wt expression levels when co-expressed in COS7 cells in a dose-dependent manner.

Methods

Using RT-PCR strategy, the full-length wildtype α ENaC transcript and the alternatively spliced form α ENaC-b were amplified, sequenced, cloned, subcloned into PCMV-sport6 expression vector, expressed and co-expressed into COS7 cells in a dose-dependent manner. A combination of denaturing and native western blotting techniques was employed to examine the expression of α ENaC-b in vitro, and to determine if an interaction between α ENaC-b and α ENaC-wt occurs in vitro, and finally to demonstrate if degradation of α ENaC-wt protein does occur.

Results

α ENaC-b is translated in COS7 cells. Co-expression of α ENaC-b together with α ENaC-wt reduced α ENaC-wt levels in a dose-dependent manner. α ENaC-wt and α ENaC-b appear to form a complex that enhances the degradation of α ENaC-wt.

Conclusions

Western blots suggest a novel mechanism in α ENaC regulation whereby α ENaC-b exerts a dominant negative effect on α ENaC-wt expression. This is potentially by sequestering α ENaC-wt, enhancing its proteolytic degradation, and possibly explaining the mechanism of salt-resistance in Dahl R rats.

Comparative analysis indicates that alternative splicing in plants has a limited role in functional expansion of the proteome

Background

Alternative splicing (AS) is a widespread phenomenon in higher eukaryotes but the extent to which it leads to functional protein isoforms and to proteome expansion at large is still a matter of debate. In contrast to animal species, for which AS has been studied extensively at the protein and functional level, protein-centered studies of AS in plant species are scarce. Here we investigate the functional impact of AS in dicot and monocot plant species using a comparative approach.

Results

Detailed comparison of AS events in alternative spliced orthologs from the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice) revealed that the vast majority of AS events in both species do not result from functional conservation. Transcript isoforms that are putative targets for the nonsense-mediated decay (NMD) pathway are as likely to contain conserved AS events as isoforms that are translated into proteins. Similar results were obtained when the same comparison was performed between the two more closely related monocot species rice and Zea mays (maize).

Genome-wide computational analysis of functional protein domains encoded in alternatively and constitutively spliced genes revealed that only the RNA recognition motif (RRM) is overrepresented in alternatively spliced genes in all species analyzed. In contrast, three domain types were overrepresented in constitutively spliced genes. AS events were found to be less frequent within than outside predicted protein domains and no domain type was found to be enriched with AS introns. Analysis of AS events that result in the removal of complete protein domains revealed that only a small number of domain types is spliced-out in all species analyzed. Finally, in a substantial fraction of cases where a domain is completely removed, this domain appeared to be a unit of a tandem repeat.

Conclusion

The results from the ortholog comparisons suggest that the ability of a gene to produce more than one functional protein through AS does not persist during evolution. Cross-species comparison of the results of the protein-domain oriented analyses indicates little correspondence between the analyzed species. Based on the premise that functional genetic features are most likely to be conserved during evolution, we conclude that AS has only a limited role in functional expansion of the proteome in plants.

Functional consequences of new exon acquisition in mammalian chromodomain Y-like (CDYL) genes

The origin of new exons is an important mechanism for proteome diversity. Here, we report the recurrent origination of new exons in mammalian chromodomain Y-like (CDYL) genes and the functional consequences associated with the acquisition of the new exons. Driven by positive selection, the newly evolved longer peptide exhibits weaker transcription repression activity and attenuates the repression activity of the old form, suggesting that the acquisition of the new exons is functionally significant.

The functional modulation of epigenetic regulators by alternative splicing

BACKGROUND

Epigenetic regulators (histone acetyltransferases, methyltransferases, chromatin-remodelling enzymes, etc) play a fundamental role in the control of gene expression by modifying the local state of chromatin. However, due to their recent discovery, little is yet known about their own regulation. This paper addresses this point, focusing on alternative splicing regulation, a mechanism already known to play an important role in other protein families, e.g. transcription factors, membrane receptors, etc.

RESULTS

To this end, we compiled the data available on the presence/absence of alternative splicing for a set of 160 different epigenetic regulators, taking advantage of the relatively large amount of unexplored data on alternative splicing available in public databases. We found that 49 % (70 % in human) of these genes express more than one transcript. We then studied their alternative splicing patterns, focusing on those changes affecting the enzyme's domain composition. In general, we found that these sequence changes correspond to different mechanisms, either repressing the enzyme's function (e.g. by creating dominant-negative inhibitors of the functional isoform) or creating isoforms with new functions.

CONCLUSION

We conclude that alternative splicing of epigenetic regulators can be an important tool for the function modulation of these enzymes. Considering that the latter control the transcriptional state of large sets of genes, we propose that epigenetic regulation of gene expression is itself strongly regulated by alternative splicing.

A procedure for identifying homologous alternative splicing events

Background

The study of the functional role of alternative splice isoforms of a gene is a very active area of research in biology. The difficulty of the experimental approach (in particular, in its high-throughput version) leaves ample room for the development of bioinformatics tools that can provide a useful first picture of the problem. Among the possible approaches, one of the simplest is to follow classical protein function annotation protocols and annotate target alternative splice events with the information available from conserved events in other species. However, the application of this protocol requires a procedure capable of recognising such events. Here we present a simple but accurate method developed for this purpose.

Results

We have developed a method for identifying homologous, or equivalent, alternative splicing events, based on the combined use of neural networks and sequence searches. The procedure comprises four steps: (i) BLAST search for homologues of the two isoforms defining the target alternative splicing event; (ii) construction of all possible candidate events; (iii) scoring of the latter with a series of neural networks; and (iv) filtering of the results. When tested in a set of 473 manually annotated pairs of homologous events, our method showed a good performance, with an accuracy of 0.99, a precision of 0.98 and a sensitivity of 0.93. When no candidates were available, the specificity of our method varied between 0.81 and 0.91.

Conclusion

The method described in this article allows the identification of homologous alternative splicing events, with a good success rate, indicating that such method could be used for the development of functional annotation of alternative splice isoforms.

The (in)dependence of alternative splicing and gene duplication

Alternative splicing (AS) and gene duplication (GD) both are processes that diversify the protein repertoire. Recent examples have shown that sequence changes introduced by AS may be comparable to those introduced by GD. In addition, the two processes are inversely correlated at the genomic scale: large gene families are depleted in splice variants and vice versa. All together, these data strongly suggest that both phenomena result in interchangeability between their effects. Here, we tested the extent to which this applies with respect to various protein characteristics. The amounts of AS and GD per gene are anticorrelated even when accounting for different gene functions or degrees of sequence divergence. In contrast, the two processes appear to be independent in their influence on variation in mRNA expression. Further, we conducted a detailed comparison of the effect of sequence changes in both alternative splice variants and gene duplicates on protein structure, in particular the size, location, and types of sequence substitutions and insertions/deletions. We find that, in general, alternative splicing affects protein sequence and structure in a more drastic way than gene duplication and subsequent divergence. Our results reveal an interesting paradox between the anticorrelation of AS and GD at the genomic level, and their impact at the protein level, which shows little or no equivalence in terms of effects on protein sequence, structure, and function. We discuss possible explanations that relate to the order of appearance of AS and GD in a gene family, and to the selection pressure imposed by the environment.

Using several pair-wise informant sequences for de novo prediction of alternatively spliced transcripts

Background

As part of the ENCODE Genome Annotation Assessment Project (EGASP), we developed the MARS extension to the Twinscan algorithm. MARS is designed to find human alternatively spliced transcripts that are conserved in only one or a limited number of extant species. MARS is able to use an arbitrary number of informant sequences and predicts a number of alternative transcripts at each gene locus.

Results

MARS uses the mouse, rat, dog, opossum, chicken, and frog genome sequences as pairwise informant sources for Twinscan and combines the resulting transcript predictions into genes based on coding (CDS) region overlap. Based on the EGASP assessment, MARS is one of the more accurate dual-genome prediction programs. Compared to the GENCODE annotation, we find that predictive sensitivity increases, while specificity decreases, as more informant species are used. MARS correctly predicts alternatively spliced transcripts for 11 of the 236 multi-exon GENCODE genes that are alternatively spliced in the coding region of their transcripts. For these genes a total of 24 correct transcripts are predicted.

Conclusion

The MARS algorithm is able to predict alternatively spliced transcripts without the use of expressed sequence information, although the number of loci in which multiple predicted transcripts match multiple alternatively spliced transcripts in the GENCODE annotation is relatively small.

Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences

Recently, it was proposed that alternative splicing may act as a mechanism for opening accelerated paths of evolution, by reducing negative selection pressure, but there has been little evidence so far that this mechanism could produce adaptive benefit. Here, we use metrics of very different types of selection pressures [e.g., against amino acid mutations (Ka/Ks), against mutations at synonymous sites (Ks), and for protein reading-frame preservation] to address this question by genomewide analyses of human, chimpanzee, mouse, and rat. These data show that alternative splicing relaxes Ka/Ks selection pressure up to 7-fold, but intriguingly this effect is accompanied by a strong increase in selection pressure against synonymous mutations, which propagates into the adjacent intron, and correlates strongly with the alternative splicing level observed for each exon. These effects are highly local to the alternatively spliced exon. Comparisons of these four genomes consistently show an increase in the density of amino acid mutations (Ka) in alternatively spliced exons and a decrease in the density of synonymous mutations (Ks). This selection pressure against synonymous mutations in alternatively spliced exons was accompanied in all four genomes by a striking increase in selection pressure for protein reading-frame preservation, and both increased markedly with increasing evolutionary age. Restricting our analysis to a subset of exons with strong evidence for biologically functional alternative splicing produced identical results. Thus alternative splicing apparently can create evolutionary "hotspots" within a protein sequence, and these events have evidently been selected for during mammalian evolution.