How did alternative splicing evolve?

Cross-kingdom patterns of alternative splicing and splice recognition

BACKGROUND

Variations in transcript splicing can reveal how eukaryotes recognize intronic splice sites. Retained introns (RIs) commonly appear when the intron definition (ID) mechanism of splice site recognition inconsistently identifies intron-exon boundaries, and cassette exons (CEs) are often caused by variable recognition of splice junctions by the exon definition (ED) mechanism. We have performed a comprehensive survey of alternative splicing across 42 eukaryotes to gain insight into how spliceosomal introns are recognized.

RESULTS

All eukaryotes we studied exhibit RIs, which appear more frequently than previously thought. CEs are also present in all kingdoms and most of the organisms in our analysis. We observe that the ratio of CEs to RIs varies substantially among kingdoms, while the ratio of competing 3' acceptor and competing 5' donor sites remains nearly constant. In addition, we find the ratio of CEs to RIs in each organism correlates with the length of its introns. In all 14 fungi we examined, as well as in most of the 9 protists, RIs far outnumber CEs. This differs from the trend seen in 13 multicellular animals, where CEs occur much more frequently than RIs. The six plants we analyzed exhibit intermediate proportions of CEs and RIs.

CONCLUSION

Our results suggest that most extant eukaryotes are capable of recognizing splice sites via both ID and ED, although ED is most common in multicellular animals and ID predominates in fungi and most protists.

Comparative analysis of transposed element insertion within human and mouse genomes reveals Alu's unique role in shaping the human transcriptome

Analysis of transposed elements in the human and mouse genomes reveals many effects on the transcriptomes, including a higher level of exonization of Alu elements than other elements.

Background

Transposed elements (TEs) have a substantial impact on mammalian evolution and are involved in numerous genetic diseases. We compared the impact of TEs on the human transcriptome and the mouse transcriptome.

Results

We compiled a dataset of all TEs in the human and mouse genomes, identifying 3,932,058 and 3,122,416 TEs, respectively. We than extracted TEs located within human and mouse genes and, surprisingly, we found that 60% of TEs in both human and mouse are located in intronic sequences, even though introns comprise only 24% of the human genome. All TE families in both human and mouse can exonize. TE families that are shared between human and mouse exhibit the same percentage of TE exonization in the two species, but the exonization level of Alu, a primate-specific retroelement, is significantly greater than that of other TEs within the human genome, leading to a higher level of TE exonization in human than in mouse (1,824 exons compared with 506 exons, respectively). We detected a primate-specific mechanism for intron gain, in which Alu insertion into an exon creates a new intron located in the 3' untranslated region (termed 'intronization'). Finally, the insertion of TEs into the first and last exons of a gene is more frequent in human than in mouse, leading to longer exons in human.

Conclusion

Our findings reveal many effects of TEs on these two transcriptomes. These effects are substantially greater in human than in mouse, which is due to the presence of Alu elements in human.

Translational control of intron splicing in eukaryotes

Most eukaryotic genes are interrupted by non-coding introns that must be accurately removed from pre-messenger RNAs to produce translatable mRNAs. Splicing is guided locally by short conserved sequences, but genes typically contain many potential splice sites, and the mechanisms specifying the correct sites remain poorly understood. In most organisms, short introns recognized by the intron definition mechanism cannot be efficiently predicted solely on the basis of sequence motifs. In multicellular eukaryotes, long introns are recognized through exon definition and most genes produce multiple mRNA variants through alternative splicing. The nonsense-mediated mRNA decay (NMD) pathway may further shape the observed sets of variants by selectively degrading those containing premature termination codons, which are frequently produced in mammals. Here we show that the tiny introns of the ciliate Paramecium tetraurelia are under strong selective pressure to cause premature termination of mRNA translation in the event of intron retention, and that the same bias is observed among the short introns of plants, fungi and animals. By knocking down the two P. tetraurelia genes encoding UPF1, a protein that is crucial in NMD, we show that the intrinsic efficiency of splicing varies widely among introns and that NMD activity can significantly reduce the fraction of unspliced mRNAs. The results suggest that, independently of alternative splicing, species with large intron numbers universally rely on NMD to compensate for suboptimal splicing efficiency and accuracy.

Alternative splicing of Alu exons--two arms are better than one

Alus, primate-specific retroelements, are the most abundant repetitive elements in the human genome. They are composed of two related but distinct monomers, left and right arms. Intronic Alu elements may acquire mutations that generate functional splice sites, a process called exonization. Most exonizations occur in right arms of antisense Alu elements, and are alternatively spliced. Here we show that without the left arm, exonization of the right arm shifts from alternative to constitutive splicing. This eliminates the evolutionary conserved isoform and may thus be selected against. We further show that insertion of the left arm downstream of a constitutively spliced non-Alu exon shifts splicing from constitutive to alternative. Although the two arms are highly similar, the left arm is characterized by weaker splicing signals and lower exonic splicing regulatory (ESR) densities. Mutations that improve these potential splice signals activate exonization and shift splicing from the right to the left arm. Collaboration between two or more putative splice signals renders the intronic left arm with a pseudo-exon function. Thus, the dimeric form of the Alu element fortuitously provides it with an evolutionary advantage, allowing enrichment of the primate transcriptome without compromising its original repertoire.

Finding exonic islands in a sea of non-coding sequence: splicing related constraints on protein composition and evolution are common in intron-rich genomes

BACKGROUND

In mammals, splice-regulatory domains impose marked trends on the relative abundance of certain amino acids near exon-intron boundaries. Is this a mammalian particularity or symptomatic of exonic splicing regulation across taxa? Are such trends more common in species that a priori have a harder time identifying exon ends, that is, those with pre-mRNA rich in intronic sequence? We address these questions surveying exon composition in a sample of phylogenetically diverse genomes.

RESULTS

Biased amino acid usage near exon-intron boundaries is common throughout the metazoa but not restricted to the metazoa. There is extensive cross-species concordance as to which amino acids are affected, and reduced/elevated abundances are well predicted by knowledge of splice enhancers. Species expected to rely on exon definition for splicing, that is, those with a higher ratio of intronic to coding sequence, more introns per gene and longer introns, exhibit more amino acid skews. Notably, this includes the intron-rich basidiomycete Cryptococcus neoformans, which, unlike intron-poor ascomycetes (Schizosaccharomyces pombe, Saccharomyces cerevisiae), exhibits compositional biases reminiscent of the metazoa. Strikingly, 5 prime ends of nematode exons deviate radically from normality: amino acids strongly preferred near boundaries are strongly avoided in other species, and vice versa. This we suggest is a measure to avoid attracting trans-splicing machinery.

CONCLUSION

Constraints on amino acid composition near exon-intron boundaries are phylogenetically widespread and characteristic of species where exon localization should be problematic. That compositional biases accord with sequence preferences of splice-regulatory proteins and are absent in ascomycetes is consistent with selection on exonic splicing regulation.

Diverse roles of hnRNP L in mammalian mRNA processing: a combined microarray and RNAi analysis

