The complete mitochondrial genome of the enigmatic bigheaded turtle (Platysternon): description of unusual genomic features and the reconciliation of phylogenetic hypotheses based on mitochondrial and nuclear DNA
© Parham et al; licensee BioMed Central Ltd. 2006
Received: 11 November 2005
Accepted: 07 February 2006
Published: 07 February 2006
The big-headed turtle (Platysternon megacephalum) from east Asia is the sole living representative of a poorly-studied turtle lineage (Platysternidae). It has no close living relatives, and its phylogenetic position within turtles is one of the outstanding controversies in turtle systematics. Platysternon was traditionally considered to be close to snapping turtles (Chelydridae) based on some studies of its morphology and mitochondrial (mt) DNA, however, other studies of morphology and nuclear (nu) DNA do not support that hypothesis.
We sequenced the complete mt genome of Platysternon and the nearly complete mt genomes of two other relevant turtles and compared them to turtle mt genomes from the literature to form the largest molecular dataset used to date to address this issue. The resulting phylogeny robustly rejects the placement of Platysternon with Chelydridae, but instead shows that it is a member of the Testudinoidea, a diverse, nearly globally-distributed group that includes pond turtles and tortoises. We also discovered that Platysternon mtDNA has large-scale gene rearrangements and possesses two, nearly identical, control regions, features that distinguish it from all other studied turtles.
Our study robustly determines the phylogenetic placement of Platysternon and provides a well-resolved outline of major turtle lineages, while demonstrating the significantly greater resolving power of comparing large amounts of mt sequence over that of short fragments. Earlier phylogenies placing Platysternon with chelydrids required a temporal gap in the fossil record that is now unnecessary. The duplicated control regions and gene rearrangements of the Platysternon mtDNA probably resulted from the duplication of part of the genome and then the subsequent loss of redundant genes. Although it is possible that having two control regions may provide some advantage, explaining why the control regions would be maintained while some of the duplicated genes were eroded, examples of this are rare. So far, duplicated control regions have been reported for mt genomes from just 12 clades of metazoans, including Platysternon.
Multiple studies have differed in the placement of Platysternon, with the results contrasted in Figure 2. Recent studies of the phylogenetic position of Platysternon using nuDNA (RAG-1 and U17 snoRNA) strongly supported testudinoid affinities [2, 3]. One of these studies  gave a detailed review of the conflicting signals from other data sets (mtDNA and morphology). These authors acknowledged the dissenting voices on the "Platysternon as a chelydrid" scheme from morphologists, but it should be noted that all such morphological hypotheses were not proposed in an explicit cladistic framework. Meanwhile, the most recent cladistic analysis of osteological characters  could not resolve the position of Platysternon beyond placing it in the same major clade (Cryptodira) that includes most extant turtle lineages including the testudinoids and chelydrids, but also softshell turtles (Trionychia), mud turtles (Kinosternoidea), and sea turtles (Chelonioidea). The combined phylogenetic analysis of short sequences of mtDNA (fragments of cob and rrnS) placed Platysternon next to chelydrids .
In order to pursue a definitive resolution of this issue, we sequenced the complete mt genome of Platysternon and nearly complete mt genome of a chelydrid and kinosternoid and compared these data to mt genomes published for other turtle lineages. In the process, we discovered several unusual mt genomic features that further distinguish this enigmatic turtle. We describe these genomic features and review the phylogenetic position of Platysternon.
Results and discussion
Phylogenetic position of Platysternon
Our phylogeny reconciles the previous conflict between mtDNA and nuDNA  by agreeing with the nuDNA data that Platysternon is a testudinoid. The fact that our large mt alignment results in a phylogenetic hypothesis that is congruent with the nuDNA rather than the analyses based on small (< 5 kb) mt sequences highlights the utility of generating large mt sequences for higher-level systematics . Because independent genetic markers (mtDNA and nuDNA) support testudinoid affinities, and there is no strong morphological argument for chelydrid affinities , the continued recognition of Platysternon as a chelydrid is no longer tenable.
The paleontological record is consistent with the "Platysternon as a testudinoid" hypothesis. The oldest fossil referred to the stem lineage of Platysternon (the Platysternidae), are from the Paleogene of Asia (55–60 mya) [12, 13], at about the same time we find the oldest testudinoids [14, 15]. Chelydrids, on the other hand, are significantly more ancient, extending back into the middle Cretaceous (~90 mya) . Consequently, the recognition of Platysternon as a testudinoid alleviates a major temporal disparity of ~30 million years. Despite this apparent congruence, it is important to realize that the reported fossil record of platysternids is poor and in need of review and confirmation . The described material from Asia is based largely on fragmentary specimens that have not been subjected to rigorous phylogenetic analysis [12, 13, 18–21]. Meanwhile, potentially relevant fossil specimens of possible platysternids in Europe  and North America ("Emydid C" ) have been mentioned in the literature, but have not been adequately described. The possibility of early platysternids in North America is especially intriguing because our study supports a sister relationship to Chrysemys, our representative of the largely North American clade Emydidae. However, until more specimens are brought to light, the paleontological perspective on platysternid origins remains highly speculative.
Genomic features of Platysternon mtDNA
One protein coding gene (nad5) and five tRNA genes are in derived positions. These have transposed from two portions of the genome (trnH, trnS, trnL, nad5 and trnT, trnP, cr) that are ancestrally near to one another, but separated by nad6, trnE, cob (Fig. 5). In Platysternon, both of these regions are inserted between trnI and trnQ, are separated by a block of non-coding sequence. This is the first true gene rearrangement reported for a turtle. In the pancake tortoise, Malacochersus, the cr and trnF are duplicated ; however, since the second cr of the pancake tortoise is highly degraded, the two trnF are essentially adjacent (i.e., no coding regions are out of sequence). The translocation of the cr and mt genes to between trnI and trnQ is interesting because this is the same position that contains a duplicated cr and rearranged tRNA genes in another reptile clade, the advanced snakes , and because this has been noted otherwise as a rearrangement "hot spot."
The arrangement of the Platysternon genome can be modelled by the "duplication-random loss" model  whereby a duplication and transposition of part of the genome occurred, then additional rearrangements resulted from the loss of supernumerary genes. Since the transposed genes in Platysternon are ancestrally separated by only a block of three genes, it may be that the originally duplicated and transposed region included the entire portion from trnH through the cr. This observation bolsters speculation that the non-coding region now found between nad5 and trnT is the degenerating vestige of what was the duplicated nad6, trnE, and/or cob and, similarly, that the non-coding region between nad4 and nad6 is the vestige of a copy of trnH, trnS, trnL, and/or nad5. The study of recent duplication events demonstrates that when parts of the genome are duplicated, redundant sequences are rapidly lost , and cr duplications have been otherwise associated with gene rearrangements .
It is unusual that there should be two similar control regions in Platysternon, and uncertain whether this indicates a very recent duplication, maintenance by selection, or some error correction mechanism resulting in their evolving in concert. Duplicated control regions have previously been reported for just 11 clades spanning the diversity of Metazoa [25, 28–30]. Some experimental data suggest that mt genomes with two crs have a selective advantage in replication over those with one cr , but there are clearly cases where one copy of a duplicated cr is degrading [24, 28, 33].
The maintenance of duplicated sequences is not restricted to crs. A recent study reported seven instances from diverse metazoans, in which reported sequences of coding regions were duplicated . To this we can add the duplication of trnK in the reptile Sphenodon . Whether all of these duplications represent cases of stable functional redundancy in coding regions or merely result from recent duplications and have not degraded into pseudogenes remains to be tested.
Platysternon is not related to chelydrids, but is instead a member of the Testudinoidea, the group that includes pond turtles and tortoises. Testudinoids diversified rapidly in Asia and North America during the Paleogene (50–60 mya) [14, 15]. Additional taxon sampling will help establish the phylogeny for extant testudinoids, including whether Platysternon is actually more closely related to emydids or testugurians. However, the best understanding of the timing and geography of this radiation will require the additional description and analysis of important, but neglected, fossil specimens.
The features of the Platysternon mitochondrial genome expand our knowledge of variation within vertebrate mitochondrial genomes, adding a new case of duplicated control regions. Moreover, the unusual mt genome of Platysternon and the pancake tortoise (Malacochersus ) are good examples of how additional sequencing of turtle mt genomes can improve our knowledge of mitochondrial variation and evolution. At the time of this writing, just ~6% of turtle diversity (18 of ~300 species) have large (> 5 kb) mt sequences reported (16 of these are complete).
Our new sequences are derived from three museum specimens: 1) Platysternon megacephalum (MVZ 230486) from Hainan Island, China; 2) Chelydra serpentina (MVZ 137436) from North Carolina, USA; 3) Kinosternon flavescens (MVZ 164999) from Texas, USA. Genomic DNA was extracted from frozen liver using the Qiagen QIAamp tissue kit. Amplification of genomic DNA was conducted using rTth long PCR enzyme (Applied Biosystems) with a denaturation at 94°C for 15 sec, annealing at 46–50°C for 20 sec, and extension at 68°C for 60 sec for a total of 38 cycles, followed by an additional extension at 72° for 12 min.
The following primers were used (listed 5' to 3'): A) TestGenPhe.f: AAAGCGTGGCATTGAAGCTG; B) 12Sa: AAACTGGGATTAGATACCCCACT; C) 16sf.2: TACGACCTCGATGTTGSATCAGG; D) TestGenCo3.f: GCTGCTTGATAYTGACACTTYGT; E) Nad4.f5: TGACTACCAAAAGCCCACGTAGA; F) 16S.r10: TCCAACATCGAGGTCGTAAACC; G) Met.r7: GCTATGGGCCCAAAAGCTT; H) Nad4.r6: TCTACGTGGGCTTTTGGTAGTCA; I) Leu.r1: TTTTACTTGGAGTTGCACCA; J) Cb.r24: CTCAGAATGATATTTGTCCTCARGG. The following primer pairs were used for each species: K. flavescens (B-I), C. serpentina (A-G, C-H, D-J), P. megacephalum (A-G, C-H, D-J, E-F).
Amplification products were sheared randomly into fragments of approximately 1.5 kb by repeated passage through a narrow aperture using a Hydroshear device. After end-repair, the sheared DNA was gel purified, ligated into a plasmid vector, and then transformed into bacterial cells to construct a library of random fragments. Automated colony pickers introduced single clones into bacterial broth with 10% glycerol in 384-well format. We sequenced 96 or 192 clones per amplification for 192–576 clones per species (192 for K. flavescens, 384 for C. serpentina, 576 for P. megacephalum). These plasmid clones were processed robotically for rolling circle amplification [34, 35]. Sequencing reactions and reaction cleanup were done using SPRI . Sequences were determined using ABI3730xl DNA sequencers and then were assembled based on overlap to form deep contigs 5X->50X).
Phylogenetic analyses and hypothesis testing
The three new DNA sequences (Platysternon [GenBank# = DQ_256377, 19,043 bp complete mt genome], Chelydra [DQ_256378, 14,567 bp sequence from rrnS to cob position 415], Kinosternon [DQ_256379, 7,288 bp sequence from trnL(taa) to trnR]) were aligned manually with those from nine other species from GenBank (Chelonia [NC_000886], Chrysemys [NC_002073], Dogania [NC_002780], Geochelone [DQ_080041], Manouria [DQ_080040], Mauremys/"Chinemys" [NC_006082], Pelodiscus [NC_006132], Pelomedusa [NC_001947], Testudo [DQ_080049]). With the exception of Platysternon, no noteworthy, unusual genomic features were found in the new sequences. However, we did note that Kinosternon lacks the "extra" nucleotide that causes a translational frameshift in nad3 in all other turtles where known [24, 37].
For our alignment, protein-coding genes were constrained to align by codon and tRNA-coding genes were constrained to align by regions of potential secondary structure . We excluded highly-variable regions that were difficult to align including the control region, 225 positions from other non-coding regions, 79 positions of rrnS, and 317 positions of rrnL. A total of 170 positions were excluded from the alignment of tRNA genes: the D-loop region was excluded from trnH, trnS, trnL(taa); and both the D- and T-loop regions were excluded from trnF, trnV, trnI, trnW, trnK, trnR, trnT, and trnP. We excluded a total of 151 positions from the protein coding gene nad5. The final alignment contains 15,289 positions and provides 4,901 parsimony informative characters.
We used maximum parsimony (MP) , maximum likelihood (ML) , and Bayesian inference (BI)  phylogenetic methods to infer phylogenetic trees. We conducted the MP and ML analyses with PAUP* 4.0b10  and BI analyses with MrBayes 3.1.1 [43, 44]. We executed MP analyses with the branch and bound search option, which guarantees an exact solution. To assess nodal support for the MP analysis, we used the bootstrap resampling method  employing 1000 pseudoreplicates of heuristic searches using TBR branch swapping and 100 random sequence additions pseudoreplication in PAUP*. We also obtained decay indices (="branch support")  for all nodes.
To determine the most appropriate model of DNA substitution for reconstructing turtle relationships under ML, we evaluated the fit of various models of molecular evolution to our data via the Akaike Information Criterion (AIC)  with the program Modeltest 3.06 . We performed ML analyses under the optimal model (GTR + I + G) with the heuristic search algorithm using TBR branch swapping with 10 random sequence additions, simultaneously estimating parameter values (with 10 Γ rate categories) and tree topology (i.e., no initial parameter estimates or starting tree). We then successively re-estimated parameter values and searched for trees until we obtained a stable topology and ML score .
We also performed ML-based BI analyses to search for additional tree topologies. Because MrBayes can perform mixed model phylogenetic analyses using different models of evolution  we assessed the best fit model of evolution for each mtDNA gene via the AIC with the program MrModeltest 2.1 . However, to avoid over-parameterization, we combined mitochondrial loci into the same data partition if they belonged to the same functional type (either rRNA, tRNA, or protein coding DNA) and conformed to the same model of evolution. This resulted in 12 partitions with the following models: (1) rrnL, rrnS= GTR+G; (2) atp6, cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad5 = GTR+I+G; (3) cob, nad6 = GTR+G; (4) atp8 = HKY+I+G; (5) nad4L = HKY+G; (6) trnA, trnD, trnG, trnQ, trnR = GTR+G; (7) trnE, trnL(nag) = GTR+I; (8) trnF = SYM+G; (9) trnM = HKY+I+G; (10) trnC, trnK, trnN, trnS(nga), trnT, trnV, trnY = HKY+G; (11) trnH, trnP, trnS(nct) = HKY+I; (12) trnI, trnL(taa), trnW = K80+G.
We then performed mixed-model BI tree searches, allowing separate parameter estimates under the chosen models of DNA substitution for each data partition. We did not specify nucleotide substitution model parameters or a topology a priori. We ran BI analyses for 3 × 106 generations using the default temperature (0.2) with four Markov chains per generation, sampling trees every 100 generations. We then computed a 50% majority rule consensus tree after excluding those trees sampled prior to the stable equilibrium (after the first 1 × 105 generations). Nodal support is given by the frequency of the recovered clade, which corresponds to the posterior probability of that clade under the assumed models of sequence evolution [43, 51].
We assessed the congruence between our hypothesized placement of Platysternon and those proposed by other molecular genetic analyses using constraint searches and subsequent topology tests in PAUP*. First, we constrained the MP and ML searches to retain only those trees with a Platysternon + Chelydra clade, consistent with previous mtDNA analyses . Second, we constrained the MP and ML searches to retain only those trees with a Platysternon + Testuguria clade, consistent with a previous nuDNA analysis . We then compared the constrained and unconstrained MP estimates of turtle phylogeny using a two-tailed Wilcoxon signed-ranks test , and compared the constrained and unconstrained ML phylogenies using a one-tailed multiple-comparisons likelihood ratio test  with 1000 RELL bootstrap pseudoreplicates.
Finally, we also performed phylogenetic analyses of a data matrix that combined our mtDNA data with nuDNA from two relevant studies [2, 3]. The combined analyses were performed using the same parameters used for the mtDNA analyses given above and with the models for the nuDNA specified in those other studies [2, 3]. Because there is non-overlapping taxonomic coverage between the three studies (ours and the two nuclear studies) we had to use Operational Taxonomic Units (OTUs) that had data from more than one species. These "chimeras" are a major problem in turtle systematics, especially in paleontological studies where the inclusion of broadly paraphyletic OTUs is a recurring phenomenon . We tried to avoid this problem by combining nuDNA and mtDNA sequences from only the most closely related taxa to ensure that our OTUs would be monophyletic with respect to one another. The one exception is the trionychids. In that case we combined the RAG-1 sequence for Apalone with the large mt sequence from Pelodiscus (no U17 snoRNA data is available for any trionychid). This combination was arbitrary since Apalone is just as closely related to Dogania as Pelodiscus, but this should not impact our results since all studies agree on the phylogentic position of these taxa within Cryptodira. The following list gives the OTU name used in Figure 5 followed by the accession numbers for the nuDNA sequences used (EMBL number for U17 snoRNA, GenBank number RAG-1): Pelomedusa (AJ306565, AY687922), Trionychidae [Dogania (no nuDNA), Pelodiscus (no U17 snoRNA, AY687901 from Apalone), the analyses were run with two separate trionychid OTUs and they were collapsed into a single terminal in Figure 5], Kinosternidae (AJ306562, AY687911 from Sternotherus), Chelydra (AJ306559, AY687906), Chelonia (AJ493419, AY687907), Emydidae [mtDNA from Chrysemys, nuDNA from Trachemys (AJ306564, AY687915)], Platysternon (AJ493418, AY687905), Geoemydidae [mtDNA and RAG-1 from Mauremys (AY687914), U17 snoRNA from Cuora (AJ493422)], Manouria (no nuDNA), Geochelone (AJ306561, AY687912), Testudo (AJ306563, no RAG-1).
Most of the suprageneric clade names used in this study are based on a recent review of phylogenetic nomenclature for turtles . We follow all of the protocols of that study with the exception of italicizing phylogenetically-defined clade names. Although most of the relevant phylogenetic definitions can be found in that study, a few names require additional discussion. For example, the first worker to hypothesize a close affinity of Platysternon and testudinoids  also coined the name Cryptoderinea to accommodate this grouping. Cryptoderinea has been phylogenetically codified, but can only be considered valid if Platysternon is sister to Testudinoidea . If Platysternon is nested within Testudinoidea, as proposed here and in other genetic studies [2, 3], then Playsternon should be considered a testudinoid and the name Cryptoderinea should not be used .
Secondly, the phylogenetically-defined name Bataguridae was proposed for the testudinoid clade that include most Asian hard-shelled turtles . However, according to a strict application of the rules of the International Congress of Zoological Nomenclature, there is an argument for the use of the name Geoemydidae for the same clade. In order to foster consensus during the transition from Linnaean taxonomy  to PhyloCode , we use the name Geoemydidae for this group.
Finally, despite the fact that a previous study  had proposed the name "Testudinoidae" for the clade that includes geoemydids and testudinids, the phylogenetic system used here  recommended using a new name, Testuguria. Testuguria was coined because Testudinoi dae was deemed too phonetically similar to clade names of the next higher and lower levels (Testudinoi dea and Testudinidae respectively). Although not explicitly listed as an objective synonym of Testuguria, Testudinoidae was given in the list of Testudo derivatives as an example of what kind of names to avoid. It is important to note that priority can not be invoked because, at the time of this writing, there is no official starting date for the validity of phylogenetically defined definitions. When the time comes to codify these names there will have to be a discussion as to which name (Testuguria or Testudinoidae) should be used. We strongly recommend the use of Testuguria for the reasons given above.
We thank Carla Cicero and David B. Wake (Museum of Vertebrate Zoology) for assistance with the loan of the tissues used in this study. Shi Haitao (Hainan Normal University) and Ted Papenfuss (Museum of Vertebrate Zoology) helped with the acquisition of the Platysternon specimen. We thank Michael A. Thomas and Luobin Yang (Idaho State University Evolutionary Genomics Group) for computational assistance. The image for Figure 1 was provided by Peter Paul van Dijk (Conservation International). Igor Danilov (Zoological Institute of St. Petersburg, Russia) helped with the literature on fossils, while Howard Hutchison (University of California) provided additional valuable comments. This work is LBNL-59073 and was performed under the auspices of the U.S. Department of Energy, Office of Biological and Environmental Research, by the University of California, Lawrence Berkeley National Laboratory, under contract No. DE-AC02-05CH11231. This is University of California Museum of Paleontology contribution # 1909.
- Shaffer HB, Meylan P, McKnight ML: Tests of turtle phylogeny: molecular, morphological, and paleontological approaches. Syst Biol. 1997, 46: 235-268.View ArticlePubMedGoogle Scholar
- Cervelli M, Oliverio M, Bellini A, Bologna M, Cecconi F, Mariottini P: Structural and sequence evolution of U17 small nucleolar RNA (snoRNA) and its phylogenetic congruence in chelonians. J Mol Evol. 2003, 57: 73-84. 10.1007/s00239-003-2453-2.View ArticlePubMedGoogle Scholar
- Krenz JG, Naylor GJP, Shaffer HB, Janzen FJ: Molecular phylogenetics and evolution of turtles. Mol Phylog Evol. 2005, 37: 178-191. 10.1016/j.ympev.2005.04.027.View ArticleGoogle Scholar
- Gaffney ES: Phylogeny of the chelydrid turtles: a study of shared derived characters in the skull. Fieldiana Geology. 1975, 33: 157-178.Google Scholar
- Brinkman DB, Wu XC: The skull of Ordosemys, an Early Cretaceous turtle from Inner Mongolia, People's Republic of China, and the interrelationships of Eucryptodira (Chelonia, Cryptodira). Paludicola. 1999, 2: 134-147.Google Scholar
- Vaillant L: Essai sur la classification générale des chéloniens. Ann Sci Natur. 1894, 16: 331-345.Google Scholar
- Williams EE: Variation and selection of the cervical central articulations of living turtles. Bull Am Mus Nat Hist. 1950, 94: 505-562.Google Scholar
- Danilov IG: Phylogenetic relationships of platysternid turtles [abstract]. Third Asian Herpetol Meetings. 1998, 14-Google Scholar
- Frair W: Taxonomic relationships among chelydrid and kinosternids turtles elucidated by serological tests. Copeia. 1972, 97-108.Google Scholar
- Danilov IG, Parham JF: A reassessment of the referral of an isolated skull from the Late Cretaceous of Uzbekistan to the stem-testudinoid turtle genus Lindholmemys. J Vert Paleo. 2005,Google Scholar
- Cummings MP, Otto SP, Wakeley J: Sampling properties of DNA sequence data in phylogenetic analysis. Mol Biol Evol. 1995, 12: 814-822.PubMedGoogle Scholar
- Nessov LA, Chkhikvadze VM: [New evidence on Paleocene turtle remains from South Kazakhstan]. Soobsh AN Gruzin SSR. 1987, 125: 177-180.Google Scholar
- Chkhikvadze VM: [Paleogene turtles of the USSR]. Izdatel 'Metsnier'. 1990, 1-96. [in Russian]Google Scholar
- Hutchison JH: Turtles across the Paleocene/Eocene epoch boundary in west-central North America. Late Paleocene-Early Eocene Climatic and Biotic Events in the Marine and Terrestrial Records. Edited by: Aubry M-P, Lucas SG, Berggren WA. 1998, Princeton: Princeton University Press, 401-408.Google Scholar
- Holroyd PA, Parham JF: The antiquity of African tortoises. J Vert Paleo. 2003, 23: 688-690.View ArticleGoogle Scholar
- Eaton JG, Cifelli RL, Hutchison JH, Kirkland JI, Parrish JM: Cretaceous vertebrate faunas from the Kaiparowits Plateau, south-central Utah. Utah Geological Survey Misc Pub. 1999, 1: 345-353.Google Scholar
- Danilov IG: Die fossilen Schildkröten Europas. Handbuch der Reptilien und Amphibien Europas. Edited by: Fritz U. 2005, Wiebelsheim: AULA-Verlag, 329-441.Google Scholar
- Chkhikvadze VM: [The first finding of a Tertiary turtle of the family Platysternidae]. Paleont Zhurn. 1971, 4: 137-149. [in Russian]Google Scholar
- Chkhikvadze VM: [New data on some fossil turtles from Mongolia, China and East Kazakhstan]. Soobsh AN Gruzin SSR. 1976, 82: 745-748. [in Russian]Google Scholar
- Chkhikvadze VM: [On the question of the origin of bigheaded turtles]. Izdatel 'Metsnier'. 1981, 131-146. [in Russian]Google Scholar
- Chkhikvadze VM: [Neogene turtles of the USSR]. Izdatel 'Metsnier'. 1989, 1-104. [in Russian]Google Scholar
- Grossens-Van Dyck MC: Les tortues du Paléocène continental du Hainin et Vinalmont (Belgique). Studia Palaeocheloniologica, 1984. 1985, 1: 133-139.Google Scholar
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.PubMed CentralView ArticlePubMedGoogle Scholar
- Parham JF, Macey JR, Papenfuss TJ, Feldman CR, Türkozan O, Polymeni RP, Boore JL: The phylogeny of Mediterranean tortoises and their close relatives based on complete mitochondrial genome sequences from museum specimens. Mol Phylog Evol.Google Scholar
- Dong S, Kumazawa Y: Complete mitochondrial DNA sequences of six snakes: phylogenetic relationships and molecular evolution of genomic features. J Mol Evol. 2005, 61: 12-22. 10.1007/s00239-004-0190-9.View ArticlePubMedGoogle Scholar
- Boore JL: The duplication/random loss model for gene rearrangement exemplified by mitochondrial genomes of deuterostome animals. Comparative Genomics. Edited by: Sankoff D, Nadeau JH. 2000, Dordrecht: Kluwer Academic Publishers, 133-147.View ArticleGoogle Scholar
- Moritz C: Evolutionary dynamics of mitochondrial DNA duplications in parthenogenetic geckos, Heteronotia binoei. Genetics. 1991, 129: 221-230.PubMed CentralPubMedGoogle Scholar
- Shao R, Barker SC, Mitani H, Aoki Y, Fukunaga M: Evolution of duplicate control regions in the mitochondrial genomes of Metazoa: a case study with Australasian Ixodes ticks. Mol Biol Evol. 2005, 23: 620-629.Google Scholar
- Rest JS, Ast JC, Austin CC, Waddell PJ, Tibbetts EA, Hay JM, Mindell DP: Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome. Mol Phylog Evol. 2003, 29: 289-297. 10.1016/S1055-7903(03)00108-8.View ArticleGoogle Scholar
- Abbot CL, Double MC, Trueman JWH, Robinson A, Cockburn A: An unusual source of apparent mitochondrial heteroplasmy: duplicate mitochondrial control regions in Thalassarche albatrosses. Mol Ecol. 2005, 14: 3605-3613. 10.1111/j.1365-294X.2005.02672.x.View ArticleGoogle Scholar
- Mueller RL, Boore JL: Molecular mechanisms of extensive mitochondrial gene rearrangement in plethodontid salamanders. Mol Biol Evol. 2005, 22: 2104-2112. 10.1093/molbev/msi204.View ArticlePubMedGoogle Scholar
- Tang YY, Manfredi G, Hirano M, Schon EA: Maintenance of human rearranged mitochondrial DNAs in long-term cultured transmitochondrial cell lines. Mol Biol Cell. 2000, 11: 2349-2358.PubMed CentralView ArticlePubMedGoogle Scholar
- Mindell DP, Sorenson MD, Dimcheff DE: Multiple independent origins of mitochondrial gene order in birds. Proc Natl Acad Sci. 1998, 95: 10693-10697. 10.1073/pnas.95.18.10693.PubMed CentralView ArticlePubMedGoogle Scholar
- Dean FB, Nelson JR, Giesler TL, Lasken RS: Rapid amplification of plasmid and phage DNA using Phi 29 DNA polymerase and multiply-primed rolling circle. Genome Res. 2001, 11: 1095-1099. 10.1101/gr.180501.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawkins TL, Detter JC, Richardson PM: Whole genome amplification-applications and advances. Curr Opin Biotechnol. 2002, 13: 65-67. 10.1016/S0958-1669(02)00286-0.View ArticlePubMedGoogle Scholar
- Elkin C, Kapur H, Smith T, Humphries D, Pollard M, Hammon N, Hawkins T: Magnetic bead purification of labeled DNA fragments for high-throughput capillary electrophoresis sequencing. Biotech. 2002, 32: 1296-1302.Google Scholar
- Mindell DP, Sorensen MD, Dimcheff DE: An extra nucleotide is not translated in mitochondrial ND3 of some birds and turtles. Mol Biol Evol. 1998, 15: 1568-1571.View ArticlePubMedGoogle Scholar
- Kumazawa Y, Nishida M: Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics. J Mol Evol. 1998, 37: 380-398.Google Scholar
- Farris JS: The logical basis of phylogenetic analysis. Advances in Cladistics. Edited by: Platnick N, Funk VA. 1983, New York: Columbia University Press, 7-36.Google Scholar
- Felsenstein J: Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981, 17: 368-376. 10.1007/BF01734359.View ArticlePubMedGoogle Scholar
- Larget B, Simon DL: Markov chain monte carlo algorithms for the bayesian analysis of phylogenetic trees. Mol Biol Evol. 1999, 16: 750-759.View ArticleGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 2002, Sunderland: Sinauer AssociatesGoogle Scholar
- Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogenetic trees. Bioinf. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinf. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticleGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evol. 1985, 39: 783-791.View ArticleGoogle Scholar
- Bremer K: Branch support and tree stability. Cladistics. 1994, 10: 295-304. 10.1111/j.1096-0031.1994.tb00179.x.View ArticleGoogle Scholar
- Akaike H: A new look at the statistical model identification. IEEE Trans Auto Cont. 1974, 19: 716-723. 10.1109/TAC.1974.1100705.View ArticleGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinf. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticleGoogle Scholar
- Wilgenbusch J, de Queiroz K: Phylogenetic relationships among the phrynosomatid sand lizards inferred from mitochondrial DNA sequences. Syst Biol. 2000, 49: 592-612. 10.1080/10635159950127411.View ArticlePubMedGoogle Scholar
- Nylander JAA: MrModeltest 2.1. 2004, Uppsala: Evolutionary Biology CentreGoogle Scholar
- Rannala B, Yang ZH: Probability distribution of molecular evolutionary trees: a new method of phylogenetic inference. J Mol Evol. 1996, 43: 304-311. 10.1007/PL00006090.View ArticlePubMedGoogle Scholar
- Templeton AR: Phylogenetic inference from restriction site endonuclease cleavage site maps with particular reference to the humans and apes. Evol. 1983, 37: 221-244.View ArticleGoogle Scholar
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116.View ArticleGoogle Scholar
- Joyce WG, Parham JF, Gauthier JA: Developing a protocol for the conversion of rank-based taxon names to phylogenetically defined clade names, as exemplified by turtles. J Paleontol. 2004, 78: 989-1013.View ArticleGoogle Scholar
- International Trust for Zoological Nomenclature: International Code of Zoological Nomenclature. 1999, London, 4Google Scholar
- PhyloCode. A phylogenetic code of biological nomenclature, version 2a. [http://www.ohio.edu/phylocode/preface.html]
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.