Association genetics of acetophenone defence against spruce budworm in mature white spruce

The Key Role of Plant Hormone Signaling Transduction and Flavonoid Biosynthesis Pathways in the Response of Chinese Pine (Pinus tabuliformis) to Feeding Stimulation by Pine Caterpillar (Dendrolimus tabulaeformis)

In nature, plants have developed a series of resistance mechanisms to face various external stresses. As understanding of the molecular mechanisms underlying plant resistance continues to deepen, exploring endogenous resistance in plants has become a hot topic in this field. Despite the multitude of studies on plant-induced resistance, how plants respond to stress under natural conditions remains relatively unclear. To address this gap, we investigated Chinese pine (Pinus tabuliformis) using pine caterpillar (Dendrolimus tabulaeformis) under natural conditions. Healthy Chinese pine trees, approximately 10 years old, were selected for studying induced resistance in Huangtuliangzi Forestry, Pingquan City, Chengde City, Hebei Province, China. Pine needles were collected at 2 h and 8 h after feeding stimulation (FS) via 10 pine caterpillars and leaf clipping control (LCC), to simulate mechanical damage caused by insect chewing for the quantification of plant hormones and transcriptome and metabolome assays. The results show that the different modes of treatments significantly influence the contents of JA and SA in time following treatment. Three types of differentially accumulated metabolites (DAMs) were found to be involved in the initial response, namely phenolic acids, lipids, and flavonoids. Weighted gene co-expression network analysis indicated that 722 differentially expressed genes (DEGs) are positively related to feeding stimulation and the specific enriched pathways are plant hormone signal transduction and flavonoid biosynthesis, among others. Two TIFY transcription factors (PtTIFY54 and PtTIFY22) and a MYB transcription factor (PtMYB26) were found to be involved in the interaction between plant hormones, mainly in the context of JA signal transduction and flavonoid biosynthesis. The results of this study provide an insight into how JA activates, serving as a reference for understanding the molecular mechanisms of resistance formation in conifers responding to mandibulate insects.

Photosynthesis, ionomics and metabolomics of the host-hemiparasite association Acacia gerrardii-Viscum schimperi

Viscum schimperi is an evergreen hemiparasitic plant that can grow on stems and branches of several tree species. It penetrates the host tissues and forms a vascular bridge (haustorium) to withdraw the nutritive resources. Its relationships with hosts remain unknown. This study aimed to investigate the physiological and biochemical attributes of the host-hemiparasite association Acacia gerrardii -Viscum schimperi . The hemiparasite exhibited 2.4- and 3.0-fold lower photosynthetic activity and water use efficiency, and 1.2- and 4.1-fold higher transpiration rate and stomatal conductance. Equally, it displayed 4.9- and 2.6-fold greater water potential and osmotic potential, and in least 3.0times more accumulated 39 K, 85 Rb and 51 V, compared to the host. Nevertheless, it had no detrimental effect on photosynthetic activity, water status and multi-element accumulations in the host. Based on metabolome profiling, V. schimperi could use xanthurenic acid and propylparaben to acquire potassium from the host, and N -1-naphthylacetamide and N -Boc-hydroxylamine to weaken or kill the distal part of the infected branch and to receive the total xylem contents. In contrast, A. gerrardii could used N -acetylserotonin, arecoline, acetophenone and 6-methoxymellein to defend against V. schimperi infection.

Genetic variation of Picea abies in response to the artificial inoculation of Heterobasidion parviporum

Norway spruce Picea abies is one of Europe's most economically important tree species. However, it is highly susceptible to the root rot fungus Heterobasidion parviporum. Climate change will benefit the pathogen as the tree host is weakened by, e.g., extended drought. Breeding can improve forest health, and several root rot genetic markers are suggested to improve the resistance of Norway spruce. This study aimed to compare the resistance (here defined as necrosis length) of Norway spruce families and genotypes against two strains of H. parviporum under different water availabilities. Our results show that the family and the genotype within the family have an impact on the necrosis length that is related to the aggressiveness of the fungal strains. Under low water conditions, the necrosis increased only in horizontal directions in phloem and sapwood. Similarly, the growth (seedling height) was not disturbed by abiotic stress (less water), indicating that the stress level (drought) was too low in this setting. The knowledge gained in this study could improve forest health in the changing climate by understanding the response of Norway spruce to pathogenic attacks under additional stress at the family level. This knowledge could be strategically used in forest breeding to improve the resistance of Norway spruce trees against root rot.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10342-023-01534-3.

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.

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.

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.

