Genomics-informed delineation of conservation units in a desert amphibian

Delineating conservation units (CUs, e.g., evolutionarily significant units, ESUs, and management units, MUs) is critical to the recovery of declining species because CUs inform both listing status and management actions. Genomic data have strengths and limitations in informing CU delineation and related management questions in natural systems. We illustrate the value of using genomic data in combination with landscape, dispersal and occupancy data to inform CU delineation in Nevada populations of the Great Basin Distinct Population Segment of the Columbia spotted frog (Rana luteiventris). R. luteiventris occupies naturally fragmented aquatic habitats in this xeric region, but beaver removal, climate change and other factors have put many of these populations at high risk of extirpation without management intervention. We addressed three objectives: (i) assessing support for ESUs within Nevada; (ii) evaluating and revising, if warranted, the current delineation of MUs; and (iii) evaluating genetic diversity, effective population size, adaptive differentiation and functional connectivity to inform ongoing management actions. We found little support for ESUs within Nevada but did identify potential revisions to MUs based on unique landscape drivers of connectivity that distinguish these desert populations from those in the northern portion of the species range. Effective sizes were uniformly small, with low genetic diversity and weak signatures of adaptive differentiation. Our findings suggest that management actions, including translocations and genetic rescue, might be warranted. Our study illustrates how a carefully planned genetic study, designed to address priority management goals that include CU delineation, can provide multiple insights to inform conservation action.

Evolving spatial conservation prioritization with intraspecific genetic data

Spatial conservation prioritization (SCP) is a planning framework used to identify new conservation areas on the basis of the spatial distribution of species, ecosystems, and their services to human societies. The ongoing accumulation of intraspecific genetic data on a variety of species offers a way to gain knowledge of intraspecific genetic diversity and to estimate several population characteristics useful in conservation, such as dispersal and population size. Here, we review how intraspecific genetic data have been integrated into SCP and highlight their potential for identifying conservation area networks that represent intraspecific genetic diversity comprehensively and that ensure the long-term persistence of biodiversity in the face of global change.

Genetic data improves niche model discrimination and alters the direction and magnitude of climate change forecasts

Ecological niche models (ENMs) have classically operated under the simplifying assumptions that there are no barriers to gene flow, species are genetically homogeneous (i.e., no population-specific local adaptation), and all individuals share the same niche. Yet, these assumptions are violated for most broadly distributed species. Here, we incorporate genetic data from the widespread riparian tree species narrowleaf cottonwood (Populus angustifolia) to examine whether including intraspecific genetic variation can alter model performance and predictions of climate change impacts. We found that (1) P. angustifolia is differentiated into six genetic groups across its range from México to Canada and (2) different populations occupy distinct climate niches representing unique ecotypes. Comparing model discriminatory power, (3) all genetically informed ecological niche models (gENMs) outperformed the standard species-level ENM (3-14% increase in AUC; 1-23% increase in pROC). Furthermore, (4) gENMs predicted large differences among ecotypes in both the direction and magnitude of responses to climate change and (5) revealed evidence of niche divergence, particularly for the Eastern Rocky Mountain ecotype. (6) Models also predicted progressively increasing fragmentation and decreasing overlap between ecotypes. Contact zones are often hotspots of diversity that are critical for supporting species' capacity to respond to present and future climate change, thus predicted reductions in connectivity among ecotypes is of conservation concern. We further examined the generality of our findings by comparing our model developed for a higher elevation Rocky Mountain species with a related desert riparian cottonwood, P. fremontii. Together our results suggest that incorporating intraspecific genetic information can improve model performance by addressing this important source of variance. gENMs bring an evolutionary perspective to niche modeling and provide a truly "adaptive management" approach to support conservation genetic management of species facing global change.

Conservation Genomics in the Sagebrush Sea: Population Divergence, Demographic History, and Local Adaptation in Sage-Grouse (Centrocercus spp.)

