Complementary roles of two classes of defense chemicals in white spruce against spruce budworm.

MAIN CONCLUSION: Monoterpenes and phenolics play distinct roles in defending white spruce trees from insect defoliators. Monoterpenes contribute to the toxicity of the foliage, deterring herbivory, whereas phenolics impede budworm growth. This study demonstrates the complex interplay between monoterpenes and phenolics and their collective influence on the defense strategy of white spruce trees against a common insect defoliator. Long-lived coniferous trees display considerable variations in their defensive chemistry. The impact of these defense phenotype variations on insect herbivores of the same species remains to be thoroughly studied, mainly due to challenges in replicating the comprehensive defense profiles of trees under controlled conditions. This study methodically examined the defensive properties of foliar monoterpenes and phenolics across 80 distinct white spruce families. These families were subsequently grouped into two chemotypes based on their foliar monoterpene concentrations. To understand the separate and combined effects of these classes on tree defenses to the eastern spruce budworm, we conducted feeding experiments using actual defense profiles from representative families. Specifically, we assessed budworm response when exposed to substrates amended with phenolics alone or monoterpenes. Our findings indicate that the ratios and amounts of monoterpenes and phenolics present in the white spruce foliage influence the survival of spruce budworms. Phenotypes associated with complete larval mortality exhibited elevated ratios (ranging from 0.4 to 0.6) and concentrations (ranging from 1143 to 1796 ng mg-1) of monoterpenes. Conversely, families characterized by higher phenolic ratios (ranging from 0.62 to 0.77) and lower monoterpene concentrations (ranging from 419 to 985 ng mg-1) were less lethal to the spruce budworm. Both classes of defense compounds contribute significantly to the overall defensive capabilities of white spruce trees. Monoterpenes appear critical in determining the general toxicity of foliage, while phenolics play a role in slowing budworm development, thereby underscoring their collective importance in white spruce defenses.

Hybridization mediated range expansion and climate change resilience in two keystone tree species of boreal forests.

Current global climate change is expected to affect biodiversity negatively at all scales leading to mass biodiversity loss. Many studies have shown that the distribution of allele frequencies across a species' range is often influenced by specific genetic loci associated with local environmental variables. This association reflects local adaptation and allele changes at those loci could thereby contribute to the evolutionary response to climate change. However, predicting how species will adapt to climate change from this type of data alone remains challenging. In the present study, we combined exome capture sequences and environmental niche reconstruction, to test multiple methods for assessing local adaptation and climate resilience in two widely distributed conifers, Norway spruce and Siberian spruce. Both species are keystone species of the boreal forest and share a vast hybrid zone. We show that local adaptation in conifers can be detected through allele frequency variation, population-level ecological preferences, and historical niche movement. Moreover, we integrated genetic and ecological information into genetic offset predictive models to show that hybridization plays a central role in expanding the niche breadth of the two conifer species and may help both species to cope better with future changing climates. This joint genetic and ecological analysis also identified spruce populations that are at risk under current climate change.

Picea wilsonii NAC31 and DREB2A Cooperatively Activate ERD1 to Modulate Drought Resistance in Transgenic Arabidopsis.

The NAC family of transcription factors (TFs) regulate plant development and abiotic stress. However, the specific function and response mechanism of NAC TFs that increase drought resistance in Picea wilsonii remain largely unknown. In this study, we functionally characterized a member of the PwNAC family known as PwNAC31. PwNAC31 is a nuclear-localized protein with transcriptional activation activity and contains an NAC domain that shows extensive homology with ANAC072 in Arabidopsis. The expression level of PwNAC31 is significantly upregulated under drought and ABA treatments. The heterologous expression of PwNAC31 in atnac072 Arabidopsis mutants enhances the seed vigor and germination rates and restores the hypersensitive phenotype of atnac072 under drought stress, accompanied by the up-regulated expression of drought-responsive genes such as DREB2A (DEHYDRATION-RESPONSIVE ELEMENT BINDING PROTEIN 2A) and ERD1 (EARLY RESPONSIVE TO DEHYDRATION STRESS 1). Yeast two-hybrid and bimolecular fluorescence complementation assays confirmed that PwNAC31 interacts with DREB2A and ABF3 (ABSCISIC ACID-RESPONSIVE ELEMENT-BINDING FACTOR 3). Yeast one-hybrid and dual-luciferase assays showed that PwNAC31, together with its interaction protein DREB2A, directly regulated the expression of ERD1 by binding to the DRE element of the ERD1 promoter. Collectively, our study provides evidence that PwNAC31 activates ERD1 by interacting with DREB2A to enhance drought tolerance in transgenic Arabidopsis.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers.

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene ßglu-1 with large variations in white spruce but not Norway spruce.

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and ßglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of ßglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of ßglu-1 and Ugt5 genes. We observed very large copy numbers of ßglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of ßglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce.

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

