Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

Microbial communities in the phyllosphere and endosphere of Norway spruce under attack by Heterobasidion

Heterobasidion annosum species complex has been regarded as the most destructive disease agent of conifer trees in boreal forests. Tree microbiome can regulate the plant-pathogen interactions by influencing both host resistance and pathogen virulence. Such information would help to improve the future health of forests and explore strategies to enhance ecosystem stability. In this study, using next-generation sequencing technology, we investigated the microbial community in different tree regions (needles, upper stem, and lower stem) of Norway spruce with and without wood decay symptoms. The primary purpose was to uncover signature characteristic microbiome harbored by asymptomatic trees compared to diseased trees. Additionally, the study was to explore the inter-kingdom and intra-kingdom interactions in microbiome (bacteria and fungi) of symptomatic versus asymptomatic trees. The results showed that in upper stem, species richness (Chao1) of fungi and bacteria were both higher in asymptomatic trees than symptomatic trees (P < 0.05). Compared to symptomatic trees, asymptomatic trees harbored a higher abundance of Actinobacteriota, bacterial genera of Methylocella, Conexibacter, Jatrophihabitans, and fungal genera of Mollisia. Fungal communities from the same anatomic region differed between the symptomatic and asymptomatic trees. Bacterial communities from the two stem regions were also distinct between the symptomatic and asymptomatic trees. The symptomatic trees possessed a less stable microbial network with more positive correlations compared to the asymptomatic trees. In the lower stem, at intra-kingdom level, the distribution of correlation numbers was more even in the bacterial network compared to the fungal network. In conclusion, the Heterobasidion attack decreased the microbial community species richness and shifted the community structure and functional structure to varying degrees. The microbial network was enlarged and became more unstable at both inter-kingdom and intra-kingdom level due to the Heterobasidion infection.

Early-Stage Infection-Specific Heterobasidion annosum (Fr.) Bref. Transcripts in H. annosum-Pinus sylvestris L. Pathosystem

Transcriptomes from stem-inoculated Scots pine saplings were analyzed to identify unique and enriched H. annosum transcripts in the early stages of infection. Comparing different time points since inoculation identified 131 differentially expressed H. annosum genes with p-values of ≤0.01. Our research supports the results of previous studies on the Norway spruce-Heterobasidion annosum s.l. pathosystem, indicating the role of carbohydrate and lignin degradation genes in pathogenesis at different time points post-inoculation and the role of lipid metabolism genes (including but not limited to the delta-12 fatty acid desaturase gene previously reported to be an important factor). The results of this study indicate that the malic enzyme could be a potential gene of interest in the context of H. annosum virulence. During this study, difficulties related to incomplete reference material of the host plant species and a low proportion of H. annosum transcripts in the RNA pool were encountered. In addition, H. annosum transcripts are currently not well annotated. Improvements in sequencing technologies (including sequencing depth) or bioinformatics focusing on small subpopulations of RNA would be welcome.

Genome sequence of Heterobasidion occidentale, a fungus that causes annosus root and butt rot among conifer trees in North America

We report an annotated draft genome of Heterobasidion occidentale, a fungus (Basidiomycota, Agaricomycetes) that has pathogenic and saprophytic lifestyles. This fungus belongs to the H. annosum (Fr.) Bref. sensu lato species complex that comprises several root rot pathogens. Heterobasidion occidentale causes annosus root and butt rot primarily in true fir (Abies spp.) and spruce (Picea spp.) species throughout western North America.

Patterns of methylation and transcriptional plasticity during thermal acclimation in a reef-building coral

