Search

Aim We quantify the degree to which co-distributed snakes across the Cochise Filter Barrier ( CFB) have a shared history of population divergence and estimate the timing of divergence for each taxon pair. Location North America. Methods A single locus dataset was collected ( n = 747 individuals) for 12 snake taxon pairs. Phylogeographical structure was estimated for each taxon. Redundancy analyses were used to assess the importance of geographical distance, climate and putative barriers to gene flow in structuring genetic diversity. Hierarchical approximate Bayesian computation was used to estimate the magnitude of synchronicity in divergence times across a well-documented phylogeographical barrier. Lastly, gene divergence and population divergence times were estimated using multiple methods. Results There is substantial phylogeographical structure in many of the snake taxa, particularly at the CFB. A model containing distance, climate and barriers explained the greatest amount of genetic variation in nearly all taxa. When each variable was examined separately, climate explained the most variation. The hABC model testing indicates that there is overwhelming support for asynchronous phylogeographical histories within these co-distributed taxa. Estimated divergence times range throughout the Quaternary and Neogene. Main Conclusions We demonstrate that the 12 snake taxon pairs studied here have diversified within the desert Southwest forming distinct Sonoran and Chihuahuan populations, illustrating the importance of this region in driving diversification in North American taxa. Although these groups exhibit the same pattern of lineage formation, there is strong support for asynchronous diversification and little concordance in divergence time estimates. [ABSTRACT FROM AUTHOR]

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 data sets. 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 2 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 data sets 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. [ Desmognathus ; integrative taxonomy; machine learning; species delimitation.] [ABSTRACT FROM AUTHOR]

Aim: A long history of isolation coupled with complex landscapes makes Madagascar ideal for exploring the historical factors that have shaped patterns of population diversity and endemism. Many species‐level studies have suggested Late Quaternary climate change may have influenced population dynamics in the tropics, but unique biomes and individual species properties may have driven idiosyncratic responses to these shifts. Here, we use community‐scale population genetic data to investigate the impact of Quaternary climate cycles on shared demographic response and investigate the contributions of both biotic and abiotic factors which shape these demographic trajectories. Location: Madagascar. Taxon Reptiles and Amphibians. Methods: Using community‐scale population genetic data, we implement a hierarchical approximate Bayesian computation approach to evaluate the degree of synchronous population expansion during glacial cycles across herpetofaunal assemblages both within and across discrete biomes and taxonomic groups. We use Bayesian model averaging to identify intrinsic and extrinsic conditions predictive of individualistic demographic change. Results: We find that demographic responses are not uniform across groups, with more than 50% of all populations showing signal of recent expansion. Our explanatory models indicate species occupying narrow elevational ranges had a higher probability of expansion, while amphibian assemblages showed higher genetic diversity and greater departures from population neutrality. Expansion events were largely asynchronous, with coexpansion found in less than half of all populations. Exceptionally, 69%–74% of all humid‐restricted populations coexpanded during the start of the Last Glacial Period at around 100 kya, supporting the hypothesis of a more extensive humid forest cover for Madagascar during this time. Main conclusions: We show that differences in life history and regional biogeography have contributed to patterns of richness and endemism found across Madagascar, and historical connectivity across humid forests has been an important factor in shaping present‐day diversity and endemism on Madagascar, particularly for amphibians. [ABSTRACT FROM AUTHOR]

Non‐native (invasive) species offer a unique opportunity to study the geographical 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 modelling 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 (TSEalien), 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 TSEalien 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 neighbouring species despite having zones of contact. We demonstrate that these native Trachemys species have distinct climatic niche suitability, probably preventing establishment of and displacement by the TSEalien. Moreover, we found major niche overlap between TSEalien 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‐centre of the native species. [ABSTRACT FROM AUTHOR]