Alternative mRNA splicing patterns are determined by the combinatorial control of regulator proteins and their target RNA sequences. We have recently characterized human hnRNP L as a global regulator of alternative splicing, binding to diverse C/A-rich elements. To systematically identify hnRNP L target genes on a genome-wide level, we have combined splice-sensitive microarray analysis and an RNAi-knockdown approach. As a result, we describe 11 target genes of hnRNP L that were validated by RT-PCR and that represent several new modes of hnRNP L-dependent splicing regulation, involving both activator and repressor functions: first, intron retention; second, inclusion or skipping of cassette-type exons; third, suppression of multiple exons; and fourth, alternative poly(A) site selection. In sum, this approach revealed a surprising diversity of splicing-regulatory processes as well as poly(A) site selection in which hnRNP L is involved.

One parameter to describe the mechanism of splice sites competition

The choice of a splice site is not only related to its own intrinsic strength, but also is influenced by its flanking competitors. Splice site competition is an important mechanism for splice site prediction, especially, it is a new insight for alternative splice site prediction. In this paper, the position weight matrix scoring function is used to represent splice site strength, and the mechanism of splice site competition is described by only one parameter: scoring function subtraction. While applying on the alternative splice site prediction, based on the only one parameter, 68.22% of donor sites and 70.86% of acceptor sites are correctly classified into alternative and constitutive. The prediction abilities are approximately equal to the recent method which is based on the mechanism of splice site competition. The results reveal that the scoring function subtraction is the best parameter to describe the mechanism of splice sites competition.

Computational analysis of splicing errors and mutations in human transcripts

Background

Most retained introns found in human cDNAs generated by high-throughput sequencing projects seem to result from underspliced transcripts, and thus they capture intermediate steps of pre-mRNA splicing. On the other hand, mutations in splice sites cause exon skipping of the respective exon or activation of pre-existing cryptic sites. Both types of events reflect properties of the splicing mechanism.

Results

The retained introns were significantly shorter than constitutive ones, and skipped exons are shorter than exons with cryptic sites. Both donor and acceptor splice sites of retained introns were weaker than splice sites of constitutive introns. The authentic acceptor sites affected by mutations were significantly weaker in exons with activated cryptic sites than in skipped exons. The distance from a mutated splice site to the nearest equivalent site is significantly shorter in cases of activated cryptic sites compared to exon skipping events. The prevalence of retained introns within genes monotonically increased in the 5'-to-3' direction (more retained introns close to the 3'-end), consistent with the model of co-transcriptional splicing. The density of exonic splicing enhancers was higher, and the density of exonic splicing silencers lower in retained introns compared to constitutive ones and in exons with cryptic sites compared to skipped exons.

Conclusion

Thus the analysis of retained introns in human cDNA, exons skipped due to mutations in splice sites and exons with cryptic sites produced results consistent with the intron definition mechanism of splicing of short introns, co-transcriptional splicing, dependence of splicing efficiency on the splice site strength and the density of candidate exonic splicing enhancers and silencers. These results are consistent with other, recently published analyses.

Alternative approach to a heavy weight problem

Obesity is reaching epidemic proportions in developed countries and represents a significant risk factor for hypertension, heart disease, diabetes, and dyslipidemia. Splicing mutations constitute at least 14% of disease-causing mutations, thus implicating polymorphisms that affect splicing as likely candidates for disease susceptibility. A recent study suggested that genes associated with obesity were significantly enriched for rare nucleotide variants. Here, we examined these variants and revealed that they are located near splice junctions and tend to affect exonic splicing regulatory sequences. We also show that the majority of the exons that harbor these SNPs are constitutively spliced, yet they exhibit weak splice sites, typical to alternatively spliced exons, and are hence suboptimal for recognition by the splicing machinery and prone to become alternatively spliced. Using ex vivo assays, we tested a few representative variants and show that they indeed affect splicing by causing a shift from a constitutive to an alternative pattern, suggesting a possible link between extreme body mass index and abnormal splicing patterns.

Insights into the connection between cancer and alternative splicing

Computational and experimental evidence has revealed that cancerous cells express transcript variants that are abnormally spliced, suggesting that mRNAs are more frequently alternatively spliced in cancerous tissues than in normal ones. We show that cancerous tissues exhibit lower levels of alternative splicing than do normal tissues. Moreover, we found that the distribution of types of alternative splicing differs between cancerous and normal tissues. We further show evidence suggesting that the lower levels of alternative splicing in cancerous tissues might be a result of disruption of splicing regulatory proteins.

Discovery of novel alternatively spliced C. elegans transcripts by computational analysis of SAGE data

Background

Alternative RNA splicing allows cells to produce multiple protein isoforms from one gene. These isoforms may have specialized functions, and may be tissue- or stage-specific. Our aim was to use computational analysis of SAGE and genomic data to predict alternatively spliced transcripts expressed in C. elegans.

Results

We predicted novel alternatively spliced variants and confirmed five of eighteen candidates selected for experimental validation by RT-PCR tests and DNA sequencing.

Conclusion

We show that SAGE data can be efficiently used to discover alternative mRNA isoforms, including those with skipped exons or retained introns. Our results also imply that C. elegans may produce a larger number of alternatively spliced transcripts than initially estimated.

Biased exonization of transposed elements in duplicated genes: A lesson from the TIF-IA gene

BACKGROUND

Gene duplication and exonization of intronic transposed elements are two mechanisms that enhance genomic diversity. We examined whether there is less selection against exonization of transposed elements in duplicated genes than in single-copy genes.

RESULTS

Genome-wide analysis of exonization of transposed elements revealed a higher rate of exonization within duplicated genes relative to single-copy genes. The gene for TIF-IA, an RNA polymerase I transcription initiation factor, underwent a humanoid-specific triplication, all three copies of the gene are active transcriptionally, although only one copy retains the ability to generate the TIF-IA protein. Prior to TIF-IA triplication, an Alu element was inserted into the first intron. In one of the non-protein coding copies, this Alu is exonized. We identified a single point mutation leading to exonization in one of the gene duplicates. When this mutation was introduced into the TIF-IA coding copy, exonization was activated and the level of the protein-coding mRNA was reduced substantially. A very low level of exonization was detected in normal human cells. However, this exonization was abundant in most leukemia cell lines evaluated, although the genomic sequence is unchanged in these cancerous cells compared to normal cells.

CONCLUSION

The definition of the Alu element within the TIF-IA gene as an exon is restricted to certain types of cancers; the element is not exonized in normal human cells. These results further our understanding of the delicate interplay between gene duplication and alternative splicing and of the molecular evolutionary mechanisms leading to genetic innovations. This implies the existence of purifying selection against exonization in single copy genes, with duplicate genes free from such constrains.

Widespread evolutionary conservation of alternatively spliced exons in Caenorhabditis

Alternative splicing (AS) contributes to increased transcriptome and proteome diversity in various eukaryotic lineages. Previous studies showed low levels of conservation of alternatively spliced (cassette) exons within mammals and within dipterans. We report a strikingly different pattern in Caenorhabditis nematodes-more than 92% of cassette exons from Caenorhabditis elegans are conserved in Caenorhabditis briggsae and/or Caenorhabditis remanei. High levels of conservation extend to minor-form exons (present in a minority of transcripts) and are particularly pronounced for exons showing complex patterns of splicing. The functionality of the vast majority of cassette exons is underscored by various other features. We suggest that differences in conservation between lineages reflect differences in levels of functionality and further suggest that these differences are due to differences in intron length and the strength of consensus boundaries across lineages. Finally, we demonstrate an inverse relationship between AS and gene duplication, suggesting that the latter may be primarily responsible for the emergence of new functional transcripts in nematodes.

Cloning and embryonic expression patterns of the chicken CELF family

The CUG-BP and ETR-3-like factor (CELF) protein family has been implicated in the regulation of pre-mRNA alternative splicing, mRNA stability, and translation. Here we discuss the evolution and radiation of the CELF protein subfamilies, and report the cloning of the chicken CELF family members. In this study, we examined the embryonic expression patterns of the CELF family in the chick by in situ hybridization. We found that the tissue specificity reported for CELF proteins in the adult is established early during embryogenesis. Members of one subfamily, CUG-BP1 and ETR-3, are broadly expressed in the early embryo, while members of the second subfamily, CELF4-6, are restricted primarily to the nervous system. Expression patterns of individual CELF genes in several tissues, including the heart, liver, eye, and neural tube, exhibit distinct, yet overlapping, expression patterns. This suggests that different members of the CELF family play distinct functional roles during embryogenesis.

Large-scale comparative analysis of splicing signals and their corresponding splicing factors in eukaryotes

Introns are among the hallmarks of eukaryotic genes. Splicing of introns is directed by three main splicing signals: the 5' splice site (5'ss), the branch site (BS), and the polypyrimdine tract/3'splice site (PPT-3'ss). To study the evolution of these splicing signals, we have conducted a systematic comparative analysis of these signals in over 1.2 million introns from 22 eukaryotes. Our analyses suggest that all these signals have dramatically evolved: The PPT is weak among most fungi, intermediate in plants and protozoans, and strongest in metazoans. Within metazoans it shows a gradual strengthening from Caenorhabditis elegans to human. The 5'ss and the BS were found to be degenerate among most organisms, but highly conserved among some fungi. A maximum parsimony-based algorithm for reconstructing ancestral position-specific scoring matrices suggested that the ancestral 5'ss and BS were degenerate, as in metazoans. To shed light on the evolutionary variation in splicing signals, we have analyzed the evolutionary changes in the factors that bind these signals. Our analysis reveals coevolution of splicing signals and their corresponding splicing factors: The strength of the PPT is correlated to changes in key residues in its corresponding splicing factor U2AF2; limited correlation was found between changes in the 5'ss and U1 snRNA that binds it; but not between the BS and U2 snRNA. Thus, although the basic ability of eukaryotes to splice introns has remained conserved throughout evolution, the splicing signals and their corresponding splicing factors have considerably evolved, uniquely shaping the splicing mechanisms of different organisms.

