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.

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.

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.

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.

Drought increases Norway spruce susceptibility to the Eurasian spruce bark beetle and its associated fungi.

Drought affects the complex interactions between Norway spruce, the bark beetle Ips typographus and associated microorganisms. We investigated the interplay of tree water status, defense and carbohydrate reserves with the incidence of bark beetle attack and infection of associated fungi in mature spruce trees. We installed roofs to induce a 2-yr moderate drought in a managed spruce stand to examine a maximum of 10 roof and 10 control trees for resin flow (RF), predawn twig water potentials, terpene, phenolic and carbohydrate bark concentrations, and bark beetle borings in field bioassays before and after inoculation with Endoconidiophora polonica and Grosmannia penicillata. Drought-stressed trees showed more attacks and significantly longer fungal lesions than controls, but maintained terpene resin defenses at predrought levels. Reduced RF and lower mono- and diterpene, but not phenolic concentrations were linked with increased host selection. Bark beetle attack and fungi stimulated chemical defenses, yet G. penicillata reduced phenolic and carbohydrate contents. Chemical defenses did not decrease under mild, prolonged drought in our simulated small-scale biotic infestations. However, during natural mass attacks, reductions in carbon fixation under drought, in combination with fungal consumption of carbohydrates, may deplete tree defenses and facilitate colonization by I. typographus.

The effects of solar radiation on daily and seasonal stem increment of canopy trees in European temperate old-growth forests.

It is well established that solar irradiance greatly influences tree metabolism and growth through photosynthesis, but its effects acting through individual climate metrics have not yet been well quantified. Understanding these effects is crucial for assessing the impacts of climate change on forest ecosystems. To describe the effects of solar irradiance on tree growth, we installed 110 automatic dendrometers in two old-growth mountain forest reserves in Central Europe, performed detailed terrestrial and aerial laser scanning to obtain precise tree profiles, and used these to simulate the sum of solar irradiance received by each tree on a daily basis. Generalized linear mixed-effect models were applied to simulate the probability of growth and the growth intensity over seven growing seasons. Our results demonstrated various contrasting effects of solar irradiance on the growth of canopy trees. On the one hand, the highest daily growth rates corresponded with the highest solar irradiance potentials (i.e. the longest photoperiod). Intense solar irradiance significantly decreased tree growth, through an increase in the vapor pressure deficit. These effects were consistent for all species but had different magnitude. Tree growth is the most effective on long rainy/cloudy days with low solar irradiance.

Accounting for photosystem I photoinhibition sheds new light on seasonal acclimation strategies of boreal conifers.

The photosynthetic acclimation of boreal evergreen conifers is controlled by regulatory and photoprotective mechanisms that allow conifers to cope with extreme environmental changes. However, the underlying dynamics of photosystem II (PSII) and photosystem I (PSI) remain unresolved. Here, we investigated the dynamics of PSII and PSI during the spring recovery of photosynthesis in Pinus sylvestris and Picea abies using a combination of chlorophyll a fluorescence, P700 difference absorbance measurements, and quantification of key thylakoid protein abundances. In particular, we derived a new set of PSI quantum yield equations, correcting for the effects of PSI photoinhibition. Using the corrected equations, we found that the seasonal dynamics of PSII and PSI photochemical yields remained largely in balance, despite substantial seasonal changes in the stoichiometry of PSII and PSI core complexes driven by PSI photoinhibition. Similarly, the previously reported seasonal up-regulation of cyclic electron flow was no longer evident, after accounting for PSI photoinhibition. Overall, our results emphasize the importance of considering the dynamics of PSII and PSI to elucidate the seasonal acclimation of photosynthesis in overwintering evergreens. Beyond the scope of conifers, our corrected PSI quantum yields expand the toolkit for future studies aimed at elucidating the dynamic regulation of PSI.

Major tree species of Central European forests differ in their proportion of positive, negative, and nonstationary growth trends.

