Ecological divergence in the yellow-bellied kingsnake (Lampropeltis calligaster) at two North American biodiversity hotspots

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.

Vicariance, dispersal, extinction and hybridization underlie the evolutionary history of Atlantic forest fire-eye antbirds (Aves: Thamnophilidae)

In order to gain insights into the biogeographic processes underlying biotic diversification in the Atlantic Forest (AF), we used a multi-locus approach to examine the evolutionary history of the White-shouldered Fire-eye (Pyriglena leucoptera) and the Fringe-backed Fire-eye (Pyriglena atra), two parapatric sister species endemic to the AF. We sequenced one mitochondrial, three Z chromosome-linked and three anonymous markers of 556 individuals from 66 localities. We recovered four lineages throughout the AF: P. atra and three populations within P. leucoptera. All populations diverged during the late Pleistocene and presented varying levels of admixture. One Z-linked locus showed the highest level of differentiation between the two species. On the other hand, a mitochondrial haplotype was shared extensively between them. Our data supported vicariance driving speciation along with extinction and dispersal as processes underlying intraspecific diversification. Furthermore, signatures of demographic expansion in most populations and areas of genetic admixture were recovered throughout the AF, suggesting that forest fragmentation was also important in differentiation. Genetic admixture areas are located between large rivers suggesting that AF rivers may diminish gene flow. Our results indicated a complex and dynamic biogeographic history of Pyriglena in the AF, with vicariance, extinction, dispersal and secondary contact followed by introgression likely influencing the current patterns of genetic distribution.

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.

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.

Taxonomic and conservation implications of population genetic admixture, mito-nuclear discordance, and male-biased dispersal of a large endangered snake, Drymarchon couperi

Accurate species delimitation and description are necessary to guide effective conservation of imperiled species, and this synergy is maximized when multiple data sources are used to delimit species. We illustrate this point by examining Drymarchon couperi (Eastern Indigo Snake), a large, federally-protected species in North America that was recently divided into two species based on gene sequence data from three loci and heuristic morphological assessment. Here, we re-evaluate the two-species hypothesis for D. couperi by evaluating both population genetic and gene sequence data. Our analyses of 14 microsatellite markers revealed 6–8 genetic population clusters with significant admixture, particularly across the contact zone between the two hypothesized species. Phylogenetic analyses of gene sequence data with maximum-likelihood methods suggested discordance between mitochondrial and nuclear markers and provided phylogenetic support for one species rather than two. For these reasons, we place Drymarchon kolpobasileus into synonymy with D. couperi. We suggest inconsistent patterns between mitochondrial and nuclear DNA are driven by high dispersal of males relative to females. We advocate for species delimitation exercises that evaluate admixture and gene flow in addition to phylogenetic analyses, particularly when the latter reveal monophyletic lineages. This is particularly important for taxa, such as squamates, that exhibit strong sex-biased dispersal. Problems associated with over-delimitation of species richness can become particularly acute for threatened and endangered species, because of high costs to conservation when taxonomy demands protection of more individual species than are supported by accumulating data.

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.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Recent advances in conservation and population genomics data analysis

New computational methods and next‐generation sequencing (NGS) approaches have enabled the use of thousands or hundreds of thousands of genetic markers to address previously intractable questions. The methods and massive marker sets present both new data analysis challenges and opportunities to visualize, understand, and apply population and conservation genomic data in novel ways. The large scale and complexity of NGS data also increases the expertise and effort required to thoroughly and thoughtfully analyze and interpret data. To aid in this endeavor, a recent workshop entitled “Population Genomic Data Analysis,” also known as “ConGen 2017,” was held at the University of Montana. The ConGen workshop brought 15 instructors together with knowledge in a wide range of topics including NGS data filtering, genome assembly, genomic monitoring of effective population size, migration modeling, detecting adaptive genomic variation, genomewide association analysis, inbreeding depression, and landscape genomics. Here, we summarize the major themes of the workshop and the important take‐home points that were offered to students throughout. We emphasize increasing participation by women in population and conservation genomics as a vital step for the advancement of science. Some important themes that emerged during the workshop included the need for data visualization and its importance in finding problematic data, the effects of data filtering choices on downstream population genomic analyses, the increasing availability of whole‐genome sequencing, and the new challenges it presents. Our goal here is to help motivate and educate a worldwide audience to improve population genomic data analysis and interpretation, and thereby advance the contribution of genomics to molecular ecology, evolutionary biology, and especially to the conservation of biodiversity.