Harnessing genomics for delineating conservation units

Genomic data have the potential to revolutionize the delineation of conservation units (CUs) by allowing the detection of adaptive genetic variation, which is otherwise difficult for rare, endangered species. In contrast to previous recommendations, we propose that the use of neutral versus adaptive markers should not be viewed as alternatives. Rather, neutral and adaptive markers provide different types of information that should be combined to make optimal management decisions. Genetic patterns at neutral markers reflect the interaction of gene flow and genetic drift that affects genome-wide variation within and among populations. This population genetic structure is what natural selection operates on to cause adaptive divergence. Here, we provide a new framework to integrate data on neutral and adaptive markers to protect biodiversity.     

Phosphoglucose isomerase (<Emphasis Type="Italic">Pgi</Emphasis>) performance and fitness effects among Arthropods and its potential role as an adaptive marker in conservation genetics

The development of conservation genomics will be greatly aided by the use of neutral as well as adaptive molecular markers. Identifying novel adaptive molecular markers that have general application across diverse taxa is challenging, especially in Arthropods where few if any examples of balanced polymorphisms exist that are shared across species. A review of literature on the Pgi gene provides strong evidence for population level fitness consequences of genetic variation in this gene, across very diverse lineages of Arthropods. While these observations demonstrate the potential of using Pgi as an adaptive molecular marker, this gene is fundamentally different from the adaptive markers MHC and SI. Rather than providing insights into individual genetic health, Pgi appears to have a role in conservation genetics by providing insights into gene by environment interactions, local adaptation and evolutionary significant units, and potentially even morphologically cryptic dispersal phenotypes. These findings argue for studying Pgi variation in more species, as it appears central to the goals of conservation genomics.     

Genome-wide SNPs reveal low effective population size within confined management units of the highly vagile Galapagos shark (<Emphasis Type="Italic">Carcharhinus galapagensis</Emphasis>)

The Galapagos shark (Carcharhinus galapagensis) is one of over thirty shark species inhabiting the Galápagos Marine Reserve (GMR), where it is a priority species for conservation. Identifying stock structure and effective population size for species-specific management and effective conservation of this top predator is important. We examined stock structure, connectivity and effective population size of Galapagos sharks among GMR locations using genome-wide neutral Single Nucleotide Polymorphism (8103 SNP) and mtDNA markers. Potential historical gene flow and/or sex-biased dispersal were also examined using the mitochondrial control region (997 bp). Cluster analyses of neutral SNPs revealed two differentiated stocks in the GMR—a western (Isabela Island) and eastern (San Cristóbal and Santa Cruz Islands) stock. Effective population size (N_e) estimates of approximately 200 suggest these populations are susceptible to ongoing natural and anthropogenic stressors and are of concern for long term resilience of populations. Mitochondrial DNA failed to identify distinct stocks, with AMOVA analyses indicating most genetic variation occurs within, rather than among locations. This pattern of genome-wide nuclear (but not mtDNA) discrimination among neighbouring islands either points to possible sex-biased dispersal by females or identifies limitations of the single organelle mtDNA marker at such small spatial scales. Regional differences across the archipelago or in behaviour may be implicated in the observed population structure. Further research focusing on a larger, Pacific wide analysis of population connectivity and effective population size at a broader spatial scale is required, to estimate the extent of discreteness and potential local adaptation. Potential adaptive units (AUs) in Galapagos sharks should ultimately be identified to leverage adaptive management and fisheries forensics applications.     

The Conservation Value of Peripheral Populations and a Relationship Between Quantitative Trait and Molecular Variation

The adaptive potential of populations and therefore their ability to cope with rapid environmental changes is a question of paramount fundamental and applied importance. However, what is still not clear is the effect of population position within the species range (i.e. core vs. edge) on population adaptive potential, and whether the adaptive potential can be predicted from extent of neutral molecular variation. In this study, we compared the extent and structure of neutral (SSR) and presumably adaptive quantitative trait genetic variation in populations of Triticum dicoccoides sampled at the species range core and two opposite edges, and related this information to multigenerational performance of plants experimentally introduced beyond the range edge. The plants from the species arid edge performed worse than plants from the more mesic core in extreme desert conditions. The core and edge populations did not differ in extent of SSR variation. In contrast to the neutral genetic variation, there was lower quantitative trait variation in the two edge as compared with the core population for many traits, and no trait in any edge population had higher variation than the core population or either of its habitats. Reduced variation in selectively important traits indicates a lower adaptive potential of the two edge as compared with the core population. Our results imply (1) that extent of variation in quantitative traits can predict plant performance in novel environments while extent of variation in molecular markers can not; and (2) caution in usage of peripheral populations in such conservation actions as relocation and creation of new populations. We also warn against usage of neutral molecular variation as a surrogate for selectively important quantitative variation in conservation decisions.     

