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

Picturing local adaptation: Spectral and structural traits from drone remote sensing reveal clinal responses to climate transfer in common-garden trials of interior spruce (Picea engelmannii × glauca)

Common-garden trials of forest trees provide phenotype data used to assess growth and local adaptation; this information is foundational to tree breeding programs, genecology, and gene conservation. As jurisdictions consider assisted migration strategies to match populations to suitable climates, in situ progeny and provenance trials provide experimental evidence of adaptive responses to climate change. We used drone technology, multispectral imaging, and digital aerial photogrammetry to quantify spectral traits related to stress, photosynthesis, and carotenoids, and structural traits describing crown height, size, and complexity at six climatically disparate common-garden trials of interior spruce (Picea engelmannii × glauca) in western Canada. Through principal component analysis, we identified key components of climate related to temperature, moisture, and elevational gradients. Phenotypic clines in remotely sensed traits were analyzed as trait correlations with provenance climate transfer distances along principal components (PCs). We used traits showing clinal variation to model best linear unbiased predictions for tree height (R2  = .98-.99, root mean square error [RMSE] = 0.06-0.10 m) and diameter at breast height (DBH, R2  = .71-.97, RMSE = 2.57-3.80 mm) and generated multivariate climate transfer functions with the model predictions. Significant (p < .05) clines were present for spectral traits at all sites along all PCs. Spectral traits showed stronger clinal variation than structural traits along temperature and elevational gradients and along moisture gradients at wet, coastal sites, but not at dry, interior sites. Spectral traits may capture patterns of local adaptation to temperature and montane growing seasons which are distinct from moisture-limited patterns in stem growth. This work demonstrates that multispectral indices improve the assessment of local adaptation and that spectral and structural traits from drone remote sensing produce reliable proxies for ground-measured height and DBH. This phenotyping framework contributes to the analysis of common-garden trials towards a mechanistic understanding of local adaptation to climate.

Building resiliency in conifer forests: Interior spruce crosses among weevil resistant and susceptible parents produce hybrids appropriate for multi-trait selection

Tree planting programs now need to consider climate change increasingly, therefore, the resistance to pests plays an essential role in enabling tree adaptation to new ranges through tree population movement. The weevil Pissodes strobi (Peck) is a major pest of spruces and substantially reduces lumber quality. We revisited a large Interior spruce provenance/progeny trial (2,964 genotypes, 42 families) of varying susceptibility, established in British Columbia. We employed multivariate mixed linear models to estimate covariances between, and genetic control of, juvenile height growth and resistance traits. We performed linear regressions and ordinal logistic regressions to test for impact of parental origin on growth and susceptibility to the pest, respectively. A significant environmental component affected the correlations between resistance and height, with outcomes dependent on families. Parents sourced from above 950 m a.s.l. elevation negatively influenced host resistance to attacks, probably due to higher P. engelmannii proportion. For the genetic contribution of parents sourced from above 1,200 m a.s.l., however, we found less attack severity, probably due to a marked mismatch in phenologies. This clearly highlights that interspecific hybrid status might be a good predictor for weevil attacks and delineates the boundaries of successful spruce population movement. Families resulting from crossing susceptible parents generally showed fast-growing trees were the most affected by weevil attacks. Such results indicate that interspecific 'hybrids' with a higher P. glauca ancestry might be genetically better equipped with an optimized resource allocation between defence and growth and might provide the solution for concurrent improvement in resistance against weevil attacks, whilst maintaining tree productivity.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Applying genomics in assisted migration under climate change: Framework, empirical applications, and case studies

The rate of global climate change is projected to outpace the ability of many natural populations and species to adapt. Assisted migration (AM), which is defined as the managed movement of climate-adapted individuals within or outside the species ranges, is a conservation option to improve species' adaptive capacity and facilitate persistence. Although conservation biologists have long been using genetic tools to increase or maintain diversity of natural populations, genomic techniques could add extra benefit in AM that include selectively neutral and adaptive regions of the genome. In this review, we first propose a framework along with detailed procedures to aid collaboration among scientists, agencies, and local and regional managers during the decision-making process of genomics-guided AM. We then summarize the genomic approaches for applying AM, followed by a literature search of existing incorporation of genomics in AM across taxa. Our literature search initially identified 729 publications, but after filtering returned only 50 empirical studies that were either directly applied or considered genomics in AM related to climate change across taxa of plants, terrestrial animals, and aquatic animals; 42 studies were in plants. This demonstrated limited application of genomic methods in AM in organisms other than plants, so we provide further case studies as two examples to demonstrate the negative impact of climate change on non-model species and how genomics could be applied in AM. With the rapidly developing sequencing technology and accumulating genomic data, we expect to see more successful applications of genomics in AM, and more broadly, in the conservation of biodiversity.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Piceaglauca [Moench] Voss), by combining genotype–environment associations, genotype–phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought‐related traits (i.e., anatomical, growth response and climate‐sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse‐controlled conditions. The interlinked phenotype–genotype–environment network revealed eight high‐confidence genes involved in white spruce adaptation to drought, of which four were drought‐responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

