Molecular phylogeography and systematics of the arid-zone members of the Egernia whitii (Lacertilia: Scincidae) species group

We assembled a molecular phylogeny for the arid-zone members of the Egernia whitii species group to test Pianka's [Zoogeography and speciation of Australian desert lizards: an ecological perspective, Copeia (1972) 127-145] hypothesis that habitat specificity to the three major arid-zone vegetation communities is the primary cause of lizard speciation within the arid interior of Australia. This hypothesis predicts that species should exhibit phylogeographic structuring concordant with the major arid-zone vegetation types. Sequence data were obtained from four of the five arid-zone members of the E. whitii species group, and from across the ranges of the ecologically generalized E. inornata and E. multiscutata and the more specialized E. striata. We targeted a fragment (696 base pair (bp)) of the mitochondrial genome comprising the 3' half of the ND4 gene. We analysed the data using parsimony, maximum likelihood and Bayesian methods. Our phylogeny confirms the monophyly of the arid-zone members of the species group, although the phylogenetic relationships among species were not fully resolved. Although our topology does not support the recognition of the existing subspecies within E. multiscutata, there is a substantial phylogeographic break between South Australian/Victorian (Clade 1) and Western Australian (Clade 2) populations. We found considerable phylogeographic structure within E. inornata, with six major clades identified. However, these clades were not concordant with the distribution of habitat types in the arid-zone. Phylogeographic structure was also observed in the more specialized E. striata, although our analysis revealed close phylogenetic affinities between the sympatric species E. striata and E. kintorei. Shimodaira-Hasegawa topology tests were equivocal in regard to whether the phylogeographic structure within E. striata was in accordance with Pianka's predictions. Although our data failed to provide strong support for the suggestion that ecological and habitat factors are responsible for the diversification of arid-zone lizards, most E. inornata and E. striata populations had similar habitats, indicating that adaptation to particular habitats may have some role in the speciation of lizards in the Australian arid-zone.     

Modeling nucleotide evolution at the mesoscale: the phylogeny of the neotropical pitvipers of the Porthidium group (viperidae: crotalinae)

We analyzed the phylogeny of the Neotropical pitvipers within the Porthidium group (including intra-specific through inter-generic relationships) using 1.4 kb of DNA sequences from two mitochondrial protein-coding genes (ND4 and cyt-b). We investigated how Bayesian Markov chain Monte-Carlo (MCMC) phylogenetic hypotheses based on this 'mesoscale' dataset were affected by analysis under various complex models of nucleotide evolution that partition models across the dataset. We develop an approach, employing three statistics (Akaike weights, Bayes factors, and relative Bayes factors), for examining the performance of complex models in order to identify the best-fit model for data analysis. Our results suggest that: (1) model choice may have important practical effects on phylogenetic conclusions even for mesoscale datasets, (2) the use of a complex partitioned model did not produce widespread increases or decreases in nodal posterior probability support, and (3) most differences in resolution resulting from model choice were concentrated at deeper nodes. Our phylogenetic estimates of relationships among members of the Porthidium group (genera: Atropoides, Cerrophidion, and Porthidium) resolve the monophyly of the three genera. Bayesian MCMC results suggest that Cerrophidion and Porthidium form a clade that is the sister taxon to Atropoides. In addition to resolving the intra-specific relationships among a majority of Porthidium group taxa, our results highlight phylogeographic patterns across Middle and South America and suggest that each of the three genera may harbor undescribed species diversity.     

Rostral horn evolution among agamid lizards of the genus Ceratophora endemic to Sri Lanka.

The first phylogenetic hypothesis for the Sri Lankan agamid lizard genus Ceratophora is presented based on 1670 aligned base positions (472 parsimony informative) of mitochondrial DNA sequences, representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Phylogenetic analysis reveals multiple origins and possibly losses of rostral horns in the evolutionary history of Ceratophora. Our data suggest a middle Miocene origin of Ceratophora with the most recent branching of recognized species occurring at the Pliocene/Pleistocene boundary. Haplotype divergence suggests that an outgroup species, Lyriocephalus scutatus, dates at least to the Pliocene. These phylogenetic results provide a framework for comparative studies of the behavioral ecological importance of horn evolution in this group.     

Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of Parsimony, Likelihood, and Bayesian approaches

