Sister Group Relationship of Turtles to the Bird-Crocodilian Clade Revealed by Nuclear DNA–Coded Proteins (original) (raw)

Journal Article

Search for other works by this author on:

Accepted:

21 December 2004

Published:

29 December 2004

Cite

Naoyuki Iwabe, Yuichiro Hara, Yoshinori Kumazawa, Kaori Shibamoto, Yumi Saito, Takashi Miyata, Kazutaka Katoh, Sister Group Relationship of Turtles to the Bird-Crocodilian Clade Revealed by Nuclear DNA–Coded Proteins, Molecular Biology and Evolution, Volume 22, Issue 4, April 2005, Pages 810–813, https://doi.org/10.1093/molbev/msi075
Close

Navbar Search Filter Mobile Enter search term Search

Abstract

The phylogenetic position of turtles is a currently controversial issue. Recent molecular studies rejected a traditional view that turtles are basal living reptiles (Hedges, S. B., and L. L. Poling. 1999. A molecular phylogeny. Science83:998–1001; Kumazawa, Y., and M. Nishida. 1999. Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink, statistical evidence for archosaurian affinity of turtles. Mol. Biol. Evol.16:784–792). Instead, these studies grouped turtles with birds and crocodiles. The relationship among turtles, birds, and crocodiles remained unclear to date. To resolve this issue, we have cloned and sequenced two nuclear genes encoding the catalytic subunit of DNA polymerase α and glycinamide ribonucleotide synthetase–aminoimidazole ribonucleotide synthetase–glycinamide ribonucleotide formyltransferase from amniotes and an amphibian. The amino acid sequences of these proteins were subjected to a phylogenetic analysis based on the maximum likelihood method. The resulting tree showed that turtles are the sister group to a monophyletic cluster of archosaurs (birds and crocodiles). All other possible tree topologies were significantly rejected.

Introduction

Unlike other reptiles and birds, turtles have no temporal holes in their skull. Most morphologists traditionally believed that turtles are the only survivors of anapsids, a primitive reptile group lacking temporal holes (e.g., Caroll 1988; Lee 1997), and that turtles are a phylogenetically isolated group in living reptiles, like figure 1_a_ (hypothesis a). A recent morphological analysis, however, did not support the traditional view and placed turtles as the sister group to lepidosaurs (squamates [lizards and snakes] and tuatara), as shown in figure 1_b_ (hypothesis b), with a moderate support value (Rieppel and deBraga 1996; deBraga and Rieppel 1997).

Four possible hypotheses about the phylogenetic position of turtles. (a) Caroll (1988) and Lee (1997). (b) Rieppel and deBraga (1996) and deBraga and Rieppel (1997). (c) Zardoya and Meyer (1998), Hedges and Poling (1999), Kumazawa and Nishida (1999), and Rest et al. (2003). (d) Hedges and Poling (1999) and Mannen and Li (1999).

Recent molecular phylogenetic studies showed the archosaurian (birds and crocodiles) affinities of turtles (Zardoya and Meyer 1998; Hedges and Poling 1999; Kumazawa and Nishida 1999; Mannen and Li 1999; Rest et al. 2003), as shown in figure 1_c_ and d (hypotheses c and d). Mitochondrial data tended to support hypothesis c, whereas nuclear data favored hypothesis d. Using both mitochondrial and nuclear data, Cao et al. (2000) conducted an extensive analysis and concluded that hypotheses a and b are significantly rejected and that more data are needed to discriminate between hypotheses c and d.

Multiple nuclear DNA–coded protein data were considered to be more reliable than mitochondrial data in recovering phylogenetic relationships among major groups of vertebrates (Hedges 2001; Kikugawa et al. 2004; Takezaki et al. 2004). Even when each single gene gives an ambiguous result, a statistically solid inference is possibly obtained by using multiple orthologous gene sequences. Considerable amount of amino acid sequence data are currently available for inferring phylogenetic position of turtles, as shown in the supplementary table. These genes are collectively referred to as “previously available” genes in this paper. Each of these genes is, however, generally short in amino acid alignment length (138–372 aa). We could not completely exclude the possibility of paralogous comparison for these genes because some of them, such as globin and lactate dehydrogenase, have multiple copies in a vertebrate genome. In order to avoid paralogy, large and single-copy genes should be used for phylogenetic inference.

