Genome-wide shifts in climate-related variation underpin responses to selective breeding in a widespread conifer.

Locally adapted temperate tree populations exhibit genetic trade-offs among climate-related traits that can be exacerbated by selective breeding and are challenging to manage under climate change. To inform climatically adaptive forest management, we investigated the genetic architecture and impacts of selective breeding on four climate-related traits in 105 natural and 20 selectively bred lodgepole pine populations from western Canada. Growth, cold injury, growth initiation, and growth cessation phenotypes were tested for associations with 18,600 single-nucleotide polymorphisms (SNPs) in natural populations to identify "positive effect alleles" (PEAs). The effects of artificial selection for faster growth on the frequency of PEAs associated with each trait were quantified in breeding populations from different climates. Substantial shifts in PEA proportions and frequencies were observed across many loci after two generations of selective breeding for height, and responses of phenology-associated PEAs differed strongly among climatic regions. Extensive genetic overlap was evident among traits. Alleles most strongly associated with greater height were often associated with greater cold injury and delayed phenology, although it is unclear whether potential trade-offs arose directly from pleiotropy or indirectly via genetic linkage. Modest variation in multilocus PEA frequencies among populations was associated with large phenotypic differences and strong climatic gradients, providing support for assisted gene flow polices. Relationships among genotypes, phenotypes, and climate in natural populations were maintained or strengthened by selective breeding. However, future adaptive phenotypes and assisted gene flow may be compromised if selective breeding further increases the PEA frequencies of SNPs involved in adaptive trade-offs among climate-related traits.

Evaluating genomic data for management of local adaptation in a changing climate: A lodgepole pine case study.

We evaluate genomic data, relative to phenotypic and climatic data, as a basis for assisted gene flow and genetic conservation. Using a seedling common garden trial of 281 lodgepole pine (Pinus contorta) populations from across western Canada, we compare genomic data to phenotypic and climatic data to assess their effectiveness in characterizing the climatic drivers and spatial scale of local adaptation in this species. We find that phenotype-associated loci are equivalent or slightly superior to climate data for describing local adaptation in seedling traits, but that climate data are superior to genomic data that have not been selected for phenotypic associations. We also find agreement between the climate variables associated with genomic variation and with 20-year heights from a long-term provenance trial, suggesting that genomic data may be a viable option for identifying climatic drivers of local adaptation where phenotypic data are unavailable. Genetic clines associated with the experimental traits occur at broad spatial scales, suggesting that standing variation of adaptive alleles for this and similar species does not require management at scales finer than those indicated by phenotypic data. This study demonstrates that genomic data are most useful when paired with phenotypic data, but can also fill some of the traditional roles of phenotypic data in management of species for which phenotypic trials are not feasible.

Growth gains from selective breeding in a spruce hybrid zone do not compromise local adaptation to climate.

Hybrid zones contain extensive standing genetic variation that facilitates rapid responses to selection. The Picea glauca × Picea engelmannii hybrid zone in western Canada is the focus of tree breeding programs that annually produce ~90 million reforestation seedlings. Understanding the direct and indirect effects of selective breeding on adaptive variation is necessary to implement assisted gene flow (AGF) polices in Alberta and British Columbia that match these seedlings with future climates. We decomposed relationships among hybrid ancestry, adaptive traits, and climate to understand the implications of selective breeding for climate adaptations and AGF strategies. The effects of selection on associations among hybrid index estimated from ~6,500 SNPs, adaptive traits, and provenance climates were assessed for ~2,400 common garden seedlings. Hybrid index differences between natural and selected seedlings within breeding zones were small in Alberta (average +2%), but larger and more variable in BC (average -7%, range -24% to +1%), slightly favoring P. glauca ancestry. The average height growth gain of selected seedlings over natural seedlings within breeding zones was 36% (range 12%-86%). Clines in growth with temperature-related variables were strong, but differed little between selected and natural populations. Seedling hybrid index and growth trait associations with evapotranspiration-related climate variables were stronger in selected than in natural seedlings, indicating possible preadaptation to drier future climates. Associations among cold hardiness, hybrid ancestry, and cold-related climate variables dominated signals of local adaptation and were preserved in breeding populations. Strong hybrid ancestry-phenotype-climate associations suggest that AGF will be necessary to match interior spruce breeding populations with shifting future climates. The absence of antagonistic selection responses among traits and maintenance of cold adaptation in selected seedlings suggests breeding populations can be safely redeployed using AGF prescriptions similar to those of natural populations.