Does size matter? An analysis of the niche width and vulnerability to climate change of fourteen species of the genus Crotalus from North America.

The niche comprises the set of abiotic and biotic environmental conditions in which a species can live. Consequently, those species that present broader niches are expected to be more tolerant to changes in climatic variations than those species that present reduced niches. In this study, we estimate the amplitude of the climatic niche of fourteen species of rattlesnakes of the genus Crotalus to evaluate whether those species that present broader niches are less susceptible to the loss of climatically suitable zones due to the projected climate change for the time period 2021-2040. Our results suggest that for the species under study, the breadth of the niche is not a factor that determines their vulnerability to climatic variations. However, 71.4% of the species will experience increasingly inadequate habitat conditions, mainly due to the increase in temperature and the contribution that this variable has in the creation of climatically suitable zones for most of these species.

Cryptic diversity across the Trans-Mexican Volcanic Belt of Mexico in the montane bunchgrass lizard Sceloporus subniger/ (Squamata: Phrynosomatidae).

Sceloporus subniger Poglaygen Smith is a montane bunchgrass lizard distributed across pine-oak forests of central Mexico. Prompted by the discovery of a new population of this lizard in far western Mexico, and by recent studies suggesting S. subniger may be a composite of several distinct species, we examined in more detail the genetic structure of S. subniger. We generated a mitochondrial DNA (mtDNA) dataset from 81 specimens and an ultraconserved elements (UCE) dataset representing thousands of genomic regions from 12 specimens to specifically evaluate the genetic distinctiveness of populations from western Michoacán and adjacent Jalisco along with the newly discovered population in the Sierra de Mascota in western Jalisco. We also recorded morphological data from 47 museum specimens to compare to our genetic data. Results from our analyses of the genetic data, augmented by specimen measurements and scale counts, support the notion that S. subniger is indeed a composite of distinct species. Montane bunchgrass lizards from western Michoacán and adjacent Jalisco, and from the Sierra de Mascota in western Jalisco, each represent distinct new species, which we describe and name here.

Dense Geographic and Genomic Sampling Reveals Paraphyly and a Cryptic Lineage in a Classic Sibling Species Complex.

Incomplete or geographically biased sampling poses significant problems for research in phylogeography, population genetics, phylogenetics, and species delimitation. Despite the power of using genome-wide genetic markers in systematics and related fields, approaches such as the multispecies coalescent remain unable to easily account for unsampled lineages. The Empidonax difficilis/Empidonax occidentalis complex of small tyrannid flycatchers (Aves: Tyrannidae) is a classic example of widely distributed species with limited phenotypic geographic variation that was broken into two largely cryptic (or "sibling") lineages following extensive study. Though the group is well-characterized north of the US Mexico border, the evolutionary distinctiveness and phylogenetic relationships of southern populations remain obscure. In this article, we use dense genomic and geographic sampling across the majority of the range of the E. difficilis/E. occidentalis complex to assess whether current taxonomy and species limits reflect underlying evolutionary patterns, or whether they are an artifact of historically biased or incomplete sampling. We find that additional samples from Mexico render the widely recognized species-level lineage E. occidentalis paraphyletic, though it retains support in the best-fit species delimitation model from clustering analyses. We further identify a highly divergent unrecognized lineage in a previously unsampled portion of the group's range, which a cline analysis suggests is more reproductively isolated than the currently recognized species E. difficilis and E. occidentalis. Our phylogeny supports a southern origin of these taxa. Our results highlight the pervasive impacts of biased geographic sampling, even in well-studied vertebrate groups like birds, and illustrate what is a common problem when attempting to define species in the face of recent divergence and reticulate evolution.

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.

Cryptic diversity in the Mexican highlands: Thousands of UCE loci help illuminate phylogenetic relationships, species limits and divergence times of montane rattlesnakes (Viperidae: Crotalus).

With the continued adoption of genome-scale data in evolutionary biology comes the challenge of adequately harnessing the information to make accurate phylogenetic inferences. Coalescent-based methods of species tree inference have become common, and concatenation has been shown in simulation to perform well, particularly when levels of incomplete lineage sorting are low. However, simulation conditions are often overly simplistic, leaving empiricists with uncertainty regarding analytical tools. We use a large ultraconserved element data set (>3,000 loci) from rattlesnakes of the Crotalus triseriatus group to delimit lineages and estimate species trees using concatenation and several coalescent-based methods. Unpartitioned and partitioned maximum likelihood and Bayesian analysis of the concatenated matrix yield a topology identical to coalescent analysis of a subset of the data in bpp. ASTRAL analysis on a subset of the more variable loci also results in a tree consistent with concatenation and bpp, whereas the SVDquartets phylogeny differs at additional nodes. The size of the concatenated matrix has a strong effect on species tree inference using SVDquartets, warranting additional investigation on optimal data characteristics for this method. Species delimitation analyses suggest up to 16 unique lineages may be present within the C. triseriatus group, with divergences occurring during the Neogene and Quaternary. Network analyses suggest hybridization within the group is relatively rare. Altogether, our results reaffirm the Mexican highlands as a biodiversity hotspot and suggest that coalescent-based species tree inference on data subsets can provide a strongly supported species tree consistent with concatenation of all loci with a large amount of missing data.