To infer the phylogenetic position of turtles based on statistically reliable data, we have cloned and sequenced two large genes encoding the catalytic subunit of DNA polymerase α (DPLA) and glycinamide ribonucleotide synthetase–aminoimidazole ribonucleotide synthetase–glycinamide ribonucleotide formyltransferase (GAG) from amniotes and an amphibian. These genes are single copy in the human, mouse, and fugu genomes. The amino acid sequences coded by these genes of human, mouse, chicken, caiman (crocodile), iguana (squamate), turtle, and axolotl (amphibian) were subjected to a phylogenetic analysis based on the maximum likelihood (ML) method.

The ML tree inferred from the concatenated alignment of DPLA and GAG (DPLA + GAG; 2,195 aa) is shown in figure 2, in which turtles are the sister group to the monophyletic cluster of chicken and caiman, corresponding to hypothesis c (fig. 1_c_). All other tree topologies, including those corresponding to hypotheses a, b, and d, were significantly rejected (P(KH) [P value based on two-sided Kishino-Hasegawa test] < 0.010, P(AU) [P value based on approximately unbiased test] < 0.044), as shown in table 1. The same tree topology was obtained independently from each of the GAG and DPLA data sets, although not statistically significant; the minimum log-likelihood differences between the ML tree and the second best tree were 12.3 ± 7.71 (P(KH) = 0.060, P(AU) = 0.072) and 8.02 ± 9.38 (P(KH) = 0.20, P(AU) = 0.39) for DPLA and GAG, respectively.

The ML tree inferred from DPLA and GAG. Branch lengths were calculated from the concatenated alignment of DPLA and GAG. RELL BP values calculated from the concatenated alignment/DPLA/GAG are shown in this order for each branch.

FIG. 2.—

Table 1

Comparison of Log-Likelihood Values Based on the DPLA and GAG Proteins