Alternative splicing, gene duplication and connectivity in the genetic interaction network of the nematode worm Caenorhabditis elegans

We examined the relationship between gene duplication, alternative splicing, and connectedness in a predicted genetic interaction network using published data from the nematode worm Caenorhabditis elegans. Similar to previous results from mammals, genes belonging to families with only one member ("singletons") were significantly more likely to lack alternative splicing than were members of large multi-gene families. Genes belonging to multi-gene families lacking alternative splicing tended to have higher connectedness in the genetic interaction network than did genes in families that included one or more alternatively spliced members. Moreover, alternatively spliced genes were significantly more likely to interact with other alternatively spliced genes. These results support the hypothesis that certain key proteins with high degrees of network connectedness are subject to selection opposing the occurrence of alternatively spliced forms.

Regulation of transcription of the RNA splicing factor hSlu7 by Elk-1 and Sp1 affects alternative splicing

Alternative splicing plays a major role in transcriptome diversity and plasticity, but it is largely unknown how tissue-specific and embryogenesis-specific alternative splicing is regulated. The highly conserved splicing factor Slu7 is involved in 3' splice site selection and also regulates alternative splicing. We show that Slu7 has a unique spatial pattern of expression among human and mouse embryonic and adult tissues. We identified several functional Ets binding sites and GC-boxes in the human Slu7 (hSlu7) promoter region. The Ets and GC-box binding transcription factors, Elk-1 and Sp1, respectively, exerted opposite effects on hSlu7 transcription: Sp1 protein enhances and Elk-1 protein represses transcription in a dose-dependent manner. Sp1 protein bound to the hSlu7 promoter in vivo, and depletion of Sp1 by RNA interference (RNAi) repressed hSlu7 expression. Elk-1 protein bound to the hSlu7 promoter in vivo, and depletion of Elk-1 by RNAi caused an increase in the endogenous level of hSlu7 mRNA. Further, depletion of either Sp1 or Elk-1 affected alternative splicing. Our results provide indications of a complex transcription regulation mechanism that controls the spatial and temporal expression of Slu7, presumably allowing regulation of tissue-specific alternative splicing events.

Association of FLT3 polymorphisms with low BMD and risk of osteoporotic fracture in postmenopausal women

The genetic effects of FLT3 polymorphisms on BMD and fracture risk in postmenopausal women were studied. We found that FLT3+13348C>T polymorphism and haplotype 2 were significantly associated with low BMD and high risk of fracture.

INTRODUCTION

FMS-related tyrosine kinase 3 (FLT3) has been shown to play a critical role in the development of myelolymphoid progenitors and in the development of osteoclasts, but any possible genetic effect of FLT3 on bone metabolism has not been studied.

MATERIALS AND METHODS

To study a possible genetic effect of FLT3, we directly sequenced the FLT3 gene in 24 Korean individuals and identified 23 sequence variants. Seven polymorphisms were selected and genotyped in Korean postmenopausal women (n = 946).

RESULTS

We found that FLT3+13348C>T was associated with low BMD at the lumbar spine (p = 0.04) and femoral neck (p = 0.04). Haplotype analysis revealed that FLT3-ht2 (TTCTT) containing the rare allele in the +13348 position also showed significant association with low BMD in the lumbar spine (p = 0.04) and femoral neck (p = 0.05). Consistent with these results, the FLT3+13348C>T polymorphism and FLT3-ht2 were also significantly associated with high risk of fracture in the vertebrae (OR = 1.44-1.58; p = 0.03-0.04 and OR = 1.45-1.59; p = 0.02-0.03, respectively) and in any sites (OR = 1.34-1.81; p = 0.02-0.03 and OR = 1.34-1.81; p = 0.02-0.03, respectively).

CONCLUSIONS

These results suggest that FLT3 polymorphisms play a role in determination of BMD and subsequent fractures in postmenopausal women.

