Search

tRNA-derived small RNA fragments (tRFs) are one class of small non-coding RNAs derived from transfer RNAs (tRNAs). tRFs play important roles in cellular processes and are involved in multiple cancers. High-throughput small RNA (sRNA) sequencing experiments can detect all the cellular expressed sRNAs, including tRFs. However, distinguishing genuine tRFs from RNA fragments generated by random degradation remains a major challenge. In this study, we developed an integrated web-based computing system, tRF2Cancer, to accurately identify tRFs from sRNA deep-sequencing data and evaluate their expression in multiple cancers. The binomial test was introduced to evaluate whether reads from a small RNA-seq data set represent tRFs or degraded fragments. A classification method was then used to annotate the types of tRFs based on their sites of origin in pre-tRNA or mature tRNA. We applied the pipeline to analyze 10 991 data sets from 32 types of cancers and identified thousands of expressed tRFs. A tool called 'tRFinCancer' was developed to facilitate the users to inspect the expression of tRFs across different types of cancers. Another tool called 'tRFBrowser' shows both the sites of origin and the distribution of chemical modification sites in tRFs on their source tRNA. The tRF2Cancer web server is available at http://rna.sysu.edu.cn/tRFfinder/.

Marek׳s disease virus (MDV) is a representative alpha herpes virus able to induce rapid-onset T-cell lymphoma in its natural host and regarded as an ideal model for the study of virus-induced tumorigenesis. Recent studies have shown that the mdv1-miR-M4-5p, a viral analog of cellular miR-155, is critical for MDV׳s oncogenicity. However, the precise mechanism whereby it was involved in MD lymphomagenesis remained unknown. We have presently identified the host mRNA targets of mdv1-miR-M4-5 and identified the latent TGF-β binding protein 1 (LTBP1) as a critical target for it. We found that during MDV infection, down-regulation of LTBP1 expression by mdv1-miR-M4-5p led to a significant decrease of the secretion and activation of TGF-β1, with suppression of TGF-β signaling and a significant activation of expression of c-Myc, a well-known oncogene which is critical for virus-induced tumorigenesis. Our findings reveal a novel and important mechanism of how mdv1-miR-M4-5p potentially contributes to MDV-induced tumorigenesis.

Inducing degradation of oncoproteins by small molecule compounds has the potential to avoid drug resistance and therefore deserves to be exploited for new therapies. Oridonin is a natural compound with promising antitumor efficacy that can trigger the degradation of oncoproteins; however, the direct cellular targets and underlying mechanisms remain unclear. Here we report that oridonin depletes BCR-ABL through chaperon-mediated proteasomal degradation in leukemia. Mechanistically, oridonin poses oxidative stress in cancer cells and directly binds to cysteines of HSF1, leading to the activation of this master regulator of the chaperone system. The resulting induction of HSP70 and ubiquitin proteins and the enhanced binding to CHIP E3 ligase hence target BCR-ABL for ubiquitin-proteasome degradation. Both wild-type and mutant forms of BCR-ABL can be efficiently degraded by oridonin, supporting its efficacy observed in cultured cells as well as mouse tumor xenograft assays with either imatinib-sensitive or -resistant cells. Collectively, our results identify a novel mechanism by which oridonin induces rapid degradation of BCR-ABL as well as a novel pharmaceutical activator of HSF1 that represents a promising treatment for leukemia.

Coordinated assembly of the ribosome is essential for proper translational activity in eukaryotic cells. It is therefore critical to coordinate the expression of components of ribosomal programs with the cell's nutritional status. However, coordinating expression of these components is poorly understood. Here, by combining experimental and computational approaches, we systematically identified box C/D snoRNAs in four fission yeasts and found that the expression of box C/D snoRNA and ribosomal protein (RP) genes were orchestrated by a common Homol-D box, thereby ensuring a constant balance of these two genetic components. Interestingly, such transcriptional coregulations could be observed in most Ascomycota species and were mediated by different cis-regulatory elements. Via the reservation of cis elements, changes in spatial configuration, the substitution of cis elements, and gain or loss of cis elements, the regulatory networks of box C/D snoRNAs evolved to correspond with those of the RP genes, maintaining transcriptional coregulation between box C/D snoRNAs and RP genes. Our results indicate that coregulation via common cis elements is an important mechanism to coordinate expression of the RP and snoRNA genes, which ensures a constant balance of these two components.

Recent studies have shown that endogenous small RNAs regulate a variety of biological processes during vertebrate development; however, little is known about the role of small RNAs in regulating developmental signaling pathways during early embryogenesis. In this study, we applied Illumina sequencing to characterize an unexpected endogenous small RNA catalog and demonstrated a dramatic transition from transposon-derived piRNA-like small RNAs (pilRNAs) to microRNAs (miRNAs) in pre- and post-gastrula chicken embryos. The comprehensive expression profile of chicken miRNAs at the pre- and post-gastrula stages revealed that most known and new miRNAs were dynamically regulated during development. In addition to embryonic stem cell-related miRNAs, Gene Ontology (GO) analysis showed that miRNAs enriched in early stage chicken embryos targeted multiple signal transduction pathways associated with the reproductive process and embryogenesis, including Wnt and TGF-β, which specifies the neural fate of blastodermal cells. Intriguingly, a large cohort of pilRNAs primarily derived from the active and most abundant transposable elements (TEs) were enriched in chicken stage X blastoderms. Within stage X blastoderms, pilRNAs were specifically localized to the primordial germ cells (PGCs), indicating their post-zygotic origin. Together, these findings imply a role for small RNAs in gastrulation in early stage chicken embryos.

microRNAs (miRNAs) represent an abundant group of small regulatory non-coding RNAs in eukaryotes. The emergence of Next-generation sequencing (NGS) technologies has allowed the systematic detection of small RNAs (sRNAs) and de novo sequencing of genomes quickly and with low cost. As a result, there is an increased need to develop fast miRNA prediction tools to annotate miRNAs from various organisms with a high level of accuracy, using the genome sequence or the NGS data. Several miRNA predictors have been proposed to achieve this purpose. However, the accuracy and fitness for multiple species of existing predictors needed to be improved. Here, we present a novel prediction tool called mirExplorer, which is based on an integrated adaptive boosting method and contains two modules. The first module named mirExplorer-genome was designed to de novo predict pre-miRNAs from genome, and the second module named mirExplorer-NGS was used to discover miRNAs from NGS data. A set of novel features of pre-miRNA secondary structure and miRNA biogenesis has been extracted to distinguish real pre-miRNAs from pseudo ones. We used outer-ten-fold cross-validation to verify the mirExplorer-genome computation, which obtained a specificity of 95.03% and a sensitivity of 93.71% on human data. This computation was made on test data from 16 species, and it achieved an overall accuracy of 95.53%. Systematic outer-ten-fold cross-validation of the mirExplorer-NGS model achieved a specificity of 98.3% and a sensitivity of 97.72%. We found that the good performance of the mirExplorer-NGS model was upheld across species from vertebrates to plants in test datasets. The mirExplorer is available as both web server and software package at http://biocenter.sysu.edu.cn/mir/.