An Axiom SNP genotyping array for Douglas-fir

BACKGROUND

In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world.

RESULTS

We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir.

CONCLUSIONS

Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

Assessing the potential for assisted gene flow using past introduction of Norway spruce in southern Sweden: Local adaptation and genetic basis of quantitative traits in trees

Norway spruce (Picea abies) is a dominant conifer species of major economic importance in northern Europe. Extensive breeding programs were established to improve phenotypic traits of economic interest. In southern Sweden, seeds used to create progeny tests were collected on about 3,000 trees of outstanding phenotype ('plus' trees) across the region. In a companion paper, we showed that some were of local origin but many were recent introductions from the rest of the natural range. The mixed origin of the trees together with partial sequencing of the exome of >1,500 of these trees and phenotypic data retrieved from the Swedish breeding program offered a unique opportunity to dissect the genetic basis of local adaptation of three quantitative traits (height, diameter and bud-burst) and assess the potential of assisted gene flow. Through a combination of multivariate analyses and genome-wide association studies, we showed that there was a very strong effect of geographical origin on growth (height and diameter) and phenology (bud-burst) with trees from southern origins outperforming local provenances. Association studies revealed that growth traits were highly polygenic and bud-burst somewhat less. Hence, our results suggest that assisted gene flow and genomic selection approaches could help to alleviate the effect of climate change on P. abies breeding programs in Sweden.

Purifying selection does not drive signatures of convergent local adaptation of lodgepole pine and interior spruce

BACKGROUND

Lodgepole pine (Pinus contorta) and interior spruce (Picea glauca, Picea engelmannii, and their hybrids) are distantly related conifer species. Previous studies identified 47 genes containing variants associated with environmental variables in both species, providing evidence of convergent local adaptation. However, if the intensity of purifying selection varies with the environment, clines in nucleotide diversity could evolve through linked (background) selection that would yield allele frequency-environment signatures resembling local adaptation. If similar geographic patterns in the strength of purifying selection occur in these species, this could result in the convergent signatures of local adaptation, especially if the landscape of recombination is conserved. In the present study, we investigated whether spatially/environmentally varying purifying selection could give rise to the convergent signatures of local adaptation that had previously reported.

RESULTS

We analyzed 86 lodgepole pine and 50 interior spruce natural populations spanning heterogeneous environments in western Canada where previous analyses had found signatures of convergent local adaptation. We estimated nucleotide diversity and Tajima's D for each gene within each population and calculated the strength of correlations between nucleotide diversity and environmental variables. Overall, these estimates in the genes with previously identified convergent local adaptation signatures had no similar pattern between pine and spruce. Clines in nucleotide diversity along environmental variables were found for interior spruce, but not for lodgepole pine. In spruce, genes with convergent adaption signatures showed a higher strength of correlations than genes without convergent adaption signatures, but there was no such disparity in pine, which suggests the pattern in spruce may have arisen due to a combination of selection and hybridization.

CONCLUSIONS

The results rule out purifying/background selection as a driver of convergent local adaption signatures in lodgepole pine and interior spruce.

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.

Ecophysiology and Growth of White Spruce Seedlings from Various Seed Sources along a Climatic Gradient Support the Need for Assisted Migration

With climate change, favorable growing conditions for tree species are shifting northwards and to higher altitudes. Therefore, local populations are becoming less adapted to their environment. Assisted migration is one of the proposed adaptive measures to reduce the vulnerability of natural populations and maintain forest productivity. It consists of moving genetic material to a territory where future climate conditions correspond to those of its current location. Eight white spruce (Picea glauca [Moench] Voss) seed sources representing as many seed orchards were planted in 2013 at three forest sites simulating a south-north climatic gradient of 1.7°C in Québec, Canada. The objectives were to (1) evaluate the morpho-physiological responses of the different seed sources and (2) determine the role of genetic adaptation and physiological plasticity on the observed variation in morpho-physiological traits. Various seedling characteristics were measured, notably height growth from nursery to the fourth year on plantation. Other traits such as biomass and carbon allocation, nutritional status, and various photosynthetic traits before bud break, were evaluated during the fourth growing season. No interaction between sites and seed sources was observed for any traits, suggesting similar plasticity between seed sources. There was no change in the rank of seed sources and sites between years for height growth. Moreover, a significant positive correlation was observed between the height from the nursery and that after 4 years in the plantation. Southern seed sources showed the best height growth, while optimum growth was observed at the central site. Juvenile height growth seems to be a good indicator of the juvenile carbon sequestration and could serve as a selection criterion for the best genetics sources for carbon sequestration. Vector analysis showed no nitrogen deficiency 4 years after planting. Neither seed sources nor planting sites had a significant effect on photosynthesis before bud break. The observed results during the establishment phase under different site conditions indicate that southern seed sources may already benefit from assisted migration to cooler climatic conditions further north. While northern seed sources are likely to benefit from anticipated local global warming, they would not match the growth performance of seedlings from southern sources.