Genetic isolation between two recently diverged populations of a symbiotic fungus

How do phytophagous insects affect phyllosphere fungi? Tracking fungi from milkweed to monarch caterpillar frass reveals communities dominated by fungal yeast

Since a significant proportion of plant matter is consumed by herbivores, a necessary adaptation for many phyllosphere microbes could be to survive through the guts of herbivores. While many studies explore the gut microbiome of herbivores by surveying the microbiome in their frass, few studies compare the phyllosphere microbiome to the gut microbiome of herbivores. High-throughput metabarcode sequencing was used to track the fungal community from milkweed (Asclepias spp.) leaves to monarch caterpillar frass. The most commonly identified fungal taxa that dominated the caterpillar frass after the consumption of leaves were yeasts, mostly belonging to the Basidiomycota phylum. While most fungal communities underwent significant bottlenecks and some yeast taxa increased in relative abundance, a consistent directional change in community structure was not identified from leaf to caterpillar frass. These results suggest that some phyllosphere fungi, especially diverse yeasts, can survive herbivory, but whether herbivory is a key stage of their life cycle remains uncertain. For exploring phyllosphere fungi and the potential coprophilous lifestyles of endophytic and epiphytic fungi, methods that target yeast and Basidiomycota fungi are recommended.

Suillus: an emerging model for the study of ectomycorrhizal ecology and evolution

Research on mycorrhizal symbiosis has been slowed by a lack of established study systems. To address this challenge, we have been developing Suillus, a widespread ecologically and economically relevant fungal genus primarily associated with the plant family Pinaceae, into a model system for studying ectomycorrhizal (ECM) associations. Over the last decade, we have compiled extensive genomic resources, culture libraries, a phenotype database, and protocols for manipulating Suillus fungi with and without their tree partners. Our efforts have already resulted in a large number of publicly available genomes, transcriptomes, and respective annotations, as well as advances in our understanding of mycorrhizal partner specificity and host communication, fungal and plant nutrition, environmental adaptation, soil nutrient cycling, interspecific competition, and biological invasions. Here, we highlight the most significant recent findings enabled by Suillus, present a suite of protocols for working with the genus, and discuss how Suillus is emerging as an important model to elucidate the ecology and evolution of ECM interactions.

Arbuscular mycorrhizal fungi-induced tolerance to chromium stress in plants

Chromium (Cr) is one of the toxic elements that harms all forms of life, including plants. Industrial discharges and mining largely contribute to Cr release into the soil environment. Excessive Cr pollution in arable land significantly reduces the yield and quality of important agricultural crops. Therefore, remediation of polluted soil is imperative not only for agricultural sustainability but also for food safety. Arbuscular mycorrhizal fungi (AMF) are widespread soil-borne endophytic fungi that form mutualistic relationships with the vast majority of land plants. In mycorrhizal symbiosis, AMF are largely dependent on the host plant-supplied carbohydrates and lipids, in return, AMF aid the host plants in acquiring water and mineral nutrients, especially phosphorus, nitrogen and sulfur from distant soils, and this distinguishing feature of the two-way exchange of resources is a functional requirement for such mutualism and ecosystem services. In addition to supplying nutrients and water to plants, the AMF symbiosis enhances plant resilience to biotic and abiotic stresses including Cr stress. Studies have revealed vital physiological and molecular mechanisms by which AMF alleviate Cr phytotoxicity and aid plants in nutrient acquisition under Cr stress. Notably, plant Cr tolerance is enhanced by both the direct effects of AMF on Cr stabilization and transformation, and the indirect effects of AMF symbiosis on plant nutrient uptake and physiological regulation. In this article, we summarized the research progress on AMF and associated mechanisms of Cr tolerance in plants. In addition, we reviewed the present understanding of AMF-assisted Cr remediation. Since AMF symbiosis can enhance plant resilience to Cr pollution, AMF may have promising prospects in agricultural production, bioremediation, and ecological restoration in Cr-polluted soils.

Population genomic analysis reveals geographic structure and climatic diversification for Macrophomina phaseolina isolated from soybean and dry bean across the United States, Puerto Rico, and Colombia

Macrophomina phaseolina causes charcoal rot, which can significantly reduce yield and seed quality of soybean and dry bean resulting from primarily environmental stressors. Although charcoal rot has been recognized as a warm climate-driven disease of increasing concern under global climate change, knowledge regarding population genetics and climatic variables contributing to the genetic diversity of M. phaseolina is limited. This study conducted genome sequencing for 95 M. phaseolina isolates from soybean and dry bean across the continental United States, Puerto Rico, and Colombia. Inference on the population structure using 76,981 single nucleotide polymorphisms (SNPs) revealed that the isolates exhibited a discrete genetic clustering at the continental level and a continuous genetic differentiation regionally. A majority of isolates from the United States (96%) grouped in a clade with a predominantly clonal genetic structure, while 88% of Puerto Rican and Colombian isolates from dry bean were assigned to a separate clade with higher genetic diversity. A redundancy analysis (RDA) was used to estimate the contributions of climate and spatial structure to genomic variation (11,421 unlinked SNPs). Climate significantly contributed to genomic variation at a continental level with temperature seasonality explaining the most variation while precipitation of warmest quarter explaining the most when spatial structure was accounted for. The loci significantly associated with multivariate climate were found closely to the genes related to fungal stress responses, including transmembrane transport, glycoside hydrolase activity and a heat-shock protein, which may mediate climatic adaptation for M. phaseolina. On the contrary, limited genome-wide differentiation among populations by hosts was observed. These findings highlight the importance of population genetics and identify candidate genes of M. phaseolina that can be used to elucidate the molecular mechanisms that underly climatic adaptation to the changing climate.

