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.

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.

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.

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.

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.

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.

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.

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.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

<i>Corrigendum to:</i> Extensive species diversification and marked geographic phylogenetic structure in the Mesoamerican genus <i>Stenopelmatus</i> (Orthoptera: Stenopelmatidae: Stenopelmatinae) revealed by mitochondrial and nuclear 3RAD data

The Jerusalem cricket subfamily Stenopelmatinae is distributed from south-western Canada through the western half of the United States to as far south as Ecuador. Recently, the generic classification of this subfamily was updated to contain two genera, the western North American <i>Ammopelmatus</i>, and the Mexican, and central and northern South American <i>Stenopelmatus</i>. The taxonomy of the latter genus was also revised, with 5, 13 and 14 species being respectively validated, declared as nomen dubium and described as new. Despite this effort, the systematics of <i>Stenopelmatus</i> is still far from complete. Here, we generated sequences of the mitochondrial DNA barcoding locus and performed two distinct DNA sequence-based approaches to assess the species&#x2019; limits among several populations of <i>Stenopelmatus</i>, with emphasis on populations from central and south-east Mexico. We reconstructed the phylogenetic relationships among representative species of the main clades within the genus using nuclear 3RAD data and carried out a molecular clock analysis to investigate its biogeographic history. The two DNA sequence-based approaches consistently recovered 34 putative species, several of which are apparently undescribed. Our estimates of phylogeny confirmed the recent generic update of Stenopelmatinae and revealed a marked phylogeographic structure within <i>Stenopelmatus</i>. Based on our results, we propose the existence of four species-groups within the genus (the <i>faulkneri</i>, <i>talpa</i>, Central America and <i>piceiventris</i> species-groups). The geographic distribution of these species-groups and our molecular clock estimates are congruent with the geological processes that took place in mountain ranges along central and southern Mexico, particularly since the Neogene. Our study emphasises the necessity to continue performing more taxonomic and phylogenetic studies on <i>Stenopelmatus</i> to clarify its actual species richness and evolutionary history in Mesoamerica.

Extensive species diversification and marked geographic phylogenetic structure in the Mesoamerican genus Stenopelmatus (Orthoptera: Stenopelmatidae: Stenopelmatinae) revealed by mitochondrial and nuclear 3RAD data

The Jerusalem cricket subfamily Stenopelmatinae is distributed from south-western Canada through the western half of the United States to as far south as Ecuador. Recently, the generic classification of this subfamily was updated to contain two genera, the western North American Ammopelmatus, and the Mexican, and central and northern South American Stenopelmatus. The taxonomy of the latter genus was also revised, with 5, 13 and 14 species being respectively validated, declared as nomen dubium and described as new. Despite this effort, the systematics of Stenopelmatus is still far from complete. Here, we generated sequences of the mitochondrial DNA barcoding locus and performed two distinct DNA sequence-based approaches to assess the species’ limits among several populations of Stenopelmatus, with emphasis on populations from central and south-east Mexico. We reconstructed the phylogenetic relationships among representative species of the main clades within the genus using nuclear 3RAD data and carried out a molecular clock analysis to investigate its biogeographic history. The two DNA sequence-based approaches consistently recovered 34 putative species, several of which are apparently undescribed. Our estimates of phylogeny confirmed the recent generic update of Stenopelmatinae and revealed a marked phylogeographic structure within Stenopelmatus. Based on our results, we propose the existence of four species-groups within the genus (the faulkneri, talpa, Central America and piceiventris species-groups). The geographic distribution of these species-groups and our molecular clock estimates are congruent with the geological processes that took place in mountain ranges along central and southern Mexico, particularly since the Neogene. Our study emphasises the necessity to continue performing more taxonomic and phylogenetic studies on Stenopelmatus to clarify its actual species richness and evolutionary history in Mesoamerica.

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.

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.

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.

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.

Metabarcoding of Fecal Samples to Determine Herbivore Diets: A Case Study of the Endangered Pacific Pocket Mouse

