Fossil-informed biogeographic analysis suggests Eurasian regionalization in crown Squamata during the early Jurassic

Background

Squamata (lizards, snakes, and amphisbaenians) is a Triassic lineage with an extensive and complex biogeographic history, yet no large-scale study has reconstructed the ancestral range of early squamate lineages. The fossil record indicates a broadly Pangaean distribution by the end- Cretaceous, though many lineages (e.g., Paramacellodidae, Mosasauria, Polyglyphanodontia) subsequently went extinct. Thus, the origin and occupancy of extant radiations is unclear and may have been localized within Pangaea to specific plates, with potential regionalization to distinct Laurasian and Gondwanan landmasses during the Mesozoic in some groups.

Methods

We used recent tectonic models to code extant and fossil squamate distributions occurring on nine discrete plates for 9,755 species, with Jurassic and Cretaceous fossil constraints from three extinct lineages. We modeled ancestral ranges for crown Squamata from an extant-only molecular phylogeny using a suite of biogeographic models accommodating different evolutionary processes and fossil-based node constraints from known Jurassic and Cretaceous localities. We hypothesized that the best-fit models would not support a full Pangaean distribution (i.e., including all areas) for the origin of crown Squamata, but would instead show regionalization to specific areas within the fragmenting supercontinent, likely in the Northern Hemisphere where most early squamate fossils have been found.

Results

Incorporating fossil data reconstructs a localized origin within Pangaea, with early regionalization of extant lineages to Eurasia and Laurasia, while Gondwanan regionalization did not occur until the middle Cretaceous for Alethinophidia, Scolecophidia, and some crown Gekkotan lineages. While the Mesozoic history of extant squamate biogeography can be summarized as a Eurasian origin with dispersal out of Laurasia into Gondwana, their Cenozoic history is complex with multiple events (including secondary and tertiary recolonizations) in several directions. As noted by previous authors, squamates have likely utilized over-land range expansion, land-bridge colonization, and trans-oceanic dispersal. Tropical Gondwana and Eurasia hold more ancient lineages than the Holarctic (Rhineuridae being a major exception), and some asymmetries in colonization (e.g., to North America from Eurasia during the Cenozoic through Beringia) deserve additional study. Future studies that incorporate fossil branches, rather than as node constraints, into the reconstruction can be used to explore this history further.

The Draft Genome Sequences of 50 Salamander species (Caudata, Amphibia)

We present partial genome sequences of 50 salamander species (Urodela) from 10 genera and 4 families. These span nearly the entire range of genome sizes in salamanders, from ~14-130GB, the latter of which is among the largest of all animal genomes. Only three salamander genomes were available to this point, from Ambystomatidae (one species) and Salamandridae (two species from two genera), to which we have added Amphiumidae (one species), Plethodontidae (45 species from 6 genera), Proteidae (one species), and Sirenidae (three species from two genera). These span ~140 million years of evolutionary divergence, leaving only Cryptobranchidae, Hynobiidae, and Rhyacotritonidae as salamander families without genome assemblies. These data should facilitate additional future work on speciation and genome evolution, both within Urodela and across Animalia.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

A giant tyrannosaur from the Campanian-Maastrichtian of southern North America and the evolution of tyrannosaurid gigantism

Tyrannosaurid dinosaurs dominated as predators in the Late Cretaceous of Laurasia, culminating in the evolution of the giant Tyrannosaurus rex, both the last and largest tyrannosaurid. Where and when Tyrannosaurini (T. rex and kin) originated remains unclear. Competing hypotheses place tyrannosaurin origins in Asia, or western North America (Laramidia). We report a new tyrannosaurin, Tyrannosaurus mcraeensis, from the Campanian-Maastrichtian Hall Lake Formation of New Mexico, based on a fossil previously referred to T. rex. T. mcraeensis predates T. rex by ~ 6-7 million years, yet rivaled it in size. Phylogenetic analysis recovers T. mcraeensis as sister to T. rex and suggests Tyrannosaurini originated in southern Laramidia. Evolution of giant tyrannosaurs in southern North America, alongside giant ceratopsians, hadrosaurs, and titanosaurs suggests large-bodied dinosaurs evolved at low latitudes in North America.

A Dataset of Amphibian Species in U.S. National Parks

National parks and other protected areas are important for preserving landscapes and biodiversity worldwide. An essential component of the mission of the United States (U.S.) National Park Service (NPS) requires understanding and maintaining accurate inventories of species on protected lands. We describe a new, national-scale synthesis of amphibian species occurrence in the NPS system. Many park units have a list of amphibian species observed within their borders compiled from various sources and available publicly through the NPSpecies platform. However, many of the observations in NPSpecies remain unverified and the lists are often outdated. We updated the amphibian dataset for each park unit by collating old and new park-level records and had them verified by regional experts. The new dataset contains occurrence records for 292 of the 424 NPS units and includes updated taxonomy, international and state conservation rankings, hyperlinks to a supporting reference for each record, specific notes, and related fields which can be used to better understand and manage amphibian biodiversity within a single park or group of parks.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Early origin and diverse phenotypic implementation of iridescent UV patterns for sexual signaling in pierid butterflies

