skip to main content


Title: Generating continuous maps of genetic diversity using moving windows
Abstract

Genetic diversity plays a key role in maintaining population viability by preventing inbreeding depression and providing the building blocks for adaptation. Understanding how genetic diversity varies across space is, therefore, of key interest in conservation and population genetics.

Here, we introducewingen, anrpackage for calculating continuous maps of genetic diversity, including nucleotide diversity, allelic richness, and heterozygosity, from standard genotypic and spatial data using a spatial moving window approach. We provide functions to account for variation in sample size across space using rarefaction, to create kriging‐interpolated maps of genetic diversity, and to mask any areas that are outside the area of interest.

Tests with simulated and empirical datasets demonstrate thatwingencan successfully capture variation in genetic diversity across landscapes from both reduced‐representation and whole genome sequencing datasets. For reduced‐representation datasets,wingen's functions can be run easily on a standard laptop computer, and we provide options for parallelization to increase the efficiency of running larger whole genome datasets.

wingenprovides novel and computationally tractable tools for creating informative maps of genetic diversity with applications for conservation prioritization as well as population and landscape genetic analyses.

 
more » « less
Award ID(s):
1845682
NSF-PAR ID:
10401324
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Methods in Ecology and Evolution
Volume:
14
Issue:
5
ISSN:
2041-210X
Page Range / eLocation ID:
p. 1175-1181
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Conservation translocation projects must carefully balance multiple, potentially competing objectives (e.g. population viability, retention of genetic diversity, delivery of key ecological services) against conflicting stakeholder values and severe time and cost constraints. Advanced decision support tools would facilitate identifying practical solutions.

    We examined how to achieve compromise across competing objectives in conservation translocations via an examination of giant tortoises in the Galapagos Islands with ancestry from the extinct Floreana Island species (Chelonoidis niger). Efforts have begun to populate Floreana Island with tortoises genetically similar to its historical inhabitants while balancing three potentially competing objectives – restoring ecosystem services (sustaining a high tortoise population size), maximizing genome representation of the extinctC. nigerspecies and maintaining a genetically diverse population – under realistic cost constraints.

    We developed a novel approach to this conservation decision problem by coupling an individual‐based simulation model with generalized additive models and global optimization. We identified several incompatibilities among programme objectives, with quasi‐optimal single‐objective solutions (sets of management actions) differing substantially in programme duration, translocation age, incubation temperature (determinant of sex ratio) and the number of individuals directly translocated from the source population.

    Quasi‐optimal single‐objective solutions were able to produce outcomes (i.e. population size and measures of genetic diversity andC. nigergenome representation) to within 75% of their highest simulated outcomes (e.g. highest population size achieved across all simulations) within a cost constraint ofc. $2m USD, but these solutions resulted in severe declines (up to 74% reduction) in outcomes for non‐focal objectives. However, when all programme objectives were equally weighted to produce a multi‐objective solution, all objectives were met to within 90% of the highest achievable mean values across all cost constraints.

    Synthesis and applications. Multi‐objective conservation translocations are likely to encounter complex trade‐offs and conflicts among programme objectives. Here, we developed a novel combination of modelling approaches to identify optimal management strategies. We found that solutions that simultaneously addressed multiple, competing objectives performed better than single‐objective solutions. Our model‐based decision support tool demonstrates that timely, cost‐effective solutions can be identified in cases where management objectives appear to be incompatible.

     
    more » « less
  2. Summary

    Global change forecasts in ecosystems require knowledge of within‐species diversity, particularly of dominant species within communities. We assessed site‐level diversity and capacity for adaptation inBouteloua gracilis, the dominant species in the Central US shortgrass steppe biome.

    We quantified genetic diversity from 17 sites across regional scales, north to south from New Mexico to South Dakota, and local scales in northern Colorado. We also quantified phenotype and plasticity within and among sites and determined the extent to which phenotypic diversity inB. graciliswas correlated with climate.

    Genome sequencing indicated pronounced population structure at the regional scale, and local differences indicated that gene flow and/or dispersal may also be limited. Within a common environment, we found evidence of genetic divergence in biomass‐related phenotypes, plasticity, and phenotypic variance, indicating functional divergence and different adaptive potential. Phenotypes were differentiated according to climate, chiefly median Palmer Hydrological Drought Index and other aridity metrics.

    Our results indicate conclusive differences in genetic variation, phenotype, and plasticity in this species and suggest a mechanism explaining variation in shortgrass steppe community responses to global change. This analysis ofB. gracilisintraspecific diversity across spatial scales will improve conservation and management of the shortgrass steppe ecosystem in the future.

     
    more » « less
  3. Abstract Background

    Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada.

    Results

    Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θW = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θW = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θW = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’sDwas positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide.

    Conclusions

    Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

     
    more » « less
  4. Abstract

    Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species‐specific genotyping methods (e.g. array‐based genotyping) and must instead use sequencing‐based approaches. Although costs are decreasing, high‐coverage whole‐genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low‐coverage WGS combined with imputation using a high‐confidence reference panel is a cost‐effective alternative, but the accuracy of genotyping using low‐coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low‐coverage sequencing (0.1–10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (medianr2 ≥ 0.92). We show that low‐coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI‐for‐non‐model‐species). Our results endorse accurate and effective genotype imputation from low‐coverage sequencing, enabling the cost‐effective generation of population‐scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation.

     
    more » « less
  5. Key points

    Induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs) capture patient‐specific genotype–phenotype relationships, as well as cell‐to‐cell variability of cardiac electrical activity

    Computational modelling and simulation provide a high throughput approach to reconcile multiple datasets describing physiological variability, and also identify vulnerable parameter regimes

    We have developed a whole‐cell model of iPSC‐CMs, composed of single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM outputs

    We have utilized experimental data across multiple laboratories to model experimental variability and investigate subcellular phenotypic mechanisms in iPSC‐CMs

    This framework links molecular mechanisms to cellular‐level outputs by revealing unique subsets of model parameters linked to known iPSC‐CM phenotypes

    Abstract

    There is a profound need to develop a strategy for predicting patient‐to‐patient vulnerability in the emergence of cardiac arrhythmia. A promisingin vitromethod to address patient‐specific proclivity to cardiac disease utilizes induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs). A major strength of this approach is that iPSC‐CMs contain donor genetic information and therefore capture patient‐specific genotype–phenotype relationships. A cited detriment of iPSC‐CMs is the cell‐to‐cell variability observed in electrical activity. We postulated, however, that cell‐to‐cell variability may constitute a strength when appropriately utilized in a computational framework to build cell populations that can be employed to identify phenotypic mechanisms and pinpoint key sensitive parameters. Thus, we have exploited variation in experimental data across multiple laboratories to develop a computational framework for investigating subcellular phenotypic mechanisms. We have developed a whole‐cell model of iPSC‐CMs composed of simple model components comprising ion channel models with single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM data for all major ionic currents. By optimizing ionic current model parameters to multiple experimental datasets, we incorporate experimentally‐observed variability in the ionic currents. The resulting population of cellular models predicts robust inter‐subject variability in iPSC‐CMs. This approach links molecular mechanisms to known cellular‐level iPSC‐CM phenotypes, as shown by comparing immature and mature subpopulations of models to analyse the contributing factors underlying each phenotype. In the future, the presented models can be readily expanded to include genetic mutations and pharmacological interventions for studying the mechanisms of rare events, such as arrhythmia triggers.

     
    more » « less