Search

Recently, attention has been shifting toward the perspective of the existence of plants and microbes as a functioning ecological unit. However, studies highlighting the impacts of the microbial community on plant health are still limited. In this study, fungal community (mycobiome) of leaf, root, and soil of symptomatic leaf-spot diseased (SS) oil palm were compared against asymptomatic (AS) trees using ITS2 rRNA gene metabarcoding. A total of 3,435,417 high-quality sequences were obtained from 29 samples investigated. Out of the 14 phyla identified, Ascomycota and Basidiomycota were the most dominant accounting for 94.2 and 4.7% of the total counts in AS, and 75 and 21.2% in SS, respectively. Neopestalotiopsis is the most abundant genus for AS representing 8.0% of the identified amplicons compared to 2.0% in SS while Peniophora is the most abundant with 8.6% of the identified amplicons for SS compared to 0.1% in AS. The biomarker discovery algorithm LEfSe revealed different taxa signatures for the sample categories, particularly soil samples from asymptomatic trees, which were the most enriched. Network analysis revealed high modularity across all groups, except in root samples. Additionally, a large proportion of the identified keystone species consisted of rare taxa, suggesting potential role in ecosystem functions. Surprisingly both AS and SS leaf samples shared taxa previously associated with oil palm leaf spot disease. The significant abundance of Trichoderma asperellum in the asymptomatic root samples could be further explored as a potential biocontrol agent against oil palm disease.

Rigidoporus microporus, the causative agent of the white root rot disease of rubber trees, poses a significant threat to natural rubber production worldwide. Understanding the molecular mechanisms facilitating its pathogenicity would be crucial for developing effective disease management strategies. The pathogen secretes effector proteins, which play pivotal roles in modulating host immune responses and infection. In this study, in silico analyses identified 357 putative secreted effector proteins from the R. microporus genome. These were then integrated into previous RNA-seq data obtained in response to rubber tree latex exposure. Annotation of putative effectors suggested the abundance of proteins in several families associated with the virulence of R. microporus, especially hydrophobin proteins and glycoside hydrolase (GH) proteins. The contribution of secreted effectors to fungal pathogenicity was discussed, particularly in response to rubber tree latex exposure. Some unknown highly expressed effectors were predicted for the protein structures, revealing their similarity to aminopeptidase, ubiquitin ligase, spherulin, and thaumatin protein. This integrative study further elucidates the molecular mechanism of R. microporus pathogenesis and offers alternative targets for developing control strategies for managing white root rot disease in rubber plantations.

In recent years, the importance of microorganisms for plant survival has been increasingly recognized. Endophytic fungi, as part of holobiont, can confer growth advantages to plants. Most studies have shown that the endophytic fungi of forest trees can promote host plant growth, increase adversity resistance, and thus improve the survival competitiveness of forest trees. However, the beneficial examples of endophytic fungi on the growth and development of woody plants have not been systematically summarized. This review is focused on various aspects of beneficial endophytic fungi in forest trees (definition, classification, colonization mechanisms, etc.), with an emphasis on their beneficial roles in woody plant growth, protection against biotic and abiotic stresses, as well as the response of forest trees to endophytic fungi. In addition, this review lists a series of experiments on screening beneficial endophytic fungi from Chinese fir (Cunninghamia lanceolata) and verifying their beneficial functions, to explore the mutualistic relationships between them. This review not only provides a theoretical basis for the study of beneficial endophytic fungi in forest trees in the future but also sheds light on the molecular perspectives for a mechanistic understanding of their potential future significance for the sustainable utilization of forest resources and ecological environment protection.

Abstract The natural world is under unprecedented and accelerating pressure. Much work on understanding resilience to local and global environmental change has, so far, focussed on ecosystems. However, understanding a system’s behaviour requires knowledge of its component parts and their interactions. Here we call for increased efforts to understand ‘biological resilience’, or the processes that enable components across biological levels, from genes to communities, to resist or recover from perturbations. Although ecologists and evolutionary biologists have the tool-boxes to examine form and function, efforts to integrate this knowledge across biological levels and take advantage of big data (e.g. ecological and genomic) are only just beginning. We argue that combining eco-evolutionary knowledge with ecosystem-level concepts of resilience will provide the mechanistic basis necessary to improve management of human, natural and agricultural ecosystems, and outline some of the challenges in achieving an understanding of biological resilience.