Iridescent ultraviolet (IUV) patterns on pierid butterfly wings are phenotypic adaptations commonly used as sexual signals, generated by scales with ultrastructural modifications. Pierid IUV patterns are sexually dichromatic, with reduced size in females, where conspicuous sexual signaling balances courtship against ecological predation. There have been no phylogenetic reconstructions of IUV within Pieridae and little morphological characterization of phenotypic diversity. Our genus-wide characterization of IUV revealed the uniform similarity of stacked lamellar ridges on the dorsal surface of cover scales. We tested a hypothesis of single versus multiple origins by reconstructing a phylogeny of 534 species (~43.2% described species), with all genera represented, and a trait matrix of 734 species (~59.4%) screened for IUV. A single, early dimorphic origin of IUV followed by several losses and gains received strong support, concluding that IUV patterns and structural coloration are old traits. Collectively, these results support the homology of IUV scales and patterns that diversified within several lineages, suggesting an interplay between female-mediated sexual selection and ecological predatory selection.

Unsupervised machine learning for species delimitation, integrative taxonomy, and biodiversity conservation

Integrative taxonomy, combining data from multiple axes of biologically relevant variation, is a major goal of systematics. Ideally, such taxonomies will derive from similarly integrative species-delimitation analyses. Yet, most current methods rely solely or primarily on molecular data, with other layers often incorporated only in a post hoc qualitative or comparative manner. A major limitation is the difficulty of devising quantitative parametric models linking different datasets in a unified ecological and evolutionary framework. Machine Learning (ML) methods offer flexibility in this arena by easily learning high-dimensional associations between observations (e.g., individual specimens) across a wide array of input features (e.g., genetics, geography, environment, and phenotype) to delimit statistically meaningful clusters. Here, I implement an unsupervised method using Self-Organizing (or "Kohonen") Maps (SOMs) for such purposes. Recent extensions called "SuperSOMs" can integrate multiple layers, each of which exerts independent influence on a two-dimensional output grid via empirically estimated weights. The grid cells are then delimited into K distinct units that can be interpreted as species or other entities. I show empirical examples in salamanders (Desmognathus) and snakes (Storeria) with layers representing alleles, space, climate, and traits. Simulations reveal that the SuperSOM approach can detect K = 1, tends not to over-split, reflects contributions from all layers, and limits large layers (e.g., genetic matrices) from overwhelming other datasets, desirable properties addressing major concerns from previous studies. Finally, I suggest that these and similar methods could integrate conservation-relevant layers such as population trends and human encroachment to delimit management units from an explicitly quantitative framework grounded in the ecology and evolution of species limits and boundaries.

Systematic revision of the Spotted and Northern Dusky Salamanders (Plethodontidae: Desmognathus conanti and D. fuscus), with six new species from the eastern United States

Spotted and Northern Dusky Salamanders (Desmognathus conanti and D. fuscus) have a long and complex taxonomic history. At least 10 other currently recognized species in the genus were either described from populations previously considered D. fuscus, described as or later considered subspecies thereof, or later considered synonyms thereof, before ultimately being recognized as distinct. Recent molecular analyses have also revealed extensive cryptic diversity within both species, which are polyphyletic assemblages of 13 distinct mitochondrial lineages with 5.7-10.3% uncorrected 'p' distances in the COI barcode locus. Based on phylogenomic data and population-clustering analyses considering admixture between lineages, 11 candidate species were circumscribed by recent authors. Those within D. conanti are also ecomorphologically variable, comprising both large, robust, keel-tailed populations, and small, gracile, round-tailed forms. Evaluating their distinctiveness based on genetic, geographic, and morphological evidence, we conclude that six of the candidates represent new species: Desmognathus anicetus sp. nov., D. bairdi sp. nov., D. campi sp. nov., D. catahoula sp. nov., D. lycos sp. nov., and D. tilleyi sp. nov. Consequently, we recognize eight total species from populations formerly associated with the nominal species D. conanti and D. fuscus, the re-delimited concepts of which also contain additional phylogeographic lineage diversity that may represent further distinct species. In addition to existing mitochondrial and nuclear phylogenetic, network, and clustering results, we present preliminary analyses of linear morphometrics to bolster diagnostic specificity based on phenotypic characteristics. These changes stabilize the previously paraphyletic taxonomy of species-level lineages within Desmognathus, though additional cryptic diversity may exist both within the species considered here, and elsewhere in the genus.

A systematic revision of the Shovel-nosed Salamander (Plethodontidae: Desmognathus marmoratus), with re-description of the related D. aureatus and D. intermedius

Shovel-nosed Salamanders, Desmognathus marmoratus (Moore, 1899), were long thought to represent a single species from the southern Appalachian Mountains of the eastern United States, ranging from northeastern Georgia to extreme southwestern Virginia. These populations have a highly derived ecomorphology, being fully aquatic with a specialized flattened and elongated phenotype adapted to rocky riffle zones in fast-flowing, high-gradient mountain streams. Because of this, they were originally described in a separate genus, Leurognathus Moore, 1899. Four additional species or subspecies were described from 1928-1956 based on regional geographic variation in phenotype before being synonymized with L. marmoratus in 1962, which was reassigned to Desmognathus in 1996. Molecular analyses subsequently revealed four distinct candidate lineages in two distantly related clades, which were recently re-delimited into three species. These are D. aureatus (Martof, 1956) from northeastern Georgia, D. intermedius (Pope, 1928) from western North Carolina, and D. marmoratus from northwestern North Carolina. We provide a systematic revision of these taxa, which do not represent a natural group but instead exhibit convergent phenotypes across multiple species, potentially driven by ancient episodes of adaptive introgression between ancestral lineages. Our recent fieldwork revealed an astonishingly disjunct and morphologically distinct population of D. marmoratus in the New River Gorge of West Virginia, which were previously confused with D. kanawha Pyron and Beamer, 2022. This locality is ~120 airline km away from the nearest populations of D. marmoratus in Virginia. No Shovel-nosed Salamanders have ever been found in the New River drainage during our extensive previous explorations or credibly reported in museum specimens or the literature. Additional cryptic populations of these taxa may remain.