Phenotypic plasticity can buffer organisms against short-term environmental fluctuations. For example, previous exposure to increased temperatures can increase thermal tolerance in many species. Prior studies have found that acclimation to higher temperature can influence the magnitude of transcriptional response to subsequent acute thermal stress (hereafter, "transcriptional response modulation"). However, mechanisms mediating this gene expression response and, ultimately, phenotypic plasticity remain largely unknown. Epigenetic modifications are good candidates for modulating transcriptional response, as they broadly correlate with gene expression. Here, we investigate changes in DNA methylation as a possible mechanism controlling shifts in gene expression plasticity and thermal acclimation in the reef-building coral Acropora nana. We find that gene expression response to acute stress is altered in corals acclimated to different temperatures, with many genes exhibiting a dampened response to heat stress in corals pre-conditioned to higher temperatures. At the same time, we observe shifts in methylation during both acclimation (11 days) and acute heat stress (24 h). We observed that the acute heat stress results in shifts in gene-level methylation and elicits an acute transcriptional response in distinct gene sets. Further, acclimation-induced shifts in gene expression plasticity and differential methylation also largely occur in separate sets of genes. Counter to our initial hypothesis no overall correlation between the magnitude of differential methylation and the change in gene expression plasticity. We do find a small but statistically significant overlap in genes exhibiting both dampened expression response and shifts in methylation (14 genes), which could be candidates for further inquiry. Overall, our results suggest transcriptional response modulation occurs independently from methylation changes induced by thermal acclimation.

Phenotypic variation and genomic variation in insect virulence traits reveal patterns of intraspecific diversity in a locust-specific fungal pathogen

Intraspecific pathogen diversity is crucial for understanding the evolution and maintenance of adaptation in host-pathogen interactions. Traits associated with virulence are often a significant source of variation directly impacted by local selection pressures. The specialist fungal entomopathogen, Metarhizium acridum, has been widely implemented as a biological control agent of locust pests in tropical regions of the world. However, few studies have accounted for natural intraspecific phenotypic and genetic variation. Here, we examine the diversity of nine isolates of M. acridum spanning the known geographic distribution, in terms of (1) virulence towards two locust species, (2) growth rates on three diverse nutrient sources, and (3) comparative genomics to uncover genomic variability. Significant variability in patterns of virulence and growth was shown among the isolates, suggesting intraspecific ecological specialization. Different patterns of virulence were shown between the two locust species, indicative of potential host preference. Additionally, a high level of diversity among M. acridum isolates was observed, revealing increased variation in subtilisin-like proteases from the Pr1 family. These results culminate in the first in-depth analysis regarding multiple facets of natural variation in M. acridum, offering opportunities to understand critical evolutionary drivers of intraspecific diversity in pathogens.

Genomic and Secretomic Analyses of the Newly Isolated Fungus Perenniporia fraxinea SS3 Identified CAZymes Potentially Related to a Serious Pathogenesis of Hardwood Trees

Perenniporia fraxinea can colonize living trees and cause severe damage to standing hardwoods by secreting a number of carbohydrate-activate enzymes (CAZymes), unlike other well-studied Polyporales. However, significant knowledge gaps exist in understanding the detailed mechanisms for this hardwood-pathogenic fungus. To address this issue, five monokaryotic P. fraxinea strains, SS1 to SS5, were isolated from the tree species Robinia pseudoacacia, and high polysaccharide-degrading activities and the fastest growth were found for P. fraxinea SS3 among the isolates. The whole genome of P. fraxinea SS3 was sequenced, and its unique CAZyme potential for tree pathogenicity was determined in comparison to the genomes of other nonpathogenic Polyporales. These CAZyme features are well conserved in a distantly related tree pathogen, Heterobasidion annosum. Furthermore, the carbon source-dependent CAZyme secretions of P. fraxinea SS3 and a nonpathogenic and strong white-rot Polyporales member, Phanerochaete chrysosporium RP78, were compared by activity measurements and proteomic analyses. As seen in the genome comparisons, P. fraxinea SS3 exhibited higher pectin-degrading activities and higher laccase activities than P. chrysosporium RP78, which were attributed to the secretion of abundant glycoside hydrolase family 28 (GH28) pectinases and auxiliary activity family 1_1 (AA1_1) laccases, respectively. These enzymes are possibly related to fungal invasion into the tree lumens and the detoxification of tree defense substances. Additionally, P. fraxinea SS3 showed secondary cell wall degradation capabilities at the same level as that of P. chrysosporium RP78. Overall, this study suggested mechanisms for how this fungus can attack the cell walls of living trees as a serious pathogen and differs from other nonpathogenic white-rot fungi. IMPORTANCE Many studies have been done to understand the mechanisms underlying the degradation of plant cell walls of dead trees by wood decay fungi. However, little is known about how some of these fungi weaken living trees as pathogens. P. fraxinea belongs to the Polyporales, a group of strong wood decayers, and is known to aggressively attack and fell standing hardwood trees all over the world. Here, we report CAZymes potentially related to plant cell wall degradation and pathogenesis factors in a newly isolated fungus, P. fraxinea SS3, by genome sequencing in conjunction with comparative genomic and secretomic analyses. The present study provides insights into the mechanisms of the degradation of standing hardwood trees by the tree pathogen, which will contribute to the prevention of this serious tree disease.

