How prevalent is functional alternative splicing in the human genome?

Comparative analyses of ESTs and cDNAs with genomic DNA predict a high frequency of alternative splicing in human genes. However, there is an ongoing debate as to how many of these predicted splice variants are functional and how many are the result of aberrant splicing (or 'noise'). To address this question, we compared alternatively spliced cassette exons that are conserved between human and mouse with EST-predicted cassette exons that are not conserved in the mouse genome. Presumably, conserved exon-skipping events represent functional alternative splicing. We show that conserved (functional) cassette exons possess unique characteristics in size, repeat content and in their influence on the protein. By contrast, most non-conserved cassette exons do not share these characteristics. We conclude that a significant portion of cassette exons evident in EST databases is not functional, and might result from aberrant rather than regulated splicing.     

Recognition of unknown conserved alternatively spliced exons

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.     

Human and mouse mitochondrial orthologs of bacterial ClpX

We have determined the cDNA sequence and exon/intron structure of the human CLPX gene encoding a human ortholog of the E. coli ClpX chaperone and protease subunit. The CLPX gene comprises 14 exons and encodes a 633-amino acid-long precursor polypeptide. The polypeptide contains an N-terminal putative mitochondrial transit peptide, and expression of a full-length ClpX cDNA tagged at its C-terminus (Myc-His) shows that the polypeptide is transported into mitochondria. FISH analysis localized the CLPX gene to human Chromosome (Chr) 15q22.1-22.32. This localization was refined by radiation hybrid mapping placing the CLPX gene 4.6 cR distal to D15S159. Murine ClpX cDNA was sequenced, and the mouse Clpx locus was mapped to a position between 31 and 42 cM offset from the centromere on mouse Chr 9. Experimental observations indicate the presence of a pseudogene in the mouse genome and sequence variability between mouse ClpX cDNAs from different strains. Alignment of the human and mouse ClpX amino acid sequences with ClpX sequences from other organisms shows that they display the typical modular organization of domains with one AAA⁺ domain common to a large group of ATPases and several other domains conserved in ClpX orthologs linked by non-conserved sequences. Notably, a C-4 zinc finger type motif is recognized in human and mouse ClpX. This motif of so far unknown function is present only in a subset of the known ClpX sequences. Received: 5 April 2000 / Accepted: 14 June 2000     

Correlation between sequence conservation and structural thermodynamics of microRNA precursors from human,mouse, and chicken genomes

Background

Previous studies have shown that microRNA precursors (pre-miRNAs) have considerably more stable secondary structures than other native RNAs (tRNA, rRNA, and mRNA) and artificial RNA sequences. However, pre-miRNAs with ultra stable secondary structures have not been investigated. It is not known if there is a tendency in pre-miRNA sequences towards or against ultra stable structures? Furthermore, the relationship between the structural thermodynamic stability of pre-miRNA and their evolution remains unclear.

Results

We investigated the correlation between pre-miRNA sequence conservation and structural stability as measured by adjusted minimum folding free energies in pre-miRNAs isolated from human, mouse, and chicken. The analysis revealed that conserved and non-conserved pre-miRNA sequences had structures with similar average stabilities. However, the relatively ultra stable and unstable pre-miRNAs were more likely to be non-conserved than pre-miRNAs with moderate stability. Non-conserved pre-miRNAs had more G+C than A+U nucleotides, while conserved pre-miRNAs contained more A+U nucleotides. Notably, the U content of conserved pre-miRNAs was especially higher than that of non-conserved pre-miRNAs. Further investigations showed that conserved and non-conserved pre-miRNAs exhibited different structural element features, even though they had comparable levels of stability.

Conclusions

We proposed that there is a correlation between structural thermodynamic stability and sequence conservation for pre-miRNAs from human, mouse, and chicken genomes. Our analyses suggested that pre-miRNAs with relatively ultra stable or unstable structures were less favoured by natural selection than those with moderately stable structures. Comparison of nucleotide compositions between non-conserved and conserved pre-miRNAs indicated the importance of U nucleotides in the pre-miRNA evolutionary process. Several characteristic structural elements were also detected in conserved pre-miRNAs.