Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae)

Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.

Specimen collection is essential for modern science

Natural history museums are vital repositories of specimens, samples and data that inform about the natural world; this Formal Comment revisits a Perspective that advocated for the adoption of compassionate collection practices, querying whether it will ever be possible to completely do away with whole animal specimen collection.

Considering admixture when producing draft genomes: an example in North American ratsnakes (Pantherophis alleghaniensis/Pantherophis obsoletus)

The number of reference genomes of snakes lags behind several other vertebrate groups (e.g. birds and mammals). However, in the last two years, a concerted effort by researchers from around the world has produced new genomes of snakes representing members from several new families. Here, we present a high-quality, annotated genome of the central ratsnake (Pantherophis alleghaniensis), a member of the most diverse snake lineage, Colubroidea. Pantherophis alleghaniensis is found in the central part of the Nearctic, east of the Mississippi River. This genome was sequenced using 10X Chromium synthetic long reads and polished using Illumina short reads. The final genome assembly had an N50 of 21.82 Mb and an L50 of 22 scaffolds with a maximum scaffold length of 82.078 Mb. The genome is composed of 49.24% repeat elements dominated by long interspersed elements. We annotated this genome using transcriptome assemblies from 14 tissue types and recovered 28,368 predicted proteins. Finally, we estimated admixture proportions between two species of ratsnakes and discovered that this specimen is an admixed individual containing genomes from the western (Pantherophis obsoletus) and central ratsnakes (P. alleghaniensis). We discuss the importance of considering interspecific admixture in downstream approaches for inferring demography and phylogeny.

Speciation rates are unrelated to the formation of population structure in Malagasy gemsnakes

Speciation rates vary substantially across the tree of life. These rates should be linked to the rate at which population structure forms if a continuum between micro and macroevolutionary patterns exists. Previous studies examining the link between speciation rates and the degree of population formation in clades have been shown to be either correlated or uncorrelated depending on the group, but no study has yet examined the relationship between speciation rates and population structure in a young group that is constrained spatially to a single-island system. We examine this correlation in 109 gemsnakes (Pseudoxyrhophiidae) endemic to Madagascar and originating in the early Miocene, which helps control for extinction variation across time and space. We find no relationship between rates of speciation and the formation rates of population structure over space in 33 species of gemsnakes. Rates of speciation show low variation, yet population structure varies widely across species, indicating that speciation rates and population structure are disconnected. We suspect this is largely due to the persistence of some lineages not susceptible to extinction. Importantly, we discuss how delimiting populations versus species may contribute to problems understanding the continuum between shallow and deep evolutionary processes.

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

Ultraconserved elements-based phylogenomic systematics of the snake superfamily Elapoidea, with the description of a new Afro-Asian family

The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov., for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes.

Widespread Reticulate Evolution in an Adaptive Radiation

A fundamental assumption of evolutionary biology is that phylogeny follows a bifurcating process. However, hybrid speciation and introgression are becoming more widely documented in many groups. Hybrid inference studies have been historically limited to small sets of taxa, while exploration of the prevalence and trends of reticulation at deep time scales remains unexplored. We study the evolutionary history of an adaptive radiation of 109 gemsnakes in Madagascar (Pseudoxyrhophiinae) to identify potential instances of introgression. Using several network inference methods, we find twelve reticulation events within the 22-million-year evolutionary history of gemsnakes, producing 28% of the diversity for the group, including one reticulation that resulted in the diversification of an 18 species radiation. These reticulations occur at nodes with high gene tree discordance. Hybridization events occurred between north-south distributed parentals which share similar ecologies. Younger hybrids occupy intermediate contact zones between the parentals, showing that post-speciation dispersal in this group has not eroded the spatial signatures of introgression. Reticulations accumulated consistently over time, despite drops in overall speciation rates during the Pleistocene. This suggests that while bifurcating speciation may decline as the result of species accumulation and environmental change, speciation by hybridization may be more robust to these processes.

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.

Empirical and philosophical problems with the subspecies rank