Temperate forests are undergoing significant transformations due to the influence of climate change, including varying responses of different tree species to increasing temperature and drought severity. To comprehensively understand the full range of growth responses, representative datasets spanning extensive site and climatic gradients are essential. This study utilizes tree-ring data from 550 sites from the temperate forests of Czechia to assess growth trends of six dominant Central European tree species (European beech, Norway spruce, Scots pine, silver fir, sessile and pedunculate oak) over 1990-2014. By modeling mean growth series for each species and site, and employing principal component analysis, we identified the predominant growth trends. Over the study period, linear growth trends were evident across most sites (56% increasing, 32% decreasing, and 10% neutral). The proportion of sites with stationary positive trends increased from low toward high elevations, whereas the opposite was true for the stationary negative trends. Notably, within the middle range of their distribution (between 500 and 700 m a.s.l.), Norway spruce and European beech exhibited a mix of positive and negative growth trends. While Scots pine growth trends showed no clear elevation-based pattern, silver fir and oaks displayed consistent positive growth trends regardless of site elevation, indicating resilience to the ongoing warming. We demonstrate divergent growth trajectories across space and among species. These findings are particularly important as recent warming has triggered a gradual shift in the elevation range of optimal growth conditions for most tree species and has also led to a decoupling of growth trends between lowlands and mountain areas. As a result, further future shifts in the elevation range and changes in species diversity of European temperate forests can be expected.

Temperature dependence of pollen germination and tube growth in conifers relates to their distribution along an elevational gradient in Washington State, USA.

BACKGROUND AND AIMS: Pollen germination and tube growth are essential processes for successful fertilization. They are among the most temperature-vulnerable stages and subsequently affect seed production and determine population persistence and species distribution under climate change. Our study aims to investigate intra- and inter-specific variations in the temperature dependence of pollen germination and tube length growth and to explore how these variations differ for pollen from elevational gradients. METHODS: We focused on three conifer species, Pinus contorta, Picea engelmannii, and Pinus ponderosa, with pollen collected from 350 to 2200m elevation in Washington State, USA. We conducted pollen viability tests at temperatures from 5 to 40°C in 5°C intervals. After testing for four days, we took images of these samples under a microscope to monitor pollen germination percentage (GP) and tube length (TL). We applied the Gamma function to describe the temperature dependence of GP and TL and estimated key parameters, including the optimal temperature for GP (Topt_GP) and TL (Topt_TL). KEY RESULTS: Results showed that pollen from three species and different elevations within a species have different GP, TL, Topt_GP, and Topt_TL. The population with a higher Topt_GP would also have a higher Topt_TL, while Topt_TL was generally higher than Topt_GP, i.e., a positive but not one-to-one relationship. However, only Pinus contorta showed that populations from higher elevations have lower Topt_GP and Topt_TL and vice versa. The variability in GP increased at extreme temperatures, whereas the variability in TL was greatest near Topt_TL. CONCLUSIONS: Our study demonstrates the temperature dependences of three conifers across a wide range of temperatures. Pollen germination and tube growth are highly sensitive to temperature conditions and vary among species and elevations, affecting their reproduction success during warming. Our findings can provide valuable insights to advance our understanding of how conifer pollen responds to rising temperatures.

Response of tree growth to drought variability in arid areas: Local hydroclimate and large-scale precipitation.

The impact of drought on terrestrial ecosystems is increasing, and the spatiotemporal heterogeneity of drought changes exacerbates the difficulty of determining ecosystem responses, especially in arid regions far from oceans. Tree rings have been widely used to understand how forest ecosystems respond to drought. However, the link between local hydroclimate variations related to tree rings and large-scale climate changes is not clear in the Qilian Mountains. Here, we used the tree ring width index to analyze the trend of Picea crassifolia growth and its relationship with climate in the middle Qilian Mountains. The results showed that the radial growth trend of Picea crassifolia is synchronized in the middle Qilian Mountains by calculating the Gleichläufigkeit index (GLK). Our analyses indicated that tree radial growth is positively correlated with drought during the growing season. Tree growth responds stably to drought (scPDSI and SPEI) and precipitation but unstably to temperature during 1950-2019. We further traced the meteorological factors that cause regional drought changes associated with radial growth. An increased total precipitation and decreased evaporation contribute to drought alleviation, favoring an increased tree radial growth. The increased total precipitation is mainly due to increased large-scale precipitation, which is related to water vapor transport changes. This study attempts to explore the influence of large-scale meteorology on regional drought change and its related tree radial growth response, which helps us to better understand the changes in forest ecosystems under climate change.

Fungicide-Resistant Botrytis in Forest Nurseries May Impact Disease Control in Norway Spruce.

