skip to main content


Title: Assessing thermal adaptation using family‐based association and F ST outlier tests in a threatened trout species
Abstract

Discovering genetic markers associated with phenotypic or ecological characteristics can improve our understanding of adaptation and guide conservation of key evolutionary traits. The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) of the northern Great Basin Desert, USA, demonstrated exceptional tolerance to high temperatures in the desert lakes where it resided historically. This trait is central to a conservation hatchery effort to protect the genetic legacy of the nearly extinct lake ecotype. We genotyped full‐sibling families from this conservation broodstock and samples from the only two remaining, thermally distinct, native lake populations at 4,644 new single nucleotide polymorphisms (SNPs). Family‐based genome‐wide association testing of the broodstock identified nine and 26 SNPs associated with thermal tolerance (p < 0.05 andp < 0.1), measured in a previous thermal challenge experiment. Genes near the associated SNPs had complex functions related to immunity, growth, metabolism and ion homeostasis. Principal component analysis using the thermotolerance‐related SNPs showed unexpected divergence between the conservation broodstock and the native lake populations at these loci.FSToutlier tests on the native lake populations identified 18 loci shared between two or more of the tests, with two SNPs identified by all three tests (p < 0.01); none overlapped with loci identified by association testing in the broodstock. A recent history of isolation and the complex genetic and demographic backgrounds of Lahontan cutthroat trout probably limited our ability to find shared thermal tolerance loci. Our study extends the still relatively rare application of genomic tools testing for markers associated with important phenotypic or environmental characteristics in species of conservation concern.

 
more » « less
NSF-PAR ID:
10375498
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
28
Issue:
10
ISSN:
0962-1083
Page Range / eLocation ID:
p. 2573-2593
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Understanding the consequences of exotic diseases on native forests is important to evolutionary ecology and conservation biology because exotic pathogens have drastically altered US eastern deciduous forests.Cornus floridaL. (flowering dogwood tree) is one such species facing heavy mortality. Characterizing the genetic structure ofC. floridapopulations and identifying the genetic signature of adaptation to dogwood anthracnose (an exotic pathogen responsible for high mortality) remain vital for conservation efforts. By integrating genetic data from genotype by sequencing (GBS) of 289 trees across the host species range and distribution of disease, we evaluated the spatial patterns of genetic variation and population genetic structure ofC. floridaand compared the pattern to the distribution of dogwood anthracnose. Using genome‐wide association study and gradient forest analysis, we identified genetic loci under selection and associated with ecological and diseased regions. The results revealed signals of weak genetic differentiation of three or more subgroups nested within two clusters—explaining up to 2%–6% of genetic variation. The groups largely corresponded to the regions within and outside the eastern Hot‐Continental ecoregion, which also overlapped with areas within and outside the main distribution of dogwood anthracnose. The fungal sequences contained in the GBS data of sampled trees bolstered visual records of disease at sampled locations and were congruent with the reported range ofDiscula destructiva, suggesting that fungal sequences within‐host genomic data were informative for detecting or predicting disease. The genetic diversity between populations at diseased vs. disease‐free sites across the range ofC. floridashowed no significant difference. We identified 72 single‐nucleotide polymorphisms (SNPs) from 68 loci putatively under selection, some of which exhibited abrupt turnover in allele frequencies along the borders of the Hot‐Continental ecoregion and the range of dogwood anthracnose. One such candidate SNP was independently identified in two prior studies as a possible L‐type lectin‐domain containing receptor kinase. Although diseased and disease‐free areas do not significantly differ in genetic diversity, overall there are slight trends to indicate marginally smaller amounts of genetic diversity in disease‐affected areas. Our results were congruent with previous studies that were based on a limited number of genetic markers in revealing high genetic variation and weak population structure inC. florida.

     
    more » « less
  2. Abstract

    Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low‐coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L,ROCK1,SIX6). Genotype–environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high‐elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high‐elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

     
    more » « less
  3. Abstract

    The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish,Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine‐tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

     
    more » « less
  4. Abstract

    Understanding how genetic diversity is distributed across spatiotemporal scales in species of conservation or management concern is critical for identifying large‐scale mechanisms affecting local conservation status and implementing large‐scale biodiversity monitoring programmes. However, cross‐scale surveys of genetic diversity are often impractical within single studies, and combining datasets to increase spatiotemporal coverage is frequently impeded by using different sets of molecular markers. Recently developed molecular tools make surveys based on standardized single‐nucleotide polymorphism (SNP) panels more feasible than ever, but require existing genomic information. Here, we conduct the first survey of genome‐wide SNPs across the native range of brook trout (Salvelinus fontinalis), a cold‐adapted species that has been the focus of considerable conservation and management effort across eastern North America. Our dataset can be leveraged to easily design SNP panels that allow datasets to be combined for large‐scale analyses. We performed restriction site‐associated DNA sequencing for wild brook trout from 82 locations spanning much of the native range and domestic brook trout from 24 hatchery strains used in stocking efforts. We identified over 24,000 SNPs distributed throughout the brook trout genome. We explored the ability of these SNPs to resolve relationships across spatial scales, including population structure and hatchery admixture. Our dataset captures a wide spectrum of genetic diversity in native brook trout, offering a valuable resource for developing SNP panels. We highlight potential applications of this resource with the goal of increasing the integration of genomic information into decision‐making for brook trout and other species of conservation or management concern.

     
    more » « less
  5. Abstract Context

    Processes that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species.Fontaineais a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species ofFontaineaare subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.

    Objectives

    Here, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatenedFontaineaspecies. Our aim was to compare spatial patterns of genetic variation between the vulnerableFontainea australis(Southern Fontainea) and critically endangeredF. oraria(Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.

    Methods

    We used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.

    Results

    Our study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.

    Conclusions

    Divergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the twoFontaineaspecies. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

     
    more » « less