Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus).

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium.

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

RESUMEN: El veneno de las serpientes puede variar entre y dentro de las especies. Mientras algunos grupos de viperidos del Nuevo Mundo­como las cascabeles­han sido bien estudiadas, muy poco se sabe acerca del veneno de las nauyacas de frío (Cerrophidion) que se encuentran en las zonas altas de Mesoamérica. Comparadas con las extensamente estudiadas cascabeles, que estan ampliamente distribuidas, las poblaciones de Cerrophidion, aisladas en montañas, pueden poseer trayectorias evolutivas y diferenciación en su veneno unicos. En el presente trabajo, describimos el transcriptoma de las glándulas de veneno de poblaciones de C. petlalcalensis, C. tzotzilorum, y C. godmani de México, y un individuo de C. sasai de Costa Rica. Exploramos la variación en la expresión de toxinas en Cerrophidion y la evolución en las secuencias geneticas en C. godmani específicamente. El transcriptoma de la glándula de veneno de Cerrophidion esta compuesto principalmente de Metaloproteinasas de Veneno de Serpiente, Fosfolipasas A[Formula: see text] (PLA[Formula: see text]s), y Serin Proteasas de Veneno de Serpiente. Cerrophidion petlalcalensis presenta poca variación intraespecífica; sin embargo, los transcriptomas de la glandula de veneno de C. godmani y C. tzotzilorum difieren significativamente entre poblaciones geográficamente aisladas. Curiosamente, la variación intraespecífica estuvo atribuida principalmente a la expresión de las toxinas ya que no encontramos señales de selección en las toxinas de C. godmani. Adicionalmente, encontramos miotoxinas similares a PLA[Formula: see text] en todas las especies excepto C. petlalcalensis, y PLA[Formula: see text]s similares a crotoxina en la población sureña de C. godmani. Nuestros resultados demuestran la presencia de variacion intraespecífica presente en el veneno de C. godmani y C. tzotzilorum. Las toxinas de Cerrophidion godmani muestran poca evidencia de selección direccional, y la variación en la secuencias de las toxinas es consistente con evolucion bajo un modelo de equilibrio de mutación-deriva. Algunos individuos de C. godmani de la población del sur potencialmente tienen un veneno neurotóxico dada la presencia de PLA[Formula: see text]s similares a la crotoxina, sin embargo, se necesita más evidencia para corroborar esta hipótesis.

Interspecific Gene Flow and Mitochondrial Genome Capture during the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae).

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

Predictors of genomic diversity within North American squamates.

Comparisons of intraspecific genetic diversity across species can reveal the roles of geography, ecology, and life history in shaping biodiversity. The wide availability of mitochondrial DNA (mtDNA) sequences in open-access databases makes this marker practical for conducting analyses across several species in a common framework, but patterns may not be representative of overall species diversity. Here, we gather new and existing mtDNA sequences and genome-wide nuclear data (genotyping-by-sequencing; GBS) for 30 North American squamate species sampled in the Southeastern and Southwestern United States. We estimated mtDNA nucleotide diversity for 2 mtDNA genes, COI (22 species alignments; average 16 sequences) and cytb (22 species; average 58 sequences), as well as nuclear heterozygosity and nucleotide diversity from GBS data for 118 individuals (30 species; 4 individuals and 6,820 to 44,309 loci per species). We showed that nuclear genomic diversity estimates were highly consistent across individuals for some species, while other species showed large differences depending on the locality sampled. Range size was positively correlated with both cytb diversity (phylogenetically independent contrasts: R2 = 0.31, P = 0.007) and GBS diversity (R2 = 0.21; P = 0.006), while other predictors differed across the top models for each dataset. Mitochondrial and nuclear diversity estimates were not correlated within species, although sampling differences in the data available made these datasets difficult to compare. Further study of mtDNA and nuclear diversity sampled across species' ranges is needed to evaluate the roles of geography and life history in structuring diversity across a variety of taxonomic groups.

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

Giant Tree Frog diversification in West and Central Africa: Isolation by physical barriers, climate, and reproductive traits.

