| Nucleic Acids Research | Pages |
Comparative structure analysis of vertebrate ribonuclease P RNA
Introduction
Materials And Methods
DNA isolation
PCR amplification of RNase P RNA genes
Cloning of PCR products
Plasmid isolation and sequencing
Sequence alignment, design of secondary structure and covariation
Results
Discussion
Acknowledgements
References
Comparative structure analysis of vertebrate ribonuclease P RNA
DDBJ/EMBL/GenBank accession nos+
ABSTRACT
INTRODUCTION
Ribonuclease P (RNase P) is the ubiquitous endonuclease that carries out maturation of the 5[prime] end of tRNA (see 1 and 2 for reviews). All cellular versions of RNase P from each of the three domains of life-Bacteria, Eucarya and Archaea-are ribonucleoproteins, consisting of one RNA and at least one protein. The RNA component of bacterial RNase P is catalytically active in the absence of protein in vitro (3), but the archaeal and eucaryal RNase P RNAs are not independently active. Although RNase P RNAs from diverse organisms exhibit extensive variation in sequence and structure, all contain five universally conserved regions (4) that are considered to reflect similar core structure in the RNAs and possibly homologous function.
The global structure of the bacterial type of RNase P RNA is increasingly well-characterized and there is considerable information on the archaeal counterpart; however, less is known about the eucaryotic version of the RNA. Accumulation of sequences for comparative structure analysis of the eucaryal RNase P RNA has focused mainly on yeasts (5,6). Animal RNAs for comparison have so far been represented by those from only three types of animals: mammals, one frog and one fish (7-10). Thus, the diversity of sequences available for comparative structure analysis has been insufficient to specify a sequence alignment and a basal structure for the vertebrate RNase P RNA.
To expand the diversity of vertebrate sequences available for comparison, we determined 13 partial sequences of RNase P RNAs from 11 of the main vertebrate lineages. We used PCR to isolate the main portions of the RNase P RNA genes from five frogs (Ascaphus truei, Alytes obstetricans, Bufo bufo, Discoglossus jeanneae and Pelobates cultripes), two salamanders (Ensatina eschscholtzii and Euproctus asper), two lizards (Phrynosoma douglasi and Uta stansburiana), one snake (Lachesis muta) and one fish (Oncorhynchus sp.). Although previously examined, the corresponding partial genes from Homo sapiens, Xenopus laevis and Danio rerio also were obtained by PCR and re-examined. We aligned the sequences consistently with the proposed universally conserved core of the RNA (4), and inferred a minimum-consensus secondary structure model of the vertebrate RNase P RNA.
MATERIALS AND METHODS
DNA isolation
DNA from frogs and salamanders was extracted as described previously (11). For the three reptiles, 20 mg of frozen tissue from each species was incubated overnight at 55°C in 1 ml of 100 mM EDTA, 10 mM Tris-HCl (pH 8), 10 mM NaCl, 1% SDS, 8 mg dithiothreitol and 2 mg proteinase K, and DNA was purified by phenol-chloroform extraction (12). Genomic DNAs were precipitated from 0.3 M NaCl with 2.5 vol ethanol and resuspended in 200 µl of 10 mM Tris-HCl (pH 8) and 1 mM EDTA. In the case of H.sapiens, DNA was isolated from 80 µl of whole blood treated with heparin (20 U/ml) as described earlier (13). Salmon testis DNA was purchased from Sigma.
PCR amplification of RNase P RNA genes
Each 100 µl PCR reaction contained 8 U of Taq DNA polymerase (Promega), 100 ng of each primer (below), 300 ng genomic DNA, 0.2 mM of each dNTP, 1.5 mM MgCl2, 50 mM KCl, 10 mM Tris-HCl (pH 9 at 25°C) and 0.1% Triton X-100. All PCR reactions started with an initial denaturation step for 2 min at 94°C, followed by 20-25 amplification cycles (below), and a final extension step for 10 min at 72°C. To avoid non-specific PCR products, we chose the highest annealing temperature that resulted in PCR fragments of the expected size. PCR products were visualized by agarose gel electrophoresis. In case of low yield or if more than one PCR product was obtained, a fragment (~5 µl) of agarose containing the desired DNA fragments was excised from the gel and used in a second PCR amplification. The following primers (nucleotide positions correspond to the individual RNase P RNA) and cycle profiles (denaturation-primer annealing-extension) were used for the respective taxa: X.laevis, primers XEN-F (positions 1-33, 5[prime]-AGATGTGGAGGGAAGCTCAGTGCTGAGGCTACA-3[prime]) and XEN-R (positions 273-241, 5[prime]-CTCAGCCATTGAGCTTGAATGTGCCAGCTGTGG-3[prime]), 20 cycles (2 min at 92°C, 45 s at 72°C). The profile for reamplification was 20 cycles (1 min at 92°C, 45 s at 67°C, 30 s at 72°C); D.rerio, primers Euk-Z1 (positions 6-30, 5[prime]-GCGGAAGGAAGCTCACTGTAGAGG-3[prime]) and Euk-Z2 (positions 254-277, 5[prime]-CTAGCATGTTATGCTAGCTCACAC-3[prime]), 25 cycles (1.5 min at 92°C, 45 s at 55°C, 45 s at 72°C). Reamplification was performed under the same conditions, except the annealing temperature was 50°C; H.sapiens, primers Euk-D2 (positions 4-36, 5[prime]-GGGCGGAGGGAAGCTCATCAGTGGGGCCACGAG-3[prime]) and HeLa-A (positions 336-304, 5[prime]-GCGGAGGAGAGTGGTCTGAATTGGGTTATGAGG-3[prime]), 25 cycles (2 min at 92°C, 45 s at 72°C); Oncorhynchus sp., primers Euk-D (7) (positions 1-26, 5[prime]-GGGCGGAGGGAAGCTCATCAGTGGGG-3[prime]) and Euk-A (7) (positions 308-278 and 312-282, 5[prime]-CGGAATTCCGGAGAGTGGTCTGAATTGGGTT-3[prime]), 20 cycles (1.5 min at 92°C, 45 s at 50°C, 1 min at 72°C). Reamplification was performed under the same conditions; A.truei, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 320-290), 25 cycles (2 min at 92°C, 1 min at 47°C, 45 s at 72°C). Reamplification was performed under the same conditions, except the annealing temperature was 49°C; A.obstetricans, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 280-310), 25 cycles (2 min at 92°C, 1 min at 49°C, 45 s at 72°C); B.bufo, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 391-361), 25 cycles (2 min at 92°C, 1 min at 48°C, 45 s at 72°C); D.jeanneae, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 406-376), 25 cycles (2 min at 92°C, 1 min at 48°C, 45 s at 72°C); E.asper, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 303-273), 25 cycles (2 min at 92°C, 1 min at 48°C, 45 s at 72°C). Reamplification was performed under the same conditions; E.eschscholtzii, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 326-296), 25 cycles (2 min at 92°C, 1 min at 50°C, 45 s at 72°C); L.muta, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 330-300 and 348-318), 25 cycles (2 min at 92°C, 1 min at 49°C, 45 s at 72°C); P.cultripes, primers Euk-D2 (positions 1-33) and Euk-A (7) (positions 310-280), 25 cycles (2 min at 92°C, 1 min at 49°C, 45 s at 72°C); P.douglasi, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 331-301), 25 cycles (2 min at 92°C, 1 min at 50°C, 45 s at 72°C); U.stansburiana, primers Euk-Z1 (positions 1-24) and Euk-A (7) (positions 328-298), 25 cycles (2 min at 92°C, 1 min at 49°C, 45 s at 72°C).
In case of sequence ambiguities or heterogeneity, the PCR reaction was repeated under the same conditions, except that 2.5 U of Pfu DNA polymerase (Stratagene) and 1× PCR buffer (Stratagene) were used.
Cloning of PCR products
PCR fragments were cloned directly using the pGEMr-T Vector System (Promega), the TA Cloningr Kit (Invitrogen) or the TOPO TA Cloningr Kit (Invitrogen), following the protocols provided by the manufacturers. PCR products derived from Pfu DNA polymerase were blunt-end ligated into pGEMr5Zf(+) vector for 1 h at 37°C in 10 µl containing 50 ng of EcoRV digested pGEMr5Zf(+), 70 ng of PCR product, 1 Weiss U T4 DNA ligase (Promega) and 1× ligase buffer (Promega). In the case of the pGEMr-T Vector and pGEMr5Zf(+) vector, each ligation reaction was purified over a Chromabond-100 TE spin column (Clontech) prior to the transformation of Escherichia coli DH5[alpha] by electroporation with an E.coli Pulser Transformation Apparatus (Bio-Rad). Preparation of the E.coli DH5[alpha] cells and electroporation were performed as described (14) using 20 ng of recombinant plasmid DNA per 50 µl of cell suspension. Recombinants were selected by blue/white color of colonies.
Plasmid isolation and sequencing
Plasmid isolation was performed using the QIAprep[trade] Plasmid Kit (Qiagen). Recombinant plasmid DNA (200 ng) was sequenced by cycle sequencing (24 cycles: 30 s at 96°C; 15 s at 45°C, 4 min at 60°C) with AmpliTaqrFS (Perkin Elmer) using 20 ng primer T7 (5[prime]-GTAATACGACTCACTATAGGG-3[prime]) and primer M13-rev (5[prime]-AGGAAACAGCTATGA-3[prime]) or primer SP6 (5[prime]-ATTTAGGTGACACTATAG-3[prime]). Sequences were analyzed on an ABI 377 or ABI 373A automated DNA sequencer (Applied Biosystems) according to the manufacturer's protocols. Sequences were deposited in GenBank under the accession numbers AF004373 (A.truei), AF004374 (A.obstetricans), AF044737 (B.bufo), AF004372 (D.jeanneae), AF004375 (E.asper), AF044328(E.eschscholtzii), AF004376 (L.muta), AF044326 (L.muta with insertion), AF004377 (Oncorhynchus sp.), AF044331 (Oncorhynchus sp. with deletion), AF004371 (P.cultripes), AF044327 (P.douglasi), AF044325 (U.stansburiana) and AF044330 (X.laevis).
Sequence alignment, design of secondary structure and covariation
All available RNase P RNA sequences from vertebrates were aligned manually using the GDE sequence editor (15). Comparative analysis of the secondary structures by covariation was performed as described earlier (16-18). The nomenclature of structural elements in eucaryal RNase P follows that described for bacterial RNase P RNA (4,19).
RESULTS
As detailed in Materials and Methods, we amplified, cloned and sequenced 13 partial (~85% of full-length) RNase P RNA genes from 11 animals, more than doubling the previously known number of vertebrate sequences for this RNA and substantially expanding the breadth of vertebrate diversity sampled for comparative analysis of the RNA. Figure
Figure 1. Diagrammatic representation of phylogenetic relationships of vertebrates with available RNase P RNA sequences. Relationships are based on morphological and molecular characters among extant vertebrates (22,29). Taxa were chosen for this study (text) to span a broad extent of vertebrate lineages and thereby delineate the core secondary structure of vertebrate RNase P RNA. References to previously determined sequences are shown as superscripts. Asterisks indicate sequences determined or confirmed in this study. Figure 2. Alignment of RNase P RNAs including 16 partial sequences investigated or confirmed in this study. Brackets above the alignment indicate helices in the corresponding RNA structures. Sequences derived from PCR primers are in lower case letters. Regions that could not be aligned with confidence are in italics. Bold letters in brackets below the alignment indicate universally conserved nucleotides that define the conserved regions (CR) I-V. The sequences in the variable region P12 could not be aligned so were omitted from the alignment and represented by the total number of nucleotides. Abbreviations are: Hum, H.sapiens (8,*); Mou, Mus musculus (7); Dis, D.jeanneae; Aly, A.obstretricans; Buf, B.bufo (*); Lac, L.muta (*); Uta, U.stansburiana (*); Phr, P.douglasi (*); Ens, E.eschscholtzii (*); Eup, E.asper (*); Zeb, D.rerio (10,*); Xen, X.laevis (9,*); Pel, P.cultripes (*); Sal, Onchorhynchus sp. (*); Asc, A.truei (*). (Asterisks indicate organisms for which sequences were determined or confirmed by this study.) The structures of the vertebrate RNase P RNAs conform to the general eucaryotic model (4) and were consistent with numerous covariations in the amplified regions of the genes. One example of the derived secondary structures of vertebrate RNase P RNA, with a bacterial version of the RNA for comparison, is shown in Figure Figure 3. Comparison of proposed secondary structures of RNase P RNAs of E.coli (bacterial) and D.rerio (vertebrate). The nomenclature for the structural elements of D.rerio RNase P RNA is based on putatively homologous bacterial counterparts (4,19). Homologs to the bacterial helices P5, P6, P13, P14, P15, P16, P17 and P18 are not identifiable in the RNase P RNAs of D.rerio or other vertebrates. In contrast, the region of P3 is more complex in vertebrate and other eucaryotic RNase P RNAs than in the bacterial version. Watson-Crick base pairs are indicated by dashes, and G/U base pairs by dots. Gray shading indicates base pairs confirmed by covariation among the vertebrate RNAs. The structure of helix P12 of the D.rerio RNA shown is based on results with the computer program Mfold (30) and is not confirmed by covariation. Highlighted bases are universally conserved components of conserved regions (CRs) in all RNase P RNAs. Figure 4. Variation in vertebrate RNase P RNA secondary structure. The proposed minimum consensus secondary structure for vertebrate RNase P RNA is shown in the center. The nomenclature for the structure is based on putatively homologous bacterial counterparts (4,19). Different possible structures for P3, P12 and P19 from the organisms identified illustrate size and structural variation among vertebrate RNase P RNAs. Structures for P12 helices, which because of sequence length variation cannot be aligned for detection of covariations, are based on results using the computer program Mfold (30). Watson-Crick base pairs are indicated by dashes, and G/U base pairs by dots. Gray lines indicate regions that bind PCR primers. Highlighted bases are universally conserved regions in all RNase P RNAs. The vertebrate RNAs display all of the conserved regions (CR I-V, Fig. We detected two types of RNase P RNA genes from three of the vertebrates: in two out of nine clones for the Oncorhynchus sp. RNA we observed a deletion of positions 117-120 in P9. We also discovered a triple-repeat of positions 187-193 in three out of 15 clones for the L.muta gene. This results in an extension of helix P12 in that type of RNA. Minor RNAs were designated Oncorhynchus sp. 2 and L.muta 2, respectively. These variant types of RNase P genes were all detected in independent PCR amplifications, and presumably represent genomic diversity in alleles or sequence repeats. The published RNase P RNA sequences of H.sapiens (8) and D.rerio (10) were confirmed. In the case of the X.laevis sequence (9), a few ambiguities were clarified. The detail of the recently proposed secondary structure of helix P3 of human and zebrafish RNase P RNA (4) was adjusted (Fig.
DISCUSSION
Earlier proposed secondary structures for vertebrate RNase P RNAs (4,7,10) were based on consideration of only a few sequences. To provide further examples for comparison, we chose organisms to cover a broad phylogenetic diversity (Fig.
Despite the fact that the vertebrate RNAs lack some structural elements present in the catalytically active bacterial RNA, the new sequences offer no new clues as to why the eucaryal RNase P RNAs are catalytically inactive without the protein components of the native holoenzyme. For instance, the eucaryal RNAs lack the bacterial helix P15/16, an internal loop of which interacts with the 3[prime]-terminal CCA of pre-tRNA (24,25). On the other hand, some structural elements of eucaryal RNase P RNA are more complex than their bacterial and archaeal counterparts. The eucaryal helix P3, for instance, consists of two short helices (P3a and P3b) connected by an internal loop (Fig.
Helix P19 is present in all vertebrate RNase P RNAs examined so far. In contrast, it is sometimes absent from bacterial and archaeal RNAs. The function of helix P19, if any, is not known. Commonly, however, the terminal loops of such ephemeral helices dock elsewhere in the RNA to bolster the superstructure (19).
The numbers of gene copies for RNase P RNA that occur in the genomes of the vertebrates for which sequences have been determined are not known. Multiple homologs of RNase P RNA genes have been detected in genomic DNA of mouse, rat, rabbit and human (26). Three genes analyzed from mouse are truncated within CRIV and therefore do not contain full length RNase P RNA genes. Although the corresponding RNAs were detected (26), it seems unlikely that these truncated genes code for functional RNase P RNAs. Some sequence degeneracies were detected in the present study in multiple PCR reactions. Twenty percent of the RNase P RNA genes obtained from Oncorhynchus sp. show a deletion of positions 117-120 in the region of P9. In the case of L.muta, 20% of the RNase P RNA genes obtained have a triple repeat of nucleotides 187-193 in P12. Those types of repeats also have been discovered recently in helix P12 of several cyanobacterial RNase P RNAs (27). Such sequence idiosyncrasies might indicate the occurrence of more than one gene-copy for RNase P RNA in these vertebrates or represent two alleles. Appropriate Southern hybridization studies have not been conducted to determine copy-numbers and organization of these genes.
Partial or complete RNase P RNA sequences from mammals, reptiles, amphibia and fish provide broad representation of this RNA among vertebrates. We were not able to obtain the corresponding gene from bird (two chickens, Gallus domesticus and Callipela californica, and two hawks, Buteo magnirostis and Zonotrichia atricapilla tested), and from turtle (Clemmysmarmorata and Chelonia mydas tested), presumably because of primer-incompatibility. Since birds are derived from deeply within the reptiles, however, we believe that the available sequences broadly represent the vertebrate RNase P RNA as summarized in the proposed minimum consensus secondary structure (Fig.
ACKNOWLEDGEMENTS
We thank A.Barth and M.Dunaway for the gifts of D.rerio and X.laevis DNAs, respectively. Many thanks to J.W.Brown, G.E.Fox, J.K.Harris and M.A.Tanner for helpful discussions. This work was supported by an NIH grant to N.R.P.
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: www-admin{at}oup.co.uk
Last modification: 2 Jul 1998
Copyright©Oxford University Press, 1998.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
A. Busch and R. Backofen
INFO-RNA--a fast approach to inverse RNA folding
Bioinformatics,
August 1, 2006;
22(15):
1823 - 1831.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
S. M. MARQUEZ, J. K. HARRIS, S. T. KELLEY, J. W. BROWN, S. C. DAWSON, E. C. ROBERTS, and N. R. PACE
Structural implications of novel diversity in eucaryal RNase P RNA
RNA,
May 1, 2005;
11(5):
739 - 751.
[Abstract]
[Full Text]
[PDF]
![]()
![]()
![]()

![]()
![]()
![]()
C. Bardeleben, R. L. Moore, and R. K. Wayne
Isolation and Molecular Evolution of the Selenocysteine tRNA (Cf TRSP) and RNase P RNA (Cf RPPH1) Genes in the Dog Family, Canidae
Mol. Biol. Evol.,
February 1, 2005;
22(2):
347 - 359.
[Abstract]
[Full Text]
[PDF]
![]()
This Article ![]()
![]()
Abstract
![]()
Print PDF (192K)
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (20)
![]()
Request Permissions ![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Pitulle, C.
![]()
Articles by Pace, N. R.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Pitulle, C.
![]()
Articles by Pace, N. R.
![]()
Social Bookmarking ![]()
![]()
What's this?