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

Transcriptional Profiling of Early Defense Response to White Pine Blister Rust Infection in Pinus albicaulis (Whitebark Pine)

Pathogen perception generates the activation of signal transduction cascades to host defense. White pine blister rust (WPBR) is caused by Cronartium ribicola J.C. Fisch and affects a number of species of Pinus. One of the most severely affected species is Pinus albicaulis Engelm (whitebark pine). WPBR resistance in the species is a polygenic and complex trait that requires an optimized immune response. We identified early responses in 2-year-old seedlings after four days of fungal inoculation and compared the underlying transcriptomic response with that of healthy non-inoculated individuals. A de novo transcriptome assembly was constructed with 56,796 high quality-annotations derived from the needles of susceptible and resistant individuals in a resistant half-sib family. Differential expression analysis identified 599 differentially expressed transcripts, from which 375 were upregulated and 224 were downregulated in the inoculated seedlings. These included components of the initial phase of active responses to abiotic factors and stress regulators, such as those involved in the first steps of flavonoid biosynthesis. Four days after the inoculation, infected individuals showed an overexpression of chitinases, reactive oxygen species (ROS) regulation signaling, and flavonoid intermediates. Our research sheds light on the first stage of infection and emergence of disease symptoms among whitebark pine seedlings. RNA sequencing (RNA-seq) data encoding hypersensitive response, cell wall modification, oxidative regulation signaling, programmed cell death, and plant innate immunity were differentially expressed during the defense response against C. ribicola.

Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Vertical and horizontal gene transfer shaped plant colonization and biomass degradation in the fungal genus Armillaria

The fungal genus Armillaria contains necrotrophic pathogens and some of the largest terrestrial organisms that cause tremendous losses in diverse ecosystems, yet how they evolved pathogenicity in a clade of dominantly non-pathogenic wood degraders remains elusive. Here we show that Armillaria species, in addition to gene duplications and de novo gene origins, acquired at least 1,025 genes via 124 horizontal gene transfer events, primarily from Ascomycota. Horizontal gene transfer might have affected plant biomass degrading and virulence abilities of Armillaria, and provides an explanation for their unusual, soft rot-like wood decay strategy. Combined multi-species expression data revealed extensive regulation of horizontally acquired and wood-decay related genes, putative virulence factors and two novel conserved pathogenicity-induced small secreted proteins, which induced necrosis in planta. Overall, this study details how evolution knitted together horizontally and vertically inherited genes in complex adaptive traits of plant biomass degradation and pathogenicity in important fungal pathogens.

Hybrid Genome Assembly of Berkeleyomyces rouxiae, an Emerging Cannabis Fungal Pathogen Causing Black Root Rot in an Aeroponic Facility

The resurged interest in cultivation of Cannabis sativa has presented an array of new challenges. Among them are the difficult-to-control pests and pathogens that infect cannabis plants. The limited methods for disease control available to cannabis growers necessitates early detection of plant pathogens, something that molecular techniques such as DNA sequencing has greatly improved. This study reports for the first time the fungal plant pathogen Berkeleyomyces rouxiae causing black root rot in high THC-containing cannabis. Aeroponically grown cannabis plants at a licenced production facility in Cranbrook BC, Canada, rapidly displayed root discoloration and rot symptoms despite testing negative for all commercially available pathogen tests. Developing sequencing-based disease diagnostics requires genomic information, so this study presents the first whole genome sequence of the multihost, widespread black root rot pathogen B. rouxiae. Hybrid genome assembly using Oxford Nanopore long-reads and Illumina short-reads yielded a genome size of 28.2 Mb represented over 404 contigs with an N50 of 267 kb. Genome annotation predicted 6,960 protein-coding genes with 59,477 functional annotations. The availability of this genome will assist in sequence-based diagnostic development, comparative genomics, and taxonomic resolution of this globally important plant pathogen.

Lignin-degrading enzymes from a pathogenic canker-rot fungus Inonotus obliquus strain IO-B2

Inonotus obliquus is a pathogenic fungus found in living trees and has been widely used as a traditional medicine for cancer therapy. Although lignocellulose-degrading enzymes are involved in the early stages of host infection, the parasitic life cycle of this fungus has not been fully understood. In this study, we aimed to investigate the activities of laccase (Lac), manganese peroxidase (MnP), and lignin peroxidase (LiP) from I. obliquus cultivated in Kirk's medium. The fungus was subjected to genome sequencing, and genes related to wood degradation were identified. The draft genome sequence of this fungus comprised 21,203 predicted protein-coding genes, of which 134 were estimated to be related to wood degradation. Among these, 47 genes associated with lignin degradation were found to have the highest number of mnp genes. Furthermore, we cloned the cDNA encoding a putative MnP, referred to as IoMnP1, and characterized its molecular structure. The results show that IoMnP1 has catalytic properties analogous to MnP. Phylogenetic analysis also confirmed that IoMnP1 was closely related to the MnPs from Pyrrhoderma noxium, Fomitiporia mediterranea, and Sanghuangporus baumii, which belong to the same family of Hymenochaetaceae. From the above results, we suggest that IoMnP1 is a member of MnPs.

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.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

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

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

SUPPLEMENTARY INFORMATION

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

Multicopper oxidases with laccase-ferroxidase activity: Classification and study of ferroxidase activity determinants in a member from Heterobasidion annosum s. l

