Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Alternative splicing and genetic diversity: silencers are more frequently modified by SNVs associated with alternative exon/intron borders

With the availability of a large amount of genomic data it is expected that the influence of single nucleotide variations (SNVs) in many biological phenomena will be elucidated. Here, we approached the problem of how SNVs affect alternative splicing. First, we observed that SNVs and exonic splicing regulators (ESRs) independently show a biased distribution in alternative exons. More importantly, SNVs map more frequently in ESRs located in alternative exons than in ESRs located in constitutive exons. By looking at SNVs associated with alternative exon/intron borders (by their common presence in the same cDNA molecule), we observed that a specific type of ESR, the exonic splicing silencers (ESSs), are more frequently modified by SNVs. Our results establish a clear association between genetic diversity and alternative splicing involving ESSs.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

Alternative exon usage in the single CPT1 gene of Drosophila generates functional diversity in the kinetic properties of the enzyme:differential expression of alternatively spliced variants in drosophila tissues

The Drosophila melanogaster genome contains only one CPT1 gene (Jackson, V. N., Cameron, J. M., Zammit, V. A., and Price, N. T. (1999) Biochem. J. 341, 483-489). We have now extended our original observation to all insect genomes that have been sequenced, suggesting that a single CPT1 gene is a universal feature of insect genomes. We hypothesized that insects may be able to generate kinetically distinct variants by alternative splicing of their single CPT1 gene. Analysis of the insect genomes revealed that (a) the single CPT1 gene in each and every insect genome contains two alternative exons and (ii) in all cases, the putative alternative splicing site occurs within a small region corresponding to 21 amino acid residues that are known to be essential for the binding of substrates and of malonyl-CoA in mammalian CPT1A.Weperformed PCR analyses of mRNA from different Drosophila tissues; both of the anticipated splice variants of CPT1mRNAwere found to be expressed in all of the tissues tested (both in larvae and adults), with the expression level for one of the splice variants being significantly different between flight muscle and the fat body of adult Drosophila. Heterologous expression of the full-length cDNAs corresponding to the two putative variants of Drosophila CPT1 in the yeast Pichia pastoris revealed two important differences between the properties of the two variants: (i) their affinity (K 0.5) for one of the substrates, palmitoyl-CoA, differed by 5-fold, and (ii) the sensitivity to inhibition by malonyl-CoA at fixed, higher palmitoyl-CoA concentrations was 2-fold different and associated with different kinetics of inhibition. These data indicate that alternative splicing that specifically affects a structurally crucial region of the protein is an important mechanism through which functional diversity of CPT1 kinetics is generated from the single gene that occurs in insects. © 2010 by The American Society for Biochemistry and Molecular Biology, Inc.

Alternatively spliced products of the human kinesin light chain 1 (KNS2) gene

Conventional kinesin is a microtubule-based molecular motor involved in the transport of membranous and non-membranous cargoes. The kinesin holoenzyme exists as a heterotetramer, consisting of two heavy chain and two light chain subunits. It is thought that one function of the light chains is to interact with the cargo. Alternative splicing of kinesin light chain pre-mRNA has been observed in lower organisms, although evidence for alternative splicing of the human gene has not been reported. We have identified 19 variants of the human KNS2 gene (KLC1) that are generated by alternative splicing of downstream exons, but calculate that KNS2 has the potential to produce 285919 spliceforms. Corresponding spliceforms of the mouse KLC1 gene were also identified. The alternative exons are all located 3' of exon 12 and the novel spliceforms produce both alternative carboxy termini and alternative 3' untranslated regions. The observation of multiple light chain isoforms is consistent with their proposed role in specific cargo attachment.

A new portal for the analysis of human splicing variants

Understanding alternative splicing is crucial to elucidate the mechanisms behind several biological phenomena, including diseases. The huge amount of expressed sequences available nowadays represents an opportunity and a challenge to catalog and display alternative splicing events (ASEs). Although several groups have faced this challenge with relative success, we still lack a computational tool that uses a simple and straightforward method to retrieve, name and present ASEs. Here we present SPLOOCE, a portal for the analysis of human splicing variants. SPLOOCE uses a method based on regular expressions for retrieval of ASEs. We propose a simple syntax that is able to capture the complexity of ASEs.