Root and stem rot caused by Heterobasidion annosum is a severe problem in boreal Scots pine. Dissecting the features of disease resistance is generally an essential step in resistance breeding in plants and forest trees. In this study, we explored inherent resistance factors of Scots pine against H. annosum. A total of 236 families consisting of 85 full-sib (FS), 35 half-sib population mix (HSpm), and 116 half-sib (HS) families of Scots pine seedlings were inoculated with a H. annosum isolate. We sampled needle tissues before inoculation for terpene measurements and RNA sequencing. Based on the lesion area, the extremes of 12 resistant and 12 susceptible families were selected for further analyses. Necrotic lesions resulting from fungal infection were in a weak to moderate relationship with the plant height. Monoterpenes were the principal terpene compounds observed in Scots pine seedlings. Concentrations of 3-carene were significantly higher in pine genotypes inherently resistant compared with susceptible seedlings. By contrast, susceptible genotypes had significantly higher proportions of α-pinene. Gene ontology analysis of differential expressed transcripts (DETs) revealed that response to biotic factors was enriched in resistant seedlings. Functional characterization of individual DETs revealed that higher expression of transcripts involved in response to abiotic stress was common in susceptible genotypes. This observation was supported by the annotation of hub genes in a key module that was significantly correlated with the lesion trait through weighted gene co-expression network analysis (WGCNA) of 16 HS and HSpm samples. These findings contribute to our understanding of constitutive resistance factors of Scots pine against Heterobasidion root and stem rot diseases.

Several economically important diseases of forest trees and agricultural crops in many parts of the world have been linked to the ascomycete fungal pathogen Thielaviopsis paradoxa. This study compared the growth rate of 41 isolates of T. paradoxa sourced from different hosts and two countries (Nigeria and Papua New Guinea (PNG)) under six temperature levels (22 °C, 25 °C, 30 °C, 32 °C, 34 °C and 35 °C). Phylogenetic relationships were obtained from the analysis of their nuclear ribosomal DNA internal transcribed sequence (ITS) data. While all the isolates from PNG and few from Nigeria grew optimally between 22 °C and 32 °C, the majority had their highest growth rate (2.9 cm/day) between 25 °C and 32 °C. Growth performances were generally low between 34 °C and 35 °C; no isolate from the sugar cane grew at these high temperatures. The oil palm isolate DA029 was the most resilient, with the highest growth rate (0.97 cm/day) at 35 °C. Phylogenetic analysis delineated five clusters: a very large clade which accommodates the majority (30 Nigerian and 3 PNG oil palm isolates) and four small clades containing two members each. To a large extent, the clustering pattern failed to address the temperature–isolate relationship observed. However, only the four small clades represent isolates with similar temperature tolerances. It is most likely that wider and robust analyses with more diverse isolates and genetic markers will provide better insight on thermal resilience of T. paradoxa. Additionally, future research to establish relationships between vegetative growth at different temperatures and of different pathogenicity and disease epidemiology merits being explored. The results might provide useful information for the formulation of effective management and control strategies against the pathogen, especially in this era of climate change.