Integrative species delimitation in practice: Revealing cryptic lineages within the short-nosed skink Plestiodon brevirostris (Squamata: Scincidae).

Integrative taxonomy has been generally considered as a goal in systematics for more than a decade. Here, we employed environmental, molecular, and morphological data to evaluate the species boundaries within the short-nosed skink Plestiodon brevirostris from south-central Mexico, one member of the morphologically conservative P. brevirostris group. Our molecular dataset includes one mitochondrial and two nuclear loci. The mitochondrial fragment includes the full length of the gene coding for the NADH dehydrogenase subunit 1 protein, a segment of the gene coding for 16S ribosomal RNA, and flanking tRNAs. The nuclear dataset includes fragments of the genes coding for the megakaryoblastic leukemia 1 and RNA fingerprint 35 proteins. We employed phylogenetic reconstruction, analyses of population structure and morphological variation, and species delimitation methods (including the integration of the three kinds of data in a unified probabilistic framework) to evaluate species limits. Our results suggest that P. brevirostris represents four distinct species. The information provided by each kind of data allowed us to discern between alternative explanations for the observed patterns of geographic structure. Two of the newly recognized lineages are poorly differentiated morphologically but apparently differ in environmental preferences and are allopatric. Additionally, one lineage is microendemic and parapatric with respect to another one. Moreover, our phylogenetic analyses suggest that other taxa within the P. brevirostris group may represent species complexes. We discuss our results in the context of integrative species delimitation.

Venom Ontogeny in the Mexican Lance-Headed Rattlesnake (Crotalus polystictus).

As trophic adaptations, rattlesnake venoms can vary in composition depending on several intrinsic and extrinsic factors. Ontogenetic changes in venom composition have been documented for numerous species, but little is known of the potential age-related changes in many rattlesnake species found in México. In the current study, venom samples collected from adult and neonate Crotalus polystictus from Estado de México were subjected to enzymatic and electrophoretic analyses, toxicity assays (LD50), and MALDI-TOF mass spectrometry, and a pooled sample of adult venom was analyzed by shotgun proteomics. Electrophoretic profiles of adult males and females were quite similar, and only minor sex-based variation was noted. However, distinct differences were observed between venoms from adult females and their neonate offspring. Several prominent bands, including P-I and P-III snake venom metalloproteinases (SVMPs) and disintegrins (confirmed by MS/MS) were present in adult venoms and absent/greatly reduced in neonate venoms. Age-dependent differences in SVMP, kallikrein-like, phospholipase A2 (PLA2), and L-amino acid oxidase (LAAO) activity levels were confirmed by enzymatic activity assays, and like many other rattlesnake species, venoms from adult snakes have higher SVMP activity than neonate venoms. Conversely, PLA2 activity was approximately 2.5 × greater in venoms from neonates, likely contributing to the increased toxicity (neonate venom LD50 = 4.5 μg/g) towards non-Swiss albino mice when compared to adult venoms (LD50 = 5.5 μg/g). Thrombin-like (TLE) and phosphodiesterase activities did not vary significantly with age. A significant effect of sex (between adult male and adult female venoms) was also observed for SVMP, TLE, and LAAO activities. Analysis of pooled adult venom by LC-MS/MS identified 14 toxin protein families, dominated by bradykinin-inhibitory peptides, SVMPs (P-I, P-II and P-III), disintegrins, PLA2s, C-type-lectins, CRiSPs, serine proteinases, and LAAOs (96% of total venom proteins). Neonate and adult C. polystictus in this population consume almost exclusively mammals, suggesting that age-based differences in composition are related to physical differences in prey (e.g., surface-to-volume ratio differences) rather than taxonomic differences between prey. Venoms from adult C. polystictus fit a Type I pattern (high SVMP activity, lower toxicity), which is characteristic of many larger-bodied rattlesnakes of North America.

Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus).

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene.

