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.

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.

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.

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.

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.

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 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.

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.

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.

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.

The dark septate endophyte Phialocephala sphaeroides suppresses conifer pathogen transcripts and promotes root growth of Norway spruce

Plant-associated microbes including dark septate endophytes (DSEs) of forest trees play diverse functional roles in host fitness including growth promotion and increased defence. However, little is known about the impact on the fungal transcriptome and metabolites during tripartite interaction involving plant host, endophyte and pathogen. To understand the transcriptional regulation of endophyte and pathogen during co-infection, Norway spruce (Picea abies) seedlings were infected with DSE Phialocephala sphaeroides, or conifer root-rot pathogen Heterobasidion parviporum, or both. Phialocephala sphaeroides showed low but stable transcripts abundance (a decrease of 40%) during interaction with Norway spruce and conifer pathogen. By contrast, H. parviporum transcripts were significantly reduced (92%) during co-infection. With RNA sequencing analysis, P. sphaeroides experienced a shift from cell growth to anti-stress and antagonistic responses, while it repressed the ability of H. parviporum to access carbohydrate nutrients by suppressing its carbohydrate/polysaccharide-degrading enzyme machinery. The pathogen on the other hand secreted cysteine peptidase to restrict free growth of P. sphaeroides. The expression of both DSE P. sphaeroides and pathogen H. parviporum genes encoding plant growth promotion products were equally detected in both dual and tripartite interaction systems. This was further supported by the presence of tryptophan-dependent indolic compound in liquid culture of P. sphaeroides. Norway spruce and Arabidopsis seedlings treated with P. sphaeroides culture filtrate exhibited auxin-like phenotypes, such as enhanced root hairs, and primary root elongation at low concentration but shortened primary root at high concentration. The results suggested that the presence of the endophyte had strong repressive or suppressive effect on H. parviporum transcripts encoding genes involved in nutrient acquisition.

Reducing Pre- and Post-Treatments in Cryopreservation Protocol and Testing Storage at -80 °C for Norway Spruce Embryogenic Cultures

Somatic embryogenesis (SE) is considered the most effective method for vegetative propagation of Norway spruce (Picea abies L. Karst). For mass propagation, a storage method that is able to handle large quantities of embryogenic tissues (ETs) reliably and at a low cost is required. The aim of the present study was to compare freezing at -80 °C in a freezer to cryopreservation using liquid nitrogen (LN) as a method for storing Norway spruce ETs. The possibility of simplifying both the pre-treatment and thawing processes in cryopreservation was also studied. The addition of abscisic acid (ABA) to the pre-treatment media and using polyethylene glycol PEG4000 instead of PEG6000 in a cryoprotectant solution were tested. Both the pre-and post-treatments on semi-solid media could be simplified by reducing the number of media, without any loss of genotype or embryo production capacity of ETs. On the contrary, the storage of ETs in a freezer at -80 °C instead of using LN was not possible, and the addition of ABA to the pre-treatment media did not provide benefits but increased costs. The lower regeneration rate after using PEG4000 instead of PEG6000 in a cryoprotectant solution in cryovials was unexpected and unwanted. The simplified pre-and post-treatment protocol will remarkably reduce the workload and costs in the mass-cryopreservation of future forest regeneration materials and in thawing the samples for mass propagations, respectively.

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.

Cone-setting in spruce is regulated by conserved elements of the age-dependent flowering pathway

Reproductive phase change is well characterized in angiosperm model species, but less studied in gymnosperms. We utilize the early cone-setting acrocona mutant to study reproductive phase change in the conifer Picea abies (Norway spruce), a gymnosperm. The acrocona mutant frequently initiates cone-like structures, called transition shoots, in positions where wild-type P. abies always produces vegetative shoots. We collect acrocona and wild-type samples, and RNA-sequence their messenger RNA (mRNA) and microRNA (miRNA) fractions. We establish gene expression patterns and then use allele-specific transcript assembly to identify mutations in acrocona. We genotype a segregating population of inbred acrocona trees. A member of the SQUAMOSA BINDING PROTEIN-LIKE (SPL) gene family, PaSPL1, is active in reproductive meristems, whereas two putative negative regulators of PaSPL1, miRNA156 and the conifer specific miRNA529, are upregulated in vegetative and transition shoot meristems. We identify a mutation in a putative miRNA156/529 binding site of the acrocona PaSPL1 allele and show that the mutation renders the acrocona allele tolerant to these miRNAs. We show co-segregation between the early cone-setting phenotype and trees homozygous for the acrocona mutation. In conclusion, we demonstrate evolutionary conservation of the age-dependent flowering pathway and involvement of this pathway in regulating reproductive phase change in the conifer P. abies.

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.

Enhanced lignin synthesis and ecotypic variation in defense-related gene expression in response to shade in Norway spruce

During the growth season, northern forests in Sweden daily receive more hours of far-red (FR)-enriched light or twilight (shade) as compared to southern forests. Norway spruce (shade-tolerant) are adapted to latitudinal variation in twilight characterized by a northward increase in FR requirement to maintain growth. Shade is a stressful condition that affects plant growth and increases plant's susceptibility to pathogen attack. Lignin plays a central role in plant defense and its metabolism is regulated by light wavelength composition (light quality). In the current work, we studied regulation of lignin synthesis and defense-related genes (growth-defense trade-offs) in response to shade in Norway spruce. In most angiosperms, light promotes lignin synthesis, whereas shade decreases lignin production leading to weaker stem, which may make plants more disease susceptible. In contrast, enhanced lignin synthesis was detected in response to shade in Norway spruce. We detected a higher number of immunity/defense-related genes up-regulated in northern populations as compared to south ones in response to shade. Enhanced lignin synthesis coupled with higher defense-related gene expression can be interpreted as an adaptive strategy for better survival in northern populations. Findings will contribute to ensuring deployment of well-adapted genetic material and identifying tree families with enhanced disease resistance.