Global analysis of alternative splicing differences between humans and chimpanzees

Alternative splicing is a powerful mechanism affording extensive proteomic and regulatory diversity from a limited repertoire of genes. However, the extent to which alternative splicing has contributed to the evolution of primate species-specific characteristics has not been assessed previously. Using comparative genomics and quantitative microarray profiling, we performed the first global analysis of alternative splicing differences between humans and chimpanzees. Surprisingly, 6%-8% of profiled orthologous exons display pronounced splicing level differences in the corresponding tissues from the two species. Little overlap is observed between the genes associated with alternative splicing differences and the genes that display steady-state transcript level differences, indicating that these layers of regulation have evolved rapidly to affect distinct subsets of genes in humans and chimpanzees. The alternative splicing differences we detected are predicted to affect diverse functions including gene expression, signal transduction, cell death, immune defense, and susceptibility to diseases. Differences in expression at the protein level of the major splice variant of Glutathione S-transferase omega-2 (GSTO2), which functions in the protection against oxidative stress and is associated with human aging-related diseases, suggests that this enzyme is less active in human cells compared with chimpanzee cells. The results of this study thus support an important role for alternative splicing in establishing differences between humans and chimpanzees.

Transcript specificity in yeast pre-mRNA splicing revealed by mutations in core spliceosomal components

Appropriate expression of most eukaryotic genes requires the removal of introns from their pre–messenger RNAs (pre-mRNAs), a process catalyzed by the spliceosome. In higher eukaryotes a large family of auxiliary factors known as SR proteins can improve the splicing efficiency of transcripts containing suboptimal splice sites by interacting with distinct sequences present in those pre-mRNAs. The yeast Saccharomyces cerevisiae lacks functional equivalents of most of these factors; thus, it has been unclear whether the spliceosome could effectively distinguish among transcripts. To address this question, we have used a microarray-based approach to examine the effects of mutations in 18 highly conserved core components of the spliceosomal machinery. The kinetic profiles reveal clear differences in the splicing defects of particular pre-mRNA substrates. Most notably, the behaviors of ribosomal protein gene transcripts are generally distinct from other intron-containing transcripts in response to several spliceosomal mutations. However, dramatically different behaviors can be seen for some pairs of transcripts encoding ribosomal protein gene paralogs, suggesting that the spliceosome can readily distinguish between otherwise highly similar pre-mRNAs. The ability of the spliceosome to distinguish among its different substrates may therefore offer an important opportunity for yeast to regulate gene expression in a transcript-dependent fashion. Given the high level of conservation of core spliceosomal components across eukaryotes, we expect that these results will significantly impact our understanding of how regulated splicing is controlled in higher eukaryotes as well.

The spliceosome is a large RNA-protein machine responsible for removing the noncoding (intron) sequences that interrupt eukaryotic genes. Nearly everything known about the behavior of this machine has been based on the analysis of only a handful of genes, despite the fact that individual introns vary greatly in both size and sequence. Here we have utilized a microarray-based platform that allows us to simultaneously examine the behavior of all intron-containing genes in the budding yeast S. cerevisiae. By systematically examining the effects of individual mutants in the spliceosome on the splicing of all substrates, we have uncovered a surprisingly complex relationship between the spliceosome and its full complement of substrates. Contrary to the idea that the spliceosome engages in “generic” interactions with all intron-containing substrates in the cell, our results show that the identity of the transcript can differentially affect splicing efficiency when the machine is subtly perturbed. We propose that the wild-type spliceosome can also distinguish among its many substrates as external conditions warrant to function as a specific regulator of gene expression.

Many eukaryotic gene transcripts are spliced; here the authors show that components of the splicing complex can distinguish between different introns in highly homologous transcripts.

Functional and evolutionary analysis of alternatively spliced genes is consistent with an early eukaryotic origin of alternative splicing

Background

Alternative splicing has been reported in various eukaryotic groups including plants, apicomplexans, diatoms, amoebae, animals and fungi. However, whether widespread alternative splicing has evolved independently in the different eukaryotic groups or was inherited from their last common ancestor, and may therefore predate multicellularity, is still unknown. To better understand the origin and evolution of alternative splicing and its usage in diverse organisms, we studied alternative splicing in 12 eukaryotic species, comparing rates of alternative splicing across genes of different functional classes, cellular locations, intron/exon structures and evolutionary origins.

Results

For each species, we find that genes from most functional categories are alternatively spliced. Ancient genes (shared between animals, fungi and plants) show high levels of alternative splicing. Genes with products expressed in the nucleus or plasma membrane are generally more alternatively spliced while those expressed in extracellular location show less alternative splicing. We find a clear correspondence between incidence of alternative splicing and intron number per gene both within and between genomes. In general, we find several similarities in patterns of alternative splicing across these diverse eukaryotes.

Conclusion

Along with previous studies indicating intron-rich genes with weak intron boundary consensus and complex spliceosomes in ancestral organisms, our results suggest that at least a simple form of alternative splicing may already have been present in the unicellular ancestor of plants, fungi and animals. A role for alternative splicing in the evolution of multicellularity then would largely have arisen by co-opting the preexisting process.

The birth of new exons: mechanisms and evolutionary consequences

A significant amount of literature was dedicated to hypotheses concerning the origin of ancient introns and exons, but accumulating evidence indicates that new exons are also constantly being added to evolving genomes. Several mechanisms contribute to the creation of novel exons in metazoan genomes, including whole gene and single exon duplications, but perhaps the most intriguing are events of exonization, where intronic sequences become exons de novo. Exonizations of intronic sequences, particularly those originating from repetitive elements, are now widely documented in many genomes including human, mouse, dog, and fish. Such de novo appearance of exons is very frequently associated with alternative splicing, with the new exon-containing variant typically being the rare one. This allows the new variant to be evolutionarily tested without compromising the original one, and provides an evolutionary strategy for generation of novel functions with minimum damage to the existing functional repertoire. This review discusses the molecular mechanisms leading to exonization, its extent in vertebrate genomes, and its evolutionary implications.

Evolution of Nova-dependent splicing regulation in the brain

A large number of alternative exons are spliced with tissue-specific patterns, but little is known about how such patterns have evolved. Here, we study the conservation of the neuron-specific splicing factors Nova1 and Nova2 and of the alternatively spliced exons they regulate in mouse brain. Whereas Nova RNA binding domains are 94% identical across vertebrate species, Nova-dependent splicing silencer and enhancer elements (YCAY clusters) show much greater divergence, as less than 50% of mouse YCAY clusters are conserved at orthologous positions in the zebrafish genome. To study the relation between the evolution of tissue-specific splicing and YCAY clusters, we compared the brain-specific splicing of Nova-regulated exons in zebrafish, chicken, and mouse. The presence of YCAY clusters in lower vertebrates invariably predicted conservation of brain-specific splicing across species, whereas their absence in lower vertebrates correlated with a loss of alternative splicing. We hypothesize that evolution of Nova-regulated splicing in higher vertebrates proceeds mainly through changes in cis-acting elements, that tissue-specific splicing might in some cases evolve in a single step corresponding to evolution of a YCAY cluster, and that the conservation level of YCAY clusters relates to the functions encoded by the regulated RNAs.

