Identification of a novel box C/D snoRNA from mouse nucleolar cDNA library

By construction and screen of mouse nucleolar cDNA library, a novel mammalian small nucleolar RNAs (snoRNA) was identified. The novel snoRNA, 70 nt in length, displays structural features typical of C/D box snoRNA family. The snoRNA possesses an 11-nt-long rRNA antisense element and is predicted to guide the 2'-O-methylation of mouse 28S rRNA at G4043, a site unknown so far to be modified in vertebrates. The comparison of functional element of snoRNA guides among eukaryotes reveals that the novel snoRNA is a mammalian counterpart of yeast snR38 despite highly divergent sequence between them. Mouse and human snR38 and other cognates in distant vertebrates were positively detected with slight length variability. As expected, the rRNA ribose-methylation site predicted by mouse snR38 was precisely mapped by specific-primer extension assay. Furthermore, our analyses show that mouse and human snR38 gene have multiple variants and are nested in the introns of different host genes with unknown function. Thus, snR38 is a phylogenetically conserved methylation guide but exhibits different genomic organization in eukaryotes.     

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.     

Evolutionary profiling of the U49 snoRNA gene

Small nucleolar RNAs (snoRNAs) are involved in precursor ribosomal RNA (pre-rRNA) processing and rRNA base modifications (2'-O-ribose methylation and pseudouridylation). Their genomic organization show great flexibility: some are individually or polycistronically transcribed, while others are encoded within introns of other genes. Here, we present an evolutionary analysis of the U49 gene in seven species. In all species analyzed, U49 contains the typical hallmarks of C and D box motifs, and a conserved 12-15 nt sequence complementary to rRNA that define them as homologs. In mouse, human, and Drosophila U49 is found encoded within introns of different genes, and in plants it is transcribed polycistronically from four different locations. In addition, U49 has two copies in two different introns of the RpL14 gene in Drosophila. The results indicate a substantial degree of duplication and translocation of the U49 gene in evolution. In light of its variable organization we discuss which of the two proposed mechanisms of rearrangement has acted upon the U49 snoRNA gene: chromosomal duplication or transposition through an RNA intermediate.     

Conserved sequence elements associated with exon skipping

One of the major forms of alternative splicing, which generates multiple mRNA isoforms differing in the precise combinations of their exon sequences, is exon skipping. While in constitutive splicing all exons are included, in the skipped pattern(s) one or more exons are skipped. The regulation of this process is still not well understood; so far, cis- regulatory elements (such as exonic splicing enhancers) were identified in individual cases. We therefore set to investigate the possibility that exon skipping is controlled by sequences in the adjacent introns. We employed a computer analysis on 54 sequences documented as undergoing exon skipping, and identified two motifs both in the upstream and downstream introns of the skipped exons. One motif is highly enriched in pyrimidines (mostly C residues), and the other motif is highly enriched in purines (mostly G residues). The two motifs differ from the known cis-elements present at the 5′ and 3′ splice site. Interestingly, the two motifs are complementary, and their relative positional order is conserved in the flanking introns. These suggest that base pairing interactions can underlie a mechanism that involves secondary structure to regulate exon skipping. Remarkably, the two motifs are conserved in mouse orthologous genes that undergo exon skipping.     

Different expression strategy: multiple intronic gene clusters of box H/ACA snoRNA in Drosophila melanogaster

The high degree of rRNA pseudouridylation in Drosophila melanogaster provides a good model for studying the genomic organization, structural and functional diversity of box H/ACA small nucleolar RNAs (snoRNAs). Accounting for both conserved sequence motifs and secondary structures, we have developed a computer-assisted method for box H/ACA snoRNA searching. Ten snoRNA clusters containing 42 box H/ACA snoRNAs were identified from D.melanogaster. Strikingly, they are located in the introns of eight protein-coding genes. In contrast to the mode of one snoRNA per intron so far observed in all animals, our results demonstrate for the first time a novel polycistronic organization that implies a different expression strategy for a box H/ACA snoRNA gene when compared to box C/D snoRNAs in D.melanogaster. Mutiple isoforms of the box H/ACA snoRNAs, from which most clusters are made up, were observed in D.melanogaster. The degree of sequence similarity between the isoforms varies from 99% to 70%, implying duplication events in different periods and a trend of enlarging the intronic snoRNA clusters. The variation in the functional elements of the isoforms could lead to partial alternation of snoRNA's function in loss or gain of rRNA complementary sequences and probably contributes to the great diversity of rRNA pseudouridylation in D.melanogaster.     

Evolutionary conservation suggests a regulatory function of AUG triplets in 5'-UTRs of eukaryotic genes

By comparing sequences of human, mouse and rat orthologous genes, we show that in 5′-untranslated regions (5′-UTRs) of mammalian cDNAs but not in 3′-UTRs or coding sequences, AUG is conserved to a significantly greater extent than any of the other 63 nt triplets. This effect is likely to reflect, primarily, bona fide evolutionary conservation, rather than cDNA annotation artifacts, because the excess of conserved upstream AUGs (uAUGs) is seen in 5′-UTRs containing stop codons in-frame with the start AUG and many of the conserved AUGs are found in different frames, consistent with the location in authentic non-coding sequences. Altogether, conserved uAUGs are present in at least 20–30% of mammalian genes. Qualitatively similar results were obtained by comparison of orthologous genes from different species of the yeast genus Saccharomyces. Together with the observation that mammalian and yeast 5′-UTRs are significantly depleted in overall AUG content, these findings suggest that AUG triplets in 5′-UTRs are subject to the pressure of purifying selection in two opposite directions: the uAUGs that have no specific function tend to be deleterious and get eliminated during evolution, whereas those uAUGs that do serve a function are conserved. Most probably, the principal role of the conserved uAUGs is attenuation of translation at the initiation stage, which is often additionally regulated by alternative splicing in the mammalian 5′-UTRs. Consistent with this hypothesis, we found that open reading frames starting from conserved uAUGs are significantly shorter than those starting from non-conserved uAUGs, possibly, owing to selection for optimization of the level of attenuation.     

Complete genomic structure of human rpS3: identification of functional U15b snoRNA in the fifth intron

Analysis of the complete genomic structure of the human ribosomal protein S3 (rpS3) gene revealed the presence of a functional U15b snoRNA gene in its intron. Human ribosomal protein S3 (rpS3) gene of 6115 bp long has been identified to contain six introns and seven exons in this study. The first and fifth introns of human S3 gene contain functional U15 snoRNA genes. Although Xenopus and Fugu counterparts also have six introns and seven exons, S3 gene of Fugu contains two functional U15 snoRNAs in the fourth and sixth introns and two pseudo genes for U15 snoRNAs in the first and fifth introns. In Xenopus S1 gene encoding ribosomal protein S3, however, three of its six introns contain U15 snoRNA gene sequence. Sequence comparison of the U15 genes from Xenopus, Fugu and human revealed that the regions involved in binding to 28S rRNA and the consensus sequence (C, D and D' boxes) for snoRNAs are highly conserved among those genes from these three species. Human U15a and U15b RNAs which are derived from the first and the fifth introns, respectively, have been identified to be functional by microinjection of human U15a and U15b snoRNAs into Xenopus oocyte. Northern blot and primer extension analyses confirm that human U15b snoRNA is expressed in vivo.     

Molecular evolution of the mammalian ribosomal protein gene, RPS14

Ribosomal protein S14 genes (RPS14) in eukaryotic species from protozoa to primates exhibit dramatically different intron-exon structures yet share homologous polypeptide-coding sequences. To recognize common features of RPS14 gene architectures in closely related mammalian species and to evaluate similarities in their noncoding DNA sequences, we isolated the intron-containing S14 locus from Chinese hamster ovary (CHO) cell DNA by using a PCR strategy and compared it with human RPS14. We found that rodent and primate S14 genes are composed of identical protein-coding exons interrupted by introns at four conserved DNA sites. However, the structures of corresponding CHO and human RPS14 introns differ significantly. Nonetheless, individual intron splice donor, splice acceptor, and upstream flanking motifs have been conserved within mammalian S14 homologues as well as within RPS14 gene fragments PCR amplified from other vertebrate genera (birds and bony fish). Our data indicate that noncoding, intronic DNA sequences within highly conserved, single-copy ribosomal protein genes are useful molecular landmarks for phylogenetic analysis of closely related vertebrate species.      

Bioinformatics analysis of plant orthologous introns: identification of an intronic tRNA-like sequence

Orthologous introns have identical positions relative to the coding sequence in orthologous genes of different species. By analyzing the complete genomes of five plants we generated a database of 40,512 orthologous intron groups of dicotyledonous plants, 28,519 orthologous intron groups of angiosperms, and 15,726 of land plants (moss and angiosperms). Multiple sequence alignments of each orthologous intron group were obtained using the Mafft algorithm. The number of conserved regions in plant introns appeared to be hundreds of times fewer than that in mammals or vertebrates. Approximately three quarters of conserved intronic regions among angiosperms and dicots, in particular, correspond to alternatively-spliced exonic sequences. We registered only a handful of conserved intronic ncRNAs of flowering plants. However, the most evolutionarily conserved intronic region, which is ubiquitous for all plants examined in this study, including moss, possessed multiple structural features of tRNAs, which caused us to classify it as a putative tRNA-like ncRNA. Intronic sequences encoding tRNA-like structures are not unique to plants. Bioinformatics examination of the presence of tRNA inside introns revealed an unusually long-term association of four glycine tRNAs inside the Vac14 gene of fish, amniotes, and mammals.     

Homologue disjunction in mouse oocytes requires proteolysis of securin and cyclin B1

Disjunction of pairs of homologous chromosomes during the first meiotic division (MI) requires anaphase-promoting complex (APC)-mediated activation of separase in budding yeast and Caenorhabditis elegans, but not Xenopus laevis. It is not clear which model best fits the mammalian system. Here we show that homologue disjunction in mouse oocytes is dependent on proteolysis of the separase inhibitor securin and the Cdk1 regulatory sub-unit cyclin B1. Proteolysis of both proteins was entirely dependent on their conserved destruction box (D-box) motifs, through which they are targeted to the APC. These data indicate that the mechanisms regulating homologue disjunction in mammalian oocytes are similar to those of budding yeast and C.elegans.     

Divergence of Mammalian Higher Order Chromatin Structure Is Associated with Developmental Loci

Several recent studies have examined different aspects of mammalian higher order chromatin structure – replication timing, lamina association and Hi-C inter-locus interactions — and have suggested that most of these features of genome organisation are conserved over evolution. However, the extent of evolutionary divergence in higher order structure has not been rigorously measured across the mammalian genome, and until now little has been known about the characteristics of any divergent loci present. Here, we generate a dataset combining multiple measurements of chromatin structure and organisation over many embryonic cell types for both human and mouse that, for the first time, allows a comprehensive assessment of the extent of structural divergence between mammalian genomes. Comparison of orthologous regions confirms that all measurable facets of higher order structure are conserved between human and mouse, across the vast majority of the detectably orthologous genome. This broad similarity is observed in spite of many loci possessing cell type specific structures. However, we also identify hundreds of regions (from 100 Kb to 2.7 Mb in size) showing consistent evidence of divergence between these species, constituting at least 10% of the orthologous mammalian genome and encompassing many hundreds of human and mouse genes. These regions show unusual shifts in human GC content, are unevenly distributed across both genomes, and are enriched in human subtelomeric regions. Divergent regions are also relatively enriched for genes showing divergent expression patterns between human and mouse ES cells, implying these regions cause divergent regulation. Particular divergent loci are strikingly enriched in genes implicated in vertebrate development, suggesting important roles for structural divergence in the evolution of mammalian developmental programmes. These data suggest that, though relatively rare in the mammalian genome, divergence in higher order chromatin structure has played important roles during evolution.