Genome-wide SNP data and morphology support the distinction of two new species of Kovarikia Soleglad, Fet & Graham, 2014 endemic to California (Scorpiones, Vaejovidae).

Morphologically conserved taxa such as scorpions represent a challenge to delimit. We recently discovered populations of scorpions in the genus Kovarikia Soleglad, Fet & Graham, 2014 on two isolated mountain ranges in southern California. We generated genome-wide single nucleotide polymorphism data and used Bayes factors species delimitation to compare alternative species delimitation scenarios which variously placed scorpions from the two localities with geographically adjacent species or into separate lineages. We also estimated a time-calibrated phylogeny of Kovarikia and examined and compared the morphology of preserved specimens from across its distribution. Genetic results strongly support the distinction of two new lineages, which we describe and name here. Morphology among the species of Kovarikia was relatively conserved, despite deep genetic divergences, consistent with recent studies of stenotopic scorpions with limited vagility. Phylogeographic structure discovered in several previously described species also suggests additional cryptic species are probably present in the genus.

Species delimitation and biogeography of the gnatcatchers and gnatwrens (Aves: Polioptilidae).

The New World avian family Polioptilidae (gnatcatchers and gnatwrens) is distributed from Argentina to Canada and includes 15 species and more than 60 subspecies. No study to date has evaluated phylogenetic relationships within this family and the historical pattern of diversification within the group remains unknown. Moreover, species limits, particularly in widespread taxa that show geographic variation, remain unclear. In this study, we delimited species and estimated phylogenetic relationships using multilocus data for the entire family. We then used the inferred diversity along with alternative taxonomic classification schemes to evaluate how lumping and splitting of both taxa and geographical areas influenced biogeographic inference. Species-tree analyses grouped Polioptilidae into four main clades: Microbates, Ramphocaenus, a Polioptila guianensis complex, and the remaining members of Polioptila. Ramphocaenus melanurus was sister to the clade containing M. cinereiventris and M. collaris, which formed a clade sister to all species within Polioptila. Polioptila was composed of two clades, the first of which included the P. guianensis complex; the other contained all remaining species in the genus. Using multispecies coalescent modeling, we inferred a more than 3-fold increase in species diversity, of which 87% represent currently recognized species or subspecies. Much of this diversity corresponded to subspecies that occur in the Neotropics. We identified three polyphyletic species, and delimited 4-6 previously undescribed candidate taxa. Probabilistic modeling of geographic ranges on the species tree indicated that the family likely had an ancestral origin in South America, with all three genera independently colonizing North America. Support for this hypothesis, however, was sensitive to the taxonomic classification scheme used and the number of geographical areas allowed. Our study proposes the first phylogenetic hypothesis for Polioptilidae and provides genealogical support for the reclassification of species limits. Species limits and the resolution of geographical areas that taxa inhabit influence the inferred spatial diversification history.

Phylogenomic insights into the diversification of salamanders in the Isthmura bellii group across the Mexican highlands.

Mountain formation in Mexico has played an important role in the diversification of many Mexican taxa. The Trans-Mexican Volcanic Belt in particular has served as both a cradle of diversification and conduit for dispersal. We investigated the evolutionary history of the Isthmura bellii group of salamanders, a widespread amphibian across the Mexican highlands, using sequence capture of ultraconserved elements. Results suggest that the I. bellii group probably originated in southeastern Mexico in the late Miocene and later dispersed across the Trans-Mexican Volcanic Belt and into the Sierra Madre Occidental. Pre-Pleistocene uplift of the Trans-Volcanic Belt likely promoted early diversification by serving as a mesic land-bridge across central Mexico. These findings highlight the importance of the Trans-Volcanic Belt in generating Mexico's rich biodiversity.

Venom phenotypes of the Rock Rattlesnake (Crotalus lepidus) and the Ridge-nosed Rattlesnake (Crotalus willardi) from México and the United States.