Phylogenetic analysis of large datasets using complex nucleotide substitution models under a maximum likelihood framework can be computationally infeasible, especially when attempting to infer confidence values by way of nonparametric bootstrapping. Recent developments in phylogenetics suggest the computational burden can be reduced by using Bayesian methods of phylogenetic inference. However, few empirical phylogenetic studies exist that explore the efficiency of Bayesian analysis of large datasets. To this end, we conducted an extensive phylogenetic analysis of the wide-ranging and geographically variable Eastern Fence Lizard (Sceloporus undulatus). Maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses were performed on a combined mitochondrial DNA dataset (12S and 16S rRNA, ND1 protein-coding gene, and associated tRNA; 3,688 bp total) for 56 populations of S. undulatus (78 total terminals including other S. undulatus group species and outgroups). Maximum parsimony analysis resulted in numerous equally parsimonious trees (82,646 from equally weighted parsimony and 335 from weighted parsimony). The majority rule consensus tree derived from the Bayesian analysis was topologically identical to the single best phylogeny inferred from the maximum likelihood analysis, but required approximately 80% less computational time. The mtDNA data provide strong support for the monophyly of the S. undulatus group and the paraphyly of "S. undulatus" with respect to S. belli, S. cautus, and S. woodi. Parallel evolution of ecomorphs within "S. undulatus" has masked the actual number of species within this group. This evidence, along with convincing patterns of phylogeographic differentiation suggests "S. undulatus" represents at least four lineages that should be recognized as evolutionary species.     

Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: evidence for a misleading mitochondrial gene tree

It has proven remarkably difficult to obtain a well-resolved and strongly supported phylogeny for horned lizards (Phrynosoma) because of incongruence between morphological and mitochondrial DNA sequence data. We infer the phylogenetic relationships among all 17 extant Phrynosoma species using >5.1 kb of mtDNA (12S rRNA, 16S rRNA, ND1, ND2, ND4, Cyt b, and associated tRNA genes), and >2.2kb from three nuclear genes (RAG-1, BDNF, and GAPD) for most taxa. We conduct separate and combined phylogenetic analyses of these data using maximum parsimony, maximum likelihood, and Bayesian methods. The phylogenetic relationships inferred from the mtDNA data are congruent with previous mtDNA analyses based on fewer characters and provide strong support for most branches. However, we detected strong incongruence between the mtDNA and nuclear data using comparisons of branch support and Shimodaira-Hasegawa tests, with the (P. platyrhinos+P. goodei) clade identified as the primary source of this conflict. Our analysis of a P. mcalliixP. goodei hybrid suggests that this incongruence is caused by reticulation via introgressive hybridization. Our preferred phylogeny based on an analysis of the combined data (excluding the introgressed mtDNA data) provides a new framework for interpreting character evolution and biogeography within Phrynosoma. In the context of this improved phylogeny we propose a phylogenetic taxonomy highlighting four clades: (1) Tapaja, containing the viviparous short-horned lizards P. ditmarsi, P. hernandesi, P. douglasii, and P. orbiculare; (2) Anota, containing species with prominent cranial horns (P. solare, P. mcallii, and the P. coronatum group); (3) Doliosaurus, containing three species lacking antipredator blood-squirting (P. modestum, P. platyrhinos, and P. goodei); and (4) Brevicauda, containing two viviparous species with extremely short tails that lack blood-squirting (P. braconnieri and P. taurus).     

Comparative phylogeography of two open forest frogs from eastern Australia

We investigated the phylogeography of two closely related Australian frog species from open forest habitats, Limnodynastes tasmaniensis and L. peronii , using mitochondrial ND4 sequence data. Comparison of our results with previous work on Litoria fallax allowed us to test the generality of phylogeographic patterns among non-rainforest anurans along the east coast of Australia. In general, there was no strong evidence for congruence between overall patterns of genetic structure in the three species. However, phylogenetic breaks congruent with the position of the Burdekin Gap were detected at some level in all species. As previously noted for closed forest taxa, this area of dry habitat appears to have been an important influence on the evolution of several open forest taxa. There were broad geographic similarities in the phylogenetic structuring of southern populations of L. peronii and L. tasmaniensis. Contrarily, although the McPherson Range has previously been noted to coincide geographically with a major mtDNA phylogenetic break in Litoria fallax this pattern is not apparent in L. peronii or L. tasmaniensis. It appears that major phylogeographic splits within L. peronii and L. tasmaniensis may predate the Quaternary. We conclude that phylogeographies of open forest frogs are complex and more difficult to predict than for rainforest taxa, mainly due to an absence of palæomodels for historical distributions of non-rainforest habitats.     