Integrated Lipidomic and Transcriptomic Analysis Reveals Phospholipid Changes in Somatic Embryos of Picea asperata in Response to Partial Desiccation

Partial desiccation treatment (PDT) is an effective technology for promoting the germination and conversion of conifer somatic embryos (SEs). PDT, as a drought stress, induces intensive physiological responses in phospholipid metabolism, which are not well understood in the conifer SEs. Here, we integrated lipidomics, transcriptomics and proteomics analyses to reveal the molecular basis of lipid remodeling under PDT in Picea asperata SEs. Among the 82 lipid molecular species determined by mass spectrometry, phosphatidic acid (PA) had a significant effect after PDT and was the most critical lipid in the response to PDT. The transcriptomics results showed that multiple transcripts in the glycerolipid and glycerophospholipid metabolism pathways were differentially expressed, and these included five PLDα1 transcripts that catalyze the conversion of phosphatidylcholine (PC) to PA. Furthermore, the enzyme activity of this phospholipase D (PLD) was significantly enhanced in response to PDT, and PDT also significantly increased the protein level of PLDα1 (MA_10436582g0020). In addition, PA is a key factor in gibberellin, abscisic acid and ethylene signal transduction. One GDI1, one DELLA, three ABI1s, two SnRK2s, one CTR and 12 ERFs showed significantly differential expression between SEs before and after PDT in this study. Our data suggest that the observed increases in the PA contents might result from the activation of PLDα by PDT. PA not only affects the physical and chemical properties of the cell membrane but also participates in plant hormone signal transduction. Our work provides novel insight into the molecular mechanism through which PDT promotes the germination of SEs of coniferous tree species and fills the gap in the understanding of the mechanism of somatic embryo lipid remodeling in response to PDT.

Genotypic variation and plasticity in climate-adaptive traits after range expansion and fragmentation of red spruce (Picea rubens Sarg.)

Shifting range limits are predicted for many species as the climate warms. However, the rapid pace of climate change will challenge the natural dispersal capacity of long-lived, sessile organisms such as forest trees. Adaptive responses of populations will, therefore, depend on levels of genetic variation and plasticity for climate-responsive traits, which likely vary across the range due to expansion history and current patterns of selection. Here, we study levels of genetic and plastic variation for phenology and growth traits in populations of red spruce (Picea rubens), from the range core to the highly fragmented trailing edge. We measured more than 5000 offspring sampled from three genetically distinct regions (core, margin and edge) grown in three common gardens replicated along a latitudinal gradient. Genetic variation in phenology and growth showed low to moderate heritability and differentiation among regions, suggesting some potential to respond to selection. Phenology traits were highly plastic, but this plasticity was generally neutral or maladaptive in the effect on growth, revealing a potential liability under warmer climates. These results suggest future climate adaptation will depend on the regional availability of genetic variation in red spruce and provide a resource for the design and management of assisted gene flow. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Molecular and Chemical Screening for Inherent Disease Resistance Factors of Norway Spruce (Picea abies) Clones Against Conifer Stem Rot Pathogen Heterobasidion parviporum

Root and stem rot of conifer trees caused by Heterobasidion annosum species complex leads to huge economic losses in Europe, yet not much is known about the molecular and chemical basis for host resistance. To identify inherent chemical or molecular markers in clones found to be either resistant or susceptible, we sampled needle tissues of all the clones before pathogen inoculation. We conducted a short-term resistance screening by using the pathogen H. parviporum to inoculate 70 Norway spruce clones. Based on lesion size, subsets of highly susceptible and resistant clones were further analyzed. Terpene detection and RNA sequencing were performed to explore inherent variations in genotypes differing in resistance to pathogenic challenge at chemical and transcriptional levels. A negative correlation emerged between resistance and growth. Terpene profiles of resistant clones showed higher content of monoterpenes and sesquiterpenes, with concomitant increased transcript abundance of genes involved in the terpenoid pathway. A set of upregulated genes relevant to flavonoid biosynthesis was observed in resistant genotypes, whereas higher transcripts of lignin biosynthetic genes were prevalent in susceptible clones. Genes involved in flavonoid and lignin biosynthesis as well as terpene content may have a role in facilitating resistance of Norway spruce against H. parviporum. Our results provide strong support on the feasibility of sampling needle tissues before pathogen inoculation, and the approach could be of value for large-scale screening of novel biomarkers for durable resistance. The additional insights could form a basis for further research on resistance screening in this pathosystem.

Effect of Drought and Methyl Jasmonate Treatment on Primary and Secondary Isoprenoid Metabolites Derived from the MEP Pathway in the White Spruce Picea glauca

