Reference genome of the California glossy snake, Arizona elegans occidentalis: A declining California Species of Special Concern.

The glossy snake (Arizona elegans) is a polytypic species broadly distributed across southwestern North America. The species occupies habitats ranging from California's coastal chaparral to the shortgrass prairies of Texas and southeastern Nebraska, to the extensive arid scrublands of central México. Three subspecies are currently recognized in California, one of which is afforded state-level protection based on the extensive loss and modification of its preferred alluvial coastal scrub and inland desert habitat. We report the first genome assembly of A. elegans occidentalis as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technologies to produce a de novo assembled genome. The assembly comprises a total of 140 scaffolds spanning 1,842,602,218 base pairs, has a contig NG50 of 61 Mb, a scaffold NG50 of 136 Mb, and a BUSCO complete score of 95.9%, and is one of the most complete snake genome assemblies. The A. e. occidentalis genome will be a key tool for understanding the genomic diversity and the basis of adaptations within this species and close relatives within the hyperdiverse snake family Colubridae.

Subspecies differentiation and range-wide genetic structure are driven by climate in the California gnatcatcher, a flagship species for coastal sage scrub conservation.

Understanding genetic structure and diversity within species can uncover associations with environmental and geographic attributes that highlight adaptive potential and inform conservation and management. The California gnatcatcher, Polioptila californica, is a small songbird found in desert and coastal scrub habitats from the southern end of Baja California Sur to Ventura County, California. Lack of congruence among morphological subspecies hypotheses and lack of measurable genetic structure found in a few genetic markers led to questions about the validity of subspecies within P. californica and the listing status of the coastal California gnatcatcher, P. c. californica. As a U.S. federally threatened subspecies, P. c. californica is recognized as a flagship for coastal sage scrub conservation throughout southern California. We used restriction site-associated DNA sequencing to develop a genomic dataset for the California gnatcatcher. We sampled throughout the species' range, examined genetic structure, gene-environment associations, and demographic history, and tested for concordance between genetic structure and morphological subspecies groups. Our data support two distinct genetic groups with evidence of restricted movement and gene flow near the U.S.- Mexico international border. We found that climate-associated outlier loci were more strongly differentiated than climate neutral loci, suggesting that local climate adaptation may have helped to drive differentiation after Holocene range expansions. Patterns of habitat loss and fragmentation are also concordant with genetic substructure throughout the southern California portion of the range. Finally, our genetic data supported the morphologically defined P. c. californica as a distinct group, but there was little evidence of genetic differentiation among other previously hypothesized subspecies in Baja California. Our data suggest that retaining and restoring connectivity, and protecting populations, particularly at the northern range edge, could help preserve existing adaptive potential to allow for future range expansion and long-term persistence of the California gnatcatcher.

Correction: Combining genetic and demographic monitoring better informs conservation of an endangered urban snake.

[This corrects the article DOI: 10.1371/journal.pone.0231744.].

Combining genetic and demographic monitoring better informs conservation of an endangered urban snake.

Conversion and fragmentation of wildlife habitat often leads to smaller and isolated populations and can reduce a species' ability to disperse across the landscape. As a consequence, genetic drift can quickly lower genetic variation and increase vulnerability to extirpation. For species of conservation concern, quantification of population size and connectivity can clarify the influence of genetic drift in local populations and provides important information for conservation management and recovery strategies. Here, we used genome-wide single nucleotide polymorphism (SNP) data and capture-mark-recapture methods to evaluate the genetic diversity and demography within seven focal sites of the endangered San Francisco gartersnake (Thamnophis sirtalis tetrataenia), a species affected by alteration and isolation of wetland habitats throughout its distribution. The primary goals were to determine the population structure and degree of genetic isolation among T. s. tetrataenia populations and estimate effective size and population abundance within sites to better understand the present and future importance of genetic drift. We also used temporally sampled datasets to examine the magnitude of genetic change over time. We found moderate population genetic structure throughout the San Francisco Peninsula that partitions sites into northern and southern regional clusters. Point estimates of both effective size and population abundance were generally small (≤ 100) for a majority of the sites, and estimates were particularly low in the northern populations. Genetic analyses of temporal datasets indicated an increase in genetic differentiation, especially for the most geographically isolated sites, and decreased genetic diversity over time in at least one site (Pacifica). Our results suggest that drift-mediated processes as a function of small population size and reduced connectivity from neighboring populations may decrease diversity and increase differentiation. Improving genetic diversity and connectivity among T. s. tetrataenia populations could promote persistence of this endangered snake.

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.

