Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression.

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Landscape genetics of a sub-alpine toad: climate change predicted to induce upward range shifts via asymmetrical migration corridors.

Climate change is expected to have a major hydrological impact on the core breeding habitat and migration corridors of many amphibians in the twenty-first century. The Yosemite toad (Anaxyrus canorus) is a species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. Despite living entirely on federal lands, it has recently faced severe extirpations, yet our understanding of climatic influences on population connectivity is limited. In this study, we used a previously published double-digest RADseq dataset along with numerous remotely sensed habitat features in a landscape genetics framework to answer two primary questions in Yosemite National Park: (1) Which fine-scale climate, topographic, soil, and vegetation features most facilitate meadow connectivity? (2) How is climate change predicted to influence both the magnitude and net asymmetry of genetic migration? We developed an approach for simultaneously modeling multiple toad migration paths, akin to circuit theory, except raw environmental features can be separately considered. Our workflow identified the most likely migration corridors between meadows and used the unique cubist machine learning approach to fit and forecast environmental models of connectivity. We identified the permuted modeling importance of numerous snowpack-related features, such as runoff and groundwater recharge. Our results highlight the importance of considering phylogeographic structure, and asymmetrical migration in landscape genetics. We predict an upward elevational shift for this already high-elevation species, as measured by the net vector of anticipated genetic movement, and a north-eastward shift in species distribution via the network of genetic migration corridors across the park.

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.

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.

Habitat fragmentation in coastal southern California disrupts genetic connectivity in the cactus wren (Campylorhynchus brunneicapillus).

Achieving long-term persistence of species in urbanized landscapes requires characterizing population genetic structure to understand and manage the effects of anthropogenic disturbance on connectivity. Urbanization over the past century in coastal southern California has caused both precipitous loss of coastal sage scrub habitat and declines in populations of the cactus wren (Campylorhynchus brunneicapillus). Using 22 microsatellite loci, we found that remnant cactus wren aggregations in coastal southern California comprised 20 populations based on strict exact tests for population differentiation, and 12 genetic clusters with hierarchical Bayesian clustering analyses. Genetic structure patterns largely mirrored underlying habitat availability, with cluster and population boundaries coinciding with fragmentation caused primarily by urbanization. Using a habitat model we developed, we detected stronger associations between habitat-based distances and genetic distances than Euclidean geographic distance. Within populations, we detected a positive association between available local habitat and allelic richness and a negative association with relatedness. Isolation-by-distance patterns varied over the study area, which we attribute to temporal differences in anthropogenic landscape development. We also found that genetic bottleneck signals were associated with wildfire frequency. These results indicate that habitat fragmentation and alterations have reduced genetic connectivity and diversity of cactus wren populations in coastal southern California. Management efforts focused on improving connectivity among remaining populations may help to ensure population persistence.

Best practices for genetic and genomic data archiving.

