|
|
|
|
Published online before print
November 21, 2007, 10.1101/gr.6539108 Genome Res. 18:104-112, 2008 ©2008 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/08 $5.00
Letter Xenopus microRNA genes are predominantly located within introns and are differentially expressed in adult frog tissues via post-transcriptional regulationDepartment of Molecular and Structural Biochemistry, North Carolina State University, Raleigh, North Carolina 27695-7622, USA
The amphibian Xenopus provides a model organism for investigating microRNA expression during vertebrate embryogenesis and development. Searching available Xenopus genome databases using known human pre-miRNAs as query sequences, more than 300 genes encoding 142 Xenopus tropicalis miRNAs were identified. Analysis of Xenopus tropicalis miRNA genes revealed a predominate positioning within introns of protein-coding and nonprotein-coding RNA Pol II-transcribed genes. MiRNA genes were also located in pre-mRNA exons and positioned intergenically between known protein-coding genes. Many miRNA species were found in multiple locations and in more than one genomic context. MiRNA genes were also clustered throughout the genome, indicating the potential for the cotranscription and coordinate expression of miRNAs located in a given cluster. Northern blot analysis confirmed the expression of many identified miRNAs in both X. tropicalis and X. laevis. Comparison of X. tropicalis and X. laevis blots revealed comparable expression profiles, although several miRNAs exhibited species-specific expression in different tissues. More detailed analysis revealed that for some miRNAs, the tissue-specific expression profile of the pri-miRNA precursor was distinctly different from that of the mature miRNA profile. Differential miRNA precursor processing in both the nucleus and cytoplasm was implicated in the observed tissue-specific differences. These observations indicated that post-transcriptional processing plays an important role in regulating miRNA expression in the amphibian Xenopus.
MicroRNAs are small noncoding RNAs of between 19 and 23 nucleotides in length that play fundamental roles in regulating gene expression in both plants and animals (Carrington and Ambros 2003 30% of human protein-coding genes are post-transcriptionally regulated by miRNAs (Lewis et al. 2005
MiRNAs are well conserved in evolution, particularly in animals, suggesting that they play conserved functions in different organisms (Pasquinelli et al. 2000
Analysis of miRNA genes has revealed a variety of genomic organizations (Lai et al. 2003
A number of approaches have been used to define miRNAs in various organisms. The cloning and sequencing of total small RNAs of appropriate size from isolated tissues or whole organisms has led to the identification of many miRNA species (for review, see Berezikov et al. 2006
The objective of this study was to computationally identify Xenopus miRNA genes, characterize their genomic organization, and then examine their expression in adult frog tissues. We present here 142 Xenopus miRNAs homologs. Like other vertebrates,
Computational identification of Xenopus microRNAs Human pre-miRNA sequences were obtained from the microRNA database (miRBase: Giffiths-Jones et al. 2006
A total of 142 putative X. tropicalis miRNAs were identified in our computational search (Table 1; Supplemental Table 1). Included in these 142 species are 28 of the 29 homologs previously shown to be differentially expressed in X. laevis during embryogenesis (Watanabe et al. 2005
Each of the identified Xenopus miRNAs is able to fold into the hairpin–loop structure characteristic of a folded pre-miRNA. Folding analysis was first accomplished using Mfold (Zuker 2003
Xenopus miRNA-coding sequences are predominantly located within introns of RNA Pol II-transcribed host genes
The X. tropicalis intronic miRNAs can be grouped into three types. Some miRNAs are located within introns of known protein-coding genes. Others are encoded in putative introns, predicted to be within hypothetical protein-coding sequences (designated in Supplemental Table 2 by scaffold or protein ID numbers). MicroRNAs were also identified within designated introns of "non-spliced EST" sequences (Mott 1997
Xenopus miRNA genes are clustered in the genome
A particularly striking example of miRNA clustering was the Xenopus xtr-mir-427 coding sequences. We have identified 80 xtr-mir-427 genes in the X. tropicalis genome (Supplemental Table 2), many of them clustered on specific DNA scaffolds as shown in Figure 2. The length of DNA scaffold containing xtr-mir-427 miRNAs varied greatly from only several hundred to more than 100,000 bp. The close proximity of some clustered miRNAs suggests cotranscription in a single pri-miRNA transcript, whereas distantly positioned miRNAs may suggest otherwise. Strikingly, analysis of available X. tropicalis EST data indicates that all but two indicated clusters could well be transcribed as a single, common precursor transcript (designated as +ESTs). Even more striking was the observation that some of the miRNAs within a cluster were absent from the EST (open triangles), indicating differential processing of the pri-miRNA transcript.
Expression of X. tropicalis and X. laevis miRNAs Total RNA was prepared from X. tropicalis and X. laevis adult tissues and miRNA expression assessed using Northern blot analysis. Of the 63 miRNAs assessed, 47 species or 77% (designated in Supplemental Table 1) were observed in one or more of the tissue-specific RNA preparations (38 of these are presented in Fig. 3). Those miRNAs tested but not detected may be expressed in other tissues not examined, expressed during specific stages of embryogenesis, or expressed at such low levels as to make detection by Northern blot analysis difficult. Most of the observed miRNAs were expressed in almost all of the adult tissues assayed. The lack of miRNA expression in oocytes is consistent with miRNA roles in regulating gene expression during embryogenesis and organism development.
The expression patterns were comparable for most of the miRNAs in both X. tropicalis and X. laevis tissues. However, apparent quantitative differences in the expression patterns of some miRNAs were observed between the two different Xenopus species. These miRNAs included xtr-mir-98, xtr-mir-133a, xtr-mir-142-5p, and xtr-mir-143 (Fig. 3). Differences in expression patterns could reflect species-specific differences in transcription or differential pri-mRNA processing. Some clustered miRNAs exhibited comparable expression patterns (e.g., xtr-mir-99/xtr-let-7c cluster), consistent with cotranscription in a common pri-miRNA as previously suggested (Lagos-Quintana et al. 2001
Evidence for the regulated, tissue-specific expression of X. laevis miRNAs
Both xtr-mir-98 and xtr-let-7f are clustered on a 427-nt DNA fragment (Xentr3_750) and can reasonably be assumed to be cotranscribed as a common pri-miRNA precursor transcript. The size of the pri-miRNA transcript observed is consistent with the cotranscription of these three miRNAs. Notably, while the patterns of xtr-let-7f and xtr-mir-98 pri-miRNA expression are identical, the patterns of tissue-specific expression for the two mature miRNAs are clearly different (Fig. 4C). Particularly prominent is the presence of xtr-let-7f miRNA in heart tissue, but the absence of xtr-mir-98 in the same tissue. This suggested the differential processing of the pri-miRNA precursor or the pre-miRNA intermediate to produce xtr-let-7f but not xtr-mir-98. Finally, both mature xtr-mir-215 and xtr-mir-220b were not observed in initial Northern blot analysis. However, additional examination revealed the presence of pri-miRNA primary transcripts for both miRNAs in all adult tissues and the presence of xtr-mir-200b pre-miRNA in some adult tissues (Fig. 4D). These observations suggested the lack of xtr-mir-215 pri-miRNA primary transcript processing to produce either pre-miRNAs or mature miRNAs in all tissues and the differential processing of xtr-pre-mir-200b in some tissues, along with the possible disruption of pre-miRNA transport and/or cytoplasmic Dicer processing in all tissues.
A computational approach using human pre-miRNA query sequences has identified more than 300 miRNA genes encoding some 142 Xenopus miRNA species. Xenopus miRNA genes are predominantly located within introns and highly clustered in the genome. Using Northern blot analysis as well as the EST databases, a substantial proportion of the identified Xenopus miRNAs have been verified as expressed RNA transcripts. The shared expression patterns of almost all X. tropicalis and X. laevis miRNAs in the various adult tissues are consistent with highly conserved miRNA function in vertebrates. However, a few probed miRNAs exhibited tissue-specific differences in these two Xenopus species, suggesting perhaps species-specific miRNA expression. Most striking however, was the observation of post-transcriptional regulation of miRNA expression at both nuclear and cytoplasmic processing steps.
Xenopus miRNAs exhibit several genomic contexts, but their predominant positioning within pre-mRNA introns is similar to other vertebrates, including human, where it has been estimated that 70% of the miRNAs are intronic (Rodriguez et al. 2004
Like other vertebrates, miRNA clustering is a prominent feature in the Xenopus genome. Approximately half of the identified Xenopus miRNAs are found in clusters in varying combinations. The clustering of Xenopus miRNA genes could also have important implications for the coordinated, tissue-specific expression of clustered miRNAs during development. Certainly, pri-miRNAs containing clustered miRNAs have been reported (Lau et al. 2001
Of particular importance will be defining promoters for pri-miRNA transcription of intergenic, intronic, and clustered miRNAs. Analysis of miRNA synthesis has revealed that the majority of pri-miRNAs are transcribed by Pol II polymerase and possess 5' caps and 3' poly (A) tails (Cai et al. 2004
Further examination of transcribed Xenopus miRNA sequences in the EST databases may suggest an even more complex picture of miRNA transcription and processing. Those miRNAs occurring in the EST databases are indicated in Supplemental Tables 2 and 3. These previously designated "spliced" miRNAs are flanked by upstream and downstream transcript corresponding to the genomic sequences flanking the pre-miRNA sequence. However, far more intriguing is the occurrence of so-called "non-spliced" RNA transcripts (Mott 1997
Finally, analysis of Xenopus miRNA expression has revealed unanticipated expression profiles for some Xenopus pri-miRNAs and their processed, mature miRNAs. Differences in tissue-specific expression between the two Xenopus species were surprising, yet other laboratories have reported similar observations (Ason et al. 2006
Computational identification of Xenopus microRNAs Human pre-miRNA sequences were retrieved from miRbase::sequences (http://microrna.sanger.ac.uk/sequences/) and used as query sequences against available X. tropicalis genomic databases (Xentr3-Oct. 2004 assembly and the Xentr4-Aug. 2005 assembly with most of our reported data obtained from the 2004 assembly and designated as Xentr3). This computational approach is summarized in Supplemental Figure 1 and consisted of homology searching, alignment and annotation, and pre-miRNA secondary structural analysis. The X. tropicalis genome database is presently under construction (Klein et al. 2002
Candidate Xenopus miRNAs were further analyzed for their folded structures with specific pre-miRNA secondary structural features being used as filters in the verification process. Identified Xenopus sequences (60–120 nt long) were folded using the Mfold program (http://www.bioinfo.rpi.edu/applications/mfold/old/rna/) and manually adjusted according to specific criteria based upon recent analysis and definition of conserved or preferred pre-miRNA structures. These features included (1) the length of the pre-miRNA core region between 60 and 80 nt (Ambros et al. 2003
Northern blot analysis
We gratefully acknowledge support from the College of Agriculture and Life Sciences at North Carolina State University, Raleigh, NC.
1 Corresponding author.
E-mail stu_maxwell{at}ncsu.edu; fax (919) 515-2047. [Supplemental material is available online at www.genome.org.] Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6539108
Altuvia, Y., Landgraf, P., Lithwick, G., Elefant, N., Pfeffer, S., Aravin, A., Brownstein, M.J., Tuschl, T., and Margalit, H. 2005. Clustering and conservation patterns of human microRNAs. Nucl. Acids Res. 33: 2697–2706. doi: 10.1093/nar/gki567. Ambros, V., Bartel, B., Bartel, D.P., Burge, C., Carrington, J., Chen, X., Dreyfuss, G., Eddy, S., Griffiths-Jones, S., Marshall, M., et al. 2003. A uniform system for microRNA annotation. RNA 9: 277–279. Ason, B., Darnell, D., Wittbrodt, B., Berezikov, E., Kloosterman, W., Wittbrodt, J., Antin, P., and Plasterk, R. 2006. Differences in vertebrate microRNA expression. Proc. Natl. Acad. Sci. 103: 14385–14389. Baskerville, S. and Bartel, D. 2005. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA 11: 241–247. Bartel, D. 2004. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 116: 281–297.[CrossRef][Medline] Bentwich, I. 2005. Prediction and validation of microRNAs and their targets. FEBS Lett. 579: 5904–5910.[CrossRef][Medline] Berezikov, E., Cuppen, E., and Plasterk, R. 2006. Approaches to microRNA discovery. Nat. Genet. 38: S2–S7 (Suppl) doi: 10.1038/ng1794.[CrossRef][Medline] Bilen, J., Liu, N., and Bonini, N. 2006. A new role for microRNA pathways: Modulation of degeneration induced by pathogenic human disease proteins. Cell Cycle 5: 2835–2838.[Medline] Bonnet, E., Wuyts, J., Rouze, P., and de Van Peer, Y. 2004. Detection of 91 potential conserved plant microRNAs in Arabidopsis thaliana and Oryza sativa identifies important target genes. Proc. Natl. Acad. Sci. 101: 11511–11516. Borchert, G., Lanier, W., and Davidson, B. 2006. RNA polymerase III transcribes human microRNAs. Nat. Struct. Mol. Biol. 13: 1097–1101.[CrossRef][Medline] Cai, X., Hagedorn, C., and Cullen, B. 2004. Human microRNAs are processed from capped, polyadenylated transcripts that can also function as mRNAs. RNA 10: 1957–1966. Carrington, J. and Ambros, V. 2003. Role of microRNAs in plant and animal development. Science 301: 336–338. Chen, K. and Rajewsky, N. 2007. The evolution of gene regulation by transcription factors and microRNAs. Nat. Rev. Genet. 8: 93–103.[Medline] Griffiths-Jones, S., Grocock, R., van Dongen, S., Bateman, A., and Enright, A. 2006. miRBase: microRNA sequences, targets and gene nomenclature. Nucl. Acids Res. 34: D140–D144. doi: 10.1093/nar/gkj112. He, L., Thomson, J., Hemann, M., Hernando-Monge, E., Mu, D., Goodson, S., Powers, S., Cordon-Cardo, C., Lowe, S., Hannon, G., et al. 2005. A microRNA polycistron as a potential human oncogene. Nature 435: 828–833.[CrossRef][Medline] John, B., Enright, A., Aravin, A., Tuschl, T., Sander, C., and Marks, D. 2005. Erratum. PLoS Biol 3: e363. doi: 10.1371/journal.pbio.0020363.[CrossRef][Medline] Kim, N. and Nam, J. 2006. Genomics of microRNA. Trends Genet. 22: 165–173.[CrossRef][Medline] Klein, S., Strausberg, R., Wagner, L., Pontius, J., Clifton, S., and Richardson, P. 2002. Genetic and genomic tools for Xenopus research: The NIH Xenopus initiative. Dev. Dyn. 225: 384–391.[CrossRef][Medline] Kloosterman, W. and Plasterk, R. 2006. The diverse functions of microRNAs in animal development and disease. Dev. Cell 11: 441–450.[CrossRef][Medline] Lagos-Quintana, M., Rauhut, R., Lendeckel, W., and Tuschl, T. 2001. Identification of novel genes coding for small expressed RNAs. Science 294: 853–858. Lai, E., Tomancak, P., Williams, R., and Rubin, G. 2003. Computational identification of Drosophila microRNA genes. Genome Biol. 4: R42. doi: 10.1186/gb-2003-4-7-r42.[CrossRef][Medline] Lau, N., Lim, L., Weinstein, E., and Bartel, D. 2001. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science 294: 858–862. Lee, Y., Kim, M., Han, J., Yeom, K., Lee, S., Baek, S., and Kim, V. 2004. MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 23: 4051–4060.[CrossRef][Medline] Lewis, B., Burge, C., and Bartel, D. 2005. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120: 15–20.[CrossRef][Medline] Li, S., Pan, C., and Lin, W. 2006. Bioinformatic discovery of microRNA precursors from human ESTs and introns. BMC Genomics 7: 164.[CrossRef][Medline] Lim, L., Lau, N., Weinstein, E., Abdelhakim, A., Yekta, S., Rhoades, M., Burge, C., and Bartel, D. 2003a. The microRNAs of Caenorhabditis elegans. Genes & Dev. 17: 991–1008. Lim, L., Glasner, M., Yekta, S., Burge, C., and Bartel, D. 2003b. Vertebrate microRNA genes. Science 299: 1540. Liu, S., Sun, Y.H., Shi, R., Clark, C., Li, L., and Chiang, V.L. 2005. Novel and mechanical stress-responsive microRNAs in Populus trichocarpa that are absent from Arabidopsis. Plant Cell 17: 2186–2203. Lund, E., Guttinger, S., Calado, A., Dahlberg, J., and Kutay, U. 2004. Nuclear export of microRNA precursors. Science 303: 95–98. Mott, R. 1997. EST_GENOME: A program to align spliced DNA sequences to unspliced genomic DNA. Comput. Appl. Biosci. 13: 477–488. Obernosterer, G., Leuschner, P., Alenius, M., and Martinez, J. 2006. Post-transcriptional regulation of microRNA expression. RNA 12: 1161–1167. Pasquinelli, A., Reinhart, B., Slack, F., Martindale, M., Kuroda, M., Maller, B., Hayward, D., Ball, E., Degnan, B., Muller, P., et al. 2000. Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature 408: 86–89.[CrossRef][Medline] Rajewsky, N. 2006. MicroRNA target predictions in animals. Nat. Genet. 38: S8–S13 (Suppl) doi: 10.1038/ng1798.[CrossRef][Medline] Rodriguez, A., Griffiths-Jones, S., Ashurst, J., and Bradley, A. 2004. Identification of mammalian microRNA host genes and transcription units. Genome Res. 14: 1902–1910. Samols, M., Hu, J., Skalsky, R., and Renne, R. 2005. Cloning and identification of a microRNA cluster within the latency-associated region of Kaposi's sarcoma-associated herpesvirus. J. Virol. 79: 9301–9315. Seitz, H. and Zamore, P. 2006. Rethinking the microprocessor. Cell 125: 827–829.[CrossRef][Medline] Smalheiser, N. 2003. EST analyses predict the existence of a population of chimeric microRNA precursor-mRNA transcripts expressed in normal human and mouse tissues. Genome Biol. 4: 403.[CrossRef][Medline] Sood, P., Krek, A., Zavolan, M., Macino, G., and Rajewsky, N. 2006. Cell-type-specific signatures of microRNAs on target mRNA expression. Proc. Natl. Acad. Sci. 103: 2746–2751. Tang, G. 2005. siRNA and miRNA: An insight into RISCs. Trends Biochem. Sci. 30: 106–114.[CrossRef][Medline] Thomson, J., Newman, M., Parker, J., Morin-Kensicki, E., Wright, T., and Hammond, S. 2007. Extensive post-transcriptional regulation of microRNAs and its implications for cancer. Genes & Dev. 20: 2202–2207.[CrossRef] Valencia-Sanchez, M., Liu, J., Hannon, G., and Parker, R. 2006. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes & Dev. 20: 515–524. van Rooij, E., Sutherland, L., Liu, N., Williams, A., McAnally, J., Gerard, R., Richardson, J., and Olson, E. 2006. A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc. Natl. Acad. Sci. 103: 18255–18260. Watanabe, T., Takeda, A., Mise, K., Okuno, T., Suzuki, T., Minami, N., and Imai, H. 2005. Stage-specific expression of microRNAs during Xenopus development. FEBS Lett. 579: 318–324.[CrossRef][Medline] Weber, M. 2005. New human and mouse microRNA genes found by homology search. FEBS J. 272: 59–73.[CrossRef][Medline] Ying, S. and Lin, S. 2006. Current perspectives in intronic microRNAs (miRNAs). J. Biomed. Sci. 13: 5–15.[CrossRef][Medline] Zeng, Y., Yi, R., and Cullen, B. 2005. Recognition and cleavage of primary microRNA precursors by the nuclear processing enzyme Drosha. EMBO J. 12: 138–148. Zuker, M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucl. Acids Res. 31: 3406–3415. doi: 10.1093/nar/gkg595.
Received March 23, 2007; accepted in revised format October 10, 2007.
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||