A new species of Dipsas (Serpentes, Dipsadidae) from central Panama

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 Dipsas temporalis (Werner, 1909), the sister species, is updated. Finally, a key to the species of Dipsas currently known from Middle America is presented.

Bald sea urchin disease shifts the surface microbiome on purple sea urchins in an aquarium

Bald sea urchin disease (BSUD) is most likely a bacterial infection that occurs in a wide range of sea urchin species and causes the loss of surface appendages. The disease has a variety of additional symptoms, which may be the result of the many bacteria that are associated with BSUD. Previous studies have investigated causative agents of BSUD, however, there are few reports on the surface microbiome associated with the infection. Here, we report changes to the surface microbiome on purple sea urchins in a closed marine aquarium that contracted and then recovered from BSUD in addition to the microbiome of healthy sea urchins in a separate aquarium. 16S rRNA gene sequencing shows that microhabitats of different aquaria are characterized by different microbial compositions, and that diseased, recovered, and healthy sea urchins have distinct microbial compositions, which indicates that there is a correlation between microbial shifts and recovery from disease.

Speciation Hypotheses from Phylogeographic Delimitation Yield an Integrative Taxonomy for Seal Salamanders (Desmognathus monticola)

Significant advances have been made in species delimitation and numerous methods can test precisely defined models of speciation, though the synthesis of phylogeography and taxonomy is still sometimes incomplete. Emerging consensus treats distinct genealogical clusters in genome-scale data as strong initial evidence of speciation in most cases; a hypothesis that must therefore be falsified under an explicit evolutionary model. We can now test speciation hypotheses linking trait differentiation to specific mechanisms of divergence with increasingly large datasets. Integrative taxonomy can therefore reflect an understanding of how each axis of variation relates to underlying speciation processes, with nomenclature for distinct evolutionary lineages. We illustrate this approach here with Seal Salamanders (Desmognathus monticola) and introduce a new unsupervised machine-learning approach for species delimitation. Plethodontid salamanders are renowned for their morphological conservatism despite extensive phylogeographic divergence. We discover two geographic genetic clusters, for which demographic and spatial models of ecology and gene flow provide robust support for ecogeographic speciation despite limited phenotypic divergence. These data are integrated under evolutionary mechanisms (e.g., spatially localized gene flow with reduced migration) and reflected in emergent properties expected under models of reinforcement (e.g., ethological isolation and selection against hybrids). Their genetic divergence is prima facie evidence for species-level distinctiveness, supported by speciation models and divergence along axes such as behavior, geography, and climate that suggest an ecological basis with subsequent reinforcement through prezygotic isolation. As datasets grow more comprehensive, species delimitation models can be tested, rejected, or corroborated as explicit speciation hypotheses, providing for reciprocal illumination of evolutionary processes and integrative taxonomies.

Evolutionary legacies in contemporary tetrapod imperilment

The Tree of Life will be irrevocably reshaped as anthropogenic extinctions continue to unfold. Theory suggests that lineage evolutionary dynamics, such as age since origination, historical extinction filters and speciation rates, have influenced ancient extinction patterns - but whether these factors also contribute to modern extinction risk is largely unknown. We examine evolutionary legacies in contemporary extinction risk for over 4000 genera, representing ~30,000 species, from the major tetrapod groups: amphibians, birds, turtles and crocodiles, squamate reptiles and mammals. We find consistent support for the hypothesis that extinction risk is elevated in lineages with higher recent speciation rates. We subsequently test, and find modest support for, a primary mechanism driving this pattern: that rapidly diversifying clades predominantly comprise range-restricted, and extinction-prone, species. These evolutionary patterns in current imperilment may have important consequences for how we manage the erosion of biological diversity across the Tree of Life.

Amphibian Speciation Rates Support a General Role of Mountains as Biodiversity Pumps

AbstractContinental mountain areas cover <15% of global land surface, yet these regions concentrate >80% of global terrestrial diversity. One prominent hypothesis to explain this pattern proposes that high mountain diversities could be explained by higher diversification rates in regions of high topographic complexity (HTC). While high speciation in mountains has been detected for particular clades and regions, the global extent to which lineages experience faster speciation in mountains remains unknown. Here we addressed this issue using amphibians as a model system (>7,000 species), and we found that families showing high speciation rates contain a high proportion of species distributed in mountains. Moreover, we found that lineages inhabiting areas of HTC speciate faster than lineages occupying areas that are topographically less complex. When comparing across regions, we identified the same pattern in five biogeographical realms where higher speciation rates are associated with higher levels of complex topography. Low-magnitude differences in speciation rates between some low and high complex topographies suggest that high mountain diversity is also affected by low extinction and/or high colonization rates. Nevertheless, our results bolster the importance of mountains as engines of speciation at different geographical scales and highlight their importance for the conservation of global biodiversity.

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Phylogenetic and spatial distribution of evolutionary diversification, isolation, and threat in turtles and crocodilians (non-avian archosauromorphs)

Background

The origin of turtles and crocodiles and their easily recognized body forms dates to the Triassic and Jurassic. Despite their long-term success, extant species diversity is low, and endangerment is extremely high compared to other terrestrial vertebrate groups, with ~ 65% of ~ 25 crocodilian and ~ 360 turtle species now threatened by exploitation and habitat loss. Here, we combine available molecular and morphological evidence with statistical and machine learning algorithms to present a phylogenetically informed, comprehensive assessment of diversification, threat status, and evolutionary distinctiveness of all extant species.

