- Research article
- Open Access
Genome-wide analysis of chicken snoRNAs provides unique implications for the evolution of vertebrate snoRNAs
BMC Genomics volume 10, Article number: 86 (2009)
Small nucleolar RNAs (snoRNAs) represent one of the largest groups of functionally diverse trans-acting non-protein-coding (npc) RNAs currently known in eukaryotic cells. Chicken snoRNAs have been very poorly characterized when compared to other vertebrate snoRNAs. A genome-wide analysis of chicken snoRNAs is therefore of great importance to further understand the functional evolution of snoRNAs in vertebrates.
Two hundred and one gene variants encoding 93 box C/D and 62 box H/ACA snoRNAs were identified in the chicken genome and are predicted to guide 86 2'-O-ribose methylations and 69 pseudouridylations of rRNAs and spliceosomal RNAs. Forty-four snoRNA clusters were grouped into four categories based on synteny characteristics of the clustered snoRNAs between chicken and human. Comparative analyses of chicken snoRNAs revealed extensive recombination and separation of guiding function, with cooperative evolution between the guiding duplexes and modification sites. The gas5-like snoRNA host gene appears to be a hotspot of snoRNA gene expansion in vertebrates. Our results suggest that the chicken is a good model for the prediction of functional snoRNAs, and that intragenic duplication and divergence might be the major driving forces responsible for expansion of novel snoRNA genes in the chicken genome.
We have provided a detailed catalog of chicken snoRNAs that aids in understanding snoRNA gene repertoire differences between avians and other vertebrates. Our genome-wide analysis of chicken snoRNAs improves annotation of the 'darkness matter' in the npcRNA world and provides a unique perspective into snoRNA evolution in vertebrates.
The term small nucleolar RNAs (snoRNAs) was originally coined to describe the nucleolar localization of this group of RNAs relative to the other small nucleoplasmic RNAs. In sharp contrast to the relatively low abundance spliceosomal nuclear RNA (snRNA) species, snoRNAs represent one of the largest groups of functionally diverse trans-acting non-protein-coding RNAs (npcRNAs) currently known in eukaryotic cells [1, 2]. On the basis of conserved sequence elements and characteristic secondary structures, snoRNAs can be divided into two major classes, box C/D and box H/ACA snoRNAs. Box C/D snoRNAs contain two conserved motifs, the 5' end box C (RUGAUGA, where R stands for any purine) and the 3' end box D (CUGA). Box H/ACA snoRNAs exhibit a common hairpin-hinge-hairpin-tail secondary structure with the H box (ANANNA, where N stands for any nucleotide) in the hinge region and the ACA motif three nucleotides from the 3' end of the molecule. During the post-transcriptional processing of diverse RNAs most members of the known C/D and H/ACA snoRNAs respectively guide 2'-O-ribose methylation and pseudouridylation (Ψ). Recently, a new class of guide RNAs has been found to accumulate in the small Cajal body  and are thus termed small Cajal body-specific RNAs (scaRNAs). scaRNAs are often composed of both C/D box and H/ACA box domains  and guide the modification of RNA-polymerase-II-transcribed snRNAs . Remarkably, an increasing number of 'orphan' snoRNAs lacking antisense to known RNA targets have been identified . Many of them exhibit a tissue-specific or restricted expression pattern [6, 7] and are linked to genomic imprinting .
Interestingly, various snoRNA gene organizations have been characterized in different organisms [5, 8]. Most snoRNAs are encoded in the introns of protein-coding or non-protein-coding genes in vertebrates . Many snoRNA paralogs are usually clustered in different introns of the same host genes (HGs) or in the introns of different HGs by intragenic or intergenic duplication (including retroposition) from existing snoRNAs [7, 10–13], respectively. The distinct character of clustering gene organizations and evolutionary conservation of vertebrate snoRNAs facilitates detection of snoRNA homologs by sequence similarity alone in the genome . However, many other snoRNAs in mammals cannot be found by simple homology search.
To date, hundreds of snoRNAs have been identified in mammals [7, 13, 15–18] by approaches including computational and experimental RNomics. Although a limited number of snoRNAs were predicted in the chicken (Gallus gallus) genome by similarity search , the nature of chicken snoRNAs is poorly understood when compared with other vertebrates and their numbers far underrepresented. Additionally, detailed information on snoRNA guiding functions, genomic organization and evolution in the chicken genome is still unavailable. As a typical amniote, the chicken has evolved separately from mammals for about 310 million years . The identification of chicken snoRNAs using conventional prediction methods such as a similarity search might be hindered by the sufficient nucleotide variation occurring in the genome. Recently, we developed an advanced computational package snoSeeker for the specific detection of guide box C/D (CDseeker) and box H/ACA (ACAseeker) snoRNAs, as well as orphan snoRNA genes in the human genome . In the present work, 93 box C/D and 62 box H/ACA snoRNAs have been identified in the chicken genome by applying the computational package and experimental methods based on RT-PCR. The characteristics of the guiding function and genomic organization of the chicken snoRNAs have been extensively compared with the human counterparts. As a result, we provide for the first time a detailed catalog of chicken snoRNAs that facilitates understanding of snoRNA gene repertoire differences between the avian and other vertebrate lineages.
Computational identification of box C/D snoRNA genes from G. gallus
The CDseeker program was applied to search the G.gallus genome for box C/D snoRNAs. In total, 132 gene variants encoding 83 box C/D snoRNAs with the ability to guide 2'-O-ribose methylation at 86 residues in rRNAs and snRNAs and 10 orphan box C/D snoRNAs were identified from the G.gallus genome (Table 1, see Additional file 1 and 2). Sixty-five box C/D snoRNAs are singleton. The other 28 snoRNAs have undergone one or more duplications in the chicken, which account for 67 paralogs. The majority of these RNAs (73 snoRNAs) have been assigned to guide only one methylation of the rRNAs or snRNAs, and are known as single-guide snoRNAs. Seventy-one snoRNAs uniquely guide methylations of the rRNAs, and 12 box C/D RNAs are predicted to guide methylation in the snRNAs or both the rRNAs and snRNAs. Interestingly, two methylation sites at 18S rRNA-C757 and U2-U47 (corresponding to human 18S rRNA-C797 and U2-U47), which have been previously reported to lack potential guide snoRNAs, were predicted to be guided by GGgCD20 and GGgCD76, respectively. Among the 83 guide RNAs, only a fraction of box C/D snoRNAs (~15%) are double-guide snoRNAs. A comparative analysis of these chicken box C/D snoRNAs and their counterparts at the corresponding genomic loci in six other vertebrate genomes (human, mouse, opossum, platypus, lizard and frog) revealed that 13 snoRNA genes appeared to be specific to the chicken or avian lineage, and the remainder had their cognate snoRNAs in at least one other vertebrate genome (see Additional file 3). In total, 55 chicken box C/D snoRNA genes are conserved in amniotic species. These snoRNAs are assigned as the core amniotic box C/D snoRNAs, out of which 31 box C/D snoRNA genes are the core vertebrate box C/D snoRNAs conserved in vertebrates (Figure 1A).
Computational identification of box H/ACA snoRNA genes from G. gallus
Based on the conserved 'hairpin-hinge-hairpin-tail' structure and the H and ACA/ATA box motifs, the ACAseeker program was performed to identify chicken box H/ACA snoRNA genes. In total, 69 RNA variants encoding 52 guide and 10 orphan box H/ACA snoRNAs were identified in the chicken genome (Table 2, see Additional file 2). Compared with the chicken box C/D snoRNAs, a higher percentage of box H/ACA snoRNAs (~90%) are singleton. Only six snoRNAs had undergone one or more duplications, yielding 13 paralogs. Fifty-two guide box H/ACA RNAs were predicted to guide 69 Ψs in rRNAs and snRNAs (see Additional file 1). The majority of these RNAs (42 snoRNAs) have been assigned to guide Ψs of the rRNAs. Nine scaRNAs have been predicted to guide Ψs of the snRNAs (GGgACA48~52) or both the rRNAs and snRNAs (GGgACA10, GGgACA14, GGgACA18 and GGgACA29). Three Ψs, 28S rRNA-Ψ3751, U2-Ψ58, U2-Ψ91 (corresponding to human 28S rRNA-Ψ4266, U2-Ψ58 and U2-Ψ91), which have been previously reported to lack a potential guide snoRNA, were predicted to be guided by GGgACA43, GGgACA18 and GGgACA10, respectively. In sharp contrast to the box C/D snoRNAs, approximately half of the H/ACA snoRNAs are capable of directing more than one Ψ which are often located on the same rRNA or snRNA. Interestingly, GGgACA5 and GGgACA6 show the potential of directing three nonadjacent Ψs by a single guide sequence and are the first to be reported in the chicken. In the case of GGgACA51, whose sequence shows high similarity to that of half of the U93 composed of two tandem arranged box H/ACA RNA domains , it should be annotated as the half U93 homolog present in humans. Comparison of the conserved regions encoding snoRNAs in the chicken and six other vertebrate genomes revealed that eight are found to be chicken-specific, and 26 are the core amniotic box H/ACA snoRNA genes which include 12 core vertebrate snoRNA genes (Figure 1B, see Additional file 3).
Experimental confirmation of the computational results
Recently, PCR-based methods have been successfully used for small RNA detection, as well as expression profiling [21, 22]. In this study, we developed an improved method for specifically detecting the expression of snoRNA candidates (see Additional file 4). As different expression levels might be detected for snoRNAs in different host genes, we first detected the expression of 14 guide snoRNAs located in different host genes to test the method we developed (Figure 2A). In most cases, unique and obvious bands within our expected sizes were detected under stringent conditions of PCR. Notably, the expression of four snoRNAs (GGgCD20, GGgCD76, GGgACA10 and GGgACA43) that are predicted to guide modification sites that have been previously reported to lack potential guide snoRNAs was detected. Two snoRNAs (GGgCD11a and GGgACA46) hosted within the ribosomal protein gene family exhibit a robust signal; whereas the other 12 snoRNAs, most of which are not located within house keeping genes, were detected only weakly. It remains to be further tested whether the different intensities of the electrophoresis bands observed under the same PCR conditions reflect the different expression levels of the snoRNAs in vivo. We next applied this method to detect the 20 orphan snoRNAs predicted in this study. All these orphan snoRNAs were amplified and the sizes of the PCR products were consistent with their predicted sizes (Figure 2B). Our experimental analysis by RT-PCR demonstrated that the snoRNA genes predicted in this study are expressed.
Clustered snoRNAs most often show conserved synteny between humans and chickens
Almost all of the snoRNAs identified in the chicken are located in introns of known genes or spliced expressed sequence tags (ESTs), which is consistent with previous reports in mammals[1, 9, 23]. The majority of chicken snoRNAs (80%) appear on chromosomes 1, 2, 4, 5, 8, 10 and Z. None of the snoRNA genes was located on chromosomes 21, 22, 26, 29, 30, 31 and W.
From this study, a total of 98 box C/D and 38 box H/ACA snoRNAs are organized into 44 clusters (see Additional file 5). Eleven clusters carry both box C/D and box H/ACA snoRNAs, and the other clusters harbor exclusively either box C/D or box H/ACA snoRNAs. Five clusters reside in the introns of non-protein-coding genes or ESTs or novel transcripts, and the remaining 39 clusters are located in protein-coding HGs. Genomic analysis was carried out on these clustered snoRNAs and their HGs. Eighty-six of 136 clustered snoRNAs (63%) in chicken are syntenic to human, and 36 protein-coding HGs are the orthologs of the human HGs. Four types of snoRNA clusters were determined depending on the synteny characteristics of the clustered snoRNAs between chicken and human (Figure 3).
The Type-1 snoRNA clusters (clusters 1–14) show perfect synteny between chicken and human (Figure 3A). The HGs of these clustered snoRNAs in chicken are orthologous to those in human, indicating that they are derived from common ancestral genes over 310 million years ago.
The Type-2 snoRNA clusters (clusters 15–33) are conserved in content but neither in snoRNA copy numbers nor in order. Fifteen HGs of the type-2 snoRNA clusters are orthologous between chicken and human, and the other snoRNA clusters reside in either non-protein-coding genes, novel transcripts or the ESTs of which their mRNA sequences show low or no similarity between the two species. Interestingly, the HG of cluster 15 is an Ensemble novel protein coding gene (ENSGALESTG00000033205) and tandemly encodes multiple chicken counterparts of human snoRNAs (ACA40, ACA18, mgh28S-2409, ACA8, ACA1, mgh28S-2411, ACA32 and ACA25) located within different introns of the same host gene JOSD3 (Figure 3B). However, these eight clustered chicken snoRNAs are arrayed in a different order when compared to those of the human and other vertebrates, which might have resulted from lineage-specific intragenic translocation (see Additional file 6). Sequence comparison revealed that the mRNA sequence of ENSGALESTG00000033205 showed low similarity (45%) to that of JOSD3, which is in sharp contrast to the high similarity of the snoRNA counterpart between the two species.
The Type-3 snoRNA clusters (clusters 34–37) have HGs that are orthologous to human genes but with the insertion of new snoRNAs. For example, four novel snoRNAs are inserted in cluster 34 when compared to that in human, and the HGs (DKC1) between chicken and human are orthologous (Figure 3C).
The Type-4 snoRNA clusters (clusters 38–44) have independently evolved in the chicken (Figure 3D). Corresponding snoRNAs could not be identified in human orthologous regions, suggesting that these snoRNA clusters might have changed their HGs during chicken and human divergence.
A gas5-like non-protein-coding snoRNA transcript found in the chicken
Whilst characterizing the synteny of the clustered snoRNAs between the chicken and human, we noted that snoRNA cluster 26 is of special interest, being located in an HG (mRNA CR387333) whose genomic organization is similar to that of human and mouse gas5 (growth arrest-specific transcript 5) , a non-protein-coding snoRNA HG. Both contain short exons (< 100 nt) and tandem array snoRNA counterparts in the corresponding intronic regions (Figure 4A). With the exception of the chicken counterpart of human U77 (in intron 4) that is replaced by two cognate snoRNAs of human U80 (GGgCD36a and GGgCD36b in introns 4 and 5 respectively), the other snoRNA counterparts array in the same order on their corresponding HGs of the chicken and human. The gene encoding mRNA CR387333 has a typical TATA box (TATATAA) and a 5'TOP sequence (TCTgCCTTTCCgCCCCT) at position -14 to +3, indicating that it is also a member of the 5' TOP gene family . Comparison of the two transcripts (mRNA CR387333 and gas5) revealed that their sequence similarity varied significantly with each pair of corresponding portions. The most highly conserved regions are the snoRNA sequences, whereas the regions not encoding snoRNA are much more discrepant (42%) between chicken and human. Furthermore, the presence of only short ORFs and numerous stop codons suggests the low probability of protein coding (Figure 4B), which is similar to gas5. Therefore, the gene encoding mRNA CR387333 is classified as a gas5-like non-protein-coding snoRNA HG, and the intron-encoded snoRNAs may be the only functional portions of the transcript.
Structural and functional evolution of chicken snoRNAs
Comparative functional analyses have revealed that many guide snoRNAs involved in rRNA posttranscriptional modification are phylogenetically conserved in mammals. However, extensive recombination and separation of guiding function are also discerned in the chicken snoRNAs when compared with the human snoRNAs. For example, human box C/D snoRNAs U36A and U36C possess the conserved function as guides for the 2'-O-ribose methylation of 18S rRNA-A668 and 28S rRNA-A3703 , respectively, whereas both of the corresponding sites in the chicken rRNAs are predicted to be guided by the single snoRNA GGgCD18a, as well as its paralog GGgCD18b (Figure 5A). A similar case of recombination can be found in the box H/ACA snoRNA GGgACA29 which is predicted to guide the Ψs at 28S rRNA-Ψ1612 and snRNA U5-Ψ43, while the two corresponding modifications are reported to be independently guided by human snoRNAs ACA56 and ACA57, respectively [16, 27, 28]. In contrast, the human snoRNA U32A/B and U69 are double guiders and potentially guide two methylations and two Ψs respectively, whereas the four corresponding modification functions are separately possessed by four single guiders in the chicken (Figure 5B).
It is worth noting that the cooperative evolution of the sequences and guiding function sites can be found in some of the chicken snoRNAs. A case in point is the GGgACA30 and its human cognate E3 snoRNA. Although high sequence similarity (~75%) was found between the two counterparts, GGgACA30 potentially guides the 28S rRNA-Ψ1612, -Ψ3875 and -Ψ3976, whereas E3 is predicted to simply guide the 28S rRNA-Ψ4390 (corresponding to chicken 28S rRNA-Ψ3875) . The analysis of functional loci has revealed that loss of the other two guiding functions for E3 results from extensive nucleotide substitutions and indels, disrupting the snoRNA/substrate RNA base-pairing potential of the 5' and 3' regions in the snoRNA. In contrast, the chicken cognate presents a perfect Ψ-guiding domain and functional motif. Another example is GGgACA14 as a guider for the Ψs of the U12-19 and 18S rRNA-775, whereas its human counterpart ACA22 is predicted to guide the Ψs at different residues of rRNA (28S rRNA-Ψ4966  and -Ψ4975 ). In addition, nucleotide variation can also destroy the guiding function and results in the conversion of guide snoRNAs to orphan snoRNAs, such as GGoACA3 versus the human cognate ACA8 snoRNA.
A model for the prediction of functional snoRNAs in the genome
Comparative analysis of genomes between closely or distantly related species might provide limited information on conserved regions . The chicken bridges the evolutionary gap between mammals and other vertebrates and represents an intermediate-level comparison for the human, making it very useful for detecting functional elements with high specificity . As a family of the most abundant and important noncoding RNAs, many snoRNAs were found to be conserved in different organisms. In this study, we have provided a detailed catalog of chicken snoRNAs to understand snoRNA gene repertoire differences between avian and other vertebrate lineages. A total of 201 gene variants encoding 93 box C/D and 62 box H/ACA snoRNAs were identified in the chicken genome. In contrast to extensive functional and/or pseudogene paralogs found in mammals [17, 31], the majority of snoRNAs in chicken remain singletons, whereas some other novel paralogs are produced by duplication and seem specific to the chicken lineage. Notably, clustered snoRNAs show a large degree of conserved synteny between chicken and human, which greatly simplifies classification of the chicken orthologs of human snoRNAs. Unlike the eutherian mammals with tandem repeated snoRNAs in the imprinted regions [32, 33], none of the imprinted snoRNA orthologs are found in the chicken. Our result might provide some clues to support that imprinting seems to have evolved in therians and hence has only been confirmed in marsupials and eutherian mammals [34–36]. Intriguingly, most of the imprinted snoRNAs in mammals are orphan snoRNAs and constitute new members of the 'dark matter' in the RNA world. With the exception of snoRNA HBII-52 that is reported to be involved in RNA editing  and alternative splicing , the function of the other orphan snoRNAs in mammals remains enigmatic. In our study, we also identified 20 orphan snoRNAs that are expressed in chicken embryos and conserved in different vertebrate species. As with microRNAs, the conservation of snoRNAs among species suggests that they bear conserved biological functions. These orphan snoRNAs might therefore be subject to purifying selection and hence are predicted to be functional in still unknown biological processes.
Intragenic duplication and divergence dominated in the expansion of novel snoRNAs in the chicken
Diverse molecular mechanisms are involved in the creation of new gene (protein-coding gene) structures, such as gene duplication and retroposition . Compared with protein-coding genes, little is known about the creation of novel npcRNA genes in the genome. To our knowledge, two main strategies are responsible for the generation of most of the novel snoRNA paralogs in vertebrates. The first strategy is that some snoRNA paralogs are generated by intragenic duplication where the snoRNAs are tandemly duplicated within the same gene, a process termed cis-duplication [12, 13]. In many cases, the sequence and the secondary structure of these snoRNA paralogs are highly conserved. The second strategy is where snoRNAs may duplicate and insert into a new host gene or a paralogous host gene in a different genomic location, and is termed trans-duplication [12, 13]. In this study, we observed that almost all of the chicken snoRNA paralogs are generated via intragenic duplication, and similar cases are also found in other vertebrates including platypus snoRNAs . There is only one case of five GGgCD61 paralogs created by the combination of intragenic and intergenic duplication during the chicken snoRNA expansion. Most of the novel snoRNA paralogs (~18%) emerge in the adjacent intron regions where the corresponding loci in other mammals lack the cognate snoRNAs, which is in contrast to the few snoRNAs (~7%) found in platypus . The sequences of many snoRNA paralogs have undergone extensive nucleotide variation, but the guiding function regions and the conserved structures are maintained, indicating that strong purifying selection is acting upon them. However, the nucleotide substitutions and insertions or deletions (indels) might also endow the snoRNAs with novel guiding functions (such as GGgACA30 and GGgACA14) or destroy the guiding function (such as GGoACA3), depending on whether or not a perfect functional domain is present within the sequence variation. Recently, hundreds of snoRNAs derived from non-autonomous retroposition have been reported in the human [10, 11] and platypus genomes , revealing a new dimension in the evolution of novel snoRNAs. However, there is no evident trace of the snoRNA-retroposon-like counterparts found in the chicken genome (data not shown), which is consistent with the paucity of functional genes formed by retroposition . Therefore, based on both our results and information obtained from other vertebrate genome analyses, snoRNAs derived from retroposition may originate from mammals. Intragenic duplication and divergence might be the major driving force responsible for expansion of novel snoRNAs in the chicken genome.
A hotspot of snoRNA gene expansion in the gas5-like non-protein-coding snoRNA HG
Non-protein-coding HGs encoding snoRNAs are unusual because they do not appear to specify protein products and snoRNAs may be the only functional portions of the transcripts . Different members of the non-protein-coding HGs encoding snoRNAs (UHG, gas5, U17HG, U19HG, U50HG) have been reported in humans and mice [24, 40–43], and even in fruit flies . Non-protein-coding HGs have not been found in Caenorhabditis elegans, which suggests a different strategy adopted in higher metazoa for regulating snoRNA expression . Among all the members of non-protein-coding HGs, gas5 is of interest because of its large snoRNA-coding capability and deviant accumulation in cells undergoing serum starvation or density arrest . In this study, we detect the gas5-like HG which also lacks protein-coding potential but instead encodes 10 box C/D snoRNAs within its introns. Unexpectedly, we could not find any other non-protein-coding HGs and corresponding snoRNAs in the chicken genome. Comparative analyses of the snoRNA content in the gas5 and other gas5-like HGs reveal diverse snoRNA numbers and gene orders in different vertebrates. In this study, nine human cognate snoRNAs were identified in the chicken gas5-like HG. However, the human U77 counterpart within the fourth intron is replaced by its cognate U80 snoRNA (GGgCD36a), and the other U80 paralog (GGgCD36b) is tandemly arrayed in the next intron, which leads to a different gene order compared to human (Figure 4A). Our comparative analyses of the snoRNAs in gas5 of vertebrates indicted that the mammalian U77 snoRNA might evolve from the other non-mammalian vertebrate cognate U80 snoRNA which has undergone mutation and lost the function of guiding adenine methylation at the corresponding site of 28S rRNA. In addition, a novel snoRNA-like sequence is also detected in intron 6 and found to be chicken-specific, suggesting the process of snoRNA expansion is actively ongoing. Intriguingly, similar cases of extensive snoRNA gene tandem duplication and intragenic transposition can be detected in the corresponding HGs of other vertebrates, such as Danio rerio and Xenopus tropicalis (data not shown). Although the mechanism by which snoRNA sequences become embedded in the introns of their HGs is still enigmatic, the gas5-like HGs appears to be a hotspot for gaining snoRNA novelties in vertebrates.
This is the first genome-wide and systematic screen for snoRNAs in the chicken by applying a computational package and experimental methods. The characteristics of the guiding function and genomic organization of the chicken snoRNAs were extensively compared with that of the human counterparts. We have provided a detailed catalog of chicken snoRNAs to understand snoRNA gene repertoire differences between avian and other vertebrate lineages. Our results improve annotation of the 'darkness matter' in the npcRNA world of the vertebrate genome and provide a unique perspective into snoRNA evolution in vertebrates.
Ten chicken-vertebrate pairwise alignments (galGal3/hg18, galGal3/mm8, galGal3/rn4, galGal3/equCab1, galGal3/monDom4, galGal3/ornAna1, galGal3/anoCar1, galGal3/xenTro2, galGal3/danRer4 and galGal3/fr2) whole-genome alignment (WGA) sequences and the Zebra finch (Taeniopygia guttata) sequence data (taeGut1) were downloaded from the UCSC Genome Bioinformatics site http://genome.ucsc.edu. The repeat families were removed by RepeatMasker http://www.repeatmasker.org/PreMaskedGenomes.html. Sequences and annotation data for known human snoRNA genes (which were used in program training) were downloaded from snoRNA-LBME-db on March 2008 as references for chicken snoRNAs http://www-snorna.biotoul.fr/. UCSC KnownGene, RefGene, Genscan and Ensembl annotation for chicken protein genes and transcript units were downloaded from the UCSC Genome Bioinformatics site and Ensembl Genome Browser http://www.ensembl.org/.
The chicken 28S and 18S rRNA gene sequences (see Additional files 7) were obtained by combining experimental and bioinformatical approaches. The 5.8S rRNA and snRNA (U1-U6, and U12) sequences were retrieved from the UCSC Genome Bioinformatics Site based on similarity search. The 2'-O-methylation and pseudouridylation sites of human rRNAs and snRNAs were mapped to the chicken rRNA and snRNA sequences.
snoRNA searching strategies and sequence analyses
We systematically searched the chicken genome for snoRNAs with 10 whole-genome alignments (WGA) using snoSeeker http://genelab.zsu.edu.cn/snoseeker/ as described previously . For box C/D snoRNAs, the 15 nt flanking sequences of known snoRNA were extracted for folding the terminal stem and training (the CDseeker program extends 5' and 3' stems of the snoRNA sequences by 15 nt). For box H/ACA snoRNAs, the maximum tail length was 3 nt downstream of the ACA box. The potential targets (2'-O-methylation and Ψ) in rRNAs and snRNAs were also determined using snoSeeker. For the comparative analyses of vertebrate snoRNAs, we used the BLAT and the convert track of the UCSC Genome Brower. The human snoRNA-HGs are used as the reference for retrieving the chicken counterparts using the UCSC convert track. Next, the introns of chicken HGs were extracted as a dataset for searching the snoRNA candidates using snoSeeker. The synteny analysis of snoRNA genes in the chicken and human genomes was also implemented on the UCSC Genome Brower.
RNA isolation and construction of snoRNA cDNA libraries
Total cellular RNA was isolated from stage HH34 chicken embryos by the guanidine thiocyanate/phenol-chloroform procedure described by Chomoczynski et al. . The construction of two families of snoRNA cDNA libraries was performed as described previously  with little revision. Briefly, 10 μg total RNA was polyadenylated using a poly(A) polymerase (Takara) at 30°C for 20 min and subsequently reverse-transcribed into the first-stranded cDNAs using [γ-32P]dATP labeled anchor primer dT16-TGT (for box H/ACA snoRNA) and dT16-TCAG (for box CD snoRNA) and MMLV reverse transcriptase (Promega) at 42°C for 1 hr. The reaction mixture was size-fractioned on a denaturing 10% polyacrylamide gel (8 M urea and 1 × TBE buffer). cDNAs with sizes ranging from 90 to 180 nt (for box CD snoRNA) and from 140 to 170 nt (for box H/ACA snoRNA) were excised and eluted from the gel in 0.3 M NaCl buffer. The selected cDNA was tailed with dGTP at the 3' end by using terminal deoxynucleotidyl transferase (Takara) at 37°C for 30 min. G-tailed cDNAs were then amplified by PCR with a forward primer dT23H2 and a reverse primer polyCM. A snoRNA-specific primer and a universal reverse primer polyCM were used for amplification of each of the snoRNAs. The routine PCR protocol (94°C for 5 min, then 30 cycles of 94°C for 1 min, 50°C for 1 min, 72°C for 1 min, and 72°C for 5 min) was used. A 5 μl sample of each PCR product was analyzed on a 2.5% agarose gel.
Oligonucleotides for PCR-based detection and primers used for PCR of the chicken 18S and 28S rRNA genes were synthesized by Invitrogen Co. (Shanghai, China) and are shown in Additional file 8. The primers used in the reverse transcription reaction were 5' end-labeled with [γ-32P]ATP (Yahui Co.) and subjected to purification according to standard laboratory protocols.
Small nucleolar RNAs
small Cajal body-specific RNA
spliceosomal nuclear RNA
reversed transcript PCR
expressed sequence tag
insertion or deletion
open reading frame
Kiss T: Small nucleolar RNAs: an abundant group of noncoding RNAs with diverse cellular functions. Cell. 2002, 109 (2): 145-148. 10.1016/S0092-8674(02)00718-3.
Matera AG, Terns RM, Terns MP: Non-coding RNAs: lessons from the small nuclear and small nucleolar RNAs. Nat Rev Mol Cell Biol. 2007, 8 (3): 209-220. 10.1038/nrm2124.
Darzacq X, Jady BE, Verheggen C, Kiss AM, Bertrand E, Kiss T: Cajal body-specific small nuclear RNAs: a novel class of 2'-O-methylation and pseudouridylation guide RNAs. Embo J. 2002, 21 (11): 2746-2756. 10.1093/emboj/21.11.2746.
Lestrade L, Weber MJ: snoRNA-LBME-db, a comprehensive database of human H/ACA and C/D box snoRNAs. Nucleic Acids Res. 2006, D158-162. 10.1093/nar/gkj002. 34 Database
Bachellerie JP, Cavaille J, Huttenhofer A: The expanding snoRNA world. Biochimie. 2002, 84 (8): 775-790. 10.1016/S0300-9084(02)01402-5.
Cavaille J, Buiting K, Kiefmann M, Lalande M, Brannan CI, Horsthemke B, Bachellerie JP, Brosius J, Huttenhofer A: Identification of brain-specific and imprinted small nucleolar RNA genes exhibiting an unusual genomic organization. Proc Natl Acad Sci USA. 2000, 97 (26): 14311-14316. 10.1073/pnas.250426397.
Yang JH, Zhang XC, Huang ZP, Zhou H, Huang MB, Zhang S, Chen YQ, Qu LH: snoSeeker: an advanced computational package for screening of guide and orphan snoRNA genes in the human genome. Nucleic Acids Res. 2006, 34 (18): 5112-5123. 10.1093/nar/gkl672.
Brown JW, Echeverria M, Qu LH: Plant snoRNAs: functional evolution and new modes of gene expression. Trends Plant Sci. 2003, 8 (1): 42-49. 10.1016/S1360-1385(02)00007-9.
Smith CM, Steitz JA: Sno storm in the nucleolus: new roles for myriad small RNPs. Cell. 1997, 89 (5): 669-672. 10.1016/S0092-8674(00)80247-0.
Weber MJ: Mammalian small nucleolar RNAs are mobile genetic elements. PLoS Genet. 2006, 2 (12): e205-10.1371/journal.pgen.0020205.
Luo Y, Li S: Genome-wide analyses of retrogenes derived from the human box H/ACA snoRNAs. Nucleic Acids Res. 2007, 35 (2): 559-571.
Zemann A, op de Bekke A, Kiefmann M, Brosius J, Schmitz J: Evolution of small nucleolar RNAs in nematodes. Nucleic Acids Res. 2006, 34 (9): 2676-2685. 10.1093/nar/gkl359.
Schmitz J, Zemann A, Churakov G, Kuhl H, Grutzner F, Reinhardt R, Brosius J: Retroposed SNOfall-a mammalian-wide comparison of platypus snoRNAs. Genome Res. 2008, 18 (6): 1005-1010. 10.1101/gr.7177908.
Griffiths-Jones S: Annotating noncoding RNA genes. Annu Rev Genomics Hum Genet. 2007, 8: 279-298. 10.1146/annurev.genom.8.080706.092419.
Huttenhofer A, Kiefmann M, Meier-Ewert S, O'Brien J, Lehrach H, Bachellerie JP, Brosius J: RNomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger RNAs in mouse. Embo J. 2001, 20 (11): 2943-2953. 10.1093/emboj/20.11.2943.
Kiss AM, Jady BE, Bertrand E, Kiss T: Human box H/ACA pseudouridylation guide RNA machinery. Mol Cell Biol. 2004, 24 (13): 5797-5807. 10.1128/MCB.24.13.5797-5807.2004.
Schattner P, Barberan-Soler S, Lowe TM: A computational screen for mammalian pseudouridylation guide H/ACA RNAs. Rna. 2006, 12 (1): 15-25. 10.1261/rna.2210406.
Gu AD, Zhou H, Yu CH, Qu LH: A novel experimental approach for systematic identification of box H/ACA snoRNAs from eukaryotes. Nucleic Acids Res. 2005, 33 (22): e194-10.1093/nar/gni185.
Consortium ICGS: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432 (7018): 695-716. 10.1038/nature03154.
Kiss AM, Jady BE, Darzacq X, Verheggen C, Bertrand E, Kiss T: A Cajal body-specific pseudouridylation guide RNA is composed of two box H/ACA snoRNA-like domains. Nucleic Acids Res. 2002, 30 (21): 4643-4649. 10.1093/nar/gkf592.
Ro S, Park C, Jin J, Sanders KM, Yan W: A PCR-based method for detection and quantification of small RNAs. Biochem Biophys Res Commun. 2006, 351 (3): 756-763. 10.1016/j.bbrc.2006.10.105.
Ro S, Park C, Sanders KM, McCarrey JR, Yan W: Cloning and expression profiling of testis-expressed microRNAs. Dev Biol. 2007, 311 (2): 592-602. 10.1016/j.ydbio.2007.09.009.
Brown JW, Marshall DF, Echeverria M: Intronic noncoding RNAs and splicing. Trends Plant Sci. 2008, 13 (7): 335-342. 10.1016/j.tplants.2008.04.010.
Smith CM, Steitz JA: Classification of gas5 as a multi-small-nucleolar-RNA (snoRNA) host gene and a member of the 5'-terminal oligopyrimidine gene family reveals common features of snoRNA host genes. Mol Cell Biol. 1998, 18 (12): 6897-6909.
Amaldi F, Pierandrei-Amaldi P: TOP genes: a translationally controlled class of genes including those coding for ribosomal proteins. Prog Mol Subcell Biol. 1997, 18: 1-17.
Nicoloso M, Qu LH, Michot B, Bachellerie JP: Intron-encoded, antisense small nucleolar RNAs: the characterization of nine novel species points to their direct role as guides for the 2'-O-ribose methylation of rRNAs. J Mol Biol. 1996, 260 (2): 178-195. 10.1006/jmbi.1996.0391.
Ofengand J, Bakin A: Mapping to nucleotide resolution of pseudouridine residues in large subunit ribosomal RNAs from representative eukaryotes, prokaryotes, archaebacteria, mitochondria and chloroplasts. J Mol Biol. 1997, 266 (2): 246-268. 10.1006/jmbi.1996.0737.
Krol A, Gallinaro H, Lazar E, Jacob M, Branlant C: The nuclear 5S RNAs from chicken, rat and man. U5 RNAs are encoded by multiple genes. Nucleic Acids Res. 1981, 9 (4): 769-787. 10.1093/nar/9.4.769.
Ganot P, Bortolin ML, Kiss T: Site-specific pseudouridine formation in preribosomal RNA is guided by small nucleolar RNAs. Cell. 1997, 89 (5): 799-809. 10.1016/S0092-8674(00)80263-9.
Furlong RF: Insights into vertebrate evolution from the chicken genome sequence. Genome Biol. 2005, 6 (2): 207-10.1186/gb-2005-6-2-207.
Schmitz J, Zemann A, Churakov G, Kuhl H, Grutzner F, Reinhardt R, Brosius J: Retroposed SNOfall – a mammalian-wide comparison of platypus snoRNAs. Genome Res. 2008, 18 (6): 1005-1010. 10.1101/gr.7177908.
Seitz H, Royo H, Lin SP, Youngson N, Ferguson-Smith AC, Cavaille J: Imprinted small RNA genes. Biol Chem. 2004, 385 (10): 905-911. 10.1515/BC.2004.118.
Cavaille J, Seitz H, Paulsen M, Ferguson-Smith AC, Bachellerie JP: Identification of tandemly-repeated C/D snoRNA genes at the imprinted human 14q32 domain reminiscent of those at the Prader-Willi/Angelman syndrome region. Hum Mol Genet. 2002, 11 (13): 1527-1538. 10.1093/hmg/11.13.1527.
Wilkins JF, Haig D: What good is genomic imprinting: the function of parent-specific gene expression. Nat Rev Genet. 2003, 4 (5): 359-368. 10.1038/nrg1062.
Edwards CA, Rens W, Clarke O, Mungall AJ, Hore T, Graves JA, Dunham I, Ferguson-Smith AC, Ferguson-Smith MA: The evolution of imprinting: chromosomal mapping of orthologues of mammalian imprinted domains in monotreme and marsupial mammals. BMC Evol Biol. 2007, 7: 157-10.1186/1471-2148-7-157.
Warren WC, Hillier LW, Marshall Graves JA, Birney E, Ponting CP, Grutzner F, Belov K, Miller W, Clarke L, Chinwalla AT, et al: Genome analysis of the platypus reveals unique signatures of evolution. Nature. 2008, 453 (7192): 175-183. 10.1038/nature06936.
Vitali P, Basyuk E, Le Meur E, Bertrand E, Muscatelli F, Cavaille J, Huttenhofer A: ADAR2-mediated editing of RNA substrates in the nucleolus is inhibited by C/D small nucleolar RNAs. J Cell Biol. 2005, 169 (5): 745-753. 10.1083/jcb.200411129.
Kishore S, Stamm S: The snoRNA HBII-52 regulates alternative splicing of the serotonin receptor 2C. Science. 2006, 311 (5758): 230-232. 10.1126/science.1118265.
Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4 (11): 865-875. 10.1038/nrg1204.
Tycowski KT, Shu MD, Steitz JA: A mammalian gene with introns instead of exons generating stable RNA products. Nature. 1996, 379 (6564): 464-466. 10.1038/379464a0.
Pelczar P, Filipowicz W: The host gene for intronic U17 small nucleolar RNAs in mammals has no protein-coding potential and is a member of the 5'-terminal oligopyrimidine gene family. Mol Cell Biol. 1998, 18 (8): 4509-4518.
Bortolin ML, Kiss T: Human U19 intron-encoded snoRNA is processed from a long primary transcript that possesses little potential for protein coding. Rna. 1998, 4 (4): 445-454.
Tanaka R, Satoh H, Moriyama M, Satoh K, Morishita Y, Yoshida S, Watanabe T, Nakamura Y, Mori S: Intronic U50 small-nucleolar-RNA (snoRNA) host gene of no protein-coding potential is mapped at the chromosome breakpoint t(3;6)(q27;q15) of human B-cell lymphoma. Genes Cells. 2000, 5 (4): 277-287. 10.1046/j.1365-2443.2000.00325.x.
Huang ZP, Zhou H, He HL, Chen CL, Liang D, Qu LH: Genome-wide analyses of two families of snoRNA genes from Drosophila melanogaster, demonstrating the extensive utilization of introns for coding of snoRNAs. Rna. 2005, 11 (8): 1303-1316. 10.1261/rna.2380905.
Huang ZP, Chen CJ, Zhou H, Li BB, Qu LH: A combined computational and experimental analysis of two families of snoRNA genes from Caenorhabditis elegans, revealing the expression and evolution pattern of snoRNAs in nematodes. Genomics. 2007, 89 (4): 490-501. 10.1016/j.ygeno.2006.12.002.
Chomczynski P, Sacchi N: The single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction: twenty-something years on. Nat Protoc. 2006, 1 (2): 581-585. 10.1038/nprot.2006.83.
This research was supported by the National Basic Research Program (No. 2005CB724600), the funds from the Ministry of Education of China and Guangdong Province (No. IRT0447, NSF05200303) and the National Natural Science Foundation of China (30570398, 30771151, 30830066). We thank three anonymous referees who made valuable suggestions that led to considerable improvements in the paper.
PS and LHQ conceived the study and contributed to manuscript writing. JHY and PS collected the data and carried out the data analyses. PS performed the experiments. HZ and DGG assisted in experimental design and the data analysis, respectively. All authors read and approved the final manuscript.
Peng Shao, Jian-Hua Yang contributed equally to this work.
Electronic supplementary material
About this article
Cite this article
Shao, P., Yang, J., Zhou, H. et al. Genome-wide analysis of chicken snoRNAs provides unique implications for the evolution of vertebrate snoRNAs. BMC Genomics 10, 86 (2009). https://doi.org/10.1186/1471-2164-10-86
- Chicken Genome
- Vertebrate Genome
- snoRNA Gene
- Intragenic Duplication
- Detailed Catalog