Search

Effective biodiversity conservation planning starts with genetic characterization within and among focal populations, in order to understand the likely impact of threats for ensuring the long-term viability of a species. The Wonder Gecko, Teratoscincus keyserlingii, is one of nine members of the genus. This species is distributed in Iran, Afghanistan, and Pakistan, with a small isolated population in the United Arab Emirates (UAE), where it is classified nationally as Critically Endangered. Within its Arabian range, anthropogenic activity is directly linked to the species’ decline, with highly localised and severely fragmented populations. Here we describe the evolutionary history of Teratoscincus, by reconstructing its phylogenetic relationships and estimating its divergence times and ancestral biogeography. For conservation implications of T. keyserlingii we evaluate the genetic structure of the Arabian population using genomic data. This study supports the monophyly of most species and reveals considerable intraspecific variability in T. microlepis and T. keyserlingii, which necessitate broad systematic revisions. The UAE population of T. keyserlingii likely arrived from southern Iran during the Pleistocene and no internal structure was recovered within, implying a single population status. Regional conservation of T. keyserlingii requires improved land management and natural habitat restoration in the species’ present distribution, and expansion of current protected areas, or establishment of new areas with suitable habitat for the species, mostly in northern Abu Dhabi Emirate.

The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing., Correction in: Nature, 585, E3, 2020DOI: 10.1038/s41586-020-2661-6

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana(1,2)-is an iconic species that is endemic to New Zealand(2,3). A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes(2,4). Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.The approximately 5-Gb tuatara (Sphenodon punctatus) genome assembly provides a resource for analysing amniote evolution, and highlights the imperative for meaningful cultural engagement with Indigenous communities in genome-sequencing endeavours.

The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing., Correction in: Nature, 585, E3, 2020DOI: 10.1038/s41586-020-2661-6

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana(1,2)-is an iconic species that is endemic to New Zealand(2,3). A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes(2,4). Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.The approximately 5-Gb tuatara (Sphenodon punctatus) genome assembly provides a resource for analysing amniote evolution, and highlights the imperative for meaningful cultural engagement with Indigenous communities in genome-sequencing endeavours.

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, that can be especially powerful. We describe here the protocols commonly used for physically isolating mtDNA, for amplifying these by PCR or RCA, for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences to date with determining and comparing complete mtDNA sequences.

Mitochondrial (mt) DNA sequences are used extensively to reconstruct evolutionary relationships among recently diverged animals, and have constituted the most widely used markers for species- and generic-level relationships for the last decade or more. However, most studies to date have employed relatively small portions of the mt-genome. In contrast, complete mt-genomes primarily have been used to investigate deep divergences, including several studies of the amount of mt sequence necessary to recover ancient relationships. We sequenced and analyzed 24 complete mt-genomes from a group of salamander species exhibiting divergences typical of those in many species-level studies. We present the first comprehensive investigation of the amount of mt sequence data necessary to consistently recover the mt-genome tree at this level, using parsimony and Bayesian methods. Both methods of phylogenetic analysis revealed extremely similar results. A surprising number of well supported, yet conflicting, relationships were found in trees based on fragments less than ~;2000 nucleotides (nt), typical of the vast majority of the thousands of mt-based studies published to date. Large amounts of data (11,500+ nt) were necessary to consistently recover the whole mt-genome tree. Some relationships consistently were recovered with fragments of all sizes, but many nodes required the majority of the mt-genome to stabilize, particularly those associated with short internal branches. Although moderate amounts of data (2000-3000 nt) were adequate to recover mt-based relationships for which most nodes were congruent with the whole mt-genome tree, many thousands of nucleotides were necessary to resolve rapid bursts of evolution. Recent advances in genomics are making collection of large amounts of sequence data highly feasible, and our results provide the basis for comparative studies of other closely related groups to optimize mt sequence sampling and phylogenetic resolution at the "tips" of th

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, that can be especially powerful. We describe here the protocols commonly used for physically isolating mtDNA, for amplifying these by PCR or RCA, for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences to date with determining and comparing complete mtDNA sequences.