The statistical rigor of species delimitation has increased dramatically over the past decade. Coalescent theory provides powerful models for population genetic inference, and is now increasingly important in phylogenetics and speciation research. By applying probabilistic models, coalescent-based species delimitation provides clear and objective testing of alternative hypotheses of evolutionary independence. As acquisition of multilocus data becomes increasingly automated, coalescent-based species delimitation will improve the discovery, resolution, consistency, and stability of the taxonomy of species. Along with other tools and data types, coalescent-based species delimitation will play an important role in an integrative taxonomy that emphasizes the identification of species limits and the processes that have promoted lineage diversification. [Copyright &y& Elsevier]

Aim: To test three different scenarios to account for the geographic distribution of genetic variation in Dermatonotus muelleri, a fossorial frog endemic of the South American diagonal of open formations (DOF) formed by Chaco, Cerrado, and Caatinga biomes: (a) The pan‐DOF hypothesis, where these biomes behave as a single biogeographical unit and DOF species lack population structure; (b) the southwest‐northeast (SW‐NE) hypothesis, where structured populations along the DOF have either asymmetric or unidirectional gene flow in the northeast direction; (c) the vicariance hypothesis, where a pattern of population structure with no or reduced gene flow is observed. Location: Caatinga, Cerrado, and Chaco in South America. Methods: We sampled 179 individuals across the species distribution and sequenced one mitochondrial and two nuclear markers. We evaluated population structure, the presence, direction, and asymmetry of gene flow, and demographic changes to test three different phylogeographic hypotheses for D. muelleri using a model‐based simulation approach. Results: We found two geographically structured lineages with no gene flow, refuting pan‐DOF hypothesis. The Southwest lineage occurs in Chaco and western Cerrado and the Northeast lineage in eastern Cerrado and Caatinga. Our results support the model with demographic change for the Northeast lineage and constant size for Southwest lineage. The vicariance model with divergence time estimated around 4 Ma was preferred. Main conclusions: Vicariance was the best explanation for the geographic distribution of genetic variation in Dermatonotus. Divergence time supports the role of Neogene orogeny in central Brazil for the splitting of Dermatonotus lineages. The influence of Pleistocene climatic fluctuations seems restricted to changes in one lineage's population size. Therefore, both ancient and recent events shaped the phylogeography of Dermatonotus. Genetic divergence, lack of gene flow, and the presence of a geographic barrier between lineages indicate the existence of a cryptic, undescribed species in northeastern Brazil. [ABSTRACT FROM AUTHOR]

[Display omitted] • UCE data resolves previously problematic relationships in Thamnophiini. • Support for reversing taxonomic changes of Haldea striatula and Liodytes pygaea. • Thamnophis is the only major genus to have crossed westward multiple times. • High level of gene tree discordance with species tree despite high nodal support. 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. [ABSTRACT FROM AUTHOR]

The genus-level taxonomy of the New World racers and whipsnakes ( Coluber and Masticophis) has long been contentious regarding whether the two genera are mutually exclusive clades. This argument is based on morphological characters and largely single-locus analyses. Herein we examine the phylogenetic history of this group using multi-locus data in a coalescent framework, where paraphyly of Masticophis would result in support for the recognition of only a single genus ( Coluber) for these species. We sample all currently recognized species and incorporate broad geographic sampling for the more widespread species groups to explore biogeographic patterns across North America. Our analyses suggest that Masticophis is monophyletic with respect to Coluber constrictor, albeit with low support. These results also demonstrate that there is undescribed cryptic diversity in this group, and we underscore additional avenues of study to further delimit unrecognized species in this clade. The biogeography of the island endemic, Masticophis anthonyi, is discussed with respect to what is known about other codistributed vertebrates. Lastly we provide an overview of the history of the arguments for or against the use of the generic name Masticophis and suggest its continued use. [ABSTRACT FROM AUTHOR]

[Display omitted] • >4500 ultraconserved elements loci resolved the phylogeny of the snake superfamily Elapoidea. • Anatomical data and phylogeny supported establishment of a new Afro-Asian family of snakes, Micrelapidae. • Different phylogeny estimation methods inferred largely congruent topologies. • Time calibrated phylogeny supported the hypothesis of a rapid elapoid radiation in the Eocene. 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. [ABSTRACT FROM AUTHOR]