Species‐level taxonomy derives from empirical sources (data and techniques) that assess the existence of spatiotemporal evolutionary lineages via various species “concepts.” These concepts determine if observed lineages are independent given a particular methodology and ontology, which relates the metaphysical species concept to what “kind” of thing a species is in reality. Often, species concepts fail to link epistemology back to ontology. This lack of coherence is in part responsible for the persistence of the subspecies rank, which in modern usage often functions as a placeholder between the evolutionary events of divergence or collapse of incipient species. Thus, prospective events like lineages merging or diverging require information from unknowable future information. This is also conditioned on evidence that the lineage already has a detectably distinct evolutionary history. Ranking these lineages as subspecies can seem attractive given that many lineages do not exhibit intrinsic reproductive isolation. We argue that using subspecies is indefensible on philosophical and empirical grounds. Ontologically, the rank of subspecies is either identical to that of species or undefined in the context of evolutionary lineages representing spatiotemporally defined individuals. Some species concepts more inclined to consider subspecies, like the Biological Species Concept, are disconnected from evolutionary ontology and do not consider genealogy. Even if ontology is ignored, methods addressing reproductive isolation are often indirect and fail to capture the range of scenarios linking gene flow to species identity over space and time. The use of subspecies and reliance on reproductive isolation as a basis for an operational species concept can also conflict with ethical issues governing the protection of species. We provide a way forward for recognizing and naming species that links theoretical and operational species concepts regardless of the magnitude of reproductive isolation.

Complex genetic patterns and distribution limits mediated by native congeners of the worldwide invasive red-eared slider turtle

Non-native (invasive) species offer a unique opportunity to study the geographic distribution and range limits of species, wherein the evolutionary change driven by interspecific interactions between native and non-native closely related species is a key component. The red-eared slider turtle, Trachemys scripta elegans (TSE), has been introduced and successfully established worldwide. It can coexist with its native congeners T. cataspila, T. venusta and T. taylori in Mexico. We performed comprehensive fieldwork, executed a battery of genetic analyses and applied a novel species distribution modeling approach to evaluate their historical lineage relationships and contemporary population genetic patterns. Our findings support the historical common ancestry between native TSE and non-native (TSE_alien ), while also highlighting the genetic differentiation of the exotic lineage. Genetic patterns are associated with their range size/endemism gradient, the microendemic T. taylori showed significant reduced genetic diversity and high differentiation, whereas TSE_alien showed the highest diversity and signals of population size expansion. Counter to our expectations, lower naturally occurring distribution overlap and little admixture patterns were found between TSE and its congeners, exhibiting reduced gene flow and clear genetic separation across neighboring species despite having zones of contact. We demonstrate that these native Trachemys species have distinct climatic niche suitability, likely preventing establishment of and displacement by the TSE_alien . Moreover, we found major niche overlap between TSE_alien and native species worldwide, supporting our prediction that sites with closer ecological optima to the invasive species have higher establishment risk than those that are closer to the niche-center of the native species.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

Morphological description of a new specimen of Herpetoreas burbrinki Guo et al 2014 (Serpentes: Colubridae)

The original description of Burbrinks Keelback, Herpetoreas burbrinki was based on a sole damaged specimen collected from Zayu County, Xizang Autonomous Region, China in September 2007. On 16 August 2019, we collected a second live adult female specimen from the type locality. The identity of the species is established based on morphological and molecular comparison with the holotype. One mitochondrial gene (Cytb) and three nuclear genes (C-mos, Rag1, NT3) of the new specimen were sequenced. The four sequences all share the same haplotypes with the holotype. We describe the coloration in life, variation with the type and expand the morphological description of this species.

A new snake species of the genus Gonyosoma Wagler, 1828 (Serpentes: Colubridae) from Hainan Island, China

A new species of the genus Gonyosoma Wagler, 1828 is described herein based on six specimens from the Diaoluoshan Mountains, Hainan Island, Hainan Province, China. The new species, Gonyosoma hainanensesp. nov., is most similar to its continental sister species, Gonyosoma boulengeri (Mocquard, 1897). Both taxa have a scaled protrusion on the anterior portion of the rostrum, distinct from other congeners. However, Gonyosoma hainanensesp. nov. can be distinguished from G. boulengeri by two significant morphological characters: (1) black orbital stripe absent in adults (vs. present in G. boulengeri); and (2) two loreals (vs. one loreal in G. boulengeri). The new species is also genetically divergent and forms a unique clade from its sister species and all other congeners based on sequences of the mitochondrial gene cytochrome b (cyt b).

Quaternary climatic fluctuations influence the demographic history of two species of sky-island endemic amphibians in the Neotropics