Multi-copper oxidases (MCO) share a common molecular architecture and the use of copper ions as cofactors to reduce O₂ to H₂O, but show high sequence heterogeneity and functional diversity. Many new emerging MCO genes are wrongly annotated as laccases, the largest group of MCOs, with the widest range of biotechnological applications (particularly those from basidiomycete fungi) due to their ability to oxidise aromatic compounds and lignin. Thus, comprehensive studies for a better classification and structure-function characterisation of MCO families are required. Laccase-ferroxidases (LAC-FOXs) constitute a separate and unexplored group of MCOs with proposed dual features between laccases and ferroxidases. We aim to better define this cluster and the structural determinants underlying putative hybrid activity. We performed a phylogenetic analysis of the LAC-FOXs from basidiomycete fungi, that resulted in two subgroups. This division seemed to correlate with the presence or absence of some of the three acidic residues responsible for ferroxidase activity in Fet3p from Saccharomyces cerevisiae. One of these LAC-FOXs (with only one of these residues) from the fungus Heterobasidion annosum s. l. (HaLF) was synthesised, heterologously expressed and characterised to evaluate its catalytic activity. HaLF oxidised typical laccase substrates (phenols, aryl amines and N-heterocycles), but no Fe (II). The enzyme was subjected to site-directed mutagenesis to determine the key residues that confer ferroxidase activity. The mutated HaLF variant with full restoration of the three acidic residues exhibited efficient ferroxidase activity, while it partially retained the wide-range oxidative activity of the native enzyme associated to laccases sensu stricto.

Taxol and β-tubulins from endophytic fungi isolated from the Himalayan Yew, Taxus wallichiana Zucc

Paclitaxel, better known as the anticancer drug Taxol®, has been isolated from several plant species and has been shown to be produced by fungi, actinomycetes, and even bacteria isolated from marine macroalgae. Given its cytostatic effect, studies conducted in the 1990's showed that paclitaxel was toxic to many pathogenic fungi and oomycetes. Further studies led to the idea that the differences in paclitaxel sensitivity exhibited by different fungi were due to differences in the β-tubulin protein sequence. With the recent isolation of endophytic fungi from the leaves and bark of the Himalayan Yew, Taxus wallichiana Zucc., and the availability of genomes from paclitaxel-producing fungi, we decided to further explore the idea that endophytic fungi isolated from Yews should be well-adapted to their environment by encoding β-tubulin proteins that are insensitive to paclitaxel. Our results found evidence of episodic positive/diversifying selection at 10 sites (default p-value threshold of 0.1) in the β-tubulin sequences, corresponding to codon positions 33, 55, 172, 218, 279, 335, 359, 362, 379, and 406. Four of these positions (i.e., 172, 279, 359, and 362) have been implicated in the binding of paclitaxel by β-tubulin or formed part of the binding pocket. As expected, all the fungal endophytes grew in different media regardless of the paclitaxel concentration tested. Furthermore, our results also showed that Taxomyces andreanae CBS 279.92, the first fungus shown to produce paclitaxel, is a Basidiomycete fungus as the two beta tubulins encoded by the fungus clustered together with other Basidiomycete fungi.

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.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

Lytic polysaccharide monooxygenases (LPMOs) producing microbes: A novel approach for rapid recycling of agricultural wastes

Out of the huge quantity of agricultural wastes produced globally, rice straw is one of the most abundant ligno-cellulosic waste. For efficient utilization of these wastes, several cost-effective biological processes are available. The practice of field level in-situ or ex-situ decomposition of rice straw is having less degree of adoption due to its poor decomposition ability within a short time span between rice harvest and sowing of the next crop. Agricultural wastes including rice straw are in general utilized by using lignocellulose degrading microbes for industrial metabolite or compost production. However, bioconversion of crystalline cellulose and lignin present in the waste, into simple molecules is a challenging task. To resolve this issue, researchers have identified a novel new generation microbial enzyme i.e., lytic polysaccharide monooxygenases (LPMOs) and reported that the combination of LPMOs with other glycolytic enzymes are found efficient. This review explains the progress made in LPMOs and their role in lignocellulose bioconversion and the possibility of exploring LPMOs producers for rapid decomposition of agricultural wastes. Also, it provides insights to identify the knowledge gaps in improving the potential of the existing ligno-cellulolytic microbial consortium for efficient utilization of agricultural wastes at industrial and field levels.

Discovery and Community Dynamics of Novel ssRNA Mycoviruses in the Conifer Pathogen Heterobasidion parviporum

Heterobasidion species are highly destructive basidiomycetous conifer pathogens of the Boreal forest region. Earlier studies have revealed dsRNA virus infections of families Curvulaviridae and Partitiviridae in Heterobasidion strains, and small RNA deep sequencing has also identified infections of Mitoviridae members in these fungi. In this study, the virome of Heterobasidion parviporum was examined for the first time by RNA-Seq using total RNA depleted of rRNA. This method successfully revealed new viruses representing two established (+)ssRNA virus families not found earlier in Heterobasidion: Narnaviridae and Botourmiaviridae. In addition, we identified the presence of a recently described virus group tentatively named “ambiviruses” in H. parviporum. The H. parviporum isolates included in the study originated from experimental forest sites located within 0.7 km range from each other, and a population analysis including 43 isolates was conducted at one of the experimental plots to establish the prevalence of the newly identified viruses in clonally spreading H. parviporum individuals. Our results indicate that viral infections are considerably more diverse and common among Heterobasidion isolates than known earlier and include ssRNA viruses with high prevalence and interspecies variation.

Mitonuclear Genetic Interactions in the Basidiomycete Heterobasidion parviporum Involve a Non-conserved Mitochondrial Open Reading Frame

The mitochondrial and nuclear genomes of Eukaryotes are inherited separately and consequently follow distinct evolutionary paths. Nevertheless, the encoding of many mitochondrial proteins by the nuclear genome shows the high level of integration they have reached, which makes mitonuclear genetic interactions all the more conceivable. For each species, natural selection has fostered the evolution of coadapted alleles in both genomes, but a population-wise divergence of such alleles could lead to important phenotypic variation, and, ultimately, to speciation. In this study in the Basidiomycete Heterobasidion parviporum, we have investigated the genetic basis of phenotypic variation among laboratory-designed heterokaryons carrying the same pair of haploid nuclei, but a different mitochondrial genome. Radial growth rate data of thirteen unrelated homokaryotic parents and of their heterokaryotic offspring were combined with SNP data extracted from parental genome sequences to identify nuclear and mitochondrial loci involved in mitonuclear interactions. Two nuclear loci encoding mitochondrial proteins appeared as best candidates to engage in a genetic interaction affecting radial growth rate with a non-conserved mitochondrial open reading frame of unknown function and not reported apart from the Russulales order of Basidiomycete fungi. We believe our approach could be useful to investigate several important traits of fungal biology where mitonuclear interactions play a role, including virulence of fungal pathogens.