White spruce (Picea glauca) emits monoterpenes that function as defensive signals and weapons after herbivore attack. We assessed the effects of drought and methyl jasmonate (MeJA) treatment, used as a proxy for herbivory, on monoterpenes and other isoprenoids in P. glauca. The emission of monoterpenes was significantly increased after MeJA treatment compared to the control, but drought suppressed the MeJA-induced increase. The composition of the emitted blend was altered strongly by stress, with drought increasing the proportion of oxygenated compounds and MeJA increasing the proportion of induced compounds such as linalool and (E)-β-ocimene. In contrast, no treatment had any significant effect on the levels of stored monoterpenes and diterpenes. Among other MEP pathway-derived isoprenoids, MeJA treatment decreased chlorophyll levels by 40%, but had no effect on carotenoids, while drought stress had no impact on either of these pigment classes. Of the three described spruce genes encoding 1-deoxy-D-xylulose-5-phosphate synthase (DXS) catalyzing the first step of the MEP pathway, the expression of only one, DXS2B, was affected by our treatments, being increased by MeJA and decreased by drought. These findings show the sensitivity of monoterpene emission to biotic and abiotic stress regimes, and the mediation of the response by DXS genes.

Integrating genomic information and productivity and climate-adaptability traits into a regional white spruce breeding program

Tree improvement programs often focus on improving productivity-related traits; however, under present climate change scenarios, climate change-related (adaptive) traits should also be incorporated into such programs. Therefore, quantifying the genetic variation and correlations among productivity and adaptability traits, and the importance of genotype by environment interactions, including defense compounds involved in biotic and abiotic resistance, is essential for selecting parents for the production of resilient and sustainable forests. Here, we estimated quantitative genetic parameters for 15 growth, wood quality, drought resilience, and monoterpene traits for Picea glauca (Moench) Voss (white spruce). We sampled 1,540 trees from three open-pollinated progeny trials, genotyped with 467,224 SNP markers using genotyping-by-sequencing (GBS). We used the pedigree and SNP information to calculate, respectively, the average numerator and genomic relationship matrices, and univariate and multivariate individual-tree models to obtain estimates of (co)variance components. With few site-specific exceptions, all traits examined were under genetic control. Overall, higher heritability estimates were derived from the genomic- than their counterpart pedigree-based relationship matrix. Selection for height, generally, improved diameter and water use efficiency, but decreased wood density, microfibril angle, and drought resistance. Genome-based correlations between traits reaffirmed the pedigree-based correlations for most trait pairs. High and positive genetic correlations between sites were observed (average 0.68), except for those pairs involving the highest elevation, warmer, and moister site, specifically for growth and microfibril angle. These results illustrate the advantage of using genomic information jointly with productivity and adaptability traits, and defense compounds to enhance tree breeding selection for changing climate.

Norway spruce deploys tissue-specific responses during acclimation to cold

Climate change in the conifer-dominated boreal forest is expected to lead to warmer but more dynamic winter air temperatures, reducing the depth and duration of snow cover and lowering winter soil temperatures. To gain insight into the mechanisms that have enabled conifers to dominate extreme cold environments, we performed genome-wide RNA-Seq analysis from needles and roots of non-dormant two-year Norway spruce (Picea abies (L.) H. Karst), and contrasted these response to herbaceous model Arabidopsis We show that the main transcriptional response of Norway spruce needles exposed to cold was delayed relative to Arabidopsis, and this delay was associated with slower development of freezing tolerance. Despite this difference in timing, Norway spruce principally utilizes early response transcription factors (TFs) belonging to the same gene families as Arabidopsis, indicating broad evolutionary conservation of cold response networks. In keeping with their different metabolic and developmental states, needles and root of Norway spruce showed contrasting results. Regulatory network analysis identified both conserved TFs with known roles in cold acclimation (e.g. homologs of ICE1, AKS3, and of the NAC and AP2/ERF superfamilies), but also a root-specific bHLH101 homolog, providing functional insights into cold stress response strategies in Norway spruce.

Cryo-Treatment Enhances the Embryogenicity of Mature Somatic Embryos via the lncRNA-miRNA-mRNA Network in White Spruce

In conifers, somatic embryogenesis is uniquely initiated from immature embryos in a narrow time window, which is considerably hindered by the difficulty to induce embryogenic tissue (ET) from other tissues, including mature somatic embryos. In this study, the embryogenic ability of newly induced ET and DNA methylation levels was detected, and whole-transcriptome sequencing analyses were carried out. The results showed that ultra-low temperature treatment significantly enhanced ET induction from mature somatic embryos, with the induction rate from 0.4% to 15.5%, but the underlying mechanisms remain unclear. The newly induced ET showed higher capability in generating mature embryos than the original ET. DNA methylation levels fluctuated during the ET induction process. Here, WGCNA analysis revealed that OPT4, TIP1-1, Chi I, GASA5, GST, LAX3, WRKY7, MYBS3, LRR-RLK, PBL7, and WIN1 genes are involved in stress response and auxin signal transduction. Through co-expression analysis, lncRNAs MSTRG.505746.1, MSTRG.1070680.1, and MSTRG.33602.1 might bind to pre-novel_miR_339 to promote the expression of WRKY7 genes for stress response; LAX3 could be protected by lncRNAs MSTRG.1070680.1 and MSTRG.33602.1 via serving as sponges for novel_miR_495 to initiate auxin signal transduction; lncRNAs MSTRG.505746.1, MSTRG.1070680.1, and MSTRG.33602.1 might serve as sponges for novel_miR_527 to enhance the expression of Chi I for early somatic embryo development. This study provides new insight into the area of stress-enhanced early somatic embryogenesis in conifers, which is also attributable to practical applications.

Population structure and the influence of microenvironment and genetic similarity on individual growth at Alaskan white spruce treelines