International audience; Plant genomes have undergone multiple rounds of duplications that contributed massively to the growth of gene families. The structure of resulting families has been studied in depth for protein-coding genes. However, little is known about the impact of duplications on noncoding RNA (ncRNA) genes. Here we perform a systematic analysis of duplicated regions in the rice genome in search of such ncRNA repeats. We observe that, just like their protein counterparts, most ncRNA genes have undergone multiple duplications that left visible sequence conservation footprints. The extent of ncRNA gene duplication in plants is such that these sequence footprints can be exploited for the discovery of novel ncRNA gene families on a large scale. We developed an SVM model that is able to retrieve likely ncRNA candidates among the 100,000+ repeat families in the rice genome, with a reasonably low false-positive discovery rate. Among the nearly 4000 ncRNA families predicted by this means, only 90 correspond to putative snoRNA or miRNA families. About half of the remaining families are classified as structured RNAs. New candidate ncRNAs are particularly enriched in UTR and intronic regions. Interestingly, 89% of the putative ncRNA families do not produce a detectable signal when their sequences are compared to another grass genome such as maize. Our results show that a large fraction of rice ncRNA genes are present in multiple copies and are species-specific or of recent origin. Intragenome comparison is a unique and potent source for the computational annotation of this major class of ncRNA.

MicroRNAs are involved in cancer-related processes. The microRNA-21(miR-21) has been identified as the only miRNA over-expressed in a wide variety of cancers, including cervical cancer. However, the function of miR-21 is unknown in cervical carcinomas. In this study, we found that the inhibition of miR-21 in HeLa cervical cancer cells caused profound suppression of cell proliferation, and up-regulated the expression of the tumor suppressor gene PDCD4. We also provide direct evidence that PDCD4-3'UTR is a functional target of miR-21 and that the 18bp putative target site can function as the sole regulatory element in HeLa cells. These results suggest that miR-21 may play an oncogenic role in the cellular processes of cervical cancer and may serve as a target for effective therapies.

Ribose methylation of eukaryotic rRNA is directed by box C/D small nucleolar RNAs (snoRNAs), which pinpoint the nucleotide to be methylated in specific position within the rRNA sequence. Here, we report the identification of a novel double-guide C/D box snoRNA termed snR88 that directs methylation of two previously undetermined sites in 25S rRNA from the fission yeast. Knockout of the predicted TATA box of the snR88 gene resulted in the complete blocking of its expression, showing that snR88 is an independently transcribed gene and dispensable for yeast viability. The depletion of snR88 abolished 25S rRNA methylation at U2304 and U2497 simultaneously. Interestingly, an unusual pause of reverse transcription at U2495 was observed, which implies an unknown structure of 25S rRNA related to ribose methylation at U2497 in the fission yeast.

Spermatogenesis in mammals is characterized by two waves of piRNA expression: one corresponds to classic piRNAs responsible for silencing retrotransponsons and the second wave is predominantly derived from nontransposon intergenic regions in pachytene spermatocytes, but the function of these pachytene piRNAs is largely unknown. Here, we report the involvement of pachytene piRNAs in instructing massive mRNA elimination in mouse elongating spermatids (ES). We demonstrate that a piRNA-induced silencing complex (pi-RISC) containing murine PIWI (MIWI) and deadenylase CAF1 is selectively assembled in ES, which is responsible for inducing mRNA deadenylation and decay via a mechanism that resembles the action of miRNAs in somatic cells. Such a highly orchestrated program appears to take full advantage of the enormous repertoire of diversified targeting capacity of pachytene piRNAs derived from nontransposon intergenic regions. These findings suggest that pachytene piRNAs are responsible for inactivating vast cellular programs in preparation for sperm production from ES.

Most, if not all, of snoRNAs in mammals are intron-encoded, implying the expressional and functional relativeness between the snoRNA and their hosts. By computational analysis of an intron database extracted from 65 known mouse imprinted genes, three novel orphan box C/D snoRNAs were identified from Irm gene which is maternally expressed and related to human disorders. The snoRNAs were positively detected and found to express in all the mouse tissues except kidney. The imprinted snoRNAs exhibit stringent structures, but quite variable in locations at their host introns, suggesting their maturation probably through a splicing independent manner. We characterized Irm as a new kind of snoRNA host gene which has no protein-coding capacity and no 5'TOP structure in its mRNA. The newly identified snoRNAs appear mouse-specific, however, their function remains to be elucidated.

The extensive post-transcriptional modification of U2 snRNA is required for spliceosome assembly and pre-mRNA splicing in vertebrates. However, the rare modification of U2 snRNA in yeast implies a different mechanism for regulating spliceosome biogenesis in single-celled eukaryotes. To understand the evolutionary pattern of U2 snRNA methylation, we determined for the first time, the 2'-O-methylations of U2 snRNA in Oryza sativa, Arabidopsis thaliana and Drosophila melanogaster, and revealed two methylations which are conserved in a crucial region of U2 snRNA in plants. Interestingly, one of the methylations, U2-Cm29 is also methylated in D. melanogaster, but not in vertebrates. According to the methylation of U2-C29, computational analysis of databases identified three canonical box C/D snoRNAs, named OsmgU2-29, AtmgU2-29 and DmmgU2-28, as small methylation guides of U2 snRNA from O. sativa, A. thaliana and D. melanogaster, respectively. Although very divergent in their sequence, the three snoRNAs exhibit in common an 11 nucleotide-long sequence complementarity to corresponding U2 snRNA, implying a functional constraint on the modification during evolution. Interestingly, a nucleotide is found to be inserted both in U2 snRNA and DmmgU2-28 and maintains a perfect match of duplex specifying the methylation of C28 in Drosophila U2 snRNA. This is the first time a new model is being provided for compensatory mutations between a small guide RNA and its target by nucleotide insertion, instead of the known nucleotide substitution. In contrast to small Cajal body-specific RNAs (scaRNAs), the snoRNAs are similar to the reported singlet guide RNAs and are known to localize in nucleolus.