| | | | | | | | | DPLA+GAG (2,195 aa) | | | | | | | | --------------------------------------------- | ------------- | ------------------------------------------------ | ----------------- | ------------ | ------------------------------------------------ | ------- | ------------------- | ------------------------------------------------ | ------- | ------------ | ------------------------------------------------ | ------- | ------- | | DPLA (1,339 aa) | GAG (856 aa) | Concatenated Alignment Analysis | Separate Analysis | | | | | | | | | | | | Treea | | log Lb | P(AU) | P(KH) | log Lb | P(AU) | P(KH) | log Lb | P(AU) | P(KH) | log Lb | P(AU) | P(KH) | | a | (((B,C),S),T) | −12.3 ± 7.71 | 0.072 | 0.060 | −10.7 ± 5.65 | 0.059 | 0.039 | −23.5 ± 9.56 | 0.012 | 0.006 | −23.1 ± 9.55 | 0.019 | 0.011 | | b | ((B,C),(T,S)) | −14.1 ± 7.26 | 0.027 | 0.032 | −10.8 ± 5.67 | 0.069 | 0.040 | −24.6 ± 9.34 | 0.001 | 0.010 | −24.8 ± 9.21 | 0.009 | 0.004 | | c | (((B,C),T),S) | 〈−9,520.7〉 | — | — | 〈−5,828.0〉 | — | — | 〈−15,401.0〉 | — | — | 〈−15,348.7〉 | — | — | | d | ((B,(T,C)),S) | −19.6 ± 7.82 | 0.011 | 0.011 | −8.02 ± 9.38 | 0.389 | 0.195 | −28.5 ± 11.4 | 0.007 | 0.007 | −27.6 ± 12.2 | 0.013 | 0.015 | | (((B,T),C),S) | −19.3 ± 7.94 | 0.007 | 0.014 | −13.8 ± 7.80 | 0.041 | 0.039 | −31.7 ± 10.7 | <0.001 | <0.001 | −33.0 ± 11.1 | 0.001 | 0.003 | | | (((B,T),S),C) | −36.2 ± 10.9 | <0.001 | <0.001 | −23.8 ± 10.3 | 0.017 | 0.017 | −60.1 ± 14.6 | <0.001 | <0.001 | −60.0 ± 15.0 | <0.001 | <0.001 | | | ((B,T),(C,S)) | −33.9 ± 11.4 | 0.005 | 0.001 | −24.2 ± 10.3 | <0.001 | 0.016 | −57.4 ± 15.1 | <0.001 | <0.001 | −58.1 ± 15.3 | 0.021 | <0.001 | | | (B,((T,C),S)) | −28.3 ± 12.8 | 0.055 | 0.017 | −18.1 ± 11.9 | 0.049 | 0.059 | −46.9 ± 17.0 | 0.013 | 0.004 | −46.4 ± 17.5 | 0.006 | 0.004 | | | ((B,S),(T,C)) | −36.7 ± 10.9 | 0.045 | <0.001 | −19.6 ± 11.5 | 0.004 | 0.046 | −56.6 ± 15.3 | 0.044 | <0.001 | −56.3 ± 15.9 | <0.001 | <0.001 | | | ((B,(C,S)),T) | −24.5 ± 11.1 | 0.026 | 0.015 | −26.8 ± 10.0 | 0.049 | 0.009 | −50.5 ± 14.8 | 0.018 | 0.001 | −51.3 ± 14.9 | 0.021 | <0.001 | | | (((B,S),C),T) | −26.1 ± 10.7 | <0.001 | 0.008 | −26.6 ± 10.1 | <0.001 | 0.011 | −53.3 ± 14.2 | 0.029 | <0.001 | −52.7 ± 14.7 | 0.013 | <0.001 | | | (B,(T,(C,S))) | −24.7 ± 12.6 | 0.063 | 0.030 | −19.5 ± 9.40 | 0.119 | 0.021 | −46.9 ± 17.2 | 0.004 | 0.003 | −49.8 ± 17.1 | 0.002 | 0.002 | | | (B,((T,S),C)) | −25.9 ± 13.2 | 0.049 | 0.028 | −23.9 ± 10.9 | 0.029 | 0.019 | −41.8 ± 16.0 | 0.023 | 0.007 | −44.2 ± 15.7 | 0.008 | 0.003 | | | (((B,S),T),C) | −37.5 ± 10.7 | <0.001 | <0.001 | −25.3 ± 10.5 | <0.001 | 0.013 | −64.3 ± 14.3 | <0.001 | <0.001 | −62.8 ± 15.0 | 0.008 | <0.001 | | | | ((B,(T,S)),C) | −33.2 ± 10.6 | <0.001 | 0.004 | −21.1 ± 8.84 | 0.019 | 0.010 | −54.5 ± 13.6 | <0.001 | <0.001 | −54.3 ± 13.8 | <0.001 | <0.001 | |

NOTE.—The taxon names are B, bird; C, crocodile; T, turtle; and S, squamate. Only 15 topologies clustering birds and reptiles relative to mammals and amphibians are shown because both the proteins strongly rejected the other tree topologies (data not shown). The highest log-likelihood values are in boldface.

Trees a–d correspond to hypotheses a_–_d in figure 1.

The log-likelihood values of the highest likelihood trees are shown in angle brackets, and the log-likelihood difference from the ML tree and standard error are shown for each of the alternative trees.

Table 1

Comparison of Log-Likelihood Values Based on the DPLA and GAG Proteins

Trees a–d correspond to hypotheses a_–_d in figure 1.

When the DPLA and GAG sequences were concatenated together with the previously available data (total length was 5,189 aa), hypothesis c was strongly supported, whereas all other tree topologies were significantly rejected (P(KH) < 0.026, P(AU) < 0.032), as shown in table 2.

Table 2

Comparison of Log-Likelihood Values Based on GAG, DPLA, and Previously Available Data