As part of an ongoing project to generate a mitochondrial database for terrestrial tortoises based on museum specimens, the complete mitochondrial genome sequences of 10 species and a ~;14 kb sequence from an eleventh species are reported. The sampling of the present study emphasizes Mediterranean tortoises (genus Testudo and their close relatives). Our new sequences are aligned, along with those of two testudinoid turtles from GenBank, Chrysemys picta and Mauremys reevesii, yielding an alignment of 14,858 positions, of which 3,238 are parsimony informative. We develop a phylogenetic taxonomy for Testudo and related species based on well-supported, diagnosable clades. Several well-supported nodes are recovered, including the monophyly of a restricted Testudo, T. kleinmanni + T. marginata (the Chersus clade), and the placement of the enigmatic African pancake tortoise (Malacochersus tornieri) within the predominantly Palearctic greater Testudo group (Testudona tax. nov.). Despite the large amount of sequence reported, there is low statistical support for some nodes within Testudona and so we do not propose names for those groups. A preliminary and conservative estimation of divergence times implies a late Miocene diversification for the testudonan clade (6-12 million years ago), matching their first appearance in the fossil record. The multi-continental distribution of testudonan turtles can be explained by the establishment of permanent connections between Europe, Africa, and Asia at this time. The arrival of testudonan turtles to Africa occurred after one or more initial tortoise invasions gave rise to the diverse (>25 species) 'Geochelone complex.' Two unusual genomic features are reported for the mtDNA of one tortoise, M. tornieri: (1) nad4 has a shift of reading frame that we suggest is resolved by translational frameshifting of the mRNA on the ribosome during protein synthesis and (2) there are two copies of the control region and trnF, with the latter havin

Phylogenetic relationships among species of the salamander family Salamandridae are investigated using nearly 3000 nucleotide bases of newly reported mitochondrial DNA sequence data from the mtDNA genic region spanning the genes tRNALeu-COI. This study uses nearly comprehensive species-level sampling to provide the first complete phylogeny for the Salamandridae. Deep phylogenetic relationships among the three most divergent lineages in the family Salamandrina terdigitata, a clade comprising the "True" salamanders, and a clade comprising all newts except S. terdigitata are difficult to resolve. However, most relationships within the latter two lineages are resolved with robust levels of branch support. The genera Euproctus and Triturus are statistically shown to be nonmonophyletic, instead each contains a diverse set of lineages positioned within the large newt clade. The genus Paramesotriton is also resolve as a nonmonophyletic group, with the newly described species P. laoensis constituting a divergent lineage placed in a sister position to clade containing all Pachytriton species and all remaining Paramesotriton species. Sequence divergences between P. laoensis and other Paramesotriton species are as great as those comparing P. laoensis and species of the genera Cynops and Pachytriton. Analyses of lineage diversification across the Salamandridae indicate that, despite its exceptional diversity, lineage accumulation appears to have been constant across time, indicating that it does not represent a true species radiation.

The analysis of interactions between lineages at varying levels of genetic divergence can provide insights into the process of speciation through the accumulation of incompatible mutations. Ring species, and especially the Ensatina eschscholtzii system exemplify this approach. The plethodontid salamanders Ensatina eschscholtzii xanthoptica and Ensatina eschscholtzii platensis hybridize in the Central Sierran foothills of California. We compared the genetic structure across two transects (southern and northern Calaveras Co.), one of which was re-sampled over 20 years, and examined diagnostic molecular markers (eight allozyme loci and mitochondrial DNA) and a diagnostic quantitative trait (color pattern). Key results across all studies were: i) cline centers for all markers were coincident and the zones were narrow, with width estimates of 730m to 2000m; ii) cline centers at the northern Calaveras transect were coincident between 1981 and 2001, demonstrating repeatability over 5 generations; iii) there are very few if any putative F1's, but a relatively high number of backcrossed individuals (57-86 percent) in the central portion of transects; iv) we found substantial linkage disequilibrium in all three studies and strong heterozygote deficit both in northern Calaveras, in 2001, and southern Calaveras. Both linkage disequilibrium and heterozygote deficit show maximum values near the center of the zones (R and Fis, approx. equal to 0.5). Using estimates of cline width and dispersal, we infer strong selection against hybrids (s* approx. equal to 46-75 percent). This is sufficient to promote accumulation of differences at loci that are neutral or under divergent selection, but would still allow for introgression of a daptive alleles. The evidence for strong, but incomplete isolation across this centrally located contact is consistent with theory suggesting a gradual increase in postzygotic incompatibility between allopatric populations subject to divergent selection and