A conserved graft formation process in Norway spruce and Arabidopsis identifies the PAT gene family as central regulators of wound healing.

The widespread use of plant grafting enables eudicots and gymnosperms to join with closely related species and grow as one. Gymnosperms have dominated forests for over 200 million years, and despite their economic and ecological relevance, we know little about how they graft. Here we developed a micrografting method in conifers using young tissues that allowed efficient grafting with closely related species and between distantly related genera. Conifer graft junctions rapidly connected vasculature and differentially expressed thousands of genes including auxin and cell-wall-related genes. By comparing these genes to those induced during Arabidopsis thaliana graft formation, we found a common activation of cambium, cell division, phloem and xylem-related genes. A gene regulatory network analysis in Norway spruce (Picea abies) predicted that PHYTOCHROME A SIGNAL TRANSDUCTION 1 (PAT1) acted as a core regulator of graft healing. This gene was strongly up-regulated during both spruce and Arabidopsis grafting, and Arabidopsis mutants lacking PAT genes failed to attach tissues or successfully graft. Complementing Arabidopsis PAT mutants with the spruce PAT1 homolog rescued tissue attachment and enhanced callus formation. Together, our data show an ability for young tissues to graft with distantly related species and identifies the PAT gene family as conserved regulators of graft healing and tissue regeneration.

Genomic signals of local adaptation in Picea crassifolia.

BACKGROUND: Global climate change poses a grave threat to biodiversity and underscores the importance of identifying the genes and corresponding environmental factors involved in the adaptation of tree species for the purposes of conservation and forestry. This holds particularly true for spruce species, given their pivotal role as key constituents of the montane, boreal, and sub-alpine forests in the Northern Hemisphere. RESULTS: Here, we used transcriptomes, species occurrence records, and environmental data to investigate the spatial genetic distribution of and the climate-associated genetic variation in Picea crassifolia. Our comprehensive analysis employing ADMIXTURE, principal component analysis (PCA) and phylogenetic methodologies showed that the species has a complex population structure with obvious differentiation among populations in different regions. Concurrently, our investigations into isolation by distance (IBD), isolation by environment (IBE), and niche differentiation among populations collectively suggests that local adaptations are driven by environmental heterogeneity. By integrating population genomics and environmental data using redundancy analysis (RDA), we identified a set of climate-associated single-nucleotide polymorphisms (SNPs) and showed that environmental isolation had a more significant impact than geographic isolation in promoting genetic differentiation. We also found that the candidate genes associated with altitude, temperature seasonality (Bio4) and precipitation in the wettest month (Bio13) may be useful for forest tree breeding. CONCLUSIONS: Our findings deepen our understanding of how species respond to climate change and highlight the importance of integrating genomic and environmental data in untangling local adaptations.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response.

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Habitat and climate influence hybridization among three genetically distinct Canada jay (Perisoreus canadensis) morphotypes in an avian hybrid zone complex.

Examining the frequency and distribution of hybrids across contact zones provide insights into the factors mediating hybridization. In this study, we examined the effect of habitat and climate on hybridization patterns for three phenotypically, genetically, and ecologically distinct groups of the Canada jay (Perisoreus canadensis) in a secondary contact zone in western North America. Additionally, we tested whether the frequency of hybridization involving the three groups (referred to as Boreal, Pacific and Rocky Mountain morphotypes) is similar across the hybrid zones or whether some pairs have hybridized more frequently than others. We reanalyzed microsatellite, mtDNA and plumage data, and new microsatellite and plumage data for 526 individuals to identify putative genetic and phenotypic hybrids. The genetically and phenotypically distinct groups are associated with different habitats and occupy distinct climate niches across the contact zone. Most putative genetic hybrids (86%) had Rocky Mountain ancestry. Hybrids were observed most commonly in intermediate climate niches and in habitats where Engelmann spruce (Picea engelmannii) overlaps broadly with boreal and subalpine tree species. Our finding that hybrids occupy intermediate climate niches relative to parental morphotypes matches patterns for other plant and animal species found in this region. This study demonstrates how habitat and climate influence hybridization patterns in areas of secondary contact and adds to the growing body of research on tri-species hybrid zones.