Evidence for a trade-off between translational efficiency and splicing regulation in determining synonymous codon usage in Drosophila melanogaster

In Drosophila melanogaster, synonymous codons corresponding to the most abundant cognate tRNAs are used more frequently, especially in highly expressed genes. Increased use of such "optimal" codons is considered an adaptation for translational efficiency. Need it always be the case that selection should favor the use of a translationally optimal codon? Here, we investigate one possible confounding factor, namely, the need to specify information in exons necessary to enable correct splicing. As expected from such a model, in Drosophila many codons show different usage near intron-exon boundaries versus exon core regions. However, this finding is in principle also consistent with Hill-Robertson effects modulating usage of translationally optimal codons. However, several results support the splice model over the translational selection model: 1) the trends in codon usage are strikingly similar to those in mammals in which codon usage near boundaries correlates with abundance in exonic splice enhancers (ESEs), 2) codons preferred near boundaries tend to be enriched for A and avoid C (conversely those avoided near boundaries prefer C rather than A), as expected were ESEs involved, and 3) codons preferred near boundaries are typically not translationally optimal. We conclude that usage of translationally optimal codons usage is compromised in the vicinity of splice junctions in intron-containing genes, to the effect that we observe higher levels of usage of translationally optimal codons at the center of exons. On the gene level, however, controlling for known correlates of codon bias, the impact on codon usage patterns is quantitatively small. These results have implications for inferring aspects of the mechanism of splicing given nothing more than a well-annotated genome.

Gene discovery and transcript analyses in the corn smut pathogen Ustilago maydis: expressed sequence tag and genome sequence comparison

Background

Ustilago maydis is the basidiomycete fungus responsible for common smut of corn and is a model organism for the study of fungal phytopathogenesis. To aid in the annotation of the genome sequence of this organism, several expressed sequence tag (EST) libraries were generated from a variety of U. maydis cell types. In addition to utility in the context of gene identification and structure annotation, the ESTs were analyzed to identify differentially abundant transcripts and to detect evidence of alternative splicing and anti-sense transcription.

Results

Four cDNA libraries were constructed using RNA isolated from U. maydis diploid teliospores (U. maydis strains 518 × 521) and haploid cells of strain 521 grown under nutrient rich, carbon starved, and nitrogen starved conditions. Using the genome sequence as a scaffold, the 15,901 ESTs were assembled into 6,101 contiguous expressed sequences (contigs); among these, 5,482 corresponded to predicted genes in the MUMDB (MIPS Ustilago maydis database), while 619 aligned to regions of the genome not yet designated as genes in MUMDB. A comparison of EST abundance identified numerous genes that may be regulated in a cell type or starvation-specific manner. The transcriptional response to nitrogen starvation was assessed using RT-qPCR. The results of this suggest that there may be cross-talk between the nitrogen and carbon signalling pathways in U. maydis. Bioinformatic analysis identified numerous examples of alternative splicing and anti-sense transcription. While intron retention was the predominant form of alternative splicing in U. maydis, other varieties were also evident (e.g. exon skipping). Selected instances of both alternative splicing and anti-sense transcription were independently confirmed using RT-PCR.

Conclusion

Through this work: 1) substantial sequence information has been provided for U. maydis genome annotation; 2) new genes were identified through the discovery of 619 contigs that had previously escaped annotation; 3) evidence is provided that suggests the regulation of nitrogen metabolism in U. maydis differs from that of other model fungi, and 4) Alternative splicing and anti-sense transcription were identified in U. maydis and, amid similar observations in other basidiomycetes, this suggests these phenomena may be widespread in this group of fungi. These advances emphasize the importance of EST analysis in genome annotation.

U7 snRNA-mediated correction of aberrant splicing caused by activation of cryptic splice sites

A considerable fraction of mutations associated with hereditary disorders and cancers affect splicing. Some of them cause exon skipping or the inclusion of an additional exon, whereas others lead to the inclusion of intronic sequences or deletion of exonic sequences through the activation of cryptic splice sites. We focused on the latter cases and have designed a series of vectors that express modified U7 small nuclear RNAs (snRNAs) containing a sequence antisense to the cryptic splice site. Three cases of such mutation were investigated in this study. In two of them, which occurred in the PTCH1 and BRCA1 genes, canonical splice donor sites had been partially impaired by mutations that activated nearby intronic cryptic splice donor sites. Another mutation found in exonic region in CYP11A created a novel splice donor site. Transient expression of the engineered U7 snRNAs in HeLa cells restored correct splicing in a sequence-specific and dose-dependent manner in the former two cases. In contrast, the third case, in which the cryptic splice donor site in the exonic sequence was activated, the expression of modified U7 snRNA resulted in exon skipping. The correction of aberrant splicing by suppressing intronic cryptic splice sites with modified U7 is expected be a promising alternative to gene replacement therapy.

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.

Molecular and phylogenetic characterization of a novel putative membrane transporter (SLC10A7), conserved in vertebrates and bacteria

The 'Solute Carrier Family SLC10' consists of six annotated members in humans, comprising two bile acid carriers (SLC10A1 and SLC10A2), one steroid sulfate transporter (SLC10A6), and three orphan carriers (SLC10A3 to SLC10A5). In this study we report molecular characterization and expression analysis of a novel member of the SLC10 family, SLC10A7, previously known as C4orf13. SLC10A7 proteins consist of 340-343 amino acids in humans, mice, rats, and frogs and show an overall amino acid sequence identity of >85%. SLC10A7 genes comprise 12 coding exons and show broad tissue expression pattern. When expressed in Xenopus laevis oocytes and HEK293 cells, SLC10A7 was detected in the plasma membrane but revealed no transport activity for bile acids and steroid sulfates. By immunofluorescence analysis of dual hemagglutinin (HA)- and FLAG-labeled SLC10A7 proteins in HEK293 cells, we established a topology of 10 transmembrane domains with an intracellular cis orientation of the N-terminal and C-terminal ends. This topology pattern is clearly different from the seven-transmembrane domain topology of the other SLC10 members but similar to hitherto uncharacterized non-vertebrate SLC10A7-related proteins. In contrast to the established SLC10 members, which are restricted to the taxonomic branch of vertebrates, SLC10A7-related proteins exist also in yeasts, plants, and bacteria, making SLC10A7 taxonomically the most widespread member of this carrier family. Vertebrate and bacterial SLC10A7 proteins exhibit >20% sequence identity, which is higher than the sequence identity of SLC10A7 to any other member of the SLC10 carrier family.

Identification of two GTP-independent alternatively spliced forms of tissue transglutaminase in human leukocytes, vascular smooth muscle, and endothelial cells