Small nucleolar RNAs (snoRNAs) are an abundant group of noncoding RNAs mainly involved in the post-transcriptional modifications of rRNAs in eukaryotes. In this study, a large-scale genome-wide analysis of the two major families of snoRNA genes in the fruit fly Drosophila melanogaster has been performed using experimental and computational RNomics methods. Two hundred and twelve gene variants, encoding 56 box H/ACA and 63 box C/D snoRNAs, were identified, of which 57 novel snoRNAs have been reported for the first time. These snoRNAs were predicted to guide a total of 147 methylations and pseudouridylations on rRNAs and snRNAs, showing a more comprehensive pattern of rRNA modification in the fruit fly. With the exception of nine, all the snoRNAs identified to date in D. melanogaster are intron encoded. Remarkably, the genomic organization of the snoRNAs is characteristic of 8 dUhg genes and 17 intronic gene clusters, demonstrating that distinct organizations dominate the expression of the two families of snoRNAs in the fruit fly. Of the 267 introns in the host genes, more than half have been identified as host introns for coding of snoRNAs. In contrast to mammals, the variation in size of the host introns is mainly due to differences in the number of snoRNAs they contain. These results demonstrate the extensive utilization of introns for coding of snoRNAs in the host genes and shed light on further research of other noncoding RNA genes in the large introns of the Drosophila genome.

From a specialized cDNA library of Giardia lamblia, 20 snoRNA-like RNAs, including 16 box C/D sRNAs and four box H/ACA sRNAs, were first identified. The sRNAs were predicted to guide a total of 11 2'-O-methylation and four pseudouridylation sites on the G. lamblia rRNAs, respectively. By using primer extension assay, seven methylation sites were precisely mapped in the G. lamblia 16S rRNA, despite its high GC content. All of the sRNA genes locate on the small intergenic regions of the G. lamblia genome and seem to be independently transcribed from their own promoters. Particularly, a cluster composed of GlsR17 and GlsR18 genes is transcribed as a dicistronic precursor, implying a mechanism of endonuclease cleavage for the maturation of the two sRNAs. The systematic identification of the sRNAs in G. lamblia has provided valuable information about the characteristics of the two major families of small guide RNAs in one of the most primitive eukaryotes and would contribute to the understanding of the evolution of small non-messenger RNA genes from prokaryotes to eukaryotes.

In this study, we have determined the complete sequence of the mitochondrial genome of an Old World tree frog Polypedates megacephalus (Anura: Rhacophoridae) by using a long polymerase chain reaction (PCR) technique and shotgun strategy of sequencing. The entire mtDNA sequence is 16,473 nt long with a novel mitogenomic gene organization in amphibians. Unlike other neobatrachian frogs, the transfer ribonucleic acid (tRNA)-LeuCUN and tRNA-Thr genes exchange their positions in P. megacephalus and form a Thr-LeuCUN-Pro-Phe tRNA gene tetrad. Moreover, we found that the ATP8 gene was replaced by a noncoding sequence of 853 nt long and that the ND5 gene was absent in the new mitogenome. These peculiar features of P. megacephalus mtDNA were further studied among related anuran species by PCR amplification. The new sequence data was used to assess the phylogenetic relationships of the three living amphibian orders using neighbor-joining, maximum likelihood, and Bayesian methods. In agreement with most morphological studies, phylogenetic analyses of a whole mitochondrial genome data set suggest a close relationship between salamanders and frogs. Moreover, using a molecular clock-independent Bayesian approach for inferring dating information from molecular phylogenies, we have provided a rough timescale for living amphibian evolution. This timescale provides a working framework for future paleontological researches on amphibian evolution and improves our understanding of the evolutionary history of modern amphibians.

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.

The mitochondrial genome of the Chinese giant salamander Andrias davidianus was isolated using the long-and-accurate polymerase chain reaction (LA PCR) method. The sequencing work adopted the shotgun strategy accompanying with seven internal primers to cover the gaps where overlapping clones were not available. The entire mtDNA sequence is 16,503 bp long, with a gene content of 13 protein-coding, two ribosomal RNA and 22 transfer RNA genes, and order identical to that observed in most other vertebrates except for an additional 318 bp non-coding sequence between tRNA-Thr and tRNA-Pro genes. In order to carry out molecular phylogenetic analyses, all 13 protein sequences deduced from whole mitochondrial genomes for eight vertebrate species (six amphibians, two lobe-finned fishes) were combined to a single data set. This data set was refined by the program Gblocks using a stringent parameter setting and then subjected to MP, ML and NJ analyses. Thus phylogenetic relationships among living amphibians were discussed.

Using a powerful computer-assisted analysis strategy, a large-scale search of small nucleolar RNA (snoRNA) genes in the recently released draft sequence of the rice genome was carried out. This analysis identified 120 different box C/D snoRNA genes with a total of 346 gene variants, which were predicted to guide 135 2'-O-ribose methylation sites in rice rRNAs. Though not exhaustive, this analysis has revealed that rice has the highest number of known box C/D snoRNAs among eukaryotes. Interestingly, although many snoRNA genes are conserved between rice and Arabidopsis, almost half of the identified snoRNA genes are rice specific, which may highlight further the differences in rRNA methylation patterns between monocotyledons and dicotyledons. In addition to 76 singletons, 70 clusters involving 270 snoRNA genes were also found in rice. The large number of the novel snoRNA polycistrons found in the introns of rice protein-coding genes is in contrast to the one-snoRNA-per-intron organization of vertebrates and yeast, and of Arabidopsis in which only a few intronic snoRNA gene clusters were identified. Furthermore, due to a high degree of gene duplication, rice snoRNA genes are clearly redundant and exhibit great sequence variation among isoforms, allowing generation of new snoRNAs for selection. Thus, the large snoRNA gene family in plants can serve as an excellent model for a rapid and functional evolution.

Small nucleolar RNAs (snoRNAs) are a well-characterized family of non-coding RNAs whose main function is rRNA modification. The diversity and complexity of this gene family continues to expand with the discovery of snoRNAs with non-rRNA or unknown targets. Plants contain more snoRNAs than other eukaryotes and have developed novel expression and processing strategies. The increased number of modifications, which will influence ribosome function, and the novel modes of expression might reflect the environmental conditions to which plants are exposed. Polyploidy and chromosomal rearrangements have generated multiple copies of snoRNA genes, allowing the generation of new snoRNAs for selection. The large snoRNA family in plants is an ideal model for investigation of mechanisms of evolution of gene families in plants.

Small RNAs (sRNAs), including microRNAs and endogenous siRNAs (endo-siRNAs), regulate most important biologic processes in eukaryotes, such as cell division and differentiation. Although sRNAs have been extensively studied in various eukaryotes, the role of sRNAs in the early emergence of eukaryotes is unclear. To address these questions, we deep sequenced the sRNA transcriptome of four different stages in the differentiation of Giardia lamblia, one of the most primitive eukaryotes. We identified a large number of endo-siRNAs in this fascinating parasitic protozoan and found that they were produced from live telomeric retrotransposons and three genomic regions (i.e., endo-siRNA generating regions [eSGRs]). eSGR-derived endo-siRNAs were proven to target mRNAs in trans. Gradual up-regulation of endo-siRNAs in the differentiation of Giardia suggested that they might be involved in the regulation of this process. This hypothesis was supported by the impairment of the differentiation ability of Giardia when GLDICER, essential for the biogenesis of endo-siRNAs, was knocked down. Endo-siRNAs are not the only sRNA regulators in Giardia differentiation, because a great number of tRNAs-derived sRNAs showed more dramatic expression changes than endo-siRNAs in this process. We totally identified five novel kinds of tRNAs-derived sRNAs and found that the biogenesis in four of them might be correlated with that of stress-induced tRNA-derived RNA (sitRNA), which was discovered in our previous studies. Our studies reveal an unexpected complex panorama of sRNA in G. lamblia and shed light on the origin and functional evolution of eukaryotic sRNAs.