Genetic and genomic data are collected for a vast array of scientific and applied purposes. Despite mandates for public archiving, data are typically used only by the generating authors. The reuse of genetic and genomic datasets remains uncommon because it is difficult, if not impossible, due to non-standard archiving practices and lack of contextual metadata. But as the new field of macrogenetics is demonstrating, if genetic data and their metadata were more accessible and FAIR (findable, accessible, interoperable and reusable) compliant, they could be reused for many additional purposes. We discuss the main challenges with existing genetic and genomic data archives, and suggest best practices for archiving genetic and genomic data. Recognizing that this is a longstanding issue due to little formal data management training within the fields of ecology and evolution, we highlight steps that research institutions and publishers could take to improve data archiving.

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus).

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future "Geminate Evolutionary Units" (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad (Anaxyrus canorus), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection (R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

Move it or lose it: Predicted effects of culverts and population density on Mojave desert tortoise (Gopherus agassizii) connectivity.

Roadways and railways can reduce wildlife movements across landscapes, negatively impacting population connectivity. Connectivity may be improved by structures that allow safe passage across linear barriers, but connectivity could be adversely influenced by low population densities. The Mojave desert tortoise is threatened by habitat loss, fragmentation, and population declines. The tortoise continues to decline as disturbance increases across the Mojave Desert in the southwestern United States. While underground crossing structures, like hydrological culverts, have begun receiving attention, population density has not been considered in tortoise connectivity. Our work asks a novel question: How do culverts and population density affect connectivity and potentially drive genetic and demographic patterns? To explore the role of culverts and population density, we used agent-based spatially explicit forward-in-time simulations of gene flow. We constructed resistance surfaces with a range of barriers to movement and representative of tortoise habitat with anthropogenic disturbance. We predicted connectivity under variable population densities. Simulations were run for 200 non-overlapping generations (3400 years) with 30 replicates using 20 microsatellite loci. We evaluated population genetic structure and diversity and found that culverts would not entirely negate the effects of linear barriers, but gene flow improved. Our results also indicated that density is important for connectivity. Low densities resulted in declines regardless of the landscape barrier scenario (> 75% population census size, > 97% effective population size). Results from our simulation using current anthropogenic disturbance predicted decreased population connectivity over time. Genetic and demographic effects were detectable within five generations (85 years) following disturbance with estimated losses in effective population size of 69%. The pronounced declines in effective population size indicate this could be a useful monitoring metric. We suggest management strategies that improve connectivity, such as roadside fencing tied to culverts, conservation areas in a connected network, and development restricted to disturbed areas.

Growth, drought response, and climate-associated genomic structure in whitebark pine in the Sierra Nevada of California.

Whitebark pine (Pinus albicaulis Engelm.) has experienced rapid population declines and is listed as threatened under the Endangered Species Act in the United States. Whitebark pine in the Sierra Nevada of California represents the southernmost end of the species' distribution and, like other portions of its range, faces threats from an introduced pathogen, native bark beetles, and a rapidly warming climate. Beyond these chronic stressors, there is also concern about how this species will respond to acute stressors, such as drought. We present patterns of stem growth from 766 large (average diameter at breast height >25 cm), disease-free whitebark pine across the Sierra Nevada before and during a recent period of drought. We contextualize growth patterns using population genomic diversity and structure from a subset of 327 trees. Sampled whitebark pine generally had positive to neutral stem growth trends from 1970 to 2011, which was positively correlated with minimum temperature and precipitation. Indices of stem growth during drought years (2012 to 2015) relative to a predrought interval were mostly positive to neutral at our sampled sites. Individual tree growth response phenotypes appeared to be linked to genotypic variation in climate-associated loci, suggesting that some genotypes can take better advantage of local climatic conditions than others. We speculate that reduced snowpack during the 2012 to 2015 drought years may have lengthened the growing season while retaining sufficient moisture to maintain growth at most study sites. Growth responses may differ under future warming, however, particularly if drought severity increases and modifies interactions with pests and pathogens.

Fecal metabarcoding of the endangered Pacific pocket mouse (Perognathus longimembris pacificus) reveals a diverse and forb rich diet that reflects local habitat availability.

Information on diet breadth and preference can assist in understanding links between food resources and population growth and inform habitat restoration for rare herbivores. We assessed the diet of the endangered Pacific pocket mouse using metabarcoding of fecal samples and compared it to plant community composition in long-term study plots in two populations on Marine Corps Base Camp Pendleton, San Diego County, CA. Fecal samples (n = 221) were collected between spring 2016 and fall 2017 during monthly live-trap surveys. Concurrently, percent cover and plant phenology were measured in plots centered on trap locations. Fecal samples were sequenced with paired-end reads of the internal transcribed spacer 2 region of the nuclear ribosomal gene, and the resulting amplicons were matched to a regionally specific database. Seventy-three plant taxa were detected, which were mostly forbs and perennial herbs (70-90%). Diet composition differed between populations, years, seasons, and plots. Overall, diet and local habitat composition in plots were significantly correlated. However, we detected some differences in above-ground seed availability and proportion in fecal samples that indicate diet preferences for some forbs, perennial herbs, and native bunch grasses over perennial shrubs and non-native grasses. This is the first study of PPM to pair plant phenology surveys with diet metabarcoding to estimate resource selection, and results suggest that managing habitat for diverse native forb communities and reducing non-native grass cover may be beneficial for this critically endangered species.

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.

Generic relationships of New World Jerusalem crickets (Orthoptera: Stenopelmatoidea:Stenopelmatinae), including all known species of Stenopelmatus.

The New World Jerusalem crickets currently consist of 4 genera: Stenopelmatus Burmeister, 1838, with 33 named entities; Ammopelmatus Tinkham, 1965, with 2 described species; Viscainopelmatus Tinkham, 1970, with 1 described species, and Stenopelmatopterus Gorochov, 1988, with 3 described species. We redefine the generic boundaries of these 4 genera, synonymize Stenopelmatopterus under Stenopelmatus, and synonymize Viscainopelmatus under Ammopelmatus. We then discuss, and illustrate, all the types of the species of Stenopelmatus, all of which only occur south of the United States' border.                We recognize as valid the following 5 described Mexican and Central American species: S. ater, S. piceiventris, S. sartorianus, S. talpa, and S. typhlops. We declare the following 13 described Mexican and Central American Stenopelmatus taxa as nomen dubium: S. calcaratus, S. erythromelus, S. guatemalae, S. histrio, S. lessonae, S. lycosoides, S. mexicanus, S. minor, S. nieti, S. sallei, S. sumichrasti, S. toltecus, and S. vicinus. We designate a neotype for S. talpa and lectotypes for S. ater, S. guatemalae, S. histrio, S. lessonae, S. mexicanus, S. minor, S. nieti, S. sallei, S. sumichrasti, and S. toltecus. We assign a type locality for S. piceiventris. We concur with the previous synonymy of S. politus under S. sartorianus. We describe 14 new species of Stenopelmatus from Mexico, Honduras and Ecuador, based on a combination of adult morphology, DNA, calling song drumming pattern, distribution, and karyotype: S. chiapas sp. nov., S. cusuco sp. nov., S. diezmilpies sp. nov., S. durango sp. nov., S. ecuadorensis sp. nov., S. faulkneri sp. nov., S. honduras sp. nov., S. hondurasito sp. nov., S. mineraldelmonte sp. nov., S. nuevoleon sp. nov., S. perote sp. nov., S. saltillo sp. nov., S. sanfelipe sp. nov., and S. zimapan sp. nov.                  We transfer the following 16 described United States taxa, plus S. cephalotes from the "west coast of North America", from Stenopelmatus to Ammopelmatus: A. cahuilaensis, A. californicus, A. cephalotes, A. fasciatus, A. fuscus, A. hydrocephalus, A. intermedius, A. irregularis, A. longispinus, A. mescaleroensis, A. monahansensis, A. navajo, A. nigrocapitatus, A. oculatus, A. pictus, and A. terrenus, along with the Mexican taxon A. comanchus: these species will be discussed in a subsequent paper (Weissman et al. in prep).                We believe that all new Jerusalem cricket species descriptions should include, at a minimum, calling drum (most important) and DNA information.

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.

Using movement to inform conservation corridor design for Mojave desert tortoise.

BACKGROUND: Preserving corridors for movement and gene flow among populations can assist in the recovery of threatened and endangered species. As human activity continues to fragment habitats, characterizing natural corridors is important in establishing and maintaining connectivity corridors within the anthropogenic development matrix. The Mojave desert tortoise (Gopherus agassizii) is a threatened species occupying a variety of habitats in the Mojave and Colorado Deserts. Desert tortoises have been referred to as corridor-dwellers, and understanding how they move within suitable habitat can be crucial to defining corridors that will sustain sufficient gene flow to maintain connections among populations amidst the increases in human development. METHODS: To elucidate how tortoises traverse available habitat and interact with potentially inhospitable terrain and human infrastructure, we used GPS dataloggers to document fine-scale movement of individuals and estimate home ranges at ten study sites along the California/Nevada border. Our sites encompass a variety of habitats, including mountain passes that serve as important natural corridors connecting neighboring valleys, and are impacted by a variety of linear anthropogenic features. We used path selection functions to quantify tortoise movements and develop resistance surfaces based on landscape characteristics including natural features, anthropogenic alterations, and estimated home ranges with autocorrelated kernel density methods. Using the best supported path selection models and estimated home ranges, we determined characteristics of known natural corridors and compared them to mitigation corridors (remnant habitat patches) that have been integrated into land management decisions in the Ivanpah Valley. RESULTS: Tortoises avoided areas of high slope and low perennial vegetation cover, avoided moving near low-density roads, and traveled along linear barriers (fences and flood control berms). CONCLUSIONS: We found that mitigation corridors designated between solar facilities should be wide enough to retain home ranges and maintain function. Differences in home range size and movement resistance between our two natural mountain pass corridors align with differences in genetic connectivity, suggesting that not all natural corridors provide the same functionality. Furthermore, creation of mitigation corridors with fences may have unintended consequences and may function differently than natural corridors. Understanding characteristics of corridors with different functionality will help future managers ensure that connectivity is maintained among Mojave desert tortoise populations.

Pleistocene glacial cycles drove lineage diversification and fusion in the Yosemite toad (Anaxyrus canorus).

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this "glacial pulse" model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed ("fused") lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double-digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low- and high-elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.

Distinguishing recent dispersal from historical genetic connectivity in the coastal California gnatcatcher.

Habitat loss and fragmentation are primary threats to biodiversity worldwide. We studied the impacts of habitat loss and fragmentation on genetic connectivity and diversity among local aggregations of the California gnatcatcher (Polioptila californica californica) across its U.S. range. With a dataset of 268 individuals genotyped at 19 microsatellite loci, we analyzed genetic structure across the range using clustering analyses, exact tests for population differentiation, and a pedigree analysis to examine the spatial distribution of first-order relatives throughout the study area. In addition, we developed a habitat suitability model and related percent suitable habitat to genetic diversity indices within aggregations at two spatial scales. We detected a single genetic cluster across the range, with weak genetic structure among recently geographically isolated aggregations in the northern part of the range. The pedigree analysis detected closely related individuals across disparate aggregations and across large geographic distances in the majority of the sampled range, demonstrating that recent long-distance dispersal has occurred within this species. Genetic diversity was independent of suitable habitat at a local 5-km scale, but increased in a non-linear fashion with habitat availability at a broader, 30-km scale. Diversity declined steeply when suitable habitat within 30-km fell below 10%. Together, our results suggest that California gnatcatchers retain genetic connectivity across the majority of the current distribution of coastal sage scrub fragments, with the exception of some outlying aggregations. Connectivity may help support long-term persistence under current conservation and management strategies. However, emerging structure among more remote aggregations and associations between available habitat and genetic diversity also suggest that continued loss of habitat could threaten diversity and connectivity in the future.

Prioritizing conserved areas threatened by wildfire and fragmentation for monitoring and management.

In many parts of the world, the combined effects of habitat fragmentation and altered disturbance regimes pose a significant threat to biodiversity. This is particularly true in Mediterranean-type ecosystems (MTEs), which tend to be fire-prone, species rich, and heavily impacted by human land use. Given the spatial complexity of overlapping threats and species' vulnerability along with limited conservation budgets, methods are needed for prioritizing areas for monitoring and management in these regions. We developed a multi-criteria Pareto ranking methodology for prioritizing spatial units for conservation and applied it to fire threat, habitat fragmentation threat, species richness, and genetic biodiversity criteria in San Diego County, California, USA. We summarized the criteria and Pareto ranking results (from west to east) within the maritime, coastal, transitional, inland climate zones within San Diego County. Fire threat increased from the maritime zone eastward to the transitional zone, then decreased in the mountainous inland climate zone. Number of fires and fire return interval departure were strongly negatively correlated. Fragmentation threats, particularly road density and development density, were highest in the maritime climate zone, declined towards the east, and were positively correlated. Species richness criteria showed distributions among climate zones similar to those of the fire threat variables. When using species richness and fire threat criteria, most lower-ranked (higher conservation priority) units occurred in the coastal and transitional zones. When considering genetic biodiversity, lower-ranked units occurred more often in the mountainous inland zone. With Pareto ranking, there is no need to select criteria weights as part of the decision-making process. However, negative correlations and larger numbers of criteria can result in more units assigned to the same rank. Pareto ranking is broadly applicable and can be used as a standalone decision analysis method or in conjunction with other methods.