Computer identification of snoRNA genes using a Mammalian Orthologous
http://www.100md.com
《核酸研究医学期刊》
Department of Medicine, Program in Bioinformatics and Proteomics/Genomics, Medical University of Ohio Toledo, OH 43614, USA 1Department of Biological Sciences, Life Sciences, Bowling Green State University Bowling Green, OH 43403, USA 2Department of Bioinformatics, Institute of Molecular Genetics, RAS Moscow 123182, Russia
*To whom correspondence should be addressed. Tel: +1 419 383 5270; Fax: +1 419 383 3102; Email: afedorov@meduohio.edu
ABSTRACT
Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.
INTRODUCTION
Small nucleolar RNA (snoRNA) is a major component of small nucleolar ribonucleoprotein (snoRNP) particles that are located inside the nucleolus of eukaryotes and participate in post-transcriptional chemical modification or processing of different RNAs, including ribosomal RNA (rRNA) and small nuclear RNA (snRNA) . SnoRNAs are ancient genes, as they are widespread through the entire domain of Eukaryota, as well as Archaea. There are two types of these RNAs, called C/D-box and H/ACA-box snoRNAs, and they are characterized by distinct 3D structures and conserved sequence elements. The C/D-box RNA is associated with 2'-O-ribose methylation, and H/ACA-box RNA with pseudouridylation of substrate RNAs. The direct role of snoRNAs is determining the site for chemical modification, via complementary pairing of its specific sequence with the segment of substrate RNA undergoing modification. The major catalytic activity does not belong to snoRNA, but rather to a fibrillarin, a protein component of the snoRNP complex (1,6). Besides pseudouridylation and 2'-O-ribose methylation, some snoRNAs perform other documented and putative functions, such as (i) cleavage of the substrate precursor in rRNA (2), (ii) facilitation of rRNA folding (7), (iii) regulation of alternative splicing (8) and (iv) possible unknown functions carried out by so-called ‘orphan’ snoRNAs (2).
The first well-known computer program for genomic prediction of C/D-box snoRNA, named SNOSCAN, was created by Lowe and Eddy (9). SNOSCAN was used by several laboratories (10) to study compact genomes of eukaryotes such as yeast, Drosophila and Arabidopsis. SNOSCAN can predict classical C/D-box snoRNAs that guide modification of rRNA or snRNA molecules, and requires that these rRNA and snRNA sequences be available for the program before its invocation. In 2003, Vitali et al. (7) used this program on a restricted set of human and mouse ribosomal protein gene sequences, and found several novel C/D snoRNAs inside their introns. Recently, two computational approaches were developed to predict H/ACA snoRNAs from genomic sequences (11,12). Again, these programs can predict snoRNAs that are involved in chemical modifications of known rRNA and snRNA sequences, and they are suitable for searching compact eukaryotic genomes. However, it is problematic to apply the existing computational algorithms for studying vertebrate genomes in their entirety, because they are dozens of times larger than the genomes of yeast, insects or Arabidopsis. In the case of vertebrates, numerous false-positive signals may arise due to the computer processing of very large non-coding genomic sequences and thus computer prediction of novel genes may not be as efficient.
We have taken advantage of comparative genomics and present here a new algorithm for computational identification of mammalian C/D-box snoRNA genes with high efficiency. We exploit two well-known characteristics of mammalian genomic structure: (i) all known snoRNA genes of mammals are located inside introns (3) and (ii) the exon–intron gene structure of mammals is highly conserved. In fact, no well-documented case of intron gain exists between fully sequenced mammalian genomes, and there are only a few intron losses reported (13,14). Therefore, introns are ‘fossilized’ in mammalian genes, and the non-coding RNAs (ncRNAs, also known as small non-messenger RNAs, non-protein coding RNAs and untranslated RNAs) located within them could well be fixed inside the same intron since the origin of this taxon. For this reason, we created a database of orthologous introns comprising human, mouse and rat sequences. We define ‘orthologous introns’ as introns from orthologous genes that have the same position relative to the coding sequence. Hence, orthologous introns should have descended from the corresponding intronic sequence of the last common ancestor for the taxon.
The evolutionary divergence of primate and rodent lineages occurred 70–90 million years ago, and during this period non-functional DNA segments representing a major portion of human and rodent orthologous introns lost almost all sequence similarity. For example, Figure 1 demonstrates an alignment of the third introns of the ribosomal protein S3a gene from human and mouse genomes that are orthologous to each other and are 1708 and 1990 nt long, respectively. It is known that these two introns contain U73b snoRNA sequences of human and mouse (NG_000961 and Z83331 , respectively). A BLAST-2 alignment with default parameters reveals 88% sequence identity of a 58 nt long fragment of these introns containing snoRNA, and no sequence similarity outside this short functional intronic fragment except one 37 nt long segment (Figure 1). This second conserved segment could also represent a functional sequence not yet known. Therefore, when conserved snoRNA structures are identified within orthologous intron triplets from all three species (human, mouse and rat), it strongly suggests that the conserved sequence is functional. Here, we present our Mammalian Orthologous Intron Database (MOID) for public usage and two programs that search for conserved C/D-box snoRNA motifs inside the entire MOID and characterize putative ASEs of snoRNAs based on comparative genomics of mammalian species.
Figure 1 Results of an online Blast2 alignment of the 1708 nt long third intron of the human ribosomal protein S3a gene, and its mouse ortholog (1990 nt long third intron of the mouse rps3a gene). (A) Dot plot figure of the intron comparison obtained by bl2seq online program (http://www.ncbi.nlm.nih.gov/blast/bl2seq/wblast2.cgi) with default parameters. (B) Alignment of human and mouse introns. The sequences of snoRNA C- and D-boxes inside these introns are boxed.
Mammalian introns bear all types of snoRNA genes (3), about three-fourth of microRNA genes (15), and also many conserved regulatory motifs (16). Our programs can be easily adapted to search for all kinds of functional sequences using MOID. We are further developing this approach by extending our Orthologous Intron Database with other mammalian and vertebrate species.
MATERIALS AND METHODS
Intron databases
Using the previously published program package for the generation of the Exon–Intron Database (EID) (17), we created three species-specific EIDs containing all known genes of mammals with entirely sequenced genomes. These EID for human, mouse and rat genomes and their documentation are freely available from our website (http://www.meduohio.edu/bioinfo/eid/index.html). We then applied the CIP.pl program (18) for the comparison of intron positions in all orthologous genes from human, mouse and rat genomes. Based on this comparison, we generated an MOID, which is also publicly available from the same web page (documentation provided). This first release of our MOID contains 100 000 orthologous introns from human, mouse and rat genomes. The next release of MOID, with a complete set of orthologous introns, is in preparation.
Programs for computations of snoRNA
We created a SNO.pl program for the identification of snoRNA-like sequences within MOID. This program was written in the PERL language and is publicly available from our web page http://www.meduohio.edu/bioinfo/eid/index.html. SNO.pl scans all human introns for the characteristic conserved C/D-box snoRNA elements defined by PERL regular expressions. Then, for those human introns with snoRNA-like structures, their orthologous mouse and rat intron sequences were extracted from MOID and were analyzed for the presence of the same snoRNA conserved elements. snoRNA-like sequences that were present in all three orthologous introns within human, mouse and rat genomes were output in a special file and are provided as a supplementary document. It takes only 4 min for a desktop computer (AMD Athlon 2200+ processor) to run this program. Further, simple modifications of the regular expressions inside this PERL script should allow the search for any type of conserved motifs within orthologous introns.
Another program written by using PERL, TARGET.pl, was designed to search for all putative targets for snoRNA ASE within databases of rRNA, snRNA and mRNA sequences. No mismatches were allowed in this ASE-target pairing, except non-Watson–Crick (non-WC) G–T pairs. The maximal number of G–T pairs were defined by the user. The execution time of this program was <2 min and it is also available from the same website.
RESULTS AND DISCUSSION
Computation of C/D-box snoRNA-like structures in human and rodents
Using our program SNO.pl, we scanned the entire set of 238 014 human introns from the EID looking for the following conserved structures characteristic for C/D-box snoRNA: (6 nt 5'-stem)-(N)-(C-box)-(loop 40–100 nt)-(D-box)-(6 nt 3'-stem). The C-box is RTGATGA, N stands for any nucleotide, the D-box is CTGA, and the loop can include any nucleotide sequence. In addition, the 5'-stem must have complementarity with the 3'-stem and a score of at least 4 points (each Watson–Crick pair adds 1 point, non-WC G–T pair adds 0.5 points; mismatch: 0 points). In total, 3693 of such snoRNA-like structures were found and selected within 3382 human introns. Then, for each selected human intron, its orthologous counterpart in mouse and rat were searched in the MOID. As a result, 1441 mouse orthologous intron sequences and 1079 rat orthologous intron sequences were obtained and were used to search for the same conserved C/D-box snoRNA elements within these orthologous introns. Finally, 224 snoRNA-like sequences were found in 193 mouse orthologous introns and 124 were found in 108 rat orthologous introns. Intersection of these data represents 57 orthologous intron triplets from human, mouse and rat genomes. Every intron from these 57 triplets contains computed snoRNA-like structures that are described in Table 1. All of the sequences and the supporting materials described in this table are available in the Supplementary File ‘57snoRNA.doc’.
Table 1 Description of computed snoRNA and snoRNA-like sequences detected inside each intron of 57 orthologous intron triplets from human, mouse and rat genomes
As shown in Table 1, we found 16 known snoRNA sequences inside 15 orthologous intron triplets (the 19th intron of the human predicted gene KIAA1731 contains two snoRNA molecules, mgU2-19/30 and Z32). Another six snoRNA-like structures found in the set of 57 orthologous intron triplets represent real snoRNA molecules because all six computer-detected snoRNA-like sequences from mouse were identical to partially characterized murine snoRNA sequences that were detected in a large-scale experimental approach (19). A thorough description of these six cases is presented in Table 2, while their sequences are illustrated in Figure 2. One snoRNA-like structure presumably represents an unknown snoRNA gene. It has 88% sequence identity with the human U53 snoRNA and its 59–79 nt fragment is identical to mouse MBII-35 snoRNA. The predicted 12 nt long ASE-2 of this novel gene has a target on the 28S rRNA with one G–T non-WC base pair (Table 2 and Figure 2). Another three snoRNA-like sequences represent putative snoRNA genes with less certainty in respect to their functionality. With one exception, we did not find convincing targets for their ASEs among rRNA, ncRNA from RNAdb and mRNA sequences (Table 2 and Figure 2). Finally, the remaining 56% of introns with snoRNA-like structures (32 out of 57 orthologous triplets) most probably manifested false-positive results because they were found in extra-long introns and in most cases they did not share significant sequence similarity that we observed in real and putative snoRNAs. All of these cases were easy to detect and could be filtered out since our SNO.pl program output the length of each intron with the snoRNA-like structure. Yet, one should keep in mind that long introns could bear novel snoRNA molecules that do not have homologs in other species. The Supplementary File 57snoRNA.doc presents complete information on the 57 predicted sequences and shows that many extra-long introns contain several snoRNA-like sequences per intron.
Table 2 Description of partially characterized, novel and putative snoRNAs, detected by the SNO.pl program
Figure 2 Sequences and conserved motifs of partially characterized, novel and putative snoRNAs that were detected by the SNO.pl program. C-, C'-, D- and D'-boxes are boxed. ASE-1s are underlined by a single line, ASE-2s are underlined by a double line. Hypothetical ASEs, which do not have strong targets, are underlined by a dotted line. All ASE targets are listed in Table 2.
Computation of snoRNA targets
As described by Huttenhofer et al. (20), the C/D-box snoRNA ASEs are 9–20 nt long sequences with 1 nt upstream of the D-boxes and are complementary to the ncRNA molecules that they are guiding for modification. Up to three G–T non-WC pairs are allowed within this ASE-target pairing and the rest of the bases should have perfect complementarity. We have generated and applied the TARGET.pl program to the search for putative targets for our predicted snoRNA-like structures in the databases of rRNA, ncRNA and mRNA sequences. For searching the putative targets within rRNA, we required that the length of ASE be at least 9 nt (L 9) and that the maximal number of non-WC G–T pairs should not exceed three (GT 3). In the search for targets among the ncRNA database, we required (L 12, GT 3), and among the mRNA database (L 16, GT 3). All calculated putative targets are described in Table 2 and their corresponding ASE sequences are underlined in Figure 2.
Simultaneously computing the snoRNA structures for three different species assists in validating the ASE targets and also sheds light on the coevolution of the ASEs and their target sequences. For instance, the tttcgactc ASE1 for the human counterpart of the MBII-316 snoRNA has a computed gagtcGggg target on the 28S rRNA (positions 1340–1348) with three non-WC G–T pairs (as shown in the first row of Table 2 and in Figure 2). In mouse and rat, this target on the 28S rRNA has a single nucleotide GA change in the middle of the corresponding segment gagtcAggg on the 28S rRNA when compared with human (this difference between human and rodents is shown in uppercase). The corresponding mouse putative ASE1 tttcgactc is identical to human and, therefore, owing to this mutation in the rRNA sequence, it most probably cannot guide modification of this site on the 28S rRNA. At the same time, there is a compensatory mutation of 2 nt for the rat ASE1 ttctgactc that restores the complementarity of the rat ASE1 with the corresponding site gagtcaggg on the 28S rRNA (positions 1252–1260). Such coevolution of the rat ASE1 and its computed target on the 28S rRNA supports the assumption that this putative ASE1 is a functional guiding sequence for chemical modification of the target on the human and rat 28S rRNA.
In contrast to the results with rRNAs, our results for possible ASE targets in mRNA, obtained from the mRNA database using our TARGET.pl program, were not as strong as the one described by Cavaille et al. (8) for ASE1 of MBII-52. All of our computed targets within mRNA corresponded to different genes in human, mouse and rat; thus, we were not inclined to presume that these targets were functional.
In conclusion, many snoRNA molecules have variations in the conserved C- and D-boxes and also in the terminal stem structures; thus, only a small fraction of snoRNA-like sequences from human and rodent genomes have been detected in this search. Modifications of the search pattern in the SNO.pl program should reveal many more putative snoRNAs. Different types of ncRNAs (C/D-box snoRNA, ACA/H-box snoRNA, microRNA and probably others) are also located inside introns. Hence, our approach can be easily adapted for searching all kinds of ncRNAs and functional motifs inside introns.
SUPPLEMENTARY MATERIAL
Supplementary Material is available at NAR Online.
ACKNOWLEDGEMENTS
We would like to thank Dr Robert Blumenthal, Medical University of Ohio, for discussion and valuable suggestions on our manuscript. Support for this work was provided by the Medical University of Ohio Foundation and the Stranahan Foundation, through the Program in Bioinformatics and Proteomics/Genomics. Funding to pay the Open Access publication charges for this article was provided by start-up fund of A.F.
REFERENCES
Fatica, A. and Tollervey, D. (2003) Insights into the structure and function of a guide RNP Nature Struct. Biol., 10, 237–239 .
Bachellerie, J.P., Cavaille, J., Huttenhofer, A. (2002) The expanding snoRNA world Biochimie, 84, 775–790 .
Huttenhofer, A., Brosius, J., Bachellerie, J.P. (2002) RNomics: identification and function of small, non-messenger RNAs Curr. Opin. Chem. Biol., 6, 835–843 .
Maxwell, E.S. and Fournier, M.J. (1995) The small nucleolar RNAs Annu. Rev. Biochem., 64, 897–934 .
Weinstein, L.B. and Steitz, J.A. (1999) Guided tours: from precursor snoRNA to functional snoRNP Curr. Opin. Cell. Biol., 11, 378–384 .
Aittaleb, M., Rashid, R., Chen, Q., Palmer, J.R., Daniels, C.J., Li, H. (2003) Structure and function of archaeal box C/D sRNP core proteins Nature Struct. Biol., 10, 256–263 .
Vitali, P., Royo, H., Seitz, H., Bachellerie, J.P., Huttenhofer, A., Cavaille, J. (2003) Identification of 13 novel human modification guide RNAs Nucleic Acids Res., 31, 6543–6551 .
Cavaille, J., Buiting, K., Kiefmann, M., Lalande, M., Brannan, C.I., Horsthemke, B., Bachellerie, J.P., Brosius, J., Huttenhofer, A. (2000) Identification of brain-specific and imprinted small nucleolar RNA genes exhibiting an unusual genomic organization Proc. Natl Acad. Sci. USA, 97, 14311–14316 .
Lowe, T.M. and Eddy, S.R. (1999) A computational screen for methylation guide snoRNAs in yeast Science, 283, 1168–1171 .
Accardo, M.C., Giordano, E., Riccardo, S., Digilio, F.A., Iazzetti, G., Calogero, R.A., Furia, M. (2004) A computational search for box C/D snoRNA genes in the Drosophila melanogaster genome Bioinformatics, 20, 3293–3301 .
Huang, Z.P., Zhou, H., Liang, D., Qu, L.H. (2004) Different expression strategy: multiple intronic gene clusters of box H/ACA snoRNA in Drosophila melanogaster J. Mol. Biol., 341, 669–683 .
Schattner, P., Decatur, W.A., Davis, C.A., Ares, M., Jr, Fournier, M.J., Lowe, T.M. (2004) Genome-wide searching for pseudouridylation guide snoRNAs: analysis of the Saccharomyces cerevisiae genome Nucleic Acids Res., 32, 4281–4296 .
Fedorov, A., Roy, S., Fedorova, L., Gilbert, W. (2003) Mystery of intron gain Genome Res., 13, 2236–2241 .
Roy, S.W., Fedorov, A., Gilbert, W. (2003) Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain Proc. Natl Acad. Sci. USA, 100, 7158–7162 .
Cullen, B.R. (2004) Transcription and processing of human microRNA precursors Mol. Cell, 16, 861–865 .
Fedorova, L. and Fedorov, A. (2003) Introns in gene evolution Genetica, 118, 123–131 .
Saxonov, S., Daizadeh, I., Fedorov, A., Gilbert, W. (2000) EID: the Exon–Intron Database—an exhaustive database of protein-coding intron-containing genes Nucleic Acids Res., 28, 185–190 .
Fedorov, A., Merican, A.F., Gilbert, W. (2002) Large-scale comparison of intron positions between plant, animal and fungal genes Proc. Natl Acad. Sci. USA, 99, 16128–16133 .
Huttenhofer, A., Kiefmann, M., Meier-Ewert, S., O'Brien, J., Lehrach, H., Bachellerie, J.P., Brosius, J. (2001) RNomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger RNAs in mouse EMBO J., 20, 2943–2953 .
Huttenhofer, A., Cavaille, J., Bachellerie, J.P. (2004) Experimental RNomics: a global approach to identifying small nuclear RNAs and their targets in different model organisms Methods Mol. Biol., 265, 409–428 .(Alexei Fedorov*, Jesse Stombaugh1, Micha)
*To whom correspondence should be addressed. Tel: +1 419 383 5270; Fax: +1 419 383 3102; Email: afedorov@meduohio.edu
ABSTRACT
Based on comparative genomics, we created a bioinformatic package for computer prediction of small nucleolar RNA (snoRNA) genes in mammalian introns. The core of our approach was the use of the Mammalian Orthologous Intron Database (MOID), which contains all known introns within the human, mouse and rat genomes. Introns from orthologous genes from these three species, that have the same position relative to the reading frame, are grouped in a special orthologous intron table. Our program SNO.pl searches for conserved snoRNA motifs within MOID and reports all cases when characteristic snoRNA-like structures are present in all three orthologous introns of human, mouse and rat sequences. Here we report an example of the SNO.pl usage for searching a particular pattern of conserved C/D-box snoRNA motifs (canonical C- and D-boxes and the 6 nt long terminal stem). In this computer analysis, we detected 57 triplets of snoRNA-like structures in three mammals. Among them were 15 triplets that represented known C/D-box snoRNA genes. Six triplets represented snoRNA genes that had only been partially characterized in the mouse genome. One case represented a novel snoRNA gene, and another three cases, putative snoRNAs. Our programs are publicly available and can be easily adapted and/or modified for searching any conserved motifs within mammalian introns.
INTRODUCTION
Small nucleolar RNA (snoRNA) is a major component of small nucleolar ribonucleoprotein (snoRNP) particles that are located inside the nucleolus of eukaryotes and participate in post-transcriptional chemical modification or processing of different RNAs, including ribosomal RNA (rRNA) and small nuclear RNA (snRNA) . SnoRNAs are ancient genes, as they are widespread through the entire domain of Eukaryota, as well as Archaea. There are two types of these RNAs, called C/D-box and H/ACA-box snoRNAs, and they are characterized by distinct 3D structures and conserved sequence elements. The C/D-box RNA is associated with 2'-O-ribose methylation, and H/ACA-box RNA with pseudouridylation of substrate RNAs. The direct role of snoRNAs is determining the site for chemical modification, via complementary pairing of its specific sequence with the segment of substrate RNA undergoing modification. The major catalytic activity does not belong to snoRNA, but rather to a fibrillarin, a protein component of the snoRNP complex (1,6). Besides pseudouridylation and 2'-O-ribose methylation, some snoRNAs perform other documented and putative functions, such as (i) cleavage of the substrate precursor in rRNA (2), (ii) facilitation of rRNA folding (7), (iii) regulation of alternative splicing (8) and (iv) possible unknown functions carried out by so-called ‘orphan’ snoRNAs (2).
The first well-known computer program for genomic prediction of C/D-box snoRNA, named SNOSCAN, was created by Lowe and Eddy (9). SNOSCAN was used by several laboratories (10) to study compact genomes of eukaryotes such as yeast, Drosophila and Arabidopsis. SNOSCAN can predict classical C/D-box snoRNAs that guide modification of rRNA or snRNA molecules, and requires that these rRNA and snRNA sequences be available for the program before its invocation. In 2003, Vitali et al. (7) used this program on a restricted set of human and mouse ribosomal protein gene sequences, and found several novel C/D snoRNAs inside their introns. Recently, two computational approaches were developed to predict H/ACA snoRNAs from genomic sequences (11,12). Again, these programs can predict snoRNAs that are involved in chemical modifications of known rRNA and snRNA sequences, and they are suitable for searching compact eukaryotic genomes. However, it is problematic to apply the existing computational algorithms for studying vertebrate genomes in their entirety, because they are dozens of times larger than the genomes of yeast, insects or Arabidopsis. In the case of vertebrates, numerous false-positive signals may arise due to the computer processing of very large non-coding genomic sequences and thus computer prediction of novel genes may not be as efficient.
We have taken advantage of comparative genomics and present here a new algorithm for computational identification of mammalian C/D-box snoRNA genes with high efficiency. We exploit two well-known characteristics of mammalian genomic structure: (i) all known snoRNA genes of mammals are located inside introns (3) and (ii) the exon–intron gene structure of mammals is highly conserved. In fact, no well-documented case of intron gain exists between fully sequenced mammalian genomes, and there are only a few intron losses reported (13,14). Therefore, introns are ‘fossilized’ in mammalian genes, and the non-coding RNAs (ncRNAs, also known as small non-messenger RNAs, non-protein coding RNAs and untranslated RNAs) located within them could well be fixed inside the same intron since the origin of this taxon. For this reason, we created a database of orthologous introns comprising human, mouse and rat sequences. We define ‘orthologous introns’ as introns from orthologous genes that have the same position relative to the coding sequence. Hence, orthologous introns should have descended from the corresponding intronic sequence of the last common ancestor for the taxon.
The evolutionary divergence of primate and rodent lineages occurred 70–90 million years ago, and during this period non-functional DNA segments representing a major portion of human and rodent orthologous introns lost almost all sequence similarity. For example, Figure 1 demonstrates an alignment of the third introns of the ribosomal protein S3a gene from human and mouse genomes that are orthologous to each other and are 1708 and 1990 nt long, respectively. It is known that these two introns contain U73b snoRNA sequences of human and mouse (NG_000961 and Z83331 , respectively). A BLAST-2 alignment with default parameters reveals 88% sequence identity of a 58 nt long fragment of these introns containing snoRNA, and no sequence similarity outside this short functional intronic fragment except one 37 nt long segment (Figure 1). This second conserved segment could also represent a functional sequence not yet known. Therefore, when conserved snoRNA structures are identified within orthologous intron triplets from all three species (human, mouse and rat), it strongly suggests that the conserved sequence is functional. Here, we present our Mammalian Orthologous Intron Database (MOID) for public usage and two programs that search for conserved C/D-box snoRNA motifs inside the entire MOID and characterize putative ASEs of snoRNAs based on comparative genomics of mammalian species.
Figure 1 Results of an online Blast2 alignment of the 1708 nt long third intron of the human ribosomal protein S3a gene, and its mouse ortholog (1990 nt long third intron of the mouse rps3a gene). (A) Dot plot figure of the intron comparison obtained by bl2seq online program (http://www.ncbi.nlm.nih.gov/blast/bl2seq/wblast2.cgi) with default parameters. (B) Alignment of human and mouse introns. The sequences of snoRNA C- and D-boxes inside these introns are boxed.
Mammalian introns bear all types of snoRNA genes (3), about three-fourth of microRNA genes (15), and also many conserved regulatory motifs (16). Our programs can be easily adapted to search for all kinds of functional sequences using MOID. We are further developing this approach by extending our Orthologous Intron Database with other mammalian and vertebrate species.
MATERIALS AND METHODS
Intron databases
Using the previously published program package for the generation of the Exon–Intron Database (EID) (17), we created three species-specific EIDs containing all known genes of mammals with entirely sequenced genomes. These EID for human, mouse and rat genomes and their documentation are freely available from our website (http://www.meduohio.edu/bioinfo/eid/index.html). We then applied the CIP.pl program (18) for the comparison of intron positions in all orthologous genes from human, mouse and rat genomes. Based on this comparison, we generated an MOID, which is also publicly available from the same web page (documentation provided). This first release of our MOID contains 100 000 orthologous introns from human, mouse and rat genomes. The next release of MOID, with a complete set of orthologous introns, is in preparation.
Programs for computations of snoRNA
We created a SNO.pl program for the identification of snoRNA-like sequences within MOID. This program was written in the PERL language and is publicly available from our web page http://www.meduohio.edu/bioinfo/eid/index.html. SNO.pl scans all human introns for the characteristic conserved C/D-box snoRNA elements defined by PERL regular expressions. Then, for those human introns with snoRNA-like structures, their orthologous mouse and rat intron sequences were extracted from MOID and were analyzed for the presence of the same snoRNA conserved elements. snoRNA-like sequences that were present in all three orthologous introns within human, mouse and rat genomes were output in a special file and are provided as a supplementary document. It takes only 4 min for a desktop computer (AMD Athlon 2200+ processor) to run this program. Further, simple modifications of the regular expressions inside this PERL script should allow the search for any type of conserved motifs within orthologous introns.
Another program written by using PERL, TARGET.pl, was designed to search for all putative targets for snoRNA ASE within databases of rRNA, snRNA and mRNA sequences. No mismatches were allowed in this ASE-target pairing, except non-Watson–Crick (non-WC) G–T pairs. The maximal number of G–T pairs were defined by the user. The execution time of this program was <2 min and it is also available from the same website.
RESULTS AND DISCUSSION
Computation of C/D-box snoRNA-like structures in human and rodents
Using our program SNO.pl, we scanned the entire set of 238 014 human introns from the EID looking for the following conserved structures characteristic for C/D-box snoRNA: (6 nt 5'-stem)-(N)-(C-box)-(loop 40–100 nt)-(D-box)-(6 nt 3'-stem). The C-box is RTGATGA, N stands for any nucleotide, the D-box is CTGA, and the loop can include any nucleotide sequence. In addition, the 5'-stem must have complementarity with the 3'-stem and a score of at least 4 points (each Watson–Crick pair adds 1 point, non-WC G–T pair adds 0.5 points; mismatch: 0 points). In total, 3693 of such snoRNA-like structures were found and selected within 3382 human introns. Then, for each selected human intron, its orthologous counterpart in mouse and rat were searched in the MOID. As a result, 1441 mouse orthologous intron sequences and 1079 rat orthologous intron sequences were obtained and were used to search for the same conserved C/D-box snoRNA elements within these orthologous introns. Finally, 224 snoRNA-like sequences were found in 193 mouse orthologous introns and 124 were found in 108 rat orthologous introns. Intersection of these data represents 57 orthologous intron triplets from human, mouse and rat genomes. Every intron from these 57 triplets contains computed snoRNA-like structures that are described in Table 1. All of the sequences and the supporting materials described in this table are available in the Supplementary File ‘57snoRNA.doc’.
Table 1 Description of computed snoRNA and snoRNA-like sequences detected inside each intron of 57 orthologous intron triplets from human, mouse and rat genomes
As shown in Table 1, we found 16 known snoRNA sequences inside 15 orthologous intron triplets (the 19th intron of the human predicted gene KIAA1731 contains two snoRNA molecules, mgU2-19/30 and Z32). Another six snoRNA-like structures found in the set of 57 orthologous intron triplets represent real snoRNA molecules because all six computer-detected snoRNA-like sequences from mouse were identical to partially characterized murine snoRNA sequences that were detected in a large-scale experimental approach (19). A thorough description of these six cases is presented in Table 2, while their sequences are illustrated in Figure 2. One snoRNA-like structure presumably represents an unknown snoRNA gene. It has 88% sequence identity with the human U53 snoRNA and its 59–79 nt fragment is identical to mouse MBII-35 snoRNA. The predicted 12 nt long ASE-2 of this novel gene has a target on the 28S rRNA with one G–T non-WC base pair (Table 2 and Figure 2). Another three snoRNA-like sequences represent putative snoRNA genes with less certainty in respect to their functionality. With one exception, we did not find convincing targets for their ASEs among rRNA, ncRNA from RNAdb and mRNA sequences (Table 2 and Figure 2). Finally, the remaining 56% of introns with snoRNA-like structures (32 out of 57 orthologous triplets) most probably manifested false-positive results because they were found in extra-long introns and in most cases they did not share significant sequence similarity that we observed in real and putative snoRNAs. All of these cases were easy to detect and could be filtered out since our SNO.pl program output the length of each intron with the snoRNA-like structure. Yet, one should keep in mind that long introns could bear novel snoRNA molecules that do not have homologs in other species. The Supplementary File 57snoRNA.doc presents complete information on the 57 predicted sequences and shows that many extra-long introns contain several snoRNA-like sequences per intron.
Table 2 Description of partially characterized, novel and putative snoRNAs, detected by the SNO.pl program
Figure 2 Sequences and conserved motifs of partially characterized, novel and putative snoRNAs that were detected by the SNO.pl program. C-, C'-, D- and D'-boxes are boxed. ASE-1s are underlined by a single line, ASE-2s are underlined by a double line. Hypothetical ASEs, which do not have strong targets, are underlined by a dotted line. All ASE targets are listed in Table 2.
Computation of snoRNA targets
As described by Huttenhofer et al. (20), the C/D-box snoRNA ASEs are 9–20 nt long sequences with 1 nt upstream of the D-boxes and are complementary to the ncRNA molecules that they are guiding for modification. Up to three G–T non-WC pairs are allowed within this ASE-target pairing and the rest of the bases should have perfect complementarity. We have generated and applied the TARGET.pl program to the search for putative targets for our predicted snoRNA-like structures in the databases of rRNA, ncRNA and mRNA sequences. For searching the putative targets within rRNA, we required that the length of ASE be at least 9 nt (L 9) and that the maximal number of non-WC G–T pairs should not exceed three (GT 3). In the search for targets among the ncRNA database, we required (L 12, GT 3), and among the mRNA database (L 16, GT 3). All calculated putative targets are described in Table 2 and their corresponding ASE sequences are underlined in Figure 2.
Simultaneously computing the snoRNA structures for three different species assists in validating the ASE targets and also sheds light on the coevolution of the ASEs and their target sequences. For instance, the tttcgactc ASE1 for the human counterpart of the MBII-316 snoRNA has a computed gagtcGggg target on the 28S rRNA (positions 1340–1348) with three non-WC G–T pairs (as shown in the first row of Table 2 and in Figure 2). In mouse and rat, this target on the 28S rRNA has a single nucleotide GA change in the middle of the corresponding segment gagtcAggg on the 28S rRNA when compared with human (this difference between human and rodents is shown in uppercase). The corresponding mouse putative ASE1 tttcgactc is identical to human and, therefore, owing to this mutation in the rRNA sequence, it most probably cannot guide modification of this site on the 28S rRNA. At the same time, there is a compensatory mutation of 2 nt for the rat ASE1 ttctgactc that restores the complementarity of the rat ASE1 with the corresponding site gagtcaggg on the 28S rRNA (positions 1252–1260). Such coevolution of the rat ASE1 and its computed target on the 28S rRNA supports the assumption that this putative ASE1 is a functional guiding sequence for chemical modification of the target on the human and rat 28S rRNA.
In contrast to the results with rRNAs, our results for possible ASE targets in mRNA, obtained from the mRNA database using our TARGET.pl program, were not as strong as the one described by Cavaille et al. (8) for ASE1 of MBII-52. All of our computed targets within mRNA corresponded to different genes in human, mouse and rat; thus, we were not inclined to presume that these targets were functional.
In conclusion, many snoRNA molecules have variations in the conserved C- and D-boxes and also in the terminal stem structures; thus, only a small fraction of snoRNA-like sequences from human and rodent genomes have been detected in this search. Modifications of the search pattern in the SNO.pl program should reveal many more putative snoRNAs. Different types of ncRNAs (C/D-box snoRNA, ACA/H-box snoRNA, microRNA and probably others) are also located inside introns. Hence, our approach can be easily adapted for searching all kinds of ncRNAs and functional motifs inside introns.
SUPPLEMENTARY MATERIAL
Supplementary Material is available at NAR Online.
ACKNOWLEDGEMENTS
We would like to thank Dr Robert Blumenthal, Medical University of Ohio, for discussion and valuable suggestions on our manuscript. Support for this work was provided by the Medical University of Ohio Foundation and the Stranahan Foundation, through the Program in Bioinformatics and Proteomics/Genomics. Funding to pay the Open Access publication charges for this article was provided by start-up fund of A.F.
REFERENCES
Fatica, A. and Tollervey, D. (2003) Insights into the structure and function of a guide RNP Nature Struct. Biol., 10, 237–239 .
Bachellerie, J.P., Cavaille, J., Huttenhofer, A. (2002) The expanding snoRNA world Biochimie, 84, 775–790 .
Huttenhofer, A., Brosius, J., Bachellerie, J.P. (2002) RNomics: identification and function of small, non-messenger RNAs Curr. Opin. Chem. Biol., 6, 835–843 .
Maxwell, E.S. and Fournier, M.J. (1995) The small nucleolar RNAs Annu. Rev. Biochem., 64, 897–934 .
Weinstein, L.B. and Steitz, J.A. (1999) Guided tours: from precursor snoRNA to functional snoRNP Curr. Opin. Cell. Biol., 11, 378–384 .
Aittaleb, M., Rashid, R., Chen, Q., Palmer, J.R., Daniels, C.J., Li, H. (2003) Structure and function of archaeal box C/D sRNP core proteins Nature Struct. Biol., 10, 256–263 .
Vitali, P., Royo, H., Seitz, H., Bachellerie, J.P., Huttenhofer, A., Cavaille, J. (2003) Identification of 13 novel human modification guide RNAs Nucleic Acids Res., 31, 6543–6551 .
Cavaille, J., Buiting, K., Kiefmann, M., Lalande, M., Brannan, C.I., Horsthemke, B., Bachellerie, J.P., Brosius, J., Huttenhofer, A. (2000) Identification of brain-specific and imprinted small nucleolar RNA genes exhibiting an unusual genomic organization Proc. Natl Acad. Sci. USA, 97, 14311–14316 .
Lowe, T.M. and Eddy, S.R. (1999) A computational screen for methylation guide snoRNAs in yeast Science, 283, 1168–1171 .
Accardo, M.C., Giordano, E., Riccardo, S., Digilio, F.A., Iazzetti, G., Calogero, R.A., Furia, M. (2004) A computational search for box C/D snoRNA genes in the Drosophila melanogaster genome Bioinformatics, 20, 3293–3301 .
Huang, Z.P., Zhou, H., Liang, D., Qu, L.H. (2004) Different expression strategy: multiple intronic gene clusters of box H/ACA snoRNA in Drosophila melanogaster J. Mol. Biol., 341, 669–683 .
Schattner, P., Decatur, W.A., Davis, C.A., Ares, M., Jr, Fournier, M.J., Lowe, T.M. (2004) Genome-wide searching for pseudouridylation guide snoRNAs: analysis of the Saccharomyces cerevisiae genome Nucleic Acids Res., 32, 4281–4296 .
Fedorov, A., Roy, S., Fedorova, L., Gilbert, W. (2003) Mystery of intron gain Genome Res., 13, 2236–2241 .
Roy, S.W., Fedorov, A., Gilbert, W. (2003) Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain Proc. Natl Acad. Sci. USA, 100, 7158–7162 .
Cullen, B.R. (2004) Transcription and processing of human microRNA precursors Mol. Cell, 16, 861–865 .
Fedorova, L. and Fedorov, A. (2003) Introns in gene evolution Genetica, 118, 123–131 .
Saxonov, S., Daizadeh, I., Fedorov, A., Gilbert, W. (2000) EID: the Exon–Intron Database—an exhaustive database of protein-coding intron-containing genes Nucleic Acids Res., 28, 185–190 .
Fedorov, A., Merican, A.F., Gilbert, W. (2002) Large-scale comparison of intron positions between plant, animal and fungal genes Proc. Natl Acad. Sci. USA, 99, 16128–16133 .
Huttenhofer, A., Kiefmann, M., Meier-Ewert, S., O'Brien, J., Lehrach, H., Bachellerie, J.P., Brosius, J. (2001) RNomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger RNAs in mouse EMBO J., 20, 2943–2953 .
Huttenhofer, A., Cavaille, J., Bachellerie, J.P. (2004) Experimental RNomics: a global approach to identifying small nuclear RNAs and their targets in different model organisms Methods Mol. Biol., 265, 409–428 .(Alexei Fedorov*, Jesse Stombaugh1, Micha)