An evaluation of LSU rDNA D1-D2 sequences for their use in species identification
© Sonnenberg et al; licensee BioMed Central Ltd. 2007
Received: 16 November 2006
Accepted: 16 February 2007
Published: 16 February 2007
Identification of species via DNA sequences is the basis for DNA taxonomy and DNA barcoding. Currently there is a strong focus on using a mitochondrial marker for this purpose, in particular a fragment from the cytochrome oxidase I gene (COI). While there is ample evidence that this marker is indeed suitable across a broad taxonomic range to delineate species, it has also become clear that a complementation by a nuclear marker system could be advantageous. Ribosomal RNA genes could be suitable for this purpose, because of their global occurrence and the possibility to design universal primers. However, it has so far been assumed that these genes are too highly conserved to allow resolution at, or even beyond the species level. On the other hand, it is known that ribosomal gene regions harbour also highly divergent parts. We explore here the information content of two adjacent divergence regions of the large subunit ribosomal gene, the D1-D2 region.
Universal primers were designed to amplify the D1-D2 region from all metazoa. We show that amplification products in the size between 800–1300 bp can be obtained across a broad range of animal taxa, provided some optimizations of the PCR procedure are implemented. Although the ribosomal genes occur in multiple copies in the genomes, we find generally very little intra-individual polymorphism (<< 0.1% on average) indicating that concerted evolution is very effective in most cases. Studies in two fish taxa (genus Cottus and genus Aphyosemion) show that the D1-D2 LSU sequence can resolve even very closely related species with the same fidelity as COI sequences. In one case we can even show that a mitochondrial transfer must have occurred, since the nuclear sequence confirms the taxonomic assignment, while the mitochondrial sequence would have led to the wrong classification. We have further explored whether hybrids between species can be detected with the nuclear sequence and we show for a test case of natural hybrids among cyprinid fish species (Alburnus alburnus and Rutilus rutilus) that this is indeed possible.
The D1-D2 LSU region is a suitable marker region for applications in DNA based species identification and should be considered to be routinely used as a marker complementing broad scale studies based on mitochondrial markers.
The use of DNA sequences as a tool for the identification of species has been widely discussed in the past years [1–6]. One important step in implementing a DNA taxonomy or DNA barcoding system is the choice of appropriate markers. There is currently a preference for the COI gene, a mitochondrial protein, as the marker of choice for animals . But other markers have also been suggested, for example the nuclear rDNA [4, 7, 8] or the mitochondrial LSU gene (16S rDNA) . It is generally expected that mitochondrial markers could provide a better resolution, because their maternal inheritance results in a smaller effective population size and hence a faster fixation of neutral mutations. However, mitochondrial markers might not always reflect the full history of a species . Nuclear markers, on the other hand, underlie recombination and are therefore less suitable to trace phylogenetic lineages within species.
Apart from the question of choosing the appropriate marker type, there are also general problems with relying on a single marker system only, most notably the possibility of hybridization among lineages. Contrary to a common perception of species being entirely independent from one another, the units that are commonly considered 'species' hybridise frequently in nature [11–14]. Several model systems for interspecific hybridisation have been studied using multiple genetic markers, which suggest that mitochondrial haplotypes and possibly also nuclear genes can be exchanged between distinct species in the course of hybridisation events [15–17]. Even in the absence of hybridisation, incomplete lineage sorting can result in shared ancestral haplotype lineages beyond speciation events . The relevance of these problems in the context of taxonomic applications has been well recognized [2, 19], but has so far been very little addressed with respect to developing appropriate comparative markers.
Here we evaluate the utility of the large subunit ribosomal rRNA (LSU) as a potential marker for species identification. The LSU is part of the rDNA gene complex which occurs in tandem repeats, arranged in ribosomal clusters in the nuclear genome . Ribosomal genes are generally considered to be highly conserved, but they are actually composed of a mixture of conserved and divergent regions. These have been called "divergence regions – D" and are numbered in 5'to 3'direction of the mature rRNA . In the present paper we report primers and protocols that permit the amplification of a highly variable part of the nuclear rDNA. We find that the variation in the LSU D1-D2 fragment permits its use in DNA barcoding approaches for a wide variety of metazoa. In addition to this general utility, we explore the use of ribosomal genes to directly identify hybrids, which would not be possible with mitochondrial markers alone.
List of taxa analysed with observed sequence length range for the D1-D2 fragment in the final alignment [see additional file for further details]
number of species
D1-D2 fragment length range (bp)
We found that the invertebrate sequences could generally be amplified with standard protocols. However, vertebrate samples were generally more difficult to amplify, which cannot only be ascribed to major size differences. Interestingly, DNA traces originating from vertebrate commensals or parasites in DNA extracts were sometimes amplified, while the host sequences were not. We found that residual ribosomal rRNA was a major inhibitor in these cases, since extensive RNAseA digestion of the DNA samples leads to a much better amplification success. Since the same treatment leads also to improved results with invertebrate samples, we have included this in the recommended standard procedure.
Apart of the removal of rRNA we found that the addition of PCR enhancer buffers (see Methods) leads to significantly improved amplification of many rDNA templates, which is in agreement with the observations of Ralser et al. . This is most likely related to high GC content and pronounced secondary structure of our target gene. Finally, single strand binding (SSB)-protein was found to further improve the quality of the PCR product but is not a necessary additive in most reactions.
Further optimizations will still be necessary for some taxa. Most notably, we found that samples from mammals, birds and amphibians are difficult to amplify, i.e. the success rate is low.
Information about the affinities of a given sequence was obtained by conducting BLAST searches (blastn)  against the GenBank database  and by using simple neighbour joining clustering based on uncorrected p-distances with MEGA 3.1 . We found that the phylogenetic signal in the data is sufficient to group most sequences into taxa that are known to be monophyletic, although deep phylogenetic relationships of this tree appeared not to be very reliable (not shown). However, given our focus on a highly variable part of the LSU, we would not have expected to retrieve reliable deep level phylogeny resolution from this anyway. Still, the information content of the D1-D2 region is sufficient to place new sequences into the right phylogenetic context.
Evaluation of hybrids
Our data provide evidence that the D1-D2 region can be readily amplified from a wide range of taxa and that it contains sufficient phylogenetic information to allow assignment to related taxa. At the same time the D1-D2 region appears to carry sufficient variability to distinguish congeneric species although we have so far shown this only for two fish groups. However, at least one of these two groups may be considered as a taxonomic test case because of their close relationships. The Cottus gobio complex has previously been considered to represent a single pan-European species. Only a recent taxonomic revision has challenged this view [28, 41]. The D1-D2 LSU sequence clusters correspond to the newly described species confirming that a high level resolution is indeed possible with this gene region.
It is noteworthy that much of the variation in the D1-D2 region is found in length variable parts. While such indels cause a problem with alignments for deep phylogenetic inferences, they should be considered as an advantage for taxonomic purposes, since they add another character state beyond transitions and transversions. Moreover, indels may be particularly suitable for fast determination assays independent of sequencing, such as microarrays .
Potential problems associated with the use of ribosomal genes
Ribosomal genes are organized in clusters that contain hundreds of copies per haploid genome. It is generally assumed that these evolve in concert [23, 24]. Evidently, intra-genomic variation among these copies could cause problems. However, we have generally observed a low level of ambiguities in our sequences.
Another potential problem in studies of rDNA genes is the occurrence of pseudogenes. Recent studies show that this is a particular problem for the ITS region in animals and plants where it may confound phylogenetic inferences [31, 32]. Similarly, multiple variants for the SSU gene were found in sturgeons [33, 34], a Platyhelminth  and a dinoflagellate . However, both the study by Keller et al.  and by Fuerst and Krieger  suggest that only the functional variants are expressed at high levels and that these are not polymorphic. Hence, in cases where evidence for pseudogenes is found, it may be useful to directly sequence the ribosomal RNA to assess which of the sequence variants is the functional one.
We note that both problems listed here apply in a similar way also to mitochondrial markers. Insufficient homogenization of rDNA repeats can be compared to heteroplasmy of different mitochondrial haplotypes in a single organism (reviewed in ) and rDNA pseudogenes can be compared to nuclear copies of mitochondrial genes (e.g. [38, 39]). Thus, for either marker system one has to be aware that complications can arise through these phenomena, although they are not so frequent that they would seriously compromise their broad applicability.
rDNA variation and taxonomic group delimitation
Any new sequence variant in rDNA has to be homogenized across at least most of the other rDNA repeats before it becomes detectable. This homogenization is most likely achieved by a succession of many unequal cross-over events within the rDNA cluster. Hence, new alleles should only become visible after some time of genetic exchange has passed in an interbreeding population. However, it is currently not clear how long it takes to fix a new variant, or even whether this occurs predominantly via intra- or interchromosomal homogenization . But in any case, a new fixed variant can be seen as a reflection of an extended history of interbreeding in a given population. Accordingly, even a single mutational difference in rDNA may be taken as an informative character that could potentially delineate a new group.
For mitochondrial DNA analysis it was suggested that a threshold level of intraspecies versus interspecies variance is used as a proxy to identify groups or species . However, species recognition (as well as the underlying concepts) is not only an issue of distances but also of distinctness (e.g. ). Taxonomic resolution in the context of alpha taxonomy is either present or not and is basically given once a single diagnostic character can be identified. In the case of Cottus, we found that the different rDNA variants depicted in Fig. 4 correspond to animals from different streams and that one usually finds only one variant per stream, at least in the cases where the respective data are available. Thus, a separate group recognition based on molecular differences as small as a single substitution is corroborated by their allopatric occurrence and their presumably independent evolution since several thousand years. Accordingly, single nucleotide differences in rDNA genes may be considered as a first indicator of genetic separation. In cases where the respective groups occur under sympatric conditions, it may even be an indicator of incipient speciation. Clearly, this issue has to be further explored, but these considerations show that the collection of sequences from the coding parts of the ribosomal genes may provide an additional value for taxon delimitation.
Our results show that the D1-D2 region of the LSU rDNA gene has the potential to be developed as a taxonomic marker. It can be amplified with truly universal primers and shows a divergence rate, which is suitable to differentiate even closely related species. Because it is bi-parentally inherited, it can also be used for detecting hybrids and their corresponding parental species. It appears therefore that it could be an ideal marker for complementing DNA barcoding studies based on mitochondrial COI sequences.
Materials and methods
Macrozoobenthos invertebrates were collected and determined by M. Hess (Munich), nematodes were obtained as cultures from E. Schierenberg (Cologne), most fish material and some invertebrates were collected and/or determined by AN and RS, the Galaxias and Brachygalaxias samples by K. Busse (Bonn). F1 Hybrids of Cottus were produced in the aquarium as described in Stemshorn et al. .
For most samples total DNA was extracted from ethanol preserved or fresh material using a standard Proteinase-K in SDS/EDTA buffer [tissue digestion in 500 μL HOM buffer (0.5% SDS, 100 mM Tris-HCl, 80 mM EDTA pH 8.0) and 5 μL Proteinase-K (20 mg/mL) for at least 3 h at 55°C; addition of 500 μL NaCl (4.5 M) and 300 μL Chloroform, gentle mixing for 15 min.; centrifugation for 10 min. at 10.000 rpm, transfer of upper phase (750 μL) without interphase in new tube; precipitation with 750 μL 99% Ethanol, gentle mixing and incubation at room temperature for 5 min., centrifugation for 10 min. at 13.000 rpm, removal of supernatant; 2× washing of the pellet with 500 μL 70% Ethanol, incubation at room temperature for 5 min., centrifugation for 10 min. at 13.000 rpm, complete removal of supernatant; airdried pellet dissolved in 100 μL TE buffer (10 mM Tris-HCl, 0.1 mM EDTA pH8.0)]. Alternatively, we used a standard CTAB buffer protocol [tissue digestion in 500 μL 2% CTAB buffer (2 g/100 mL CTAB, 1.4 M NaCl, 100 mM Tris-HCl, 20 mM EDTA, pH 8.0) and 15 μL Proteinase-K (20 mg/mL) for 1 h – overnight at 64°C; 2× extraction with Chloroform/Isoamylalcohol (24:1), gentle mixing for 10 min., centrifugation for 10 min. at 13.000 rpm, transfer of upper phase without interphase in new tube; precipitation with 400 μL 98% Ethanol, incubation at room temperature for 1 h, centrifugation for 20 min. at 13.000 rpm, removal of supernatant; 2× washing of the pellet with 500 μL 75% Ethanol, centrifugation for 10 min. at 13.000 rpm, complete removal of supernatant; airdried pellet solved in 50–100 μL TE]. For some samples we used also DNA isolated with a commercial kit (procedure according to the manufacturer; Qiagen, Düsseldorf) or released the DNA with a Chelex/Proteinase-K protocol (500 μL 5% Chelex suspension, 10 μl Proteinase K (20 mg/ml]) incubation 1 h to overnight at 64°C, 15 min. 95°C denaturation of Proteinase-K (important for the following RNAseA treatment). RNA was digested for all samples with RNAseA (10 mg/ml, Fermentas) before PCR reactions. We added to 50 μL DNA solution 2 μl Fermentas RNAseA and incubated for 1–3 h at room temperature.
Primers were designed according to partial or complete LSU rDNA sequences from GenBank for a variety of taxa, ranging from plathelminths, nematodes and arthropods to vertebrates. Primer sequences are listed in Figure 2.
PCR conditions were tested with temperature and MgCl2 gradients. Amplification of LSU fragments were considerably enhanced by the addition of Q-solution (Qiagen, Düsseldorf) and single strand binding Protein (SSB, Sigma Aldrich). Final concentration of Q-solution is 1× (from 5× stock) and 1 μg SSB protein in a 20 μl PCR reaction mix.
The following PCR program was used to amplify the D1-D2 fragments: 4 min. at 94°C for initial denaturation; 45 cycles with 20 sec. 94°C, 20 sec. 52,5°C and 90 sec. 72°C, followed by 8 min. at 72°C for final extension. For most fragments we used the primer combination fw1 and rev1, some invertebrates amplified better with the combination fw1 and rev2.
In addition we amplified and sequenced a COI fragment currently applied in DNA barcoding applications for a sample of Cottus and the species of the Aphyosemion calliurum group (Cyprinodontiformes: Nothobranchiidae) for comparison with the LSU fragment. PCR primer for the amplification of the COI fragment for these taxa were taken from the literature (HCO-2198 [1, 43]) or designed according to published complete mitochondrial sequences and own data for the Cottus samples and the A. calliurum group. The following primers were used : Cottus COI forward: 5'-TTC TCG ACT AAT CAC AAA GAC ATT-3, Cottus COI reverse: 5'-TAG ACT TCA GGG TGA CCA AAG AAT CA-3, Aphyosemion forward: 5'-TAA GAA AAG GAT TTA AAC CT-3': "universal" reverse: 5'-TAA ACT TCA GGG TGA CCA AAA AAT CA-3'.
All Aphyosemion and Cottus COI PCR reactions are done in 15 microliter reactions with the Qiagen Multiplex PCR Kit, including 3 microliter Q-Solution and 0.5 microliter of a 10 pmol solution of each primer. The following PCR program was used to amplify the COI fragments: 15 min. 95°C for initial denaturation and activation of the polymerase enzyme; 45 cycles with 20 sec. 94°C, 90 sec. 52°C and 90 sec. 72°C, followed by 8 min. at 72°C for final extension.
PCR products were checked on a 1.5 % agarose gel with ethidium bromide staining (130 V, 30–40 min.), cleaned with Millipore PCR cleaning plates and sequenced according to the manual with ABI BigDye Terminator ver.3.1 in both directions on an ABI 3700. Sequencing was done with the same primers as used in the PCR reaction for both gene fragments. Very long LSU sequences (> 1100 bp), especially if they contain GC rich stretches, were sequenced in addition with the internal primers fw2 and rev4. For sequencing it sometimes turned out to be helpful to increase the amount of template DNA to get better reads in difficult sequence regions.
Contigs were assembled with Lasergene SeqMan II (DNA-Star) and resulting sequences checked against GenBank for contamination. All contigs were checked by eye for ambiguous nucleotides in the regions sequenced for both strands. We counted positions with double peaks from one third up to same height in both strands to estimate the occurrence of different alleles or copies in the rDNA cluster.
The COI and LSU sequences for Cottus and Aphyosemion were aligned with Clustal X  and checked by eye with BioEdit 5.0.9 . Aligned protein sequences were checked for a functional coding sequence to test against non-functional nuclear copies. A cluster analysis was done with the neighbour joining algorithm (NJ) as implemented in MEGA 3.1 . We employed no model of sequence evolution and used p-distances to compare only the raw data without any assumptions on sequence evolution. Missing data or gaps were deleted in the pairwise comparison.
All sequences are deposited in Genbank under the accession numbers EF416965 – EF417284).
We thank Monika Hess (ÖkoKart, Munich) and Einhard Schierenberg (Inst. for Zoology, University of Cologne) and Klaus Busse (Zoologisches Forschungsmuseum Alexander Koenig, Bonn) for their help in providing reference specimens for DNA analyses. We thank Henri Persat, Maurice Kottelat, Eric Taylor, David Neely, Joerg Freyhof, Filip Volckaert, Bernd Haenfling and Vendula Slechtova for providing DNA samples of various Cottus species. Andreas Hartl (Dorfen) has been so kind to provide the Alburnus × Rutilus hybrid as well as photographs.
- Hebert PDN, Cywinska A, Ball SL, deWaard JR: Biological identifications through DNA barcodes. Proceedings of the Royal Society of London, Series B. 2003, 270: 313-321. 10.1098/rspb.2002.2218.View ArticleGoogle Scholar
- Tautz D, Arctander P, Minelli A, Thomas RH, Vogler AP: A plea for DNA taxonomy. Trends in Ecology & Evolution. 2003, 18: 70-74. 10.1016/S0169-5347(02)00041-1.View ArticleGoogle Scholar
- Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM: Identification of Birds through DNA Barcodes. PLoS Biology. 2004, 2: 7-10.1371/journal.pbio.0020312.View ArticleGoogle Scholar
- Blaxter M: The promise of a DNA taxonomy. Philosophical Transactions of the Royal Society, London, Series B, Biological Sciences. 2004, 359: 669-679. 10.1098/rstb.2003.1447.View ArticleGoogle Scholar
- Savolainen V, Cowan RS, Vogler AP, Roderick GK, Lane R: Towards writing the encyclopedia of life: an introduction to DNA barcoding. Philosophical Transactions of the Royal Society, London, Series B, Biological Science. 2005, 360: 1805-1811. 10.1098/rstb.2005.1730.View ArticleGoogle Scholar
- Vogler AP, Monaghan MT: Recent advances in DNA taxonomy. Journal of Zoological Systematics and Evolutionary Research.
- Markmann M, Tautz D: Reverse taxonomy: an approach towards determining the diversity of meiobenthic organisms based on ribosomal RNA signature sequences. Philosophical Transactions of the Royal Society of London, Series B, Biological Sciences. 2005, 360: 1917-1924. 10.1098/rstb.2005.1723.View ArticleGoogle Scholar
- Monaghan MT, Balke M, Gregory TR, Vogler AP: DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers. Philosophical Transactions of the Royal Society of London, Series B, Biological Sciences. 2005, 360: 1925-1933. 10.1098/rstb.2005.1724.View ArticleGoogle Scholar
- Vences M, Thomas M, van der Meijden A, Chiari Y, Vieites DR: Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Front Zool. 2005, 2 (1): 5-10.1186/1742-9994-2-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballard JWO, Whitlock MC: The incomplete natural history of mitochondria. Molecular Ecology. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.View ArticlePubMedGoogle Scholar
- Grant PR, Grant BR: Hybridization of bird species. Science. 1992, 256: 193-197. 10.1126/science.256.5054.193.View ArticlePubMedGoogle Scholar
- Scribner KT, Page KS, Bartron ML: Hybridization in freshwater fishes: a review of case studies and cytonuclear methods of biological inference. Reviews in Fish Biology and Fisheries. 2001, 10: 293-323. 10.1023/A:1016642723238.View ArticleGoogle Scholar
- Seehausen O: Hybridization and adaptive radiation. Trends in Ecology and Evolution. 2004, 19: 198-207. 10.1016/j.tree.2004.01.003.View ArticlePubMedGoogle Scholar
- Gross BL, Rieseberg LH: The ecological genetics of homoploid hybrid speciation. J Hered. 2005, 96: 241-252. 10.1093/jhered/esi026.PubMed CentralView ArticlePubMedGoogle Scholar
- Nolte AW, Freyhof J, Stemshorn K, Tautz D: An invasive lineage of sculpins, Cottus sp. (Pisces, Teleostei) in the Rhine with new habitat adaptations has originated by hybridization between old phylogeographic groups. Proc Biol Sci. 2005, 272 (1579): 2379-87.PubMed CentralView ArticlePubMedGoogle Scholar
- Herder F, Nolte AW, Pfaender J, Schwarzer J, Hadiaty RK, Schliewen UK: Adaptive radiation and hybridization in Wallace's Dreamponds: evidence from sailfin silversides in the Malili Lakes of Sulawesi. Proc Biol Sci. 2006, 273 (1598): 2209-17. 10.1098/rspb.2006.3558.PubMed CentralView ArticlePubMedGoogle Scholar
- Bachtrog D, Thornton K, Clark A, Andolfatto P: Extensive introgression of mitochondrial DNA relative to nuclear genes in the Drosophila yakuba species group. Evolution. 2006, 60: 292-302. 10.1554/05-337.1.View ArticlePubMedGoogle Scholar
- Avise JC, Wollenberger K: Phylogenetics and the origin of species. Proceedings of the National Academy of Scienes USA. 1997, 94: 7748-7755. 10.1073/pnas.94.15.7748.View ArticleGoogle Scholar
- Moritz C, Cicero C: DNA Barcoding: Promise and Pitfalls. PLoS Biology. 2004, 2 (10): e354-10.1371/journal.pbio.0020354.PubMed CentralView ArticlePubMedGoogle Scholar
- Long EO, Dawid IB: Repeated genes in eukaryotes. Annual Reviews in Biochemistry. 1980, 49: 727-764. 10.1146/annurev.bi.49.070180.003455.View ArticleGoogle Scholar
- Hassouna N, Michot B, Bachellerie JP: The complete nucleotide sequence of mouse 28S rRNA gene. Implications for the process of size increase of the large subunit rRNA in higher eukaryotes. Nucleic Acids Research. 1984, 12: 3563-3583. 10.1093/nar/12.8.3563.PubMed CentralView ArticlePubMedGoogle Scholar
- Ralser M, Querfurth R, Warnatz H-J, Lehrach H, Yaspo M-L, Krobitsch S: An efficient and economic enhancer mix for PCR. Biochemical and Biophysical Research Communications. 2006, 347: 747-751. 10.1016/j.bbrc.2006.06.151.View ArticlePubMedGoogle Scholar
- Elder JF, Turner BJ: Concerted evolution of repetitive DNA sequences in eukaryotes. Quarterly Review in Biology. 1995, 70: 297-320. 10.1086/419073.View ArticleGoogle Scholar
- Liao D: Concerted evolution: molecular mechanism and biological implications. American Journal of Human Genetics. 1999, 64: 24-30. 10.1086/302221.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. Journal of Molecular Biology. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Research. 2006, 34: D16-D20. 10.1093/nar/gkj157.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
- Freyhof J, Kottelat M, Nolte AW: Taxonomic diversity of European Cottus with description of eight new species (Teleostei: Cottidae). Ichthyological Exploration of Freshwaters. 2005, 16: 107-172.Google Scholar
- Huber JH: Killi – Data 2000. 2000, Paris: Société Française d'Ichthyologie;, germanGoogle Scholar
- Busse K, Freyhof J, Zarske A: Serum-electrophoretic patterns of various cyprinid hybrids from the Rivers Elbe and Sieg (Teleostei: Cyprinidae). Faunistische Abhandlungen, Staatliches Museum für Tierkunde Dresden. 1999, 21: 265-272.Google Scholar
- Keller I, Chintauan-Marquier IC, Veltsos P, Nichols RA: Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution. Genetics. 2006, 174: 863-874. 10.1534/genetics.106.061341.PubMed CentralView ArticlePubMedGoogle Scholar
- Harpke D, Peterson A: Non-concerted ITS evolution in Mammillaria (Cactaceae). Molecular Phylogenetics and Evolution. 2006,Google Scholar
- Krieger J, Fuerst PA: Characterization of nuclear 18S rRNA gene sequence diversity and expression in an individual lake sturgeon (Acipenser fulvescens). Journal of Applied Ichthyology. 2004, 20: 433-439. 10.1111/j.1439-0426.2004.00610.x.View ArticleGoogle Scholar
- Krieger J, Hett AK, Fuerst PA, Birstein VJ, Ludwig A: Unusual intraindividual variation of the nuclear 18S rRNA gene is widespread within the acipenseridae. Journal of Heredity. 2006, 97: 218-225. 10.1093/jhered/esj035.View ArticlePubMedGoogle Scholar
- Carranza S, Giribet G, Ribera C, Barguna Riutort M: Evidence that two types of 18S rDNA coexist in the genome of Dugesia (Schmidtea) mediterranea (Platyhelminthes, Turbellaria, Tricladida). Molecular Biology and Evolution. 1996, 13: 824-832.View ArticlePubMedGoogle Scholar
- Santos SR, Kinzie RA, Sakai K, Coffroth MA: Molecular characterization of nuclear small subunit (18S)-rDNA pseudogenes in a symbiotic dinoflagellate (Symbiodinium, Dinophyta). Journal of Eukaryotic Microbiology. 2003, 50: 417-421. 10.1111/j.1550-7408.2003.tb00264.x.View ArticlePubMedGoogle Scholar
- Kmiec B, Woloszynska M, Janska H: Heteroplasmy as a common state of mitochondrial genetic information in plants and animals. Current Genetics. 2006, 50: 149-159. 10.1007/s00294-006-0082-1.View ArticlePubMedGoogle Scholar
- Pons J, Vogler AP: Complex pattern of coalescence and fast evolution of a mitochondrial rRNA pseudogene in a recent radiation of tiger beetles. Molecular Biology and Evolution. 2005, 22: 991-1000. 10.1093/molbev/msi085.View ArticlePubMedGoogle Scholar
- Hazkani-Covo E, Graur D: A Comparative Analysis of Numt Evolution in Human and Chimpanzee. Molecular Biology and Evolution. 2007, 24: 13-18. 10.1093/molbev/msl149.View ArticlePubMedGoogle Scholar
- Schlötterer C, Tautz D: Chromosomal homogeneity of Drosophila ribosomal DNA arrays suggests intrachromosomal exchanges drive concerted evolution. Current Biology. 1994, 4: 777-783. 10.1016/S0960-9822(00)00175-5.View ArticlePubMedGoogle Scholar
- Kottelat M: European freshwater fishes. An heuristic checklist of the freshwater fishes of Europe (exclusive of former USSR), with an introduction for non-systematists and comments on nomenclature and conservation. Biologia, Bratislava, Section Zoology. 1997, 52 (Supplement 5): 1-271.Google Scholar
- Stemshorn KC, Nolte AW, Tautz D: A Genetic Map of Cottus gobio (Pisces, Teleostei) based on microsatellites can be linked to the Physical Map of Tetraodon nigroviridis . Journal Evolutionary Biology. 2005, 18: 1619-1624.View ArticleGoogle Scholar
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology. 1994, 3: 294-299.PubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research. 1997, 24: 4876-4882. 10.1093/nar/25.24.4876.View ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 41: 95-98.Google 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.