We evaluated the role of Quaternary climatic fluctuations on the demographic history and population structure of amphibian species endemic to the 'campo rupestre' in the Neotropics, evaluating their distributional shifts, demographic changes, and lineage formation from the end of Pleistocene to present. We chose two anurans endemic to the high-elevation 'campo rupestre' in the Espinhaço Range (ER) in northeastern and southeastern Brazil (Bokermannohyla alvarengai and Bokermannohyla oxente), as models to test the role of Quaternary climatic fluctuations over their distribution range in this region. We collected tissue samples throughout their distribution range and used statistical phylogeography to examine processes of divergence and population demography. We generated spatial-temporal reconstructions using Bayesian inference in a coalescent framework in combination with hind-cast projections of species distribution models (SDMs). We also used the results and literature information to test alternative diversification scenarios via approximate Bayesian computation (ABC). Our results show that Quaternary climatic fluctuations influenced the geographic ranges of both species showing population expansion during the last glacial maximum (LGM) and range contraction during interglacial periods, as inferred from selected ABC models and from past projections of SDMs. We recovered Pleistocene diversification for both species occuring in distinctly unique periods for each taxon. An older and range-restricted lineage was recovered in a geographically isolated geological massif, deserving conservation and further taxonomic study. The diversification and distribution of these amphibian species endemic to the Neotropical 'campo rupestre' were influenced by Quaternary climatic fluctuations. The expansion of cold adapted species restricted to higher elevations during glacial periods and their concomitant retraction during interglacial periods may have been crucial for producing patterns of species richness and endemism along elevation gradients in tropical and subtropical domains. Such processes may influence the evolution of the biota distributed in heterogeneous landscapes with varied topography.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

Phylogenetic relationships and biogeographic range evolution in cat-eyed snakes, Boiga (Serpentes: Colubridae)

The genus Boiga includes 35, primarily arboreal snake species distributed from the Middle East to Australia and many islands in the western Pacific, with particularly high species diversity in South-East Asia. Despite including the iconic mangrove snakes (Boiga dendrophila complex) and the brown tree snake (Boiga irregularis; infamous for avian extinctions on small islands of the Pacific), species-level phylogenetic relationships and the biogeographic history of this ecologically and morphologically distinct clade are poorly understood. In this study, we sequenced mitochondrial and nuclear DNA for 24 Boiga species and used these data to estimate a robust phylogenetic inference, in order to (1) test the hypothesis that Boiga is monophyletic, (2) evaluate the validity of current species-level taxonomy and (3) examine whether geographic range evolution in Boiga is consistent with expectations concerning dispersal and colonization of vertebrates between continents and islands. Our results support the prevailing view that most dispersal events are downstream – from continents to oceanic islands – but we also identify a role for upstream dispersal from oceanic islands to continents. Additionally, the novel phylogeny of Boiga presented here is informative for updating species-level taxonomy within the genus.

Biogeographic barriers, Pleistocene refugia, and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species-specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.

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).

Female-biased gape and body-size dimorphism in the New World watersnakes (tribe: Thamnophiini) oppose predictions from Rensch's rule

ABSTRACT

Sexual-size dimorphism (SSD) is ubiquitous across animals and often biased in the direction of larger females in snakes and other ectothermic organisms. To understand how SSD evolves across species, Rensch's rule predicts that in taxa where males are larger, SSD increases with body size. In contrast, where females are larger, SSD decreases with body size. While this rule holds for many taxa, it may be ambiguous for others, particularly ectothermic vertebrates. Importantly, this rule suggests that the outcomes of SSD over phylogenetic time scales depend on the direction of dimorphism predicated on the difference in reproductive efforts between males and females. Here, we examine SSD in the context of Rensch's rule in Thamnophiini, the gartersnakes and watersnakes, a prominent group that in many areas comprises the majority of the North American snake biota. Using a dated phylogeny, measurements of gape, body, and tail size, we show that these snakes do not follow Rensch's rule, but rather female-biased SSD increases with body size. We in turn find that this allometry is most pronounced with gape and is correlated with both neonate and litter size, suggesting that acquiring prey of increased size may be directly related to fecundity selection. These changes in SSD are not constrained to any particular clade; we find no evidence of phylogenetic shifts in those traits showing SSD. We suggest several ways forward to better understand the anatomical units of selection for SSD and modularity.

OPEN RESEARCH BADGES

This article has been awarded Open Data and Open Materials Badges. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.5061/dryad.3pn57h0.

The Origins and Diversification of the Exceptionally Rich Gemsnakes (Colubroidea: Lamprophiidae: Pseudoxyrhophiinae) in Madagascar