Analysis of the transcriptome of the needles and bark of Pinus radiata induced by bark stripping and methyl jasmonate

BACKGROUND

Plants are attacked by diverse insect and mammalian herbivores and respond with different physical and chemical defences. Transcriptional changes underlie these phenotypic changes. Simulated herbivory has been used to study the transcriptional and other early regulation events of these plant responses. In this study, constitutive and induced transcriptional responses to artificial bark stripping are compared in the needles and the bark of Pinus radiata to the responses from application of the plant stressor, methyl jasmonate. The time progression of the responses was assessed over a 4-week period.

RESULTS

Of the 6312 unique transcripts studied, 86.6% were differentially expressed between the needles and the bark prior to treatment. The most abundant constitutive transcripts were related to defence and photosynthesis and their expression did not differ between the needles and the bark. While no differential expression of transcripts were detected in the needles following bark stripping, in the bark this treatment caused an up-regulation and down-regulation of genes associated with primary and secondary metabolism. Methyl jasmonate treatment caused differential expression of transcripts in both the bark and the needles, with individual genes related to primary metabolism more responsive than those associated with secondary metabolism. The up-regulation of genes related to sugar break-down and the repression of genes related with photosynthesis, following both treatments was consistent with the strong down-regulation of sugars that has been observed in the same population. Relative to the control, the treatments caused a differential expression of genes involved in signalling, photosynthesis, carbohydrate and lipid metabolism as well as defence and water stress. However, non-overlapping transcripts were detected between the needles and the bark, between treatments and at different times of assessment. Methyl jasmonate induced more transcriptional responses in the bark than bark stripping, although the peak of expression following both treatments was detected 7 days post treatment application. The effects of bark stripping were localised, and no systemic changes were detected in the needles.

CONCLUSION

There are constitutive and induced differences in the needle and bark transcriptome of Pinus radiata. Some expression responses to bark stripping may differ from other biotic and abiotic stresses, which contributes to the understanding of plant molecular responses to diverse stresses. Whether the gene expression changes are heritable and how they differ between resistant and susceptible families identified in earlier studies needs further investigation.

Additive genetic variation in Pinus radiata bark chemistry and the chemical traits associated with variation in mammalian bark stripping

Secondary metabolites are suggested as a major mechanism explaining genetic variation in herbivory levels in Pinus radiata. The potential to incorporate these chemical traits into breeding/deployment programmes partly depends on the presence of additive genetic variation for the relevant chemical traits. In this study, near-infrared spectroscopy was used to quantify the constitutive and induced levels of 54 compounds in the bark of trees from 74 P. radiata full-sib families. The trees sampled for chemistry were protected from browsing and induced levels were obtained by subjecting half of the trees to artificial bark stripping. The treatment effect on bark chemistry was assessed along with narrow-sense heritability, the significance of non-additive genetic effects and the additive genetic correlations of compounds with bark stripping by mammalian herbivores that was observed in unprotected replicates of the field trial. The results indicated: (i) significant additive genetic variation, with low-moderate narrow-sense heritability estimates for most compounds; (ii) while significant induced effects were detected for some chemicals, no significant genetic variation in inducibility was detected; and (iii) sugars, fatty acids and a diterpenoid positively genetically correlated while a sesquiterpenoid negatively genetically correlated with bark stripping by the mammalian herbivore, the Bennett's wallaby (Macropus rufogriseus). At the onset of browsing, a trade-off with height was detected for selecting higher amounts of this sesquiterpenoid. However, overall, results showed potential to incorporate chemical traits into breeding/deployment programmes. The quantitative genetic analyses of the near infrared predicted chemical traits produced associations with mammalian bark stripping that mostly conform with those obtained using standard wet chemistry.

Evolution and Adaptation of Forest and Crop Pathogens in the Anthropocene

Anthropocene marks the era when human activity is making a significant impact on earth, its ecological and biogeographical systems. The domestication and intensification of agricultural and forest production systems have had a large impact on plant and tree health. Some pathogens benefitted from these human activities and have evolved and adapted in response to the expansion of crop and forest systems, resulting in global outbreaks. Global pathogen genomics data including population genomics and high-quality reference assemblies are crucial for understanding the evolution and adaptation of pathogens. Crops and forest trees have remarkably different characteristics, such as reproductive time and the level of domestication. They also have different production systems for disease management with more intensive management in crops than forest trees. By comparing and contrasting results from pathogen population genomic studies done on widely different agricultural and forest production systems, we can improve our understanding of pathogen evolution and adaptation to different selection pressures. We find that in spite of these differences, similar processes such as hybridization, host jumps, selection, specialization, and clonal expansion are shaping the pathogen populations in both crops and forest trees. We propose some solutions to reduce these impacts and lower the probability of global pathogen outbreaks so that we can envision better management strategies to sustain global food production as well as ecosystem services.