| | | Previously Available Data (2,994 aa) | DPLA + GAG + Previously Available Data (5,189 aa) | | | | | | | | | | ------------------------------- | ------------------------------------ | ------------------------------------------------- | ------- | ------------ | --------------- | ------- | ------------ | --------------- | ------- | ------- | | Concatenated Alignment Analysis | Concatenated Alignment Analysis | Separate Analysis | | | | | | | | | | Treea | | log _L_b | P(AU) | P(KH) | log _L_b | P(AU) | P(KH) | log _L_b | P(AU) | P(KH) | | a | (((B,C),S),T) | −1.69 ± 9.03 | 0.567 | 0.563 | −23.6 ± 12.9 | 0.032 | 0.026 | −30.8 ± 14.9 | 0.021 | 0.025 | | b | ((B,C),(T,S)) | −12.3 ± 6.23 | <0.001 | 0.031 | −34.9 ± 10.9 | 0.001 | <0.001 | −43.2 ± 12.9 | <0.001 | 0.001 | | c | (((B,C),T),S) | 〈−19,591.2〉 | — | — | 〈−35,028.5〉 | — | — | 〈−34,342.5〉 | — | — | | d | ((B,(T,C)),S) | −7.26 ± 10.9 | 0.342 | 0.249 | −34.6 ± 15.5 | 0.022 | 0.021 | −34.5 ± 18.4 | 0.031 | 0.029 | | (((B,T),C),S) | −9.70 ± 10.6 | 0.347 | 0.174 | −39.4 ± 14.9 | 0.012 | 0.004 | −46.5 ± 17.0 | 0.011 | 0.004 | | | (((B,T),S),C) | −18.4 ± 13.8 | 0.118 | 0.085 | −75.4 ± 19.7 | <0.001 | <0.001 | −88.3 ± 22.9 | 0.015 | <0.001 | | | ((B,T),(C,S)) | −17.2 ± 14.0 | 0.133 | 0.116 | −69.6 ± 20.3 | <0.001 | <0.001 | −83.5 ± 23.9 | 0.012 | <0.001 | | | (B,((T,C),S)) | −17.6 ± 14.1 | 0.031 | 0.100 | −62.6 ± 21.4 | 0.007 | 0.001 | −70.8 ± 25.1 | <0.001 | 0.001 | | | ((B,S),(T,C)) | −12.4 ± 14.7 | 0.173 | 0.195 | −68.2 ± 20.6 | 0.015 | <0.001 | −78.6 ± 24.0 | 0.012 | 0.001 | | | ((B,(C,S)),T) | −7.40 ± 13.5 | 0.286 | 0.302 | −52.9 ± 19.7 | 0.009 | 0.003 | −65.4 ± 23.0 | 0.007 | 0.001 | | | (((B,S),C),T) | −3.42 ± 14.2 | 0.626 | 0.427 | −53.5 ± 19.7 | 0.004 | 0.003 | −64.3 ± 22.5 | 0.013 | 0.003 | | | (B,(T,(C,S))) | −21.7 ± 13.9 | 0.016 | 0.058 | −65.8 ± 21.5 | <0.001 | <0.001 | −77.4 ± 25.3 | 0.008 | 0.001 | | | (B,((T,S),C)) | −29.0 ± 12.3 | 0.015 | 0.004 | −69.3 ± 19.7 | <0.001 | <0.001 | −82.2 ± 22.2 | 0.007 | <0.001 | | | (((B,S),T),C) | −15.5 ± 14.7 | 0.133 | 0.139 | −78.6 ± 20.0 | 0.004 | <0.001 | −87.9 ± 23.3 | 0.009 | <0.001 | | | | ((B,(T,S)),C) | −28.4 ± 12.4 | 0.021 | 0.009 | −82.4 ± 17.9 | <0.001 | <0.001 | −94.1 ± 20.9 | 0.007 | <0.001 | |

NOTE.—See the footnotes of Table 1.

Table 2

Comparison of Log-Likelihood Values Based on GAG, DPLA, and Previously Available Data

NOTE.—See the footnotes of Table 1.

We also carried out separate analyses, in which the log-likelihood values were separately calculated for 15 proteins and then summed, and obtained similar results to those from concatenated alignment analyses (the right three columns of tables 1 and 2).