Genomic and Experimental Investigations of Auriscalpium and Strobilurus Fungi Reveal New Insights into Pinecone Decomposition

Saprophytic fungi (SPF) play vital roles in ecosystem dynamics and decomposition. However, because of the complexity of living systems, our understanding of how SPF interact with each other to decompose organic matter is very limited. Here we studied their roles and interactions in the decomposition of highly specialized substrates between the two genera Auriscalpium and Strobilurus fungi-colonized fallen pinecones of the same plant sequentially. We obtained the genome sequences from seven fungal species with three pairs: A. orientale-S. luchuensis, A. vulgare-S. stephanocystis and A. microsporum-S. pachcystidiatus/S. orientalis on cones of Pinus yunnanensis, P. sylvestris and P. armandii, respectively, and the organic profiles of substrate during decomposition. Our analyses revealed evidence for both competition and cooperation between the two groups of fungi during decomposition, enabling efficient utilization of substrates with complementary profiles of carbohydrate active enzymes (CAZymes). The Auriscalpium fungi are highly effective at utilizing the primary organic carbon, such as lignin, and hemicellulose in freshly fallen cones, facilitated the invasion and colonization by Strobilurus fungi. The Strobilurus fungi have genes coding for abundant CAZymes to utilize the remaining organic compounds and for producing an arsenal of secondary metabolites such as strobilurins that can inhibit other fungi from colonizing the pinecones.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

Draft Genome Sequence of the Ectomycorrhizal Fungus Astraeus odoratus from Northern Thailand

We report the draft genome sequence of Astraeus odoratus, an edible ectomycorrhizal fungus from northern Thailand. The assembled genome has a size of 45.1 Mb and 13,403 annotated protein-coding genes. This reference genome will provide a better understanding of the biology of mushroom-forming ectomycorrhizal fungi in the family Diplocystidiaceae.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Basidiomycetes Are Particularly Sensitive to Bacterial Volatile Compounds: Mechanistic Insight Into the Case Study of Pseudomonas protegens Volatilome Against Heterobasidion abietinum

Volatile organic compounds (VOCs) play an important role in the communication among organisms, including plants, beneficial or pathogenic microbes, and pests. In vitro, we observed that the growth of seven out of eight Basidiomycete species tested was inhibited by the VOCs of the biocontrol agent Pseudomonas protegens strain CHA0. In the Ascomycota phylum, only some species were sensitive (e.g., Sclerotinia sclerotiorum, Botrytis cinerea, etc.) but others were resistant (e.g., Fusarium oxysporum f. sp. cubense, Verticillium dahliae, etc.). We further discovered that CHA0 as well as other ten beneficial or phytopathogenic bacterial strains were all able to inhibit Heterobasidion abietinum, which was used in this research as a model species. Moreover, such an inhibition occurred only when bacteria grew on media containing digested proteins like peptone or tryptone (e.g., Luria-Bertani agar or LBA). Also, the inhibition co-occurred with a pH increase of the agar medium where the fungus grew. Therefore, biogenic ammonia originating from protein degradation by bacteria was hypothesized to play a major role in fungus inhibition. Indeed, when tested as a synthetic compound, it was highly toxic to H. abietinum (effective concentration 50% or EC₅₀ = 1.18 M; minimum inhibitory concentration or MIC = 2.14 M). Using gas chromatography coupled to mass spectrometry (GC/MS), eight VOCs were found specifically emitted by CHA0 grown on LBA compared to the bacterium grown on potato dextrose agar (PDA). Among them, two compounds were even more toxic than ammonia against H. abietinum: dimethyl trisulfide had EC₅₀ = 0.02 M and MIC = 0.2 M, and 2-ethylhexanol had EC₅₀ = 0.33 M and MIC = 0.77 M. The fungus growth inhibition was the result of severe cellular and sub-cellular alterations of hyphae occurring as early as 15 min of exposure to VOCs, as evidenced by transmission and scanning electron microscopy observations. Transcriptome reprogramming of H. abietinum induced by CHA0’s VOCs pointed out that detrimental effects occurred on ribosomes and protein synthesis while the cells tried to react by activating defense mechanisms, which required a lot of energy diverted from the growth and development (fitness cost).

Genic introgression from an invasive exotic fungal forest pathogen increases the establishment potential of a sibling native pathogen

Significant hybridization between the invasive North American fungal plant pathogen Heterobasidion irregulare and its Eurasian sister species H. annosum is ongoing in Italy. Whole genomes of nine natural hybrids were sequenced, assembled and compared with those of three genotypes each of the two parental species. Genetic relationships among hybrids and their level of admixture were determined. A multi-approach pipeline was used to assign introgressed genomic blocks to each of the two species. Alleles that introgressed from H. irregulare to H. annosum were associated with pathways putatively related to saprobic processes, while alleles that introgressed from the native to the invasive species were mainly linked to gene regulation. There was no overlap of allele categories introgressed in the two directions. Phenotypic experiments documented a fitness increase in H. annosum genotypes characterized by introgression of alleles from the invasive species, supporting the hypothesis that hybridization results in putatively adaptive introgression. Conversely, introgression from the native into the exotic species appeared to be driven by selection on genes favoring genome stability. Since the introgression of specific alleles from the exotic H. irregulare into the native H. annosum increased the invasiveness of the latter species, we propose that two invasions may be co-occurring: the first one by genotypes of the exotic species, and the second one by alleles belonging to the exotic species. Given that H. irregulare represents a threat to European forests, monitoring programs need to track not only exotic genotypes in native forest stands, but also exotic alleles introgressed in native genotypes.