Intraspecific Comparative Analysis Reveals Genomic Variation of Didymella arachidicola and Pathogenicity Factors Potentially Related to Lesion Phenotype

Didymella arachidicola is one of the most important fungal pathogens, causing foliar disease and leading to severe yield losses of peanuts (Arachis hypogaea L.) in China. Two main lesion phenotypes of peanut web blotch have been identified as reticulation type (R type) and blotch type (B type). As no satisfactory reference genome is available, the genomic variations and pathogenicity factors of D. arachidicola remain to be revealed. In the present study, we collected 41 D. arachidicola isolates from 26 geographic locations across China (33 for R type and 8 for B type). The chromosome-scale genome of the most virulent isolate (YY187) was assembled as a reference using PacBio and Hi-C technologies. In addition, we re-sequenced 40 isolates from different sampling sites. Genome-wide alignments showed high similarity among the genomic sequences from the 40 isolates, with an average mapping rate of 97.38%. An average of 3242 SNPs and 315 InDels were identified in the genomic variation analysis, which revealed an intraspecific polymorphism in D. arachidicola. The comparative analysis of the most and least virulent isolates generated an integrated gene set containing 512 differential genes. Moreover, 225 genes individually or simultaneously harbored hits in CAZy-base, PHI-base, DFVF, etc. Compared with the R type reference, the differential gene sets from all B type isolates identified 13 shared genes potentially related to lesion phenotype. Our results reveal the intraspecific genomic variation of D. arachidicola isolates and pathogenicity factors potentially related to different lesion phenotypes. This work sets a genomic foundation for understanding the mechanisms behind genomic diversity driving different pathogenic phenotypes of D. arachidicola.

Genomic Comparisons of Two Armillaria Species with Different Ecological Behaviors and Their Associated Soil Microbial Communities

Armillaria species show considerable variation in ecological roles and virulence, from mycorrhizae and saprophytes to important root pathogens of trees and horticultural crops. We studied two Armillaria species that can be found in coniferous forests of northwestern USA and southwestern Canada. Armillaria altimontana not only is considered as a weak, opportunistic pathogen of coniferous trees, but it also appears to exhibit in situ biological control against A. solidipes, formerly North American A. ostoyae, which is considered a virulent pathogen of coniferous trees. Here, we describe their genome assemblies and present a functional annotation of the predicted genes and proteins for the two Armillaria species that exhibit contrasting ecological roles. In addition, the soil microbial communities were examined in association with the two Armillaria species within a 45-year-old plantation of western white pine (Pinus monticola) in northern Idaho, USA, where A. altimontana was associated with improved tree growth and survival, while A. solidipes was associated with reduced growth and survival. The results from this study reveal a high similarity between the genomes of the beneficial/non-pathogenic A. altimontana and pathogenic A. solidipes; however, many relatively small differences in gene content were identified that could contribute to differences in ecological lifestyles and interactions with woody hosts and soil microbial communities.

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.

The dark septate endophyte Phialocephala sphaeroides confers growth fitness benefits and mitigates pathogenic effects of Heterobasidion on Norway spruce

Forest trees frequently interact with a diverse range of microorganisms including dark septate endophytes (DSEs) and fungal pathogens. Plant defense responses to either individual pathogens or endophytes have been widely studied, but very little is known on the effect of coinfection on host defenses. To study the impact of coinfection or tripartite interaction on plant growth and host defenses, Norway spruce (Picea abies (L.) Karst) seedlings were inoculated with a DSE Phialocephala sphaeroides or with a root pathogen Heterobasidion parviporum Niemela & Korhonen or coinfected with both fungi. The results showed that the DSE promoted the root growth of spruce seedlings. Control seedlings without any inoculum were subjected to sequencing and used as a baseline for identification of differentially expressed genes (DEGs). RNA-seq analysis of seedlings inoculated with P. sphaeroides, infected with H. parviporum or coinfected with both fungi resulted in a total of 5269 DEGs. The majority of DEGs were found in P. sphaeroides-inoculated seedlings. Lignin biosynthesis pathways were generally activated during fungal infections. The pattern was distinct with endophyte inoculation. The majority of the genes in the flavonoid biosynthesis pathway were generally suppressed during fungal infections. A specific transcriptional response to P. sphaeroides inoculation was the increased transcripts of genes involved in jasmonic acid biosynthesis, mitogen-activated protein kinases signaling pathway, plant hormone signal transduction and calcium-mediated signaling. This may have potentially contributed to promoting the root growth of seedlings. Although the coinfection suppressed the induction of numerous genes, no negative effect on the growth of the spruce seedlings occurred. We conclude that the subsequent H. parviporum infection triggered reprogramming of host metabolism. Conversely, the endophyte (P. sphaeroides), on the other hand, counteracted the negative effects of H. parviporum on the growth of the spruce seedlings.

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.