A molecular phylogenetic study of ecological diversification in the Australian lizard genus Ctenophorus.

We present phylogenetic analyses of the lizard genus Ctenophorus using 1,639 aligned positions of mitochondrial DNA sequences containing 799 parsimony-informative characters for samples of 22 species of Ctenophorus and 12 additional Australian agamid genera. Sequences from three protein-coding genes (ND1, ND2, and COI) and eight intervening tRNA genes are examined using both parsimony and maximum-likelihood analyses. Species of Ctenophorus form a monophyletic group with Rankinia adelaidensis, which we suggest placing in Ctenophorus. Ecological differentiation among species of Ctenophorus is most evident in the kinds of habitats used for shelter. Phylogenetic analyses suggest that the ancestral condition is to use burrows for shelter, and that habits of sheltering in rocks and shrubs/hummock grasses represent separately derived conditions. Ctenophorus appears to have undergone extensive cladogenesis approximately 10-12 million years ago, with all three major ecological modes being established at that time.     

Molecular phylogeny of the lemur family cheirogaleidae (primates) based on mitochondrial DNA sequences

Cheirogaleidae currently comprises five genera whose relationships remain contentious. The taxonomic status and phylogenetic position of both Mirza coquereli and Allocebus trichotis are still unclear. The taxonomic status of the recently discovered Microcebus ravelobensis (a sympatric sibling species of Microcebus murinus) and its phylogenetic position also require further examination. A approximately 2.4-kb mitochondrial DNA sequence including part of the COIII gene, complete ND3, ND4L, and ND4 genes, and 5 tRNAs was used to clarify relationships among cheirogaleids. Mirza and Microcebus form a clade representing the sister group of Allocebus, with a clade containing Cheirogaleus major and Cheirogaleus medius diverging first. M. ravelobensis and Microcebus rufus form a subclade within Microcebus, with M. murinus as its sister group. The molecular data support the generic status of Mirza coquereli and species-level divergence of M. ravelobensis. Furthermore, "M. rufus" may well represent more than one species.     

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.     

Molecular phylogenetics and evolution of turtles

Turtles are one of Earth's most instantly recognizable life forms, distinguished for over 200 million years in the fossil record. Even so, key nodes in the phylogeny of turtles remain uncertain. To address this issue, we sequenced >90% of the nuclear recombination activase gene 1 (RAG-1) for 24 species representing all modern turtle families. RAG-1 exhibited negligible saturation and base composition bias, and extensive base composition homogeneity. Most of the relationships suggested by prior phylogenetic analyses were also supported by RAG-1 and, for at least two critical nodes, with a much higher level of support. RAG-1 also indicates that the enigmatic Platysternidae and Chelydridae, often considered sister taxa based on morphological evidence, are not closely related, although their precise phylogenetic placement in the turtle tree is still unresolved. Although RAG-1 is phylogenetically informative, our research revealed fundamental conflicts among analytical methods for estimating phylogenetic hypotheses. Maximum parsimony analyses of RAG-1 alone and in combination with two mitochondrial genes suggest the earliest phylogenetic splits separating into three basal branches, the pig-nosed turtles (Carettochelyidae), the softshell turtles (Trionychidae), and a clade comprising all remaining extant turtles. Maximum likelihood and Bayesian analyses group Carettochelyidae and Trionychidae (=Trionychoidae) in their more traditional location as the sister taxon to all other hidden-necked turtles, collectively forming the Cryptodira. Our research highlights the utility of molecular data in identifying issues of character homology in morphological datasets, while shedding valuable light on the biodiversity of a globally imperiled taxon.     

A phylogenetic analysis of the frog genera Vibrissaphora and Leptobrachium, and the correlated evolution of nuptial spine and reversed sexual size dimorphism