Results

In contrast to other terrestrial vertebrates and their own diversity in the fossil record, the recent extant lineages of turtles and crocodilians have not experienced any global mass extinctions or lineage-wide shifts in diversification rate or body-size evolution over time. We predict threat statuses for 114 as-yet unassessed or data-deficient species and identify a concentration of threatened turtles and crocodilians in South and Southeast Asia, western Africa, and the eastern Amazon. We find that unlike other terrestrial vertebrate groups, extinction risk increases with evolutionary distinctiveness: a disproportionate amount of phylogenetic diversity is concentrated in evolutionarily isolated, at-risk taxa, particularly those with small geographic ranges. Our findings highlight the important role of geographic determinants of extinction risk, particularly those resulting from anthropogenic habitat-disturbance, which affect species across body sizes and ecologies.

Conclusions

Extant turtles and crocodilians maintain unique, conserved morphologies which make them globally recognizable. Many species are threatened due to exploitation and global change. We use taxonomically complete, dated molecular phylogenies and various approaches to produce a comprehensive assessment of threat status and evolutionary distinctiveness of both groups. Neither group exhibits significant overall shifts in diversification rate or body-size evolution, or any signature of global mass extinctions in recent, extant lineages. However, the most evolutionarily distinct species tend to be the most threatened, and species richness and extinction risk are centered in areas of high anthropogenic disturbance, particularly South and Southeast Asia. Range size is the strongest predictor of threat, and a disproportionate amount of evolutionary diversity is at risk of imminent extinction.

Specialized breeding in plants affects diversification trajectories in Neotropical frogs

Many animals breed exclusively in plants that accumulate rainwater (phytotelma; e.g., bromeliad, bamboo, fruit husk, and tree hole), to which they are either physiologically or behaviorally specialized for this microhabitat. Of the numerous life-history modes observed in frogs, few are as striking or potentially consequential as the transition from pond- or stream-breeding to the deposition of eggs or larvae in phytotelmata. Such specialization can increase offspring survivorship due to reduced competition and predation, representing potential ecological opportunity for adaptive radiation, yet few lineages of phytotelma-breeding frogs appear to have diversified extensively after such a transition, at least in the New World. We use a phylogeny of Neotropical frogs and data on breeding microhabitat to understand the evolutionary transitions involved with specialized phytotelma-breeding. First, we find that phytotelma-breeding is present in at least 168 species in 10 families of frogs. Across the phylogeny, we estimate 14 origins of phytotelma-breeding and 115 reversals, showing that phytotelma-breeding is a highly labile character. Second, phytotelma-breeding lineages overall have higher net diversification than nonphytotelma-breeding ones. This specialization represents an ecological opportunity resulting in increased diversification in most families with phytotelma-breeding lineages, whereas phytotelma-breeding toads have restricted diversification histories.

Availability of Recently Described Ptychadena (Anura: Ptychadenidae) Nomina from Ethiopia

Frogs of the genus Ptychadena (Boulenger, 1917) have long been identified as harboring cryptic diversity and comprising numerous species-complexes (Largen 1997; Zimkus et al. 2017), and many authors have recognized the particularly high hidden richness in the Ethiopian highlands (Largen 1997; see Largen Spawls 2010 and refs. within). This cryptic diversity was confirmed by recent molecular studies (Freilich et al. 2014; Smith et al. 2017a, Reyes-Velasco et al. 2018). Those authors identified a congruent set of evolutionarily distinct candidate species using both mitochondrial and nuclear DNA, and described the geographic and ecological isolation of these species in detail (Freilich et al. 2014; Smith et al. 2017a).

A new species of African Snake-Eyed Skink (Scincidae: Panaspis) from Ethiopia

A recent molecular phylogenetic revision of the snake-eyed skinks (genus Panaspis Cope, 1868) uncovered extensive cryptic diversity, including several new species from throughout sub-Saharan Africa. Here, we describe one of these from Ethiopia as Panaspis annettesabinae sp. nov. This description is based on a previous molecular phylogeny and morphological, scalation, and coloration data collected from the type specimen. Phylogenetic analyses place the species alone in what we term the P. annettesabinae species group from Ethiopia. This group forms the sister lineage to a large southern African radiation and suggests a potential northern origin for much of the extant diversity of Panaspis. Many new taxa have recently been discovered in the genus and region, and there are several historical records of Panaspis from elsewhere in Ethiopia. Thus, we suggest that the range of this new species (known only from a single specimen at present) may be much larger, and that additional undescribed species may exist in northern sub-Saharan Africa.

Phylogenetic and Trait-Based Prediction of Extinction Risk for Data-Deficient Amphibians

