A resource of identified and annotated lincRNAs expressed during somatic embryogenesis development in Norway spruce.

Long non-coding RNAs (lncRNAs) have emerged as important regulators of many biological processes, although their regulatory roles remain poorly characterized in woody plants, especially in gymnosperms. A major challenge of working with lncRNAs is to assign functional annotations, since they have a low coding potential and low cross-species conservation. We utilised an existing RNA-Sequencing resource and performed short RNA sequencing of somatic embryogenesis developmental stages in Norway spruce (Picea abies L. Karst). We implemented a pipeline to identify lncRNAs located within the intergenic space (lincRNAs) and generated a co-expression network including protein coding, lincRNA and miRNA genes. To assign putative functional annotation, we employed a guilt-by-association approach using the co-expression network and integrated these results with annotation assigned using semantic similarity and co-expression. Moreover, we evaluated the relationship between lincRNAs and miRNAs, and identified which lincRNAs are conserved in other species. We identified lincRNAs with clear evidence of differential expression during somatic embryogenesis and used network connectivity to identify those with the greatest regulatory potential. This work provides the most comprehensive view of lincRNAs in Norway spruce and is the first study to perform global identification of lincRNAs during somatic embryogenesis in conifers. The data have been integrated into the expression visualisation tools at the PlantGenIE.org web resource to enable easy access to the community. This will facilitate the use of the data to address novel questions about the role of lincRNAs in the regulation of embryogenesis and facilitate future comparative genomics studies.

Insights into the Detoxification of Spruce Monoterpenes by the Eurasian Spruce Bark Beetle.

Plant defence mechanisms, including physical barriers like toughened bark and chemical defences like allelochemicals, are essential for protecting them against pests. Trees allocate non-structural carbohydrates (NSCs) to produce secondary metabolites like monoterpenes, which increase during biotic stress to fend off pests like the Eurasian spruce bark beetle, ESBB (Ips typographus). Despite these defences, the ESBB infests Norway spruce, causing significant ecological damage by exploiting weakened trees and using pheromones for aggregation. However, the mechanism of sensing and resistance towards host allelochemicals in ESBB is poorly understood. We hypothesised that the exposure of ESBB to spruce allelochemicals, especially monoterpenes, leads to an upsurge in the important detoxification genes like P450s, GSTs, UGTs, and transporters, and at the same time, genes responsible for development must be compromised. The current study demonstrates that exposure to monoterpenes like R-limonene and sabiene effectively elevated detoxification enzyme activities. The differential gene expression (DGE) analysis revealed 294 differentially expressed (DE) detoxification genes in response to R-limonene and 426 DE detoxification genes in response to sabiene treatments, with 209 common genes between the treatments. Amongst these, genes from the cytochrome P450 family 4 and 6 genes (CP4 and CP6), esterases, glutathione S-transferases family 1 (GSTT1), UDP-glucuronosyltransferase 2B genes (UDB), and glucose synthesis-related dehydrogenases were highly upregulated. We further validated 19 genes using RT-qPCR. Additionally, we observed similar high expression levels of detoxification genes across different monoterpene treatments, including myrcene and α-pinene, suggesting a conserved detoxification mechanism in ESBB, which demands further investigation. These findings highlight the potential for molecular target-based beetle management strategies targeting these key detoxification genes.

Evolutionary drivers of reproductive fitness in two endangered forest trees.

Population genetics theory predicts a relationship between fitness, genetic diversity (H0) and effective population size (Ne), which is often tested through heterozygosity-fitness correlations (HFCs). We tested whether population and individual fertility and heterozygosity are correlated in two endangered Mexican spruces (Picea martinezii and Picea mexicana) by combining genomic, demographic and reproductive data (seed development and germination traits). For both species, there was a positive correlation between population size and seed development traits, but not germination rate. Individual genome-wide heterozygosity and seed traits were only correlated in P. martinezii (general-effects HFC), and none of the candidate single nucleotide polymorphisms (SNPs) associated with individual fertility showed heterozygote advantage in any species (no local-effects HFC). We observed a single and recent (c. 30 thousand years ago (ka)) population decline for P. martinezii; the collapse of P. mexicana occurred in two phases separated by a long period of stability (c. 800 ka). Recruitment always contributed more to total population census than adult trees in P. mexicana, while this was only the case in the largest populations of P. martinezii. Equating fitness to either H0 or Ne, as traditionally proposed in conservation biology, might not always be adequate, as species-specific evolutionary factors can decouple the expected correlation between these parameters.

Phenotypic variation in 1,100 provenances of Picea abies measured over 50 years on 33 German trial sites.