Quantifying adaptive evolution and the effects of natural selection across the Norway spruce genome.

Detecting natural selection is one of the major goals of evolutionary genomics. Here, we sequenced the whole genome of 25 Picea abies individuals and quantified the amount of selection across the genome. Using an estimate of the distribution of fitness effects, we showed that both negative selection and the rate of positively selected substitutions are very limited in coding regions. We found a positive correlation between the rate of adaptive substitutions and recombination rate and a negative correlation between the rate of adaptive substitutions and gene density, suggesting a widespread influence from Hill-Robertson interference on the efficiency of protein adaptation in P. abies. Finally, the distinct population statistics between genomic regions under either positive or balancing selection with that under neutral regions indicated the impact of natural selection on the genomic architecture of Norway spruce. Further gene ontology enrichment analysis for genes located in regions identified as undergoing either positive or long-term balancing selection also highlighted the specific molecular functions and biological processes that appear to be targets of selection in Norway spruce.

Effects of plant growth regulators and sucrose on proliferation and quality of embryogenic tissue in Picea pungens.

Embryogenic tissue (ET) is important for genetic modification and plant re-generation. The proliferation ability and vigor of ET are crucial for plant propagation via somatic embryogenesis. In this study, ET was induced from mature zygotic embryos in blue spruce (Picea pungens Engelm.). There were significant differences in ET induction between two provenances, i.e. 78.8 ± 12.5% and 62.50 ± 12.8% respectively. Effects of 2,4-Dichlorophenoxy acetic acid (2,4-D), 6-Benzyl amino-purine (6-BA) and/or sucrose on ET proliferation and somatic embryo (SE) maturation were further investigated with four cell lines. The highest ET proliferation rate reached 1473.7 ± 556.0% biweekly. Concentrations of 2,4-D or 6-BA applied at tissue proliferation stage impacted SE maturation among the cell lines, whereas sucrose showed less effects. The highest rate, 408 ± 230 mature SEs/g FW, was achieved in SE maturation cultures. This research demonstrated that the culture conditions, i.e. the specific concentrations of 2,4-D and BA, at ET proliferation stage affected not only ET growth, but also the quality of ET for SE maturation. This study revealed the necessity and benefit in developing both the general and the genotype-specific protocols for efficient production of mature SEs, or somatic plants in blue spruce.

Transcriptome features of stone cell development in weevil-resistant and susceptible Sitka spruce.

Stone cells are a specialized, highly lignified cell type found in both angiosperms and gymnosperms. In conifers, abundance of stone cells in the cortex provides a robust constitutive physical defense against stem feeding insects. Stone cells are a major insect-resistance trait in Sitka spruce (Picea sitchensis), occurring in dense clusters in apical shoots of trees resistant (R) to spruce weevil (Pissodes strobi) but being rare in susceptible (S) trees. To learn more about molecular mechanisms of stone cell formation in conifers, we used laser microdissection and RNA sequencing to develop cell-type-specific transcriptomes of developing stone cells from R and S trees. Using light, immunohistochemical, and fluorescence microscopy, we also visualized the deposition of cellulose, xylan, and lignin associated with stone cell development. A total of 1293 genes were differentially expressed at higher levels in developing stone cells relative to cortical parenchyma. Genes with potential roles in stone cell secondary cell wall formation (SCW) were identified and their expression evaluated over a time course of stone cell formation in R and S trees. The expression of several transcriptional regulators was associated with stone cell formation, including a NAC family transcription factor and several genes annotated as MYB transcription factors with known roles in SCW formation.

Clinal variation in PHY (PAS domain) and CRY (CCT domain)-Signs of local adaptation to light quality in Norway spruce.