Amphibians are among the most highly threatened lineages, with at least 2,000 species estimated to be in danger of extinction [1, 2]. Alarmingly, another ∼2,200 species (∼25% of all ∼7,900 known species) are data deficient or not evaluated (hereinafter termed data deficient) by the International Union for Conservation of Nature (IUCN) [1]. Without an estimate of their status, data-deficient species are usually overlooked in conservation planning and resource allocation [3]. Amphibians have the highest proportion of data-deficient species of any vertebrate group [1, 4], which highlights the need to estimate their threat status considering potentially imminent extinctions. We apply a trait-based spatio-phylogenetic statistical framework [5] to predict threat status for data-deficient species. Because ecological, geographical, and evolutionary attributes increase extinction risk [6, 7], we used geographic distribution data [1, 8], phylogenetically imputed ecological traits, and an amphibian phylogeny [9] to provide initial baseline predictions. We estimate that half of the ∼2,200 data-deficient species are threatened with extinction (vulnerable, endangered, or critically endangered), primarily in the Neotropics and Southeast Asia. This increases the number of amphibian species estimated to be threatened with extinction by ∼50%. Of these, we predict that ∼500 species are endangered or critically endangered, and three may be extinct already. We highlight families that are most at risk and suggest where urgent conservation is needed to avert their loss. We show that some of the most vulnerable species may also be the most poorly known and offer an analytical framework for preliminary analysis of their threat status in the face of deficient empirical data.

Additions to the phylogeny of colubrine snakes in Southwestern Asia, with description of a new genus and species (Serpentes: Colubridae: Colubrinae)

Reptiles are still being described worldwide at a pace of hundreds of species a year. While many discoveries are from remote tropical areas, biodiverse arid regions still harbor many novel taxa. Here, we present an updated phylogeny of colubrid snakes from the Western Palearctic by analyzing a supermatrix of all available global snake species with molecular data and report on the discovery of a new genus and species of colubrine snake from southeastern Iran. The new taxon, named Persiophis fahimii Gen. et sp. nov., is nested within a clade containing Middle Eastern and South Asian ground racers (Lytorhynchus, Rhynchocalamus, Wallaceophis, and Wallophis). This species has a derived morphology including an edentulous pterygoid and occurrence of short and blunt teeth on the palatine, maxillae and dentary bones, an elongated snout and a relatively trihedral first supralabial scale that is slightly bigger than the second, and elongated toward the tip of rostral. We also report on the osteology and phylogenetic placement of several poorly studied colubrines: Hierophis andreanus (reassigned to Dolichophis) and Muhtarophis barani.

Interrogating Genomic-Scale Data for Squamata (Lizards, Snakes, and Amphisbaenians) Shows no Support for Key Traditional Morphological Relationships

Genomics is narrowing uncertainty in the phylogenetic structure for many amniote groups. For one of the most diverse and species-rich groups, the squamate reptiles (lizards, snakes, and amphisbaenians), an inverse correlation between the number of taxa and loci sampled still persists across all publications using DNA sequence data and reaching a consensus on the relationships among them has been highly problematic. In this study, we use high-throughput sequence data from 289 samples covering 75 families of squamates to address phylogenetic affinities, estimate divergence times, and characterize residual topological uncertainty in the presence of genome-scale data. Importantly, we address genomic support for the traditional taxonomic groupings Scleroglossa and Macrostomata using novel machine-learning techniques. We interrogate genes using various metrics inherent to these loci, including parsimony-informative sites (PIS), phylogenetic informativeness, length, gaps, number of substitutions, and site concordance to understand why certain loci fail to find previously well-supported molecular clades and how they fail to support species-tree estimates. We show that both incomplete lineage sorting and poor gene-tree estimation (due to a few undesirable gene properties, such as an insufficient number of PIS), may account for most gene and species-tree discordance. We find overwhelming signal for Toxicofera, and also show that none of the loci included in this study supports Scleroglossa or Macrostomata. We comment on the origins and diversification of Squamata throughout the Mesozoic and underscore remaining uncertainties that persist in both deeper parts of the tree (e.g., relationships between Dibamia, Gekkota, and remaining squamates; among the three toxicoferan clades Iguania, Serpentes, and Anguiformes) and within specific clades (e.g., affinities among gekkotan, pleurodont iguanians, and colubroid families).

Systematics of South American snail-eating snakes (Serpentes, Dipsadini), with the description of five new species from Ecuador and Peru

A molecular phylogeny of the Neotropical snail-eating snakes (tribe Dipsadini) is presented including 43 (24 for the first time) of the 77 species, sampled for both nuclear and mitochondrial genes. Morphological and phylogenetic support was found for four new species of Dipsas and one of Sibon, which are described here based on their unique combination of molecular, meristic, and color pattern characteristics. Sibynomorphus is designated as a junior subjective synonym of Dipsas. Dipsas latifrontalis and D. palmeri are resurrected from the synonymy of D. peruana. Dipsas latifasciata is transferred from the synonymy of D. peruana to the synonymy of D. palmeri. A new name, D. jamespetersi, is erected for the taxon currently known as Sibynomorphus petersi. Re-descriptions of D. latifrontalis and D. peruana are presented, as well as the first photographic voucher of an adult specimen of D. latifrontalis, along with photographs of all known Ecuadorian Dipsadini species. The first country record of D. variegata in Ecuador is provided and D. oligozonata removed from the list of Peruvian herpetofauna. With these changes, the number of Dipsadini reported in Ecuador increases to 22, 18 species of Dipsas and four of Sibon.

Novel Approaches for Phylogenetic Inference from Morphological Data and Total-Evidence Dating in Squamate Reptiles (Lizards, Snakes, and Amphisbaenians)

