Killing two enemies with one stone? Genomics of resistance to two sympatric pathogens in Norway spruce.

Trees must cope with the attack of multiple pathogens, often simultaneously during their long lifespan. Ironically, the genetic and molecular mechanisms controlling this process are poorly understood. The objective of this study was to compare the genetic component of resistance in Norway spruce to Heterobasidion annosum s.s. and its sympatric congener Heterobasidion parviporum. Heterobasidion root- and stem-rot is a major disease of Norway spruce caused by members of the Heterobasidion annosum species complex. Resistance to both pathogens was measured using artificial inoculations in half-sib families of Norway spruce trees originating from central to northern Europe. The genetic component of resistance was analysed using 63,760 genome-wide exome-capture sequenced SNPs and multitrait genome-wide associations. No correlation was found for resistance to the two pathogens; however, associations were found between genomic variants and resistance traits with synergic or antagonist pleiotropic effects to both pathogens. Additionally, a latitudinal cline in resistance in the bark to H. annosum s.s. was found; trees from southern latitudes, with a later bud-set and thicker stem diameter, allowed longer lesions, but this was not the case for H. parviporum. In summary, this study detects genomic variants with pleiotropic effects which explain multiple disease resistance from a genic level and could be useful for selection of resistant trees to both pathogens. Furthermore, it highlights the need for additional research to understand the evolution of resistance traits to multiple pathogens in trees.

Annulohypoxylon sp. strain MUS1, an endophytic fungus isolated from Taxus wallichiana Zucc., produces taxol and other bioactive metabolites.

The current study focuses on the isolation and in vitro characterization of bioactive metabolites produced by endophytic fungi isolated from the Himalayan yew (Taxus wallichiana Zucc.). The endophytic fungi were isolated on artificial media from inner tissues of bark and needles. Antimicrobial and antioxidant activity, along with total phenolic- and flavonoid-content assays were used in the evaluation of bioactivity of the fermented crude extracts. The ability of the endophytes to produce the anticancer compound Taxol was also analyzed using thin-layer chromatography (TLC) and reverse-phase high-performance liquid chromatography (RP-HPLC). A total of 16 fungal morphotypes were obtained from asymptomatic inner tissues of the bark and needles of T. wallichiana. Among the 16 isolates, the ethyl acetate (EA) fraction of isolate MUS1, showed antibacterial and antifungal activity against all test-pathogens used (Streptococcus faecalis ATCC 19433, Staphylococcus aureus ATCC 12600, Bacillus subtilis ATCC 6633, Escherichia coli ATCC 25922, Salmonella enterica ATCC 13076, Pseudomonas aeruginosa ATCC 27853, and Candida albicans). MUS1 showed significant inhibition against Pseudomonas aeruginosa ATCC 27853 (minimum inhibitory concentration (MIC): 250 µg/ml) and the pathogenic yeast, Candida albicans (MIC: 125 µg/ml). Antioxidant activity, total phenolic, and total flavonoid content as well as in vitro Taxol production were evaluated for EA fraction of isolate MUS1. Antioxidant activity was evaluated using 2,2-diphenyl-1-picrylhydrazyl (DPPH) assay. At a concentration of 100 µg/ml, the % DPPH radical scavenging activity was 83.15 ± 0.40, 81.62 ± 0.11, and 62.36 ± 0.29, for ascorbic acid, butylated hydroxytoluene (BHT), and the EA fraction of MUS1, respectively. The DPPH-Half maximal inhibitory concentration (DPPH-IC50) value for the EA fraction was 81.52 ± 0.23 µg/ml, compared to BHT (62.87 ± 0.08 µg/ml) and ascorbic acid (56.15 ± 0.19 µg/ml). The total phenolic and flavonoid content in the EA fraction were 16.90 ± 0.075 µg gallic acid equivalent (GAE) and 11.59 ± 0.148 µg rutin equivalent (RE), per mg of dry crude extract, respectively. TLC and RP-HPLC analysis showed that the isolate MUS1 also produces Taxol (282.05 µg/l of fermentation broth). Isolate MUS1 was identified as Annulohypoxylon sp. by internal transcribed spacer (ITS) sequencing. Having the ability to produce antimicrobial and antioxidant metabolites, as well as the anticancer compound Taxol, makes Annulohypoxylon sp. strain MUS1, a promising candidate for further study of naturally occurring bioactive metabolites. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02693-z.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce.

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model 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.

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.