A phylogeny of the frog genus Vibrissaphora and seven Leptobrachium species was reconstructed based on 3808 base pairs of mitochondrial and nuclear DNA sequence data. Maximum parsimony, Bayesian, and statistical parsimony approaches were employed to reveal the historical relationships among the recovered haplotypes. The species of the genera Vibrissaphora and Leptobrachium failed to form their respective monophyletic groups. Therefore, Vibrissaphora should be part of the genus Leptobrachium. Our analysis also further delineated several species boundaries; Leptobrachium chapense is not a single species, rather, a species complex. On the other hand, although L. liui demonstrated substantial morphological variation, all populations belong to one species. The RAG-1 data revealed a case of species-level non-monophyly; when studying slowly evolving genes, population coalescence may not have been reached within species. The phylogenetic comparative analysis revealed a positive correlation between evolution of male nuptial spines and reversed sexual size dimorphism in this group. Resource defense polygyny and male parental care are possible explanations for the presence of large male size.     

中华雏蝗（Chorthippus chinensis Tarb）线粒体基因组分析

采用Lon-PCR扩增线粒体全基因组和保守引物步移法结合克隆方法测定并拼接获得了中华雏蝗（Chorthippus chinensis Tarb）线粒体基因组全序列．序列的注释和分析结果表 明,中华雏蝗线粒体基因组序列全长15 599 bp,共有13个编码蛋白质基因、22个tRNA基因、2个rRNA基因和1个A+T富集区．基因顺序与非洲飞蝗（Locusta migratoria）相同,也发生了2个 tRNA Asp(D)和tRNALys(K)的倒置．13个编码蛋白质基因都使用了ATN作为起始密码子．除ND1以TAG和ND5的终止密码子为不完全的T外,其余11个编码蛋白质基因的终止密码子都为完整的TAA．6种直翅类昆虫13个蛋白质的氨基酸序列的联合数据集构建的系统树与形态分类系统一致,中华雏蝗与非洲飞蝗为姐妹群,并与东方蝼蛄构成一单系群．     

Phylogeny and phylogenetic classification of the antbirds, ovenbirds, woodcreepers, and allies (Aves: Passeriformes: infraorder Furnariides)

The infraorder Furnariides is a diverse group of suboscine passerine birds comprising a substantial component of the Neotropical avifauna. The included species encompass a broad array of morphologies and behaviours, making them appealing for evolutionary studies, but the size of the group (ca. 600 species) has limited well-sampled higher-level phylogenetic studies. Using DNA sequence data from the nuclear RAG-1 and RAG-2 exons, we undertook a phylogenetic analysis of the Furnariides sampling 124 (more than 88%) of the genera. Basal relationships among family-level taxa differed depending on phylogenetic method, but all topologies had little nodal support, mirroring the results from earlier studies in which discerning relationships at the base of the radiation was also difficult. In contrast, branch support for family-rank taxa and for many relationships within those clades was generally high. Our results support the Melanopareidae and Grallariidae as distinct from the Rhinocryptidae and Formicariidae, respectively. Within the Furnariides our data contradict some recent phylogenetic hypotheses and suggest that further study is needed to resolve these discrepancies. Of the few genera represented by multiple species, several were not monophyletic, indicating that additional systematic work remains within furnariine families and must include dense taxon sampling. We use this study as a basis for proposing a new phylogenetic classification for the group and in the process erect new family-group names for clades having high branch support across methods.     

Out of Africa? Phylogenetic relationships between Falco biarmicus and the other hierofalcons (Aves: Falconidae)

The phylogeographic history of the lanner falcon ( Falco biarmicus ) and the phylogenetic relationships among hierofalcons ( F. biarmicus , Falco cherrug , Falco jugger and Falco rusticolus ) were investigated using mitochondrial (mt) DNA sequences. Of the two non-coding mt sections tested, the control region (CR) appeared more suitable as phylogenetic marker sequence compared with the pseudo control region (ΨCR). For the comprehensive analysis samples from a broad geographic range representing all four hierofalcon species and their currently recognized subspecies were included. Moreover, samples of Falco mexicanus were analysed to elucidate its phylogenetic relationships to the hierofalcons. The sequence data indicate that this species is more closely related to Falco peregrinus than to the hierofalcons. In the DNA-based trees and in the maximum parsimony network all hierofalcons appear closely related and none of the species represents a monophyletic group. The close relationships among haplotypes suggest that the hierofalcon complex is an assemblage of morphospecies not yet differentiated in the genetic markers used in the present study and that the radiation of the four hierofalcon species took place rather recently. Based on the high intraspecific diversity found within F. biarmicus we assume an African origin of the hierofalcon complex. The observed pattern of haplotype distribution in the extant species may be due to incomplete lineage sorting of ancestral polymorphisms, and interspecific gene flow through hybridization.     