Both gene-tree discordance and unrecognized diversity are sources of error for accurate estimation of species trees, and can affect downstream diversification analyses by obscuring the correct number of nodes, their density, and the lengths of the branches subtending them. Although the theoretical impact of gene-tree discordance on evolutionary analyses has been examined previously, the effect of unsampled and cryptic diversity has not. Here, we examine how delimitation of previously unrecognized diversity in the milksnake (Lampropeltis triangulum) and use of a species-tree approach affects both estimation of the Lampropeltis phylogeny and comparative analyses with respect to the timing of diversification. Coalescent species delimitation indicates that L. triangulum is not monophyletic and that there are multiple species of milksnake, which increases the known species diversity in the genus Lampropeltis by 40%. Both genealogical and temporal discordance occurs between gene trees and the species tree, with evidence that mitochondrial DNA (mtDNA) introgression is a main factor. This discordance is further manifested in the preferred models of diversification, where the concatenated gene tree strongly supports an early burst of speciation during the Miocene, in contrast to species-tree estimates where diversification follows a birth–death model and speciation occurs mostly in the Pliocene and Pleistocene. This study highlights the crucial interaction among coalescent-based phylogeography and species delimitation, systematics, and species diversification analyses. [ABSTRACT FROM AUTHOR]

Spiny-rayed fishes, or acanthomorphs, comprise nearly one-third of all living vertebrates. Despite their dominant role in aquatic ecosystems, the evolutionary history and tempo of acanthomorph diversification is poorly understood. We investigate the pattern of lineage diversification in acanthomorphs by using a well-resolved time-calibrated phylogeny inferred from a nuclear gene super-matrix that includes 520 acanthomorph species and 37 fossil age constraints. This phylogeny provides resolution for what has been classically referred to as the "bush at the top" of the teleost tree, and indicates acanthomorphs originated in the Early Cretaceous. Pale-ontological evidence suggests acanthomorphs exhibit a pulse of morphological diversification following the end Cretaceous mass extinction; however, the role of this event on the accumulation of living acanthomorph diversity remains unclear. Lineage diversification rates through time exhibit no shifts associated with the end Cretaceous mass extinction, but there is a global decrease in lineage diversification rates 50 Ma that occurs during a period when morphological disparity among fossil acanthomorphs increases sharply. Analysis of clade-specific shifts in diversification rates reveal that the hyperdiversity of living acanthomorphs is highlighted by several rapidly radiating lineages including tunas, gobies, blennies, snail-fishes, and Afro-American cichlids. These lineages with high diversification rates are not associated with a single habitat type, such as coral reefs, indicating there is no single explanation for the success of acanthomorphs, as exceptional bouts of diversification have occurred across a wide array of marine and freshwater habitats. [ABSTRACT FROM AUTHOR]

[Display omitted] • Quaternary climatic fluctuations influenced in the demographic history of amphibians. • Population expansions took place during periods of climate cooling. • A microendemic lineage was recovered in a geographic isolated region. • The results found are in accordance with other studies with codistributed species. 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. [ABSTRACT FROM AUTHOR]

Convergence offers an opportunity to explore to what extent evolution can be predictable when genomic composition and environmental triggers are similar. Here, we present an emergent model system to study convergent evolution in nature in a mammalian group, the bat genus Myotis. Three foraging strategies-gleaning, trawling, and aerial hawking, each characterized by different sets of phenotypic features-have evolved independently multiple times in different biogeographic regions in isolation for millions of years. To investigate the genomic basis of convergence and explore the functional genomic changes linked to ecomorphological convergence, we sequenced and annotated 17 new genomes and screened 16,426 genes for positive selection and associations between relative evolutionary rates and foraging strategies across 30 bat species representing all Myotis ecomorphs across geographic regions as well as among sister groups. We identify genomic changes that describe both phylogenetic and ecomorphological trends. We infer that colonization of new environments may have first required changes in genes linked to hearing sensory perception, followed by changes linked to fecundity and development, metabolism of carbohydrates, and heme degradation. These changes may be linked to prey acquisition and digestion and match phylogenetic trends. Our findings also suggest that the repeated evolution of ecomorphs does not always involve changes in the same genes but rather in genes with the same molecular functions such as developmental and cellular processes., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe., Competing Interests: We declare we have no competing interests., (© 2024 The Author(s).)