Sage-grouse are two closely related iconic species of the North American West, with historically broad distributions across sagebrush-steppe habitat. Both species are dietary specialists on sagebrush during winter, with presumed adaptations to tolerate the high concentrations of toxic secondary metabolites that function as plant chemical defenses. Marked range contraction and declining population sizes since European settlement have motivated efforts to identify distinct population genetic variation, particularly that which might be associated with local genetic adaptation and dietary specialization of sage-grouse. We assembled a reference genome and performed whole-genome sequencing across sage-grouse from six populations, encompassing both species and including several populations on the periphery of the species ranges. Population genomic analyses reaffirmed genome-wide differentiation between greater and Gunnison sage-grouse, revealed pronounced intraspecific population structure, and highlighted important differentiation of a small isolated population of greater sage-grouse in the northwest of the range. Patterns of genome-wide differentiation were largely consistent with a hypothesized role of genetic drift due to limited gene flow among populations. Inferred ancient population demography suggested persistent declines in effective population sizes that have likely contributed to differentiation within and among species. Several genomic regions with single-nucleotide polymorphisms exhibiting extreme population differentiation were associated with candidate genes linked to metabolism of xenobiotic compounds. In vitro activity of enzymes isolated from sage-grouse livers supported a role for these genes in detoxification of sagebrush, suggesting that the observed interpopulation variation may underlie important local dietary adaptations, warranting close consideration for conservation strategies that link sage-grouse to the chemistry of local sagebrush.

Fish conservation in the land of steppe and sky: Evolutionarily significant units of threatened salmonid species in Mongolia mirror major river basins

Mongolia's salmonids are suffering extensive population declines; thus, more comprehensive fisheries management and conservation strategies are required. To assist with their development, a better understanding of the genetic structure and diversity of these threatened species would allow a more targeted approach for preserving genetic variation and ultimately improve long-term species recoveries. It is hypothesized that the unfragmented river basins that have persisted across Mongolia provide unobstructed connectivity for resident salmonid species. Thus, genetic structure is expected to be primarily segregated between major river basins. We tested this hypothesis by investigating the population structure for three salmonid genera (Hucho, Brachymystax and Thymallus) using different genetic markers to identify evolutionarily significant units (ESUs) and priority rivers to focus conservation efforts. Fish were assigned to separate ESUs when the combined evidence of mitochondrial and nuclear data indicated genetic isolation. Hucho taimen exhibited a dichotomous population structure forming two ESUs, with five priority rivers. Within the Brachymystax genus, there were three B. lenokESUs and one B. tumensisESU, along with six priority rivers. While B. tumensiswas confirmed to display divergent mtDNA haplotypes, haplotype sharing between these two congeneric species was also identified. For T. baicalensis,only a single ESU was assigned, with five priority rivers identified plus Lake Hovsgol. Additionally, we confirmed that T. nigrescens from Lake Hovsgol is a synonym of T. baicalensis. Across all species, the most prominent pattern was strong differentiation among major river basins with low differentiation and weak patterns of isolation by distance within river basins, which corroborated our hypothesis of high within-basin connectivity across Mongolia. This new genetic information provides authorities the opportunity to distribute resources for management between ESUs while assigning additional protection for the more genetically valuable salmonid rivers so that the greatest adaptive potential within each species can be preserved.

Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors

Climate change and increasing habitat loss greatly impact species survival, requiring range shifts, phenotypic plasticity and/or evolutionary change for long-term persistence, which may not readily occur unaided in threatened species. Therefore, defining conservation actions requires a detailed assessment of evolutionary factors. Existing genetic diversity needs to be thoroughly evaluated and spatially mapped to define conservation units (CUs) in an evolutionary context, and we address that here. We also propose a multidisciplinary approach to determine corridors and functional connectivity between CUs by including genetic diversity in the modelling while controlling for isolation by distance and phylogeographic history. We evaluate our approach on a Near Threatened Iberian endemic rodent by analysing genotyping-by-sequencing (GBS) genomic data from 107 Cabrera voles (Microtus cabrerae), screening the entire species distribution to define categories of CUs and their connectivity: We defined six management units (MUs) which can be grouped into four evolutionarily significant units (ESUs) and three (putatively) adaptive units (AUs). We demonstrate that the three different categories of CU can be objectively defined using genomic data, and their characteristics and connectivity can inform conservation decision-making. In particular, we show that connectivity of the Cabrera vole is very limited in eastern Iberia and that the pre-Pyrenean and part of the Betic geographic nuclei contribute the most to the species genetic diversity. We argue that a multidisciplinary framework for CU definition is essential and that this framework needs a strong evolutionary basis.