Genomic determination of breeding systems and trans-specific evolution of HD MAT genes in suilloid fungi

Studying the signatures of evolution can help to understand genetic processes. Here, we demonstrate how the existence of balancing selection can be used to identify the breeding systems of fungi from genomic data. The breeding systems of fungi are controlled by self-incompatibility loci that determine mating types between potential mating partners, resulting in strong balancing selection at the loci. Within the fungal phylum Basidiomycota, two such self-incompatibility loci, namely HD MAT locus and P/R MAT locus, control mating types of gametes. Loss of function at one or both MAT loci results in different breeding systems and relaxes the MAT locus from balancing selection. By investigating the signatures of balancing selection at MAT loci, one can infer a species' breeding system without culture-based studies. Nevertheless, the extreme sequence divergence among MAT alleles imposes challenges for retrieving full variants from both alleles when using the conventional read-mapping method. Therefore, we employed a combination of read-mapping and local de novo assembly to construct haplotypes of HD MAT alleles from genomes in suilloid fungi (genera Suillus and Rhizopogon). Genealogy and pairwise divergence of HD MAT alleles showed that the origins of mating types predate the split between these two closely related genera. High sequence divergence, trans-specific polymorphism, and the deeply diverging genealogy confirm the long-term functionality and multiallelic status of HD MAT locus in suilloid fungi. This work highlights a genomics approach to studying breeding systems regardless of the culturability of organisms based on the interplay between evolution and genetics.

Advancing research on ectomycorrhizal fungal adaptation with landscape genomics

Ectomycorrhizal (ECM) fungi serve key functions in forest ecosystems by supplying water and nutrients to tree hosts, yet mutualistic plant-fungi interactions are jeopardised by environmental alterations. Here, we discuss the great potential and current limitations of landscape genomics in investigating signatures of local adaptation in natural populations of ECM fungi.

Local adaptation in fungi

In this review, I explore the pervasive but underappreciated role of local adaptation in fungi. It has been difficult historically to study local adaptation in fungi because of the limited understanding of fungal species and their traits, but new hope has been offered with technological advances in sequencing. The filamentous nature of fungi invalidates some assumptions made in evolution because of their ability to exist as multinucleate entities with genetically different nuclei sharing the same cytoplasm. Many insights on local adaptation have come from studying fungi, and much of the empirical evidence gathered about local adaptation in the context of host-pathogen interactions comes from studying fungal virulence genes, drug resistance, and environmental adaptation. Together, these insights paint a picture of the variety of processes involved in fungal local adaptation and their connections to the unusual cell biology of Fungi (multinucleate, filamentous habit), but there is much that remains unknown, with major gaps in our knowledge of fungal species, their phenotypes, and the ways by which they adapt to local conditions.

Mechanisms of stress tolerance and their effects on the ecology and evolution of mycorrhizal fungi

Stress is ubiquitous and disrupts homeostasis, leading to damage, decreased fitness, and even death. Like other organisms, mycorrhizal fungi evolved mechanisms for stress tolerance that allow them to persist or even thrive under environmental stress. Such mechanisms can also protect their obligate plant partners, contributing to their health and survival under hostile conditions. Here we review the effects of stress and mechanisms of stress response in mycorrhizal fungi. We cover molecular and cellular aspects of stress and how stress impacts individual fitness, physiology, growth, reproduction, and interactions with plant partners, along with how some fungi evolved to tolerate hostile environmental conditions. We also address how stress and stress tolerance can lead to adaptation and have cascading effects on population- and community-level diversity. We argue that mycorrhizal fungal stress tolerance can strongly shape not only fungal and plant physiology, but also their ecology and evolution. We conclude by pointing out knowledge gaps and important future research directions required for both fully understanding stress tolerance in the mycorrhizal context and addressing ongoing environmental change.

Physiological and transcriptional responses of the ectomycorrhizal fungus Cenococcum geophilum to salt stress

Ectomycorrhizal (ECM) fungi improve the host plant's tolerance to abiotic and biotic stresses. Cenococcum geophilum (Cg) is among the most common ECM fungi worldwide and often grows in saline environments. However, the physiological and molecular mechanisms of salt tolerance in this fungus are largely unknown. In the present study, 12 isolates collected from different ecogeographic regions were used to investigate the mechanism of salt tolerance of Cg. The isolates were classified into four groups (salt-sensitive, moderately salt-tolerant, salt-tolerant, and halophilic) based on their in vitro mycelial growth under 0, 50, 125, 250, and 500 mM NaCl concentrations. Hence, the Na, Ca, P, and K concentrations of mycelia and the pH of the culture solution were determined. Compared with salt-tolerant isolates, treatment with 250 mM NaCl significantly increased the sodium concentration and decreased the potassium concentration of salt-sensitive isolates. RNA-sequencing and qRT-PCR analysis were conducted to identify differentially expressed genes (DEGs) involved in transmembrane transport and oxidoreductase activity pathways. The hydrogen peroxide concentration and activities of peroxidase and superoxide dismutase in mycelia were determined, and the accumulation and scavenging of reactive oxygen species in the salt-sensitive isolates were more active than those in the salt-tolerant isolates. The results supply functional validations to RNA-seq and qRT-PCR analysis. This study provides novel insights into the salt-stress response of Cg isolates and provides a foundation for elucidation of the salt-tolerance mechanism of ECM fungi.

Effect of sodium chloride on basidiospore germination and vegetative mycelial growth of the ectomycorrhizal fungus <i>Rhizopogon roseolus</i>