Species richness and endemism in the Western Australian flora

Aim Estimates of endemic and non‐endemic native vascular plant species in each of the three Western Australian Botanical Provinces were made by East in 1912 and Beard in 1969. The present paper contains an updated assessment of species endemism in the State. Location Western Australia comprises one third of the continental Australian land mass. It extends from 13° to 35° S and 113° to 129° W. Methods Western Australia is recognized as having three Botanical Provinces (Northern, Eremaean and South‐West) each divided into a number of Botanical Districts. Updated statistics for number of species and species endemism in each Province are based on the Census of Western Australian Plants data base at the Western Australian Herbarium ( Western Australian Herbarium, 1998 onwards). Results The number of known species in Western Australia has risen steadily over the years but reputed endemism has declined in the Northern and Eremaean Provinces where cross‐continental floras are common. Only the isolated South‐West Province retains high rates of endemism (79%). Main conclusions With 5710 native species, the South‐West Province contains about the same number as the California Floristic Province which has a similar area. The Italian mediterranean zone also contains about this number but in a smaller area, while the much smaller Cape Floristic Region has almost twice as many native species. The percentage of endemic species is highest at the Cape, somewhat less in south‐western Australia and less again in California. Italy, at 12.5%, has the lowest value. Apart from Italy, it is usual for endemism to reach high values in the largest plant families. In Western Australia, these mainly include woody sclerophyll shrubs and herbaceous perennials with special adaptations to environmental conditions. While those life forms are prominent in the Cape, that region differs in the great importance of herbaceous families and succulents, both of which are virtually absent from Western Australia. In California and Italy, most endemics are in families of annual, herbaceous perennial and soft shrub plants. It is suggested that the dominant factor shaping the South‐West Province flora is the extreme poverty of the area’s soils, a feature that emphasizes sclerophylly, favours habitat specialization and ensures relatively many local endemic species.     

Genetic elucidation of cryptic and ancient diversity in a group of Australian diplodactyline geckos: the Diplodactylus vittatus complex

We examine species boundaries and phylogenetic relationships within the Australian diplodactyline geckos currently assigned to Diplodactylus granariensis and Diplodactylus vittatus using data from multiple allozyme loci, the mitochondrial ND2 gene and published karyotypic studies. These analyses uncover at least six morphologically cryptic but genetically distinctive species. The taxa identified correspond poorly with the existing taxonomy and represent a paraphyletic assemblage with respect to a number of other Diplodactylus species included in analyses. Molecular dating indicates that the species identified evolved considerably before the Pleistocene, and that the basal node of the D. vittatus species group (as redefined here) is of comparable age to entire radiations of other Australian squamate families. The antiquity of this one small group within the diplodactyline geckos suggests that further study of patterns of diversification in the diplodactylines will provide both a valuable insight into long-term patterns of environmental change in the Australian continent and a useful contrast for analysis of patterns of radiation in other Australian lizard groups.     

Evolution of Galapagos Island Lava Lizards (Iguania: Tropiduridae: Microlophus)

Nucleotide sequences of mitochondrial genes (ND1, ND2, COI, and tRNAs) were determined for 38 samples representing 15 taxa of tropidurid lizards from the Galapagos Islands and mainland South America. Phylogenetically informative characters (759 of 1,956) were analyzed under Bayesian, maximum likelihood, and parsimony frameworks. This study supports the hypothesis that tropidurid lizards dispersed to the Galapagos on at least two separate occasions. One dispersal event involved an eastern Galapagos clade (Microlophus habelii and M. bivittatus, on Marchena and San Cristobal islands, respectively) the sister taxon of which is M. occipitalis from coastal Ecuador and Peru; the closest mainland relative of the western Galapagos clade was not unambiguously identified. The wide-ranging M. albemarlensis is revealed to be a complex of weakly divergent lineages that is paraphyletic with respect to the insular species M. duncanensis, M. grayii, and M. pacificus.     