Transcriptome studies have revealed that many non-coding RNAs (ncRNAs) are located near the 3' sense terminus of protein-coding genes. However, the transcription and function of these RNAs remain elusive. Here, we identify a 3' sense termini-associated sRNA (TASR) downstream of rpl26 in Schizosaccharomyces pombe (S. pombe). Structure and function assays indicate that the TASR is an H/ACA box snoRNA required for 18S rRNA pseudouridylation at U121 and U305 sites and is therefore a cognate of snR49 from the budding yeast. Transcriptional studies show that pre-snR49 overlaps most of the coding sequence (CDS) of rpl26. Using scanning deletion analysis within promoter region, we show that the rpl26 promoter is required for the 3' TASR transcription. Interestingly, chromosomal synteny of rpl26-snR49 is found in the Schizosaccharomyces groups. Taken together, we have revealed a new transcriptional mechanism for 3' sense TASRs, which are transcribed by the same promoter as their upstream protein genes. These results further suggest that the origin and function of 3' sense ncRNAs are associated with upstream genes in higher eukaryotes.

Trypanosoma brucei, a pathogen of human and domestic animals, is an early evolved parasitic protozoan with a complex life cycle. Most genes of this parasite are post-transcriptionally regulated. However, the mechanisms and the molecules involved remain largely unknown. We have deep-sequenced the small RNAs of two life stages of this parasite—the bloodstream form and the procyclic form. Our results show that the small RNAs of T. brucei could derive from multiple sources, including NATs (natural antisense transcripts), tRNAs, and rRNAs. Most of these small RNAs in the two stages were found to share uniform characteristics. However, our results demonstrate that their variety and expression show significant differences between different stages, indicating possible functional differentiation. Dicer-knockdown evidence further proved that some of the small interfering RNAs (siRNAs) could regulate the expression of genes. Based on the genome-wide analysis of the small RNAs in the two stages of T. brucei, our results not only provide evidence to study their differentiation but also shed light on questions regarding the origins and evolution of small RNA-based mechanisms in early eukaryotes.

MicroRNAs (miRNAs) are endogenous non-coding RNAs (ncRNAs) of approximately 22 nt that regulate the expression of a large fraction of genes by targeting messenger RNAs (mRNAs). However, determining the biologically significant targets of miRNAs is an ongoing challenge. In this chapter, we describe how to identify miRNA-target interactions and miRNA regulatory networks from high-throughput deep sequencing, CLIP-Seq (HITS-CLIP, PAR-CLIP) and degradome sequencing data using starBase platforms. In starBase, several web-based and stand-alone computational tools were developed to discover Argonaute (Ago) binding and cleavage sites, miRNA-target interactions, perform enrichment analysis of miRNA target genes in Gene Ontology (GO) categories and biological pathways, and identify combinatorial effects between Ago and other RNA-binding proteins (RBPs). Investigating target pathways of miRNAs in human CLIP-Seq data, we found that many cancer-associated miRNAs modulate cancer pathways. Performing an enrichment analysis of genes targeted by highly expressed miRNAs in the mouse brain showed that many miRNAs are involved in cancer-associated MAPK signaling and glioma pathways, as well as neuron-associated neurotrophin signaling and axon guidance pathways. Moreover, thousands of combinatorial binding sites between Ago and RBPs were identified from CLIP-Seq data suggesting RBPs and miRNAs coordinately regulate mRNA transcripts. As a means of comprehensively integrating CLIP-Seq and Degradome-Seq data, the starBase platform is expected to identify clinically relevant miRNA-target regulatory relationships, and reveal multi-dimensional post-transcriptional regulatory networks involving miRNAs and RBPs. starBase is available at http://starbase.sysu.edu.cn/ .

U21 is an intron-encoded snoRNA in vertebrates which contains a 13-nt tract of complementarity to an invariant sequence in eukaryotic 28S rRNA. Here, we report the characterization of its Drosophila melanogaster homolog which is the first case of an intron-encoded snoRNA in an invertebrate metazoan. In D. melanogaster , U21 is encoded within the ARF-1 (ADP ribosylation factor 1) gene, whereas in chicken and mammals it is found in the ribosomal protein L5 gene. In D. melanogaster , like in vertebrates, U21 is devoid of a 5′ trimethylguanosine cap, thus, likely resulting from processing of intronic RNA. The only portion of U21 sequence preserved between D. melanogaster and vertebrates, in addition to the hallmark box C and box D motifs, corresponds precisely to the 13-nt complementary to rRNA, pointing to an important role of the pairing of U21 to pre-rRNA.

cyclin D3 (CCND3) is one of the three D-type cyclins that regulate the G1/S phase transition of the cell cycle. Expression of CCND3 is observed in nearly all proliferating cells; however, the presence of high levels of CCND3 has been linked to a poor prognosis for several types of cancer. Therefore, further mechanistic studies on the regulation of CCND3 expression are urgently needed to provide therapeutic implications. In this study, we report that a conserved RNA G-quadruplex-forming sequence (hereafter CRQ), located in the 5′ UTR of mammalian CCND3 mRNA, is able to fold into an extremely stable, intramolecular, parallel G-quadruplex in vitro. The CRQ G-quadruplex dramatically reduces the activity of a reporter gene in human cell lines, but it has little impact on its mRNA level, indicating a translational repression. Moreover, the CRQ sequence in its natural context inhibits translation of CCND3. Disruption of the G-quadruplex structure by G/U-mutation or deletion results in an elevated expression of CCND3 and an increased phosphorylation of Rb, a downstream target of CCND3, which promotes progression of cells through the G1 phase. Our results add to the growing understanding of the regulation of CCND3 expression and provide a potential therapeutic target for cancer treatment.

Background The Siberian salamander (Ranodon sibiricus), distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations. Methodology/Principal Findings We used two genetic markers (mtDNA and microsatellites) to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC) method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines. Conclusion/Significance Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.

