- Research article
- Open Access
Phylogenetic footprinting of non-coding RNA: hammerhead ribozyme sequences in a satellite DNA family of Dolichopodacave crickets (Orthoptera, Rhaphidophoridae)
© Martinsen et al; licensee BioMed Central Ltd. 2010
- Received: 21 July 2009
- Accepted: 4 January 2010
- Published: 4 January 2010
The great variety in sequence, length, complexity, and abundance of satellite DNA has made it difficult to ascribe any function to this genome component. Recent studies have shown that satellite DNA can be transcribed and be involved in regulation of chromatin structure and gene expression. Some satellite DNAs, such as the pDo500 sequence family in Dolichopoda cave crickets, have a catalytic hammerhead (HH) ribozyme structure and activity embedded within each repeat.
We assessed the phylogenetic footprints of the HH ribozyme within the pDo500 sequences from 38 different populations representing 12 species of Dolichopoda. The HH region was significantly more conserved than the non-hammerhead (NHH) region of the pDo500 repeat. In addition, stems were more conserved than loops. In stems, several compensatory mutations were detected that maintain base pairing. The core region of the HH ribozyme was affected by very few nucleotide substitutions and the cleavage position was altered only once among 198 sequences. RNA folding of the HH sequences revealed that a potentially active HH ribozyme can be found in most of the Dolichopoda populations and species.
The phylogenetic footprints suggest that the HH region of the pDo500 sequence family is selected for function in Dolichopoda cave crickets. However, the functional role of HH ribozymes in eukaryotic organisms is unclear. The possible functions have been related to trans cleavage of an RNA target by a ribonucleoprotein and regulation of gene expression. Whether the HH ribozyme in Dolichopoda is involved in similar functions remains to be investigated. Future studies need to demonstrate how the observed nucleotide changes and evolutionary constraint have affected the catalytic efficiency of the hammerhead.
- Shannon Entropy
- Compensatory Mutation
- Phylogenetic Footprint
- Ribozyme Activity
- satDNA Family
Noncoding tandem-repetitive satellite DNA (satDNA) has long been known to constitute a large portion of any eukaryotic genome . Recent studies in Drosophila through underreplication documented a strong positive correlation between genome size and the amount of satDNA . SatDNA is usually located in the heterochromatic parts of the chromosomes close to the centromeres and telomeres. The high evolutionary turn-over of satDNA and the frequent species-specific patterns have puzzled researchers for many years. Such data challenge the hypothesis that there is a general biological function of satDNA. The enormous diversity of satDNA in nucleotide sequence, length of repeats, complexity, and genomic abundance, instead suggest that a specific function cannot be ascribed to satDNA.
However, transcription of satDNA and other non-coding DNA has been found in many species (e.g. [3–10]). The function of the vast majority of these transcripts is unclear, but recent studies using RNA interference have documented the role of non-coding RNAs in regulation of gene expression, chromatin organization, and genome functioning (for a review, see ). For example, it has been demonstrated that small non-coding RNAs can play a crucial role in regulating heterochromatin formation  and it has been found that transcripts of tandem repeated non-coding DNA can give rise to dsRNAs .
The transcription of satDNA was in some instances also linked to hammerhead (HH) ribozyme activity which may affect certain regulatory mechanisms in the cell. SatDNA-derived HH ribozyme structures have been detected in several organisms as different as schistosomes, cave crickets, and salamanders [14–16]. The HH ribozyme is one of the smallest catalytic RNAs and was first identified in viroid and viroid-like satellite RNA sequences where they catalyze a specific phosphodiester bond isomerization reaction in the course of rolling-circle replication[17, 18]. All HH ribozymes detected in animal satDNA so far have been shown to self-cleave in cis long multimeric transcripts into monomers [6, 14–16, 19, 20]. The HH ribozyme from Schistosoma mansoni is certainly the best characterized of natural HH ribozymes. This particular ribozyme has also been shown to perform efficient ligation  and trans-cleavage of an RNA target . Recently, an active HH ribozyme was also found in a mammalian messenger RNA that self-cleaved both in vitro and in vivo and reduced protein expression of C-type lectin type II (Clec2) genes in mouse cells . This study showed that ribozymes might also act in post-transcriptional regulation of protein coding gene expression in a way that is similar to destabilizing protein factors.
HH ribozymes are characterized by a conserved central core which is mandatory for cleavage activity. The consensus sequence for a HH ribozyme consists of three paired stems that branch from the core as defined by conserved nucleotides (reviewed in ). In addition, the natural HH ribozymes require additional sequence elements outside of the conserved catalytic core to facilitate intracellular activity. These elements are variable and are therefore not phylogenetically conserved . The correlation between structure and function of the HH ribozyme became clear only recently [25–27]. For many years studies of HH ribozymes were restricted to a minimal construct. However, there have been major discrepancies between the crystal structure of this particular construct and the results of biochemical experiments [26, 28, 29]. The biochemical data could only be explained if a large-scale rearrangement of the core region was assumed during catalysis [28, 30, 31]. Recently, the HH ribozymes found in schistosomes were characterized in more detail and loop-loop interactions were shown to have a large impact on HH ribozyme activity [25, 26, 32].
Previously, database searches identified a potential HH ribozyme (49 base pairs) in Dolichopoda cave crickets embedded within the pDo500satDNA family and self-cleavage of satDNA transcripts has been shown in vitro for RNA from D. baccettii . The roughly 500 bp repeats of the pDo500 have been found in all Dolichopoda species studied so far  and constitute approximately 5% of the D. schiavazzii genome . Phylogenetic analyses of the pDo500 sequences from 12 Dolichopoda species showed that this satDNA family evolved gradually at a rate of 3.4% per one million year . Thus, pDo500-based phylogenetic hypotheses were mainly congruent with the mitochondrial DNA phylogeny [33, 35]. The HH ribozyme in Dolichopoda performed self-cleavage in vitro and this was associated with processing of long multimer transcripts into monomers in vivo . For the HH ribozyme of Dolichopoda two cleavage mechanisms were suggested: 1) cis cleavage by a single-HH (sHH) ribozyme which is the folding of the 49 bp HH sequence, and 2) trans cleavage by a double-hammerhead (dHH) which is a hybrid structure between two extended HH sequences of 68 bp.
We took a phylogenetic footprint approach to analyze in detail the pDo500 satDNA family from 38 different populations representing 12 species of Dolichopoda cave crickets with particular emphasis on the HH region. We assessed whether the HH motif in the pDo500 satDNA is more conserved than other parts of the pDo500 repeat unit. HH stems should be more conserved than loops, due to their importance for stabilizing secondary structures. Furthermore, the core region was expected invariant. RNA folding patterns of the HH regions of pDo500 sequences from different populations and species may indicate which Dolichopoda species are expected to utilize active HH ribozymes. The results indicated that the putative HH motif in the pDo500 satDNA is likely to be under selective constraints and may be functional in Dolichopoda.
Mann-Whitney U Tests of the bootstrapped Shannon-Wiener index and the Homozygosity index (Approach 1 see Methods).
Number of sites compared
Shannon-Wiener index, p-values
Honmozygosity index, p-values
The HH region versus the NHH region
p < 0.0001
Z = -6.287
p < 0.0001
Z = -21.664
The stems of the HH region versus the NHH region
p < 0.0001
Z = -13.033
p < 0.0001
Z = -19.622
The loops of the HH region versus the NHH region
p = 0.019
Z = -2.347
p < 0.0001
Z = -5,542
The stems of the HH region versus the loops of the HH region
p < 0.0001
Z = -12.328
p < 0.0001
Z = -12.199
The double HH region versus the NHH region
p < 0.0001
Z = -12.328
p < 0.0001
Z =- 28.221
The stems of the double-HH region versus NHH region
p < 0.0001
Z =- 30.723
p < 0.0001
Z =- 31,211
The loops of the double-HH region versus the NHH region
p < 0.0001
Z =- 5.670
p < 0.0001
Z =- 10.690
The stems of the double-HH region versus the loops of the HH region
p < 0.0001
Z =- 29.081
p < 0.0001
Z =- 25.513
Mann-Whitney U Tests of the most frequent nucleotide in the HH region versus NHH region (Approach 2 see Methods).
Double-hammerhead (68 bp)
Single-hammerhead (49 bp)
The HH region versus the NHH region
p = 0.001 n1 = 429 n2 = 68 Z = -3.214
p = 0.142 n1 = 448 n2 = 49 Z =- 1.469
The stems of the HH region versus the NHH region
p = 0.006 n1 = 429 n2 = 35 Z = -2.762
p = 0.026 n1 = 28 n2 = 448 Z = -2.226
The loops of the HH region versus the NHH region
p = 0.801 n1 = 429 n2 = 33 Z = -0.252
p = 0.281 n1 = 21 n2 = 448 Z = -1.078
The stems of the HH region versus the loops of the HH region
p = 0.064 N1 = 35 n2 = 33 Z = -1.854
p = 0.457 n1 = 28 n2 = 21 Z = -0.743
Figure 2 illustrates the sequence variation in both the single and the double HH region. The variation affects mainly those bases that are not involved in either the hypothetical base pairing in the stems or the core region in the catalytic centre of the ribozyme. The hypothetical cleavage site is altered in only 2 clones. In the uridine turn most mutations affect position U7, which is the only position in the core region that may vary [26, 36]. The most common mutation in U7 is the C→T transition, which is species specific for the two closely related D. linderi and D. bolivari. The positions C3 and G8 form a Watson-Crick base pair in the three-dimensional structure of the HH ribozyme, and these positions are highly conserved and mutated in only four and two sequences, respectively. G12 which interacts with the cleavage site, is mutated in only one sequence. Of the two other core positions A13 and A14, only A13 is altered in two sequences.
There is a relatively high number of mutations in the stem I region, although the two positions closest to the core , i.e. C1.1 and G2.2, are almost invariant. The positions G1.2 and C2.2 are altered in a number of sequences, but frequently these changes are compensatory mutations that sustain the Watson-Crick base pairing. Compensatory mutations were also found in stem II where they affect the conserved positions G10.1 and C11.1. Stem III contains a non-canonical base pair which in most sequences is G-U, which is the most common wobble found in RNA secondary structures [37, 38]. The non-canonical base pair U-U was also found in this stem and was species specific to D. schiavazzii.
Altogether, there are 19 compensatory mutations distributed in 13 sequences that can potentially restore the base pairing stems (in the total alignment of 198 sequences). Six of them fold into an active HH structure (see below). The remaining seven sequences cannot fold into an active HH structure due to other mutations.
The number of pDo500 sequences per species that contains or lacks the proposed functional HH ribozyme sequence motif.
Total number of pDo500sequences per species
Number of sequences with the Pst3-like structure
Number of sequences with the For6-2-like structure
Total number of sequences with a potentially active hammerhead structure (A+B)
Number of sequences without a potentially active hammerhead structure
Folding analyses were also performed for population specific consensus sequences as determined previously . Again, the vast majority of sequences, i.e. 34 out of 39, folded into active HH ribozyme structures. The consensus sequences of D. schiavazzii and the PRA and TUS populations of D. geniculata formed the For6-2 structure, whereas all other consensus sequences folded into the Pst3 structure.
The HH ribozyme is the best studied ribozyme structure to date . However, HH ribozymes have still not been found in a large number of species. Of the non-viroid organisms, the natural HH ribozyme has primarily been studied in Schistosoma mansoni. In the current study, we investigated the conservation of the HH ribozyme sequence, embedded in the pDo500 satDNA repeats, in 12 species of Dolichopoda cave crickets. This is currently the most comprehensive dataset of HH ribozyme sequences in terms of the number of species involved.
Here we applied a phylogenetic footprinting approach [40, 41] in order to assess the phylogenetic conservation of the HH ribozyme region. Using three different measures of sequence variation, we found a significantly higher level of sequence conservation for the HH region than for the NHH region of the pDo500 repeats. This clearly suggests an evolutionary constraint on the HH structure in the pDo500 satellite family, which has been hypothetized earlier on the basis of self cleavage of pDo500 transcripts from D. baccettii . In particular pairing stem regions as well as the core of the HH ribozyme are highly conserved in the pDo500 satDNA family.
SatDNA is known for its high evolutionary turnover; it may thus be surprising to find highly conserved HH ribozyme sequences within the pDo500 satDNA family. However, the pDo500 satDNA was characterized earlier as relatively homogeneous. Bachmann et al.  noted that the pDo500 satDNA was relatively conserved and low in copy number in D. schiavazzi, at least compared to the two species-specific satDNA families pDoP102  and pDsPv400 . Two trends for the mode of evolution of the three satDNA families were detected: 1) a positive correlation of sequence variability and copy number, and 2) a negative correlation between sequence variability and length of repeat. While the first trend was considered in line with the theory on molecular evolution of satDNA, the second trend was not (e.g. . High sequence homogeneity is expected to correlate to high recombination rate, which in turn should reduce repeat length. The results presented here add a further aspect to the mode of evolution of the pDo500 satDNA in Dolichopoda - the evolutionary constraint brought upon the pDo500 repeats by the conserved HH ribozyme may explain the unexpected negative correlation between sequence variability and repeat length described by Bachmann et al. .
In the pDo500 data set we find compensatory mutations in at least 13 sequences. The compensatory mutations were observed in stems I and II and can restore the hypothetical base pairing in the single-HH structure. If we consider the double-HH, the number of potential compensatory mutations may even be higher because a double-HH consists of two different repeats that need to hybridize. In stem I we found mutations in position 2.2 which in some instances were compensated for by mutations in position 1.2 in the same sequence, but in other instances were not. The compensatory mutations in the stem regions do, to some extent, covary and this is considered further evidence for a functional constraint on pDo500 satellite DNA. Although some HH sequences probably have lost the ability to form an active structure, the majority of our single-HH sequences could fold into a potentially active HH structure as previously described .
Further functional analyses are certainly necessary to address the biological significance of the pDo500 HH ribozyme in Dolichopoda cave crickets. Future experiments may focus on the transcription pattern of pDo500 sequences in the different Dolichopoda species and on the self-cleavage of pDo500 repeats that carry nucleotide substitutions in the 49 bp and 68 bp regions described above. However, the data presented herein strongly suggest that the HH ribozyme in Dolichopoda cave crickets is an example of sequence conservation in non-coding DNA due to functional importance to the organism. Nevertheless, this is certainly a special case similar to that described earlier in amphibians [19, 45]. Coding for functional ribozymes, at present, cannot be considered a general function of satDNA. There are two major known functions related to satDNA: 1) the proper functioning of the centromere [48–51], and 2) processes related to non-coding RNAs. The HH ribozyme in Dolichopoda is a special case of category two. The encoding of the HH ribozyme may explain the slow evolutionary rate for the pDo500 satDNA family.
PDo500 satDNA sequences from 38 populations representing 12 Dolichopoda cave cricket species from Italy and Spain that were previously described were included in the analyses [33, 34]. The 198 sequences were downloaded from GenBank; accession numbers GU322143-GU322341 . Details on the cloning and PCR amplification strategies, as well as the sequence alignment of the pDo500 satDNA sequences, were described previously . An alignment of the HH region is provided in Additional File 2.
The region comprising the potential HH structure was compared with the detailed descriptions of recently described active HH ribozymes [22, 26, 27, 36, 52]. Numbering of specific nucleotides comprising the HH ribozyme was according to Hertel et al. .
Two mechanisms of cleavage were suggested in a previous study of the HH ribozymes in Dolichopoda : 1) a single cis-cleavage mechanism carried out by a 49 bp sequence folding into a regular HH ribozyme, and 2) a trans-cleavage mechanism carried out by a dHH structure where two copies of an extended HH sequence of 68 bp pair with each other. The latter was thought to explain the efficient processing of long primary transcripts of the pDo500 satDNA in vivo. The current study focused on both the 49 bp sequence (Figure 1) and the extended 68 bp region (Figure 2).
Genetic variation in the potential HH region was tested against the genetic variation in the remainder of the pDo500 sequences following two different approaches: 1) Two different measures of variation were estimated for each position in the alignment and subsequently used in bootstrap approaches. The Shannon entropy  was calculated as , where p is the proportion of each character state in each position and S is the number of characters. For calculating the Shannon entropy, gaps were treated as a fifth character state; i.e. S = 4 at positions without gaps and S = 5 at positions with gaps. A variability estimate was calculated in analogy to the homozygosity index  as follows H = [A]2+[T]2 [C]2+[G]2. In positions where gaps amount to more than 3% the equation was extended to H = [A]2+[T]2 [C]2+[G]+[-]2. One thousand bootstrap replicates were run for both the Shannon entropy and the homozygosity index. 2) The frequency of the nucleotide at each position of the alignment was determined for both the HH region and the non-HH (NHH) region of the pDo500 sequence and illustrated by histograms (see Additional File 1). For all three statistical procedures (Shannon entropy, homozygosity index and tests of nucleotide frequency), each of the following comparisons were made: the HH region versus the NHH region of the pDo500 repeats; the stems of the HH region versus the NHH region of the pDo500 monomer; the loops of the HH region versus the NHH region of the pDo500 monomer; the stems of the HH region versus the loops of the HH region. Each comparison was tested in a Mann-Whitney U test using the program SPSS . Both approaches were used on both the 49 bp and the 68 bp HH regions.
Secondary structure of 198 HH sequences was predicted and illustrated with the programs MFOLD - using the web based program Quickfold that allows folding of many sequences simultaneously on the DINAMelt server [57, 58] - and RNAdraw. The proportion of sequences per species that may fold into potentially active HH ribozyme structures was determined. As criteria indicating potentially active HH ribozyme structures, we used secondary structure similarity of the sequences to the structure of either the Pst3 or For6-2 sequences of the pDo500 satDNA family which were previously identified as active HH ribozyme sequences in Dolichopoda . In addition, we used invariance of the core region since it was shown that mutations in the core region abolish or significantly reduce HH cleavage . These analyses were done only on the 49 bp HH sequence.
We want to thank Eirik Rindal and Hallvard Elven for help with analyses, Tore Schweder and Øyvind Hammer for advice on statistical analyses, and Fahri Saatcioglu for advice in writing the manuscript. This work was supported by the 'National Centre for Biosystematics' (Project no. 146515/420), co-funded by the Norwegian Research Council and the Natural History Museum, University of Oslo, Norway.
- Beridze T: Satellite DNA. 1986, Berlin: Springer VerlagView ArticleGoogle Scholar
- Bosco G, Campbell P, Leiva-Neto JT, Markow TA: Analysis of Drosophila species genome size and satellite DNA content reveals significant differences among strains as well as between species. Genetics. 2007, 177 (3): 1277-1290. 10.1534/genetics.107.075069.PubMed CentralView ArticlePubMedGoogle Scholar
- Deryusheva S, Krasikova A, Kulikova T, Gaginskaya E: Tandem 41-bp repeats in chicken and Japanese quail genomes: FISH mapping and transcription analysis on lampbrush chromosomes. Chromosoma. 2007, 116 (6): 519-530. 10.1007/s00412-007-0117-5.View ArticlePubMedGoogle Scholar
- Rouleux-Bonnin F, Bigot S, Bigot Y: Structural and transcriptional features of Bombus terrestris satellite DNA and their potential involvement in the differentiation process. Genome. 2004, 47 (5): 877-888. 10.1139/g04-053.View ArticlePubMedGoogle Scholar
- Renault S, Rouleux-Bonnin F, Periquet G, Bigot Y: Satellite DNA transcription in Diadromus pulchellus (Hymenoptera). Insect Biochem Molec Biol. 1999, 29 (2): 103-111. 10.1016/S0965-1748(98)00113-1.View ArticleGoogle Scholar
- Epstein L, Mahon K, Gall J: Transcription of a satellite DNA in the newt. J Cell Biol. 1986, 103 (4): 1137-1144. 10.1083/jcb.103.4.1137.View ArticlePubMedGoogle Scholar
- Komine Y, Tanaka NK, Yano R, Takai S, Yuasa S, Shiroishi T, Tsuchiya K, Yamamori T: A novel type of non-coding RNA expressed in the rat brain. Mol Brain Res. 1999, 66 (1-2): 1-13. 10.1016/S0169-328X(98)00343-X.View ArticlePubMedGoogle Scholar
- Carchilan M, Delgado M, Ribeiro T, Costa-Nunes P, Caperta A, Morais-Cecilio L, Jones RN, Viegas W, Houben A: Transcriptionally active heterochromatin in rye B chromosomes. Plant Cell. 2007, 19 (6): 1738-1749. 10.1105/tpc.106.046946.PubMed CentralView ArticlePubMedGoogle Scholar
- Martins C, Baptista CS, Ienne S, Cerqueira GC, Bartholomeu DC, Zingales B: Genomic organization and transcription analysis of the 195-bp satellite DNA in Trypanosoma cruzi. Mol Biochem Parasitol. 2008, 160 (1): 60-64. 10.1016/j.molbiopara.2008.03.004.View ArticlePubMedGoogle Scholar
- Pezer Z, Ugarkovic D: RNA pol II promotes transcription of centromeric satellite DNA in beetles. PLoS ONE. 2008, 3 (2): e1594-10.1371/journal.pone.0001594.PubMed CentralView ArticlePubMedGoogle Scholar
- Carthew RW, Sontheimer EJ: Origins and mechanisms of miRNAs and siRNAs. Cell. 2009, 136 (4): 642-655. 10.1016/j.cell.2009.01.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein E, Allis CD: RNA meets chromatin. Genes Dev. 2005, 19 (14): 1635-1655. 10.1101/gad.1324305.View ArticlePubMedGoogle Scholar
- Martens JHA, O'Sullivan RJ, Braunschweig U, Opravil S, Radolf M, Steinlein P, Jenuwein T: The profile of repeat-associated histone lysine methylation states in the mouse epigenome. EMBO Journal. 2005, 24 (4): 800-812. 10.1038/sj.emboj.7600545.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferbeyre G, Smith JM, Cedergren R: Schistosome satellite DNA encodes active hammerhead ribozymes. Mol Cell Biol. 1998, 18 (7): 3880-3888.PubMed CentralView ArticlePubMedGoogle Scholar
- Rojas AA, Vazquez-Tello A, Ferbeyre G, Venanzetti F, Bachmann L, Paquin B, Sbordoni V, Cedergren R: Hammerhead-mediated processing of satellite pDo500 family transcripts from Dolichopoda cave crickets. Nucl Acids Res. 2000, 28 (20): 4037-4043. 10.1093/nar/28.20.4037.PubMed CentralView ArticlePubMedGoogle Scholar
- Epstein LM, Gall JG: Self-cleavage transcripts of satellite DNA from the newt. Cell. 1987, 48 (3): 535-543. 10.1016/0092-8674(87)90204-2.View ArticlePubMedGoogle Scholar
- Prody GA, Bakos JT, Buzayan JM, Schneider IR, Bruening G: Autolytic processing of dimeric plant virus satellite RNA. Science. 1986, 231 (4745): 1577-1580. 10.1126/science.231.4745.1577.View ArticlePubMedGoogle Scholar
- Hutchins CJ, Rathjen PD, Forster AC, Symons RH: Self-cleavage of plus and minus RNA transcripts of avocado sunblotch viroid. Nucl Acids Res. 1986, 14 (9): 3627-3640. 10.1093/nar/14.9.3627.PubMed CentralView ArticlePubMedGoogle Scholar
- Denti MA, Alba AEMd, Sagesser R, Tsagris M, Tabler M: A novel RNA-binding protein from Triturus carnifex identified by RNA-ligand screening with the newt hammerhead ribozyme. Nucl Acids Res. 2000, 28 (5): 1045-1052. 10.1093/nar/28.5.1045.PubMed CentralView ArticlePubMedGoogle Scholar
- Vazquez-Tello A, Castan P, Moreno R, Smith JM, Berenguer J, Cedergren R: Efficient trans-cleavage by the Schistosoma mansoni SM alpha 1 hammerhead ribozyme in the extreme thermophile Thermus thermophilus. Nucl Acids Res. 2002, 30 (7): 1606-1612. 10.1093/nar/30.7.1606.PubMed CentralView ArticlePubMedGoogle Scholar
- Canny MD, Jucker FM, Pardi A: Efficient ligation of the Schistosoma hammerhead ribozyme. Biochemistry. 2007, 46 (12): 3826-3834. 10.1021/bi062077r.PubMed CentralView ArticlePubMedGoogle Scholar
- Martick M, Horan LH, Noller HF, Scott WG: A discontinuous hammerhead ribozyme embedded in a mammalian messenger RNA. Nature. 2008, 454 (7206): 899-902. 10.1038/nature07117.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammann C, Westhof E: The unforeseeable hammerhead ribozyme. Biology Reports Ltd. 2009, 1 (6):Google Scholar
- Nelson JA, Uhlenbeck OC: When to Believe What You See. Mol Cell. 2006, 23 (4): 447-450. 10.1016/j.molcel.2006.08.001.View ArticlePubMedGoogle Scholar
- Osborne EM, Schaak JE, DeRose VJ: Characterization of a native hammerhead ribozyme derived from schistosomes. RNA. 2005, 11 (2): 187-196. 10.1261/rna.7950605.PubMed CentralView ArticlePubMedGoogle Scholar
- Martick M, Scott WG: Tertiary contacts distant from the active site prime a ribozyme for catalysis. Cell. 2006, 126 (2): 309-320. 10.1016/j.cell.2006.06.036.PubMed CentralView ArticlePubMedGoogle Scholar
- Khvorova A, Lescoute A, Westhof E, Jayasena SD: Sequence elements outside the hammerhead ribozyme catalytic core enable intracellular activity. Nat Struct Biol. 2003, 10 (9): 708-10.1038/nsb959.View ArticlePubMedGoogle Scholar
- Blount KF, Uhlenbeck OC: The structure-function dilemma of the hammerhead ribozyme. Annu Rev Biophys Biomol Struct. 2005, 34 (1): 415-440. 10.1146/annurev.biophys.34.122004.184428.View ArticlePubMedGoogle Scholar
- Nelson JA, Uhlenbeck OC: Hammerhead redux: Does the new structure fit the old biochemical data?. RNA. 2008, 14 (4): 605-615. 10.1261/rna.912608.PubMed CentralView ArticlePubMedGoogle Scholar
- Scott WG: Ribozymes. Curr Opin Struct Biol. 2007, 17 (3): 280-286. 10.1016/j.sbi.2007.05.003.View ArticlePubMedGoogle Scholar
- Murray JB, Terwey DP, Maloney L, Karpeisky A, Usman N, Beigelman L, Scott WG: The structural basis of hammerhead ribozyme self-cleavage. Cell. 1998, 92 (5): 665-673. 10.1016/S0092-8674(00)81134-4.View ArticlePubMedGoogle Scholar
- Westhof E: A tale in molecular recognition: the hammerhead ribozyme. J Mol Recogn. 2007, 20 (1): 1-3. 10.1002/jmr.806.View ArticleGoogle Scholar
- Martinsen L, Venanzetti F, Johnsen A, Bachmann L: Molecular evolution of the pDo500 family in Dolichopoda cave crickets (Rhaphidophoridae). BMC Evol Biol. 2009, 9: 301-10.1186/1471-2148-9-301.PubMed CentralView ArticlePubMedGoogle Scholar
- Bachmann L, Venanzetti F, Sbordoni V: Tandemly repeated satellite DNA of Dolichopoda schiavazzii: a test for models on the evolution of highly repetitive DNA. J Mol Evol. 1996, 43 (2): 135-144. 10.1007/BF02337358.View ArticlePubMedGoogle Scholar
- Allegrucci G, Todisco V, Sbordoni V: Molecular phylogeography of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae): A scenario suggested by mitochondrial DNA. Mol Phylogenet Evol. 2005, 37 (1): 153-164. 10.1016/j.ympev.2005.04.022.View ArticlePubMedGoogle Scholar
- Przybilski R, Hammann C: The tolerance to exchanges of the Watson Crick base pair in the hammerhead ribozyme core is determined by surrounding elements. RNA. 2007, 13 (10): 1625-1630. 10.1261/rna.631207.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagaswamy U, Voss N, Zhang Z, Fox GE: Database of non-canonical base pairs found in known RNA structures. Nucl Acids Res. 2000, 28 (1): 375-376. 10.1093/nar/28.1.375.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu D, Landon T, Greenbaum NL, Fenley MO: The electrostatic characteristics of G·U wobble base pairs. Nucl Acids Res. 2007, 35 (11): 3836-3847. 10.1093/nar/gkm274.PubMed CentralView ArticlePubMedGoogle Scholar
- Hean J, Weinberg MS: The hammerhead ribozyme revisited: New biological insights for the development of therapeutic agents and for reverse genomics applications. RNA and the Regulation of Gene Expression: A Hidden Layer of Complexity. Edited by: Morris KV. 2008, Caister Academic Press, 1-18.Google Scholar
- Allende ML, Manzanares M, Tena JJ, Feijóo CG, Gómez-Skarmeta JL: Cracking the genome's second code: Enhancer detection by combined phylogenetic footprinting and transgenic fish and frog embryos. Methods. 2006, 39 (3): 212-10.1016/j.ymeth.2005.12.005.View ArticlePubMedGoogle Scholar
- Tagle DA, Koop BF, Goodman M, Slightom JL, Hess DL, Jones RT: Embryonic e and g globin genes of a prosimian primate (Galago crassicaudatus): Nucleotide and amino acid sequences, developmental regulation and phylogenetic footprints. J Mol Biol. 1988, 203 (2): 439-455. 10.1016/0022-2836(88)90011-3.View ArticlePubMedGoogle Scholar
- Bachmann L, Venanzetti F, Sbordoni V: Characterization of a species-specific satellite DNA family of Dolichopoda schiavazzii (Orthoptera, Rhaphidophoridae) cave crickets. J Mol Evol. 1994, 39 (3): 274-281. 10.1007/BF00160151.View ArticlePubMedGoogle Scholar
- Stephan W, Cho S: Possible role of natural selection in the formation of tandem-repetitive noncoding DNA. Genetics. 1994, 136 (1): 333-341.PubMed CentralPubMedGoogle Scholar
- Flores R, Navarro JA, de la Pena M, Navarro B, Ambros S, Vera A: Viroids with hammerhead ribozymes: Some unique structural and functional aspects with respect to other members of the group. Biol Chem. 1999, 380 (7-8): 849-854. 10.1515/BC.1999.104.View ArticlePubMedGoogle Scholar
- Luzi E, Eckstein F, Barsacchi G: The newt ribozyme is part of a riboprotein complex. Proc Natl Acad Sci USA. 1997, 94 (18): 9711-9716. 10.1073/pnas.94.18.9711.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Lloyd ME: Cloning and characterization of extended hammerheads from a diverse set of caudate amphibians. Gene. 1996, 172 (2): 183-190. 10.1016/0378-1119(96)00126-6.View ArticlePubMedGoogle Scholar
- Garrett TA, Pabón-Peña LM, Gokaldas N, Epstein LM: Novel requirements in peripheral structures of the extended satellite 2 hammerhead. RNA. 1996, 2 (7): 699-706.PubMed CentralPubMedGoogle Scholar
- Koch J: Neocentromeres and alpha satellite: a proposed structural code for functional human centromere DNA. Hum Mol Genet. 2000, 9 (2): 149-154. 10.1093/hmg/9.2.149.View ArticlePubMedGoogle Scholar
- Rudd MK, Mays RW, Schwartz S, Willard HF: Human artificial chromosomes with alpha satellite-based de novo centromeres show increased frequency of nondisjunction and anaphase lag. Mol Cell Biol. 2003, 23 (21): 7689-7697. 10.1128/MCB.23.21.7689-7697.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Yunis JJ, Yasmineh WG: Heterochromatin, satellite DNA, and cell function. Science. 1971, 174 (4015): 1200-1209. 10.1126/science.174.4015.1200.View ArticlePubMedGoogle Scholar
- Vagnarelli P, Ribeiro SA, Earnshaw WC: Centromeres: Old tales and new tools. FEBS Lett. 2008, 582 (14): 1950-1959. 10.1016/j.febslet.2008.04.014.View ArticlePubMedGoogle Scholar
- Przybilski R, Hammann C: The hammerhead ribozyme structure brought in line. ChemBioChem. 2006, 7 (11): 1641-1644. 10.1002/cbic.200600312.View ArticlePubMedGoogle Scholar
- Hertel KJ, Pardi A, Uhlenbeck OC, Koizumi M, Ohtsuka E, Uesugi S, Cedergren R, Eckstein F, Gerlach WL, Hodgson R, Symons RH: Numbering system for the hammerhead. Nucl Acids Res. 1992, 20 (12): 3252-10.1093/nar/20.12.3252.PubMed CentralView ArticlePubMedGoogle Scholar
- Durbin R, Eddy S, Krogh A, Mitchison G: Biological sequence analysis: probabilistic models of proteins and nucleic acids. 1998, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University PressGoogle Scholar
- SPSS: SPSS for Windows: Version: Rel. 16.0.2. Version: Rel. 126.96.36.1998 edition. 2008, City: SPSS IncGoogle Scholar
- Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction. Nucl Acids Res. 2005, W577-W581. 10.1093/nar/gki591. 33 Web serverGoogle Scholar
- Markham NR, Zuker M: UNAFold: software for nucleic acid folding and hybriziation in Methods in Molecular Biology, chapter 1. Bioinformatics, Volume II Structure, Functions and Applications. Edited by: Keith JM. 2008, Totowa: Humana Press, 3-31.Google Scholar
- Hammann C, Lilley DMJ: Folding and activity of the hammerhead ribozyme. ChemBioChem. 2002, 3 (8): 690-700. 10.1002/1439-7633(20020802)3:8<690::AID-CBIC690>3.0.CO;2-C.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.