Processes leading to spectacular diversity of both form and species on islands have been well-documented under island biogeography theory, where distance from source and island size are key factors determining immigration and extinction resistance. But far less understood are the processes governing in situ diversification on the world’s mega islands, where large and isolated land masses produced morphologically distinct radiations from related taxa on continental regions. Madagascar has long been recognized as a natural laboratory due to its isolation, lack of influence from adjacent continents, and diversification of spectacular vertebrate radiations. However, only a handful of studies have examined rate shifts of in situ diversification for this island. Here, we examine rates of diversification in the Malagasy snakes of the family Pseudoxyrhophiinae (gemsnakes) to understand if rates of speciation were initially high, enhanced by diversification into distinct biomes, and associated with key dentition traits. Using a genomic sequence-capture data set for 366 samples, we determine that all previously described and newly discovered species are delimitable and therefore useful candidates for understanding diversification trajectories through time. Our analysis detected no shifts in diversification rate between clades or changes in biome or dentition type. Remarkably, we demonstrate that rates of diversification of the gemsnake radiation, which originated in Madagascar during the early Miocene, remained steady throughout the Neogene. However, we do detect a significant slowdown in diversification during the Pleistocene. We also comment on the apparent paradox where most living species originated in the Pleistocene, despite diversification rates being substantially higher during the earlier 15 myr.

Hibernation in bats (Mammalia: Chiroptera) did not evolve through positive selection of leptin

Temperature regulation is an indispensable physiological activity critical for animal survival. However, relatively little is known about the origin of thermoregulatory regimes in a phylogenetic context, or the genetic mechanisms driving the evolution of these regimes. Using bats as a study system, we examined the evolution of three thermoregulatory regimes (hibernation, daily heterothermy, and homeothermy) in relation to the evolution of leptin, a protein implicated in regulation of torpor bouts in mammals, including bats. A threshold model was used to test for a correlation between lineages with positively selected lep, the gene encoding leptin, and the thermoregulatory regimes of those lineages. Although evidence for episodic positive selection of lep was found, positive selection was not correlated with lineages of heterothermic bats, a finding that contradicts results from previous studies. Evidence from our ancestral state reconstructions suggests that the most recent common ancestor of bats used daily heterothermy and that the presence of hibernation is highly unlikely at this node. Hibernation likely evolved independently at least four times in bats-once in the common ancestor of Vespertilionidae and Molossidae, once in the clade containing Rhinolophidae and Rhinopomatidae, and again independently in the lineages leading to Taphozous melanopogon and Mystacina tuberculata. Our reconstructions revealed that thermoregulatory regimes never transitioned directly from hibernation to homeothermy, or the reverse, in the evolutionary history of bats. This, in addition to recent evidence that heterothermy is best described along a continuum, suggests that thermoregulatory regimes in mammals are best represented as an ordered continuous trait (homeothermy ← → daily torpor ← → hibernation) rather than as the three discrete regimes that evolve in an unordered fashion. These results have important implications for methodological approaches in future physiological and evolutionary research.

The Biogeography of Deep Time Phylogenetic Reticulation

Most phylogenies are typically represented as purely bifurcating. However, as genomic data have become more common in phylogenetic studies, it is not unusual to find reticulation among terminal lineages or among internal nodes (deep time reticulation; DTR). In these situations, gene flow must have happened in the same or adjacent geographic areas for these DTRs to have occurred and therefore biogeographic reconstruction should provide similar area estimates for parental nodes, provided extinction or dispersal has not eroded these patterns. We examine the phylogeny of the widely distributed New World kingsnakes (Lampropeltis), determine if DTR is present in this group, and estimate the ancestral area for reticulation. Importantly, we develop a new method that uses coalescent simulations in a machine learning framework to show conclusively that this phylogeny is best represented as reticulating at deeper time. Using joint probabilities of ancestral area reconstructions on the bifurcating parental lineages from the reticulating node, we show that this reticulation likely occurred in northwestern Mexico/southwestern US, and subsequently, led to the diversification of the Mexican kingsnakes. This region has been previously identified as an area important for understanding speciation and secondary contact with gene flow in snakes and other squamates. This research shows that phylogenetic reticulation is common, even in well-studied groups, and that the geographic scope of ancient hybridization is recoverable.

Host susceptibility to snake fungal disease is highly dispersed across phylogenetic and functional trait space