Knowledge on the adaptation of trees to rapid environmental changes is essential to preserve forests and their ecosystem services under climate change. Treeline populations are particularly suitable for studying adaptation processes in trees, as environmental stress together with reduced gene flow can enhance local adaptation. We investigated white spruce (Picea glauca) populations in Alaska on one moisture-limited and two cold-limited treeline sites with a paired plot design of one forest and one treeline population each, resulting in six plots. Additionally, one forest plot in the middle of the distribution range complements the study design. We combined spatial, climatic and dendrochronological data with neutral genetic marker of 2203 trees to investigate population genetic structure and drivers of tree growth. We used several individual-based approaches including random slope mixed-effects models to test the influence of genetic similarity and microenvironment on growth performance. A high degree of genetic diversity was found within each of the seven plots associated with high rates of gene flow. We discovered a low genetic differentiation between the three sites which was better explained by geographic distances than by environmental differences, indicating genetic drift as the main driver of population differentiation. Our findings indicated that microenvironmental features had an overall larger influence on growth performances than genetic similarity among individuals. The effects of climate on growth differed between sites but were smaller than the effect of tree size. Overall, our results suggest that the high genetic diversity of white spruce may result in a wider range of phenotypes which enhances the efficiency of selection when the species is facing rapid climatic changes. In addition, the large intra-individual variability in growth responses may indicate the high phenotypic plasticity of white spruce which can buffer short-term environmental changes and, thus, allow enduring the present changing climate conditions.

An atlas of the Norway spruce needle seasonal transcriptome

Boreal conifers possess a tremendous ability to survive and remain evergreen during harsh winter conditions and resume growth during summer. This is enabled by coordinated regulation of major cellular functions at the level of gene expression, metabolism, and physiology. Here we present a comprehensive characterization of the annual changes in the global transcriptome of Norway spruce (Picea abies) needles as a resource to understand needle development and acclimation processes throughout the year. In young, growing needles (May 15 until June 30), cell walls, organelles, etc., were formed, and this developmental program heavily influenced the transcriptome, explained by over-represented Gene Ontology (GO) categories. Later changes in gene expression were smaller but four phases were recognized: summer (July-August), autumn (September-October), winter (November-February), and spring (March-April), where over-represented GO categories demonstrated how the needles acclimated to the various seasons. Changes in the seasonal global transcriptome profile were accompanied by differential expression of members of the major transcription factor families. We present a tentative model of how cellular activities are regulated over the year in needles of Norway spruce, which demonstrates the value of mining this dataset, accessible in ConGenIE together with advanced visualization tools.

Candidate regulators and target genes of drought stress in needles and roots of Norway spruce

Drought stress impacts seedling establishment, survival and whole-plant productivity. Molecular responses to drought stress have been most extensively studied in herbaceous species, mostly considering only aboveground tissues. Coniferous tree species dominate boreal forests, which are predicted to be exposed to more frequent and acute drought as a result of ongoing climate change. The associated impact at all stages of the forest tree life cycle is expected to have large-scale ecological and economic impacts. However, the molecular response to drought has not been comprehensively profiled for coniferous species. We assayed the physiological and transcriptional response of Picea abies (L.) H. Karst seedling needles and roots after exposure to mild and severe drought. Shoots and needles showed an extensive reversible plasticity for physiological measures indicative of drought-response mechanisms, including changes in stomatal conductance (gs), shoot water potential and abscisic acid (ABA). In both tissues, the most commonly observed expression profiles in response to drought were highly correlated with the ABA levels. Still, root and needle transcriptional responses contrasted, with extensive root-specific down-regulation of growth. Comparison between previously characterized Arabidopsis thaliana L. drought-response genes and P. abies revealed both conservation and divergence of transcriptional response to drought. In P. abies, transcription factors belonging to the ABA responsive element(ABRE) binding/ABRE binding factors ABA-dependent pathway had a more limited role. These results highlight the importance of profiling both above- and belowground tissues, and provide a comprehensive framework to advance the understanding of the drought response of P. abies. The results demonstrate that a short-term, severe drought induces severe physiological responses coupled to extensive transcriptome modulation and highlight the susceptibility of Norway spruce seedlings to such drought events.

Leveraging breeding programs and genomic data in Norway spruce (Picea abies L. Karst) for GWAS analysis

BACKGROUND

Genome-wide association studies (GWAS) identify loci underlying the variation of complex traits. One of the main limitations of GWAS is the availability of reliable phenotypic data, particularly for long-lived tree species. Although an extensive amount of phenotypic data already exists in breeding programs, accounting for its high heterogeneity is a great challenge. We combine spatial and factor-analytics analyses to standardize the heterogeneous data from 120 field experiments of 483,424 progenies of Norway spruce to implement the largest reported GWAS for trees using 134 605 SNPs from exome sequencing of 5056 parental trees.

RESULTS

We identify 55 novel quantitative trait loci (QTLs) that are associated with phenotypic variation. The largest number of QTLs is associated with the budburst stage, followed by diameter at breast height, wood quality, and frost damage. Two QTLs with the largest effect have a pleiotropic effect for budburst stage, frost damage, and diameter and are associated with MAP3K genes. Genotype data called from exome capture, recently developed SNP array and gene expression data indirectly support this discovery.

CONCLUSION

Several important QTLs associated with growth and frost damage have been verified in several southern and northern progeny plantations, indicating that these loci can be used in QTL-assisted genomic selection. Our study also demonstrates that existing heterogeneous phenotypic data from breeding programs, collected over several decades, is an important source for GWAS and that such integration into GWAS should be a major area of inquiry in the future.

Physiochemical and molecular features of the necrotic lesion in the Heterobasidion-Norway spruce pathosystem