Here, I combine previously underutilized models and priors to perform more biologically realistic phylogenetic inference from morphological data, with an example from squamate reptiles. When coding morphological characters, it is often possible to denote ordered states with explicit reference to observed or hypothetical ancestral conditions. Using this logic, we can integrate across character-state labels and estimate meaningful rates of forward and backward transitions from plesiomorphy to apomorphy. I refer to this approach as MkA, for “asymmetric.” The MkA model incorporates the biological reality of limited reversal for many phylogenetically informative characters, and significantly increases likelihoods in the empirical data sets. Despite this, the phylogeny of Squamata remains contentious. Total-evidence analyses using combined morphological and molecular data and the MkA approach tend toward recent consensus estimates supporting a nested Iguania. However, support for this topology is not unambiguous across data sets or analyses, and no mechanism has been proposed to explain the widespread incongruence between partitions, or the hidden support for various topologies in those partitions. Furthermore, different morphological data sets produced by different authors contain both different characters and different states for the same or similar characters, resulting in drastically different placements for many important fossil lineages. Effort is needed to standardize ontology for morphology, resolve incongruence, and estimate a robust phylogeny. The MkA approach provides a preliminary avenue for investigating morphological evolution while accounting for temporal evidence and asymmetry in character-state changes.

Patterns, biases and prospects in the distribution and diversity of Neotropical snakes

MOTIVATION