Emerging infectious diseases (EIDs) reduce host population sizes, cause extinction, disassemble communities, and have indirect negative effects on human well-being. Fungal EIDs have reduced population abundances in amphibians and bats across many species over large areas. The recent emergence of snake fungal disease (SFD) may have caused declines in some snake populations in the Eastern United States (EUS), which is home to a phylogenetically and ecologically diverse assembly of 98 taxa. SFD has been documented in only 23 naturally occuring species, although this is likely an underestimate of the number of susceptible taxa. Using several novel methods, including artificial neural networks, we combine phylogenetic and trait-based community estimates from all taxa in this region to show that SFD hosts are both phylogenetically and ecologically randomly dispersed. This might indicate that other species of snakes in the EUS could be currently infected or susceptible to SFD. Our models also indicate that information about key traits that enhance susceptiblity is lacking. Surveillance should consider that all snake species and habitats likely harbor this pathogen.

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.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

Comparing species tree estimation with large anchored phylogenomic and small Sanger-sequenced molecular datasets: an empirical study on Malagasy pseudoxyrhophiine snakes

BACKGROUND

Using molecular data generated by high throughput next generation sequencing (NGS) platforms to infer phylogeny is becoming common as costs go down and the ability to capture loci from across the genome goes up. While there is a general consensus that greater numbers of independent loci should result in more robust phylogenetic estimates, few studies have compared phylogenies resulting from smaller datasets for commonly used genetic markers with the large datasets captured using NGS. Here, we determine how a 5-locus Sanger dataset compares with a 377-locus anchored genomics dataset for understanding the evolutionary history of the pseudoxyrhophiine snake radiation centered in Madagascar. The Pseudoxyrhophiinae comprise ~86 % of Madagascar's serpent diversity, yet they are poorly known with respect to ecology, behavior, and systematics. Using the 377-locus NGS dataset and the summary statistics species-tree methods STAR and MP-EST, we estimated a well-supported species tree that provides new insights concerning intergeneric relationships for the pseudoxyrhophiines. We also compared how these and other methods performed with respect to estimating tree topology using datasets with varying numbers of loci.

METHODS

Using Sanger sequencing and an anchored phylogenomics approach, we sequenced datasets comprised of 5 and 377 loci, respectively, for 23 pseudoxyrhophiine taxa. For each dataset, we estimated phylogenies using both gene-tree (concatenation) and species-tree (STAR, MP-EST) approaches. We determined the similarity of resulting tree topologies from the different datasets using Robinson-Foulds distances. In addition, we examined how subsets of these data performed compared to the complete Sanger and anchored datasets for phylogenetic accuracy using the same tree inference methodologies, as well as the program *BEAST to determine if a full coalescent model for species tree estimation could generate robust results with fewer loci compared to the summary statistics species tree approaches. We also examined the individual gene trees in comparison to the 377-locus species tree using the program MetaTree.

RESULTS

Using the full anchored dataset under a variety of methods gave us the same, well-supported phylogeny for pseudoxyrhophiines. The African pseudoxyrhophiine Duberria is the sister taxon to the Malagasy pseudoxyrhophiines genera, providing evidence for a monophyletic radiation in Madagascar. In addition, within Madagascar, the two major clades inferred correspond largely to the aglyphous and opisthoglyphous genera, suggesting that feeding specializations associated with tooth venom delivery may have played a major role in the early diversification of this radiation. The comparison of tree topologies from the concatenated and species-tree methods using different datasets indicated the 5-locus dataset cannot beused to infer a correct phylogeny for the pseudoxyrhophiines under any method tested here and that summary statistics methods require 50 or more loci to consistently recover the species-tree inferred using the complete anchored dataset. However, as few as 15 loci may infer the correct topology when using the full coalescent species tree method *BEAST. MetaTree analyses of each gene tree from the Sanger and anchored datasets found that none of the individual gene trees matched the 377-locus species tree, and that no gene trees were identical with respect to topology.

CONCLUSIONS

Our results suggest that ≥50 loci may be necessary to confidently infer phylogenies when using summaryspecies-tree methods, but that the coalescent-based method *BEAST consistently recovers the same topology using only 15 loci. These results reinforce that datasets with small numbers of markers may result in misleading topologies, and further, that the method of inference used to generate a phylogeny also has a major influence on the number of loci necessary to infer robust species trees.

Independent Demographic Responses to Climate Change among Temperate and Tropical Milksnakes (Colubridae: Genus Lampropeltis)

