skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Ecological genomics of local adaptation in Cornus florida L. by genotyping by sequencing

Discovering local adaptation, its genetic underpinnings, and environmental drivers is important for conserving forest species. Ecological genomic approaches coupled with next‐generation sequencing are useful means to detect local adaptation and uncover its underlying genetic basis in nonmodel species. We report results from a study on flowering dogwood trees (Cornus florida L.) using genotyping by sequencing (GBS). This species is ecologically important to easternUSforests but is severely threatened by fungal diseases. We analyzed subpopulations in divergent ecological habitats within North Carolina to uncover loci under local selection and associated with environmental–functional traits or disease infection. At this scale, we tested the effect of incorporating additional sequencing before scaling for a broader examination of the entire range. To test for biases ofGBS, we sequenced two similarly sampled libraries independently from six populations of three ecological habitats. We obtained environmental–functional traits for each subpopulation to identify associations with genotypes via latent factor mixed modeling (LFMM) and gradient forests analysis. To test whether heterogeneity of abiotic pressures resulted in genetic differentiation indicative of local adaptation, we evaluatedFstper locus while accounting for genetic differentiation between coastal subpopulations and Piedmont‐Mountain subpopulations. Of the 54 candidate loci with sufficient evidence of being under selection among both libraries, 28–39 were Arlequin–BayeScanFstoutliers. ForLFMM, 45 candidates were associated with climate (of 54), 30 were associated with soil properties, and four were associated with plant health. Reanalysis of combined libraries showed that 42 candidate loci still showed evidence of being under selection. We conclude environment‐driven selection on specific loci has resulted in local adaptation in response to potassium deficiencies, temperature, precipitation, and (to a marginal extent) disease. High allele turnover along ecological gradients further supports the adaptive significance of loci speculated to be under selection.

more » « less
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecology and Evolution
Page Range / eLocation ID:
p. 441-465
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Despite its economic importance as a bioenergy crop and key role in riparian ecosystems, little is known about genetic diversity and adaptation of the eastern cottonwood (Populus deltoides). Here, we report the first population genomics study for this species, conducted on a sample of 425 unrelated individuals collected in 13 states of the southeastern United States. The trees were genotyped by targeted resequencing of 18,153 genes and 23,835 intergenic regions, followed by the identification of single nucleotide polymorphisms (SNPs). This naturalP. deltoidespopulation showed low levels of subpopulation differentiation (FST = 0.022–0.106), high genetic diversity (θW = 0.00100, π = 0.00170), a large effective population size (Ne ≈ 32,900), and low to moderate levels of linkage disequilibrium. Additionally, genomewide scans for selection (Tajima'sD), subpopulation differentiation (XTX), and environmental association analyses with eleven climate variables carried out with two different methods (LFMMandBAYENV2) identified genes putatively involved in local adaptation. Interestingly, many of these genes were also identified as adaptation candidates in another poplar species,Populus trichocarpa, indicating possible convergent evolution. This study constitutes the first assessment of genetic diversity and local adaptation inP. deltoidesthroughout the southern part of its range, information we expect to be of use to guide management and breeding strategies for this species in future, especially in the face of climate change.

    more » « less
  2. Abstract

    Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation‐by‐distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses inO. lurida, 13,424 single nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235SNPs across 129GBSloci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine‐scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open‐sourced research on the population structure ofO. lurida.

    more » « less
  3. Abstract

    Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwiseFST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin,USA. Using 151 putatively neutral and 29 candidate adaptiveSNPloci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables andFSTacross all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables andFSTwithinthree metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (= 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.

    more » « less
  4. Abstract

    Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (FST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured byFST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.

    more » « less
  5. Abstract

    Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identifySNPmarkers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein‐coding and nearby 5′ and 3′ untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis) exon capture data and directly from the domestic sheep genome (Ovis ariesv. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR‐basedSNPchip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositanandbayescan), we detected 28SNPloci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease‐regulating functions (e.g. Ovar‐DRA,APC,BATF2,MAGEB18), cell regulation signalling pathways (e.g.KRIT1,PI3K,ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene‐targetedSNPdiscovery and subsequentSNPchip genotyping using low‐quality samples in a nonmodel species.

    more » « less