Overexpression of PaNAC03, a stress induced NAC gene family transcription factor in Norway spruce leads to reduced flavonol biosynthesis and aberrant embryo development.

BACKGROUND: The NAC family of transcription factors is one of the largest gene families of transcription factors in plants and the conifer NAC gene family is at least as large, or possibly larger, as in Arabidopsis. These transcription factors control both developmental and stress induced processes in plants. Yet, conifer NACs controlling stress induced processes has received relatively little attention. This study investigates NAC family transcription factors involved in the responses to the pathogen Heterobasidion annosum (Fr.) Bref. sensu lato. RESULTS: The phylogeny and domain structure in the NAC proteins can be used to organize functional specificities, several well characterized stress-related NAC proteins are found in III-3 in Arabidopsis (Jensen et al. Biochem J 426:183-196, 2010). The Norway spruce genome contain seven genes with similarity to subgroup III-3 NACs. Based on the expression pattern PaNAC03 was selected for detailed analyses. Norway spruce lines overexpressing PaNAC03 exhibited aberrant embryo development in response to maturation initiation and 482 misregulated genes were identified in proliferating cultures. Three key genes in the flavonoid biosynthesis pathway: a CHS, a F3'H and PaLAR3 were consistently down regulated in the overexpression lines. In accordance, the overexpression lines showed reduced levels of specific flavonoids, suggesting that PaNAC03 act as a repressor of this pathway, possibly by directly interacting with the promoter of the repressed genes. However, transactivation studies of PaNAC03 and PaLAR3 in Nicotiana benthamiana showed that PaNAC03 activated PaLAR3A, suggesting that PaNAC03 does not act as an independent negative regulator of flavan-3-ol production through direct interaction with the target flavonoid biosynthetic genes. CONCLUSIONS: PaNAC03 and its orthologs form a sister group to well characterized stress-related angiosperm NAC genes and at least PaNAC03 is responsive to biotic stress and appear to act in the control of defence associated secondary metabolite production.

Transcriptional Responses Associated with Virulence and Defence in the Interaction between Heterobasidion annosum s.s. and Norway Spruce.

Heterobasidion annosum sensu lato is a serious pathogen causing root and stem rot to conifers in the northern hemisphere and rendering the timber defective for sawing and pulping. In this study we applied next-generation sequencing to i) identify transcriptional responses unique to Heterobasidion-inoculated Norway spruce and ii) investigate the H. annosum transcripts to identify putative virulence factors. To address these objectives we wounded or inoculated 30-year-old Norway spruce clones with H. annosum and 454-sequenced the transcriptome of the interaction at 0, 5 and 15 days post inoculation. The 491,860 high-quality reads were de novo assembled and the relative expression was analysed. Overall, very few H. annosum transcripts were represented in our dataset. Three delta-12 fatty acid desaturase transcripts and one Clavaminate synthase-like transcript, both associated with virulence in other pathosystems, were found among the significantly induced transcripts. The analysis of the Norway spruce transcriptional responses produced a handful of differentially expressed transcripts. Most of these transcripts originated from genes known to respond to H. annosum. However, three genes that had not previously been reported to respond to H. annosum showed specific induction to inoculation: an oxophytodienoic acid-reductase (OPR), a beta-glucosidase and a germin-like protein (GLP2) gene. Even in a small data set like ours, five novel highly expressed Norway spruce transcripts without significant alignment to any previously annotated protein in Genbank but present in the P. abies (v1.0) gene catalogue were identified. Their expression pattern suggests a role in defence. Therefore a more complete survey of the transcriptional responses in the interactions between Norway spruce and its major pathogen H. annosum would probably provide a better understanding of gymnosperm defence than accumulated until now.

The primary module in Norway spruce defence signalling against H. annosum s.l. seems to be jasmonate-mediated signalling without antagonism of salicylate-mediated signalling.