Mitochondrial (mt) DNA sequences are used extensively to reconstruct evolutionary relationships among recently diverged animals, and have constituted the most widely used markers for species- and generic-level relationships for the last decade or more. However, most studies to date have employed relatively small portions of the mt-genome. In contrast, complete mt-genomes primarily have been used to investigate deep divergences, including several studies of the amount of mt sequence necessary to recover ancient relationships. We sequenced and analyzed 24 complete mt-genomes from a group of salamander species exhibiting divergences typical of those in many species-level studies. We present the first comprehensive investigation of the amount of mt sequence data necessary to consistently recover the mt-genome tree at this level, using parsimony and Bayesian methods. Both methods of phylogenetic analysis revealed extremely similar results. A surprising number of well supported, yet conflicting, relationships were found in trees based on fragments less than ~;2000 nucleotides (nt), typical of the vast majority of the thousands of mt-based studies published to date. Large amounts of data (11,500+ nt) were necessary to consistently recover the whole mt-genome tree. Some relationships consistently were recovered with fragments of all sizes, but many nodes required the majority of the mt-genome to stabilize, particularly those associated with short internal branches. Although moderate amounts of data (2000-3000 nt) were adequate to recover mt-based relationships for which most nodes were congruent with the whole mt-genome tree, many thousands of nucleotides were necessary to resolve rapid bursts of evolution. Recent advances in genomics are making collection of large amounts of sequence data highly feasible, and our results provide the basis for comparative studies of other closely related groups to optimize mt sequence sampling and phylogenetic resolution at the "tips" of th

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, that can be especially powerful. We describe here the protocols commonly used for physically isolating mtDNA, for amplifying these by PCR or RCA, for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences to date with determining and comparing complete mtDNA sequences.

As part of an ongoing project to generate a mitochondrial database for terrestrial tortoises based on museum specimens, the complete mitochondrial genome sequences of 10 species and a ~;14 kb sequence from an eleventh species are reported. The sampling of the present study emphasizes Mediterranean tortoises (genus Testudo and their close relatives). Our new sequences are aligned, along with those of two testudinoid turtles from GenBank, Chrysemys picta and Mauremys reevesii, yielding an alignment of 14,858 positions, of which 3,238 are parsimony informative. We develop a phylogenetic taxonomy for Testudo and related species based on well-supported, diagnosable clades. Several well-supported nodes are recovered, including the monophyly of a restricted Testudo, T. kleinmanni + T. marginata (the Chersus clade), and the placement of the enigmatic African pancake tortoise (Malacochersus tornieri) within the predominantly Palearctic greater Testudo group (Testudona tax. nov.). Despite the large amount of sequence reported, there is low statistical support for some nodes within Testudona and so we do not propose names for those groups. A preliminary and conservative estimation of divergence times implies a late Miocene diversification for the testudonan clade (6-12 million years ago), matching their first appearance in the fossil record. The multi-continental distribution of testudonan turtles can be explained by the establishment of permanent connections between Europe, Africa, and Asia at this time. The arrival of testudonan turtles to Africa occurred after one or more initial tortoise invasions gave rise to the diverse (>25 species) 'Geochelone complex.' Two unusual genomic features are reported for the mtDNA of one tortoise, M. tornieri: (1) nad4 has a shift of reading frame that we suggest is resolved by translational frameshifting of the mRNA on the ribosome during protein synthesis and (2) there are two copies of the control region and trnF, with the latter havin

The analysis of interactions between lineages at varying levels of genetic divergence can provide insights into the process of speciation through the accumulation of incompatible mutations. Ring species, and especially the Ensatina eschscholtzii system exemplify this approach. The plethodontid salamanders Ensatina eschscholtzii xanthoptica and Ensatina eschscholtzii platensis hybridize in the Central Sierran foothills of California. We compared the genetic structure across two transects (southern and northern Calaveras Co.), one of which was re-sampled over 20 years, and examined diagnostic molecular markers (eight allozyme loci and mitochondrial DNA) and a diagnostic quantitative trait (color pattern). Key results across all studies were: i) cline centers for all markers were coincident and the zones were narrow, with width estimates of 730m to 2000m; ii) cline centers at the northern Calaveras transect were coincident between 1981 and 2001, demonstrating repeatability over 5 generations; iii) there are very few if any putative F1's, but a relatively high number of backcrossed individuals (57-86 percent) in the central portion of transects; iv) we found substantial linkage disequilibrium in all three studies and strong heterozygote deficit both in northern Calaveras, in 2001, and southern Calaveras. Both linkage disequilibrium and heterozygote deficit show maximum values near the center of the zones (R and Fis, approx. equal to 0.5). Using estimates of cline width and dispersal, we infer strong selection against hybrids (s* approx. equal to 46-75 percent). This is sufficient to promote accumulation of differences at loci that are neutral or under divergent selection, but would still allow for introgression of a daptive alleles. The evidence for strong, but incomplete isolation across this centrally located contact is consistent with theory suggesting a gradual increase in postzygotic incompatibility between allopatric populations subject to divergent selection and