Although the Mexican Highlands has the highest diversity of small-bodied rattlesnakes in the world, studies on the species found throughout this region have been relatively scarce. This has led to challenges with examining venom phenotypic characteristics, as well as species misidentifications and misclassifications. In the current study we investigated venom variation among four subspecies of Crotalus lepidus (C. l. klaluberi, C. l. lepidus, C. l. maculosus, C. l. morulus) and four subspecies of C. willardi (C. w. amabilis, C. w. obscurus, C. w. silus, and C. w. willardi) that inhabit regions of southwestern United States and central México. SDS-PAGE patterns show the presence of many of the major compounds found in other rattlesnake venoms, although minor variations in protein banding patterns and intensity are recognizable. Most notably, PI-metalloproteinase (SVMP) bands appear to be very faint to absent in northern C. l. lepidus and C. l. klauberi subspecies, but are fairly prominent in all other C. lepidus and C. willardi subspecies. Enzyme activity assays revealed that C. lepidus subspecies exhibit higher SVMP and thrombin-like activities when compared to C. willardi subspecies. Significant differences between subspecies were also observed for kallikrein-like serine protease, L-amino acid oxidase, and phosphodiesterase activities, although these differences appear to be random and fail to follow a geographical or phylogenetic trend. The same relationship was also observed for fibrinogenolytic and coagulation assays. Toxicity assays conducted on lab mice (Mus musculus), house geckos (Hemidactylus frenatus), and house crickets (Acheta domestica) revealed varying toxicities between subspecies, with C. l klauberi being the most toxic towards mice (LD50 = 1.36 µg/g) and house geckos (LD50 = 0.17 µg/g), and C. w. silus being most toxic to house crickets (LD50 = 1.94 µg/g). These results provide additional evidence that geographical isolation, natural selection, and adaptive evolution in response to diets may be driving forces contributing to population-level variation in venom composition.

Cryptic diversity and discordance in single-locus species delimitation methods within horned lizards (Phrynosomatidae: Phrynosoma).

Biodiversity reduction and loss continues to progress at an alarming rate, and thus, there is widespread interest in utilizing rapid and efficient methods for quantifying and delimiting taxonomic diversity. Single-locus species delimitation methods have become popular, in part due to the adoption of the DNA barcoding paradigm. These techniques can be broadly classified into tree-based and distance-based methods depending on whether species are delimited based on a constructed genealogy. Although the relative performance of these methods has been tested repeatedly with simulations, additional studies are needed to assess congruence with empirical data. We compiled a large data set of mitochondrial ND4 sequences from horned lizards (Phrynosoma) to elucidate congruence using four tree-based (single-threshold GMYC, multiple-threshold GMYC, bPTP, mPTP) and one distance-based (ABGD) species delimitation models. We were particularly interested in cases with highly uneven sampling and/or large differences in intraspecific diversity. Results showed a high degree of discordance among methods, with multiple-threshold GMYC and bPTP suggesting an unrealistically high number of species (29 and 26 species within the P. douglasii complex alone). The single-threshold GMYC model was the most conservative, likely a result of difficulty in locating the inflection point in the genealogies. mPTP and ABGD appeared to be the most stable across sampling regimes and suggested the presence of additional cryptic species that warrant further investigation. These results suggest that the mPTP model may be preferable in empirical data sets with highly uneven sampling or large differences in effective population sizes of species.

High phylogenetic utility of an ultraconserved element probe set designed for Arachnida.

Arachnida is an ancient, diverse and ecologically important animal group that contains a number of species of interest for medical, agricultural and engineering applications. Despite their importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative approaches to arachnid biology. Biologists have made considerable efforts to resolve the arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of genetic resources, have hindered progress. Here, we present a genomic toolkit for arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that allow targeted sequencing of any species in the arachnid tree of life. We used recently developed capture probes designed from conserved regions of available arachnid genomes to enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs produced a highly resolved arachnid tree with relationships largely consistent with recent transcriptome-based phylogenies. We also tested the phylogenetic informativeness of UCE probes within the spider, scorpion and harvestman orders, demonstrating the utility of these markers at shallower taxonomic scales and suggesting that these loci will be useful for species-level differences. This probe set will open the door to phylogenomic and population genomic studies across the arachnid tree of life, enabling systematics, species delimitation, species discovery and conservation of these diverse arthropods.

Hidden histories of gene flow in highland birds revealed with genomic markers.

Genomic studies are revealing that divergence and speciation are marked by gene flow, but it is not clear whether gene flow has played a prominent role during the generation of biodiversity in species-rich regions of the world where vicariance is assumed to be the principal mode by which new species form. We revisit a well-studied organismal system in the Mexican Highlands, Aphelocoma jays, to test for gene flow among Mexican sierras. Prior results from mitochondrial DNA (mtDNA) largely conformed to the standard model of allopatric divergence, although there was also evidence for more obscure histories of gene flow in a small sample of nuclear markers. We tested for these 'hidden histories' using genomic markers known as ultraconserved elements (UCEs) in concert with phylogenies, clustering algorithms and newer introgression tests specifically designed to detect ancient gene flow (e.g. ABBA/BABA tests). Results based on 4303 UCE loci and 2500 informative SNPs are consistent with varying degrees of gene flow among highland areas. In some cases, gene flow has been extensive and recent (although perhaps not ongoing today), whereas in other cases there is only a trace signature of ancient gene flow among species that diverged as long as 5 million years ago. These results show how a species complex thought to be a model for vicariance can reveal a more reticulate history when a broader portion of the genome is queried. As more organisms are studied with genomic data, we predict that speciation-with-bouts-of-gene-flow will turn out to be a common mode of speciation.