Rhizopogon roseolus is a basidiomycetous ectomycorrhizal fungus that inhabits mainly coastal areas. Understanding the response of this fungus to salinity at each stage of its life cycle will lead to elucidation of the strategies for its propagation. This study examined the effect of sodium chloride (NaCl) on basidiospore germination and mycelial growth of both homokaryotic and heterokaryotic strains of R. roseolus, on nutrient agar media with varying concentrations of NaCl (0, 50, 150, and 300 mM). Regardless of the presence of NaCl, R. roseolus basidiospores germinated and the germlings grew, forming compatible fusions. In addition, all multispore strains, including homokaryons and heterokaryons, grew under these NaCl conditions. Most of these strains had an effective concentration inhibiting mycelial growth by 50% value greater than 300 mM of NaCl. These results indicate that R. roseolus can germinate, grow, and mate in the presence of NaCl, allowing it to propagate in saline habitats.

Population genomics of a forest fungus reveals high gene flow and climate adaptation signatures

Genome sequencing of spatially distributed individuals sheds light on how evolution structures genetic variation. Populations of Phellopilus nigrolimitatus, a red-listed wood-inhabiting fungus associated with old-growth coniferous forests, have decreased in size over the last century due to a loss of suitable habitats. We assessed the population genetic structure and investigated local adaptation in P. nigrolimitatus, by establishing a reference genome and genotyping 327 individuals sampled from 24 locations in Northern Europe by RAD sequencing. We revealed a shallow population genetic structure, indicating large historical population sizes and high levels of gene flow. Despite this weak sub-structuring, two genetic groups were recognized; a western group distributed mostly in Norway and an eastern group covering most of Finland, Poland and Russia. This sub-structuring may reflect co-immigration with the main host, Norway spruce (Picea abies), into Northern Europe after the last ice age. We found evidence of low levels of genetic diversity in southwestern Finland, which has a long history of intensive forestry and urbanization. Numerous loci were significantly associated with one or more environmental factors, indicating adaptation to specific environments. These loci clustered into two groups with different associations with temperature and precipitation. Overall, our findings indicate that the current population genetic structure of P. nigrolimitatus results from a combination of gene flow, genetic drift and selection. The acquisition of similar knowledge especially over broad geographic scales, linking signatures of adaptive genetic variation to evolutionary processes and environmental variation, for other fungal species will undoubtedly be useful for assessment of the combined effects of habitat fragmentation and climate change on fungi strongly bound to old-growth forests.

Deep sequencing across multiple host species tests pine-endophyte specificity

PREMISE

Foliar fungal endophytes vary in their distributions across landscapes or plant host taxa, indicative of specialized ecologies and host specific adaptations. Accounts of specialization, however, depend on the taxonomic breadth and geographic range of the host plants included in each study. A broad region-scale study or deep sampling of diverse potential host species still remains relatively rare but is becoming increasingly possible with high-throughput sequencing.

METHODS

Amplicon sequencing was used to rapidly identify the fungal endophytic community among six pine (Pinus, Pinaceae) species co-occurring across northeastern United States and to test for site and host specialization. We focused on the endophytic genus Lophodermium (Rhytismataceae), whose species members are thought to specialize on different pine species, to test if amplicon sequencing could rapidly verify previously implied or discover new patterns of host specificity.

RESULTS

While amplicon sequencing could analyze more samples at greater depths and recover greater numbers of unique Lophodermium taxa than when endophyte communities were surveyed with traditional culturing methods, patterns of specialization were not better supported. This may be because amplicon sequencing can indiscriminately capture non-host specific organisms found incidentally from plant tissues or because we have overestimated host-specificity in the past with biased culturing techniques.

CONCLUSIONS

Amplicon sequencing can quickly identify patterns of host specificity by allowing large-scale surveys but has limitations in quantifying the level of intimacy of these relationships.

Predictors of taxonomic and functional composition of black spruce seedling ectomycorrhizal fungal communities along peatland drainage gradients

Many trees depend on symbiotic ectomycorrhizal fungi for nutrients in exchange for photosynthetically derived carbohydrates. Trees growing in peatlands, which cover 3% of the earth's terrestrial surface area yet hold approximately one-third of organic soil carbon stocks, may benefit from ectomycorrhizal fungi that can efficiently forage for nutrients and degrade organic matter using oxidative enzymes such as class II peroxidases. However, such traits may place a higher carbon cost on both the fungi and host tree. To investigate these trade-offs that might structure peatland ectomycorrhizal fungal communities, we sampled black spruce (Picea mariana (Mill.)) seedlings along 100-year-old peatland drainage gradients in Minnesota, USA, that had resulted in higher soil nitrogen and canopy density. Structural equation models revealed that the relative abundance of the dominant ectomycorrhizal fungal genus, Cortinarius, which is known for relatively high fungal biomass coupled with elevated class II peroxidase potential, was negatively linked to site fertility but more positively affected by recent host stem radial growth, suggesting carbon limitation. In contrast, Cenococcum, known for comparatively lower fungal biomass and less class II peroxidase potential, was negatively linked to host stem radial growth and unrelated to site fertility. Like Cortinarius, the estimated relative abundance of class II peroxidase genes in the ectomycorrhizal community was more related to host stem radial growth than site fertility. Our findings indicate a trade-off between symbiont foraging traits and associated carbon costs that consequently structure seedling ectomycorrhizal fungal communities in peatlands.

Population genomics provides insights into the genetic basis of adaptive evolution in the mushroom-forming fungus Lentinula edodes

•