Phenotypic Recovery of a Heterobasidion Isolate Infected by a Debilitation-Associated Virus Is Related to Altered Host Gene Expression and Reduced Virus Titer

The fungal genus Heterobasidion includes forest pathogenic species hosting a diverse group of partitiviruses. They include the host debilitating Heterobasidion partitivirus 13 strain an1 (HetPV13-an1), which was originally observed in a slowly growing H. annosum strain 94233. In this study, a relatively fast-growing sector strain 94233-RC3 was isolated from a highly debilitated mycelial culture of 94233, and its gene expression and virus transcript quantities as well as the genomic sequence of HetPV13-an1 were examined. The sequence of HetPV13-an1 genome in 94233-RC3 was identical to that in the original 94233, and thus not the reason for the partial phenotypic recovery. According to RNA-seq analysis, the HetPV13-an1 infected 94233-RC3 transcribed eight genes differently from the partitivirus-free 94233-32D. Three of these genes were downregulated and five upregulated. The number of differentially expressed genes was considerably lower and the changes in their expression were small compared to those of the highly debilitated original strain 94233 with the exception of the most highly upregulated ones, and therefore viral effects on the host transcriptome correlated with the degree of the virus-caused debilitation. The amounts of RdRp and CP transcripts of HetPV13-an1 were considerably lower in 94233-RC3 and also in 94233 strain infected by a closely related mildly debilitating virus HetPV13-an2, suggesting that the virus titer would have a role in determining the effect of HetPV13 viruses on their hosts.

Fungi Occurring in Norway Spruce Wood Decayed by Heterobasidion parviporum in Puszcza Borecka Stands (Northeastern Poland)

In many spruce stands, trees are frequently attacked by the pathogen Heterobasidion parviporum, albeit without visible symptoms in the crown. In the present work, the results of the presence of stem rot, assessed by PICUS Sonic Tomography, and the fungal biota on trees and stumps in eight plots in the Puszcza Borecka Forest are described. The plots were located in stands on original forest soil (4) and on post-agricultural soil (4), where around a stump with H. parviporum symptoms (signs of internal rot and basidiocarps), 30 trees were selected and examined for internal rot. Wood samples were collected from two selected trees for fungal molecular analysis. A total of 79 fungal taxa were found, including 57 taxa in plots on post-agricultural soil and 45 on forest soil. There were 395 fungal records on stumps and 22 records on trees, therein, from the inner parts of felled trunks. Significant differences in the Chao-1 diversity index indicate that the origin of the soil—post-agricultural or forest soil—influenced the alpha diversity of the fungal communities in the forests studied. The values of the Shannon and Simpson indices show that the two communities were similar in terms of species numbers. The presence of basidiomata of H. parviporum and two species of Armillaria (mainly A. cepistipes) in samples on all plots is striking, although Armillaria spp. was detected more frequently. Most of the species identified were typical saprotrophs, although rare species were also found, such as Entoloma byssisedum, Onnia tomentosa, Physisporinus vitreus, Postia ptychogaster, and Ramaria apiculata. The presence of H. parviporum in the inner woody parts was confirmed by PCR analysis, and decay was detected even up to a stem height of 6 m. Armillaria was the dominant genus in the studied stands, and H. parviporum plays a significant and underestimated role in heartwood decay of old spruce trees in Puszcza Borecka Forest.