In the forest of Northern Hemisphere, the fungi Heterobasidion annosum (Fr.) Bref. s.l. causes severe root and stem rot diseases, dramatically reducing the wood quality of conifer trees. The hallmark of the host response during the infection process is the formation of necrotic lesions and reaction zones. To characterize physiochemical and molecular features of the necrotic lesion, we conducted artificial inoculations on Norway spruce plants at different developmental stages: seedlings, young and mature trees. The results were further compared against data available on the formation of reaction zones. Strong necrosis browning or enlarged necrotic lesions were observed in infected tissues. This was accompanied by elevated pH. However, the increased pH, around 6.0 in necrotic lesions, was not as high as that documented in reaction zones, above 7.0 as marked by the intensity of the blue colour in response to 2,6-dichlorophenol indophenol dye. Peroxidase activity increased in infected plants and RNA-seq analysis of necrotic lesions showed marked upregulation of defence-related genes. Our findings highlight similarities and differences between the reaction zone and necrotic lesion formation in response of conifer trees to biotic stress.

Quantitative analysis of differential dehydrin regulation in pine and spruce seedlings under water deficit

Dehydrins are well-known components of plant responses to different stresses that cause dehydration, including drought, freezing, salinity, etc. In conifers, the dehydrin gene family is very large, implying that the members of this family have important physiological functions in conifer stress tolerance. However, dehydrin gene expression displays a wide range of responses to stress, from thousand-fold increased expression to decreased expression, and it is generally unknown how regulatory systems are connected at the mRNA and protein levels. Therefore, we studied these aspects of dehydrin regulation in Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies (L.) H. Karst) seedlings under polyethylene glycol 6000-induced osmotic stress ranging from relatively low (culture medium water potential of -0.15 MPa) to very high (-1.0 MPa) intensities. In pine, the major dehydrin protein was Dhn1 in both the roots and needles, and in spruce, two isoforms of the Dhn4 protein were the major dehydrins; both of these proteins are AESK-type dehydrins. The genes encoding these major proteins were highly expressed even under control conditions; surprisingly, we also observed several highly expressed dehydrin genes that were not abundantly translated. Under osmotic stress, the most prominent expression changes were observed for the dehydrin genes with low basal expression levels, whereas highly expressed genes generally demonstrated rather modest changes in expression. We report proposed constitutive physiological functions of the AESK-type dehydrins in Pinaceae plants.

Integrated biological and chemical characterisation of a pair of leonardesque canal lock gates

The Museo Nazionale della Scienza e della Tecnologia “Leonardo da Vinci” in Milan is exposing two pairs of canal lock gates, used to control the water flow in Milan canal system, whose design appears in the Leonardo’s Codex Atlanticus. The wood present in the gates has been deeply characterised by mean of a multidisciplinary investigation involving i) DNA barcoding of wood fragments; ii) microbial community characterisation, and iii) chemical analyses. DNA barcoding revealed that two fragments of the gates belonged to wood species widely used in the middle age: Fagus sylvatica and Picea abies. The chemical characterisations were based on the use of ionic liquid as dissolving medium in order to analyse the entire cell wall material by means of Gel Permeation Chromatography (GPC) and 2D-NMR-HSQC techniques. This multidisciplinary analytical approach was able to highlight the complex nature of the degradation occurred during the gate operation (XVI-XVIII centuries): an intricate interplay between microbial populations (i.e. Shewanella), inorganic factors (i.e. iron from nails), physical factors and the lignocellulosic material.

Effect of elevated gibberellic acid application on growth and gene expression patterns in white spruce families from a tree improvement program in Alberta, Canada

Nine open-pollinated families of Picea glauca (Moench) Voss from the Region D1 Controlled Parentage Program (Alberta, Canada) were systematically chosen from fast, medium and slow-growth rankings based on breeding values for height from field progeny tests at age 30 years. Seeds from these families were sown and grown to age 3 years to analyze the performance and correlations of growth, physiological traits and expression of gibberellin-related genes, with and without elevated gibberellic acid 3 (GA3) application, under greenhouse conditions. We observed a significant interaction effect between families and growth groups subjected to 50 μg μl-1 of GA3 treatment, causing a decrease in apical internode length, diameter, volume and absolute transcript level for fast-growing families but an increase for families in the slow-growth group for the same traits. We also observed that in the apical internode, the gene PgGA20ox1 had significantly more relative expression under the elevated GA3 treatment than the control trees. In the stem, PgGA3ox1 showed a significantly higher relative expression under elevated GA3 treatment compared with control trees. Also, the slow-growth group showed more relative expression of PgGA20ox1 (in the apical internode) and PgGA3ox1 (in the stem) than the fast-growth group. The apical internode length and diameter significantly increased by 24% and 16%, respectively, with the hormone treatment in the slow growing group. In general, the PgGID1 and PgDELLA1 genes were upregulated and downregulated respectively, in spruce shoots under the GA3 treatment, meaning a positive feedback regulation by those genes were influencing PgGA20ox1 and PgGA3ox1 expression in that tissue type. Moreover, there was a significant correlation between absolute transcript levels of PgGA20ox1 in the apical internode and apical internode length, and absolute transcript levels of PgGA3ox1 in the stem and the diameter, in the fast-growth group families. This study shows that expression of GA genes is a limiting factor for growth in certain white spruce families with a complex feedback mechanism. Finally, absolute transcript levels of endogenous GA relative to growth parameters in juvenile seedlings could potentially be used to accelerate the early selection of families with inherently rapid apical and radial growth expansion.

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

Brassinosteroids synthesised by CYP85A/A1 but not CYP85A2 function via a BRI1-like receptor but not via BRI1 in Picea abies