Transcriptional profile of oil palm pathogen, Ganoderma boninense, reveals activation of lignin degradation machinery and possible evasion of host immune response

BACKGROUND

The white-rot fungi in the genus Ganoderma interact with both living and dead angiosperm tree hosts. Two Ganoderma species, a North American taxon, G. zonatum and an Asian taxon, G. boninense, have primarily been found associated with live palm hosts. During the host plant colonization process, a massive transcriptional reorganization helps the fungus evade the host immune response and utilize plant cell wall polysaccharides.

RESULTS

A publicly available transcriptome of G. boninense - oil palm interaction was surveyed to profile transcripts that were differentially expressed in planta. Ten percent of the G. boninense transcript loci had altered expression as it colonized oil palm plants one-month post inoculation. Carbohydrate active enzymes (CAZymes), particularly those with a role in lignin degradation, and auxiliary enzymes that facilitate lignin modification, like cytochrome P450s and haloacid dehalogenases, were up-regulated in planta. Several lineage specific proteins and secreted proteins that lack known functional domains were also up-regulated in planta, but their role in the interaction could not be established. A slowdown in G. boninense respiration during the interaction can be inferred from the down-regulation of proteins involved in electron transport chain and mitochondrial biogenesis. Additionally, pathogenicity related genes and chitin degradation machinery were down-regulated during the interaction indicating G. boninense may be evading detection by the host immune system.

CONCLUSIONS

This analysis offers an overview of the dynamic processes at play in G. boninense - oil palm interaction and provides a framework to investigate biology of Ganoderma fungi across plantations and landscape.

Transcript Dynamics in Wounded and Inoculated Scots Pine

Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.

Root trait-microbial relationships across tundra plant species

Fine roots, and their functional traits, influence associated rhizosphere microorganisms via root exudation and root litter quality. However, little information is known about their relationship with rhizosphere microbial taxa and functional guilds. We investigated the relationships of 11 fine root traits of 20 sub-arctic tundra meadow plant species and soil microbial community composition, using phospholipid fatty acids (PLFAs) and high-throughput sequencing. We primarily focused on the root economics spectrum, as it provides a useful framework to examine plant strategies by integrating the co-ordination of belowground root traits along a resource acquisition-conservation trade-off axis. We found that the chemical axis of the fine root economics spectrum was positively related to fungal to bacterial ratios, but negatively to Gram-positive to Gram-negative bacterial ratios. However, this spectrum was unrelated to the relative abundance of functional guilds of soil fungi. Nevertheless, the relative abundance of arbuscular mycorrhizal fungi was positively correlated to root carbon content, but negatively to the numbers of root forks per root length. Our results suggest that the fine root economics spectrum is important for predicting broader groups of soil microorganisms (i.e. fungi and bacteria), while individual root traits may be more important for predicting soil microbial taxa and functional guilds.

Genome analysis of Plectus murrayi, a nematode from continental Antarctica

Plectus murrayi is one of the most common and locally abundant invertebrates of continental Antarctic ecosystems. Because it is readily cultured on artificial medium in the laboratory and highly tolerant to an extremely harsh environment, P. murrayi is emerging as a model organism for understanding the evolutionary origin and maintenance of adaptive responses to multiple environmental stressors, including freezing and desiccation. The de novo assembled genome of P. murrayi contains 225.741 million base pairs and a total of 14,689 predicted genes. Compared to Caenorhabditis elegans, the architectural components of P. murrayi are characterized by a lower number of protein-coding genes, fewer transposable elements, but more exons, than closely related taxa from less harsh environments. We compared the transcriptomes of lab-reared P. murrayi with wild-caught P. murrayi and found genes involved in growth and cellular processing were up-regulated in lab-cultured P. murrayi, while a few genes associated with cellular metabolism and freeze tolerance were expressed at relatively lower levels. Preliminary comparative genomic and transcriptomic analyses suggest that the observed constraints on P. murrayi genome architecture and functional gene expression, including genome decay and intron retention, may be an adaptive response to persisting in a biotically simplified, yet consistently physically harsh environment.

Hallmarks of Basidiomycete Soft- and White-Rot in Wood-Decay -Omics Data of Two Armillaria Species

Wood-decaying Basidiomycetes are among the most efficient degraders of plant cell walls, making them key players in forest ecosystems, global carbon cycle, and in bio-based industries. Recent insights from -omics data revealed a high functional diversity of wood-decay strategies, especially among the traditional white-rot and brown-rot dichotomy. We examined the mechanistic bases of wood-decay in the conifer-specialists Armillaria ostoyae and Armillaria cepistipes using transcriptomic and proteomic approaches. Armillaria spp. (Fungi, Basidiomycota) include devastating pathogens of temperate forests and saprotrophs that decay wood. They have been discussed as white-rot species, though their response to wood deviates from typical white-rotters. While we observed an upregulation of a diverse suite of plant cell wall degrading enzymes, unlike white-rotters, they possess and express an atypical wood-decay repertoire in which pectinases and expansins are enriched, whereas lignin-decaying enzymes (LDEs) are generally downregulated. This combination of wood decay genes resembles the soft-rot of Ascomycota and appears widespread among Basidiomycota that produce a superficial white rot-like decay. These observations are consistent with ancestral soft-rot decay machinery conserved across asco- and Basidiomycota, a gain of efficient lignin-degrading ability in white-rot fungi and repeated, complete, or partial losses of LDE encoding gene repertoires in brown- and secondarily soft-rot fungi.

Optimized metabarcoding with Pacific biosciences enables semi-quantitative analysis of fungal communities