An integrated framework to identify wildlife populations under threat from climate change

Climate change is a major threat to global biodiversity that will produce a range of new selection pressures. Understanding species responses to climate change requires an interdisciplinary perspective, combining ecological, molecular and environmental approaches. We propose an applied integrated framework to identify populations under threat from climate change based on their extent of exposure, inherent sensitivity due to adaptive and neutral genetic variation and range shift potential. We consider intraspecific vulnerability and population‐level responses, an important but often neglected conservation research priority. We demonstrate how this framework can be applied to vertebrates with limited dispersal abilities using empirical data for the bat Plecotus austriacus. We use ecological niche modelling and environmental dissimilarity analysis to locate areas at high risk of exposure to future changes. Combining outlier tests with genotype–environment association analysis, we identify potential climate‐adaptive SNPs in our genomic data set and differences in the frequency of adaptive and neutral variation between populations. We assess landscape connectivity and show that changing environmental suitability may limit the future movement of individuals, thus affecting both the ability of populations to shift their distribution to climatically suitable areas and the probability of evolutionary rescue through the spread of adaptive genetic variation among populations. Therefore, a better understanding of movement ecology and landscape connectivity is needed for predicting population persistence under climate change. Our study highlights the importance of incorporating genomic data to determine sensitivity, adaptive potential and range shift potential, instead of relying solely on exposure to guide species vulnerability assessments and conservation planning.     

Genetic structure in peripheral populations of the natterjack toad, <Emphasis Type="Italic">Bufo calamita</Emphasis>, as revealed by AFLP

Decreased fitness due to loss of genetic variation is a well recognised issue in conservation biology. Along the Swedish west coast, the endangered natterjack toad (Bufo calamita) occurs on, for the species, highly unusual habitat of rocky islands. Although the toads inhabit a restricted geographical area (maximum distance between the populations is 71 km), the fragmented nature of the landscape makes the genetic properties of the populations of conservation interest. However, lack of genetic variation found using conventional methods (microsatellites) has impeded genetic studies within these peripheral populations so far. In this study we assess population structure and genetic variation among seven of these fringe populations using 105 polymorphic Amplified Fragment Length Polymorphism (AFLP) loci. We found a well-defined population structure without evidence for isolation by distance, implying restricted gene flow between populations. Additionally, the populations differed in their amount of genetic variation, emphasizing the need to monitor genetically impoverished populations for possible declines mediated by inbreeding depression and reduced adaptive potential. Conservation implications for these unique populations are discussed in the light of our results.     

Conserving marine biodiversity: insights from life-history trait candidate genes in Atlantic cod (Gadus morhua)

Recent technological developments have facilitated an increased focus on identifying genomic regions underlying adaptive trait variation in natural populations, and it has been advocated that this information should be important for designating population units for conservation. In marine fishes, phenotypic studies have suggested adaptation through divergence of life-history traits among natural populations, but the distribution of adaptive genetic variation in these species is still relatively poorly known. In this study, we extract information about the geographical distribution of genetic variation for 33 single nucleotide polymorphisms (SNPs) associated with life-history trait candidate genes, and compare this to variation in 70 putatively neutral SNPs in Atlantic cod (Gadus morhua). We analyse samples covering the major population complexes in the eastern Atlantic and find strong evidence for non-neutral levels and patterns of population structuring for several of the candidate gene-associated markers, including two SNPs in the growth hormone 1 gene. Thus, this study aligns with findings from phenotypic studies, providing molecular data strongly suggesting that these or closely linked genes are under selection in natural populations of Atlantic cod. Furthermore, we find that patterns of variation in outlier markers do not align with those observed at selectively neutral markers, and that outlier markers identify conservation units on finer geographical scales than those revealed when analysing only neutral markers. Accordingly, results also suggest that information about adaptive genetic variation will be useful for targeted conservation and management in this and other marine species.     

Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Low neutral and immunogenetic diversity in northern fringe populations of the green toad Bufotes viridis: implications for conservation

Genetic variation is often lower at high latitudes, which may compromise the adaptability and hence survival of organisms. Here we show that genetic variability is negatively correlated with northern latitude in European green toads (Bufotes viridis). The result holds true for both putatively neutral microsatellite variation and supposedly adaptive MHC Class IIB variation. In particular, our findings have bearing on the conservation status of this species in Sweden, on the northern limit of its distribution where local populations are small and fragmented. These genetically impoverished populations are closely related to other populations found around the Baltic Sea basin. The low neutral and adaptive variation in these fringe populations compared to population at central ranges confirms a pattern shared across all other amphibians so far studied. In Sweden, the situation of green toads is of concern as the remaining populations may not have the evolutionary potential to cope with present and future environmental challenges.

     