The microbiome of Heterobasidion-induced wood decay of living trees has been previously studied; however, less is known about the bacteria biota of its perennial fruiting body and the adhering wood tissue. In this study, we investigated the bacteria biota of the Heterobasidion fruiting body and its adhering deadwood. Out of 7,462 operational taxonomic units (OTUs), about 5,918 OTUs were obtained from the fruiting body and 5,469 OTUs were obtained from the associated dead wood. Interestingly, an average of 52.6% of bacteria biota in the fruiting body was shared with the associated dead wood. The overall and unique OTUs had trends of decreasing from decay classes 1 to 3 but increasing in decay class 4. The fruiting body had the highest overall and unique OTUs number in the fourth decay class, whereas wood had the highest OTU in decay class 1. Sphingomonas spp. was significantly higher in the fruiting body, and phylum Firmicutes was more dominant in wood tissue. The FAPROTAX functional structure analysis revealed nutrition, energy, degradation, and plant-pathogen-related functions of the communities. Our results also showed that bacteria communities in both substrates experienced a process of a new community reconstruction through the various decay stages. The process was not synchronic in the two substrates, but the community structures and functions were well-differentiated in the final decay class. The bacteria community was highly dynamic; the microbiota activeness, community stability, and functions changed with the decay process. The third decay class was an important turning point for community restructuring. Host properties, environmental factors, and microbial interactions jointly influenced the final community structure. Bacteria community in the fruiting body attached to the living standing tree was suppressed compared with those associated with dead wood. Bacteria appear to spread from wood tissue of the standing living tree to the fruiting body, but after the tree is killed, bacteria moved from fruiting body to wood. It is most likely that some of the resident endophytic bacteria within the fruiting body are either parasitic, depending on it for their nutrition, or are mutualistic symbionts.