We analyzed the concatenated alignment of previously available amino acid sequences (2,994 aa), although the possibility of paralogy could not be completely excluded. This data also supported hypothesis c, but the difference from the second best tree was only 1.69 ± 9.03 (P(KH) = 0.56, P(AU) = 0.57). This result disagreed with Hedges and Poling (1999), which supported hypothesis d based on the nuclear DNA data available at that time. This discrepancy was probably caused by the amount of available data and/or unrecognized paralogous comparisons.

In summary, we conducted a phylogenetic analysis of amniotes based on large amount of nuclear DNA–coded protein sequence data and obtained a single tree topology supporting hypothesis c, in which turtles are the sister group to a monophyletic cluster of archosaurs (birds and crocodiles). All other tree topologies, including the traditional one, were significantly rejected.

Materials and Methods

Total RNAs were extracted from liver of Caiman crocodilus (spectacled caiman), tail of Trachemys scripta (red-eared slider), tail of Iguana iguana (green iguana), tail of Ambystoma mexicanum (Mexican axolotl), and embryo of Gallus gallus (chicken) using TRIZOL reagent (Invitrogen, Carlsbad, Calif.). These total RNAs were reverse-transcribed to cDNAs using SMART RACE cDNA Amplification Kit (Clontech, Palo Alto, Calif.). These cDNAs were used as templates for PCR amplification with Expand High-Fidelity PCR System (Roche, Basel, Switzerland). The sense and antisense degenerate primers were designed from conserved amino acid residues of each gene. As fragments of chicken DPLA gene were found in the dbEST database of the National Center for Biotechnology Information, gene specific primers were used for this gene. The PCR products were purified and cloned into the pT7Blue vector (Novagen, Darmstadt, Germany). More than three independent clones were isolated and sequenced using ABI 3100 DNA Sequencer (Applied Biosystems, Foster City, Calif.). The 3′ ends of cDNAs were amplified using 3′RACE (GIBCO BRL, Invitrogen) and sequenced in the same way as above. The following sequence data were taken from the GenBank database: human and mouse DPLA; human, mouse, and chicken GAG; and previously available data listed in supplementary table.

Tuatara, one of the four major groups of reptiles (crocodiles, turtles, squamates, and tuatara), is thought to be closely related to squamates based on morphological data (Caroll 1988) and mitochondrial sequence data (Rest et al. 2003). Thus, tuatara was excluded from the present analysis.

Chicken has a pair of duplicated GAG genes (GAG-A and -B) (Smith et al. 2000). Chicken GAG-B was estimated to be more closely related to chicken GAG-A than to any of the reptile GAG genes, in distance measured by synonymous substitutions (Miyata and Yasunaga 1980). According to a comparison of amino acid substitution rates of amniote GAG genes using amphibian sequences as an out-group, the chicken GAG-B gene was estimated to evolve approximately 3.5–4.5 times faster than the others. These observations suggest that the GAG gene was duplicated on the avian lineage and that one (GAG-B) of the duplicated genes accumulated amino acid changes at an extremely rapid rate. Thus, the chicken GAG-B gene was excluded from the present analysis, in order to avoid the long-branch attraction artifact (Felsenstein 1978).

Each protein data set was multiply aligned by MAFFT (Katoh et al. 2002) and manually inspected. Unambiguously aligned amino acid positions were subjected to phylogenetic tree analyses based on the ML method (Felsenstein 1981; Kishino, Miyata, and Hasegawa 1990). For each of 945 possible topologies consisting of seven taxonomic groups (six amniote groups and an out-group), the log-likelihood values based on GAG (856 aa), DPLA (1,339 aa), each of the previously available 13 proteins, and three types of concatenated alignments (DPLA + GAG, 2,195 aa; DPLA + GAG + previously available data, 5,189 aa; and previously available data only, 2,994 aa) were calculated using the GAMT (Katoh, Kuma, and Miyata 2001) program, assuming the JTT-F model (Jones, Taylor, and Thornton 1992; Cao et al. 1994; Adachi and Hasegawa 1996). Heterogeneity of evolutionary rates among sites was modeled by a discrete Γ distribution (Yang 1994) with the shape parameter α optimized for each data set.