Conservation Investment for Rare Plants in Urban Environments

Budgets for species conservation limit actions. Expending resources in areas of high human density is costly and generally considered less likely to succeed. Yet, coastal California contains both a large fraction of narrowly endemic at-risk plant species as well as the state''s three largest metropolitan regions. Hence understanding the capacity to protect species along the highly urbanized coast is a conservation priority. We examine at-risk plant populations along California''s coastline from San Diego to north of San Francisco to better understand whether there is a relationship between human population density and: i) performance of at-risk plant populations; and ii) conservation spending. Answering these questions can help focus appropriate strategic conservation investment. Rare plant performance was measured using the annualized growth rate estimate between census periods using the California Natural Diversity Database. Human density was estimated using Census Bureau statistics from the year 2000. We found strong evidence for a lack of a relationship between human population density and plant population performance in California''s coastal counties. Analyzing US Endangered Species expenditure reports, we found large differences in expenditures among counties, with plants in San Diego County receiving much higher expenditures than other locations. We found a slight positive relationship between expenditures on behalf of endangered species and human density. Together these data support the argument that conservation efforts by protecting habitats within urban environments are not less likely to be successful than in rural areas. Expenditures on behalf of federally listed endangered and threatened plants do not appear to be related to proximity to human populations. Given the evidence of sufficient performance in urban environments, along with a high potential to leverage public support for nature in urban environments, expenditures in these areas appear to be an appropriate use of conservation funds.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

Consequences of fragmentation for the ability to adapt to novel environments in experimental <Emphasis Type="Italic">Drosophila</Emphasis> metapopulations

We used experimental populations of Drosophila melanogaster, which had either been subdivided (metapopulations) or kept undivided for 40 generations, to study the consequences of population subdivision for the tolerance and adaptive response after six generations of exposure to novel environmental factors (high temperature, medium with ethanol or salt added) for traits with different genetic architectures. In this setup, we attempted to separate the effects of the loss of fitness due to inbreeding (i.e., the survival upon first exposure to stress) from the loss of adaptive potential due to the lack of genetic variation. To place our experimental results in a more general perspective, we used individual-based simulations combining different options of levels of gene flow, intensity of selection and genetic architecture to derive quantitative hypotheses of the effects of these factors on the adaptive response to stress. We observed that population subdivision resulted in substantial inter-deme variation in tolerance due to redistribution of genetic variation from within demes to among demes. In line with the simulation results, the adaptive response was generally lower in the subdivided than in the undivided populations, particularly so for high temperature. We observed pronounced differences between stress factors that are likely related to the different genetic architectures involved in resistance to these factors. From a conservation genetics viewpoint, our results have two important implications: (i) Long-term fragmentation in combination with restricted gene flow will limit the adaptive potential of individual subpopulations because adaptive variation will become distributed among populations rather than within populations. (ii) The genetic architecture of the trait(s) under selection is of great significance to understand the possible responses to novel stresses that may be expected.     

Genomics and conservation genetics

In large part, the relevance of genetics to conservation rests on the premise that neutral marker variation in populations reflects levels of detrimental and adaptive genetic variation. Despite its prominence, this tenet has been difficult to evaluate, until now. As we discuss here, genome sequence information and new technological and bioinformatics platforms now enable comprehensive surveys of neutral variation and more direct inferences of detrimental and adaptive variation in species with sequenced genomes and in 'genome-enabled' endangered taxa. Moreover, conservation schemes could begin to consider specific pathological genetic variants. A new conservation genetic agenda would utilize data from enhanced surveys of genomic variation in endangered species to better manage functional genetic variation.     

Genetic diversity of Allium munzii (Amaryllidaceae), a rare southern California species and implication for its conservation

An understanding of genetic diversity within and among populations of rare plant species is a prerequisite to develop effective conservation management strategies and reintroduction programs. Allium munzii is a narrow endemic species distributed in western Riverside County, California, USA and known from 18 extant element occurrences. We sampled 119 individuals from 11 element occurrences and investigated within and among population genetic diversity using two variable chloroplast markers (rpL32–trnL intergenic spacer and rpoC1 intron). Of the total genetic variation detected in A. munzii, 87.65% was due to differences among occurrences. Furthermore, our results revealed that most of the element occurrences are strongly genetically differentiated. There are low levels of gene flow between occurrences, not due to isolation by distance but possibly resulting from habitat fragmentation. Non-significant values of Tajima's D and Fu's Fs were found in all occurrences suggesting no demographic expansion in A. munzii. Ex situ seed and bulb conservation is recommended to enable introduction of individuals to occurrences with low abundance and genetic diversity.     