Recent studies have questioned the use of high-throughput sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region to derive a semi-quantitative representation of fungal community composition. However, comprehensive studies that quantify biases occurring during PCR and sequencing of ITS amplicons are still lacking. We used artificially assembled communities consisting of 10 ITS-like fragments of varying lengths and guanine-cytosine (GC) contents to evaluate and quantify biases during PCR and sequencing with Illumina MiSeq, PacBio RS II and PacBio Sequel I technologies. Fragment length variation was the main source of bias in observed community composition relative to the template, with longer fragments generally being under-represented for all sequencing platforms. This bias was three times higher for Illumina MiSeq than for PacBio RS II and Sequel I. All 10 fragments in the artificial community were recovered when sequenced with PacBio technologies, whereas the three longest fragments (> 447 bases) were lost when sequenced with Illumina MiSeq. Fragment length bias also increased linearly with increasing number of PCR cycles but could be mitigated by optimization of the PCR setup. No significant biases related to GC content were observed. Despite lower sequencing output, PacBio sequencing was better able to reflect the community composition of the template than Illumina MiSeq sequencing.

Comparative Genomics Analyses of Lifestyle Transitions at the Origin of an Invasive Fungal Pathogen in the Genus Cryphonectria

Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

Emerging fungal pathogens are a threat to forest and agroecosystems, as well as animal and human health. How pathogens evolve from nonpathogenic ancestors is still poorly understood, making the prediction of future outbreaks challenging. Most pathogens have evolved lifestyle adaptations, which were enabled by specific changes in the gene content of the species. Hence, understanding transitions in the functions encoded by genomes gives valuable insight into the evolution of pathogenicity. Here, we studied lifestyle evolution in the genus Cryphonectria, including the prominent invasive pathogen Cryphonectria parasitica, the causal agent of chestnut blight on Castanea species. We assembled and compared the genomes of pathogenic and putatively nonpathogenic Cryphonectria species, as well as sister group pathogens in the family Cryphonectriaceae (Diaporthales, Ascomycetes), to investigate the evolution of genome size and gene content. We found a striking loss of genes associated with carbohydrate metabolism (CAZymes) in C. parasitica compared to other Cryphonectriaceae. Despite substantial CAZyme gene loss, experimental data suggest that C. parasitica has retained wood colonization abilities shared with other Cryphonectria species. Putative effectors substantially varied in number, cysteine content, and protein length among species. In contrast, secondary metabolite gene clusters show a high degree of conservation within the genus. Overall, our results underpin the recent lifestyle transition of C. parasitica toward a more pathogenic lifestyle. Our findings suggest that a CAZyme loss may have promoted pathogenicity of C. parasitica on Castanea species. Analyzing gene complements underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

IMPORTANCE Forest and agroecosystems, as well as animal and human health, are threatened by emerging pathogens. Following decimation of chestnuts in the United States, the fungal pathogen Cryphonectria parasitica colonized Europe. After establishment, the pathogen population gave rise to a highly successful lineage that spread rapidly across the continent. Core to our understanding of what makes a successful pathogen is the genetic repertoire enabling the colonization and exploitation of host species. Here, we have assembled >100 genomes across two related genera to identify key genomic determinants leading to the emergence of chestnut blight. We found subtle yet highly specific changes in the transition from saprotrophy to latent pathogenicity mostly determined by enzymes involved in carbohydrate metabolism. Large-scale genomic analyses of genes underlying key nutrition modes can facilitate the detection of species with the potential to emerge as pathogens.

De novo sequencing, assembly and functional annotation of Armillaria borealis genome

BACKGROUND

Massive forest decline has been observed almost everywhere as a result of negative anthropogenic and climatic effects, which can interact with pests, fungi and other phytopathogens and aggravate their effects. Climatic changes can weaken trees and make fungi, such as Armillaria more destructive. Armillaria borealis (Marxm. & Korhonen) is a fungus from the Physalacriaceae family (Basidiomycota) widely distributed in Eurasia, including Siberia and the Far East. Species from this genus cause the root white rot disease that weakens and often kills woody plants. However, little is known about ecological behavior and genetics of A. borealis. According to field research data, A. borealis is less pathogenic than A. ostoyae, and its aggressive behavior is quite rare. Mainly A. borealis behaves as a secondary pathogen killing trees already weakened by other factors. However, changing environment might cause unpredictable effects in fungus behavior.

RESULTS

The de novo genome assembly and annotation were performed for the A. borealis species for the first time and presented in this study. The A. borealis genome assembly contained ~ 68 Mbp and was comparable with ~ 60 and ~ 79.5 Mbp for the A. ostoyae and A. mellea genomes, respectively. The N50 for contigs equaled 50,544 bp. Functional annotation analysis revealed 21,969 protein coding genes and provided data for further comparative analysis. Repetitive sequences were also identified. The main focus for further study and comparative analysis will be on the enzymes and regulatory factors associated with pathogenicity.

CONCLUSIONS

Pathogenic fungi such as Armillaria are currently one of the main problems in forest conservation. A comprehensive study of these species and their pathogenicity is of great importance and needs good genomic resources. The assembled genome of A. borealis presented in this study is of sufficiently good quality for further detailed comparative study on the composition of enzymes in other Armillaria species. There is also a fundamental problem with the identification and classification of species of the Armillaria genus, where the study of repetitive sequences in the genomes of basidiomycetes and their comparative analysis will help us identify more accurately taxonomy of these species and reveal their evolutionary relationships.

Degradative Capacity of Two Strains of Rhodonia placenta: From Phenotype to Genotype