Secondary sympatry amongst sister lineages is strongly associated with genetic and ecological divergence. This pattern suggests that for closely related species to coexist in secondary sympatry, they must accumulate differences in traits that mediate ecological and/or reproductive isolation. Here, we characterized inter- and intraspecific divergence in three giant tree frog species whose distributions stretch across West and Central Africa. Using genome-wide single-nucleotide polymorphism data, we demonstrated that species-level divergence coincides temporally and geographically with a period of large-scale forest fragmentation during the late Pliocene. Our environmental niche models further supported a dynamic history of climatic suitability and stability, and indicated that all three species occupy distinct environmental niches. We found modest morphological differentiation amongst the species with significant divergence in tympanum diameter and male advertisement call. In addition, we confirmed that two species occur in secondary sympatry in Central Africa but found no evidence of hybridization. These patterns support the hypothesis that cycles of genetic exchange and isolation across West and Central Africa have contributed to globally significant biodiversity. Furthermore, divergence in both ecology and reproductive traits appear to have played important roles in maintaining distinct lineages. At the intraspecific level, we found that climatic refugia, precipitation gradients, marine incursions, and potentially riverine barriers generated phylogeographic structure throughout the Pleistocene and into the Holocene. Further studies examining phenotypic divergence and secondary contact amongst these geographically structured populations may demonstrate how smaller scale and more recent biogeographic barriers contribute to regional diversification.

La sympatrie secondaire parmi les espèces sœurs est fortement associée à la divergence génétique et écologique. Ce modèle suggère que pour que des espèces étroitement liées coexistent en sympatrie secondaire, elles doivent accumuler des différences dans les traits qui contribuent à l'isolement écologique ou reproductif. Ici, nous avons caractérisé la divergence inter- et intra-spécifique chez trois espèces de grenouilles arboricoles géantes dont les distributions s'étendent à travers l'Afrique de l'Ouest et Centrale. Avec des données génétiques, nous avons démontré que la divergence au niveau des espèces coïncide temporellement et géographiquement avec une période de fragmentation forestière à la fin du Pliocène. Nos modèles de niches environnementales ont soutenu une histoire dynamique de stabilité climatique, et ont indiqué que les trois espèces occupent des niches environnementales distinctes. Nous avons trouvé une différenciation morphologique modeste parmi les trois espèces mais une divergence significative dans le diamètre du tympan et les cris des mâles. De plus, nous avons confirmé que deux espèces sont présentes en sympatrie secondaire en Afrique Centrale mais n'avons trouvé aucune preuve d'hybridation. Ces résultats soutiennent l'hypothèse que les cycles d'échange génétique et d'isolement à travers l'Afrique de l'Ouest et Centrale ont contribué à une profonde concentration de biodiversité dans la région. De plus, la divergence des traits écologiques et reproducteurs semble avoir joué un rôle important dans le maintien de lignées distinctes. Au niveau intra-spécifique, nous avons constaté que les refuges climatiques, les gradients de précipitation, les incursions marines et potentiellement les barrières fluviales ont généré une structure phylogéographique pendant le Pléistocène et jusqu'à l'Holocène. Des études examinant la divergence phénotypique et le contact secondaire entre ces populations géographiquement structurées pourraient démontrer comment des barrières biogéographiques à échelle plus petite et plus récentes contribuent à la diversification régionale.

Genome-wide data reveal extensive gene flow during the diversification of the western rattlesnakes (Viperidae: Crotalinae: Crotalus).

Hybridization and introgression are important, but often overlooked processes when inferring phylogenies. When these processes are not accounted for and a strictly diverging phylogenetic model is applied to groups with a history of hybridization, phylogenetic inference and parameter estimation can be inaccurate. Recent developments in phylogenetic network approaches coupled with the increasing availability of genomic data allow inferences of reticulate evolutionary histories across the tree of life. The western rattlesnake species group (C. viridis species complex, C. mitchellii species complex, C. scutulutas, and C. tigris) is an iconic snake lineage that is widespread across western North America. This group is composed of several species complexes with unclear species limits, likely the result of ongoing gene flow among nascent lineages. Here I generate reduced representation genomic data and test for a history of reticulation within this group. I demonstrate that all species have undergone hybridization with at least one other lineage, suggesting introgression is widespread in this group. Topologies differ between phylogenies estimated under the multispecies coalescent and multispecies network coalescent methods, indicating that gene flow has obscured phylogenetic relationships within this group. These past introgression events are predominantly restricted to species that co-occur geographically. However, within species that have a history of introgression, this signature is detected regardless of specimen sampling across geography. Overall, my results suggest the accumulation of reproductive isolating barriers occurs slowly in rattlesnakes which likely leads to the difficulty in delimiting species, furthermore, the results of this study have implications for trait evolution in this group.

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.

Environmental heterogeneity and not vicariant biogeographic barriers generate community-wide population structure in desert-adapted snakes.

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community-wide divergence should be explored before assuming the role of biogeographic barriers.

Exploring Chihuahuan Desert diversification in the gray-banded kingsnake, Lampropeltis alterna (Serpentes: Colubridae).