International audience; Small RNAs constitute a new and unanticipated layer of gene regulation present in the three domains of life. In plants, all organs are ultimately derived from a few pluripotent stem cells localized in specialized structures called apical meristems. The development of meristems involves a coordinated balance between undifferentiated growth and differentiation, a phenomenon requiring a tight regulation of gene expression. We used in vitro cultured embryogenic calli as a model to investigate the roles of meristem-associated small RNAs. Using high throughput sequencing, we sequenced 20 million short reads with size of 18-30nt from rice undifferentiated and differentiated calli. We confirmed 50 known microRNA families, representing one third of annotated rice microRNAs. Using a specific computational pipeline for plant microRNA identification, we identified 24 novel microRNA families. Among them, 53 microRNA or microRNA* sequences appear to vary in expression between differentiated and undifferentiated calli, suggesting a role in meristem development. Our analysis also revealed a new class of plant small RNAs derived from 5' or 3' ends of mature tRNA analogous to the tRFs in human cancer cell. We independently verified the expression of these small RNAs from 5' end of mature tRNA using qRT-PCR.

BACKGROUND: The insulin-like growth factor (IGF) signaling pathway has long been established as playing critical roles in skeletal muscle development. However, the underlying regulatory mechanism is poorly understood. Recently, a large family of small RNAs, named microRNAs (miRNAs), has been identified as key regulators for many developmental processes. Because miRNAs participate in the regulation of various signaling pathways, we hypothesized that miRNAs may be involved in the regulation of IGF signaling in skeletal myogenesis. METHODOLOGY/PRINCIPAL FINDINGS: In the present study, we determined that the cell-surface receptor IGF-1R is directly regulated by a muscle-specific miRNA, microRNA-133 (miR-133). A conserved and functional binding site for miR-133 was identified in the 3'untranslated region (3'UTR) of IGF-1R. During differentiation of C2C12 myoblasts, IGF-1R protein, but not messenger RNA (mRNA) expression, was gradually reduced, concurrent with the upregulation of miR-133. Overexpression of miR-133 in C2C12 cells significantly suppressed IGF-1R expression at the posttranscriptional level. We also demonstrated that both overexpression of miR-133 and knockdown of IGF-1R downregulated the phosphorylation of Akt, the central mediator of the PI3K/Akt signaling pathway. Furthermore, upregulation of miR-133 during C2C12 differentiation was significantly accelerated by the addition of IGF-1. Mechanistically, we found that the expression of myogenin, a myogenic transcription factor reported to transactivate miR-133, was increased by IGF-1 stimulation. CONCLUSION/SIGNIFICANCE: Our results elucidate a negative feedback circuit in which IGF-1-stimulated miR-133 in turn represses IGF-1R expression to modulate the IGF-1R signaling pathway during skeletal myogenesis. These findings also suggest that miR-133 may be a potential therapeutic target in muscle diseases.

Imprinted noncoding RNAs (ncRNAs) are expressed mono-allelically in a parent-of-origin-dependent manner, which is mainly evident in mammals. Lying at a crossroad between imprinted genes and ncRNAs, imprinted ncRNAs show distinct features. They are likely to function in nontraditional ways compared to non-imprinted ncRNAs, and are much more responsible for the mechanism of genomic imprinting compared to imprinted protein-coding genes. An increasing number of human diseases have been shown to be related to abnormalities in imprinted ncRNAs. Due to their functional importance, many studies focusing on imprinted ncRNAs have been published in recent years; however, there is no systematic collection or description of imprinted ncRNAs and the rapidly growing knowledge is scattered in various places. Here, we describe a new database, ncRNAimprint, which serves as a comprehensive resource center for mammalian imprinted ncRNAs. A catalog of imprinted ncRNAs, including snoRNAs, microRNAs, piRNAs, siRNAs, antisense ncRNAs, and mRNA-like ncRNAs, was annotated in detail using information extracted from relevant literature and databases. Comprehensive collections of imprinted ncRNA-related diseases, imprinting control regions (ICRs), and imprinted regions were manually compiled to provide resources for current research focusing on imprinted ncRNAs. Small RNA deep sequencing reads that fully matched within imprinted regions were also included to offer useful evidence in detecting novel imprinted ncRNAs and to aid in analyzing expression patterns of known imprinted ncRNAs. A search page including four effective search forms and two graphical browsers was created for rapid retrieval and visualization of these data. The imprinted ncRNA database is freely accessible at http://rnaqueen.sysu.edu.cn/ncRNAimprint.

BACKGROUND: MicroRNAs (miRNAs) are approximately 22-nt small non-coding regulatory RNAs that have generally been considered to regulate gene expression at the post-transcriptional level in the cytoplasm. However, recent studies have reported that some miRNAs localize to and function in the nucleus. METHODOLOGY/PRINCIPAL FINDINGS: To determine the number of miRNAs localized to the nucleus, we systematically investigated the subcellular distribution of small RNAs (sRNAs) by independent deep sequencing sequenced of the nuclear and cytoplasmic pools of 18- to 30-nucleotide sRNAs from human cells. We identified 339 nuclear and 324 cytoplasmic known miRNAs, 300 of which overlap, suggesting that the majority of miRNAs are imported into the nucleus. With the exception of a few miRNAs evidently enriched in the nuclear pool, such as the mir-29b, the ratio of miRNA abundances in the nuclear fraction versus in the cytoplasmic fraction vary to some extent. Moreover, our results revealed that a large number of tRNA 3' trailers are exported from the nucleus and accumulate in the cytoplasm. These tRNA 3' trailers accumulate in a variety of cell types, implying that the biogenesis of tRNA 3' trailers is conserved and that they have a potential functional role in vertebrate cells. CONCLUSION/SIGNIFICANCE: Our results provide the first comprehensive view of the subcellular distribution of diverse sRNAs and new insights into the roles of miRNAs and tRNA 3' trailers in the cell.

The RNA minihelix, a proposed tRNA precursor, exhibits tRNA-like properties. Sequence-specific RNA minihelices can inhibit cell growth probably due to their binding to the cognate tRNA and naturally occurring non-tRNA substrates for aminoacylation may serve a similar purpose. Thus far, no natural RNA minihelices have been found. In the present study, we found a novel small RNA of 32 nucleotides, which is expressed abundantly in all rat tissues tested. Distinct from all of known endogenous small RNAs, this small RNA (temporarily named as tpl-sRNA) can form an RNA minihelix containing a stem-loop domain followed by ACCA. tpl-sRNA is encoded by the light-strand replication origin (Ori L) of the rat mitochondrial genome, and the 3′-terminal CCA of tpl-sRNA is post-transcriptionally added. Moreover, tpl-sRNA is chargeable in vivo. Our study demonstrates for the first time an endogenous small RNA that resembles an ancestral tRNA and exhibits some tRNA-like properties in mammals.