Brown rot fungi, such as Rhodonia placenta (previously Postia placenta), occur naturally in northern coniferous forest ecosystems and are known to be the most destructive group of decay fungi, degrading wood faster and more effectively than other wood-degrading organisms. It has been shown that brown rot fungi not only rely on enzymatic degradation of lignocellulose, but also use low molecular weight oxidative agents in a non-enzymatic degradation step prior to the enzymatic degradation. R. placenta is used in standardized decay tests in both Europe and North America. However, two different strains are employed (FPRL280 and MAD-698, respectively) for which differences in colonization-rate, mass loss, as well as in gene expression have been observed, limiting the comparability of results. To elucidate the divergence between both strains, we investigated the phenotypes in more detail and compared their genomes. Significant phenotypic differences were found between the two strains, and no fusion was possible. MAD-698 degraded scots pine more aggressively, had a more constant growth rate and produced mycelia faster than FPRL280. After sequencing the genome of FPRL280 and comparing it with the published MAD-698 genome we found 660,566 SNPs, resulting in 98.4% genome identity. Specific analysis of the carbohydrate-active enzymes, encoded by the genome (CAZome) identified differences in many families related to plant biomass degradation, including SNPs, indels, gaps or insertions within structural domains. Four genes belonging to the AA3_2 family could not be found in or amplified from FPRL280 gDNA, suggesting the absence of these genes. Differences in other CAZy encoding genes that could potentially affect the lignocellulolytic activity of the strains were also predicted by comparison of genome assemblies (e.g., GH2, GH3, GH5, GH10, GH16, GH78, GT2, GT15, and CBM13). Overall, these mutations help to explain the phenotypic differences observed between both strains as they could interfere with the enzymatic activities, substrate binding ability or protein folding. The investigation of the molecular reasons that make these two strains distinct contributes to the understanding of the development of this important brown rot reference species and will help to put the data obtained from standardized decay tests across the globe into a better biological context.

Using the CODIT model to explain secondary metabolites of xylem in defence systems of temperate trees against decay fungi

BACKGROUND

In trees, secondary metabolites (SMs) are essential for determining the effectiveness of defence systems against fungi and why defences are sometimes breached. Using the CODIT model (Compartmentalization of Damage/Dysfunction in Trees), we explain defence processes at the cellular level. CODIT is a highly compartmented defence system that relies on the signalling, synthesis and transport of defence compounds through a three-dimensional lattice of parenchyma against the spread of decay fungi in xylem.

SCOPE

The model conceptualizes 'walls' that are pre-formed, formed during and formed after wounding events. For sapwood, SMs range in molecular size, which directly affects performance and the response times in which they can be produced. When triggered, high-molecular weight SMs such as suberin and lignin are synthesized slowly (phytoalexins), but can also be in place at the time of wounding (phytoanticipins). In contrast, low-molecular weight phenolic compounds such as flavonoids can be manufactured de novo (phytoalexins) rapidly in response to fungal colonization. De novo production of SMs can be regulated in response to fungal pathogenicity levels. The protective nature of heartwood is partly based on the level of accumulated antimicrobial SMs (phytoanticipins) during the transitionary stage into a normally dead substance. Effectiveness against fungal colonization in heartwood is largely determined by the genetics of the host.

CONCLUSION

Here we review recent advances in our understanding of the role of SMs in trees in the context of CODIT, with emphasis on the relationship between defence, carbohydrate availability and the hydraulic system.We also raise the limitations of the CODIT model and suggest its modification, encompassing other defence theory concepts. We envisage the development of a new defence system that is modular based and incorporates all components (and organs) of the tree from micro- to macro-scales.

The conifer root rot pathogens Heterobasidion irregulare and Heterobasidion occidentale employ different strategies to infect Norway spruce

Heterobasidion irregulare and H. occidentale are two closely related conifer root rot pathogens in the H. annosum sensu lato (s.l.) species complex. The two species H. irregulare and H. occidentale have different host preference with pine and non-pine tree species favored, respectively. The comparison of transcriptomes of H. irregulare and H. occidentale growing in Norway spruce bark, a susceptible host non-native to North America, showed large differences in gene expression. Heterobasidion irregulare induced more genes involved in detoxification of host compounds and in production of secondary metabolites, while the transcriptome induced in H. occidentale was more oriented towards carbohydrate degradation. Along with their separated evolutionary history, the difference might be driven by their host preferences as indicated by the differentially expressed genes enriched in particular Gene Ontology terms.

Fast and reliable molecular methods to detect fungal pathogens in woody plants

Plant diseases caused by pathogenic microorganisms represent a serious threat to plant productivity, food security, and natural ecosystems. An effective framework for early warning and rapid response is a crucial element to mitigate or prevent the impacts of biological invasions of plant pathogens. For these reasons, detection tools play an important role in monitoring plant health, surveillance, and quantitative pathogen risk assessment, thus improving best practices to mitigate and prevent microbial threats. The need to reduce the time of diagnosis has prompted plant pathologists to move towards more sensitive and rapid methods such as molecular techniques. Considering prevention to be the best strategy to protect plants from diseases, this review focuses on fast and reliable molecular methods to detect the presence of woody plant pathogens at early stage of disease development before symptoms occur in the host. A harmonized pool of novel technical, methodological, and conceptual solutions is needed to prevent entry and establishment of new diseases in a country and mitigate the impact of both invasive and indigenous organisms to agricultural and forest ecosystem biodiversity and productivity.

Genomic biosurveillance of forest invasive alien enemies: A story written in code

The world's forests face unprecedented threats from invasive insects and pathogens that can cause large irreversible damage to the ecosystems. This threatens the world's capacity to provide long-term fiber supply and ecosystem services that range from carbon storage, nutrient cycling, and water and air purification, to soil preservation and maintenance of wildlife habitat. Reducing the threat of forest invasive alien species requires vigilant biosurveillance, the process of gathering, integrating, interpreting, and communicating essential information about pest and pathogen threats to achieve early detection and warning and to enable better decision-making. This process is challenging due to the diversity of invasive pests and pathogens that need to be identified, the diverse pathways of introduction, and the difficulty in assessing the risk of establishment. Genomics can provide powerful new solutions to biosurveillance. The process of invasion is a story written in four chapters: transport, introduction, establishment, and spread. The series of processes that lead to a successful invasion can leave behind a DNA signature that tells the story of an invasion. This signature can help us understand the dynamic, multistep process of invasion and inform management of current and future introductions. This review describes current and future application of genomic tools and pipelines that will provide accurate identification of pests and pathogens, assign outbreak or survey samples to putative sources to identify pathways of spread, and assess risk based on traits that impact the outbreak outcome.

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.

