Abstract Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush (Catharus guttatus), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.
more »
« less
Landscape Genomics Provides Evidence of Ecotypic Adaptation and a Barrier to Gene Flow at Treeline for the Arctic Foundation Species Eriophorum vaginatum
Global climate change has resulted in geographic range shifts of flora and fauna at a global scale. Extreme environments, like the Arctic, are seeing some of the most pronounced changes. This region covers 14% of the Earth’s land area, and while many arctic species are widespread, understanding ecotypic variation at the genomic level will be important for elucidating how range shifts will affect ecological processes. Tussock cottongrass ( Eriophorum vaginatum L.) is a foundation species of the moist acidic tundra, whose potential decline due to competition from shrubs may affect ecosystem stability in the Arctic. We used double-digest Restriction Site-Associated DNA sequencing to identify genomic variation in 273 individuals of E. vaginatum from 17 sites along a latitudinal gradient in north central Alaska. These sites have been part of 30 + years of ecological research and are inclusive of a region that was part of the Beringian refugium. The data analyses included genomic population structure, demographic models, and genotype by environment association. Genome-wide SNP investigation revealed environmentally associated variation and population structure across the sampled range of E. vaginatum , including a genetic break between populations north and south of treeline. This structure is likely the result of subrefugial isolation, contemporary isolation by resistance, and adaptation. Forty-five candidate loci were identified with genotype-environment association (GEA) analyses, with most identified genes related to abiotic stress. Our results support a hypothesis of limited gene flow based on spatial and environmental factors for E. vaginatum , which in combination with life history traits could limit range expansion of southern ecotypes northward as the tundra warms. This has implications for lower competitive attributes of northern plants of this foundation species likely resulting in changes in ecosystem productivity.
more »
« less
- Award ID(s):
- 2220863
- PAR ID:
- 10400870
- Date Published:
- Journal Name:
- Frontiers in Plant Science
- Volume:
- 13
- ISSN:
- 1664-462X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The phenology of Arctic plants is an important determinant of the pattern of carbon uptake and may be highly sensitive to continued rapid climate change. Eriophorum vaginatum L. (Cyperaceae) has a disproportionate influence over ecosystem processes in moist acidic tundra, but it is unclear whether its growth and phenology will remain competitive in the future. We investigated whether northern tundra ecotypes of E. vaginatum could extend their growing season in response to direct warming and transplanting into southern ecosystems. At the same time, we examined whether southern ecotypes could adjust their growth patterns in order to thrive further north, should they disperse quickly enough. Detailed phenology measurements across three reciprocal transplant gardens over a 2-year period showed that some northern ecotypes were capable of growing for longer when conditions were favourable, but their biomass and growing season length was still shorter than those of the southern ecotype. Southern ecotypes retained large leaf length when transplanted north and mirrored the growing season length better than the others, mainly owing to immediate green-up after snowmelt. All ecotypes retained the same senescence timing, regardless of environment, indicating a strong genetic control. Eriophorum vaginatum may remain competitive in a warming world if southern ecotypes can migrate north.more » « less
-
Abstract Global change is impacting biodiversity across all habitats on earth. New selection pressures from changing climatic conditions and other anthropogenic activities are creating heterogeneous ecological and evolutionary responses across many species' geographic ranges. Yet we currently lack standardised and reproducible tools to effectively predict the resulting patterns in species vulnerability to declines or range changes.We developed an informatic toolbox that integrates ecological, environmental and genomic data and analyses (environmental dissimilarity, species distribution models, landscape connectivity, neutral and adaptive genetic diversity, genotype‐environment associations and genomic offset) to estimate population vulnerability. In our toolbox, functions and data structures are coded in a standardised way so that it is applicable to any species or geographic region where appropriate data are available, for example individual or population sampling and genomic datasets (e.g. RAD‐seq, ddRAD‐seq, whole genome sequencing data) representing environmental variation across the species geographic range.To demonstrate multi‐species applicability, we apply our toolbox to three georeferenced genomic datasets for co‐occurring East African spiny reed frogs (Afrixalus fornasini, A. delicatusandA. sylvaticus) to predict their population vulnerability, as well as demonstrating that range loss projections based on adaptive variation can be accurately reproduced from a previous study using data for two European bat species (Myotis escaleraiandM. crypticus).Our framework sets the stage for large scale, multi‐species genomic datasets to be leveraged in a novel climate change vulnerability framework to quantify intraspecific differences in genetic diversity, local adaptation, range shifts and population vulnerability based on exposure, sensitivity and landscape barriers.more » « less
-
Identifying the evolutionary and ecological mechanisms that drive lineage diversification in the species-rich tropics is of broad interest to evolutionary biologists. Here, we use phylogeographic and demographic analyses of genomic scale RADseq data to assess the impact of a large geographic feature, the Amazon River, on lineage formation in a venomous pitviper, Bothrops atrox. We compared genetic differentiation in samples from four sites near Santarem, Brazil that spanned the Amazon and represented major habitat types. A species delimitation analysis identified each population as a distinct evolutionary lineage while a species tree analysis with populations as taxa revealed a phylogenetic tree consistent with dispersal across the Amazon from north to south. Phylogenetic analyses of mtDNA variation confirmed this pattern and suggest that all lineages originated during the mid- to late-Pleistocene. Historical demographic analyses support a population model of lineage formation through isolation between lineages with low ongoing migration between large populations and reject a model of differentiation through isolation by distance alone. Our results provide a rare example of a phylogeographic pattern demonstrating dispersal over evolutionary time scales across a large tropical river and suggest a role for the Amazon River as a driver of in-situ divergence by both impeding (but not preventing) gene flow and through parapatric differentiation along an ecological gradient.more » « less
-
Range shifts in a foundation sedge potentially induce large Arctic ecosystem carbon losses and gainsAbstract Foundation species have disproportionately large impacts on ecosystem structure and function. As a result, future changes to their distribution may be important determinants of ecosystem carbon (C) cycling in a warmer world. We assessed the role of a foundation tussock sedge (Eriophorum vaginatum) as a climatically vulnerable C stock using field data, a machine learning ecological niche model, and an ensemble of terrestrial biosphere models (TBMs). Field data indicated that tussock density has decreased by ∼0.97 tussocks per m2over the past ∼38 years on Alaska’s North Slope from ∼1981 to 2019. This declining trend is concerning because tussocks are a large Arctic C stock, which enhances soil organic layer C stocks by 6.9% on average and represents 745 Tg C across our study area. By 2100, we project that changes in tussock density may decrease the tussock C stock by 41% in regions where tussocks are currently abundant (e.g. −0.8 tussocks per m2and −85 Tg C on the North Slope) and may increase the tussock C stock by 46% in regions where tussocks are currently scarce (e.g. +0.9 tussocks per m2and +81 Tg C on Victoria Island). These climate-induced changes to the tussock C stock were comparable to, but sometimes opposite in sign, to vegetation C stock changes predicted by an ensemble of TBMs. Our results illustrate the important role of tussocks as a foundation species in determining future Arctic C stocks and highlight the need for better representation of this species in TBMs.more » « less