Chlamydomonas reinhardtii is a unicellular green alga, the lineage of which diverged from that of land plants >1 billion years ago. Using the powerful small nucleolar RNA (snoRNA) mining platform to screen the C. reinhardtii genome, we identified 322 snoRNA genes grouped into 118 families. The 74 box C/D families can potentially guide methylation at 96 sites of ribosomal RNAs (rRNAs) and snRNAs, and the 44 box H/ACA families can potentially guide pseudouridylation at 62 sites. Remarkably, 242 of the snoRNA genes are arranged into 76 clusters, of which 77% consist of homologous genes produced by small local tandem duplications. At least 70 snoRNA gene clusters are found within introns of protein-coding genes. Although not exhaustive, this analysis reveals that C. reinhardtii has the highest number of intronic snoRNA gene clusters among eukaryotes. The prevalence of intronic snoRNA gene clusters in C. reinhardtii is similar to that of rice but in contrast with the one-snoRNA-per-intron organization of vertebrates and fungi and with that of Arabidopsis thaliana in which only a few intronic snoRNA gene clusters were identified. This analysis of C. reinhardtii snoRNA gene organization shows the functional importance of introns in a single-celled organism and provides evolutionary insight into the origin of intron-encoded RNAs in the plant lineage.

During evolution, the potential for sequence (and length) variation of large-subunit rRNA has been mostly restricted over 12 divergent domains (termed D1-D12) interspersed along the molecule. Here, we have focused our attention onto the D3 divergent domain, through a detailed analysis of its pattern of variation in the phylogeny, both in terms of primary and secondary structures. We have systematically compared all the procaryotic and eucaryotic sequences published so far (i.e. 36 species), together with a series of 10 additional eucaryotic specimens, which were determined by direct RNA sequencing. Secondary structures supported by comparative evidence have been derived for archaebacteria, eubacteria and eucaryotes respectively, which shows that the D3 domain contains a subset of universally conserved structural features interspersed with four variable subdomains. Within the four portions where a structural diversification has taken place, elementary structures specific of large phylogenetic groups can be identified. Remarkably such diversified structures appear to be preserved despite sequence divergence, suggesting they correspond to functionally important structures. Accordingly, the mode of sequence variation of the D3 domain suggests this region of the molecule may encode elementary functions of rRNA which could have significantly diversified during the evolution of the major groups of organisms.

In vitro cultured embryogenic callus was employed as a model to investigate microRNAs (miRNAs) associated with embryogenesis and post-embryonic development. Thirty-one miRNAs including 16 novel species were identified from a large number of small RNAs which were cloned from both differentiated and undifferentiated rice embryogenic calli. Four target genes of the miRNAs were further validated. A set of the miRNAs, including miR397 and miR156, exhibited intriguing expression patterns during the transition from undifferentiated to differentiated calli. By exploiting the correlations between the differential expression patterns of these miRNAs and their targets, the regulatory roles of the miRNAs on meristem maintenance and embryogenesis were indicated.

Establishing the relationships among modern amphibians (lissamphibians) and their ancient relatives is necessary for our understanding of early tetrapod evolution. However, the phylogeny is still intractable because of the highly specialized anatomy and poor fossil record of lissamphibians. Paleobiologists are still not sure whether lissamphibians are monophyletic or polyphyletic, and which ancient group (temnospondyls or lepospondyls) is most closely related to them. In an attempt to address these problems, eight mitochondrial genomes of living amphibians were determined and compared with previously published amphibian sequences. A comprehensive molecular phylogenetic analysis of nucleotide sequences yields a highly resolved tree congruent with the traditional hypotheses (Batrachia). By using a molecular clock-independent approach for inferring dating information from molecular phylogenies, we present here the first molecular timescale for lissamphibian evolution, which suggests that lissamphibians first emerged about 330 million years ago. By observing the fit between molecular and fossil times, we suggest that the temnospondyl-origin hypothesis for lissamphibians is more credible than other hypotheses. Moreover, under this timescale, the potential geographic origins of the main living amphibian groups are discussed: (i) advanced frogs (neobatrachians) may possess an Africa-India origin; (ii) salamanders may have originated in east Asia; (iii) the tropic forest of the Triassic Pangaea may be the place of origin for the ancient caecilians. An accurate phylogeny with divergence times can be also helpful to direct the search for "missing" fossils, and can benefit comparative studies of amphibian evolution.

Heterosmilax japonica Kunth is well recognized for its diuretic effects in China. However, little is known about its endophytic fungi. In this study, microbial communities inhabiting the stems of H. japonica in spring and summer were investigated by light microscopy and cultivation-independent approaches, such as RFLP analysis and sequencing of rDNA ITS library. Molecular phylogenetic analysis showed that a broad spectrum of fungi, including Mycosphaerella, Phomopsis, Aureobasidium, Cladosporium, Glomerella, Botryosphaeria, Guignardia, is able to colonize the plants internally. Particularly, several rDNA sequences determined in this study like YJ4-61 are not specifically affiliated with any currently documented fungal sequences in the public database. Several sequence types, such as YJ4-9 and YJ4-70, are significantly similar to some uncultured environmental samples. Furthermore, our result also showed that the samples collected in spring harbored more abundant endophytic populations than that in summer, implying a seasonal fluctuation for the endophytes in H. japonica.

MicroRNAs (miRNAs) are present in both plant and animal kingdoms and represents a growing family of non-coding RNAs. These tiny RNAs act as small guides and direct negative regulations usually in the process of development through sequence complementarity to target mRNAs. Although a large number of miRNAs have been identified from various animals, so far plant miRNA studies have focused mainly on Arabidopsis. Here we describe the identification of 20 miRNAs from a rice cDNA library. All the miRNAs were presumably processed from precursors with stem-loop structures and were positively detected in rice cells from at least one tissue, some of which showed tissue-specific expression. Twenty-three unique rice genes were identified to be feasible targets for seven rice miRNAs, including four members of Scarecrow-like transcription factor, the targets of miR-39 that had been characterized in Arabidopsis. Lacking long complementarity, the regulatory targets of 13 miRNAs remain to be further investigated. A possible mechanism of translational repressor for plant miRNAs that lack perfect complementarity to target mRNAs is discussed.

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 .

The Amphibia were the dominant land Vertebrates in the Carboniferous and were certainly the stock from which the reptiles and in turn the mammals and birds evolved. Although two major groups of extinct amphibians, the Labyrinthodontia and the Lepospondyli did not last beyond the Triassic, amphibian lines evidently continued and are represented today by the three distinctly different groups of modern Amphibia (Lissamphibia), the Caudata, Anura, and Gymnophiona. The monophyly of the Lissamphibia and their phylogenetic relationships are still debated. The most widely accepted hypothesis, based on morphological data, supports the monophyletic origin in the Late Paleozoic (300myr) of the three living amphibian orders and a sister-group relationship between Anura and Caudata (the Batrachia hypothesis) (Laurin and Reisz, 1997; Trueb and Cloutier, 1991). However, Jarvik (1980) has suggested that Lissamphibia is not a natural group and hypothesized a polyphyletic origin of the three living orders. Carroll and coworkers (Carroll, 1988; Carroll and Holmes, 1980) have advocated a close phylogenetic relationship between salamanders and caecilians. According to these authors, frogs would have evolved from temnospondyl amphibians (Apsidospondyli) whereas salamanders and caecilians originated from different lineages of microsaurs (Lepospondyli). The affinity between salamanders and caecilians is also supported by several morphological (Bolt, 1991) and molecular (Feller and Hedges, 1998; Hedges et al., 1990) phylogenetic