Persistence of historical population structure in an endangered species despite near-complete biome conversion in California's San Joaquin Desert.

Genomic responses to habitat conversion can be rapid, providing wildlife managers with time-limited opportunities to enact recovery efforts that use population connectivity information that reflects predisturbance landscapes. Despite near-complete biome conversion, such opportunities may still exist for the endemic fauna and flora of California's San Joaquin Desert, but comprehensive genetic data sets are lacking for nearly all species in the region. To fill this knowledge gap, we studied the rangewide population structure of the endangered blunt-nosed leopard lizard Gambelia sila, a San Joaquin Desert endemic, using restriction site-associated DNA (RAD), microsatellite and mtDNA data to test whether admixture patterns and estimates of effective migration surfaces (EEMS) can identify land areas with high population connectivity prior to the conversion of native xeric habitats. Clustering and phylogenetic analyses indicate a recent shared history between numerous isolated populations and EEMS reveals latent signals of corridors and barriers to gene flow over areas now replaced by agriculture and urbanization. Conflicting histories between the mtDNA and nuclear genomes are consistent with hybridization with the sister species G. wislizenii, raising important questions about where legal protection should end at the southern range limit of G. sila. Comparative analysis of different data sets also adds to a growing list of advantages in using RAD loci for genetic studies of rare species. We demonstrate how the results of this work can serve as an evolutionary guidance tool for managing endemic, arid-adapted taxa in one of the world's most compromised landscapes.

Lineage diversification of fringe-toed lizards (Phrynosomatidae: Uma notata complex) in the Colorado Desert: Delimiting species in the presence of gene flow.

Multi-locus nuclear DNA data were used to delimit species of fringe-toed lizards of the Uma notata complex, which are specialized for living in wind-blown sand habitats in the deserts of southwestern North America, and to infer whether Quaternary glacial cycles or Tertiary geological events were important in shaping the historical biogeography of this group. We analyzed ten nuclear loci collected using Sanger sequencing and genome-wide sequence/single-nucleotide polymorphism (SNP) data collected using restriction-associated DNA (RAD) sequencing. A combination of species discovery methods (concatenated phylogenies, parametric and non-parametric clustering algorithms) and species validation approaches (coalescent-based species tree/isolation-with-migration models) were used to delimit species, infer phylogenetic relationships, and to estimate effective population sizes, migration rates, and speciation times. Uma notata, U. inornata, U. cowlesi, and an undescribed species from Mohawk Dunes, Arizona (U. sp.) were supported as distinct in the concatenated analyses and by clustering algorithms, and all operational taxonomic units were decisively supported as distinct species by ranking hierarchical nested speciation models with Bayes factors based on coalescent-based species tree methods. However, significant unidirectional gene flow (2NM>1) from U. cowlesi and U. notata into U. rufopunctata was detected under the isolation-with-migration model. Therefore, we conservatively delimit four species-level lineages within this complex (U. inornata, U. notata, U. cowlesi, and U. sp.), treating U. rufopunctata as a hybrid population (U. notata×cowlesi). Both concatenated and coalescent-based estimates of speciation times support the hypotheses that speciation within the complex occurred during the late Pleistocene, and that the geological evolution of the Colorado River delta during this period was an important process shaping the observed phylogeographic patterns.

Fuzzy boundaries: color and gene flow patterns among parapatric lineages of the western shovel-nosed snake and taxonomic implication.