A key tree species for the forest industry in Europe is Norway spruce [Picea abies (L.) Karst.]. One of its major diseases is stem and butt rot caused by Heterobasidion parviporum (Fr.) Niemelä & Korhonen, which causes extensive revenue losses every year. In this study, we investigated the parallel induction of Norway spruce genes presumably associated with salicylic acid- and jasmonic acid/ethylene-mediated signalling pathways previously observed in response to H. parviporum. Relative gene expression levels in bark samples of genes involved in the salicylic acid- and jasmonic acid/ethylene-mediated signalling pathways after wounding and inoculation with either the saprotrophic biocontrol fungus Phlebiopsis gigantea or with H. parviporum were analysed with quantitative PCR at the site of the wound and at two distal locations from the wound/inoculation site to evaluate their roles in the induced defence response to H. parviporum in Norway spruce. Treatment of Norway spruce seedlings with methylsalicylate, methyljasmonate and inhibitors of the jasmonic acid/ethylene signalling pathway, as well as the Phenylalanine ammonia lyase inhibitor 2-aminoindan-2-phosphonic acid were conducted to determine the responsiveness of genes characteristic of the different pathways to different hormonal stimuli. The data suggest that jasmonic acid-mediated signalling plays a central role in the induction of the genes analysed in this study irrespective of their responsiveness to salicylic acid. This may suggest that jasmonic acid-mediated signalling is the prioritized module in the Norway spruce defence signalling network against H. parviporum and that there seems to be no immediate antagonism between the modules in this interaction.

Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen.

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.

Expression analysis of Clavata1-like and Nodulin21-like genes from Pinus sylvestris during ectomycorrhiza formation.

The ecology and physiology of ectomycorrhizal (EcM) symbiosis with conifer trees are well documented. In comparison, however, very little is known about the molecular regulation of these associations. In an earlier study, we identified three EcM-regulated Pinus expressed sequence tags (EST), two of which were identified as homologous to the Medicago truncatula nodulin MtN21. The third EST was a homologue to the receptor-like kinase Clavata1. We have characterized the expression patterns of these genes and of auxin- and mycorrhiza-regulated genes after induction with indole-3-butyric acid in Pinus sylvestris and in a time course experiment during ectomycorrhizal initiation with the co-inoculation of 2,3,5-triiodobenzoic acid, an auxin transport inhibitor. Our results suggest that different P. sylvestris nodulin homologues are associated with diverse processes in the root. The results also suggest a potential role of the Clv1-like gene in lateral root initiation by the ectomycorrhizal fungus.

Chemical and transcriptional responses of Norway spruce genotypes with different susceptibility to Heterobasidion spp. infection.

BACKGROUND: Norway spruce [Picea abies (L.) Karst.] is one of the most important conifer species in Europe. The wood is economically important and infections by wood-rotting fungi cause substantial losses to the industry.The first line of defence in a Norway spruce tree is the bark. It is a very efficient barrier against infection based on its mechanical and chemical properties. Once an injury or an infection is recognized by the tree, induced defences are activated. In this study we examined transcriptional response, using 454-sequencing, and chemical profiles in bark of Norway spruce trees with different susceptibility to Heterobasidion annosum s.l. infection. The aim was to find associations between the transcriptome and chemical profiles to the level of susceptibility to Heterobasidion spp. in Norway spruce genotypes. RESULTS: Both terpene and phenol compositions were analysed and at 28 days post inoculation (dpi) high levels of 3-carene was produced in response to H. annosum. However, significant patterns relating to inoculation or to genotypes with higher or lower susceptibility could only be found in the phenol fraction. The levels of the flavonoid catechin, which is polymerized into proanthocyanidins (PA), showed a temporal variation; it accumulated between 5 and 15 dpi in response to H. annosum infection in the less susceptible genotypes. The transcriptome data suggested that the accumulation of free catechin was preceded by an induction of genes in the flavonoid and PA biosynthesis pathway such as leucoanthocyanidin reductase. Quantitative PCR analyses verified the induction of genes in the phenylpropanoid and flavonoid pathway. The qPCR data also highlighted genotype-dependent differences in the transcriptional regulation of these pathways. CONCLUSIONS: The varying dynamics in transcriptional and chemical patterns displayed by the less susceptible genotypes suggest that there is a genotypic variation in successful spruce defence strategies against Heterobasidion. However, both high levels of piceasides and flavonoids in the less susceptible genotypes suggested the importance of the phenolic compounds in the defence. Clearly an extended comparison of the transcriptional responses in the interaction with Heterobasidion between several independent genotypes exhibiting reduced susceptibility is needed to catalogue mechanisms of successful host defence strategies.