Genome Sequencing and Analysis of the Fungal Symbiont of Sirex noctilio, Amylostereum areolatum: Revealing the Biology of Fungus-Insect Mutualism

Sirex noctilio (F.), together with Amylostereum areolatum, a wood-decaying symbiotic fungus, causes severe damage to Pinus species worldwide. In China, it causes extensive death of Mongolian pine (Pinus sylvestris var. mongolica). There is an obligate dependency mutualism between the woodwasp and its fungus. Studies have suggested that the fungal growth rate affected the size of the wasps: larger adults emerged from sites with a higher fungus growth rate. This genome is the first reported genome sequence of a woodwasp symbiotic fungus. Genome sequence analysis of this fungus would prove the possibility of A. areolatum volatiles affecting the host selection of S. noctilio on a molecular basis. We further clarified that A. areolatum was a strict obligate symbiotic fungus and that it would provide S. noctilio with a suitable environment and with nutrients for the larval growth. These results would lay a foundation for our understanding of the mechanism of this entomogenous symbiosis.

Amylostereum areolatum is the symbiotic fungus of the Eurasian woodwasp, Sirex noctilio, a globally invasive species. The mutualistic symbiont is associated with the woodwasp, assisting the damage process and providing nutrition for its insect partners. Colonization and growth of A. areolatum have essential impacts on the development and spread of S. noctilio, though the mechanism of interaction between the two has been poorly described. In this study, the first genome of this symbiotic fungus was sequenced, assembled, and annotated. The assembled A. areolatum genome was 57.5 Mb (54.51% GC content) with 15,611 protein-coding genes. We identified 580 carbohydrate-active enzymes (CAZymes), 661 genes associated with pathogen-host interactions, and 318 genes encoding transport proteins in total. The genome annotation revealed 10 terpene/phytoene synthases responsible for terpenoid biosynthesis, which could be classified into three clades. Terpene synthase gene clusters in clade II were conserved well across Russulales. In this cluster, genes encoding mevalonate kinase (MK), EGR12 (COG1557), and nonplant terpene cyclases (cd00687) were the known biosynthesis and regulatory genes. Genome sequence analysis of this fungus would prove the possibility of A. areolatum volatiles affecting the host selection of S. noctilio on a molecular basis. We further clarified that A. areolatum was a strict obligate symbiotic fungus. The wasps might protect the fungus before it was introduced into a suitable host substrate by oviposition, while the fungus would provide S. noctilio with a suitable environment and nutrients for the larval growth. These results would lay a foundation for our understanding of the mechanism of this entomogenous symbiosis.

IMPORTANCESirex noctilio (F.), together with Amylostereum areolatum, a wood-decaying symbiotic fungus, causes severe damage to Pinus species worldwide. In China, it causes extensive death of Mongolian pine (Pinus sylvestris var. mongolica). There is an obligate dependency mutualism between the woodwasp and its fungus. Studies have suggested that the fungal growth rate affected the size of the wasps: larger adults emerged from sites with a higher fungus growth rate. This genome is the first reported genome sequence of a woodwasp symbiotic fungus. Genome sequence analysis of this fungus would prove the possibility of A. areolatum volatiles affecting the host selection of S. noctilio on a molecular basis. We further clarified that A. areolatum was a strict obligate symbiotic fungus and that it would provide S. noctilio with a suitable environment and with nutrients for the larval growth. These results would lay a foundation for our understanding of the mechanism of this entomogenous symbiosis.

The Pattern and Function of DNA Methylation in Fungal Plant Pathogens

To successfully infect plants and trigger disease, fungal plant pathogens use various strategies that are dependent on characteristics of their biology and genomes. Although pathogenic fungi are different from animals and plants in the genomic heritability, sequence feature, and epigenetic modification, an increasing number of phytopathogenic fungi have been demonstrated to share DNA methyltransferases (MTases) responsible for DNA methylation with animals and plants. Fungal plant pathogens predominantly possess four types of DNA MTase homologs, including DIM-2, DNMT1, DNMT5, and RID. Numerous studies have indicated that DNA methylation in phytopathogenic fungi mainly distributes in transposable elements (TEs), gene promoter regions, and the repetitive DNA sequences. As an important and heritable epigenetic modification, DNA methylation is associated with silencing of gene expression and transposon, and it is responsible for a wide range of biological phenomena in fungi. This review highlights the relevant reports and insights into the important roles of DNA methylation in the modulation of development, pathogenicity, and secondary metabolism of fungal plant pathogens. Recent evidences prove that there are massive links between DNA and histone methylation in fungi, and they commonly regulate fungal development and mycotoxin biosynthesis.