Two-sided KH test and the AU test were carried out using the CONSEL package (Shimodaira and Hasegawa 2001). Bootstrap probability (BP) based on the resampling of estimated log-likelihoods (RELL) approximation (Kishino, Miyata, and Hasegawa 1990) was also computed. The bootstrap probability for a cluster is calculated by totaling the RELL BP values of the tree topologies having the cluster.

Present address: Division of Bioengineering and Bioinformatics, Graduate School of Information Science and Technology, Hokkaido University, Sapporo, Japan.

Present address: Division of Morphogenesis, Department of Developmental Biology, National Institute for Basic Biology, Okazaki, Japan.

Present address: Evolutionary Regeneration Biology Group, RIKEN Center for Developmental Biology, Kobe, Japan.

Present addresses: Biohistory Research Hall, Takatsuki, Osaka, Japan, and Department of Electrical Engineering and Bioscience, Science and Engineering, Waseda University, Tokyo, Japan.

Present address: Bioinformatics Center, Institute for Chemical Research, Kyoto University, Uji, Kyoto, Japan.

Takashi Gojobori, Associate Editor

We thank K. Kuma for comments and suggestions. This work was supported in part by a Grant-in-Aid for Creative Scientific Research and a Grant for the Biodiversity Research of the 21st Century COE (A14) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.

References

Adachi, J., and M. Hasegawa.

1996

. MOLPHY Version 2.3: programs for molecular phylogenetics based on maximum likelihood.

Comput. Sci. Monogr.

–150.

Cao, Y., J. Adachi, A. Janke, S. Paabo, and M. Hasegawa.

1994

. Phylogenetic relationships among eutherian orders estimated from inferred sequences of mitochondrial proteins: instability of a tree based on a single gene.

J. Mol. Evol.

519

–527.

Cao, Y., M. D. Sorenson, Y. Kumazawa, D. P. Mindell, and M. Hasegawa.

2000

. Phylogenetic position of turtles among amniotes: evidence from mitochondrial and nuclear genes.

Gene

259

139

–148.

Caroll, R. L.

1988

. Vertebrate paleontology and evolution. Freeman, New York.

deBraga, M., and O. Rieppel.

1997

. Reptile phylogeny and the interrelationships of turtles.

Zool. J. Linn. Soc.

120

281

–354.

Felsenstein, J.

1978

. Cases in which parsimony or compatibility methods will be positively misleading.

Syst. Zool.

401

–410.

———.

1981

. Evolutionary trees from DNA sequences: a maximum likelihood approach.

J. Mol. Evol.

368

–376.

Hedges, S. B.

2001

. Molecular evidence for the early history of living vertebrates. Pp. 119–134 in P. E. Ahlberg, ed. Major events in early vertebrate evolution. Taylor & Francis, London.

Hedges, S. B., and L. L. Poling.

1999

. A molecular phylogeny of reptiles.

Science

283

998

–1001.

Jones, D. T., W. R. Taylor, and J. M. Thornton.

1992

. The rapid generation of mutation data matrices from protein sequences.

Comput. Appl. Biosci.

275

–282.

Katoh, K., K. Kuma, and T. Miyata.

2001

. Genetic algorithm-based maximum-likelihood analysis for molecular phylogeny.

J. Mol. Evol.

477

–484.

Katoh, K., K. Misawa, K. Kuma, and T. Miyata.

2002

. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform.

Nucleic Acids Res.

3059

–3066.

Kikugawa, K., K. Katoh, S. Kuraku, H. Sakurai, O. Ishida, N. Iwabe, and T. Miyata.

2004

. Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes.

BMC Biol.

Kishino, H., T. Miyata, and M. Hasegawa.

1990

. Maximum likelihood inference of protein phylogeny and the origin of chloroplasts.

J. Mol. Evol.

151

–160.

Kumazawa, Y., and M. Nishida.

1999

. Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink: statistical evidence for archosaurian affinity of turtles.

Mol. Biol. Evol.

784

–792.

Lee, M. S. Y.

1997

. Pareiasaur phylogeny and the origin of turtles.

Zool. J. Linn. Soc.

120

197

–280.

Mannen, H., and S. S. Li.

1999

. Molecular evidence for a clade of turtles.

Mol. Phylogenet. Evol.

144

–148.