We detected three subgroups of L. edodes with robust phenotypic differentiation.

•

The three subgroups were diverged 36,871 generations ago.

•

L. edodes cultivars in China might originate from the vicinity of Northeast China.

•

We dissected the genetic basis of adaptive evolution in L. edodes.

•

Genes related to fruiting body development are involved in adaptive evolution.

Introduction

Mushroom-forming fungi comprise diverse species that develop complex multicellular structures. In cultivated species, both ecological adaptation and artificial selection have driven genome evolution. However, little is known about the connections among genotype, phenotype and adaptation in mushroom-forming fungi.

Objectives

This study aimed to (1) uncover the population structure and demographic history of Lentinula edodes, (2) dissect the genetic basis of adaptive evolution in L. edodes, and (3) determine if genes related to fruiting body development are involved in adaptive evolution.

Methods

We analyzed genomes and fruiting body-related traits (FBRTs) in 133 L. edodes strains and conducted RNA-seq analysis of fruiting body development in the YS69 strain. Combined methods of genomic scan for divergence, genome-wide association studies (GWAS), and RNA-seq were used to dissect the genetic basis of adaptive evolution.

Results

We detected three distinct subgroups of L. edodes via single nucleotide polymorphisms, which showed robust phenotypic and temperature response differentiation and correlation with geographical distribution. Demographic history inference suggests that the subgroups diverged 36,871 generations ago. Moreover, L. edodes cultivars in China may have originated from the vicinity of Northeast China. A total of 942 genes were found to be related to genetic divergence by genomic scan, and 719 genes were identified to be candidates underlying FBRTs by GWAS. Integrating results of genomic scan and GWAS, 80 genes were detected to be related to phenotypic differentiation. A total of 364 genes related to fruiting body development were involved in genetic divergence and phenotypic differentiation.

Conclusion

Adaptation to the local environment, especially temperature, triggered genetic divergence and phenotypic differentiation of L. edodes. A general model for genetic divergence and phenotypic differentiation during adaptive evolution in L. edodes, which involves in signal perception and transduction, transcriptional regulation, and fruiting body morphogenesis, was also integrated here.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates—namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose–methanol–choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases—we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

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.

Gene Copy Number Variation Does Not Reflect Structure or Environmental Selection in Two Recently Diverged California Populations of Suillus brevipes

Gene copy number variation across individuals has been shown to track population structure and be a source of adaptive genetic variation with significant fitness impacts. In this study, we report opposite results for both predictions based on the analysis of gene copy number variants (CNVs) of Suillus brevipes, a mycorrhizal fungus adapted to coastal and montane habitats in California. In order to assess whether gene copy number variation mirrored population structure and selection in this species, we investigated two previously studied locally adapted populations showing a highly differentiated genomic region encompassing a gene predicted to confer salt tolerance. In addition, we examined whether copy number in the genes related to salt homeostasis was differentiated between the two populations. Although we found many instances of CNV regions across the genomes of S. brevipes individuals, we also found CNVs did not recover population structure and known salt-tolerance-related genes were not under selection across the coastal population. Our results contrast with predictions of CNVs matching single-nucleotide polymorphism divergence and showed CNVs of genes for salt homeostasis are not under selection in S. brevipes.

Contrasting genetic structuring in the closely related basidiomycetes Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales)

Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales, Basidiomycota) are closely related saprotrophic fungi, widely distributed on coniferous wood in temperate regions worldwide. Three intersterility groups have previously been detected in T. abietinum, while no prezygotic barriers have been proven within T. fuscoviolaceum. The aim of this study was to reveal the phylogeography and genetic relationship between these two closely related species and to explore whether the previously observed intersterility groups in T. abietinum are reflected in the genetic data. We assembled worldwide fruit body collections of both species (N = 314) and generated DNA sequences from three nuclear (ITS2, LSU, IGS) and one mitochondrial rDNA region (mtLSU). The two species are genetically well separated in all analyses. In correspondence with observations from earlier mating studies, our results revealed that T. fuscoviolaceum is genetically more uniform than T. abietinum. Multiple genetic sub-groups exist in T. abietinum that may correspond to the previously observed intersterility groups. However, there is low consistency across the investigated loci in delimiting the different sub-groups, except for a consistent North American group. As for many other widespread fungi, a complex phylogeographic pattern is found in T. abietinum which may have been formed by geographic, as well as multiple genetic intersterility barriers.

Intra-species genetic variability drives carbon metabolism and symbiotic host interactions in the ectomycorrhizal fungus Pisolithus microcarpus

Ectomycorrhizal (ECM) fungi are integral to boreal and temperate forest ecosystem functioning and nutrient cycling. ECM fungi, however, originate from diverse saprotrophic lineages and the impacts of genetic variation across species, and especially within a given ECM species, on function and interactions with the environment is not well understood. Here, we explore the extent of intra-species variation between four isolates of the ECM fungus Pisolithus microcarpus, in terms of gene regulation, carbon metabolism and growth, and interactions with a host, Eucalyptus grandis. We demonstrate that, while a core response to the host is maintained by all of the isolates tested, they have distinct patterns of gene expression and carbon metabolism, resulting in the differential expression of isolate-specific response pathways in the host plant. Together, these results highlight the importance of using a wider range of individuals within a species to understand the broader ecological roles of ECM fungi and their host interactions.

High levels of endemism and local differentiation in the fungal and algal symbionts of saxicolous lecideoid lichens along a latitudinal gradient in southern South America