Laccases (EC 1.10.3.2), a group of multi-copper oxidases (MCOs), play multiple biological functions and widely exist in many species. Fungal laccases have been extensively studied for their industrial applications, however, there was no database specially focused on fungal laccases. To provide a comparative genomics platform for fungal laccases, we have developed a comparative genomics platform for laccases and MCOs (http://laccase.riceblast.snu.ac.kr/). Based on protein domain profiles of characterized sequences, 3,571 laccases were predicted from 690 genomes including 253 fungi. The number of putative laccases and their properties exhibited dynamic distribution across the taxonomy. A total of 505 laccases from 68 genomes were selected and subjected to phylogenetic analysis. As a result, four clades comprised of nine subclades were phylogenetically grouped by their putative functions and analyzed at the sequence level. Our work would provide a workbench for putative laccases mainly focused on the fungal kingdom as well as a new perspective in the identification and classification of putative laccases and MCOs.

Abstract Background The white rot fungus Phlebia radiata, a type species of the genus Phlebia, is an efficient decomposer of plant cell wall polysaccharides, modifier of softwood and hardwood lignin, and is able to produce ethanol from various waste lignocellulose substrates. Thus, P. radiata is a promising organism for biotechnological applications aiming at sustainable utilization of plant biomass. Here we report the genome sequence of P. radiata isolate 79 originally isolated from decayed alder wood in South Finland. To better understand the evolution of wood decay mechanisms in this fungus and the Polyporales phlebioid clade, gene content and clustering of genes encoding specific carbohydrate-active enzymes (CAZymes) in seven closely related fungal species was investigated. In addition, other genes encoding proteins reflecting the fungal lifestyle including peptidases, transporters, small secreted proteins and genes involved in secondary metabolism were identified in the genome assembly of P. radiata. Results The PACBio sequenced nuclear genome of P. radiata was assembled to 93 contigs with 72X sequencing coverage and annotated, revealing a dense genome of 40.4 Mbp with approximately 14 082 predicted protein-coding genes. According to functional annotation, the genome harbors 209 glycoside hydrolase, 27 carbohydrate esterase, 8 polysaccharide lyase, and over 70 auxiliary redox enzyme-encoding genes. Comparisons with the genomes of other phlebioid fungi revealed shared and specific properties among the species with seemingly similar saprobic wood-decay lifestyles. Clustering of especially GH10 and AA9 enzyme-encoding genes according to genomic localization was discovered to be conserved among the phlebioid species. In P. radiata genome, a rich repertoire of genes involved in the production of secondary metabolites was recognized. In addition, 49 genes encoding predicted ABC proteins were identified in P. radiata genome together with 336 genes encoding peptidases, and 430 genes encoding small secreted proteins. Conclusions The genome assembly of P. radiata contains wide array of carbohydrate polymer attacking CAZyme and oxidoreductase genes in a composition identifiable for phlebioid white rot lifestyle in wood decomposition, and may thus serve as reference for further studies. Comparative genomics also contributed to enlightening fungal decay mechanisms in conversion and cycling of recalcitrant organic carbon in the forest ecosystems.

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.

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

Plant-beneficial microbes have drawn wide attention due to their potential application as bio-control agents and bio-fertilizers. Moso bamboo, which is among the monocots with the highest growth rate, lives perennially with abundant microbes that may benefit annually growing crops. Genome information of moso bamboo associated bacteria remains underexplored. We isolated and identified a novel Paraburkholderia strain Suichang626 from moso bamboo roots. Growth promoting effects of Suichang626 on both moso bamboo and seedlings of the model dicot Arabidopsis thaliana were documented in laboratory conditions. To gain insight into the genetic basis of this growth promotion effect, we sequenced the genome of Suichang626. Evidenced by genome-wide phylogeny data, we propose that Suichang626 is a novel strain of Paraburkholderia sacchari. Gene homologs encoding biosynthesis of the plant growth-promoting chemicals, acetoin and 2,3-butanediol, were identified in the genome of Suichang626. Comparative genomics was further performed with plant-beneficial and plant/animal pathogenic species of Paraburkholderia and Burkholderia. Genes related to volatile organic compounds, nitrogen fixation, and auxin biosynthesis were discovered specifically in the plant growth-promoting species of both genera.

Abstract Background Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species. Results To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates. Conclusion Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

The causal agent of root and butt rot of conifer trees, Heterobasidion annosum, is widespread in boreal forests and economically responsible for annual loss of approximately 50 million euros to forest industries in Finland alone and much more at European level. In order to further understand the pathobiology of this fungus at the genome level, a Finnish isolate of H. annosum sensu stricto (isolate 03012) was sequenced and analyzed with the genome sequences of 23 white-rot and 13 brown-rot fungi. The draft genome assembly of H. annosum has a size of 31.01 Mb, containing 11,453 predicted genes. Whole genome alignment showed that 84.38% of H. annosum genome sequences were aligned with those of previously sequenced H. irregulare TC 32-1 counterparts. The result is further supported by the protein sequence clustering analysis which revealed that the two genomes share 6719 out of 8647 clusters. When sequencing reads of H. annosum were aligned against the genome sequences of H. irregulare, six single nucleotide polymorphisms were found in every 1 kb, on average. In addition, 98.68% of SNPs were found to be homo-variants, suggesting that the two species have long evolved from different niches. Gene family analysis revealed that most of the white-rot fungi investigated had more gene families involved in lignin degradation or modification, including laccases and peroxidase. Comparative analysis of the two Heterobasidion spp. as well as white-/brown-rot fungi would provide new insights for understanding the pathobiology of the conifer tree pathogen.

Abstract The basidiomycete Heterobasidion annosum sensu lato (s.l.) is considered to be one of the most destructive conifer pathogens in the temperate forests of the northern hemisphere. H. annosum is characterized by a dual fungal lifestyle. The fungus grows necrotrophically on living plant cells and saprotrophically on dead wood material. In this study, we screened the H. annosum genome for small secreted proteins (HaSSPs) that could potentially be involved in promoting necrotrophic growth during the fungal infection process. The final list included 58 HaSSPs that lacked predictable protein domains. The transient expression of HaSSP encoding genes revealed the ability of 8 HaSSPs to induce cell chlorosis and cell death in Nicotiana benthamiana. In particular, one protein (HaSSP30) could induce a rapid, strong, and consistent cell death within 2 days post-infiltration. HaSSP30 also increased the transcription of host-defence-related genes in N. benthamiana, which suggested a necrotrophic-specific immune response. This is the first line of evidence demonstrating that the H. annosum genome encodes HaSSPs with the capability to induce plant cell death in a non-host plant.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pine wilt disease (PWD) caused by the nematode Bursaphelenchus xylophilus is a devastating disease in conifer forests in Eurasia. However, information on the effect of PWD on the host microbial community is limited. In this study, the bacterial community structure and potential function in the needles, roots, and soil of diseased pine were studied under field conditions using Illumina MiSeq coupled with Phylogenetic Investigation of Communities by Reconstruction of Unobserved states (PICRUSt) software. The results showed that the community and functional structure of healthy and diseased trees differed only in the roots and needles, respectively (p < 0.05). The needles, roots, and soil formed unique bacterial community and functional structures. The abundant phyla across all samples were Proteobacteria (41.9% of total sequence), Actinobacteria (29.0%), Acidobacteria (12.2%), Bacteroidetes (4.8%), and Planctomycetes (2.1%). The bacterial community in the healthy roots was dominated by Acidobacteria, Planctomycetes, and Rhizobiales, whereas in the diseased roots, Proteobacteria, Firmicutes, and Burkholderiales were dominant. Functionally, groups involved in the cell process and genetic information processing had a higher abundance in the diseased needles, which contributed to the difference in functional structure. The results indicate that PWD can only affect the host bacteria community structure and function in certain anatomical regions of the host tree.

Root and butt rot caused by members of the Heterobasidion annosum species complex is the most economically important disease of conifer trees in boreal forests. Wood decay in the infected trees dramatically decreases their value and causes considerable losses to forest owners. Trees vary in their susceptibility to Heterobasidion infection, but the genetic determinants underlying the variation in the susceptibility are not well-understood. We performed the identification of Norway spruce genes associated with the resistance to Heterobasidion parviporum infection using genome-wide exon-capture approach. Sixty-four clonal Norway spruce lines were phenotyped, and their responses to H. parviporum inoculation were determined by lesion length measurements. Afterwards, the spruce lines were genotyped by targeted resequencing and identification of genetic variants (SNPs). Genome-wide association analysis identified 10 SNPs located within 8 genes as significantly associated with the larger necrotic lesions in response to H. parviporum inoculation. The genetic variants identified in our analysis are potential marker candidates for future screening programs aiming at the differentiation of disease-susceptible and resistant trees.

Climate changes, exemplified by increased temperatures and CO2 concentration, pose a global threat to forest health. Of particular concern are pests and pathogens, with a warming climate altering their distributions and evolutionary capacity, while impairing the ability of some plants to respond to infections. Progress in understanding and mitigating such effects is currently hindered by a lack of empirical research. Norway spruce (Picea abies) is one of the most economically important tree species in northern Europe, and is considered highly vulnerable to changes in climate. It is commonly infected by the fungus Endoconidiophora polonica, and we hypothesized that damage caused to trees will increase under future climate change predictions. To test this hypothesis an in vivo greenhouse experiment was conducted to evaluate the effects of a changed growing environment on E. polonica infected Norway spruce seedlings, comparing ambient conditions to predicted temperatures and CO2 levels in Finland for the years 2030 and 2100. In total, 450 seedlings were randomized amongst the three treatments, with 25 seedlings from each allocated to inoculation with one of five different fungal strains or mock-inoculation. Seedlings were monitored throughout the thermal growing season for mortality, and lesion length and depth indices were measured at the experiment conclusion. Disease severity (mortality and lesions) was consistently greater in fungal-inoculated than mock-inoculated seedlings. However, substantial differences were observed among fungal strains in response to climate scenarios. For example, although overall seedling mortality was highest under the most distant (and severe) climate change expectations, of the two fungal strains with the highest mortality counts (referred to as F4 and F5), one produced greater mortality under the 2030 and 2100 scenarios than ambient conditions, whereas climate scenario had no effect on the other. This study contributes to a limited body of empirical research on the effects of projected climate changes on forestry pathosystems, and is the first to investigate interactions between Norway spruce and E. polonica. The results indicate the potential for future climate changes to alter the impact of forest pathogens with implications for productivity, while highlighting the need for a strain-specific level of understanding of the disease agents.