Phylogenetic relationships among species of the salamander family Salamandridae are investigated using nearly 3000 nucleotide bases of newly reported mitochondrial DNA sequence data from the mtDNA genic region spanning the genes tRNALeu-COI. This study uses nearly comprehensive species-level sampling to provide the first complete phylogeny for the Salamandridae. Deep phylogenetic relationships among the three most divergent lineages in the family Salamandrina terdigitata, a clade comprising the "True" salamanders, and a clade comprising all newts except S. terdigitata are difficult to resolve. However, most relationships within the latter two lineages are resolved with robust levels of branch support. The genera Euproctus and Triturus are statistically shown to be nonmonophyletic, instead each contains a diverse set of lineages positioned within the large newt clade. The genus Paramesotriton is also resolve as a nonmonophyletic group, with the newly described species P. laoensis constituting a divergent lineage placed in a sister position to clade containing all Pachytriton species and all remaining Paramesotriton species. Sequence divergences between P. laoensis and other Paramesotriton species are as great as those comparing P. laoensis and species of the genera Cynops and Pachytriton. Analyses of lineage diversification across the Salamandridae indicate that, despite its exceptional diversity, lineage accumulation appears to have been constant across time, indicating that it does not represent a true species radiation.

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, that can be especially powerful. We describe here the protocols commonly used for physically isolating mtDNA, for amplifying these by PCR or RCA, for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences to date with determining and comparing complete mtDNA sequences.

The evolutionary history of the largest salamander family (Plethodontidae) is characterized by extreme morphological homoplasy. Analysis of the mechanisms generating such homoplasy requires an independent, molecular phylogeny. To this end, we sequenced 24 complete mitochondrial genomes (22 plethodontids and two outgroup taxa), added data for three species from GenBank, and performed partitioned and unpartitioned Bayesian, ML, and MP phylogenetic analyses. We explored four dataset partitioning strategies to account for evolutionary process heterogeneity among genes and codon positions, all of which yielded increased model likelihoods and decreased numbers of supported nodes in the topologies (PP > 0.95) relative to the unpartitioned analysis. Our phylogenetic analyses yielded congruent trees that contrast with the traditional morphology-based taxonomy; the monophyly of three out of four major groups is rejected. Reanalysis of current hypotheses in light of these new evolutionary relationships suggests that 1) a larval life history stage re-evolved from a direct-developing ancestor multiple times, 2) there is no phylogenetic support for the "Out of Appalachia" hypothesis of plethodontid origins, and 3) novel scenarios must be reconstructed for the convergent evolution of projectile tongues, reduction in toe number, and specialization for defensive tail loss. Some of these novel scenarios imply morphological transformation series that proceed in the opposite direction than was previously thought. In addition, they suggest surprising evolutionary lability in traits previously interpreted to be conservative.

The evolutionary history of the largest salamander family (Plethodontidae) is characterized by extreme morphological homoplasy. Analysis of the mechanisms generating such homoplasy requires an independent, molecular phylogeny. To this end, we sequenced 24 complete mitochondrial genomes (22 plethodontids and two outgroup taxa), added data for three species from GenBank, and performed partitioned and unpartitioned Bayesian, ML, and MP phylogenetic analyses. We explored four dataset partitioning strategies to account for evolutionary process heterogeneity among genes and codon positions, all of which yielded increased model likelihoods and decreased numbers of supported nodes in the topologies (PP > 0.95) relative to the unpartitioned analysis. Our phylogenetic analyses yielded congruent trees that contrast with the traditional morphology-based taxonomy; the monophyly of three out of four major groups is rejected. Reanalysis of current hypotheses in light of these new evolutionary relationships suggests that 1) a larval life history stage re-evolved from a direct-developing ancestor multiple times, 2) there is no phylogenetic support for the "Out of Appalachia" hypothesis of plethodontid origins, and 3) novel scenarios must be reconstructed for the convergent evolution of projectile tongues, reduction in toe number, and specialization for defensive tail loss. Some of these novel scenarios imply morphological transformation series that proceed in the opposite direction than was previously thought. In addition, they suggest surprising evolutionary lability in traits previously interpreted to be conservative.