Saxicolous, lecideoid lichenized fungi have a cosmopolitan distribution but, being mostly cold adapted, are especially abundant in polar and high-mountain regions. To date, little is known of their origin or the extent of their trans-equatorial dispersal. Several mycobiont genera and species are thought to be restricted to either the Northern or the Southern Hemisphere, whereas others are thought to be widely distributed and occur in both hemispheres. However, these assumptions often rely on morphological analyses and lack supporting molecular genetic data. Also unknown is the extent of regional differentiation in the southern polar regions. An extensive set of lecideoid lichens (185 samples) was collected along a latitudinal gradient at the southern end of South America. Subantarctic climate conditions were maintained by increasing the elevation of the collecting sites with decreasing latitude. The investigated specimens were placed in a global context by including Antarctic and cosmopolitan sequences from other studies. For each symbiont three markers were used to identify intraspecific variation (mycobiont: ITS, mtSSU, RPB1; photobiont: ITS, psbJ-L, COX2). For the mycobiont, the saxicolous genera Lecidea, Porpidia, Poeltidea and Lecidella were phylogenetically re-evaluated, along with their photobionts Asterochloris and Trebouxia. For several globally distributed species groups, the results show geographically highly differentiated subclades, classified as operational taxonomical units (OTUs), which were assigned to the different regions of southern South America (sSA). Furthermore, several small endemic and well-supported clades apparently restricted to sSA were detected at the species level for both symbionts.

A meta-analysis of global fungal distribution reveals climate-driven patterns

The evolutionary and environmental factors that shape fungal biogeography are incompletely understood. Here, we assemble a large dataset consisting of previously generated mycobiome data linked to specific geographical locations across the world. We use this dataset to describe the distribution of fungal taxa and to look for correlations with different environmental factors such as climate, soil and vegetation variables. Our meta-study identifies climate as an important driver of different aspects of fungal biogeography, including the global distribution of common fungi as well as the composition and diversity of fungal communities. In our analysis, fungal diversity is concentrated at high latitudes, in contrast with the opposite pattern previously shown for plants and other organisms. Mycorrhizal fungi appear to have narrower climatic tolerances than pathogenic fungi. We speculate that climate change could affect ecosystem functioning because of the narrow climatic tolerances of key fungal taxa.

The authors assemble and analyse previously generated mycobiome data linked to geographical locations across the world. They describe the distribution of fungal taxa and show that climate is an important driver of fungal biogeography and that fungal diversity appears to be concentrated at high latitudes.

Maintenance of Sympatric and Allopatric Populations in Free-Living Terrestrial Bacteria

Due to the promiscuous exchange of genetic material and asexual reproduction, delineating microbial species (and, by extension, populations) remains challenging. Because of this, the vast majority of microbial studies assessing population structure often compare divergent strains from disparate environments under varied selective pressures. Here, we investigated the population structure within a single bacterial ecotype, a unit equivalent to a eukaryotic species, defined as highly clustered genotypic and phenotypic strains with the same ecological niche. Using a combination of genomic and computational analyses, we assessed the phylogenetic structure, extent of recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. To our knowledge, this study is the first to do so for a dominant soil bacterium. Our results indicate that bacterial soil populations, similarly to those in other environments, are structured by gene flow discontinuities and exhibit distributional patterns consistent with both isolation by distance and isolation by environment. Thus, both dispersal limitation and local environments contribute to the divergence among closely related soil bacteria as observed in macroorganisms.

For free-living bacteria and archaea, the equivalent of the biological species concept does not exist, creating several obstacles to the study of the processes contributing to microbial diversification. These obstacles are particularly high in soil, where high bacterial diversity inhibits the study of closely related genotypes and therefore the factors structuring microbial populations. Here, we isolated strains within a single Curtobacterium ecotype from surface soil (leaf litter) across a regional climate gradient and investigated the phylogenetic structure, recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. Our results indicate that microbial populations are delineated by gene flow discontinuities, with distinct populations cooccurring at multiple sites. Bacterial population structure was further delineated by genomic features allowing for the identification of candidate genes possibly contributing to local adaptation. These results suggest that the genetic structure within this bacterium is maintained both by ecological specialization in localized microenvironments (isolation by environment) and by dispersal limitation between geographic locations (isolation by distance).

Fifty Aureobasidium pullulans genomes reveal a recombining polyextremotolerant generalist

The black yeast Aureobasidium pullulans is a textbook example of a generalistic and ubiquitous fungus thriving in a wide variety of environments. To investigate whether A. pullulans is a true generalist, or alternatively, whether part of its versatility can be attributed to intraspecific specialization masked by cryptic diversification undetectable by traditional phylogenetic analyses, we sequenced and analysed the genomes of 50 strains of A. pullulans from different habitats and geographic locations. No population structure was observed in the sequenced strains. Decay of linkage disequilibrium over shorter physical distances (<100 bp) than in many sexually reproducing fungi indicates a high level of recombination in the species. A homothallic mating locus was found in all of the sequenced genomes. Aureobasidium pullulans appears to have a homogeneous population genetics structure, which is best explained by good dispersal and high levels of recombination. This means that A. pullulans is a true generalist that can inhabit different habitats without substantial specialization to any of these habitats at the genomic level. Furthermore, in the future, the high level of A. pullulans recombination can be exploited for the identification of genomic loci that are involved in the many biotechnologically useful traits of this black yeast.

Population Genomics of an Obligately Halophilic Basidiomycete Wallemia ichthyophaga

Background

Wallemia ichthyophaga is a highly specialized basidiomycetous fungus. It is one of the most halophilic fungi ever described, only able to grow at low water activity. This specialization is thought to explain why it is only rarely isolated from nature.

Results