Understanding the diet of an endangered species illuminates the animal’s ecology, habitat requirements, and conservation needs. However, direct observation of diet can be difficult, particularly for small, nocturnal animals such as the Pacific pocket mouse (Heteromyidae: Perognathus longimembris pacificus). Very little is known of the dietary habits of this federally endangered rodent, hindering management and restoration efforts. We used a metabarcoding approach to identify source plants in fecal samples (N = 52) from the three remaining populations known. The internal transcribed spacers (ITS) of the nuclear ribosomal loci were sequenced following the Illumina MiSeq amplicon strategy and processed reads were mapped to reference databases. We evaluated a range of threshold mapping criteria and found the best-performing setting generally recovered two distinct mock communities in proportions similar to expectation. We tested our method on captive animals fed a known diet and recovered almost all plant sources, but found substantial heterogeneity among fecal pellets collected from the same individual at the same time. Observed richness did not increase with pooling of pellets from the same individual. In field-collected samples, we identified 4–14 plant genera in individual samples and 74 genera overall, but over 50 percent of reads mapped to just six species in five genera. We simulated the effects of sequencing error, variable read length, and chimera formation to infer taxon-specific rates of misassignment for the local flora, which were generally low with some exceptions. Richness at the species and genus levels did not reach a clear asymptote, suggesting that diet breadth remained underestimated in the current pool of samples. Large numbers of scat samples are therefore needed to make inferences about diet and resource selection in future studies of the Pacific pocket mouse. We conclude that our minimally invasive method is promising for determining herbivore diets given a library of sequences from local plants.

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.

Evolutionary dynamics of a rapidly receding southern range boundary in the threatened California Red-Legged Frog (Rana draytonii)

Populations forming the edge of a species range are often imperiled by isolation and low genetic diversity, with proximity to human population centers being a major determinant of edge stability in modern landscapes. Since the 1960s, the California red‐legged frog (Rana draytonii) has undergone extensive declines in heavily urbanized southern California, where the range edge has rapidly contracted northward while shifting its cardinal orientation to an east‐west trending axis. We studied the genetic structure and diversity of these frontline populations, tested for signatures of contemporary disturbance, specifically fire, and attempted to disentangle these signals from demographic events extending deeper into the past. Consistent with the genetic expectations of the ‘abundant‐center’ model, we found that diversity, admixture, and opportunity for random mating increases in populations sampled successively further away from the range boundary. Demographic simulations indicate that bottlenecks in peripheral isolates are associated with processes extending tens to a few hundred generations in the past, despite the demographic collapse of some due to recent fire‐flood events. While the effects of recent disturbance have left little genetic imprint on these populations, they likely contribute to an extinction debt that will lead to continued range contraction unless management intervenes to stall or reverse the process.

Species Differentiation on a Dynamic Landscape: Shifts in Metapopulation Genetic Structure Using the Chronology of the Hawaiian Archipelago

Species formation during adaptive radiation often occurs in the context of a changing environment. The establishment and arrangement of populations, in space and time, sets up ecological and genetic processes that dictate the rate and pattern of differentiation. Here, we focus on how a dynamic habitat can affect genetic structure, and ultimately, differentiation among populations. We make use of the chronology and geographical history provided by the Hawaiian archipelago to examine the initial stages of population establishment and genetic divergence. We use data from a set of 6 spider lineages that differ in habitat affinities, some preferring low elevation habitats with a longer history of connection, others being more specialized for high elevation and/or wet forest, some with more general habitat affinities. We show that habitat preferences associated with lineages are important in ecological and genetic structuring. Lineages that have more restricted habitat preferences are subject to repeated episodes of isolation and fragmentation as a result of lava flows and vegetation succession. The initial dynamic set up by the landscape translates over time into discrete lineages. Further work is needed to understand how genetic changes interact with a changing set of ecological interactions amongst a shifting mosaic of landscapes to achieve species formation.

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.

Identification of co-occurring Branchinecta fairy shrimp species from encysted embryos using multiplex polymerase chain reaction

Morphological identification of many fairy shrimp species is difficult because distinguishing characters are restricted to adults. We developed two multiplex polymerase chain reaction assays that differentiate among three Branchinecta fairy shrimp with distributional overlap in southern California vernal pools. Two of the species are federally listed as threatened. Molecular identification of Branchinecta from cysts allows for species surveys to be conducted during the dry season, expanding the timeframe for population assessment and providing a less intrusive method of sampling sensitive vernal pool habitats.