Tissue transglutaminase (tTG) is a multifunctional enzyme with transglutaminase crosslinking (TGase), GTP binding, and hydrolysis activities that play a role in many different disorders. We identified, characterized, and investigated the function and stability of two alternatively spliced forms of tTG using biochemical, cellular, and molecular biological approaches. Using a human aortic vascular smooth muscle cells (VSMC) cDNA library, we identified two cDNAs encoding C-terminal truncated forms, tTG(V1) and tTG(V2). tTG(V1,2) mRNAs were synthesized by a rare splicing event using alternate splice sites within exons 12 and 13 of the tTG gene, respectively. Quantitative PCR and immunoblotting demonstrated that there was unique expression and localization of tTG(V1,2) compared with tTG in human umbilical vein endothelial cells (HUVECs), VSMC, and leukocytes. The loss of C-terminal 52 amino acid residues (AAs) in tTG(V1,2) altered GTP binding, enhanced GTP hydrolysis, rendered the variants insensitive to GTP inhibition, and resulted in <10% residual Ca(+2)-dependent TGase activity. Transfection in HEK293 demonstrated a 28- and 5-fold reduction in the expression of tTG(V1) and tTG(V2), respectively, demonstrating that the C-terminal GTP-binding domain is important in stabilizing and promoting the half-life of tTG. The altered affinity for GTP allowed tTG(V1,2) to exhibit enhanced TGase activity when there is a transient increase in Ca(+2) levels. The abundance of tTG(V1,2) and its distinct intracellular expression patterns in human vascular cells and leukocytes indicate these isoforms likely have unique physiological functions.

The emergence of alternative 3' and 5' splice site exons from constitutive exons

Alternative 3' and 5' splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3'ss and 5'ss exons. The results revealed that alternative 3'ss and 5'ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3'ss and 5'ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

Progress in understanding the biology of the human mutagen LINE-1

Long interspersed nucleotide element (LINE)-1 retrotransposon (L1) has emerged as the largest contributor to mammalian genome mass, responsible for over 35% of the human genome. Differences in the number and activity levels of L1s contribute to interindividual variation in humans, both by affecting an individual's likelihood of acquiring new L1-mediated mutations, as well as by differentially modifying gene expression. Here, we report on recent progress in understanding L1 biology, with a focus on mechanisms of L1-mediated disease. We discuss known details of L1 life cycle, including L1 structure, transcriptional regulation, and the mechanisms of translation and retrotransposition. Current views on cell type specificity, timing, and control of retrotransposition are put forth. Finally, we discuss the role of L1 as a mutagen, using the latest findings in L1 biology to illuminate molecular mechanisms of L1-mediated gene disruption.

Genomic checkpoints for exon 10 usage in the luteinizing hormone receptor type 1 and type 2

Alternative splicing is a hallmark of glycoprotein hormone receptor gene regulation, but its molecular mechanism is unknown. The LH receptor (LHR) gene possesses 11 exons, but exon 10 is constitutively skipped in the New World monkey lineage (LHR type 2), whereas it is constitutively spliced in the human (LHR type 1). This study identifies the regulatory elements of exon 10 usage. Sequencing of genomic marmoset DNA revealed that the cryptic LHR exon 10 is highly homologous to exon 10 from other species and displays intact splice sites. Functional studies using a minigene approach excluded the contribution of intronic, marmoset-specific long interspersed nucleotide-1 elements to exon 10 skipping. Sequencing of the genomic regions surrounding exon 10 from several primate lineages, sequence comparisons including the human and mouse LHR gene, revealed the presence of unique nucleotides at 3'-intronic position -19 and -10 and at position +26 within exon 10 of the marmoset LHR. Exon trap experiments and in vitro mutagenesis of these nucleotides resulted in the identification of a composite regulatory element of splicing consisting of cis-acting elements represented by two polypyrimidine tracts and a trans-acting element within exon 10, which affect the secondary RNA structure. Changes within this complex resulted either in constitutive exon inclusion, constitutive skipping, or alternative splicing of exon 10. This work delineates the molecular pathway leading to intronization of exon 10 in the LHR type 2 and reveals, for the first time, the essential function of regulatory and structural elements involved in glycoprotein hormone receptor splicing.

Association of KIT gene polymorphisms with bone mineral density in postmenopausal Korean women

Bone mineral density (BMD) is a major factor for determining bone strength and osteoporotic fracture risk, and is determined by environmental and multiple genetic factors. KIT, which encodes a transmembrane receptor with tyrosine kinase activity, plays an important role in the differentiation of osteoclasts. We examined the associations between KIT gene polymorphisms and BMD in postmenopausal Korean women. All exons, their boundaries, and the promoter region (approximately 1.5 kb) from 24 individuals were directly sequenced. Eighteen polymorphisms were identified, and three single-nucleotide polymorphisms (SNPs) were genotyped in all study participants (n=946). BMD at the lumbar spine and femoral neck was measured using dual-energy X-ray absorptiometry. The mean age of the study subjects was 58.9+/-7.5 years, and the mean number of years since menopause was 9.6+/-7.9 years. None of the three SNPs (-1694G>T, +41894A>G, and +49512G>A) was significantly associated with BMD value. However, multivariate analysis showed that the ht3 (-1694T-+41894A-+49512G) was significantly associated with lower BMD at the femoral neck (P=0.007 in the recessive model). These findings indicate that KIT-ht3 may be a useful genetic marker for osteoporosis and that KIT may have a role on bone metabolism in humans.

Coevolution of genomic intron number and splice sites

Spliceosomal intron numbers and boundary sequences vary dramatically in eukaryotes. We found a striking correspondence between low intron number and strong sequence conservation of 5' splice sites (5'ss) across eukaryotic genomes. The phylogenetic pattern suggests that ancestral 5'ss were relatively weakly conserved, but that some lineages independently underwent both major intron loss and 5'ss strengthening. It seems that eukaryotic ancestors had relatively large intron numbers and 'weak' 5'ss, a pattern associated with frequent alternative splicing in modern organisms.

Characterization of a novel alternatively spliced human transcript encoding an N-terminally truncated Vps24 protein that suppresses the effects of Bax in an ESCRT independent manner in yeast

Elucidating novel anti-apoptotic regulatory pathways is central to further understanding the molecular basis of several pathologies, including cancer. We have previously reported the identification of several mammalian cDNAs effective in preventing the lethal effects of heterologous expression of a pro-apoptotic BAX cDNA in yeast [Yang, Z., Khoury, C., Jean-Baptiste, G., Greenwood, M.T., 2006. Identification of mouse sphingomyelin synthase 1 (SMS1) as a suppressor of Bax mediated cell death in yeast. FEMS Yeast Res. 6, 751-762]. Here we report that one of the Bax suppressors encodes a novel 156 amino acid variant of the human Vps24 protein, Vps24beta, that lacks the N-terminal lipid binding domain of the well characterized 222 residue Vps24 (Vps24alpha). We demonstrate that the VPS24beta cDNA represents an expressed transcript that is likely produced by alternative splicing of the human VPS24 gene. Vps24alpha, but not Vps24beta, prevented the temperature and salt sensitive growth defects observed in a yeast mutant lacking a functional VPS24 gene. In contrast, Vps24beta, but not Vps24alpha, suppressed the inhibitory effects of Bax on yeast growth. Vps24beta protein also suppressed the effects of Bax in mutants lacking other VPS genes suggesting that a functional ESCRT pathway, of which the yeast Vps24p is an essential component, is not required for Vps24beta function. Taken together, we demonstrate that the human VPS24 gene gives rise to two functionally distinct proteins, one of which is involved in the ESCRT pathway and another novel protein that serves an anti-apoptotic role.