Based on the analysis of structural features and conserved elements, 27 novel snoRNA genes have been identified from rice. All of them belong to the C/D box-containing snoRNA family except for one that belongs to the H/ACA box type. The newly found genes fall into six clusters that comprise at least three snoRNA genes, and in one case as many as nine genes. Interestingly, four of the six clusters are located within the largest intron of a protein coding gene. The majority of intronic snoRNA gene clusters are simply formed by multiple copies of the same species of snoRNA gene that possess the identical functional elements. This implies a possible mechanism of duplication for the origin of repeating snoRNA coding regions in one intron. However, a few intronic snoRNA gene clusters consisting of different snoRNAs species were also observed. Polycistronic precursors from two independently transcribed clusters were demonstrated by RT-PCR and individual snoRNAs processed from the polycistronic precursors were positively determined by reverse transcription assay. Analyses of the intergenic spacers in the clusters showed that, in addition to a very high AT content, the processing signals in rice snoRNA polycistronic transcripts might be different from those of yeast. Our results demonstrate that, in both plants and mammals, numerous snoRNAs can be produced simultaneously from an mRNA precursor of a host gene despite the different arrangements. The intronic snoRNA gene cluster is a novel gene organization, which is so far unique to plants. The conservation of intronic snoRNA gene clusters in plants was further demonstrated by the study of a similar snoRNA gene organization in the first intron of a Hsp70 gene from wild rice and Zizania caduciflora.

Through a computer search of DNA databases, we have identified the homologs of the mgU6-47 snoRNA gene from the yeast Schizosaccharomyces pombe, the fly Drosophila melanogaster and human. The three box C/D-containing snoRNA genes showed no significant similarity in their sequences except for an 11 nt long complementarity to U6 snRNA, suggesting that the mechanism of snoRNA guided snRNA methylation is conserved from mammals to yeast. The corresponding snoRNAs have been positively detected by reverse transcription and northern blotting. Taking advantage of the fission yeast system, we have disrupted the yeast mgU6-47 gene and demonstrated that it is absolutely required for site-specific 2'-O-methylation of U6 at position A41. No growth differences between mgU6-47 gene-disrupted and wild-type cells were observed, suggesting that the mgU6-47 gene, as for most rRNA methylation guides, is dispensable in yeast. Nevertheless, it was revealed by temperature shift assay that abolition of A41 methylation in yeast U6 snRNA might cause a small decrease in mRNA splicing efficiency. The timing of S.pombe U6 pre-RNA transport in the nucleus for splicing and methylation was also analyzed and is described.

Through a computer search of the genome of the yeast Saccharomyces cerevisiae, the coding sequences of seven different box C/D antisense small nucleolar RNAs (snoRNAs) with the structural hallmarks of guides for rRNA ribose methylation have been detected clustered over a 1.4-kb tract in an inter-open reading frame region of chromosome XIII. The corresponding snoRNAs have been positively identified in yeast cells. Disruption of the nonessential snoRNA gene cluster specifically suppressed the seven cognate rRNA ribose methylations but did not result in any growth delay under the conditions of yeast culture tested. The seven snoRNAs are processed from a common polycistronic transcript synthesized from an independent promoter, similar to some plant snoRNAs but in marked contrast with their vertebrate functional homologues processed from pre-mRNA introns containing a single snoRNA. Processing of the polycistronic precursor requires nucleases also involved in rRNA processing, i.e., Rnt1p and Rat1p. After disruption of the RNT1 gene, the yeast ortholog of bacterial RNase III, production of the seven mature snoRNAs was abolished, while the polycistronic snoRNA precursor accumulated. In cells lacking functional Rat1p, an exonuclease involved in the processing of both pre-rRNA and intron-encoded snoRNAs, several processing intermediates of the polycistronic precursor accumulated. This allowed for the mapping in the precursor of the presumptive Rnt1p endonucleolytic cuts which provide entry sites for subsequent exonucleolytic trimming of the pre-snoRNAs. In line with known properties of double-stranded RNA-specific RNase III, pairs of Rnt1p cuts map next to each other on opposite strands of long double-helical stems in the secondary structure predicted for the polycistronic snoRNA precursor.

A growing number of small nucleolar RNAs (snoRNAs) are intron-encoded, contain the characteristic box C (UGAUGA) and box D (CUGA) motifs and exhibit long complementarities to conserved sequences in mature rRNAs. We have identified nine additional members of this family, U32 to U40. All but one are encoded in introns of ribosomal protein genes in vertebrates: U32 to U35 in rpL13a, U36 in rpL7a and U38 to U40 in rpS8. By contrast, U37 is encoded in elongation factor 2 gene. Interestingly, U32 and U36 each contain two complementarities (one to 18 S and the other to 28 S rRNA). U32 to U40 are fibrillarin-associated, devoid of a 5'-trimethyl-cap and display an exclusively nucleolar localization. They are all metabolically stable and roughly as abundant as previously reported members of this family. Characterization of their homologs in distant species shows that their 10 to 14 nt long rRNA complementarities are conserved. A clue on the function of this snoRNA family is provided by the comparative analysis of the largely expanded collection of their conserved duplexes with rRNA. Not only does each duplex span at least one site of 2'-O-ribose methylation in the rRNA but the modification site is always at the same position in the duplex, paired to the fifth nucleotide upstream from a box D motif in the snoRNA. Consistent with the notion that each snoRNA of this family guides one particular methylation along the rRNA sequence, we have detected several pairs of snoRNAs with overlapping complementarities to rRNA tracts with vicinal sites of ribose methylations. In each case, the two overlapping complementarities are shifted from each other by a distance equal to the spacing between the methylated sites which are thus found at the same position within each of the mutually exclusive duplexes. Finally, we have also identified, within three previously known snoRNAs, novel antisense elements able to form a canonical duplex around ribose-methylated sites in rRNA, which further supports the conclusion that the duplex structure provides the 2'-O-methyltransferase with the appropriate site-specificity on the substrate.