An Italian Research Culture Collection of Wood Decay Fungi

One of the main aims of the University of Pavia mycology laboratory was to collect wood decay fungal (WDF) strains in order to deepen taxonomic studies, species distribution, officinal properties or to investigate potential applications such as biocomposite material production based on fungi. The Italian Alps, Apennines and wood plains were investigated to collect Basidiomycota basidiomata from living or dead trees. The purpose of this study was to investigate the wood decay strains of the Mediterranean area, selecting sampling sites in North and Central Italy, including forests near the Ligurian and Adriatic seas, or near the Lombardy lakes. The isolation of mycelia in pure culture was performed according to the current methodology and the identity of the strains was confirmed by molecular analyses. The strains are maintained in the Research Culture Collection MicUNIPV of Pavia University (Italy). Among the 500 WDF strains in the collection, the most interesting isolates from the Mediterranean area are: Dichomitus squalens (basidioma collected from Pinus pinea), Hericium erinaceus (medicinal mushroom), Inocutis tamaricis (white-rot agent on Tamarix trees), Perenniporia meridionalis (wood degrader through Mn peroxidase) and P. ochroleuca. In addition, strains of species related to the Mediterranean climate (e.g., Fomitiporia mediterranea and Cellulariella warnieri) were obtained from sites with a continental-temperate climate.

Potential Interactions between Invasive Fusarium circinatum and Other Pine Pathogens in Europe

Pines are major components of native forests and plantations in Europe, where they have both economic significance and an important ecological role. Diseases of pines are mainly caused by fungal and oomycete pathogens, and can significantly reduce the survival, vigor, and yield of both individual trees and entire stands or plantations. Pine pitch canker (PPC), caused by Fusarium circinatum (Nirenberg and O’Donnell), is among the most devastating pine diseases in the world, and is an example of an emergent invasive disease in Europe. The effects of microbial interactions on plant health, as well as the possible roles plant microbiomes may have in disease expression, have been the focus of several recent studies. Here, we describe the possible effects of co-infection with pathogenic fungi and oomycetes with F. circinatum on the health of pine seedlings and mature plants, in an attempt to expand our understanding of the role that biotic interactions may play in the future of PPC disease in European nurseries and forests. The available information on pine pathogens that are able to co-occur with F. circinatum in Europe is here reviewed and interpreted to theoretically predict the effects of such co-occurrences on pine survival, growth, and yield. Beside the awareness that F. circinatum may co-occurr on pines with other pathogens, an additional outcome from this review is an updating of the literature, including the so-called grey literature, to document the geographical distribution of the relevant pathogens and to facilitate differential diagnoses, particularly in nurseries, where some of them may cause symptoms similar to those induced by F. circinatum. An early and accurate diagnosis of F. circinatum, a pathogen that has been recently introduced and that is currently regulated in Europe, is essential to prevent its introduction and spread in plantings and forests.

The Conifer Root and Stem Rot Pathogen (Heterobasidion parviporum): Effectome Analysis and Roles in Interspecific Fungal Interactions

Heterobasidion parviporum Niemelä & Korhonen is an economically important basidiomycete, causing root and stem rot disease of Norway spruce (Picea abies (L.) Karst) in Northern Europe. The H. parviporum genome encodes numerous small secreted proteins, which might be of importance for interacting with mycorrhiza symbionts, endophytes, and other saprotrophs. We hypothesized that small secreted proteins from H. parviporum (HpSSPs) are involved in interspecific fungal interaction. To identify HpSSP-coding genes potentially involved, we screened the H. parviporum effectome and compared their transcriptomic profiles during fungal development and in planta tree infection. We further conducted phylogenetic analysis, and identified a subset of hypothetical proteins with nonpredicted domain or unknown function as HpSSPs candidates for further characterization. The HpSSPs candidates were selected based on high-quality sequence, cysteine residue frequency, protein size, and in planta expression. We subsequently explored their roles during in vitro interaction in paired cultures of H. parviporum with ectomycorrhizal Cortinarius gentilis, endophytic Phialocephala sphaeroides, saprotrophs (Mycena sp., Phlebiopsis gigantea, and Phanerochaete chrysosporium), respectively. The transcriptomic profile revealed that a large proportion of effector candidates was either barely expressed or highly expressed under all growth conditions. In vitro dual-culture test showed that P. sphaeroides and C. gentilis were overgrown by H. parviporum. The barrage zone formation or no physical contact observed in paired cultures with the saprotrophs suggest they had either combative interaction or antibiosis effect with H. parviporum. Several HpSSPs individuals were up- or downregulated during the nonself interactions. The results of HpSSPs gene expression patterns provide additional insights into the diverse roles of SSPs in tree infection and interspecific fungal interactions.