Miyata, T., and T. Yasunaga.

1980

. Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application.

J. Mol. Evol.

–36.

Rest, J. S., J. C. Ast, C. C. Austin, P. J. Waddell, E. A. Tibbetts, J. M. Hay, and D. P. Mindell.

2003

. Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome.

Mol. Phylogenet. Evol.

289

–297.

Rieppel, O., and M. deBraga.

1996

. Turtles as diapsid reptiles.

Nature

384

453

–455.

Shimodaira, H., and M. Hasegawa.

2001

. CONSEL: for assessing the confidence of phylogenetic tree selection.

Bioinformatics

1246

–1247.

Smith, J., J. Reboul, G. Lutfalla, and D. W. Burt.

2000

. Human chromosomes 3 and 21 are the products of an ancestral gene arrangement that is at least 300 million years old.

Mamm. Genome

806

–807.

Takezaki, N., F. Figueroa, Z. Zaleska-Rutczynska, N. Takahata, and J. Klein.

2004

. The phylogenetic relationship of tetrapod, coelacanth, and lungfish revealed by the sequences of forty-four nuclear genes.

Mol. Biol. Evol.

1512

–1524.

Yang, Z.

1994

. Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods.

J. Mol. Evol.

306

–314.

Zardoya, R., and A. Meyer.

1998

. Complete mitochondrial genome suggests diapsid affinities of turtles.

Proc. Natl. Acad. Sci. USA

14226

–14231.

Author notes

*Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan; and †Division of Material Science, Graduate School of Science, Nagoya University, Nagoya, Japan

Supplementary data

Citations

Views

Altmetric

Metrics

Total Views 6,529

5,880 Pageviews

649 PDF Downloads

Since 11/1/2016

Month:	Total Views:
November 2016	1
January 2017	3
February 2017	12
March 2017	24
April 2017	20
May 2017	16
June 2017	4
July 2017	3
August 2017	6
September 2017	4
October 2017	21
November 2017	94
December 2017	134
January 2018	135
February 2018	107
March 2018	129
April 2018	143
May 2018	137
June 2018	133
July 2018	113
August 2018	127
September 2018	88
October 2018	115
November 2018	132
December 2018	127
January 2019	90
February 2019	106
March 2019	122
April 2019	139
May 2019	70
June 2019	69
July 2019	68
August 2019	71
September 2019	83
October 2019	95
November 2019	74
December 2019	54
January 2020	51
February 2020	66
March 2020	90
April 2020	174
May 2020	74
June 2020	125
July 2020	70
August 2020	84
September 2020	90
October 2020	119
November 2020	92
December 2020	56
January 2021	47
February 2021	82
March 2021	97
April 2021	87
May 2021	67
June 2021	62
July 2021	32
August 2021	70
September 2021	37
October 2021	28
November 2021	71
December 2021	48
January 2022	33
February 2022	43
March 2022	59
April 2022	71
May 2022	72
June 2022	49
July 2022	30
August 2022	46
September 2022	55
October 2022	69
November 2022	59
December 2022	47
January 2023	40
February 2023	68
March 2023	57
April 2023	77
May 2023	53
June 2023	28
July 2023	40
August 2023	39
September 2023	42
October 2023	84
November 2023	54
December 2023	41
January 2024	65
February 2024	59
March 2024	66
April 2024	84
May 2024	51
June 2024	31
July 2024	50
August 2024	49
September 2024	57
October 2024	65
November 2024	8

Citations

137 Web of Science

Sister Group Relationship of Turtles to the Bird-Crocodilian Clade Revealed by Nuclear DNA–Coded Proteins (original) (raw)

Cite

Abstract

Introduction

Materials and Methods

References

Author notes

Supplementary data

Citations

Views

Altmetric

Email alerts

Email alerts

Citing articles via

Latest

Most Cited

Sister Group Relationship of Turtles to the Bird-Crocodilian Clade Revealed by Nuclear DNA–Coded Proteins (original) (raw)

Cite

Abstract

Introduction

Materials and Methods

References

Author notes

Supplementary data

Citations

Views

Altmetric

Email alerts

Email alerts

Citing articles via

Latest

Most Read

Most Cited