Albinism is an uncommon phenomenon and inherited condition in animals characterized by a partial or complete lack of melanin. The family Xenodermidae Gray, 1849, is a group of caenophidian snakes widely distributed in South, East, and Southeast Asia, including five recognized genera and 36 species. However, there are currently no reports of albinism in any species in Xenodermidae. Achalinussheni Ma, Xu, Qi, Wang, Tang, Huang & Jiang, 2023 was first described based on five male specimens from Loudi City and Nanyue District, Hunan Province, China. At the time, there were no descriptions on female individuals. In this study, we report in detail a collected albinistic specimen of A.sheni , which is the first discovery of wild albinism in the family Xenodermidae. We also provide photographs and descriptions of the first three female specimens of A.sheni and extend the diagnosis of this species., Competing Interests: The authors have declared that no competing interests exist., (Yu-Hao Xu, Shuai Wang, Shun Ma, Frank T. Burbrink, Li-Fang Peng, Song Huang.)

The South American Dry Diagonal, also called the Diagonal of Open Formations, is a large region of seasonally dry vegetation extending from northeastern Brazil to northern Argentina, comprising the Caatinga, Cerrado, and Chaco subregions. A growing body of phylogeography literature has determined that a complex history of climatic changes coupled with more ancient geological events has produced a diverse and endemic-rich Dry Diagonal biota. However, the exact drivers are still under investigation, and their relative strengths and effects are controversial. Pleistocene climatic fluctuations structured lineages via vegetation shifts, refugium formation, and corridors between the Amazon and Atlantic forests. In some taxa, older geological events, such as the reconfiguration of the São Francisco River, uplift of the Central Brazilian Plateau, or the Miocene inundation of the Chaco by marine incursions, were more important. Here, we review the Dry Diagonal phylogeography literature, discussing each hypothesized driver of diversification and assessing degree of support. Few studies statistically test these hypotheses, with most support drawn from associating encountered phylogeographic patterns such as population structure with the timing of ancient geoclimatic events. Across statistical studies, most hypotheses are well supported, with the exception of the Pleistocene Arc Hypothesis. However, taxonomic and regional biases persist, such as a proportional overabundance of herpetofauna studies, and the under-representation of Chaco studies. Overall, both Pleistocene climate change and Neogene geological events shaped the evolution of the Dry Diagonal biota, though the precise effects are regionally and taxonomically varied. We encourage further use of model-based analyses to test evolutionary scenarios, as well as interdisciplinary collaborations to progress the field beyond its current focus on the traditional set of geoclimatic hypotheses., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

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 ., Competing Interests: The authors declare no competing interests., (© 2024 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.)

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., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

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., Competing Interests: We have no competing interests., (© 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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 12 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 are found at nodes with high gene tree discordance and occurred among parental lineages distributed along a north-south axis that share similar ecologies. Younger hybrids occupy intermediate contact zones between the parent lineages 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 rates may decline as the result of species accumulation and environmental change, speciation by hybridization may be more robust to these processes., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

It was not known that phylogenetic analysis of population genetic data would show geographic structure when Avise et al. (1979) introduced the unnamed field of phylogeography. Since then, the field has grown from using single gene fragments to whole genomes resulting in more than 22,000 publications. This seemingly simple relationship between geography and genetic variation has provided the foundation for studying speciation, species delimitation, hybrid zone dynamics, adaptation, conservation genetics, community assembly, and historical demography to name a few (see Burbrink & Ruane, 2021; Edwards et al., 2016). As with many burgeoning fields, there is often little consideration of how to make data sets accessible for future researchers addressing more comprehensive questions under a common framework. For example, it is important in biogeography to examine how shared environments or barriers affect population structure across communities of species, or test if range size or latitude are correlated with genetic diversity across taxa (Myers et al., 2019; Smith et al., 2017). However, addressing these types of questions using existing data requires researchers to assemble large databases manually. Genetic and geographic databases used to store this information like NCBI Genbank (National Center for Biotechnology Information) and GBIF (Global Biodiversity Information Facility) are disconnected, and often of limited general use for conducting multitaxon studies. In this issue of Molecular Ecology Resources, Pelletier et al. (2022) have automated the process of connecting geography to DNA sequences via the phylogatR (phylogeographic data aggregation and repurposing) database., (© 2022 John Wiley & Sons Ltd.)

A new species of odd-scaled snake in the genus Achalinus is described from Ningshan County, Shaanxi Province, China, based on seven female specimens. Morphologically, the new species can be distinguished from its congeners by the following combination of characters: (1) dorsum iridescent and uniformly dark brown in preservative, longitudinal vertebral line absent; (2) light brown beneath; (3) dotted black streak in the middle of the subcaudals absent; (4) tail length relatively short, TaL/ToL 12-16%, in females; (5) fewer subcaudals, 41-46 in females; (6) dorsal scales 23 rows throughout, strongly keeled; the outer-most rows on both sides are also keeled and slightly enlarged; (7) one loreal; (8) internasal not fused to prefrontal; (9) suture between internasals is similar in size when compared to the suture between prefrontals; (10) preocular and postocular absent; (11) 6 supralabials; (12) 5 infralabials, the first 3 (rarely 2) touching the first pair of chin shields; (13) 3 pairs of chin shields. A phylogenetic analysis using mitochondrial gene, cytochrome c oxidase subunit 1 (CO1), showed that this new species forms a monophyletic group with strong support. In addition, the uncorrected p-distances between the new species and other known congeners ranged from 3.4% to 13.0%. The recognition of the new species increases the number of described Achalinus species to 21.

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., Competing Interests: None declared., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., (© 2020 The Authors. Evolution © 2020 The Society for the Study of Evolution.)

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)., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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.

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., (© 2015 The Author(s).)

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., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Adaptive radiation (AR) theory predicts that groups sharing the same source of ecological opportunity (EO) will experience deterministic species diversification and morphological evolution. Thus, deterministic ecological and morphological evolution should be correlated with deterministic patterns in the tempo and mode of speciation for groups in similar habitats and time periods. We test this hypothesis using well-sampled phylogenies of four squamate groups that colonized the New World (NW) in the Late Oligocene. We use both standard and coalescent models to assess species diversification, as well as likelihood models to examine morphological evolution. All squamate groups show similar early pulses of speciation, as well as diversity-dependent ecological limits on clade size at a continental scale. In contrast, processes of morphological evolution are not easily predictable and do not show similar pulses of early and rapid change. Patterns of morphological and species diversification thus appear uncoupled across these groups. This indicates that the processes that drive diversification and disparification are not mechanistically linked, even among similar groups of taxa experiencing the same sources of EO. It also suggests that processes of phenotypic diversification cannot be predicted solely from the existence of an AR or knowledge of the process of diversification.