optomotor-blind and the horizontal and vertical system cells of the drosophila optic lobes

The horizontal and vertical system neurons (HS and VS cells) are part of a conserved set of lobula plate giant neurons (LPGNs) in the optic lobes of the adult brain. Structure and physiology of these cells are well known, predominantly from studies in larger Dipteran flies. Our knowledge about the ontogeny of these cells is limited and stems predominantly from laser ablation studies in larvae of the house fly Musca domestica. These studies suggested that the HS and VS cells stem from a single precursor, which, at least in Musca, has not yet divided in the second larval instar. A regulatory mutation (In(1)omb[H31]) in the Drosophila gene optomotor-blind (omb) leads to the selective loss of the adult HS and VS cells. This mutation causes a transient reduction in omb expression in what appears to be the entire optic lobe anlage (OLA) late in embryogenesis. Here, I have reinitiated the laser approach with the goal of identifying the presumptive embryonic HS/VS precursor cell in Drosophila. The usefulness of the laser ablation approach which has not been applied, so far, to cells lying deep within the Drosophila embryo, was first tested on two well defined embryonic sensory structures, the olfactory antenno-maxillary complex (AMC) and the light-sensitive Bolwing´s organ (BO). In the case of the AMC, the efficiency of the ablation procedure was demonstrated with a behavioral assay. When both AMCs were ablated, the response to an attractive odour (n-butanol) was clearly reduced. Interestingly, the larvae were not completely unresponsive but had a delayed response kinetics, indicating the existence of a second odour system. BO will be a useful test system for the selectivity of laser ablation when used at higher spatial resolution. An omb-Gal4 enhancer trap line was used to visualize the embryonic OLA by GFP fluorescence. This fluorescence allowed to guide the laser beam to the relevant structure within the embryo. The success of the ablations was monitored in the adult brain via the enhancer trap insertion A122 which selectively visualizes the HS and VS cell bodies. Due to their tight clustering, individual cells could not be identified in the embryonic OLA by conventional fluorescence microscopy. Nonetheless, systematic ablation of subdomains of the OLA allowed to localize the presumptive HS/VS precursor to a small area within the OLA, encompassing around 10 cells. Future studies at higher resolution should be able to identify the precursor as (an) individual cell(s). Most known lethal omb alleles do not complement the HS/VS phenotype of the In(1)omb[H31] allele. This is the expected behaviour of null alleles. Two lethal omb alleles that had been isolated previously by non-complementation of the omb hypomorphic allele bifid, have been reported, however, to complement In(1)omb[H31]. This report was based on low resolution paraffin histology of adult heads. Four mutations from this mutagenesis were characterized here in more detail (l(1)omb[11], l(1)omb[12], l(1)omb[13], and l(1)omb[15]). Using A122 as marker for the adult HS and VS cells, I could show, that only l(1)omb[11] can partly complement the HS/VS cell phenotype of In(1)omb[H31]. In order to identify the molecular lesions in these mutants, the exons and exon/intron junctions were sequenced in PCR-amplified material from heterozygous flies. Only in two mutants could the molecular cause for loss of omb function be identified: in l(1)omb[13]), a missense mutation causes the exchange of a highly conserved residue within the DNA-binding T-domain; in l(1)omb[15]), a nonsense mutation causes a C-terminal truncation. In the other two mutants apparently regulatory regions or not yet identified alternative exons are affected. To see whether mutant OMB protein in the missense mutant l(1)omb[13] is affected in DNA binding, electrophoretic shift assays on wildtype and mutant T-domains were performed. They revealed that the mutant no longer is able to bind the consensus palindromic T-box element.

Novel Gq alpha isoform is a candidate transducer of rhodopsin signaling in a Drosophila testes-autonomous pacemaker.