Brassinosteroids (BRs) are essential plant hormones. In angiosperms, brassinolide and castasterone, the first and second most active BRs, respectively, are synthesised by CYP85A2 and CYP85A/A1, respectively. BRs in angiosperms function through an essential receptor, BR Insensitive 1 (BRI1). In addition, some angiosperms also have non-essential BRI1-like 1/3 (BRL1/3). In conifers, BRs promote seed germination under drought stress; however, how BRs function in gymnosperms is unknown. In this study, we performed functional complementation of BR biosynthesis and receptor genes from Picea abies with respective Arabidopsis mutants. We found that P. abies possessed functional PaCYP85A and PaBRL1 but not PaCYP85A2 or PaBRI1, and this results in weak BR signaling, and both PaCYP85A and PaBRL1 were abundantly expressed. However, neither BR treatment of P. abies seedlings nor expression of PaBRL1 in the Arabidopsis Atbri1 mutant promoted plant height, despite the fact that BR-responsive genes were activated. Importantly, chimeric AtBRI1 replaced with the BR-binding domain of PaBRL1 complemented the Atbri1 phenotypes. Furthermore, PaBRL1 had less kinase activity than BRI1 in vitro. Overall, P. abies had weak but still active BR signaling, explaining aspects of its slow growth and high stress tolerance. Our study sheds light on the functional and evolutionary significance of distinct BR signaling that is independent of BRI1 and brassinolide.

Molecular signatures of local adaptation to light in Norway spruce

MAIN CONCLUSION

Transcriptomic and exome capture analysis reveal an adaptive cline for shade tolerance in Norway spruce. Genes involved in the lignin pathway and immunity seem to play a potential role in contributing towards local adaptation to light. The study of natural variation is an efficient method to elucidate how plants adapt to local climatic conditions, a key process for the evolution of a species. Norway spruce is a shade-tolerant conifer in which the requirement of far-red light for growth increases latitudinally northwards. The objective of the study is to characterize the genetic control of local adaptation to light enriched in far-red in Norway spruce, motivated by a latitudinal gradient for the Red:Far-red (R:FR) ratio to which Norway spruce has been proven to be genetically adapted. We have established the genomic signatures of local adaptation by conducting transcriptomic (total RNA-sequencing) and genomic analyses (exome capture), for the identification of genes differentially regulated along the cline. RNA-sequencing revealed 274 differentially expressed genes in response to SHADE (low R:FR light), between the southern and northern natural populations in Sweden. Exome capture included analysis of a uniquely large data set (1654 trees) that revealed missense variations in coding regions of nine differentially expressed candidate genes, which followed a latitudinal cline in allele and genotype frequencies. These genes included five transcription factors involved in vital processes like bud-set/bud-flush, lignin pathway, and cold acclimation and other genes that take part in cell-wall remodeling, secondary cell-wall thickening, response to starvation, and immunity. Based on these results, we suggest that the northern populations might not only be able to adjust their growing season in response to low R:FR light, but they may also be better adapted towards disease resistance by up-regulation of the lignin pathway that is linked to immunity. This forms a concrete basis for local adaptation to light quality in Norway spruce, one of the most economically important conifer tree species in Sweden.

Genetic control of tracheid properties in Norway spruce wood

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

Polycomb Repressive Complex 2-mediated histone modification H3K27me3 is associated with embryogenic potential in Norway spruce

Epigenetic reprogramming during germ cell formation is essential to gain pluripotency and thus embryogenic potential. The histone modification H3K27me3, which is catalysed by the Polycomb repressive complex 2 (PRC2), regulates important developmental processes in both plants and animals, and defects in PRC2 components cause pleiotropic developmental abnormalities. Nevertheless, the role of H3K27me3 in determining embryogenic potential in gymnosperms is still elusive. To address this, we generated H3K27me3 profiles of Norway spruce (Picea abies) embryonic callus and non-embryogenic callus using CUT&RUN, which is a powerful method for chromatin profiling. Here, we show that H3K27me3 mainly accumulated in genic regions in the Norway spruce genome, similarly to what is observed in other plant species. Interestingly, H3K27me3 levels in embryonic callus were much lower than those in the other examined tissues, but markedly increased upon embryo induction. These results show that H3K27me3 levels are associated with the embryogenic potential of a given tissue, and that the early phase of somatic embryogenesis is accompanied by changes in H3K27me3 levels. Thus, our study provides novel insights into the role of this epigenetic mark in spruce embryogenesis and reinforces the importance of PRC2 as a key regulator of cell fate determination across different plant species.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Individual reproductive success in Norway spruce natural populations depends on growth rate, age and sensitivity to temperature

Quantifying the individual reproductive success and understanding its determinants is a central issue in evolutionary research for the major consequences that the transmission of genetic variation from parents to offspring has on the adaptive potential of populations. Here, we propose to distil the myriad of information embedded in tree-ring time series into a set of tree-ring-based phenotypic traits to be investigated as potential drivers of reproductive success in forest trees. By using a cross-disciplinary approach that combines parentage analysis and a thorough dendrophenotypic characterisation of putative parents, we assessed sex-specific relationships between such dendrophenotypic traits (i.e., age, growth rate and parameters describing sensitivity to climate and to extreme climatic events) and reproductive success in Norway spruce. We applied a full probability method for reconstructing parent-offspring relationships between 604 seedlings and 518 adult trees sampled within five populations from southern and central Europe. We found that individual female and male reproductive success was positively associated with tree growth rate and age. Female reproductive success was also positively influenced by the correlation between growth and the mean temperature of the previous vegetative season. Overall, our results showed that Norway spruce individuals with the highest fitness are those who are able to keep high-growth rates despite potential growth limitations caused by reproductive costs and climatic limiting conditions. Identifying such functional links between the individual ecophysiological behaviour and its evolutionary gain would increase our understanding on how natural selection shapes the genetic composition of forest tree populations over time.