Snake diversity varies by at least two orders of magnitude among extant lineages, with numerous groups containing only one or two species, and several young clades exhibiting exceptional richness (>700 taxa). With a phylogeny containing all known families and subfamilies, we find that these patterns cannot be explained by background rates of speciation and extinction. The majority of diversity appears to derive from a radiation within the superfamily Colubroidea, potentially stemming from the colonization of new areas and the evolution of advanced venom-delivery systems. In contrast, negative relationships between clade age, clade size, and diversification rate suggest the potential for possible bias in estimated diversification rates, interpreted by some recent authors as support for ecologically mediated limits on diversity. However, evidence from the fossil record indicates that numerous lineages were far more diverse in the past, and that extinction has had an important impact on extant diversity patterns. Thus, failure to adequately account for extinction appears to prevent both rate- and diversity-limited models from fully characterizing richness dynamics in snakes. We suggest that clade-level extinction may provide a key mechanism for explaining negative or hump-shaped relationships between clade age and diversity, and the prevalence of ancient, species-poor lineages in numerous groups., (© 2011 The Author(s). Evolution © 2011 The Society for the Study of Evolution.)

Molecular phylogenies are often used to test hypotheses about the tempo and mode of speciation and extinction. One commonly used statistic is Pybus and Harvey's γ, which measures the density of ordered internode distances on an ultrametric tree to infer earlier (negative γ) or later (positive γ) bursts of diversification. However, coalescent theory predicts that γ might be biased toward negative values (inferring early bursts of diversification) when using gene trees rather than species trees. Gene divergences predate species divergences, increasingly so at higher effective population sizes (N(e)), and proportionally more so toward the tips of the tree. Thus, gene trees will have a higher density of older nodes in many cases (particularly at higher N(e)), due to the disproportionate lengthening of terminal branches. This will yield an artifactual signature of early bursts of diversification when estimating γ from gene trees. We simulate gene trees within species trees under both Yule (pure-birth) and birth-death processes, and demonstrate support for these predictions. However, for most realistic estimates of θ in natural populations, gene trees provide relatively good estimates of γ, despite the disproportionate overestimation of younger node ages. This is corroborated with an empirical dataset of North American fence lizards (Sceloporus)., (© 2011 The Author(s). Evolution© 2011 The Society for the Study of Evolution.)