Accurate delineation of lineage diversity is increasingly important, as species distributions are becoming more reduced and threatened. During the last century, the subspecies category was often used to denote phenotypic variation within a species range and to provide a framework for understanding lineage differentiation, often considered incipient speciation. While this category has largely fallen into disuse, previously recognized subspecies often serve as important units for conservation policy and management when other information is lacking. In this study, we evaluated phenotypic subspecies hypotheses within shovel-nosed snakes on the basis of genetic data and considered how evolutionary processes such as gene flow influenced possible incongruence between phenotypic and genetic patterns. We used both traditional phylogenetic and Bayesian clustering analyses to infer range-wide genetic structure and spatially explicit analyses to detect possible boundary locations of lineage contact. Multilocus analyses supported three historically isolated groups with low to moderate levels of contemporary gene exchange. Genetic data did not support phenotypic subspecies as exclusive groups, and we detected patterns of discordance in areas where three subspecies are presumed to be in contact. Based on genetic and phenotypic evidence, we suggested that species-level diversity is underestimated in this group and we proposed that two species be recognized, Chionactis occipitalis and C. annulata. In addition, we recommend retention of two subspecific designations within C. annulata (C. a. annulata and C. a. klauberi) that reflect regional shifts in both genetic and phenotypic variation within the species. Our results highlight the difficultly in validating taxonomic boundaries within lineages that are evolving under a time-dependent, continuous process.

Refugial isolation and divergence in the Narrowheaded Gartersnake species complex (Thamnophis rufipunctatus) as revealed by multilocus DNA sequence data.

Glacial-interglacial cycles of the Pleistocene are hypothesized as one of the foremost contributors to biological diversification. This is especially true for cold-adapted montane species, where range shifts have had a pronounced effect on population-level divergence. Gartersnakes of the Thamnophis rufipunctatus species complex are restricted to cold headwater streams in the highlands of the Sierra Madre Occidental and southwestern USA. We used coalescent and multilocus phylogenetic approaches to test whether genetic diversification of this montane-restricted species complex is consistent with two prevailing models of range fluctuation for species affected by Pleistocene climate changes. Our concatenated nuDNA and multilocus species analyses recovered evidence for the persistence of multiple lineages that are restricted geographically, despite a mtDNA signature consistent with either more recent connectivity (and introgression) or recent expansion (and incomplete lineage sorting). Divergence times estimated using a relaxed molecular clock and fossil calibrations fall within the Late Pleistocene, and zero gene flow scenarios among current geographically isolated lineages could not be rejected. These results suggest that increased climate shifts in the Late Pleistocene have driven diversification and current range retraction patterns and that the differences between markers reflect the stochasticity of gene lineages (i.e. ancestral polymorphism) rather than gene flow and introgression. These results have important implications for the conservation of T. rufipunctatus (sensu novo), which is restricted to two drainage systems in the southwestern US and has undergone a recent and dramatic decline.

Novel patterns of historical isolation, dispersal, and secondary contact across Baja California in the Rosy Boa (Lichanura trivirgata).

Mitochondrial DNA (mtDNA) sequence variation was examined in 131 individuals of the Rosy Boa (Lichanura trivirgata) from across the species range in southwestern North America. Bayesian inference and nested clade phylogeographic analyses (NCPA) were used to estimate relationships and infer evolutionary processes. These patterns were evaluated as they relate to previously hypothesized vicariant events and new insights are provided into the biogeographic and evolutionary processes important in Baja California and surrounding North American deserts. Three major lineages (Lineages A, B, and C) are revealed with very little overlap. Lineage A and B are predominately separated along the Colorado River and are found primarily within California and Arizona (respectively), while Lineage C consists of disjunct groups distributed along the Baja California peninsula as well as south-central Arizona, southward along the coastal regions of Sonora, Mexico. Estimated divergence time points (using a Bayesian relaxed molecular clock) and geographic congruence with postulated vicariant events suggest early extensions of the Gulf of California and subsequent development of the Colorado River during the Late Miocene-Pliocene led to the formation of these mtDNA lineages. Our results also suggest that vicariance hypotheses alone do not fully explain patterns of genetic variation. Therefore, we highlight the importance of dispersal to explain these patterns and current distribution of populations. We also compare the mtDNA lineages with those based on morphological variation and evaluate their implications for taxonomy.