Profiles of endogenous phytohormones and expression of some hormone-related genes in Scots pine and Norway spruce seedlings under water deficit

Different plant hormones are involved in plant adaptation to water deficit. In comparison to angiosperms, little is known about the impact of drought on the pool of phytohormones in gymnosperms. Therefore, we studied the effect of polyethylene glycol-induced water deficit on the changes in content of different phytohormones in Scots pine and Norway spruce seedlings, which are known for their different strategies of adaptation to water deficit. The following hormone classes were analysed: cytokinins, auxins, jasmonates, salicylic and benzoic acids, and 1-aminocyclopropane-1-carboxylic acid (an ethylene precursor). No consistent reaction to water stress was observed for the content of well-known stress-related hormones - salicylic acid and jasmonates. In contrast, drought induced a dose-dependent accumulation of cytokinins in pine needles, with less profound changes in spruce needles. The most prominent changes were observed for 1-aminocyclopropane-1-carboxylic acid content, which increased several-fold in spruce roots and pine needles under water deficit. Water-deficit-induced changes in the contents of cytokinins and 1-aminocyclopropane-1-carboxylic acid were accompanied by the differential regulation of genes involved in the metabolism of these hormones. Possible links between changes in hormone pools and the adaptation of seedlings to water deficit are discussed.

RNA-Seq and secondary metabolite analyses reveal a putative defence-transcriptome in Norway spruce (Picea abies) against needle bladder rust (Chrysomyxa rhododendri) infection

BACKGROUND

Norway spruce trees in subalpine forests frequently face infections by the needle rust fungus Chrysomyxa rhododendri, which causes significant growth decline and increased mortality of young trees. Yet, it is unknown whether trees actively respond to fungal attack by activating molecular defence responses and/or respective gene expression.

RESULTS

Here, we report results from an infection experiment, in which the transcriptomes (via RNA-Seq analysis) and phenolic profiles (via UHPLC-MS) of control and infected trees were compared over a period of 39 days. Gene expression between infected and uninfected ramets significantly differed after 21 days of infection and revealed already known, but also novel candidate genes involved in spruce molecular defence against pathogens.

CONCLUSIONS

Combined RNA-Seq and biochemical data suggest that Norway spruce response to infection by C. rhododendri is restricted locally and primarily activated between 9 and 21 days after infestation, involving a potential isolation of the fungus by a hypersensitive response (HR) associated with an activation of phenolic pathways. Identified key regulatory genes represent a solid basis for further specific analyses in spruce varieties with varying susceptibility, to better characterise resistant clones and to elucidate the resistance mechanism.

Two dominant boreal conifers use contrasting mechanisms to reactivate photosynthesis in the spring

Boreal forests are dominated by evergreen conifers that show strongly regulated seasonal photosynthetic activity. Understanding the mechanisms behind seasonal modulation of photosynthesis is crucial for predicting how these forests will respond to changes in seasonal patterns and how this will affect their role in the terrestrial carbon cycle. We demonstrate that the two co-occurring dominant boreal conifers, Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies), use contrasting mechanisms to reactivate photosynthesis in the spring. Scots pine downregulates its capacity for CO2 assimilation during winter and activates alternative electron sinks through accumulation of PGR5 and PGRL1 during early spring until the capacity for CO2 assimilation is recovered. In contrast, Norway spruce lacks this ability to actively switch between different electron sinks over the year and as a consequence suffers severe photooxidative damage during the critical spring period.

Regulation of Somatic Embryo Development in Norway Spruce

Somatic embryogenesis in Norway spruce combined with reverse genetics can be used as a model to study the regulation of embryo development in conifers. The somatic embryo system includes a sequence of developmental stages, which are similar in morphology to their zygotic counterparts. The system can be sufficiently synchronized to enable the collection and study of a large number of somatic embryos at each developmental stage.Here we describe a protocol for establishing transgenic cell lines in which genes of interest are upregulated or downregulated. Furthermore, we present methods for comparing embryo morphology and development in transgenic and control cell lines, including phenotyping the embryos, histological analysis, and tracking embryo development. The expression pattern of different genes is determined by GUS reporter assays.

Genotypic variation in Norway spruce correlates to fungal communities in vegetative buds

The taxonomically diverse phyllosphere fungi inhabit leaves of plants. Thus, apart from the fungi's dispersal capacities and environmental factors, the assembly of the phyllosphere community associated with a given host plant depends on factors encoded by the host's genome. The host genetic factors and their influence on the assembly of phyllosphere communities under natural conditions are poorly understood, especially in trees. Recent work indicates that Norway spruce (Picea abies) vegetative buds harbour active fungal communities, but these are hitherto largely uncharacterized. This study combines internal transcribed spacer sequencing of the fungal communities associated with dormant vegetative buds with a genome-wide association study (GWAS) in 478 unrelated Norway spruce trees. The aim was to detect host loci associated with variation in the fungal communities across the population, and to identify loci correlating with the presence of specific, latent, pathogens. The fungal communities were dominated by known Norway spruce phyllosphere endophytes and pathogens. We identified six quantitative trait loci (QTLs) associated with the relative abundance of the dominating taxa (i.e., top 1% most abundant taxa). Three additional QTLs associated with colonization by the spruce needle cast pathogen Lirula macrospora or the cherry spruce rust (Thekopsora areolata) in asymptomatic tissues were detected. The identification of the nine QTLs shows that the genetic variation in Norway spruce influences the fungal community in dormant buds and that mechanisms underlying the assembly of the communities and the colonization of latent pathogens in trees may be uncovered by combining molecular identification of fungi with GWAS.