The effects of Late Quaternary climate change have been examined for many temperate New World taxa, but the impact of Pleistocene glacial cycles on Neotropical taxa is less well understood, specifically with respect to changes in population demography. Here, we examine historical demographic trends for six species of milksnake with representatives in both the temperate and tropical Americas to determine if species share responses to climate change as a taxon or by area (i.e., temperate versus tropical environments). Using a multilocus dataset, we test for the demographic signature of population expansion and decline using non-genealogical summary statistics, as well as coalescent-based methods. In addition, we determine whether range sizes are correlated with effective population sizes for milksnakes. Results indicate that there are no identifiable trends with respect to demographic response based on location, and that species responded to changing climates independently, with tropical taxa showing greater instability. There is also no correlation between range size and effective population size, with the largest population size belonging to the species with the smallest geographic distribution. Our study highlights the importance of not generalizing the demographic histories of taxa by region and further illustrates that the New World tropics may not have been a stable refuge during the Pleistocene.

The phylogenetic position and taxonomic status of the Rainbow Tree Snake Gonyophis margaritatus (Peters, 1871) (Squamata: Colubridae)

Molecular phylogenies have provided strong evidence for clarifying the taxonomy of groups with ambiguous morphological traits, thus avoiding potentially misleading conclusions based on evolutionary convergence of these traits. For snakes, established molecular databases along with new sequences from rare species allows us to estimate phylogenies, to clarify the phylogenetic relationships and test the monophyly of most taxonomic groups. Using one mitochondrial gene and five nuclear loci, we evaluate the taxonomic status of a rare Southeast Asian serpent, the Rainbow Tree Snake Gonyophis margaritatus (Squamata: Colubridae) by inferring a molecular phylogeny of 101 snake species. Both maximum likelihood and time- calibrated Bayesian inference phylogenies demonstrate that G. margaritatus is sister to Rhadinophis prasinus, previously considered to be part of a radiation of Old World ratsnakes. This group is in turn sister to a group containing Rhadinophis frenatus and Rhynchophis boulengeri with the entire clade originating in the mid-Miocene (~16 Ma) in Southeast Asia. This group is sister to the genus Gonyosoma and together originated in the early Miocene (~20 Ma). We discuss three potential solutions towards eliminating polyphyly of the genus Rhadinophis, but recommend using the genus name Gonyosoma for all species within this clade, which currently contains all of the species within the genera Gonyosoma, Gonyophis, Rhadinophis, and Rhynchophis. 

A hybrid phylogenetic-phylogenomic approach for species tree estimation in African Agama lizards with applications to biogeography, character evolution, and diversification

Africa is renowned for its biodiversity and endemicity, yet little is known about the factors shaping them across the continent. African Agama lizards (45 species) have a pan-continental distribution, making them an ideal model for investigating biogeography. Many species have evolved conspicuous sexually dimorphic traits, including extravagant breeding coloration in adult males, large adult male body sizes, and variability in social systems among colorful versus drab species. We present a comprehensive time-calibrated species tree for Agama, and their close relatives, using a hybrid phylogenetic-phylogenomic approach that combines traditional Sanger sequence data from five loci for 57 species (146 samples) with anchored phylogenomic data from 215 nuclear genes for 23 species. The Sanger data are analyzed using coalescent-based species tree inference using (*)BEAST, and the resulting posterior distribution of species trees is attenuated using the phylogenomic tree as a backbone constraint. The result is a time-calibrated species tree for Agama that includes 95% of all species, multiple samples for most species, strong support for the major clades, and strong support for most of the initial divergence events. Diversification within Agama began approximately 23 million years ago (Ma), and separate radiations in Southern, East, West, and Northern Africa have been diversifying for >10Myr. A suite of traits (morphological, coloration, and sociality) are tightly correlated and show a strong signal of high morphological disparity within clades, whereby the subsequent evolution of convergent phenotypes has accompanied diversification into new biogeographic areas.

Phylogenetic estimates of speciation and extinction rates for testing ecological and evolutionary hypotheses

Phylogenies are used to estimate rates of speciation and extinction, reconstruct historical diversification scenarios, and link these to ecological and evolutionary factors, such as climate or organismal traits. Recent models can now estimate the effects of binary, multistate, continuous, and biogeographic characters on diversification rates. Others test for diversity dependence (DD) in speciation and extinction, which has become recognized as an important process in numerous clades. A third class incorporates flexible time-dependent functions, enabling reconstruction of major periods of both expanding and contracting diversity. Although there are some potential problems (particularly for estimating extinction), these methods hold promise for answering many classic questions in ecology and evolution, such as the origin of adaptive radiations, and the latitudinal gradient in species richness.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.