Detection of the genomic basis of local adaptation to environmental conditions is challenging in forest trees. Phytochromes (PHY) and cryptochromes (CRY) perceive the red (R)/far-red (FR) and blue light respectively, thus playing a fundamental role in regulating plant growth and development. PHYO and PHYP from conifers are the equivalents of PHYA/PHYC and PHYB in angiosperms, respectively. Norway spruce shows an adaptive latitudinal cline for shade (low R:FR or FR-enriched light) tolerance and requirement of FR light for its growth. We analyzed the exome capture data that included a uniquely large data set of 1654 Norway spruce trees sampled across many latitudes in Sweden to capture the natural clines for photoperiod and FR light exposure during the growth season. Statistically significant clinal variation was detected in allele and genotype frequencies of missense mutations in coding regions belonging to well-defined functional domains of PHYO (PAS-B), PHYP2 (PAS fold-2), CRY1 (CCT1) and CRY2 (CCT2) that strongly correlates with the latitudinal gradient in response to variable light quality in Norway spruce. The missense SNP in PHYO resulting in Asn835Ser, displayed the steepest cline among all other polymorphisms. We propose that these variations in the photoreceptors represent signs of local adaptation to light quality.

Backcrossing to different parents produced two distinct hybrid species.

Repeated homoploid hybrid speciation (HHS) events with the same parental species have rarely been reported. In this study, we used population transcriptome data to test paraphyly and HHS events in the conifer Picea brachytyla. Our analyses revealed non-sister relationships for two lineages of P. brachytyla, with the southern lineage being placed within the re-circumscribed P. likiangensis species complex (PLSC) and P. brachytyla sensu stricto (s.s.) consisted solely of the northern lineage, forming a distinct clade that is paratactic to both the PLSC and P. wilsonii. Our phylogenetic and coalescent analyses suggested that P. brachytyla s.s. arose from HHS between the ancestor of the PLSC before its diversification and P. wilsonii through an intermediate hybrid lineage at an early stage and backcrossing to the ancestral PLSC. Additionally, P. purpurea shares the same parents and an extinct lineage with P. brachytyla s.s. but backcrossing to the other parent, P. wilsonii at a later stage. We reveal the first case that backcrossing to different parents of the same extinct hybrid lineage produced two different hybrid species. Our results highlight the existence of more reticulate evolution during species diversification in the spruce genus and more complex homoploid hybrid events than previously identified.

Analysis of polycomb repressive complex 2 (PRC2) subunits in Picea abies with a focus on embryo development.

BACKGROUND: Conserved polycomb repressive complex 2 (PRC2) mediates H3K27me3 to direct transcriptional repression and has a key role in cell fate determination and cell differentiation in both animals and plants. PRC2 subunits have undergone independent multiplication and functional divergence in higher plants. However, relevant information is still absent in gymnosperms. RESULTS: To launch gymnosperm PRC2 research, we identified and cloned the PRC2 core component genes in the conifer model species Picea abies, including one Esc/FIE homolog PaFIE, two p55/MSI homologs PaMSI1a and PaMSI1b, two E(z) homologs PaKMT6A2 and PaKMT6A4, a Su(z)12 homolog PaEMF2 and a PaEMF2-like fragment. Phylogenetic and protein domain analyses were conducted. The Esc/FIE homologs were highly conserved in the land plant, except the monocots. The other gymnospermous PRC2 subunits underwent independent evolution with angiospermous species to different extents. The relative transcript levels of these genes were measured in endosperm and zygotic and somatic embryos at different developmental stages. The obtained results proposed the involvement of PaMSI1b and PaKMT6A4 in embryogenesis and PaKMT6A2 and PaEMF2 in the transition from embryos to seedlings. The PaEMF2-like fragment was predominantly expressed in the endosperm but not in the embryo. In addition, immunohistochemistry assay showed that H3K27me3 deposits were generally enriched at meristem regions during seed development in P. abies. CONCLUSIONS: This study reports the first characterization of the PRC2 core component genes in the coniferous species P. abies. Our work may enable a deeper understanding of the cell reprogramming process during seed and embryo development and may guide further research on embryonic potential and development in conifers.

Effects of combined drought and pathogen stress on growth, resistance and gene expression in young Norway spruce trees.

Drought-induced mortality is a major direct effect of climate change on tree health, but drought can also affect trees indirectly by altering their susceptibility to pathogens. Here, we report how a combination of mild or severe drought and pathogen infection affected the growth, pathogen resistance and gene expression in potted 5-year-old Norway spruce trees [Picea abies (L.) Karst.]. After 5 weeks of drought, trees were inoculated with the fungal pathogen Endoconidiophora polonica. Combined drought-pathogen stress over the next 8 weeks led to significant reductions in the growth of drought-treated trees relative to well-watered trees and more so in trees subjected to severe drought. Belowground, growth of the smallest fine roots was most affected. Aboveground, shoot diameter change was most sensitive to the combined stress, followed by shoot length growth and twig biomass. Both drought-related and some resistance-related genes were upregulated in bark samples collected after 5 weeks of drought (but before pathogen infection), and gene expression levels scaled with the intensity of drought stress. Trees subjected to severe drought were much more susceptible to pathogen infection than well-watered trees or trees subjected to mild drought. Overall, our results show that mild drought stress may increase the tree resistance to pathogen infection by upregulating resistance-related genes. Severe drought stress, on the other hand, decreased tree resistance. Because drought episodes are expected to become more frequent with climate change, combined effects of drought and pathogen stress should be studied in more detail to understand how these stressors interactively influence tree susceptibility to pests and pathogens.