Genomes of 21 W. ichthyophaga strains were sequenced with PE150 reads on BGISEQ500 platform. The genomes shared high similarity with the reference genome of the species, they were all smaller than 10 Mbp and had a low number of predicted genes. Groups of strains isolated in the same location encompassed clones as well as very divergent strains. There was little concordance between phylogenies of predicted genes. Linkage disequilibrium of pairs of polymorphic loci decayed relatively quickly as a function of distance between the loci (LD decay distance 1270 bp). For the first time a putative mating-type locus was identified in the genomes of W. ichthyophaga.

Conclusion

Based on the comparison of W. ichthyophaga genomes it appears that some phylogenetic lineages of the species can persist in the same location over at least several years. Apart from this, the differences between the strains do not reflect the isolation habitat or geographic location. Together with results supporting the existence of (sexual) recombination in W. ichthyophaga, the presented results indicate that strains of W. ichthyophaga can form a single recombining population even between different habitats and over large geographical distances.

Pleurotus eryngii Genomes Reveal Evolution and Adaptation to the Gobi Desert Environment

Pleurotus eryngii (King Oyster) is one of the most highly prized edible mushrooms. Among the diverse varieties within P. eryngii, P. eryngii var. eryngii is the commonest one, with a worldwide distribution, while P. eryngii var. ferulae is only distributed in Europe and China, and is especially adapted to the Gobi Desert in Xinjiang Autonomous Region of China. However, little is known about the genome-wide pattern of evolution and adaptation to the divergent environments of P. eryngii. Here, we present the high-quality genome sequences of P. eryngii var. eryngii strain PEE81 originating from Europe and P. eryngii var. ferulae strain PEF12 originating from the Gobi Desert of China. The assembled genome sizes of PEE81 and PEF12 were 53.6 and 48.0 Mbp, respectively, which are larger than other reported genomes in the genus Pleurotus. We propose that the selective amplification of long terminal repeat (LTR) retrotransposons increases the genome size of the genus Pleurotus, and may play a key role in driving their rapid species diversification. Molecular clock analyses of five Pleurotus species, namely PEE81, PEF12, P. tuoliensis, P. ostreatus and P. cf. floridanus suggest that the divergence estimates of the genus Pleurotus over time scales ranged from ∼4 to ∼38 million years ago (Mya), and PEE81 and PEF12 diverged at ∼13 Mya. The whole genome resequencing of 33 geographically diverse strains of P. eryngii var. eryngii and var. ferulae was then performed and the genome variation among and within these two populations were investigated. Comparative analyses of these two populations detected several candidate genes related to stress responses and DNA repair that are putatively involved in adaptation to the Gobi Desert environment. These findings offer insights into genome evolution of the genus Pleurotus and provide valuable genomic resources for King Oyster mushroom breeding.

Phosphorus Transport in Mycorrhiza: How Far Are We?

Mycorrhizal fungi considerably improve plant nutrition and help them to cope with changing environments. Particularly, these fungi express proteins to transfer inorganic phosphate (P_i) from the soil to colonized roots through symbiotic interfaces. The mechanisms involved in P_i transfer from fungal to plant cells are still largely unknown. Here, we discuss the recent progress made on the description of these mechanisms and we propose the most promising hypotheses and alternative mechanisms for this process. Specifically, we present a phylogenetic survey of candidate P_i transporters of mycorrhizal fungi that might ensure P_i unload into the symbiotic interfaces. Gathering additional knowledge on mycorrhizal P_i transport will improve the P_i-useefficiency in agroecological systems and will guide towards addressing future research challenges.

Ectomycorrhizal symbiosis helps plants to challenge salt stress conditions

Soil salinity is an environmental condition that is currently increasing worldwide. Plant growth under salinity induces osmotic stress and ion toxicity impairing root water and nutrient absorption, but the association with beneficial soil microorganisms has been linked to an improved adaptation to this constraint. The ectomycorrhizal (ECM) symbiosis has been proposed as a key factor for a better tolerance of woody species to salt stress, thanks to the reduction of sodium (Na⁺) uptake towards photosynthetic organs. Although no precise mechanisms for this enhanced plant salt tolerance have been described yet, in this review, we summarize the knowledge accumulated so far on the role of ECM symbiosis. Moreover, we propose several strategies by which ECM fungi might help plants, including restriction of Na⁺ entrance into plant tissues and improvement of mineral nutrition and water balances. This positive effect of ECM fungi has been proven in field assays and the results obtained point to a promising application in forestry cultures and reforestation.

Role of urban ectomycorrhizal fungi in improving the tolerance of lodgepole pine (Pinus contorta) seedlings to salt stress

With large forested urban areas, the city of Edmonton, Alberta, Canada, faces high annual costs of replacing trees injured by deicing salts that are commonly used for winter road maintenance. Ectomycorrhizal fungi form symbiotic associations with tree roots that allow trees to tolerate the detrimental effects of polluted soils. Here, we examined mycorrhizal colonization of Pinus contorta by germinating seeds in soils collected from different locations: (1) two urban areas within the city of Edmonton, and (2) an intact pine forest just outside Edmonton. We then tested the responses of seedlings to 0-, 60-, and 90-mM NaCl. Our results showed lower abundance and diversity of ectomycorrhizal fungi in seedlings colonized with the urban soils compared to those from the pine forest soil. However, when subsequently exposed to NaCl treatments, only seedlings inoculated with one of the urban soils containing fungi from the genera Tuber, Suillus, and Wilcoxina, showed reduced shoot Na accumulation and higher growth rates. Our results indicate that local ectomycorrhizal fungi that are adapted to challenging urban sites may offer a potential suitable source for inoculum for conifer trees designated for plating in polluted urban environments.