Phylogeny and evolution of the Drosophila nasuta subgroup based on mitochondrial ND4 and ND4L gene sequences

The sequences of the mitochondrial ND4 gene (1339 bp) and the ND4L gene (290 bp) were determined for all the 14 extant taxa of the Drosophila nasuta subgroup. The average A + T content of ND4 genes is 76.5% and that of ND4L genes is 83.5%. A total of 114 variable sites were scored. The ND4 gene sequence divergence ranged from 0 to 5.4% within the subgroup. The substitution rate of the ND4 gene is about 1.25% per million years. The base substitution of the genes is strongly transition biased. Neighbor-joining and parsimony were used to construct a phylogeny based on the resultant sequence data set. According to these trees, five distinct mtDNA clades can be identified. D. niveifrons represents the most diverged lineage. D. sulfurigaster bilimbata and D. kepulauana form two independent lineages. The other two clades are the kohkoa complex and the albomicans complex. The kohkoa complex consists of D. sulfurigaster sulfurigaster, D. pulaua, D. kohkoa, and Taxon-F. The albomicans complex can be divided into two groups: D. nasuta, D. sulfurigaster neonasuta, D. sulfurigaster albostrigata, and D. albomicans from Chiangmai form one group; and D. pallidifrons, Taxon-I, Taxon-J, and D. albomicans from China form the other group. High genetic differentiation was found among D. albomicans populations. Based on our phylogenetic results, we hypothesize that D. niveifrons diverged first from the D. nasuta subgroup in Papua New Guinea about 3.5 Mya. The ancestral population spread to the north and when it reached Borneo, it diversified sequentially into the kohkoa complex, D. s. bilimbata, and D. kepulauana. About 1 Mya, another radiation occurred when the ancestral populations reached the Indo-China Peninsula, forming the albomicans complex. Discrepancy between morphological groupings and phylogenetic results suggests that the male morphological traits may not be orthologous.     

Phylogenetic relationships within the genus Pimephales as inferred from ND4 and ND4L nucleotide sequences

In phylogenetic analyses, conducted on the ND4L gene and part of the ND4 gene from species of the genus Pimephales , maximum parsimony yielded four trees, with the strict consensus providing no resolution of relationships among species. Maximum likelihood and minimum evolution methods yielded identical tree topologies, which differed from previous hypotheses of relationships for these species. If this topology is correct, it implies independent evolution of morphological characters, possibly associated with convergent trophic specialization.     

Molecular phylogenetics of emydine turtles: taxonomic revision and the evolution of shell kinesis

The 10 extant species of emydine turtles represent an array of morphological and ecological forms recognizable and popular among scientists and hobbyists. Nevertheless, the phylogenetic affinities of most emydines remain contentious. Here, we examine the evolutionary relationships of emydine turtles using 2092 bp of DNA encoding the mitochondrial genes cyt b, ND4, and adjacent tRNAs. These data contain 339 parsimony informative characters that we use to erect hypotheses of relationships for the Emydinae. Both maximum parsimony and maximum likelihood methods yield a monophyletic Emydinae in which all but three nodes are well resolved. Emys orbicularis, Emydoidea blandingii, and Clemmys marmorata form a monophyletic clade, as do the species of Terrapene. Clemmys muhlenbergii and Clemmys insculpta form a third monophyletic group that may be sister to all other emydines. Clemmys guttata is problematic and probably related to Terrapene. Based on this phylogeny, and previous molecular work on the group, we suggest the following taxonomic revisions: (1) Clemmys should be restricted to a single species, C. guttata. (2) Calemys should be resurrected for C. muhlenbergii and C. insculpta. (3) Emys should be expanded to include three species: E. orbicularis, E. blandingii, and E. marmorata. Furthermore, our analyses show that neither kinetic-shelled nor akinetic-shelled emydines form monophyletic groups. Therefore, shell kinesis was either independently gained in Emys and Terrapene or secondarily lost in E. marmorata and C. guttata. Parsimony, paleontological evidence, and the multiple origins of shell kinesis in related turtle lineages (especially geoemydines) support the independent origin of plastral kinesis.