We present a database of Norway spruce phenotypic traits measured over a period of 51 years in provenance trials planted across western Germany. These trials are part of the 'Inventory Provenance Test 1964/68 with Norway spruce' (IPTNS) that tested a total of 1,100 provenances collected across Europe, both from the species natural and artificial range. Phenotypic traits (first height, later diameter at breast height, survival) were recorded at 33 trial sites of 100 provenances, each a subsample from the entire collection area that represents an approximately comparable climate space. The dataset contains 424,781 records of the same 65,518 trees ranging from age 7 to 51. Overall, it captures the considerable genetic and phenotypic intraspecific variation present in Norway spruce and should be of interest to various fields including quantitative genetics, ecology, biogeography and the adaptive management of forests.

Recurrent hybridization and gene flow shaped Norway and Siberian spruce evolutionary history over multiple glacial cycles.

Most tree species underwent cycles of contraction and expansion during the Quaternary. These cycles led to an ancient and complex genetic structure that has since been affected by extensive gene flow and by strong local adaptation. The extent to which hybridization played a role in this multi-layered genetic structure is important to be investigated. To study the effect of hybridization on the joint population genetic structure of two dominant species of the Eurasian boreal forest, Picea abies and P. obovata, we used targeted resequencing and obtained around 480 K nuclear SNPs and 87 chloroplast SNPs in 542 individuals sampled across most of their distribution ranges. Despite extensive gene flow and a clear pattern of Isolation-by-Distance, distinct genetic clusters emerged, indicating the presence of barriers and corridors to migration. Two cryptic refugia located in the large hybrid zone between the two species played a critical role in shaping their current distributions. The two species repeatedly hybridized during the Pleistocene and the direction of introgression depended on latitude. Our study suggests that hybridization helped both species to overcome main shifts in their distribution ranges during glacial cycles and highlights the importance of considering whole species complex instead of separate entities to retrieve complex demographic histories.

PaMYB11 promotes suberin deposition in Norway spruce embryogenic tissue during cryopreservation: A novel resistance mechanism against osmosis.

The osmotic resistance mechanism has been extensively studied in whole plants or plant tissues. However, little is known about it in embryogenic tissue (ET) which is widely used in plant-based biotechnological systems. Suberin, a cell wall aliphatic and aromatic heteropolymer, plays a critical role in plant cells against osmosis stress. The suberin regulatory biosynthesis has rarely been studied in gymnosperms. Here, PaMYB11, a subgroup 11 R2R3-MYB transcription factor, plays a key role in the osmotic resistance of Norway spruce (Picea abies) ETs during cryoprotectant pretreatment. Thus, RNA-seq, histological, and analytical chemical analyses are performed on the stable transformations of PaMYB11-OE and PaMYB11-SRDX in Norway spruce ETs. DAP-seq, Y1H, and LUC are further combined to explore the PaMYB11 targets. Activation of PaMYB11 is necessary and sufficient for suberin lamellae deposition on Norway spruce embryogenic cell walls, which plays a decisive role in ET survival under osmotic stress. Transcriptome analysis shows that PaMYB11 enhances suberin lamellae monomer synthesis by promoting very long-chain fatty acid (VLCFA) synthesis. PaPOP, PaADH1, and PaTET8L, the first two (PaADH1 and PaPOP, included) involved in VLCFA synthesis, are proved to be the direct targets of PaMYB11. Our study identified a novel osmotic response directed by PaMYB11 in Norway spruce ET, which provides a new understanding of the resistance mechanism against osmosis in gymnosperms.

Future carbon sequestration potential in a widespread transcontinental boreal tree species: Standing genetic variation matters!

Climate change (CC) necessitates reforestation/afforestation programs to mitigate its impacts and maximize carbon sequestration. But comprehending how tree growth, a proxy for fitness and resilience, responds to CC is critical to maximize these programs' effectiveness. Variability in tree response to CC across populations can notably be influenced by the standing genetic variation encompassing both neutral and adaptive genetic diversity. Here, a framework is proposed to assess tree growth potential at the population scale while accounting for standing genetic variation. We applied this framework to black spruce (BS, Picea mariana [Mill] B.S.P.), with the objectives to (1) determine the key climate variables having impacted BS growth response from 1974 to 2019, (2) examine the relative roles of local adaptation and the phylogeographic structure in this response, and (3) project BS growth under two Shared Socioeconomic Pathways while taking standing genetic variation into account. We modeled growth using a machine learning algorithm trained with dendroecological and genetic data obtained from over 2600 trees (62 populations divided in three genetic clusters) in four 48-year-old common gardens, and simulated growth until year 2100 at the common garden locations. Our study revealed that high summer and autumn temperatures negatively impacted BS growth. As a consequence of warming, this species is projected to experience a decline in growth by the end of the century, suggesting maladaptation to anticipated CC and a potential threat to its carbon sequestration capacity. This being said, we observed a clear difference in response to CC within and among genetic clusters, with the western cluster being more impacted than the central and eastern clusters. Our results show that intraspecific genetic variation, notably associated with the phylogeographic structure, must be considered when estimating the response of widespread species to CC.