Insight into plant cell wall degradation and pathogenesis of Ganoderma boninense via comparative genome analysis

Background

G. boninense is a hemibiotrophic fungus that infects oil palms (Elaeis guineensis Jacq.) causing basal stem rot (BSR) disease and consequent massive economic losses to the oil palm industry. The pathogenicity of this white-rot fungus has been associated with cell wall degrading enzymes (CWDEs) released during saprophytic and necrotrophic stage of infection of the oil palm host. However, there is a lack of information available on the essentiality of CWDEs in wood-decaying process and pathogenesis of this oil palm pathogen especially at molecular and genome levels.

Methods

In this study, comparative genome analysis was carried out using the G. boninense NJ3 genome to identify and characterize carbohydrate-active enzyme (CAZymes) including CWDE in the fungal genome. Augustus pipeline was employed for gene identification in G. boninense NJ3 and the produced protein sequences were analyzed via dbCAN pipeline and PhiBase 4.5 database annotation for CAZymes and plant-host interaction (PHI) gene analysis, respectively. Comparison of CAZymes from G. boninense NJ3 was made against G. lucidum, a well-studied model Ganoderma sp. and five selected pathogenic fungi for CAZymes characterization. Functional annotation of PHI genes was carried out using Web Gene Ontology Annotation Plot (WEGO) and was used for selecting candidate PHI genes related to cell wall degradation of G. boninense NJ3.

Results

G. boninense was enriched with CAZymes and CWDEs in a similar fashion to G. lucidum that corroborate with the lignocellulolytic abilities of both closely-related fungal strains. The role of polysaccharide and cell wall degrading enzymes in the hemibiotrophic mode of infection of G. boninense was investigated by analyzing the fungal CAZymes with necrotrophic Armillaria solidipes, A. mellea, biotrophic Ustilago maydis, Melampsora larici-populina and hemibiotrophic Moniliophthora perniciosa. Profiles of the selected pathogenic fungi demonstrated that necrotizing pathogens including G. boninense NJ3 exhibited an extensive set of CAZymes as compared to the more CAZymes-limited biotrophic pathogens. Following PHI analysis, several candidate genes including polygalacturonase, endo β-1,3-xylanase, β-glucanase and laccase were identified as potential CWDEs that contribute to the plant host interaction and pathogenesis.

Discussion

This study employed bioinformatics tools for providing a greater understanding of the biological mechanisms underlying the production of CAZymes in G. boninense NJ3. Identification and profiling of the fungal polysaccharide- and lignocellulosic-degrading enzymes would further facilitate in elucidating the infection mechanisms through the production of CWDEs by G. boninense. Identification of CAZymes and CWDE-related PHI genes in G. boninense would serve as the basis for functional studies of genes associated with the fungal virulence and pathogenicity using systems biology and genetic engineering approaches.

Notes, outline and divergence times of Basidiomycota

The Basidiomycota constitutes a major phylum of the kingdom Fungi and is second in species numbers to the Ascomycota. The present work provides an overview of all validly published, currently used basidiomycete genera to date in a single document. An outline of all genera of Basidiomycota is provided, which includes 1928 currently used genera names, with 1263 synonyms, which are distributed in 241 families, 68 orders, 18 classes and four subphyla. We provide brief notes for each accepted genus including information on classification, number of accepted species, type species, life mode, habitat, distribution, and sequence information. Furthermore, three phylogenetic analyses with combined LSU, SSU, 5.8s, rpb1, rpb2, and ef1 datasets for the subphyla Agaricomycotina, Pucciniomycotina and Ustilaginomycotina are conducted, respectively. Divergence time estimates are provided to the family level with 632 species from 62 orders, 168 families and 605 genera. Our study indicates that the divergence times of the subphyla in Basidiomycota are 406–430 Mya, classes are 211–383 Mya, and orders are 99–323 Mya, which are largely consistent with previous studies. In this study, all phylogenetically supported families were dated, with the families of Agaricomycotina diverging from 27–178 Mya, Pucciniomycotina from 85–222 Mya, and Ustilaginomycotina from 79–177 Mya. Divergence times as additional criterion in ranking provide additional evidence to resolve taxonomic problems in the Basidiomycota taxonomic system, and also provide a better understanding of their phylogeny and evolution.

Fungal functional ecology: bringing a trait-based approach to plant-associated fungi

Fungi play many essential roles in ecosystems. They facilitate plant access to nutrients and water, serve as decay agents that cycle carbon and nutrients through the soil, water and atmosphere, and are major regulators of macro-organismal populations. Although technological advances are improving the detection and identification of fungi, there still exist key gaps in our ecological knowledge of this kingdom, especially related to function. Trait-based approaches have been instrumental in strengthening our understanding of plant functional ecology and, as such, provide excellent models for deepening our understanding of fungal functional ecology in ways that complement insights gained from traditional and -omics-based techniques. In this review, we synthesize current knowledge of fungal functional ecology, taxonomy and systematics and introduce a novel database of fungal functional traits (Fun^Fun ). Fun^Fun is built to interface with other databases to explore and predict how fungal functional diversity varies by taxonomy, guild, and other evolutionary or ecological grouping variables. To highlight how a quantitative trait-based approach can provide new insights, we describe multiple targeted examples and end by suggesting next steps in the rapidly growing field of fungal functional ecology.

Comparative genomics reveals unique wood-decay strategies and fruiting body development in the Schizophyllaceae