Fourier-transform infrared (FT-IR) spectroscopy analysis discriminates asymptomatic and symptomatic Norway spruce trees

Conifer trees, including Norway spruce, are threatened by fungi of the Heterobasidion annosum species complex, which severely affect timber quality and cause economic losses to forest owners. The timely detection of infected trees is complicated, as the pathogen resides within the heartwood and sapwood of infected trees. The presence of the disease and the extent of the wood decay often becomes evident only after tree felling. Fourier-transform infrared (FT-IR) spectroscopy is a potential method for non-destructive sample analysis that may be useful for identifying infected trees in this pathosystem. We performed FT-IR analysis of 18 phloem, 18 xylem, and 18 needle samples from asymptomatic and symptomatic Norway spruce trees. FT-IR spectra from 1066 - 912 cm^-1 could be used to distinguish phloem, xylem, and needle tissue extracts. FT-IR spectra collected from xylem and needle extracts could also be used to discriminate between asymptomatic and symptomatic trees using spectral bands from 1657 - 994 cm^-1 and 1104 - 994 cm^-1, respectively. A partial least squares regression model predicted the concentration of condensed tannins, a defense-related compound, in phloem of asymptomatic and symptomatic trees. This work is the first to show that FT-IR spectroscopy can be used for the identification of Norway spruce trees naturally infected with Heterobasidion spp.

Genome-Wide Signatures of Selection in Colletotrichum kahawae Reveal Candidate Genes Potentially Involved in Pathogenicity and Aggressiveness

Plants and their pathogens are engaged in continuous evolutionary battles, with pathogens evolving to circumvent plant defense mechanisms and plants responding through enhanced protection to prevent or mitigate damage induced by pathogen attack. Managed ecosystems are composed of genetically identical populations of crop plants with few changes from year to year. These environments are highly conducive to the emergence and dissemination of pathogens and they exert selective pressure for both qualitative virulence factors responsible for fungal pathogenicity, and quantitative traits linked to pathogen fitness, such as aggressiveness. In this study, we used a comparative genome-wide approach to investigate the genomic basis underlying the pathogenicity and aggressiveness of the fungal coffee pathogen Colletotrichum kahawae infecting green coffee berries. The pathogenicity was investigated by comparing genomic variation between C. kahawae and its non-pathogenic sibling species, while the aggressiveness was studied by a genome-wide association approach with groups of isolates with different phenotypic profiles. High genetic differentiation was observed between C. kahawae and the most closely related species with 5,560 diagnostic SNPs identified, in which a significant enrichment of non-synonymous mutations was detected. Functional annotation of these non-synonymous mutations revealed a significant enrichment mainly in two gene ontology categories, "oxidation-reduction process" and "integral component of membrane." Finally, the annotation of several genes potentially under-selection revealed that C. kahawae's pathogenicity may be a complex biological process, in which important biological functions, such as, detoxification and transport, regulation of host and pathogen gene expression, and signaling are involved. On the other hand, the genome-wide association analyses for aggressiveness were able to identify 10 SNPs and 15 SNPs of small effect in single and multi-association analysis, respectively, from which 7 were common, giving in total 18 SNPs potentially associated. The annotation of these genomic regions allowed the identification of four candidate genes encoding F-box domain-containing, nitrosoguanidine resistance, Fungal specific transcription factor domain-containing and C6 transcription factor that could be associated with aggressiveness. This study shed light, for the first time, on the genetic mechanisms of C. kahawae host specialization.

Chlorophyll fluorescence imaging for monitoring effects of Heterobasidion parviporum small secreted protein induced cell death and in planta defense gene expression