Implications of accounting for marker-based population structure in the quantitative genetic evaluation of genetic parameters related to growth and wood properties in Norway spruce.

BACKGROUND: Forest geneticists typically use provenances to account for population differences in their improvement schemes; however, the historical records of the imported materials might not be very precise or well-aligned with the genetic clusters derived from advanced molecular techniques. The main objective of this study was to assess the impact of marker-based population structure on genetic parameter estimates related to growth and wood properties and their trade-offs in Norway spruce, by either incorporating it as a fixed effect (model-A) or excluding it entirely from the analysis (model-B). RESULTS: Our results indicate that models incorporating population structure significantly reduce estimates of additive genetic variance, resulting in substantial reduction of narrow-sense heritability. However, these models considerably improve prediction accuracies. This was particularly significant for growth and solid-wood properties, which showed to have the highest population genetic differentiation (QST) among the studied traits. Additionally, although the pattern of correlations remained similar across the models, their magnitude was slightly lower for models that included population structure as a fixed effect. This suggests that selection, consistently performed within populations, might be less affected by unfavourable genetic correlations compared to mass selection conducted without pedigree restrictions. CONCLUSION: We conclude that the results of models properly accounting for population structure are more accurate and less biased compared to those neglecting this effect. This might have practical implications for breeders and forest managers where, decisions based on imprecise selections can pose a high risk to economic efficiency.

Extracellular vesicles of Norway spruce contain precursors and enzymes for lignin formation and salicylic acid.

Lignin is a phenolic polymer in plants that rigidifies the cell walls of water-conducting tracheary elements and support-providing fibers and stone cells. Different mechanisms have been suggested for the transport of lignin precursors to the site of lignification in the cell wall. Extracellular vesicle (EV)-enriched samples isolated from a lignin-forming cell suspension culture of Norway spruce (Picea abies L. Karst.) contained both phenolic metabolites and enzymes related to lignin biosynthesis. Metabolomic analysis revealed mono-, di-, and oligolignols in the EV isolates, as well as carbohydrates and amino acids. In addition, salicylic acid (SA) and some proteins involved in SA signaling were detected in the EV-enriched samples. A proteomic analysis detected several laccases, peroxidases, ß-glucosidases, putative dirigent proteins, and cell wall-modifying enzymes, such as glycosyl hydrolases, transglucosylase/hydrolases, and expansins in EVs. Our findings suggest that EVs are involved in transporting enzymes required for lignin polymerization in Norway spruce, and radical coupling of monolignols can occur in these vesicles.

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level.

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Circadian rhythm and redox homeostasis candidate genes showed association with shallow elevation in Norway spruce.

The analysis of genetic variation underlying local adaptation in natural populations, together with the response to different external stimuli, is currently a hot topic in forest sciences, with the aim of identifying genetic markers controlling key phenotypic traits of interest for their inclusion in restoration and breeding programs. In Europe, one of the main tree species is Norway spruce (Picea abies (L.) H.Karst.). Using the MassARRAY® platform, 568 trees from North Rhine-Westphalia (Germany) were genotyped with 94 single nucleotide polymorphisms (SNPs) related to circadian and growth rhythms, and to stress response. The association analysis of the selected markers with health status and elevation was performed using three different methods, and those identified by at least two of these were considered as high confidence associated SNPs. While just five markers showed a weak association with health condition, 32 SNPs were correlated with elevation, six of which were considered as high confidence associated SNPs, as indicated by at least two different association methods. Among these genes, thioredoxin and pseudo response regulator 1 (PRR1) are involved in redox homeostasis and ROS detoxification, APETALA2-like 3 (AP2L3), a transcription factor, is involved in seasonal apical growth, and a RPS2-like is a disease resistance gene. The function of some of these genes in controlling light-dependent reactions and metabolic processes suggests signatures of adaptation to local photoperiod and the synchronization of the circadian rhythm. This work provides new insights into the genetic basis of local adaptation over a shallow elevation gradient in Norway spruce.

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.

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.

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.

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.

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.

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.

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.