Functional anthology of intrinsic disorder. 2. Cellular components, domains, technical terms, developmental processes, and coding sequence diversities correlated with long disordered regions

Biologically active proteins without stable ordered structure (i.e., intrinsically disordered proteins) are attracting increased attention. Functional repertoires of ordered and disordered proteins are very different, and the ability to differentiate whether a given function is associated with intrinsic disorder or with a well-folded protein is crucial for modern protein science. However, there is a large gap between the number of proteins experimentally confirmed to be disordered and their actual number in nature. As a result, studies of functional properties of confirmed disordered proteins, while helpful in revealing the functional diversity of protein disorder, provide only a limited view. To overcome this problem, a bioinformatics approach for comprehensive study of functional roles of protein disorder was proposed in the first paper of this series (Xie, H.; Vucetic, S.; Iakoucheva, L. M.; Oldfield, C. J.; Dunker, A. K.; Obradovic, Z.; Uversky, V. N. Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long disordered regions. J. Proteome Res. 2007, 5, 1882-1898). Applying this novel approach to Swiss-Prot sequences and functional keywords, we found over 238 and 302 keywords to be strongly positively or negatively correlated, respectively, with long intrinsically disordered regions. This paper describes approximately 90 Swiss-Prot keywords attributed to the cellular components, domains, technical terms, developmental processes, and coding sequence diversities possessing strong positive and negative correlation with long disordered regions.

SpliceMiner: a high-throughput database implementation of the NCBI Evidence Viewer for microarray splice variant analysis

Background

There are many fewer genes in the human genome than there are expressed transcripts. Alternative splicing is the reason. Alternatively spliced transcripts are often specific to tissue type, developmental stage, environmental condition, or disease state. Accurate analysis of microarray expression data and design of new arrays for alternative splicing require assessment of probes at the sequence and exon levels.

Description

SpliceMiner is a web interface for querying Evidence Viewer Database (EVDB). EVDB is a comprehensive, non-redundant compendium of splice variant data for human genes. We constructed EVDB as a queryable implementation of the NCBI Evidence Viewer (EV). EVDB is based on data obtained from NCBI Entrez Gene and EV. The automated EVDB build process uses only complete coding sequences, which may or may not include partial or complete 5' and 3' UTRs, and filters redundant splice variants. Unlike EV, which supports only one-at-a-time queries, SpliceMiner supports high-throughput batch queries and provides results in an easily parsable format. SpliceMiner maps probes to splice variants, effectively delineating the variants identified by a probe.

Conclusion

EVDB can be queried by gene symbol, genomic coordinates, or probe sequence via a user-friendly web-based tool we call SpliceMiner (). The EVDB/SpliceMiner combination provides an interface with human splice variant information and, going beyond the very valuable NCBI Evidence Viewer, supports fluent, high-throughput analysis. Integration of EVDB information into microarray analysis and design pipelines has the potential to improve the analysis and bioinformatic interpretation of gene expression data, for both batch and interactive processing. For example, whenever a gene expression value is recognized as important or appears anomalous in a microarray experiment, the interactive mode of SpliceMiner can be used quickly and easily to check for possible splice variant issues.

Identification of novel RANK polymorphisms and their putative association with low BMD among postmenopausal women

INTRODUCTION

Bone mineral density (BMD) is the major factor for determining bone strength, which is closely correlated to osteoporotic fracture risk and is largely determined by multiple genetic factors. The RANK (TNFRSF11A), receptor for RANKL, is a member of the tumor necrosis factor receptor (TNFR) superfamily and plays a central role in osteoclast development.

METHODS

In order to investigate the effects of RANK polymorphism on BMD and osteoporosis, we directly sequenced the RANK gene in 24 Korean individuals and identified 25 sequence variants. Eleven of these polymorphisms were selected and genotyped in a larger-scale study of postmenopausal women (n = 560). Areal BMD (g/cm(2)) of the anterior-posterior lumbar spine and the nondominant proximal femur were measured using dual-energy X-ray absorptiometry.

RESULTS

We found that two intronic polymorphisms in the RANK gene [RANK + 34863G > A (rs12458117) and RANK + 35928insdelC (new polymorphism found in this study) in intron 6] were significantly associated with the BMD of the lumbar spine, i.e., rare alleles were significantly associated with low BMD of the lumbar spine among Korean postmenopausal women (p = 0.04 and 0.02, respectively). These polymorphisms were also associated with low BMD of proximal femur sites, including Ward's triangle, trochanter, and total femur. Our results suggest that +34863G > A and +35928insdelC polymorphisms in RANK are possible genetic factors for low BMD in postmenopausal women.

mRNA maturation by two-step trans-splicing/polyadenylation processing in trypanosomes

Trypanosomes are unique eukaryotic cells, in that they virtually lack mechanisms to control gene expression at the transcriptional level. These microorganisms mostly control protein synthesis by posttranscriptional regulation processes, like mRNA stabilization and degradation. Transcription in these cells is polycistronic. Tens to hundreds of protein-coding genes of unrelated function are arrayed in long clusters on the same DNA strand. Polycistrons are cotranscriptionally processed by trans-splicing at the 5' end and polyadenylation at the 3' end, generating monocistronic units ready for degradation or translation. In this work, we show that some trans-splicing/polyadenylation sites may be skipped during normal polycistronic processing. As a consequence, dicistronic units or monocistronic transcripts having long 3' UTRs are produced. Interestingly, these unspliced transcripts can be processed into mature mRNAs by the conventional trans-splicing/polyadenylation events leading to translation. To our knowledge, this is a previously undescribed mRNA maturation by trans-splicing uncoupled from transcription. We identified an RNA-recognition motif-type protein, homologous to the mammalian polypyrimidine tract-binding protein, interacting with one of the partially processed RNAs analyzed here that might be involved in exon skipping. We propose that splice-site skipping might be part of a posttranscriptional mechanism to regulate gene expression in trypanosomes, through the generation of premature nontranslatable RNA molecules.

High qualitative and quantitative conservation of alternative splicing in Caenorhabditis elegans and Caenorhabditis briggsae