Ecological adaptive radiation theory predicts an increase in both morphological and specific diversification when organisms colonize new environments. Accordingly, bursts of morphological diversification, characterized by low within-subclade morphological disparity, may be associated with these increases in speciation rates. Conversely, increasing species density, reduction in available habitat, or increasing extinction rates are expected to cause rates of diversification to decline. We test these hypotheses by examining the tempo and mode of speciation in the lampropeltinine snakes, a morphologically variable group that colonized the New World approximately 24 million years ago and radiated throughout the Miocene. We show that specific diversification increased early in the history of the group, and that most morphological variation is partitioned among, rather than within subclades. These patterns provide further evidence for the hypothesis that morphological variation tends to be strongly partitioned among lineages when clades undergo early bursts of species diversification. A reduction in speciation rates may be indicative of density dependent effects due to a saturation of available ecological opportunity, rather than increases in extinction rates at the onset of the Pleistocene/Pliocene glacial cycles. This evidence runs counter to the general Pleistocene species pump model.

Niche conservatism and niche divergence are both important ecological mechanisms associated with promoting allopatric speciation across geographical barriers. However, the potential for variable responses in widely distributed organisms has not been fully investigated. For allopatric sister lineages, three patterns for the interaction of ecological niche preference and geographical barriers are possible: (i) niche conservatism at a physical barrier; (ii) niche divergence at a physical barrier; and (iii) niche divergence in the absence of a physical barrier. We test for the presence of these patterns in a transcontinentally distributed snake species, the common kingsnake (Lampropeltis getula), to determine the relative frequency of niche conservatism or divergence in a single species complex inhabiting multiple distinct ecoregions. We infer the phylogeographic structure of the kingsnake using a range-wide data set sampled for the mitochondrial gene cytochrome b. We use coalescent simulation methods to test for the presence of structured lineage formation vs. fragmentation of a widespread ancestor. Finally, we use statistical techniques for creating and evaluating ecological niche models to test for conservatism of ecological niche preferences. Significant geographical structure is present in the kingsnake, for which coalescent tests indicate structured population division. Surprisingly, we find evidence for all three patterns of conservatism and divergence. This suggests that ecological niche preferences may be labile on recent phylogenetic timescales, and that lineage formation in widespread species can result from an interaction between inertial tendencies of niche conservatism and natural selection on populations in ecologically divergent habitats.