Genomic selection for resistance to spruce budworm in white spruce and relationships with growth and wood quality traits

With climate change, the pressure on tree breeding to provide varieties with improved resilience to biotic and abiotic stress is increasing. As such, pest resistance is of high priority but has been neglected in most tree breeding programs, given the complexity of phenotyping for these traits and delays to assess mature trees. In addition, the existing genetic variation of resistance and its relationship with productivity should be better understood for their consideration in multitrait breeding. In this study, we evaluated the prospects for genetic improvement of the levels of acetophenone aglycones (AAs) in white spruce needles, which have been shown to be tightly linked to resistance to spruce budworm. Furthermore, we estimated the accuracy of genomic selection (GS) for these traits, allowing selection at a very early stage to accelerate breeding. A total of 1,516 progeny trees established on five sites and belonging to 136 full-sib families from a mature breeding population in New Brunswick were measured for height growth and genotyped for 4,148 high-quality SNPs belonging to as many genes along the white spruce genome. In addition, 598 trees were assessed for levels of AAs piceol and pungenol in needles, and 578 for wood stiffness. GS models were developed with the phenotyped trees and then applied to predict the trait values of unphenotyped trees. AAs were under moderate-to-high genetic control (h 2: 0.43-0.57) with null or marginally negative genetic correlations with other traits. The prediction accuracy of GS models (GBLUP) for AAs was high (PAAC: 0.63-0.67) and comparable or slightly higher than pedigree-based (ABLUP) or BayesCπ models. We show that AA traits can be improved and that GS speeds up the selection of improved trees for insect resistance and for growth and wood quality traits. Various selection strategies were tested to optimize multitrait gains.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

Continent-wide population genomic structure and phylogeography of North America's most destructive conifer defoliator, the spruce budworm (Choristoneura fumiferana)

The spruce budworm, Choristoneura fumiferana, is presumed to be panmictic across vast regions of North America. We examined the extent of panmixia by genotyping 3,650 single nucleotide polymorphism (SNP) loci in 1975 individuals from 128 collections across the continent. We found three spatially structured subpopulations: Western (Alaska, Yukon), Central (southeastern Yukon to the Manitoba-Ontario border), and Eastern (Manitoba-Ontario border to the Atlantic). Additionally, the most diagnostic genetic differentiation between the Central and Eastern subpopulations was chromosomally restricted to a single block of SNPs that may constitute an island of differentiation within the species. Geographic differentiation in the spruce budworm parallels that of its principal larval host, white spruce (Picea glauca), providing evidence that spruce budworm and spruce trees survived in the Beringian refugium through the Last Glacial Maximum and that at least two isolated spruce budworm populations diverged with spruce/fir south of the ice sheets. Gene flow in the spruce budworm may also be affected by mountains in western North America, habitat isolation in West Virginia, regional adaptations, factors related to dispersal, and proximity of other species in the spruce budworm species complex. The central and eastern geographic regions contain individuals that assign to Eastern and Central subpopulations, respectively, indicating that these barriers are not complete. Our discovery of previously undetected geographic and genomic structure in the spruce budworm suggests that further population modelling of this ecologically important insect should consider regional differentiation, potentially co-adapted blocks of genes, and gene flow between subpopulations.

Hydroxyacetophenone defenses in white spruce against spruce budworm

We review a recently discovered white spruce (Picea glauca) chemical defense against spruce budworm (Choristoneura fumiferana) involving hydroxyacetophenones. These defense metabolites detected in the foliage accumulate variably as the aglycons, piceol and pungenol, or the corresponding glucosides, picein and pungenin. We summarize current knowledge of the genetic, genomic, molecular, and biochemical underpinnings of this defense and its effects on C. fumiferana. We present an update with new results on the ontogenic variation and the phenological window of this defense, including analysis of transcript responses in P. glauca to C. fumiferana herbivory. We also discuss this chemical defense from an evolutionary and a breeding context.