Preselection of QTL markers enhances accuracy of genomic selection in Norway spruce.

Genomic prediction (GP) or genomic selection is a method to predict the accumulative effect of all quantitative trait loci (QTLs) in a population by estimating the realized genomic relationships between the individuals and by capturing the linkage disequilibrium between markers and QTLs. Thus, marker preselection is considered a promising method to capture Mendelian segregation effects. Using QTLs detected in a genome-wide association study (GWAS) may improve GP. Here, we performed GWAS and GP in a population with 904 clones from 32 full-sib families using a newly developed 50 k SNP Norway spruce array. Through GWAS we identified 41 SNPs associated with budburst stage (BB) and the largest effect association explained 5.1% of the phenotypic variation (PVE). For the other five traits such as growth and wood quality traits, only 2 - 13 associations were observed and the PVE of the strongest effects ranged from 1.2% to 2.0%. GP using approximately 100 preselected SNPs, based on the smallest p-values from GWAS showed the greatest predictive ability (PA) for the trait BB. For the other traits, a preselection of 2000-4000 SNPs, was found to offer the best model fit according to the Akaike information criterion being minimized. But PA-magnitudes from GP using such selections were still similar to that of GP using all markers. Analyses on both real-life and simulated data also showed that the inclusion of a large QTL SNP in the model as a fixed effect could improve PA and accuracy of GP provided that the PVE of the QTL was ≥ 2.5%.

Transcriptome-based analyses of adaptive divergence between two closely related spruce species on the Qinghai-Tibet plateau and adjacent regions.

Speciation among populations connected by gene flow is driven by adaptation to different environments, but underlying gene-environment associations remain largely unknown. Here, 162 individuals from 32 populations were sampled to obtain 191,648 independent single nucleotide polymorphisms (SNPs) across the genomes of two closely related spruce species, Picea asperata and Picea crassifolia, which occur on the Qinghai-Tibet Plateau and in surrounding regions. Using the SNP data set, genotype-environment associations and demographic modelling were used to examine local adaptation and genetic divergence between these two species. While morphologically similar, the two Picea species were genetically differentiated in multiple analyses. These species diverged despite continuous gene flow, and their initial divergence was dated back to the late Quaternary. The effective population sizes of both species have expanded since their divergence, as confirmed by niche distribution simulations. A total of 6365 genes were associated with the tested environmental variables; of these, 41 were positively selected in P. asperata and were mainly associated with temperature, while 83 were positively selected in P. crassifolia and were primarily associated with precipitation. These results deepen our understanding of the adaptive divergence and demographic histories of these two spruce species and highlight the importance of genomic data in deciphering the environmental selection underlying Quaternary interspecific divergence.

From common gardens to candidate genes: exploring local adaptation to climate in red spruce.

Local adaptation to climate is common in plant species and has been studied in a range of contexts, from improving crop yields to predicting population maladaptation to future conditions. The genomic era has brought new tools to study this process, which was historically explored through common garden experiments. In this study, we combine genomic methods and common gardens to investigate local adaptation in red spruce and identify environmental gradients and loci involved in climate adaptation. We first use climate transfer functions to estimate the impact of climate change on seedling performance in three common gardens. We then explore the use of multivariate gene-environment association methods to identify genes underlying climate adaptation, with particular attention to the implications of conducting genome scans with and without correction for neutral population structure. This integrative approach uncovered phenotypic evidence of local adaptation to climate and identified a set of putatively adaptive genes, some of which are involved in three main adaptive pathways found in other temperate and boreal coniferous species: drought tolerance, cold hardiness, and phenology. These putatively adaptive genes segregated into two 'modules' associated with different environmental gradients. This study nicely exemplifies the multivariate dimension of adaptation to climate in trees.