Mitogenomic sequences support a north-south subspecies subdivision within Solenodon paradoxus

Solenodons are insectivores found only in Hispaniola and Cuba, with a Mesozoic divergence date versus extant mainland mammals. Solenodons are the oldest lineage of living eutherian mammal for which a mitogenome sequence has not been reported. We determined complete mitogenome sequences for six Hispaniolan solenodons (Solenodon paradoxus) using next-generation sequencing. The solenodon mitogenomes were 16,454-16,457 bp long and carried the expected repertoire of genes. A mitogenomic phylogeny confirmed the basal position of solenodons relative to shrews and moles, with solenodon mitogenomes estimated to have diverged from those of other mammals ca. 78 Mya. Control region sequences of solenodons from the northern (n = 3) and southern (n = 5) Dominican Republic grouped separately in a network, with F_ST = 0.72 (p = 0.036) between north and south. This regional genetic divergence supports previous morphological and genetic reports recognizing northern (S. p. paradoxus) and southern (S. p. woodi) subspecies in need of separate conservation plans.

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range-wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid-latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black-tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.

Phylogeography of the pallid kangaroo mouse, Microdipodops pallidus: a sand-obligate endemic of the Great Basin, western North America

AIM: Kangaroo mice, genus Microdipodops Merriam, are endemic to the Great Basin and include two species: M. pallidus Merriam and M. megacephalus Merriam. The pallid kangaroo mouse, M. pallidus, is a sand-obligate desert rodent. Our principal intent is to identify its current geographical distribution and to formulate a phylogeographical hypothesis for this taxon. In addition, we test for orientation patterns in haplotype sharing for evidence of past episodes of movement and gene flow. LOCATION: The Great Basin Desert region of western North America, especially the sandy habitats of the Lahontan Trough and those in south-central Nevada. METHODS: Mitochondrial DNA sequence data from portions of three genes (16S ribosomal RNA, cytochrome b, and transfer RNA for glutamic acid) were obtained from 98 individuals of M. pallidus representing 27 general localities sampled throughout its geographical range. Molecular sequence data were analysed using neighbour-joining, maximum-parsimony, maximum-likelihood and Bayesian methods of phylogenetic inference. Directional analysis of phylogeographical patterns, a novel method, was used to examine angular measurements of haplotype sharing between pairs of localities to detect and quantify historical events pertaining to movement patterns and gene flow. RESULTS: Collecting activities showed that M. pallidus is a rather rare rodent (mean trapping success was 2.88%), and its distribution has changed little from that determined three-quarters of a century ago. Two principal phylogroups, distributed as eastern and western moieties, are evident from the phylogenetic analyses (mean sequence divergence for cytochrome b is c. 8%). The western clade shows little phylogenetic structure and seems to represent a large polytomy. In the eastern clade, however, three subgroups are recognized. Nine of the 42 unique composite haplotypes are present at two or more localities and are used for the orientation analyses. Axial data from haplotype sharing between pairwise localities show significant, non-random angular patterns: a north-west to south-east orientation in the western clade, and a north-east to south-west directional pattern in the eastern clade. MAIN CONCLUSIONS: The geographical range of M. pallidus seems to be remarkably stable in historical times and does not show a northward (or elevationally upward) movement trend, as has been reported for some other kinds of organism in response to global climate change. The eastern and western clades are likely to represent morphologically cryptic species. Estimated times of divergence of the principal clades of M. pallidus (4.38 Ma) and between M. pallidus and M. megacephalus (8.1 Ma; data from a related study) indicate that kangaroo mice diverged much earlier than thought previously. The phylogeographical patterns described here may serve as a model for other sand-obligate members of the Great Basin Desert biota.