Alternative splicing (AS) is an important contributor to proteome diversity and is regarded as an explanatory factor for the relatively low number of human genes compared with less complex animals. To assess the evolutionary conservation of AS and its developmental regulation, we have investigated the qualitative and quantitative expression of 21 orthologous alternative splice events through the development of 2 nematode species separated by 85-110 Myr of evolutionary time. We demonstrate that most of these alternative splice events present in Caenorhabditis elegans are conserved in Caenorhabditis briggsae. Moreover, we find that relative isoform expression levels vary significantly during development for 78% of the AS events and that this quantitative variation is highly conserved between the 2 species. Our results suggest that AS is generally tightly regulated through development and that the regulatory mechanisms controlling AS are to a large extent conserved during the evolution of Caenorhabditis. This strong conservation indicates that both major and minor splice forms have important functional roles and that the relative quantities in which they are expressed are crucial. Our results therefore suggest that the quantitative regulation of isoform expression levels is an intrinsic part of most AS events. Moreover, our results indicate that AS contributes little to transcript variation in Caenorhabditis genes and that gene duplication may be the major evolutionary mechanism for the origin of novel transcripts in these 2 species.

High-density yeast-tiling array reveals previously undiscovered introns and extensive regulation of meiotic splicing

Knowing gene structure is vital to understanding gene function, and accurate genome annotation is essential for understanding cellular function. To this end, we have developed a genome-wide assay for mapping introns in Saccharomyces cerevisiae. Using high-density tiling arrays, we compared wild-type yeast to a mutant deficient for intron degradation. Our method identified 76% of the known introns, confirmed 18 previously predicted introns, and revealed 9 formerly undiscovered introns. Furthermore, we discovered that all 13 meiosis-specific intronic yeast genes undergo regulated splicing, which provides posttranscriptional regulation of the genes involved in yeast cell differentiation. Moreover, we found that approximately 16% of intronic genes in yeast are incompletely spliced during exponential growth in rich medium, which suggests that meiosis is not the only biological process regulated by splicing. Our tiling-array assay provides a snapshot of the spliced transcriptome in yeast. This robust methodology can be used to explore environmentally distinct splicing responses and should be readily adaptable to the study of other organisms, including humans.

Cloning, expression, and characterization of an alternatively spliced variant of human heparanase

Heparanase is an endoglycosidase that cleaves heparan sulfate in the extracellular matrix (ECM) and hence participates in ECM degradation and remodeling. Heparanase is involved in fundamental biological processes such as cancer metastasis, angiogenesis, and inflammation. Alternative splicing in the coding region of human heparanase was not reported. Here, we report the cloning of a splice variant of human heparanase that lacks exon 5 and is missing 174 bp compared to the wild-type cDNA. Splice 5 is expressed as a 55 kDa protein compared to the 65 and 50 kDa latent and active wild-type enzyme. Splice 5 was not detected in the incubation medium of tumor cells as opposed to the wild-type latent heparanase. Splice 5 escaped proteolytic cleavage, was devoid of HS degradation activity and exhibited diffused rather than granular cellular localization.

Alternative splicing of pre-messenger RNAs in plants in the genomic era

Primary transcripts (precursor-mRNAs) with introns can undergo alternative splicing to produce multiple transcripts from a single gene by differential use of splice sites, thereby increasing the transcriptome and proteome complexity within and between cells and tissues. Alternative splicing in plants is largely an unexplored area of gene expression, as this phenomenon used to be considered rare. However, recent genome-wide computational analyses have revealed that alternative splicing in flowering plants is far more prevalent than previously thought. Interestingly, pre-mRNAs of many spliceosomal proteins, especially serine/arginine-rich (SR) proteins, are extensively alternatively spliced. Furthermore, stresses have a dramatic effect on alternative splicing of pre-mRNAs including those that encode many spliceosomal proteins. Although the mechanisms that regulate alternative splicing in plants are largely unknown, several reports strongly suggest a key role for SR proteins in spliceosome assembly and regulated splicing. Recent studies suggest that alternative splicing in plants is an important posttranscriptional regulatory mechanism in modulating gene expression and eventually plant form and function.

Detecting over-represented motifs in alternatively spliced exons using Gibbs sampling

Alternative pre-mRNA splicing is a biological mechanism with significant prevalence in complex organisms and experimentally verified association with numerous disease-causing factors. Splicing-related proteins play a significant regulatory role during this process. In this study, we applied a stochastic analysis of alternatively spliced human genes based on Gibbs sampling, in order to identify short consensus sequences that are over-represented compared to a reference Markov model describing constitutive exons of the same genes. The analysis resulted in a set of statistically significant over-represented motifs. The biological importance of these motifs was assessed by estimating the likelihood of being identified by cis-acting elements that correspond to the binding domains of splicing enhancers/silencers. The results indicate that the identified over-represented sequences are often similar to those recognized by known regulatory splicing elements.

Identification and Partial Characterization of a Novel CYP2C9 Splicing Variant Encoding a Protein Lacking Eight Amino Acid Residues

CYP2C9 is known as an enzyme responsible for the metabolism of various clinically important drugs. Recently, we cloned a cDNA corresponding to a CYP2C9 splicing variant (SV), which seemed to have an open reading frame of a protein with 482 amino acid residues. To investigate whether or not the SV can be translated as a functionally active protein, we expressed the CYP2C9SV in insect cells, and spectrophotometric and enzymatic properties were characterized. The CYP2C9SV protein showed a typical reduced CO-difference spectrum, indicating that the translated protein binds a heme moiety. However, CYP2C9SV did not metabolize tolbutamide or diclofenac at all, suggesting that the SV protein appeared to lack the ability to catalyze reactions mediated by CYP2C9. Although the CYP2C9SV mRNA was detected in all human liver samples examined in this study by real-time PCR, the level was generally low, ranging between 0.7 and 9.6% of the normal CYP2C9 mRNA. These results suggest that the CYP2C9SV protein is unlikely to contribute to CYP2C9 activities, although it appears to be expressed in most individuals.

RNA-based analysis of BRCA1 and BRCA2 gene alterations

Alterations in BRCA1 and BRCA2 genes account for a large proportion of hereditary breast and ovarian cancers. Mutations and variants of unknown pathological significance have been identified in both genes; however, most of them have been studied only at the genomic level, and their effect on mRNA expression remains unknown. We identified two BRCA1 and six BRCA2 splice site variants, and one BRCA2 alteration at exon 14. Our aim was to ascertain the effect on RNA processing of the variants still unclassified. We found that BRCA1 c.IVS11 + 1G>A, BRCA2 c.7252_7272delinsTG, BRCA2 c.IVS2 + 1G>A, BRCA2 c.IVS13-2A>G, BRCA2 c.IVS21 + 4A>G, and BRCA2 c.9345G>A lead to aberrant transcripts in lymphocytes. Five of these six splice site variants caused a complete inactivation of the mutant allele because they produced frameshift similar to previously described deleterious exonic variants. Therefore, we consider them to be true deleterious mutations, possibly associated with an increased lifetime risk of breast or ovarian cancer. BRCA1 c.IVS17 + 6C>G, BRCA2 c.IVS12-9del4, and BRCA2 IVS1-9del3 represent rare variants, not disrupting normal mRNA processing. The last two BRCA2 genetic variants had not been reported in the Breast Cancer Information Core BIC database.