Understanding the genetic effects of recent habitat fragmentation in the context of evolutionary history: phylogeography and landscape genetics of a southern California endemic Jerusalem cricket (Orthoptera: Stenopelmatidae: Stenopelmatus)

Habitat loss and fragmentation due to urbanization are the most pervasive threats to biodiversity in southern California. Loss of habitat and fragmentation can lower migration rates and genetic connectivity among remaining populations of native species, reducing genetic variability and increasing extinction risk. However, it may be difficult to separate the effects of recent anthropogenic fragmentation from the genetic signature of prehistoric fragmentation due to previous natural geological and climatic changes. To address these challenges, we examined the phylogenetic and population genetic structure of a flightless insect endemic to cismontane southern California, Stenopelmatus'mahogani' (Orthoptera: Stenopelmatidae). Analyses of mitochondrial DNA sequence data suggest that diversification across southern California began during the Pleistocene, with most haplotypes currently restricted to a single population. Patterns of genetic divergence correlate with contemporary urbanization, even after correcting for (geographical information system) GIS-based reconstructions of fragmentation during the Pleistocene. Theoretical simulations confirm that contemporary patterns of genetic structure could be produced by recent urban fragmentation using biologically reasonable assumptions about model parameters. Diversity within populations was positively correlated with current fragment size, but not prehistoric fragment size, suggesting that the effects of increased drift following anthropogenic fragmentation are already being seen. Loss of genetic connectivity and diversity can hinder a population's ability to adapt to ecological perturbations commonly associated with urbanization, such as habitat degradation, climatic changes and introduced species. Consequently, our results underscore the importance of preserving and restoring landscape connectivity for long-term persistence of low vagility native species.

Influence of volcanic activity on the population genetic structure of Hawaiian Tetragnatha spiders: fragmentation, rapid population growth and the potential for accelerated evolution

Volcanic activity on the island of Hawaii results in a cyclical pattern of habitat destruction and fragmentation by lava, followed by habitat regeneration on newly formed substrates. While this pattern has been hypothesized to promote the diversification of Hawaiian lineages, there have been few attempts to link geological processes to measurable changes in population structure. We investigated the genetic structure of three species of Hawaiian spiders in forests fragmented by a 150-year-old lava flow on Mauna Loa Volcano, island of Hawaii: Tetragnatha quasimodo (forest and lava flow generalist), T. anuenue and T. brevignatha (forest specialists). To estimate fragmentation effects on population subdivision in each species, we examined variation in mitochondrial and nuclear genomes (DNA sequences and allozymes, respectively). Population subdivision was higher for forest specialists than for the generalist in fragments separated by lava. Patterns of mtDNA sequence evolution also revealed that forest specialists have undergone rapid expansion, while the generalist has experienced more gradual population growth. Results confirm that patterns of neutral genetic variation reflect patterns of volcanic activity in some Tetragnatha species. Our study further suggests that population subdivision and expansion can occur across small spatial and temporal scales, which may facilitate the rapid spread of new character states, leading to speciation as hypothesized by H. L. Carson 30 years ago.

Mermithid parasitism of Hawaiian Tetragnatha spiders in a fragmented landscape

Hawaiian Tetragnatha spiders inhabiting small forest fragments on the Big Island of Hawaii are parasitized by mermithid nematodes. This is the first report of mermithid nematodes infecting spiders in Hawaii, and an initial attempt to characterize this host-parasite interaction. Because immature mermithids were not morphologically identifiable, a molecular identification was performed. A phylogenetic analysis based on 18S small ribosomal subunit nuclear gene sequences suggested that Hawaiian spider mermithids are more closely related to a mainland presumptive Aranimemis species that infects spiders, than to an insect-infecting mermithid collected on Oahu, HI, or to Mermis nigrescens, also a parasite of insects. Measured infection prevalence was low (ranging from 0 to 4%) but differed significantly among forest fragments. Infection prevalence was associated significantly with fragment area, but not with spider density nor spider species richness. Results suggest that mermithid populations are sensitive to habitat fragmentation, but that changes in infection prevalence do not appear to affect spider community structure.