Recently, several new snoRNAs encoded in introns of genes coding for ribosomal, ribosome-associated, or nucleolar proteins have been discovered. We are presently studying four of these intronic snoRNAs. Three of them, U20, U21, and U24, are closely related to each other on a structural basis. They are included in genes encoding nucleolin and ribosomal proteins L5 and L7a, respectively, in warm-blooded vertebrates. These three metabolically stable snoRNAs interact with nucleolar protein fibrillarin. In addition, they display common features that make them strikingly related to snoRNA U14. U14 contains two tracts of complementarity to 18S rRNA, which are required for the production of 18S rRNA. U20 displays a 21 nucleotide (nt) long complementarity to 18S rRNA. U21 contains a 13 nt complementarity to an invariant sequence in eukaryotic 28S rRNA. U24 has two separate 12 nt long complementarities to a highly conserved tract of 28S rRNA. Phylogenetic evidences support the fundamental importance of the pairings of these three snoRNAs to pre-rRNA, which could be involved in a control of pre-rRNA folding during preribosome assembly. By transfection of mouse cells, we have also analyzed the processing of U20 and found that the -cis acting signals for its processing from intronic RNA are restricted to the mature snoRNA sequence. Finally, we have documented changes of host genes for these three intronic snoRNAs during the evolution of eukaryotes.Key words: snoRNA, pre-rRNA, folding, genes, introns.

Following computer searches of sequence banks, we have positively identified a novel intronic snoRNA, U24, encoded in the ribosomal protein L7a gene in humans and chicken. Like previously reported intronic snoRNAs, U24 is devoid of a 5'-trimethyl-cap. U24 is immunoprecipitated by an antifibrillarin antibody and displays an exclusively nucleolar localization by fluorescence microscopy after in situ hybridization with antisense oligonucleotides. In vertebrates, U24 is a 76 nt long conserved RNA which is metabolically stable, present at approximately 14,000 molecules per human HeLa cell. U24 exhibits a 5'-3' terminal stem-box C-box D structure, typical for several snoRNAs, and contains two 12 nt long conserved sequences complementary to 28S rRNA. It is, therefore, strikingly related to U14, U20 and U21 snoRNAs which also possess long sequences complementary to conserved sequences of mature 18S or 28S rRNAs. In 28S rRNA the two tracts complementary to U24 are adjacent to each other, they involve several methylated nucleotides and are surprisingly close, within the rRNA secondary structure, to complementarities to snoRNAs U18 and U21. Identification of the yeast Saccharomyces cerevisiae U24 gene directly confirms the outstanding conservation of the complementarity to 28S rRNA during evolution, suggesting a key role of U24 pairing to pre-rRNA during ribosome biogenesis, possible in the control of pre-rRNA folding. Yeast S.cerevisiae U24 is also intron-encoded but not in the same host-gene as in humans or chicken.

U3 small nucleolar RNA, which participates in eukaryotic rRNA processing, is encoded by a small multigene family in mammals. In humans, the four to six gene copies code for an identical U3 RNA molecule; rodents, however, have two variant forms of these genes, U3A and U3B. We show that all U3 genes in humans map to a single chromosomal locus, 17p12→p11, which corresponds exactly to the region of mouse Chromosome 11 where the four U3B genes are clustered. By contrast, in mouse the unique U3A gene copy is not linked to the U3B gene cluster but maps to another chromosome (the B2-B4 region of Chromosome 10).

In humans, cellular 28S rRNA displays a sequence dimorphism within an evolutionarily conserved motif, with the presence, at position +60, of either a A (like the metazoan consensus) or a G. The relative abundance of the two forms of variant genes in the genome exhibit large differences among individuals. The two variant forms are generally represented in cellular 28S rRNA in proportion of their relative abundance in the genome, at least for leucocytes. However, in some cases, one form of variant may be markedly underexpressed as compared to the other. Thus, in HeLa cells, A-form genes contribute to only 1% of the cellular content in mature 28S rRNA although amounting to 15% of the ribosomal genes. The differential expression seems to result from different transcriptional activities rather than from differences in pre-rRNA processing efficiency or in stabilities of mature rRNAs. G-form ribosomal genes were not detected in other mammals, including chimpanzee, which suggests that the fixation of this variant type is a rather recent event in primate evolution.

A U3 RNA variant has been identified in mouse, the abundance of which relative to the previously characterized major form (U3B) appears to vary to a large extent depending upon the cell origin. Its partial sequence analysis shows that it is clearly related to the U3A form previously described in rat. Sequence comparisons suggest that the separation of the two forms of U3 genes now found in rat and mouse represent a relatively ancient event in rodent evolution. While mouse U3B RNA is encoded by four clustered genes, the U3A variant is encoded by a unique gene. Both mouse U3 RNAs differ substantially in primary structure (more than 10% divergence). Although rodent U3 RNAs exhibit a largely similar secondary structure, a specific difference between the A and B form can nevertheless be observed.

Using a rapid rRNA sequencing technique, we have determined the sequence of the 400 nucleotides located at the 5' end of the large subunit rRNA molecule from eight species of unicellular eukaryotes (protists). This region contains a pair of conservative domains well-suited for long-range phylogenetic evaluations among eukaryotes, due both to their substantia, length and to their intrinsic rate of sequence variation during evolution. It also comprises a central more rapidly evolving portion, which allows for a fine tuning of distance evaluation between closely related species. Molecular distances were computed between the aligned nucleotides of all presently available protist sequences and were used to derive a tentative dendrogram. Within the limitations inherent to this approach, a number of interesting observations emerge: The various protist groups appear to have separated very early from each other. The most deeply divergent protists belong to a number of orders of flagellates (mastigotes), suggesting a very ancient origin for organelles containing a 9 + 2 microtubular arrangement. Ciliates emerged late among eukaryotes, suggesting that their peculiar genetic code was derived secondarily. Moreover, a dinoflagellate clusters with ciliates, thus making it likely that the unusual features of nuclear organization and mitosis of this group are not primitive but derived characters. Finally, within groups, taxonomic and evolutionary inferences appear to be feasible using this portion of the rRNA.

We have developed a combined approach for probing native structures in large RNAs. In the first method, after digestion with a structure specific nuclease, accessible sites are mapped at sequence resolution along the entire RNA molecule which is used as a template for the reverse transcriptase elongation of a 5' end labelled selected primer (coding strand of a small restriction fragment of the cloned gene). This method circumvents any prior end-labelling of RNA, a technique with major limitations for large RNAs. In the second approach, a rapid "heterologous" sequencing can be easily applied to definite domains of an RNA molecule in a variety of species (or individuals), without additional DNA cloning nor end-labelling of RNA. By taking advantage of the presence of evolutionary conserved tracts within an RNA sequence, it allows a rapid analysis of RNA folding patterns in terms of phylogenetic comparisons : when located within such a conserved tract, selected restriction fragments from a cloned gene can be used as heterologous primers for sequencing the upstream divergent region in RNAs of other species by currently available technology, i.e. reverse transcriptase elongation in the presence of chain terminator dideoxynucleotides.