Genomic Evidence of Recombination in the Basidiomycete Wallemia mellicola

One of the most commonly encountered species in the small basidiomycetous sub-phylum Wallemiomycotina is Wallemia mellicola, a xerotolerant fungus with a widespread distribution. To investigate the population characteristics of the species, whole genomes of twenty-five strains were sequenced. Apart from identification of four strains of clonal origin, the distances between the genomes failed to reflect either the isolation habitat of the strains or their geographical origin. Strains from different parts of the world appeared to represent a relatively homogenous and widespread population. The lack of concordance between individual gene phylogenies and the decay of linkage disequilibrium indicated that W. mellicola is at least occasionally recombining. Two versions of a putative mating-type locus have been found in all sequenced genomes, each present in approximately half of the strains. W. mellicola thus appears to be capable of (sexual) recombination and shows no signs of allopatric speciation or specialization to specific habitats.

A few Ascomycota taxa dominate soil fungal communities worldwide

Despite having key functions in terrestrial ecosystems, information on the dominant soil fungi and their ecological preferences at the global scale is lacking. To fill this knowledge gap, we surveyed 235 soils from across the globe. Our findings indicate that 83 phylotypes (<0.1% of the retrieved fungi), mostly belonging to wind dispersed, generalist Ascomycota, dominate soils globally. We identify patterns and ecological drivers of dominant soil fungal taxa occurrence, and present a map of their distribution in soils worldwide. Whole-genome comparisons with less dominant, generalist fungi point at a significantly higher number of genes related to stress-tolerance and resource uptake in the dominant fungi, suggesting that they might be better in colonising a wide range of environments. Our findings constitute a major advance in our understanding of the ecology of fungi, and have implications for the development of strategies to preserve them and the ecosystem functions they provide.

Soil fungi play essential roles in ecosystems worldwide. Here, the authors sequence and analyze 235 soil samples collected from across the globe, and identify dominant fungal taxa and their associated environmental attributes.

Comparative genomics of 40 edible and medicinal mushrooms provide an insight into the evolution of lignocellulose decomposition mechanisms

Diversity comparison and phylogenetic analyses of carbohydrate-active enzymes (CAZymes), auxiliary activities (AAs) and cytochromes P450 among 40 fungi, which are based on different nutritional pathways, help clarify and explain their divergence and improvement of various life-styles. Molecular clock analyses allow us to understand the evolutionary and developmental rules in decomposition gene families. Our results suggested that fungi in different ecological types acquired an obvious preference on specific decomposing gene families during evolutionary selection. White rot and litter saprotrophic fungi possessed more complete types of varied degradation gene families and were superior in quantities. With evolution and development of lignocellulose decomposition mechanism, certain families (like CBM1, GH6, GH7, GH10, and CYP53) disappeared in brown rot fungi and symbiotic fungi. In addition, the earlier time of phylogenetic divergence determined the more integrated and larger decomposition families. And various gains and losses in gene quantity of varied decomposition families led in particularly phylogenetic clades or nodes, then accelerated in forming varied ecotypes of species.

Plant-mediated partner discrimination in ectomycorrhizal mutualisms

Although ectomycorrhizal fungi have well-recognized effects on ecological processes ranging from plant community dynamics to carbon cycling rates, it is unclear if plants are able to actively influence the structure of these fungal communities. To address this knowledge gap, we performed two complementary experiments to determine (1) whether ectomycorrhizal plants can discriminate among potential fungal partners, and (2) to what extent the plants might reward better mutualists. In experiment 1, split-root Larix occidentalis seedlings were inoculated with spores from three Suillus species (S. clintonianus, S. grisellus, and S. spectabilis). In experiment 2, we manipulated the symbiotic quality of Suillus brevipes isolates on split-root Pinus muricata seedlings by changing the nitrogen resources available, and used carbon-13 labeling to track host investment in fungi. In experiment 1, we found that hosts can discriminate in multi-species settings. The split-root seedlings inhibited colonization by S. spectabilis whenever another fungus was available, despite similar benefits from all three fungi. In experiment 2, we found that roots and fungi with greater nitrogen supplies received more plant carbon. Our results suggest that plants may be able to regulate this symbiosis at a relatively fine scale, and that this regulation can be integrated across spatially separated portions of a root system.

The Cedrus-associated truffle Trappeindia himalayensis is a morphologically unique and phylogenetically divergent species of Rhizopogon

In the northwestern Himalayan mountains of India, the hypogeous sequestrate fungus Trappeindia himalayensis is harvested from forests dominated by the ectomycorrhizal tree Cedrus deodara (Himalayan cedar). This truffle has basidiospores that are ornamented with raised reticulation. The original description of Trappeindia himalayensis suggested that the gleba of this species is similar to young specimens of Scleroderma (Boletales), whereas its basidiospores are ornamented with raised reticulation, suggesting a morphological affinity to Leucogaster (Russulales) or Strobilomyces (Boletales). Given this systematic ambiguity, we have generated DNA sequence data from type material and other herbarium specimens and present the first molecular phylogenetic analysis of this unusual Cedrus-associated truffle. Despite the irregular ornamented basidiospore morphology, T. himalayensis is resolved within the genus Rhizopogon (Suillineae, Boletales) and represents a unique lineage that has not been previously detected. All known Rhizopogon species possess an ectomycorrhizal trophic mode, and because of its placement in this lineage, it is likely that Trappeindia himalayensis is an ectomycorrhizal partner of Cedrus deodara. This study highlights the importance of generating sequence data from herbarium specimens in order to identify fungal biodiversity and clarify the systematic relationships of poorly documented fungi.