Genomic data reveal ancient microendemism in forest scorpions across the California Floristic Province.

The California Floristic Province (CFP) in western North America is a globally significant biodiversity hotspot. Elucidating patterns of endemism and the historical drivers of this diversity has been an important challenge of comparative phylogeography for over two decades. We generated phylogenomic data using ddRADseq to examine genetic structure in Uroctonus forest scorpions, an ecologically restricted and dispersal-limited organism widely distributed across the CFP north to the Columbia River. We coupled our genetic data with species distribution models (SDMs) to determine climatically suitable areas for Uroctonus both now and during the Last Glacial Maximum. Based on our analyses, Uroctonus is composed of two major genetic groups that likely diverged over 2 million years ago. Each of these groups itself contains numerous genetic groups that reveal a pattern of vicariance and microendemism across the CFP. Migration rates among these populations are low. SDMs suggest forest scorpion habitat has remained relatively stable over the last 21 000 years, consistent with the genetic data. Our results suggest tectonic plate rafting, mountain uplift, river drainage formation and climate-induced habitat fragmentation have all likely played a role in the diversification of Uroctonus. The intricate pattern of genetic fragmentation revealed across a temporal continuum highlights the potential of low-dispersing species to shed light on small-scale patterns of biodiversity and the underlying processes that have generated this diversity in biodiversity hotspots.

Pseudouroctonus maidu, a new species of scorpion from northern California (Scorpiones, Vaejovidae).

A new species of vaejovid scorpion from northern California, Pseudouroctonus maidu sp. n., is named and described. This new species appears to be most similar to Pseudouroctonus iviei (Gertsch & Soleglad, 1972) and Pseudouroctonus glimmei (Hjelle, 1972).

A new species of the "mexicanus" group of the genus Vaejovis C. L. Koch, 1836 from the Mexican state of Aguascalientes (Scorpiones: Vaejovidae).

A new species of Vaejovis is described from the Mexican state of Aguascalientes. It is assigned to the "mexicanus" group and compared with similar species from Jalisco, Guanajuato, and San Luis Potosí. A map with their known distributions is provided.

A multilocus phylogeny of a major New World avian radiation: the Vireonidae.

The family Vireonidae represents one of the most widespread and well-known New World avian radiations, but a robust species-level phylogeny of the group is lacking. Here, we infer a phylogeny of Vireonidae using multilocus data obtained from 221 individuals from 46 of 52 vireonid species (representing all four genera) and five "core Corvoidea" outgroups. Our results show Vireonidae to be monophyletic, consistent with a single colonization of the New World by an Asian ancestor. Cyclarhis and Vireolanius are monophyletic genera that diverged early from the rest of Vireonidae. Hylophilus is polyphyletic, represented by three distinct clades concordant with differences in morphology, habitat, and voice. The poorly known South American species Hylophilus sclateri is embedded within the genus Vireo. Vireo, in turn, consists of several well-supported intrageneric clades. Overall, tropical vireonid species show much higher levels of intraspecific genetic structure than temperate species and several currently recognized species are probably comprised of multiple cryptic species.

Multilocus species delimitation in the Crotalus triseriatus species group (Serpentes: Viperidae: Crotalinae), with the description of two new species.

Members of the Crotalus triseriatus species group of montane rattlesnakes are widely distributed across the highlands of Mexico and southwestern USA. Although five species are currently recognized within the group, species limits remain to be tested. Genetic studies suggest that species may be paraphyletic and that at least one cryptic species may be present. We generate 3,346 base pairs of DNA sequence data from seven nuclear loci to test competing models of species delimitation in the C. triseriatus group using Bayes factor delimitation. We also examine museum specimens from the Trans-Mexican Volcanic Belt for evidence of cryptic species. We find strong support for a nine-species model and genetic and morphological evidence for recognizing two new species within the group, which we formally describe here. Our results suggest that the current taxonomy of the C. triseriatus species group does not reflect evolutionary history. We suggest several conservative taxonomic changes to the group, but future studies are needed to better clarify relationships among species and examine genetic patterns and structure within wide-ranging lineages.