We generated a novel database of Neotropical snakes (one of the world's richest herpetofauna) combining the most comprehensive, manually compiled distribution dataset with publicly available data. We assess, for the first time, the diversity patterns for all Neotropical snakes as well as sampling density and sampling biases.

MAIN TYPES OF VARIABLES CONTAINED

We compiled three databases of species occurrences: a dataset downloaded from the Global Biodiversity Information Facility (GBIF), a verified dataset built through taxonomic work and specialized literature, and a combined dataset comprising a cleaned version of the GBIF dataset merged with the verified dataset.

SPATIAL LOCATION AND GRAIN

Neotropics, Behrmann projection equivalent to 1° × 1°.

TIME PERIOD

Specimens housed in museums during the last 150 years.

MAJOR TAXA STUDIED

Squamata: Serpentes.

SOFTWARE FORMAT

Geographical information system (GIS).

RESULTS

The combined dataset provides the most comprehensive distribution database for Neotropical snakes to date. It contains 147,515 records for 886 species across 12 families, representing 74% of all species of snakes, spanning 27 countries in the Americas. Species richness and phylogenetic diversity show overall similar patterns. Amazonia is the least sampled Neotropical region, whereas most well-sampled sites are located near large universities and scientific collections. We provide a list and updated maps of geographical distribution of all snake species surveyed.

MAIN CONCLUSIONS

The biodiversity metrics of Neotropical snakes reflect patterns previously documented for other vertebrates, suggesting that similar factors may determine the diversity of both ectothermic and endothermic animals. We suggest conservation strategies for high-diversity areas and sampling efforts be directed towards Amazonia and poorly known species.

Diversification of the rainfrog Pristimantis ornatissimus in the lowlands and Andean foothills of Ecuador

Geographic barriers and elevational gradients have long been recognized as important in species diversification. Here, we illustrate an example where both mechanisms have shaped the genetic structure of the Neotropical rainfrog, Pristimantis ornatissimus, which has also resulted in speciation. This species was thought to be a single evolutionary lineage distributed throughout the Ecuadorian Chocó and the adjacent foothills of the Andes. Based on recent sampling of P. ornatissimus sensu lato, we provide molecular and morphological evidence that support the validity of a new species, which we name Pristimantis ecuadorensis sp. nov. The sister species are elevational replacements of each other; the distribution of Pristimantis ornatissimus sensu stricto is limited to the Ecuadorian Chocó ecoregion (< 1100 m), whereas the new species has only been found at Andean localities between 1450–1480 m. Given the results of the Multiple Matrix Regression with Randomization analysis, the genetic difference between P. ecuadorensis and P. ornatissimus is not explained by geographic distance nor environment, although environmental variables at a finer scale need to be tested. Therefore this speciation event might be the byproduct of stochastic historic extinction of connected populations or biogeographic events caused by barriers to dispersal such as rivers. Within P. ornatissimus sensu stricto, morphological patterns and genetic structure seem to be related to geographic isolation (e.g., rivers). Finally, we provide an updated phylogeny for the genus, including the new species, as well as other Ecuadorian Pristimantis.

A revision and key for the tribe Diaphorolepidini (Serpentes: Dipsadidae) and checklist for the genus Synophis

The genus Synophis contains a number of enigmatic species, distributed primarily in the Andean highlands of northern South America. Their extreme crypsis and rarity has precluded detailed study of most species. A recent flurry of collection activity resulted in the accession of many new specimens, and the description of 4 new species in 2015, doubling the number of described taxa. However, lingering questions remain regarding the assignment of many new and historical specimens, the morphological limits and geographical ranges of the species, and their phylogenetic relationships. We analyze new and existing morphological and molecular data to produce a new molecular phylogeny and revised morphological descriptions. We validate the previously unavailable tribe name Diaphorolepidini Jenner, Pyron, Arteaga, Echevarría, & Torres-Carvajal tribe nov., describe a 9th species Synophis niceforomariae Pyron, Arteaga, Echevarría, & Torres-Carvajal sp. nov., and offer new Standard Names in English and Spanish for the group: Andean Shadow Snakes and Culebras Andinas de la Sombra, respectively. A variety of features such as vertebrae and hemipenes show an interesting range of variation in the group, which should be evaluated in future studies, to refine species limits and diagnoses further. Cryptic and undiscovered diversity undoubtedly remains, and we hope this summary provides a robust basis for future work.

Predicting community structure in snakes on Eastern Nearctic islands using ecological neutral theory and phylogenetic methods

Predicting species presence and richness on islands is important for understanding the origins of communities and how likely it is that species will disperse and resist extinction. The equilibrium theory of island biogeography (ETIB) and, as a simple model of sampling abundances, the unified neutral theory of biodiversity (UNTB), predict that in situations where mainland to island migration is high, species-abundance relationships explain the presence of taxa on islands. Thus, more abundant mainland species should have a higher probability of occurring on adjacent islands. In contrast to UNTB, if certain groups have traits that permit them to disperse to islands better than other taxa, then phylogeny may be more predictive of which taxa will occur on islands. Taking surveys of 54 island snake communities in the Eastern Nearctic along with mainland communities that have abundance data for each species, we use phylogenetic assembly methods and UNTB estimates to predict island communities. Species richness is predicted by island area, whereas turnover from the mainland to island communities is random with respect to phylogeny. Community structure appears to be ecologically neutral and abundance on the mainland is the best predictor of presence on islands. With regard to young and proximate islands, where allopatric or cladogenetic speciation is not a factor, we find that simple neutral models following UNTB and ETIB predict the structure of island communities.

Methodological congruence in phylogenomic analyses with morphological support for teiid lizards (Sauria: Teiidae)

A well-known issue in phylogenetics is discordance among gene trees, species trees, morphology, and other data types. Gene-tree discordance is often caused by incomplete lineage sorting, lateral gene transfer, and gene duplication. Multispecies-coalescent methods can account for incomplete lineage sorting and are believed by many to be more accurate than concatenation. However, simulation studies and empirical data have demonstrated that concatenation and species tree methods often recover similar topologies. We use three popular methods of phylogenetic reconstruction (one concatenation, two species tree) to evaluate relationships within Teiidae. These lizards are distributed across the United States to Argentina and the West Indies, and their classification has been controversial due to incomplete sampling and the discordance among various character types (chromosomes, DNA, musculature, osteology, etc.) used to reconstruct phylogenetic relationships. Recent morphological and molecular analyses of the group resurrected three genera and created five new genera to resolve non-monophyly in three historically ill-defined genera: Ameiva, Cnemidophorus, and Tupinambis. Here, we assess the phylogenetic relationships of the Teiidae using "next-generation" anchored-phylogenomics sequencing. Our final alignment includes 316 loci (488,656bp DNA) for 244 individuals (56 species of teiids, representing all currently recognized genera) and all three methods (ExaML, MP-EST, and ASTRAL-II) recovered essentially identical topologies. Our results are basically in agreement with recent results from morphology and smaller molecular datasets, showing support for monophyly of the eight new genera. Interestingly, even with hundreds of loci, the relationships among some genera in Tupinambinae remain ambiguous (i.e. low nodal support for the position of Salvator and Dracaena).

Comparative Phylogeography Reveals Cryptic Diversity and Repeated Patterns of Cladogenesis for Amphibians and Reptiles in Northwestern Ecuador

Comparative phylogeography allow us to understand how shared historical circumstances have shaped the formation of lineages, by examining a broad spectrum of co-distributed populations of different taxa. However, these types of studies are scarce in the Neotropics, a region that is characterized by high diversity, complex geology, and poorly understood biogeography. Here, we investigate the diversification patterns of five lineages of amphibians and reptiles, co-distributed across the Choco and Andes ecoregions in northwestern Ecuador. Mitochondrial DNA and occurrence records were used to determine the degree of geographic genetic divergence within species. Our results highlight congruent patterns of parapatric speciation and common geographical barriers for distantly related taxa. These comparisons indicate similar biological and demographic characteristics for the included clades, and reveal the existence of two new species of Pristimantis previously subsumed under P. walkeri, which we describe herein. Our data supports the hypothesis that widely distributed Chocoan taxa may generally experience their greatest opportunities for isolation and parapatric speciation across thermal elevational gradients. Finally, our study provides critical information to predict which unstudied lineages may harbor cryptic diversity, and how geology and climate are likely to have shaped their evolutionary history.

Biogeography of worm lizards (Amphisbaenia) driven by end-Cretaceous mass extinction

Worm lizards (Amphisbaenia) are burrowing squamates that live as subterranean predators. Their underground existence should limit dispersal, yet they are widespread throughout the Americas, Europe and Africa. This pattern was traditionally explained by continental drift, but molecular clocks suggest a Cenozoic diversification, long after the break-up of Pangaea, implying dispersal. Here, we describe primitive amphisbaenians from the North American Palaeocene, including the oldest known amphisbaenian, and provide new and older molecular divergence estimates for the clade, showing that worm lizards originated in North America, then radiated and dispersed in the Palaeogene following the Cretaceous-Palaeogene (K-Pg) extinction. This scenario implies at least three trans-oceanic dispersals: from North America to Europe, from North America to Africa and from Africa to South America. Amphisbaenians provide a striking case study in biogeography, suggesting that the role of continental drift in biogeography may be overstated. Instead, these patterns support Darwin and Wallace's hypothesis that the geographical ranges of modern clades result from dispersal, including oceanic rafting. Mass extinctions may facilitate dispersal events by eliminating competitors and predators that would otherwise hinder establishment of dispersing populations, removing biotic barriers to dispersal.

Systematics of Nothopsini (Serpentes, Dipsadidae), with a new species of Synophis from the Pacific Andean slopes of southwestern Ecuador

Within Dipsadinae, some recent authors have recognized a tribe Nothopsini containing the genera Diaphorolepis, Emmochliophis, Nothopsis, Synophis, and Xenopholis, on the basis of a number of putative morphological synapomorphies. However, molecular results suggest that Nothopsis, Synophis, and Xenopholis do not form a monophyletic group, while the remaining taxa are unsampled in recent molecular phylogenies. Here, DNA-sequence data for some Diaphorolepis and Synophis species are provided for the first time, as well as additional new sequences for Nothopsis and some Synophis species. Including these and other existing data for nothopsine species, previous studies showing that Nothopsini is not a natural group are corroborated. Nothopsini Cope, 1871 is restricted to Nothopsis. Diaphorolepidini Jenner, 1981 is resurrected and re-delimited to include only Diaphorolepis, Emmochliophis, and Synophis. Finally, Xenopholis remains Dipsadinae incertae sedis. Known material of Diaphorolepidini is reviewed to generate revised and expanded descriptions and diagnoses at the tribe, genus, and species level. Numerous cryptic species are likely present in Synophis bicolor and Synophis lasallei. Finally, a new population from the low-elevation cloud forests of SW Ecuador is reported upon, which is genetically and morphologically distinct from all other species, that is here named Synophis zaheri sp. n.

Two novel genera and one new species of treefrog (Anura: Rhacophoridae) highlight cryptic diversity in the Western Ghats of India

Amphibian diversity in the Western Ghats-Sri Lanka biodiversity hotspot is extremely high, especially for such a geographically restricted area. Frogs in particular dominate these assemblages, and the family Rhacophoridae is chief among these, with hundreds of endemic species. These taxa continue to be described at a rapid pace, and several groups have recently been found to represent unique evolutionary clades at the genus level. Here, we report DNA sequences, larval and breeding data for two species of rhacophorid treefrog (Polypedates bijui and a new, hitherto undescribed species). Remarkably, they represent unique, independent clades which form successive sister groups to the Pseudophilautus (Sri Lanka) + Raorchestes (India, China & Indochina) clades. We place these species into two new genera (Beddomixalus gen. nov. and Mercurana gen. nov.). Both of these genera exhibit a distinct reproductive mode among Rhacophoridae of peninsular India and Sri Lanka, with explosive breeding and semiterrestrial, unprotected, non-pigmented eggs oviposited in seasonal swamp pools, which hatch into exotrophic, free-living aquatic tadpoles. Relationships and representation of reproductive modes in sister taxa within the larger clade into which these novel genera are placed, is also discussed. These results suggest that more undescribed taxa may remain to be discovered in South Asia, and the crucial importance of conserving remaining viable habitats.

An evaluation of fossil tip-dating versus node-age calibrations in tetraodontiform fishes (Teleostei: Percomorphaceae)

Time-calibrated phylogenies based on molecular data provide a framework for comparative studies. Calibration methods to combine fossil information with molecular phylogenies are, however, under active development, often generating disagreement about the best way to incorporate paleontological data into these analyses. This study provides an empirical comparison of the most widely used approach based on node-dating priors for relaxed clocks implemented in the programs BEAST and MrBayes, with two recently proposed improvements: one using a new fossilized birth-death process model for node dating (implemented in the program DPPDiv), and the other using a total-evidence or tip-dating method (implemented in MrBayes and BEAST). These methods are applied herein to tetraodontiform fishes, a diverse group of living and extinct taxa that features one of the most extensive fossil records among teleosts. Previous estimates of time-calibrated phylogenies of tetraodontiforms using node-dating methods reported disparate estimates for their age of origin, ranging from the late Jurassic to the early Paleocene (ca. 150-59Ma). We analyzed a comprehensive dataset with 16 loci and 210 morphological characters, including 131 taxa (95 extant and 36 fossil species) representing all families of fossil and extant tetraodontiforms, under different molecular clock calibration approaches. Results from node-dating methods produced consistently younger ages than the tip-dating approaches. The older ages inferred by tip dating imply an unlikely early-late Jurassic (ca. 185-119Ma) origin for this order and the existence of extended ghost lineages in their fossil record. Node-based methods, by contrast, produce time estimates that are more consistent with the stratigraphic record, suggesting a late Cretaceous (ca. 86-96Ma) origin. We show that the precision of clade age estimates using tip dating increases with the number of fossils analyzed and with the proximity of fossil taxa to the node under assessment. This study suggests that current implementations of tip dating may overestimate ages of divergence in calibrated phylogenies. It also provides a comprehensive phylogenetic framework for tetraodontiform systematics and future comparative studies.

A taxonomic revision of the Asian keelback snakes, genus Amphiesma (Serpentes: Colubridae: Natricinae), with description of a new species

The Asian keelback snakes (genus Amphiesma) are a widely distributed group of Old World natricines, inhabiting a variety of niches and exhibiting significant morphological variation. Recent molecular phylogenies suggest that this genus is not monophyletic, and that additional cryptic diversity is also likely present. We conducted a phylogenetic analysis of the group based on 3162 bp of one mitochondrial gene (Cyt. b) and three nuclear genes (C-mos, Rag1, NT3), sampling 18 species in addition to those sequenced in previous works. All analyses consistently show that Amphiesma consists of three distinct, monophyletic lineages with strong support. We divide Amphiesma into three genera, Amphiesma, Hebius, and Herpetoreas. The genus Amphiesma is monotypic, Herpetoreas contains three species, and Hebius comprises the remaining 39 species. On the basis of a combination of molecular analyses and external morphological comparisons, we describe a new species in the Herpetoreas group from China as H. burbrinki sp. nov. Several other species are shown to be non-monophyletic or contain significant levels of intraspecific genetic diversity. Another Old World natricine genera, Xenochrophis is also found to be non-monophyletic. Our results indicate that further taxonomic revisions are needed in Natricinae, at multiple levels.