DGq is the alpha subunit of the heterotrimeric GTPase (G alpha), which couples rhodopsin to phospholipase C in Drosophila vision. We have uncovered three duplicated exons in dgq by scanning the GenBank data base for unrecognized coding sequences. These alternative exons encode sites involved in GTPase activity and G beta-binding, NorpA (phospholipase C)-binding, and rhodopsin-binding. We examined the in vivo splicing of dgq in adult flies and find that, in all but the male gonads, only two isoforms are expressed. One, dgqA, is the original visual isoform and is expressed in eyes, ocelli, brain, and male gonads. The other, dgqB, has the three novel exons and is widely expressed. Remarkably, all three nonvisual B exons are highly similar (82% identity at the amino acid level) to the Gq alpha family consensus, from Caenorhabditis elegans to human, but all three visual A exons are divergent (61% identity). Intriguingly, we have found a third isoform, dgqC, which is specifically and abundantly expressed in male gonads, and shares the divergent rhodopsin-binding exon of dgqA. We suggest that DGqC is a candidate for the light-signal transducer of a testes-autonomous photosensory clock. This proposal is supported by the finding that rhodopsin 2 and arrestin 1, two photoreceptor-cell-specific genes, are also expressed in male gonads.

Hybrid 'Sinta' papaya exhibits unique ACC synthase 1 cDNA isoforms

Five ripening-related ACC synthase cDNA isoforms were cloned from 80% ripe papaya cv. 'Sinta' by reverse transcription-PCR using gene-specific primers. Clone 2 had the longest transcript and contained all common exons and three alternative exons. Clones 3 and 4 contained common exons and one alternative exon each, while clone 1, the most common transcript, contained only the common exons. Clone 5 could be due to cloning artifacts and might not be a unique cDNA fragment. Thus, there are only four isoforms of ACC synthase mRNA. Southern blot analysis indicates that all five clones came from only one gene existing as a single copy in the 'Sinta' papaya genome. Multiple sequence alignment indicates that the four isoforms arise from a single gene, possibly through alternative splicing mechanisms. All the putative alternative exons were present at the 5'-end of the gene comprising the N-terminal region of the protein. 'Sinta' ACC synthase cDNAs were of the capacs 1 type and are most closely related to a 1.4 kb capacs 1-type DNA (AJ277160) from Eksotika papaya. No capacs 2-type cDNAs were cloned from 'Sinta' by RT-PCR. This is the first report of possible alternative splicing mechanism in ripening-related ACC synthase genes in hybrid papaya, possibly to modulate or fine-tune gene expression relevant to fruit ripening.

Acute intermittent porphyria: alternative splicing of hydroxymethylbilane synthase mRNA excludes exons 3 and 12

The hydroxymethylbilane synthase (HMBS) mRNAs from 44 control individuals and 30 patients suffering from acute intermittent porphyria (AIP), were screened for length differences by reverse transcriptase polymerase chain reaction (RT-PCR) and any abnormalities were characterized by direct sequencing. Examination of the mRNAs extracted from the peripheral blood lymphocytes of the samples revealed varying degrees of alternative splicing, involving the removal of exons 3 and 12. Approximately 10-50% of the mRNA molecules were affected, despite the absence of genomic splice site mutations or any major deviance from consensus splice sequence values. The preliminary data obtained from this study suggest that this event is a normal occurrence in peripheral blood lymphocytes, and may not be associated with the molecular pathology responsible for AIP. (C) 1998 Academic Press Limited.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions

It has been previously observed that the intrinsically weak variant GC donor sites, in order to be recognized by the U2-type spliceosome, possess strong consensus sequences maximized for base pair formation with U1 and U5/U6 snRNAs. However, variability in signal strength is a fundamental mechanism for splice site selection in alternative splicing. Here we report human alternative GC-AG introns (for the first time from any species), and show that while constitutive GC-AG introns do possess strong signals at their donor sites, a large subset of alternative GC-AG introns possess weak consensus sequences at their donor sites. Surprisingly, this subset of alternative isoforms shows strong consensus at acceptor exon positions 1 and 2. The improved consensus at the acceptor exon can facilitate a strong interaction with U5 snRNA, which tethers the two exons for ligation during the second step of splicing. Further, these isoforms nearly always possess alternative acceptor sites and always possess alternative acceptor sites and exhibit particularly weak polypyrimidine tracts characteristic of AG-dependent introns. The acceptor exon nucleotides are part of the consensus required for the U2AF(35)-mediated recognition of AG in such introns. Such improved consensus at acceptor exons is not found in either normal or alternative GT-AG introns having weak donor sites or weak polypyrimidine,tracts. The changes probably reflect mechanisms that allow GC-AG alternative intron isoforms to cope with two conflicting requirements, namely an apparent need for differential splice strength to direct the choice of alternative sites and a need for improved donor signals to compensate for the central mismatch base pair (C-A) in the RNA duplex of U1 snRNA and the pre-mRNA. The other important findings include (i) one in every twenty alternative introns is a GC-AG intron, and (ii) three of every five observed GC-AG introns are alternative isoforms.

Post-transcriptional regulation of the GLI1 oncogene by the expression of alternative 5 ' untranslated regions

The oncogene GLI1 is involved in the formation of basal cell carcinoma and other tumor types as a result of the aberrant signaling of the Sonic hedgehog-Patched pathway. In this study, we have identified alternative GLI1 transcripts that differ in their 5' untranslated regions (UTRs) and are generated by exon skipping. These are denoted (alpha -UTR, beta -UTR, and gamma -UTR according to the number of noncoding exons possessed (three, two, and one, respectively). The alpha- and beta -UTR forms represent the major Gli1 transcripts expressed in mouse tissues, whereas the gamma -UTR is present at relatively low levels but is markedly induced in mouse skin treated with 12-O-tetradecanoylphorbol 13-acetate, Transcripts corresponding to the murine beta and gamma forms were identified in human tissues, but significantly, only the gamma -UTR form was present in basal cell carcinomas and in proliferating cultures of a keratinocyte cell line. Flow cytometry analysis determined that the gamma -UTR variant expresses a heterologous reporter gene 14-23-fold higher than the alpha -UTR and 5-13-fold higher than the beta -UTR in a variety of cell types. Because expression of the gamma -UTR variant correlates with proliferation, consistent with a role for GLI1 in growth promotion, up-regulation of GLI1 expression through skipping of 5' noncoding exons may be an important tumorigenic mechanism.

Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human

By spliced alignment of human DNA and transcript sequence data we constructed a data set of transcript-confirmed exons and introns from 2793 genes, 796 of which (28%) were seen to have multiple isoforms. We find that over one-third of human exons can translate in more than one frame, and that this is highly correlated with G+C content. Introns containing adenosine at donor site position +3 (A3), rather than guanosine (G3), are more common in low G+C regions, while the converse is true in high G+C regions. These two classes of introns are shown to have distinct lengths, consensus sequences and correlations among splice signals, leading to the hypothesis that A3 donor sites are associated with exon definition, and G3 donor sites with intron definition. Minor classes of introns, including GC-AG, U12-type GT-AG, weak, and putative AG-dependant introns are identified and characterized. Cassette exons are more prevalent in low G+C regions, while exon isoforms are more prevalent in high G+C regions. Cassette exon events outnumber other alternative events, while exon isoform events involve truncation twice as often as extension, and occur at acceptor sites twice as often as at donor sites. Alternative splicing is usually associated with weak splice signals, and in a majority of cases, preserves the coding frame. The reported characteristics of constitutive and alternative splice signals, and the hypotheses offered regarding alternative splicing and genome organization, have important implications for experimental research into RNA processing. The 'AltExtron' data sets are available at http://www.bit.uq.edu.au/altExtron/ and http://www.ebi.ac.uk/similar tothanaraj/altExtron/.