Search

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict., (© The Author(s) 2023. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

In a genomic study by Chan and colleagues, pit-vipers of the Trimeresuruserythrurus - purpureomaculatus complex from the Ayeyarwady and Yangon regions in Myanmar were demonstrated to be a distinct species based on robust population genetic and species delimitation analyses. Here, we provide morphological characterizations and a formal description of those populations as a new species. The new species, Trimeresurusayeyarwadyensis sp. nov. , is most closely related to T.erythrurus and T.purpureomaculatus and shares morphological characteristics with both of those species. Some specimens of T.ayeyarwadyensis sp. nov. have green dorsal coloration and no distinct dorsal blotches (a trait shared with T.erythrurus but not T.purpureomaculatus ), while others have dark dorsal blotches (a trait shared with T.purpureomaculatus but not T.erythrurus ). The distinct evolutionary trajectory of the new species, coupled with the lack of obvious morphological differentiation, represents a classic example of the cryptic nature of species commonly found in the Trimeresurus group of Asian pit-vipers and underscores the need for data-rich analyses to verify species' boundaries more broadly within this genus., Competing Interests: The authors have declared that no competing interests exist., (Kin Onn Chan, Shahrul Anuar, Ananthanarayanan Sankar, Ingg Thong Law, Ing Sind Law, Rasu Shivaram, Ching Christian, Daniel G. Mulcahy, Anita Malhotra.)

A new species of Dipsas Laurenti, 1768, from Central Panama is described based on molecular analyses, hemipenial morphology, and external characters. This is the sixth species of Dipsas to be described for the country; the snake has been suspected to exist since 1977 and has not been thoroughly studied until now. Additionally, morphological comparations including scale counts are done with other species within the genus, and the current geographic distribution of Dipsastemporalis (Werner, 1909), the sister species, is updated. Finally, a key to the species of Dipsas currently known from Middle America is presented., (Julie M. Ray, Paola Sánchez-Martínez, Abel Batista, Daniel G. Mulcahy, Coleman M. Sheehy III, Eric N. Smith, R. Alexander Pyron, Alejandro Arteaga.)

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, Anolis opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting (ILS). Here, we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene tree, species tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene tree recovers the same relationships previously inferred for this group, which is strikingly different from the species tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model (MSCM), so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group; however, coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent, adaptive radiations. [Adaptive radiation; hybridization; introgression; multispecies network coalescent; posterior predictive simulation.]., (© The Authors 2021. Published by Oxford University Press on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Natural history collections are essential to a wide variety of studies in biology because they maintain large collections of specimens and associated data, including genetic material (e.g., tissues) for DNA sequence data, yet they are currently under-funded and collection staff have high workloads. With the advent of aggregate databases and advances in sequencing technologies, there is an increased demand on collection staff for access to tissue samples and associated data. Scientists are rapidly developing large DNA barcode libraries, DNA sequences of specific genes for species across the tree of life, in order to document and conserve biodiversity. In doing so, mistakes are made. For instance, inconsistent taxonomic information is commonly taken from different lending institutions and deposited in data repositories, such as the Barcode of Life Database (BOLD) and GenBank, despite explicit disclaimers regarding the need for taxonomic verification by the lending institutions. Such errors can have profound effects on subsequent research based on these mis-labelled sequences in data repositories. Here, we present the production of a large DNA barcode library of reptiles from the National Museum of Natural History tissue holdings. The library contains 2,758 sequences (2,205 COI and 553 16S) from 2260 specimens (four crocodilians, 37 turtles, and 2,219 lizards, including snakes), representing 583 named species, from 52 countries. In generating this library, we noticed several common mistakes made by scientists depositing DNA barcode data in public repositories (e.g., BOLD and GenBank). Our goal is to raise awareness of these concerns and offer advice to avoid such mistakes in the future to maintain accurate DNA barcode libraries to properly document Earth's biodiversity., Competing Interests: The authors have declared that no competing interests exist.

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

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 North American nightsnakes in the genus Hypsiglena is composed of nine named and at least two unnamed species. Here, we provide the first mt-genome of H. affinis , an additional mt-genome for H. sp. nov. 1, and four additional mt-genomes from the widespread H. jani . These mtDNA genomes were sequenced using both Illumina and Ion Torrent sequencing technologies. The resulting genomes contained the expected 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions typical of colubroid snakes. Two of the H . jani samples had partial tRNA Ile genes upstream of CR2 which has not been previously documented in colubroid snakes. A maximum likelihood gene-tree based on these data combined with previously published sequence data recovers a well-supported phylogeny and is in concordance with previous estimates of evolutionary relationships in this group., Competing Interests: No potential conflict of interest was reported by the author(s)., (© 2020 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

Fringe-toed lizards (Uma) are among the most specialized lizards in North America, adapted to insular windblown sand habitats in the hyper-arid southwestern deserts, with allopatric distributions, subtle morphological variation, and an unstable taxonomic history. We analyzed a morphological dataset of 40 characters for 65 specimens and a molecular dataset of 2,286 bases from three mitochondrial loci for 92 individuals and interpreted these data alongside published analyses of multi-locus genetic data with the goal of revising the taxonomy of the Uma notata (Baird 1858) species complex. We confirmed that fringe-toed lizards from the Mohawk Dunes in southwestern Arizona (U. sp.) constitute a cryptic species sister to the rest of the complex that can be diagnosed with DNA barcoding and geography, so we describe and name this species Uma thurmanae sp. nov. We also confirmed the evolutionary distinctiveness of U. inornata (Cope 1895), an endangered species endemic to Coachella Valley in southern California. We designate a lectotype for the taxon U. "rufopunctata", but we put its name in quotation marks to reflect its uncertain taxonomic status with respect to its neighboring species U. cowlesi and U. notata.

We report the first occurrence of the Asian Rock Gecko genus Cnemaspis Strauch, 1887 from mainland Myanmar based on a series of specimens recently collected from the Tanintharyi Region. These records come from several localities in the Tanintharyi Region and fill a significant sampling gap for the genus. Molecular phylogenetic analyses using the mitochondrial gene ND2 identify two distinct clades, the first includes Cnemaspis siamensis of the Southeast Asian Cnemaspis group and the second includes two new species belonging to the South Asian Cnemaspis kandiana group. These two species are morphologically distinct and are distinguished from all other members in the C. kandiana group by a combination of character states. The first species, Cnemaspis tanintharyi sp. nov. occurs on the mainland in southern Tanintharyi and is distinguished from all congeners by the possession of keeled pectoral scales; smooth ventral scales and abdominal scales; 2-4 precloacal pores; 4-5 femoral pores on each leg; smooth subcaudal scales with the median row enlarged; coloration of the gular region beige, dark gray-brown with dark blue hueing towards throat; 15-18 subdigital lamellae on the 4th toe; and 21-23 ventral scales at midbody. The other, Cnemaspis thayawthadangyi sp. nov. is known only from the island group of Thayawthadangyi, in the Myeik Archipelago and is distinguished from all congeners by the possession of keeled pectoral scales; keeled ventral scales and abdominal scales; three precloacal pores; four femoral pores on each leg; smooth subcaudal scales, scales on median row enlarged; coloration of the gular region silver with dark-gray irregularly shaped streaks; 16-18 subdigital lamellae on the 4th toe; and 18-20 ventral scales at midbody. In addition, we address the taxonomic status of populations referred to as C. kandiana (a species now restricted to Sri Lanka) found in peninsular Thailand and provisionally assign them as C. cf. tanintharyi sp. nov. until further analysis can be conducted. Finally, we briefly discuss the biogeography of the South Asian clade of Cnemaspis.

An integrative taxonomic analysis of the Ptychozoon lionotum group across its range in Indochina and Sundaland recovers P. lionotum sensu lato Annandale, 1905 as paraphyletic with respect to P. popaense Grismer, Wood, Thura, Grismer, Brown, Stuart, 2018a and composed of four allopatric, genetically divergent, ND2 mitochondrial lineages. Multivariate and univariate analyses of continuous and discrete morphological and color pattern characters statistically and discretely diagnose each lineage from one another and together, with maximum likelihood and Bayesian inference analyses, provide the foundation for the recognition of each lineage as a new species-hypotheses corroborated with a Generalized Mixed Yule Coalescent species delimitation analysis. Ptychozoon cicakterbang sp. nov. ranges throughout Peninsular Malaysia to Pulau Natuna Besar, Indonesia; P. kabkaebin sp. nov. is endemic to northern and central Laos; and P. tokehos sp. nov. ranges from southern Thailand south of the Isthmus of Kra northward to Chiang Mai, fringing the Chao Phraya Basin and ranging southward through Cambodia to southern Vietnam. Ptychozoon lionotum sensu stricto ranges from northwestern Laos through southern Myanmar to eastern India. The phylogeographic structure within each species varies considerably with P. lionotum s.s. showing no genetic divergence across its 1,100 km range compared to P. cicakterbang sp. nov. showing upwards of 8.2% sequence divergence between syntopic individuals. Significant phylogeographic structure exists within P. tokehos sp. nov. and increased sampling throughout Thailand may require additional taxonomic changes within this species.

Despite threats of species extinctions, taxonomic crises, and technological advances in genomics and natural history database informatics, we are still distant from cataloguing all of the species of life on earth. Amphibians and reptiles are no exceptions; in fact new species are described nearly every day and many species face possible extinction. The number of described species continues to climb as new areas of the world are explored and as species complexes are examined more thoroughly. The use of DNA barcoding provides a mechanism for rapidly estimating the number of species at a given site and has the potential to record all of the species of life on Earth. Though DNA barcoding has its caveats, it can be useful to estimate the number of species in a more systematic and efficient manner, to be followed in combination with more traditional, morphology-based identifications and species descriptions. Herein, we report the results of a voucher-based herpetological expedition to the Tanintharyi (Tenasserim) Region of Myanmar, enhanced with DNA barcode data. Our main surveys took place in the currently proposed Tanintharyi National Park. We combine our results with photographs and observational data from the Chaung-nauk-pyan forest reserve. Additionally, we provide the first checklist of amphibians and reptiles of the region, with species based on the literature and museum. Amphibians, anurans in particular, are one of the most poorly known groups of vertebrates in terms of taxonomy and the number of known species, particularly in Southeast Asia. Our rapid-assessment program combined with DNA barcoding and use of Barcode Index Numbers (BINs) of voucher specimens reveals the depth of taxonomic diversity in the southern Tanintharyi herpetofauna even though only a third of the potential amphibians and reptiles were seen. A total of 51 putative species (one caecilian, 25 frogs, 13 lizards, 10 snakes, and two turtles) were detected, several of which represent potentially undescribed species. Several of these species were detected by DNA barcode data alone. Furthermore, five species were recorded for the first time in Myanmar, two amphibians (Ichthyophis cf. kohtaoensis and Chalcorana eschatia ) and three snakes ( Ahaetulla mycterizans , Boiga dendrophila , and Boiga drapiezii ).

Genome sequencing initiatives like the Arthropod i5k project and other biodiversity genomics research rely on access to high quality DNA and/or tissue. Global collection initiatives such as the Smithsonian Global Genome Initiative (GGI) and its partner network, the Global Genome Biodiversity Network (GGBN) aim to provide access to these resources at high-quality standards. Here, we review progress toward providing genomic resources (tissues, DNA, genomes) for terrestrial arthropods, a megadiverse animal group, and compare progress in genome sequencing to all other animals., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

The taxonomic identity of the Trimeresurus (Popeia) popeiorum complex from the Isthmus of Kra and to the north was investigated. Several studies over the last decade have produced several specimens and associated mtDNA sequence data for a variety of individuals of the T. popeiorum and "T. sabahi" complexes. Here, we combine four mitochondrial genes (12S, 16S, ND4, and CytB) from all available specimens in GenBank with the addition of five new specimens collected from the mainland, Tanintharyi Region of Myanmar. Maximum Likelihood and Bayesian analyses identified that T. popeiorum sensu lato is paraphyletic with two geographically distinct clades: a northern clade representing populations from northern Myanmar, Laos and northern Thailand and a southern clade representing samples from the Tanintharyi Region and adjacent west Thailand. While the two clades have considerable genetic distance, they appear to be morphologically identical, leading to the hypothesis that the southern clade represents a cryptic, undescribed species. Because they appear to be cryptic species and the limitation of only five specimens from the southern lineage, this does not permit us to formally describe the new species. In accordance to past molecular studies, we uncovered paraphyly and lack of genetic support for the validity of taxa within the T. sabahi complex. However, we suggest recognizing these populations as subspecies within T. sabahi.

Development projects in west Central Africa are proceeding at an unprecedented rate, often with little concern for their effects on biodiversity. In an attempt to better understand potential impacts of a road development project on the anuran amphibian community, we conducted a biodiversity assessment employing multiple methodologies (visual encounter transects, auditory surveys, leaf litter plots and pitfall traps) to inventory species prior to construction of a new road within the buffer zone of Moukalaba-Doudou National Park, Gabon. Because of difficulties in morphological identification and taxonomic uncertainty of amphibian species observed in the area, we integrated a DNA barcoding analysis into the project to improve the overall quality and accuracy of the species inventory. Based on morphology alone, 48 species were recognized in the field and voucher specimens of each were collected. We used tissue samples from specimens collected at our field site, material available from amphibians collected in other parts of Gabon and the Republic of Congo to initiate a DNA barcode library for west Central African amphibians. We then compared our sequences with material in GenBank for the genera recorded at the study site to assist in identifications. The resulting COI and 16S barcode library allowed us to update the number of species documented at the study site to 28, thereby providing a more accurate assessment of diversity and distributions. We caution that because sequence data maintained in GenBank are often poorly curated by the original submitters and cannot be amended by third-parties, these data have limited utility for identification purposes. Nevertheless, the use of DNA barcoding is likely to benefit biodiversity inventories and long-term monitoring, particularly for taxa that can be difficult to identify based on morphology alone; likewise, inventory and monitoring programs can contribute invaluable data to the DNA barcode library and the taxonomy of complex groups. Our methods provide an example of how non-taxonomists and parataxonomists working in understudied parts of the world with limited geographic sampling and comparative morphological material can use DNA barcoding and publicly available sequence data (GenBank) to rapidly identify the number of species and assign tentative names to aid in urgent conservation management actions and contribute to taxonomic resolution.

The Global Genome Initiative has sequenced and released 1961 DNA barcodes for genetic samples obtained as part of the Global Genome Initiative for Gardens Program. The dataset includes barcodes for 29 plant families and 309 genera that did not have sequences flagged as barcodes in GenBank and sequences from officially recognized barcoding genetic markers meet the data standard of the Consortium for the Barcode of Life. The genetic samples were deposited in the Smithsonian Institution's National Museum of Natural History Biorepository and their records were made public through the Global Genome Biodiversity Network's portal. The DNA barcodes are now available on GenBank.

Myanmar's recent transition from military rule towards a more democratic government has largely ended decades of political and economic isolation. Although Myanmar remains heavily forested, increased development in recent years has been accompanied by exceptionally high rates of forest loss. In this study, we document the rapid progression of deforestation in and around the proposed Lenya National Park, which includes some of the largest remaining areas of lowland evergreen rainforest in mainland Southeast Asia. The globally unique forests in this area are rich in biodiversity and remain a critical stronghold for many threatened and endangered species, including large charismatic fauna such as tiger and Asian elephant. We also conducted a rapid assessment survey of the herpetofauna of the proposed national park, which resulted in the discovery of two new species of bent-toed geckos, genus Cyrtodactylus. We describe these new species, C. lenya sp. nov. and C. payarhtanensis sp. nov., which were found in association with karst (i.e., limestone) rock formations within mature lowland wet evergreen forest. The two species were discovered less than 35 km apart and are each known from only a single locality. Because of the isolated nature of the karst formations in the proposed Lenya National Park, these geckos likely have geographical ranges restricted to the proposed protected area and are threatened by approaching deforestation. Although lowland evergreen rainforest has vanished from most of continental Southeast Asia, Myanmar can still take decisive action to preserve one of the most biodiverse places on Earth.

Recent fieldwork in southern Tanintharyi revealed the presence of a small Green Crested Lizard in the wet evergreen forest. We generated mtDNA sequence data (ND2) that demonstrates that this population's nearest relative is Bronchocela rayaensis Grismer et al., 2015 of Pulau Langkawi, northwestern Peninsular Malaysia and Phuket Island. Morphologically the Burmese Bronchocela shares many features with Bronchocela rayaensis , which potentially would make this recently described Thai-Malay species a synonym of Bronchocela burmana Blanford, 1878; however, we interpret the genetic and morphological differences to reflect evolutionary divergence and recommend the recognition of both species.

Advances in biodiversity genomic sequencing will increasingly depend on the availability of DNA samples-and their quantifiable metadata-preserved in large institutional biorepositories that are discoverable to the scientific community. Improvements in sequencing technology constantly provide longer reads, such that longer fragment length, higher molecular weight, and overall "genome-quality" DNA (gDNA) will be desirable. Ideally, biorepositories should publish numerical scale measurements of DNA quality useful to the user community. However, the most widely used technique to evaluate DNA quality, the classic agarose gel, has yet to be quantified. Here we propose a simple and economical method using open source image analysis software to make gDNA gel images quantifiable, and propose percentage of gDNA "greater than X kb" as a standard of comparison, where X is a band from any widely used DNA ladder with desirably large band sizes. We employ two metadata standards ("DNA Threshold" and "Percent above Threshold") introduced as part of the Global Genome Biodiversity Network (GGBN) Darwin Core extension. We illustrate the method using the traditionally used Hin dIII ladder and the 9,416 base-pair (bp) band as a standard. We also present data, for two taxa, a vertebrate (fish) and an invertebrate (crab), on how gDNA quality varies with seven tissue preservation methods, time since death, preservation method (i.e. buffers vs. cold temperatures), and storage temperature of various buffers over time. Our results suggest that putting tissue into a buffer prior to freezing may be better than directly into ultra-cold conditions., Competing Interests: Sean G. Brady is an Academic Editor for PeerJ.

During climate change, species are often assumed to shift their geographic distributions (geographic ranges) in order to track environmental conditions - niches - to which they are adapted. Recent work, however, suggests that the niches do not always remain conserved during climate change but shift instead, allowing populations to persist in place or expand into new areas. We assessed the extent of range and niche shifts in response to the warming climate after the Last Glacial Maximum (LGM) in the desert horned lizard ( Phrynosoma platyrhinos ), a species occupying the western deserts of North America. We used a phylogeographic approach with mitochondrial DNA sequences to approximate the species range during the LGM by identifying populations that exhibit a genetic signal of population stability versus those that exhibit a signal of a recent (likely post-LGM) geographic expansion. We then compared the climatic niche that the species occupies today with the niche it occupied during the LGM using two models of simulated LGM climate. The genetic analyses indicated that P. platyrhinos persisted within the southern Mojave and Sonoran deserts throughout the latest glacial period and expanded from these deserts northwards, into the western and eastern Great Basin, after the LGM. The climatic niche comparisons revealed that P. platyrhinos expanded its climatic niche after the LGM towards novel, warmer and drier climates that allowed it to persist within the southern deserts. Simultaneously, the species shifted its climatic niche towards greater temperature and precipitation fluctuations after the LGM. We concluded that climatic changes at the end of the LGM promoted both range and niche shifts in this lizard. The mechanism that allowed the species to shift its niche remains unknown, but phenotypic plasticity likely contributes to the species ability to adjust to climate change.

Phylogenetic analysis of extinction threat is an emerging tool in the field of conservation. However, there are problems with the methods and data as commonly used. Phylogenetic sampling usually extends to the level of family or genus, but International Union for Conservation of Nature (IUCN) rankings are available only for individual species, and, although different species within a taxonomic group may have the same IUCN rank, the species may have been ranked as such for different reasons. Therefore, IUCN rank may not reflect evolutionary history and thus may not be appropriate for use in a phylogenetic context. To be used appropriately, threat-risk data should reflect the cause of extinction threat rather than the IUCN threat ranking. In a case study of the toad genus Incilius, with phylogenetic sampling at the species level (so that the resolution of the phylogeny matches character data from the IUCN Red List), we analyzed causes of decline and IUCN threat rankings by calculating metrics of phylogenetic signal (such as Fritz and Purvis' D). We also analyzed the extent to which cause of decline and threat ranking overlap by calculating phylogenetic correlation between these 2 types of character data. Incilius species varied greatly in both threat ranking and cause of decline; this variability would be lost at a coarser taxonomic resolution. We found far more phylogenetic signal, likely correlated with evolutionary history, for causes of decline than for IUCN threat ranking. Individual causes of decline and IUCN threat rankings were largely uncorrelated on the phylogeny. Our results demonstrate the importance of character selection and taxonomic resolution when extinction threat is analyzed in a phylogenetic context., (© 2015 Society for Conservation Biology.)

Squamate reptiles (lizards and snakes) are a pivotal group whose relationships have become increasingly controversial. Squamates include >9000 species, making them the second largest group of terrestrial vertebrates. They are important medicinally and as model systems for ecological and evolutionary research. However, studies of squamate biology are hindered by uncertainty over their relationships, and some consider squamate phylogeny unresolved, given recent conflicts between molecular and morphological results. To resolve these conflicts, we expand existing morphological and molecular datasets for squamates (691 morphological characters and 46 genes, for 161 living and 49 fossil taxa, including a new set of 81 morphological characters and adding two genes from published studies) and perform integrated analyses. Our results resolve higher-level relationships as indicated by molecular analyses, and reveal hidden morphological support for the molecular hypothesis (but not vice-versa). Furthermore, we find that integrating molecular, morphological, and paleontological data leads to surprising placements for two major fossil clades (Mosasauria and Polyglyphanodontia). These results further demonstrate the importance of combining fossil and molecular information, and the potential problems of estimating the placement of fossil taxa from morphological data alone. Thus, our results caution against estimating fossil relationships without considering relevant molecular data, and against placing fossils into molecular trees (e.g. for dating analyses) without considering the possible impact of molecular data on their placement.

Vertebrates are currently going extinct at an alarming rate, largely because of habitat loss, global warming, infectious diseases, and human introductions. Island ecosystems are particularly vulnerable to invasive species and other ecological disturbances. Properly documenting historic and current species distributions is critical for quantifying extinction events. Museum specimens, field notes, and other archived materials from historical expeditions are essential for documenting recent changes in biodiversity. The Islas Revillagigedo are a remote group of four islands, 700-1100 km off the western coast of mainland México. The islands are home to many endemic plants and animals recognized at the specific- and subspecific-levels, several of which are currently threatened or have already gone extinct. Here, we recount the initial discovery of an endemic snake Hypsiglena ochrorhyncha unaocularus Tanner on Isla Clarión, the later dismissal of its existence, its absence from decades of field surveys, our recent rediscovery, and recognition of it as a distinct species. We collected two novel complete mitochondrial (mt) DNA genomes and up to 2800 base-pairs of mtDNA from several other individuals, aligned these with previously published mt-genome data from samples throughout the range of Hypsiglena, and conducted phylogenetic analyses to infer the biogeographic origin and taxonomic status of this population. We found the Isla Clarión population to be most closely related to populations in the Sonora-Sinaloa state border area of mainland México and Isla Santa Catalina, in the Gulf of California. Based on genetics, morphology, and geographic distributions, we also recognize these two other lineages as distinct species. Our study shows the importance of museum specimens, field notes, and careful surveys to accurately document biodiversity and brings these island endemics (Clarión and Santa Catalina nightsnakes) and mainland population near the Sonora-Sinaloa state border to the attention of conservation biologists currently monitoring biodiversity in these fragile subtropical ecosystems.

The salamander family Hynobiidae contains over 50 species and has been the subject of a number of molecular phylogenetic investigations aimed at reconstructing branches across the entire family. In general, studies using the greatest amount of sequence data have used reduced taxon sampling, while the study with the greatest taxon sampling has used a limited sequence data set. Here, we provide insights into the phylogenetic history of the Hynobiidae using both dense taxon sampling and a large mitochondrial DNA sequence data set. We report exclusive new mitochondrial DNA data of 2566 aligned bases (with 151 excluded sites, of included sites 1157 are variable with 957 parsimony informative). This is sampled from two genic regions encoding a 12S-16S region (the 3' end of 12S rRNA, tRNA(VAI), and the 5' end of 16S rRNA), and a ND2-COI region (ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), the origin for light strand replication--O(L), tRNA(Cys), tRNAT(Tyr), and the 5' end of COI). Analyses using parsimony, Bayesian, and maximum likelihood optimality criteria produce similar phylogenetic trees, with discordant branches generally receiving low levels of branch support. Monophyly of the Hynobiidae is strongly supported across all analyses, as is the sister relationship and deep divergence between the genus Onychodactylus with all remaining hynobiids. Within this latter grouping our phylogenetic results identify six clades that are relatively divergent from one another, but for which there is minimal support for their phylogenetic placement. This includes the genus Batrachuperus, the genus Hynobius, the genus Pachyhynobius, the genus Salamandrella, a clade containing the genera Ranodon and Paradactylodon, and a clade containing the genera Liua and Pseudohynobius. This latter clade receives low bootstrap support in the parsimony analysis, but is consistent across all three analytical methods. Our results also clarify a number of well-supported relationships within the larger Batrachuperus and Hynobius clades. While the relationships identified in this study do much to clarify the phylogenetic history of the Hynobiidae, the poor resolution among major hynobiid clades, and the contrast of mtDNA-derived relationships with recent phylogenetic results from a small number of nuclear genes, highlights the need for continued phylogenetic study with larger numbers of nuclear loci.

Squamate reptiles (lizards and snakes) are one of the most diverse groups of terrestrial vertebrates. Recent molecular analyses have suggested a very different squamate phylogeny relative to morphological hypotheses, but many aspects remain uncertain from molecular data. Here, we analyse higher-level squamate phylogeny with a molecular dataset of unprecedented size, including 161 squamate species for up to 44 nuclear genes each (33 717 base pairs), using both concatenated and species-tree methods for the first time. Our results strongly resolve most squamate relationships and reveal some surprising results. In contrast to most other recent studies, we find that dibamids and gekkotans are together the sister group to all other squamates. Remarkably, we find that the distinctive scolecophidians (blind snakes) are paraphyletic with respect to other snakes, suggesting that snakes were primitively burrowers and subsequently re-invaded surface habitats. Finally, we find that some clades remain poorly supported, despite our extensive data. Our analyses show that weakly supported clades are associated with relatively short branches for which individual genes often show conflicting relationships. These latter results have important implications for all studies that attempt to resolve phylogenies with large-scale phylogenomic datasets.

Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1,000 bp per gene; 19,020bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1 kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0 Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0 Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82 Myr older (mean=73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50 Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates., (Copyright © 2012 Elsevier Inc. All rights reserved.)

We investigated the evolutionary relationships among populations of two threatened Red Data Book fossorial skinks, Scelotes gronovii and Scelotes kasneri, along the Western Cape Coast of South Africa. The genus Scelotes shows considerable variation in limb and digit reduction. We sampled four localities purported to contain S. gronovii and seven of S. kasneri, encompassing all of each species' limited distribution. Each of these species lack forelimbs, and differ by the number of digits on the hind limbs, among other morphological characters; S. gronovii bears a single digit and S. kasneri bears two digits on the hind limbs. Sequence data obtained from three mtDNA (16S ribosomal RNA, cytochrome b, and nicotinamide adenine dinucleotide dehydrogenase 1 unit; 2035 bp ttl.) and two nuclear (dynein axonemal heavy chain 3 and the natural killer tumor recognition; 1848 bp ttl.) gene regions were used to reconstruct the evolutionary relationships among the two focal species and several other co-distributed species (Scelotes bipes, Scelotes montispectus, and Scelotes sexlineatus). Phylogenetic results (Bayesian and parsimony) revealed that several populations previously considered S. kasneri actually belong to other species, and others are paraphyletic with respect to one another. Additionally, populations of S. gronovii were also found to be paraphyletic, with populations south of the Berg River supported as sister to S. bipes, and populations north of the Berg River sister the remaining sampled species. Our results require a redefinition of S. sexlineatus to encompass populations morphologically convergent with S. kasneri and restrict the ranges of the already threatened S. kasneri and S. gronovii even further. The paraphyly of S. gronovii and the placement of each clade as sister to clades of species bearing two digits on the hind limbs suggests that digit loss has occurred at least twice in this group., (Copyright © 2011. Published by Elsevier Inc.)

Iguanian lizards form a diverse clade whose members have been the focus of many comparative studies of ecology, behavior, and evolution. Despite the importance of phylogeny to such studies, interrelationships among many iguanian clades remain uncertain. Within the Old World clade Acrodonta, Agamidae is sometimes found to be paraphyletic with respect to Chamaeleonidae, and recent molecular studies have produced conflicting results for many major clades. Within the largely New World clade Pleurodonta, relationships among the 12 currently recognized major subclades (mostly ranked as families) have been largely unresolved or poorly supported in previous studies. To clarify iguanian evolutionary history, we first infer phylogenies using concatenated maximum-likelihood (ML) and Bayesian analyses of DNA sequence data from 29 nuclear protein-coding genes for 47 iguanian and 29 outgroup taxa. We then estimate a relaxed-clock Bayesian chronogram for iguanians using BEAST. All three methods produce identical topologies. Within Acrodonta, we find strong support for monophyly of Agamidae with respect to Chamaeleonidae, and for almost all relationships within agamids. Within Pleurodonta, we find strong Bayesian support for almost all relationships, and strong ML support for some interfamilial relationships and for monophyly of almost all families (excepting Polychrotidae). Our phylogenetic results suggest a non-traditional biogeographic scenario in which pleurodonts originated in the Northern Hemisphere and subsequently spread southward into South America. The pleurodont portion of the tree is characterized by several very short, deep branches, raising the possibility of deep coalescences that may confound concatenated analyses. We therefore also use 27 of these genes to implement a coalescent-based species-tree approach for pleurodonts. Although this analysis strongly supports monophyly of the pleurodont families, interfamilial relationships are generally different from those in the concatenated tree, and support is uniformly poor. However, a species-tree analysis using only the seven most variable loci yields higher support and more congruence with the concatenated tree. This suggests that low support in the 27-gene species-tree analysis may be an artifact of the many loci that are uninformative for very short branches. This may be a general problem for the application of species-tree methods to rapid radiations, even with phylogenomic data sets. Finally, we correct the non-monophyly of Polychrotidae by recognizing the pleurodont genus Anolis (sensu lato) as a separate family (Dactyloidae), and we correct the non-monophyly of the agamid genus Physignathus by resurrection of the genus Istiurus for the former Physignathus lesueurii., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has caused dramatic amphibian population declines and extinctions in Australia, Central and North America, and Europe. Bd is associated with >200 species extinctions of amphibians, but not all species that become infected are susceptible to the disease. Specifically, Bd has rapidly emerged in some areas of the world, such as in Australia, USA, and throughout Central and South America, causing population and species collapse. The mechanism behind the rapid global emergence of the disease is poorly understood, in part due to an incomplete picture of the global distribution of Bd. At present, there is a considerable amount of geographic bias in survey effort for Bd, with Asia being the most neglected continent. To date, Bd surveys have been published for few Asian countries, and infected amphibians have been reported only from Indonesia, South Korea, China and Japan. Thus far, there have been no substantiated reports of enigmatic or suspected disease-caused population declines of the kind that has been attributed to Bd in other areas. In order to gain a more detailed picture of the distribution of Bd in Asia, we undertook a widespread, opportunistic survey of over 3,000 amphibians for Bd throughout Asia and adjoining Papua New Guinea. Survey sites spanned 15 countries, approximately 36° latitude, 111° longitude, and over 2000 m in elevation. Bd prevalence was very low throughout our survey area (2.35% overall) and infected animals were not clumped as would be expected in epizootic events. This suggests that Bd is either newly emerging in Asia, endemic at low prevalence, or that some other ecological factor is preventing Bd from fully invading Asian amphibians. The current observed pattern in Asia differs from that in many other parts of the world.

Molecular data offer great potential to resolve the phylogeny of living taxa but can molecular data improve our understanding of relationships of fossil taxa? Simulations suggest that this is possible, but few empirical examples have demonstrated the ability of molecular data to change the placement of fossil taxa. We offer such an example here. We analyze the placement of snakes among squamate reptiles, combining published morphological data (363 characters) and new DNA sequence data (15,794 characters, 22 nuclear loci) for 45 living and 19 fossil taxa. We find several intriguing results. First, some fossil taxa undergo major changes in their phylogenetic position when molecular data are added. Second, most fossil taxa are placed with strong support in the expected clades by the combined data Bayesian analyses, despite each having >98% missing cells and despite recent suggestions that extensive missing data are problematic for Bayesian phylogenetics. Third, morphological data can change the placement of living taxa in combined analyses, even when there is an overwhelming majority of molecular characters. Finally, we find strong but apparently misleading signal in the morphological data, seemingly associated with a burrowing lifestyle in snakes, amphisbaenians, and dibamids. Overall, our results suggest promise for an integrated and comprehensive Tree of Life by combining molecular and morphological data for living and fossil taxa.

Plate tectonics can have profound effects on organismal distribution and is often the driving force in speciation. Through geologic processes, the Baja California Peninsula depicts two faunal patterns: one through southern vicariance with Cape separation, and the other through dispersal onto the northern peninsula, referred to as a 'dual-peninsular effect.' Here we apply a hierarchical sampling strategy that combines population-level sequence data ( approximately 800bp, nad4 region) with complete mt-genome data (aligned 15,549bp) and 5 nuclear protein encoding loci (3315bp), to test whether both patterns have occurred in one group of nightsnakes (Hypsiglena). The geologic formation of the peninsula is thought to have occurred in three stages: (1) Cape separation from mainland Mexico; (2) the northern peninsula separated, forming the northern Gulf of California; and (3) the peninsula was united through volcanic activity, while moving northward causing collision with southern California. However, the timing of events is debated. We explore phylogenetic relationships and estimate dates of divergence for nightsnakes using our hierarchical sampling strategy. Our data support both 'southern-vicariance' and 'northern-dispersal' have occurred in nightsnakes, forming a ring distribution around the Gulf of California. Two divergent forms are sympatric on the southern half of the peninsula with no indication of hybridization. Nightsnakes represent the first group to depict the 'dual-peninsular effect' with extensive overlap on the Baja California Peninsula.

Many authors have claimed that short branches in the Tree of Life will be very difficult to resolve with strong support, even with the large multilocus data sets now made possible by genomic resources. Short branches may be especially problematic because the underlying gene trees are expected to have discordant phylogenetic histories when the time between branching events is very short. Although there are many examples of short branches that are difficult to resolve, surprisingly, no empirical studies have systematically examined the relationships between branch lengths, branch support, and congruence among genes. Here, we examine these fundamental relationships quantitatively using a data set of 20 nuclear loci for 50 species of snakes (representing most traditionally recognized families). A combined maximum likelihood analysis of the 20 loci gives strong support for 69% of the nodes, but many remain weakly supported, with bootstrap values for 20% ranging from 21% to 66%. For the combined-data tree, we find significant correlations between the length of a branch, levels of bootstrap support, and the proportion of genes that are congruent with that branch in the separate analyses of each gene. We also find that strongly supported conflicts between gene trees over the resolution of individual branches are common (roughly 35% of clades), especially for shorter branches. Overall, our results support the hypothesis that short branches may be very difficult to confidently resolve, even with large, multilocus data sets. Nevertheless, our study provides strong support for many clades, including several that were controversial or poorly resolved in previous studies of snake phylogeny.

The subspecies concept has received considerable debate throughout the past century. Subspecies were originally used to delineate potential incipient species, but were later employed to simply capture geographical variation. There is a recent trend to eliminate the trinomial in light of new evidence. Discrete, diagnosable lineages are elevated to specific status, while those that show clinal variation and/or appear to represent ecological pattern classes are placed in synonymy with the parent species and the subspecific epithets are disregarded. Here, I examine the species boundaries of nightsnakes (Hypsiglena torquata) using standard phylogeographic methods and mtDNA data from 178 individuals. Previously, seventeen subspecies of H. torquata were described. In this study, I recognize six species in what was previously considered H. torquata: one is novel, two were previously recognized subspecies, while the remaining three are wide-spread, polymorphic lineages, composed of multiple subspecies. I make the case to maintain the subspecific lineages in these wide-ranging species because they are geographically cohesive, morphologically discrete, and may represent incipient species within each complex, which have not yet achieved speciation. These subspecies are maintained, not only pending future investigations, but because they provide a useful identity for the taxonomy of this diverse lineage.

The broad distribution of the Sceloporus magister species group (squamata: phrynosomatidae) throughout western North America provides an appropriate model for testing biogeographical hypotheses explaining the timing and origins of diversity across mainland deserts and the Baja California Peninsula. We inferred concordant phylogenetic trees describing the higher-level relationships within the magister group using 1.6 kb of mitochondrial DNA (mtDNA) and 1.7 kb of nuclear DNA data. These data provide strong support for the parallel divergence of lineages endemic to the Baja California Peninsula (S. zosteromus and the orcutti complex) in the form of two sequential divergence events at the base of the magister group phylogeny. A relaxed phylogenetic analysis of the mtDNA data using one fossil and one biogeographical constraint provides a chronology of these divergence events and evidence that further diversification within the Baja California clades occurred simultaneously, although patterns of geographical variation and speciation between clades differ. We resolved four major phylogeographical clades within S. magister that (i) provide a novel phylogenetic placement of the Chihuahuan Desert populations sister to the Mojave Desert; (ii) illustrate a mixed history for the Colorado Plateau that includes Mojave and Sonoran Desert components; and (iii) identify an area of overlap between the Mojave and Sonoran Desert clades near Yuma, Arizona. Estimates of bidirectional migration rates among populations of S. magister using four nuclear loci support strong asymmetries in gene flow among the major mtDNA clades. Based on the nonexclusivity of mtDNA haplotypes, nuclear gene flow among populations and wide zones of phenotypic intergradation, S. magister appears to represent a single geographically variable and widespread species.

Two species of horned lizards are sympatric along the periphery of the Salton Trough. Phrynosoma mcallii, endemic to the trough, is of conservation concern because its limited habitat has been fragmented by human activities. A more common and widespread species, Phrynosoma platyrhinos, occurs around the periphery of the trough and much further to the North. The two species are syntopic at a few localities, where morphological intermediates have also been found. Here, we used nested clade phylogeographical analysis (NCPA) and analysis of molecular variance (amova), to examine 781 bp of mitochondrial DNA (nad4 and two tRNAs) from 82 individuals of P. mcallii. We tested whether populations of this species were recently connected, or if they were historically isolated prior to human modification of the region. Our NCPA results indicated significant population structure associated with the Colorado River, suggesting limited gene flow and potential isolation across this barrier. Populations west of the Colorado River, currently isolated from one another by human development, show less genetic differentiation. We also collected homologous sequence data from 34 individuals of P. platyrhinos and seven specimens morphologically intermediate between P. mcallii and P. platyrhinos, as a preliminary investigation of hybridization between these two species. From phylogenetic results of these data, we identified a species (Phrynosoma goodei) previously recognized as a subspecies of P. platyrhinos. Six of the morphologically intermediate specimens shared mtDNA haplotypes with P. goodei, while one was nested among P. mcallii haplotypes.

Parallel "arms races" involving the same or similar phenotypic interfaces allow inference about selective forces driving coevolution, as well as the importance of phylogenetic and phenotypic constraints in coevolution. Here, we report the existence of apparent parallel arms races between species pairs of garter snakes and their toxic newt prey that indicate independent evolutionary origins of a key phenotype in the interface. In at least one area of sympatry, the aquatic garter snake, Thamnophis couchii, has evolved elevated resistance to the neurotoxin tetrodotoxin (TTX), present in the newt Taricha torosa. Previous studies have shown that a distantly related garter snake, Thamnophis sirtalis, has coevolved with another newt species that possesses TTX, Taricha granulosa. Patterns of within population variation and phenotypic tradeoffs between TTX resistance and sprint speed suggest that the mechanism of resistance is similar in both species of snake, yet phylogenetic evidence indicates the independent origins of elevated resistance to TTX.

The widespread, lowland toad Bufo valliceps has an unusual distribution in North and Middle America that straddles two major biogeographic areas; previous morphological studies of this species suggested the existence of two species. We used mitochondrial DNA sequences to examine the phylogeography of this species and discovered the existence of two distinct clades. We recognize these as two species: B. valliceps and B. nebulifer. These molecular data support morphological data from previous studies. Our results show low levels of molecular variation in a morphologically uniform temperate species (B. nebulifer) and high levels of molecular variation in a morphologically variable tropical species (B. valliceps), providing an example of molecules matching morphology. Two biogeographic hypotheses are tested to explain the current distribution of these species, based on a calibrated rate of evolution and the percent sequence divergence between the two species. A more recent Pleistocene dispersal event, followed by vicariance associated with rising sea level, is rejected in favor of an earlier Miocene-Pliocene vicariant hypothesis associated with the formation of the Trans-Mexican Neovolcanic Belt., (Copyright 2000 Academic Press.)