Agaricomycetes are fruiting body-forming fungi that produce some of the most efficient enzyme systems to degrade wood. Despite decades-long interest in their biology, the evolution and functional diversity of both wood-decay and fruiting body formation are incompletely known. We performed comparative genomic and transcriptomic analyses of wood-decay and fruiting body development in Auriculariopsis ampla and Schizophyllum commune (Schizophyllaceae), species with secondarily simplified morphologies, an enigmatic wood-decay strategy and weak pathogenicity to woody plants. The plant cell wall-degrading enzyme repertoires of Schizophyllaceae are transitional between those of white rot species and less efficient wood-degraders such as brown rot or mycorrhizal fungi. Rich repertoires of suberinase and tannase genes were found in both species, with tannases restricted to Agaricomycetes that preferentially colonize bark-covered wood, suggesting potential complementation of their weaker wood-decaying abilities and adaptations to wood colonization through the bark. Fruiting body transcriptomes revealed a high rate of divergence in developmental gene expression, but also several genes with conserved expression patterns, including novel transcription factors and small-secreted proteins, some of the latter which might represent fruiting body effectors. Taken together, our analyses highlighted novel aspects of wood-decay and fruiting body development in an important family of mushroom-forming fungi.

Screening of the High-Rhizosphere Competent Limoniastrum monopetalum' Culturable Endophyte Microbiota Allows the Recovery of Multifaceted and Versatile Biocontrol Agents

Halophyte Limoniastrum monopetalum, an evergreen shrub inhabiting the Mediterranean region, has well-documented phytoremediation potential for metal removal from polluted sites. It is also considered to be a medicinal halophyte with potent activity against plant pathogens. Therefore, L. monopetalum may be a suitable candidate for isolating endophytic microbiota members that provide plant growth promotion (PGP) and resistance to abiotic stresses. Selected for biocontrol abilities, these endophytes may represent multifaceted and versatile biocontrol agents, combining pathogen biocontrol in addition to PGP and plant protection against abiotic stresses. In this study 117 root culturable bacterial endophytes, including Gram-positive (Bacillus and Brevibacillus), Gram-negative (Proteus, Providencia, Serratia, Pantoea, Klebsiella, Enterobacter and Pectobacterium) and actinomycete Nocardiopsis genera have been recovered from L. monopetalum. The collection exhibited high levels of biocontrol abilities against bacterial (Agrobacterium tumefaciens MAT2 and Pectobacterium carotovorum MAT3) and fungal (Alternaria alternata XSZJY-1, Rhizoctonia bataticola MAT1 and Fusarium oxysporum f. sp. radicis lycopersici FORL) pathogens. Several bacteria also showed PGP capacity and resistance to antibiotics and metals. A highly promising candidate Bacillus licheniformis LMRE 36 with high PGP, biocontrol, metal and antibiotic, resistance was subsequently tested in planta (potato and olive trees) for biocontrol of a collection of 14 highly damaging Fusarium species. LMRE 36 proved very effective against the collection in both species and against an emerging Fusarium sp. threatening olive trees culture in nurseries. These findings provide a demonstration of our pyramiding strategy. Our strategy was effective in combining desirable traits in biocontrol agents towards broad-spectrum resistance against pathogens and protection of crops from abiotic stresses. Stacking multiple desirable traits into a single biocontrol agent is achieved by first, careful selection of a host for endophytic microbiota recovery; second, stringent in vitro selection of candidates from the collection; and third, application of the selected biocontrol agents in planta experiments. That pyramiding strategy could be successfully used to mitigate effects of diverse biotic and abiotic stresses on plant growth and productivity. It is anticipated that the strategy will provide a new generation of biocontrol agents by targeting the microbiota of plants in hostile environments.

Fungal evolution: major ecological adaptations and evolutionary transitions

Fungi are a highly diverse group of heterotrophic eukaryotes characterized by the absence of phagotrophy and the presence of a chitinous cell wall. While unicellular fungi are far from rare, part of the evolutionary success of the group resides in their ability to grow indefinitely as a cylindrical multinucleated cell (hypha). Armed with these morphological traits and with an extremely high metabolical diversity, fungi have conquered numerous ecological niches and have shaped a whole world of interactions with other living organisms. Herein we survey the main evolutionary and ecological processes that have guided fungal diversity. We will first review the ecology and evolution of the zoosporic lineages and the process of terrestrialization, as one of the major evolutionary transitions in this kingdom. Several plausible scenarios have been proposed for fungal terrestralization and we here propose a new scenario, which considers icy environments as a transitory niche between water and emerged land. We then focus on exploring the main ecological relationships of Fungi with other organisms (other fungi, protozoans, animals and plants), as well as the origin of adaptations to certain specialized ecological niches within the group (lichens, black fungi and yeasts). Throughout this review we use an evolutionary and comparative-genomics perspective to understand fungal ecological diversity. Finally, we highlight the importance of genome-enabled inferences to envision plausible narratives and scenarios for important transitions.

Global population genomics of the forest pathogen Dothistroma septosporum reveal chromosome duplications in high dothistromin-producing strains

Dothistroma needle blight is one of the most devastating pine tree diseases worldwide. New and emerging epidemics have been frequent over the last 25 years, particularly in the Northern Hemisphere, where they are in part associated with changing weather patterns. One of the main Dothistroma needle blight pathogens, Dothistroma septosporum, has a global distribution but most molecular plant pathology research has been confined to Southern Hemisphere populations that have limited genetic diversity. Extensive genomic and transcriptomic data are available for a D. septosporum reference strain from New Zealand, where an introduced clonal population of the pathogen predominates. Due to the global importance of this pathogen, we determined whether the genome of this reference strain is representative of the species worldwide by sequencing the genomes of 18 strains sampled globally from different pine hosts. Genomic polymorphism shows substantial variation within the species, clustered into two distinct groups of strains with centres of diversity in Central and South America. A reciprocal chromosome translocation uniquely identifies the New Zealand strains. Globally, strains differ in their production of the virulence factor dothistromin, with extremely high production levels in strain ALP3 from Germany. Comparisons with the New Zealand reference revealed that several strains are aneuploids; for example, ALP3 has duplications of three chromosomes. Increased gene copy numbers therefore appear to contribute to increased production of dothistromin, emphasizing that studies of population structure are a necessary adjunct to functional analyses of genetic polymorphisms to identify the molecular basis of virulence in this important forest pathogen.