Ray Parenchymal Cells Contribute to Lignification of Tracheids in Developing Xylem of Norway Spruce

A comparative transcriptomic study and a single-cell metabolome analysis were combined to determine whether parenchymal ray cells contribute to the biosynthesis of monolignols in the lignifying xylem of Norway spruce (Picea abies). Ray parenchymal cells may function in the lignification of upright tracheids by supplying monolignols. To test this hypothesis, parenchymal ray cells and upright tracheids were dissected with laser-capture microdissection from tangential cryosections of developing xylem of spruce trees. The transcriptome analysis revealed that among the genes involved in processes typical for vascular tissues, genes encoding cell wall biogenesis-related enzymes were highly expressed in both developing tracheids and ray cells. Interestingly, most of the shikimate and monolignol biosynthesis pathway-related genes were equally expressed in both cell types. Nonetheless, 1,073 differentially expressed genes were detected between developing ray cells and tracheids, among which a set of genes expressed only in ray cells was identified. In situ single cell metabolomics of semi-intact plants by picoliter pressure probe-electrospray ionization-mass spectrometry detected monolignols and their glycoconjugates in both cell types, indicating that the biosynthetic route for monolignols is active in both upright tracheids and parenchymal ray cells. The data strongly support the hypothesis that in developing xylem, ray cells produce monolignols that contribute to lignification of tracheid cell walls.

Testing of the narrow crowned Norway spruce ideotype (Picea abies f. pendula) and the hybrids with normal crown form (pyramidalis) in multisite comparative trials

The study aims to analyse the stability of the narrow crowned Norway spruce (pendula form) compared to the normal spruce form (pyramidalis form) and the hybrids of the two forms, in 5 field trials (Comandău, Lepşa 1&2, Ilva Mică and Voineasa) located in the Romanian Carpathians. Trees height (Th), breast height diameter (Dbh), height growth of the last year, crown diameter (Cd), number of branches per whorl (Nbw) and dominant branch diameter (Dbd) traits were measured and survival rate (Sr) was determined, at 20 years old. Also, branches finesse (Bf), trees volume (Tv) and trees slenderness (Ts) were calculated. In order to compare the wood density (Cwd) there were collected cores. In all trials ANOVA revealed significant (p < 0.05) differences between the two forms of spruce and the hybrids (mainly between those that have a different crown form mother), especially for the stability and quality traits. Factorial ANOVA revealed a high influence (p < 0.001) of the locality and also a significant influence (p < 0.05) of the locality × spruce form interaction. The factor "form" was significant for some traits involved in Norway spruce stability (Ts, Cd, Nbw). The pendula trees present higher values for Sr, Dbh and Tv, and lower values for Ts, Cd, Nbw, Dbd and Bf, compared to pyramidalis spruce form, which showed a higher stability. Heritability was in generally low (h² < 0.4), with exceptions of Ts which presents a medium rate of heritage. For the same trait, different heritability was registered in different environmental conditions. The Cwd was higher only with 2% for the pendula form in Lepşa trial, while in Comandău trial the pyramidalis registered a higher value (7%). In the new breeding programme, the selection strategy may be pursued with the pendula trees selection based on Ts and branches traits.

Comparative sensitivity to gamma radiation at the organismal, cell and DNA level in young plants of Norway spruce, Scots pine and Arabidopsis thaliana

Persistent DNA damage in gamma-exposed Norway spruce, Scots pine and Arabidopsis thaliana, but persistent adverse effects at the organismal and cellular level in the conifers only. Gamma radiation emitted from natural and anthropogenic sources may have strong negative impact on plants, especially at high dose rates. Although previous studies implied different sensitivity among species, information from comparative studies under standardized conditions is scarce. In this study, sensitivity to gamma radiation was compared in young seedlings of the conifers Scots pine and Norway spruce and the herbaceous Arabidopsis thaliana by exposure to ⁶⁰Co gamma dose rates of 1-540 mGy h^-1 for 144 h, as well as 360 h for A. thaliana. Consistent with slightly less prominent shoot apical meristem, in the conifers growth was significantly inhibited with increasing dose rate ≥ 40 mGy h^-1. Post-irradiation, the conifers showed dose-rate-dependent inhibition of needle and root development consistent with increasingly disorganized apical meristems with increasing dose rate, visible damage and mortality after exposure to ≥ 40 mGy h^-1. Regardless of gamma duration, A. thaliana showed no visible or histological damage or mortality, only delayed lateral root development after ≥ 100 mGy h^-1 and slightly, but transiently delayed post-irradiation reproductive development after ≥ 400 mGy h^-1. In all species dose-rate-dependent DNA damage occurred following ≥ 1-10 mGy h^-1 and was still at a similar level at day 44 post-irradiation. In conclusion, the persistent DNA damage (possible genomic instability) following gamma exposure in all species may suggest that DNA repair is not necessarily mobilized more extensively in A. thaliana than in Norway spruce and Scots pine, and the far higher sensitivity at the organismal and cellular level in the conifers indicates lower tolerance to DNA damage than in A. thaliana.