Within many biomes, the cause of phylogeographic structure remains unknown even across regions throughout North America, including within the biodiverse Chihuahuan Desert. For example, little is known about population structure or the timing of diversification of Chihuahuan endemics. This is due largely to the lack of population genomic studies within this region. We generated ultra-conserved element data for the gray-banded kingsnake (Lampropeltis alterna) to investigate lineage divergence and historical demography across the Chihuahuan Desert. We found three unique lineages corresponding to the Trans-Pecos and Mapimian biogeographic regions of the Chihuahuan Desert, and a distinct population in the Sierra Madre Occidental. Using several mutation rates to calibrate the timing of divergence among these lineages, we show that lineage divergence likely occurred during the Pleistocene, which indicates that careful consideration needs to be used when applying mutation rates to ultra-conserved elements. We suggest that biogeographic provinces within the Chihuahuan Desert may have served as allopatric refugia during climatic fluctuations of the Quaternary. This work serves as an important template for further testing biogeographic hypotheses within the region.

Asynchronous demographic responses to Pleistocene climate change in Eastern Nearctic vertebrates.

Pleistocene climatic cycles altered species distributions in the Eastern Nearctic of North America, yet the degree of congruent demographic response to the Pleistocene among codistributed taxa remains unknown. We use a hierarchical approximate Bayesian computational approach to test if population sizes across lineages of snakes, lizards, turtles, mammals, birds, salamanders and frogs in this region expanded synchronously to Late Pleistocene climate changes. Expansion occurred in 75% of 74 lineages, and of these, population size trajectories across the community were partially synchronous, with coexpansion found in at least 50% of lineages in each taxonomic group. For those taxa expanding outside of these synchronous pulses, factors related to when they entered the community, ecological thresholds or biotic interactions likely condition their timing of response to Pleistocene climate change. Unified timing of population size change across communities in response to Pleistocene climate cycles is likely rare in North America.

Complex longitudinal diversification across South China and Vietnam in Stejneger's pit viper, Viridovipera stejnegeri (Schmidt, 1925) (Reptilia: Serpentes: Viperidae).

Viridovipera stejnegeri is one of the most common pit vipers in Asia, with a wide distribution in southern China and Vietnam. We investigated historical demography and explored how the environment and climatic factors have shaped genetic diversity and the evolutionary history of this venomous snake. A total of 171 samples from 47 localities were sequenced and analysed for two mitochondrial gene fragments and three nuclear genes. Gene trees reveal the existence of two well-supported clades (Southwest China and Southeast China) with seven distinct and strongly supported, geographically structured subclades within V. stejnegeri. Estimation of divergence time and ancestral area suggests that V. stejnegeri originated at ~6.0 Ma in the late Miocene on the Yunnan-Guizhou Plateau. The estimated date of origin and divergence of the island populations of Taiwan and Hainan closely matches the geological origin of the both islands. The mtDNA gene tree reveals the presence of west-east diversification in V. stejnegeri populations. Complex orogenesis and heterogeneous habitats, as well as climate-mediated habitat differentiation including glacial cycles, all have influenced population structure and the distribution of this taxon. The validity of V. stejnegeri chenbihuii is questionable, and this subspecies most probably represents an invalid taxon.

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.

Assessing species boundaries and the phylogenetic position of the rare Szechwan ratsnake, Euprepiophis perlaceus (Serpentes: Colubridae), using coalescent-based methods.

Delimiting species and clarifying phylogenetic relationships are the main goals of systematics. For species with questionable taxonomic status, species delimitation approaches using multi-species coalescent models with multiple loci are recommended if morphological data are unavailable or unhelpful. Moreover, these methods will also reduce subjectivity based on genetic distance or requirement of monophyletic genetic lineages. We determine the validity and phylogenetic position of a rare and long controversial species of Chinese reptile, the Szechwan ratsnake (Euprepiophis perlaceus), using multi-locus data from multiple individuals and coalescent-based approaches. Species were first delimited using Bayesian Phylogenetics & Phylogeography (BP&P), Brownie and Bayes Factor model comparison approaches, while relationships among species were estimated using species tree inference in (*)BEAST. Results indicate that Euprepiophis perlaceus is a distinct species sister to Euprepiophis mandarinus. Despite gene tree discrepancy, the coalescent model-based approaches used here demonstrate the taxonomic validity and the phylogenetic position of Euprepiophis perlaceus. These approaches objectively test the validity of questionable species diagnoses based on morphological characters and determine their phylogenetic position.

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.

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

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

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.

Evidence for determinism in species diversification and contingency in phenotypic evolution during adaptive radiation.

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.

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.