First genome-wide analysis of the endangered, endemic lichen Cetradonia linearis reveals isolation by distance and strong population structure

PREMISE OF THE STUDY

Lichenized fungi are evolutionarily diverse and ecologically important, but little is known about the processes that drive their diversification and genetic differentiation. Distributions are often assumed to be wholly shaped by ecological requirements rather than dispersal limitations. Furthermore, although asexual and sexual reproductive structures are observable, the lack of information about recombination rates makes inferences about reproductive strategies difficult. We investigated the population genomics of Cetradonia linearis, a federally endangered lichen in the southern Appalachians of eastern North America, to test the relative contributions of environmental and geographic distance in shaping genetic structure, and to characterize the mating system and genome-wide recombination.

METHODS

Whole-genome shotgun sequencing was conducted to generate data for 32 individuals of C. linearis. A reference genome was assembled, and reads from all samples were aligned to generate a set of single-nucleotide polymorphisms for further analyses.

KEY RESULTS

We found evidence for low rates of recombination and for isolation by distance, but not for isolation by environment. The species is putatively unisexual, given that only one mating-type locus was found. Hindcast species distribution models and the distribution of genetic diversity support C. linearis having a larger range during the Last Glacial Maximum in the southern portion of its current extent.

CONCLUSIONS

Our findings contribute to the understanding of factors that shape genetic diversity in C. linearis and in fungi more broadly. Because all populations are highly genetically differentiated, the extirpation of any population would mean the loss of unique genetic diversity; therefore, our results support the continued conservation of this species.

The origin and adaptive evolution of domesticated populations of yeast from Far East Asia

The yeast Saccharomyces cerevisiae has been an essential component of human civilization because of its long global history of use in food and beverage fermentation. However, the diversity and evolutionary history of the domesticated populations of the yeast remain elusive. We show here that China/Far East Asia is likely the center of origin of the domesticated populations of the species. The domesticated populations form two major groups associated with solid- and liquid-state fermentation and appear to have originated from heterozygous ancestors, which were likely formed by outcrossing between diverse wild isolates primitively for adaptation to maltose-rich niches. We found consistent gene expansion and contraction in the whole domesticated population, as well as lineage-specific genome variations leading to adaptation to different environments. We show a nearly panoramic view of the diversity and life history of S. cerevisiae and provide new insights into the origin and evolution of the species.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

Population genomics of Fusarium graminearum reveals signatures of divergent evolution within a major cereal pathogen

The cereal pathogen Fusarium graminearum is the primary cause of Fusarium head blight (FHB) and a significant threat to food safety and crop production. To elucidate population structure and identify genomic targets of selection within major FHB pathogen populations in North America we sequenced the genomes of 60 diverse F. graminearum isolates. We also assembled the first pan-genome for F. graminearum to clarify population-level differences in gene content potentially contributing to pathogen diversity. Bayesian and phylogenomic analyses revealed genetic structure associated with isolates that produce the novel NX-2 mycotoxin, suggesting a North American population that has remained genetically distinct from other endemic and introduced cereal-infecting populations. Genome scans uncovered distinct signatures of selection within populations, focused in high diversity, frequently recombining regions. These patterns suggested selection for genomic divergence at the trichothecene toxin gene cluster and thirteen additional regions containing genes potentially involved in pathogen specialization. Gene content differences further distinguished populations, in that 121 genes showed population-specific patterns of conservation. Genes that differentiated populations had predicted functions related to pathogenesis, secondary metabolism and antagonistic interactions, though a subset had unique roles in temperature and light sensitivity. Our results indicated that F. graminearum populations are distinguished by dozens of genes with signatures of selection and an array of dispensable accessory genes, suggesting that FHB pathogen populations may be equipped with different traits to exploit the agroecosystem. These findings provide insights into the evolutionary processes and genomic features contributing to population divergence in plant pathogens, and highlight candidate genes for future functional studies of pathogen specialization across evolutionarily and ecologically diverse fungi.

Dissimilar pigment regulation in Serpula lacrymans and Paxillus involutus during inter-kingdom interactions

Production of basidiomycete atromentin-derived pigments like variegatic acid (pulvinic acid-type) and involutin (diarylcyclopentenone) from the brown-rotter Serpula lacrymans and the ectomycorrhiza-forming Paxillus involutus, respectively, is induced by complex nutrition, and in the case of S. lacrymans, bacteria. Pigmentation in S. lacrymans was stimulated by 13 different bacteria and cell-wall-damaging enzymes (lytic enzymes and proteases), but not by lysozyme or mechanical damage. The use of protease inhibitors with Bacillus subtilis or heat-killed bacteria during co-culturing with S. lacrymans significantly reduced pigmentation indicating that enzymatic hyphal damage and/or released peptides, rather than mechanical injury, was the major cause of systemic pigment induction. Conversely, no significant pigmentation by bacteria was observed from P. involutus. We found additional putative transcriptional composite elements of atromentin synthetase genes in P. involutus and other ectomycorrhiza-forming species that were absent from S. lacrymans and other brown-rotters. Variegatic and its precursor xerocomic acid, but not involutin, in return inhibited swarming and colony biofilm spreading of Bacillus subtilis, but did not kill B. subtilis. We suggest that dissimilar pigment regulation by fungal lifestyle was a consequence of pigment bioactivity and additional promoter motifs. The focus on basidiomycete natural product gene induction and regulation will assist in future studies to determine global regulators, signalling pathways and associated transcription factors of basidiomycetes.