Heterobasidion parviporum Niemelä & Korhonen is a necrotrophic fungal pathogen of Norway spruce (Picea abies). The H. parviporum genome encodes numerous necrotrophic small secreted proteins (SSP) which might be important for promoting and sustaining the disease development. However, their transcriptional dynamics and plant defense response during infection are largely unknown. In this study, we identified a necrotrophic SSP named HpSSP35.8 and its coding gene was highly expressed in the pre-symptomatic phase of the host (Norway spruce) infection. We explored the impact of HpSSP35.8 on non-host Nicotiana benthamiana using Agrobacterium-mediated transient expression system under visible spectrum RGB imaging and chlorophyll fluorescence imaging. The results showed that HpSSP35.8 triggered a form of SSP-associated programmed cell death, accompanied by a decrease in the plant photosynthetic activity. Defense-related genes including WRKY12, ethylene response factor (ERF1α) and a chitinase gene PR4 were up-regulated in both HpSSP35.8-N. benthamiana interaction and H. parviporum-Norway spruce pathosystem. This study also highlighted the potential to use the chlorophyll fluorescence imaging approach to monitor both the indirect effects of SSP and also for the selection of other potential effector-like protein candidates.

Dual RNA-seq analysis provides new insights into interactions between Norway spruce and necrotrophic pathogen Heterobasidion annosum s.l

BACKGROUND

Root and butt rot of conifer trees caused by fungi belonging to the Heterobasidion annosum species complex is one of the most economically important fungal diseases in commercial conifer plantations throughout the Northern hemisphere. We investigated the interactions between Heterobasidion fungi and their host by conducting dual RNA-seq and chemical analysis on Norway spruce trees naturally infected by Heterobasidion spp. We analyzed host and pathogen transcriptome and phenolic and terpenoid contents of the spruce trees.

RESULTS

Presented results emphasize the role of the phenylpropanoid and flavonoid pathways in the chemical defense of Norway spruce trees. Accumulation of lignans was observed in trees displaying symptoms of wood decay. A number of candidate genes with a predicted role in the higher level regulation of spruce defense responses were identified. Our data indicate a possible role of abscisic acid (ABA) signaling in the spruce defense against Heterobasidion infection. Fungal transcripts corresponding to genes encoding carbohydrate- and lignin-degrading enzymes, secondary metabolism genes and effector-like genes were expressed during the host colonization.

CONCLUSIONS

Our results provide additional insight into defense strategies employed by Norway spruce trees against Heterobasidion infection. The potential applications of the identified candidate genes as markers for higher resistance against root and butt rot deserve further evaluation.

Genome-wide DNA methylation and transcriptomic profiles in the lifestyle strategies and asexual development of the forest fungal pathogen Heterobasidion parviporum

Heterobasidion parviporum is the most devastating fungal pathogen of conifer forests in Northern Europe. The fungus has dual life strategies, necrotrophy on living trees and saprotrophy on dead woods. DNA cytosine methylation is an important epigenetic modification in eukaryotic organisms. Our presumption is that the lifestyle transition and asexual development in H. parviporum could be driven by epigenetic effects. Involvements of DNA methylation in the regulation of aforementioned processes have never been studied thus far. RNA-seq identified lists of highly induced genes enriched in carbohydrate-active enzymes during necrotrophic interaction with host trees and saprotrophic sawdust growth. It also highlighted signaling- and transcription factor-related genes potentially associated with the transition of saprotrophic to necrotrophic lifestyle and groups of primary cellular activities throughout asexual development. Whole-genome bisulfite sequencing revealed that DNA methylation displayed pronounced preference in CpG dinucleotide context across the genome and mostly targeted transposable element (TE)-rich regions. TE methylation level demonstrated a strong negative correlation with TE expression, reinforcing the protective function of DNA methylation in fungal genome stability. Small groups of genes putatively subject to methylation transcriptional regulation in response to saprotrophic and necrotrophic growth in comparison with free-living mycelia were also explored. Our study reported on the first methylome map of a forest pathogen. Analysis of transcriptome and methylome variations associated with asexual development and different lifestyle strategies provided further understanding of basic biological processes in H. parviporum. More importantly, our work raised additional potential roles of DNA methylation in fungi apart from controlling the proliferation of TEs.