Genetic differentiation across eastern Pacific oceanographic barriers in the threatened seahorse <Emphasis Type="Italic">Hippocampus ingens</Emphasis>

Understanding the population structure and evolutionary history of the eastern Pacific seahorse Hippocampus ingens is critical for the effective management of this threatened species. Life history characteristics of H. ingens (site fidelity and brooding of young) may limit gene flow and lead to population differentiation. A recent study analyzing conserved fragments of the mitochondrial cyt b and control region found no population structure. We re-assess this conclusion with a phylogeographic analysis of relationships among 115 individuals of H. ingens over a broader geographic range (San Diego Bay in California, Gulf of California, Central America, Ecuador, and Peru) based on a more variable 428 base pair fragment of the control region. This expanded analysis affirms low overall nucleotide diversity relative to other seahorses (θ_π = 0.004), and shows evidence of a recent bottleneck and population expansion since the middle Pleistocene. AMOVA analysis shows moderate overall population structure (Φ_ST = 0.10, P val = 0.00), and pairwise Φ_ST estimates indicate structure between the Gulf of California and all Pacific coast localities. Knowledge of population structure in H. ingens may improve conservation efforts by identifying evolutionarily important management units, and could determine source regions in the continuing trade of seahorses for traditional Chinese medicine. The level of genetic divergence observed between the Gulf of California and all other localities sampled may distinguish the Gulf as a separate management unit. Additional phylogeographic research with more quickly evolving genetic markers and targeted sampling at the mouth of the Gulf of California is warranted to inform strategies for conservation of this threatened seahorse.     

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus,Osmeridae), an anadromous Pacific smelt

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction‐site‐associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single‐nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F_ST outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N_e) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome‐wide neutral and adaptive genetic variation for the applications of conservation and management.     

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included F_ST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.     

Seascape genomics provides evidence for thermal adaptation and current‐mediated population structure in American lobster (Homarus americanus)

Investigating how environmental features shape the genetic structure of populations is crucial for understanding how they are potentially adapted to their habitats, as well as for sound management. In this study, we assessed the relative importance of spatial distribution, ocean currents and sea surface temperature (SST) on patterns of putatively neutral and adaptive genetic variation among American lobster from 19 locations using population differentiation (PD) approaches combined with environmental association (EA) analyses. First, PD approaches (using bayescan , arlequin and outflank ) found 28 outlier SNPs putatively under divergent selection and 9770 neutral SNPs in common. Redundancy analysis revealed that spatial distribution, ocean current‐mediated larval connectivity and SST explained 31.7% of the neutral genetic differentiation, with ocean currents driving the majority of this relationship (21.0%). After removing the influence of spatial distribution, no SST were significant for putatively neutral genetic variation whereas minimum annual SST still had a significant impact and explained 8.1% of the putatively adaptive genetic variation. Second, EA analyses (using Pearson correlation tests, bayescenv and lfmm ) jointly identified seven SNPs as candidates for thermal adaptation. Covariation at these SNPs was assessed with a spatial multivariate analysis that highlighted a significant temperature association, after accounting for the influence of spatial distribution. Among the 505 candidate SNPs detected by at least one of the three approaches, we discovered three polymorphisms located in genes previously shown to play a role in thermal adaptation. Our results have implications for the management of the American lobster and provide a foundation on which to predict how this species will cope with climate change.     

Near panmixia at the distribution‐wide scale but evidence of genetic differentiation in a geographically isolated population of the Terek Sandpiper Xenus cinereus

Populations from different parts of a species range may vary in their genetic structure, variation and dynamics. Geographically isolated populations or those located at the periphery of the range may differ from those located in the core of the range. Such peripheral populations may harbour genetic variation important for the adaptive potential of the species. We studied the distribution‐wide population genetic structure of the Terek Sandpiper Xenus cinereus using 13 microsatellite loci and the mitochondrial DNA (mtDNA) control region. In addition, we estimated whether genetic variation changes from the core towards the edge of the breeding range. We used the results to evaluate the management needs of the sampled populations. Distribution‐wide genetic structure was negligible; the only population that showed significant genetic differentiation was the geographically isolated Dnieper River basin population in Eastern Europe. The genetic variation of microsatellites decreased towards the edge of the distribution, supporting the abundant‐centre hypotheses in which the core area of the distribution preserves the most genetic variation; however, no such trend could be seen with mtDNA. Overall genetic variation was low and there were signs of past population contractions followed by expansion; this pattern is found in most northern waders. The current effective population size (N_e) is large, and therefore global conservation measures are not necessary. However, the marginal Dnieper River population needs to be considered its own management unit. In addition, the Finnish population warrants conservation actions due to its extremely small size and degree of isolation from the main range, which makes it vulnerable to genetic depletion.