Gray mold, caused by Botrytis spp., is a serious problem in Norway spruce seedling production in forest nurseries. From 2013 to 2019, 125 isolates of Botrytis were obtained from eight forest nurseries in Norway: 53 from Norway spruce seedlings, 16 from indoor air, 52 from indoor surfaces, and four from weeds growing close to seedlings. The majority of isolates were identified as B. cinerea, and over 60% of these were characterized as Botrytis group S. B. pseudocinerea isolates were obtained along with isolates with DNA sequence similarities to B. prunorum. Fungicide resistance was assessed with a mycelial growth assay, and resistance was found for the following: boscalid (8.8%), fenhexamid (33.6%), fludioxonil (17.6%), pyraclostrobin (36.0%), pyrimethanil (13.6%), and thiophanate-methyl (50.4%). Many isolates (38.4%) were resistant to two to six different fungicides. A selection of isolates was analyzed for the presence of known resistance-conferring mutations in the cytb, erg27, mrr1, sdhB, and tubA genes, and mutations leading to G143A, F412S, ΔL497, H272R, and E198A/F200Y were detected, respectively. Detection of fungicide resistance in Botrytis from Norway spruce and forest nursery facilities reinforces the necessity of employing resistance management strategies to improve control and delay development of fungicide resistance in the gray mold pathogens.[Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A pilot study of antimicrobial effects and ototoxicity of a Norway spruce (Picea abies) resin-based canine otic rinse product.

BACKGROUND: Norway spruce (Picea abies) resin-based products are used in human medicine. A resin-based otic rinse also could be useful in supportive care of canine otitis externa (COE), yet information on its antimicrobial effect against canine pathogens or ototoxicity is lacking. OBJECTIVES: To investigate the antimicrobial properties and ototoxicity of a commercial resin-based otic product. MATERIALS AND METHODS: Antimicrobial effect was evaluated using a standardised challenge test on Staphylococcus pseudintermedius, Corynebacterium auriscanis, Pseudomonas aeruginosa, Escherichia coli, Malassezia pachydermatis, and Streptococcus halichoeri strains to measure reduction in growth after 24 h exposure to the product. Effect on cell morphology was investigated by exposing S. pseudintermedius, C. auriscanis, P. aeruginosa and M. pachydermatis to the product in 20% and 100% (v/v) concentrations for 6, 24 and 48 h, and evaluating cells by transmission (TEM) and scanning (SEM) electron microscopy. An in vitro microbial kill-rate assay also was performed. Auditory brain stem response test, clinical evaluation and postmortem histological evaluation of ear canals were undertaken on experimental guinea pigs treated with the test product or saline controls. RESULTS: The product showed >log 5 growth reduction for all strains in the challenge test. TEM and SEM images showed clear changes in the cells' inner structures and deterioration of cells, and 100% (v/v) test product exposure induced microbial killing in 1-2 h. Ototoxicity was not detected in guinea pigs. CONCLUSIONS AND CLINICAL RELEVANCE: The product may be an option in supportive care of COE because of antimicrobial effects and lack of ototoxic properties in a guinea pig model.

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.

Phytochemical Characterization, Antioxidant, and Antimicrobial Activity of the Vegetative Buds from Romanian Spruce, Picea abies (L.) H. Karst.

This study aims to investigate the vegetative buds from Picea abies (spruce), naturally found in a central region of Romania, through a comprehensive analysis of the chemical composition to identify bioactive compounds responsible for pharmacological properties. Using HPLC/derivatization technique of GC-MS and quantitative spectrophotometric assays, the phenolic profile, and main components of an ethanolic extract from the buds were investigated. The essential oil was characterized by GC-MS. Moreover, the antioxidant activity with the DPPH method, and the antimicrobial activity were tested. Heavy metal detection was performed by graphite furnace atomic absorption spectrometry. The main components of the alcoholic extract were astragalin, quercetin, kaempferol, shikimic acid, and quinic acid. A total content of 25.32 ± 2.65 mg gallic acid equivalent per gram of dry plant (mg GAE/g DW) and of 10.54 ± 0.083 mg rutin equivalents/g of dry plant (mg RE/g DW) were found. The essential oil had D-limonene, α-cadinol, δ-cadinene, 13-epimanool, and δ-3-carene as predominant components. The spruce vegetative buds exhibited significant antioxidant activity (IC50 of 53 µg/mL) and antimicrobial effects against Staphylococcus aureus. Furthermore, concentrations of heavy metals Pb and Cd were below detection limits, suggesting that the material was free from potentially harmful contaminants. The results confirmed the potential of this indigenous species to be used as a source of compounds with pharmacological utilities.

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.

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

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

Transcriptome Analysis of Embryogenic and Non-Embryogenic Callus of Picea Mongolica.

Picea mongolica is a rare tree species in China, which is of great significance in combating desertification and improving the harsh ecological environment. Due to the low rate of natural regeneration, high mortality, and susceptibility to pests and cold springs, Picea mongolica has gradually become extinct. At present, somatic embryogenesis (SE) is the most effective method of micro-proliferation in conifers, but the induction rate of embryogenic callus (EC) is low, and EC is difficult to differentiate from non-embryonic callus (NEC). Therefore, the EC and NEC of Picea mongolica were compared from the morphology, histological, physiological, and transcriptional levels, respectively. Morphological observation showed that the EC was white and transparent filamentous, while the NEC was compact and brownish-brown lumpy. Histological analyses showed that the NEC cells were large and loosely arranged; the nuclei attached to the edge of the cells were small; the cytoplasm was low; and the cell gap was large and irregular. In the EC, small cells, closely arranged cells, and a large nucleus and nucleolus were observed. Physiological studies showed significant differences in ROS-scavenging enzymes between the EC and NEC. Transcriptome profiling revealed that 13,267 differentially expressed genes (DEGs) were identified, 3682 were up-regulated, and 9585 were down-regulated. In total, 63 GO terms had significant enrichment, 32 DEGs in plant hormone signal transduction pathway were identified, and 502 different transcription factors (TFs) were characterized into 38 TF families. Meanwhile, we identified significant gene expression trends associated with somatic embryo development in plant hormones (AUX/IAA, YUCCA, LEA, etc.), stress (GST, HSP, GLP, etc.), phenylpropanoid metabolism (4CL, HCT, PAL, etc.), and transcription factors (AP2/ERF, MYB, WOX, etc.). In addition, nine genes were chosen for RT-qPCR, and the results were consistent with RNA-Seq data. This study revealed the changes in morphology, histology, physiology, and gene expression in the differentiation of NEC into EC and laid the foundation for finding the key genes to promote EC formation.

Beyond the surface: exploring the mycobiome of Norway spruce under drought stress and with Heterobasidion parviporum.

The mycobiome, comprising fungi inhabiting plants, potentially plays a crucial role in tree health and survival amidst environmental stressors like climate change and pathogenic fungi. Understanding the intricate relationships between trees and their microbial communities is essential for developing effective strategies to bolster the resilience and well-being of forest ecosystems as we adopt more sustainable forest management practices. The mycobiome can be considered an integral aspect of a tree's biology, closely linked to its genotype. To explore the influence of host genetics and environmental factors on fungal composition, we examined the mycobiome associated with phloem and roots of Norway spruce (Picea abies (L.) Karst.) cuttings under varying watering conditions. To test the "mycobiome-associated-fitness" hypothesis, we compared seedlings artificially inoculated with Heterobasidion parviporum and control plants to evaluate mycobiome interaction on necrosis development. We aimed to 1) identify specific mycobiome species for the Norway spruce genotypes/families within the phloem and root tissues and their interactions with H. parviporum and 2) assess stability in the mycobiome species composition under abiotic disturbances (reduced water availability). The mycobiome was analyzed by sequencing the ribosomal ITS2 region. Our results revealed significant variations in the diversity and prevalence of the phloem mycobiome among different Norway spruce genotypes, highlighting the considerable impact of genetic variation on the composition and diversity of the phloem mycobiome. Additionally, specific mycobiome genera in the phloem showed variations in response to water availability, indicating the influence of environmental conditions on the relative proportion of certain fungal genera in Norway spruce trees. In the root mycobiome, key fungi such as Phialocephala fortinii and Paraphaeosphaeria neglecta were identified as conferring inhibitory effects against H. parviporum growth in Norway spruce genotypes. Furthermore, certain endophytes demonstrated greater stability in root ecosystems under low water conditions than ectomycorrhizal fungi. This knowledge can contribute to developing sustainable forest management practices that enhance the well-being of trees and their ecosystems, ultimately bolstering forest resilience.

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.

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.