Continental-level population differentiation and environmental adaptation in the mushroom Suillus brevipes

Ectomycorrhizal fungi enhance pine growth by stimulating iron-dependent mechanisms with trade-offs in symbiotic performance

Iron (Fe) is crucial for metabolic functions of living organisms. Plants access occluded Fe through interactions with rhizosphere microorganisms and symbionts. Yet, the interplay between Fe addition and plant-mycorrhizal interactions, especially the molecular mechanisms underlying mycorrhiza-assisted Fe processing in plants, remains largely unexplored. We conducted mesocosms in Pinus plants inoculated with different ectomycorrhizal fungi (EMF) Suillus species under conditions with and without Fe coatings. Meta-transcriptomic, biogeochemical, and X-ray fluorescence imaging analyses were applied to investigate early-stage mycorrhizal roots. While Fe addition promoted Pinus growth, it concurrently reduced mycorrhiza formation rate, symbiosis-related metabolites in plant roots, and aboveground plant carbon and macronutrient content. This suggested potential trade-offs between Fe-enhanced plant growth and symbiotic performance. However, the extent of this trade-off may depend on interactions between host plants and EMF species. Interestingly, dual EMF species were more effective at facilitating plant Fe uptake by inducing diverse Fe-related functions than single-EMF species. This subsequently triggered various Fe-dependent physiological and biochemical processes in Pinus roots, significantly contributing to Pinus growth. However, this resulted in a greater carbon allocation to roots, relatively reducing the aboveground plant carbon content. Our study offers critical insights into how EMF communities rebalance benefits of Fe-induced effects on symbiotic partners.

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.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

Population genetic structure of Hymenopellis radicata germplasm resources based on genome re-sequencing

Through whole-genome re-sequencing of 18 Hymenopellis radicata germplasm resources collected from diverse regions in China, we identified significant variations in the form of Single Nucleotide Polymorphisms (SNPs) and Insertions and Deletions (InDels). These variations were comprehensively annotated, shedding light on the mutation types present in the entire genome of the H. radicata germplasm. This analysis revealed the number and position information of each mutation and provided insights into the overall genomic landscape of H. radicata germplasm. Utilizing SNP data, we delved into the population structure of the 18 H. radicata germplasm resources. The results indicated the presence of 2,335,179 Indel sites and 12,050,448 SNP sites. The population structure analysis unveiled two distinct subgroups among the H. radicata germplasm resources. Phenotypic statistics, principal component analysis, and phylogenetic tree results echoed the findings of the population structure analysis. Different strains of H. radicata from various regions in China exhibited notable differences in genetic diversity, mycelial growth rate, yield, and fruiting body characteristics. Significant disparities were observed between the two subgroups, while strains within each subgroup shared common characteristics. This research establishes a solid foundation for integrating H. radicata into diverse breeding programs. The data underscore the potential of H. radicata for genetic improvement and exploitation in breeding initiatives, paving the way for future advancements in this field.

Whole-genome resequencing reveals genetic diversity and selection signals in warm temperate and subtropical Sillago sinica populations

BACKGROUND

Genetic diversity and heterogeneous genomic signatures in marine fish populations may result from selection pressures driven by the strong effects of environmental change. Nearshore fishes are often exposed to complex environments and human activities, especially those with small ranges. However, studies on genetic diversity and population selection signals in these species have mostly been based on a relatively small number of genetic markers. As a newly recorded species of Sillaginidae, the population genetics and genomic selection signals of Sillago sinica are fragmented or even absent.

RESULTS

To address this theoretical gap, we performed whole-genome resequencing of 43 S. sinica individuals from Dongying (DY), Qingdao (QD) and Wenzhou (WZ) populations and obtained 4,878,771 high-quality SNPs. Population genetic analysis showed that the genetic diversity of S. sinica populations was low, but the genetic diversity of the WZ population was higher than that of the other two populations. Interestingly, the three populations were not strictly clustered within the group defined by their sampling location but showed an obvious geographic structure signal from the warm temperate to the subtropics. With further analysis, warm-temperate populations exhibited strong selection signals in genomic regions related to nervous system development, sensory function and immune function. However, subtropical populations showed more selective signalling for environmental tolerance and stress signal transduction.

CONCLUSIONS

Genome-wide SNPs provide high-quality data to support genetic studies and localization of selection signals in S. sinica populations. The reduction in genetic diversity may be related to the bottleneck effect. Considering that low genetic diversity leads to reduced environmental adaptability, conservation efforts and genetic diversity monitoring of this species should be increased in the future. Differences in genomic selection signals between warm temperate and subtropical populations may be related to human activities and changes in environmental complexity. This study deepened the understanding of population genetics and genomic selection signatures in nearshore fishes and provided a theoretical basis for exploring the potential mechanisms of genomic variation in marine fishes driven by environmental selection pressures.

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.

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.

Contrasting continental patterns of adaptive population divergence in the holarctic ectomycorrhizal fungus Boletus edulis

In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome-wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis, a globally distributed species complex with broad ecological amplitude. By whole-genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphisms to characterize patterns of genome-wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62-2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene-Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome-wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation.

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.

A call to reconceptualize lichen symbioses

Several decades of research across disciplines have overturned historical perspectives of symbioses dominated by binary characterizations of highly specific species-species interactions. This paradigm shift has unlocked the previously underappreciated and overlooked dynamism of fungal mutualisms such as mycorrhizae. Lichens are another example of important fungal mutualisms where reconceptualization is urgently needed to realize their potential as model systems. This reconceptualization requires both an objective synthesis of new data and envisioning a revised integrative approach that unifies the spectrum of ecology and evolution. We propose a ten-theme framework that if pursued would propel lichens to the vanguard of symbiotic theory.

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.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.

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.

Contrasting demographic histories revealed in two invasive populations of the dry rot fungus Serpula lacrymans

Globalization and international trade have impacted organisms around the world leading to a considerable number of species establishing in new geographic areas. Many organisms have taken advantage of human-made environments, including buildings. One such species is the dry rot fungus Serpula lacrymans, which is the most aggressive wood-decay fungus in indoor environments in temperate regions. Using population genomic analyses of 36 full genome sequenced isolates, we demonstrated that European and Japanese isolates are highly divergent and the populations split 3000-19,000 generations ago, probably predating human influence. Approximately 250 generations ago, the European population went through a tight bottleneck, probably corresponding to the fungus colonization of the built environment in Europe. The demographic history of these populations, probably lead to low adaptive potential. Only two loci under selection were identified using a F_st outlier approach, and selective sweep analyses identified three loci with extended haplotype homozygosity. The selective sweep analyses found signals in genes possibly related to decay of various substrates in Japan and in genes involved DNA replication and protein modification in Europe. Our results suggest that the dry rot fungus independently established in indoor environments in Europe and Japan and that invasive species can potentially establish large populations in new habitats based on a few colonizing individuals.

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.

Population structure and adaptation of a bacterial pathogen in California grapevines

Xylella fastidiosa subsp. fastidiosa causes Pierce's disease of grapevine (PD) and has been present in California for over a century. A singly introduced genotype spread across the state causing large outbreaks and damaging the grapevine industry. This study presents 122 X. fastidiosa subsp. fastidiosa genomes from symptomatic grapevines, and explores pathogen genetic diversity associated with PD in California. A total of 5218 single-nucleotide polymorphisms (SNPs) were found in the dataset. Strong population genetic structure was found; isolates split into five genetic clusters divided into two lineages. The core/soft-core genome constituted 41.2% of the total genome, emphasizing the high genetic variability of X. fastidiosa genomes. An ecological niche model was performed to estimate the environmental niche of the pathogen within California and to identify key climatic factors involved in dispersal. A landscape genomic approach was undertaken aiming to link local adaptation to climatic factors. A total of 18 non-synonymous polymorphisms found to be under selective pressures were correlated with at least one environmental variable highlighting the role of temperature, precipitation and elevation on X. fastidiosa adaptation to grapevines in California. Finally, the contribution to virulence of three of the genes under positive selective pressure and of one recombinant gene was studied by reverse genetics.

Climate change influences mycorrhizal fungal-plant interactions, but conclusions are limited by geographical study bias

Climate change is altering the interactions among plants and soil organisms in ways that will alter the structure and function of ecosystems. We reviewed the literature and developed a map of studies focused on how the three most common types of mycorrhizal fungi (arbuscular mycorrhizal [AM], ectomycorrhizal [EcM], and ericoid mycorrhizal [ErM] fungi) respond to elevated atmospheric carbon dioxide concentrations (eCO₂ ), climatic warming, and changes in the distribution of precipitation. Broadly, we ask how do mycorrhizal fungi respond to climate change, how do these responses vary by fungal type, and how do mycorrhizal traits influence plant adaptation, movement, or extinction in response to climatic change? First, we found that 92% of studies were conducted in the northern hemisphere, and plant host, ecosystem type and study location were only correlated with each other in the northern hemisphere because studies across all mycorrhizal fungal types were only common in the northern hemisphere. Second, we show that temperature and rainfall variability had more variable effects than eCO₂ on mycorrhizal fungal structures, but these effects were context dependent. Third, while mycorrhizal fungal types vary in their responses to climate change, it appears that warming leads to more variable responses in ectomycorrhizal than in arbuscular mycorrhizal fungi. Finally, we discuss common traits of mycorrhizal fungi that could aid in fungal and plant adaption to climate change. We posit that mycorrhizal fungi can buffer plant hosts against extinction risk, they can facilitate or retard the dispersal success of plants moving away from poor environments, and, by buffering host plants, they can enable host plant adaptation to new climates. All of these influences are, however, context dependent a finding that reflects the complex traits of mycorrhizal fungi as a group, the diversity of plant species they associate with and the variation in ecosystems in which they reside. Overall, while we point out many gaps in our understanding of the influence of climate changes on mycorrhizal fungi, we also highlight the large number of opportunities for researching plant and mycorrhizal fungal responses to and mitigation of climate changes.

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.

Genomic Analyses Reveal Evidence of Independent Evolution, Demographic History, and Extreme Environment Adaptation of Tibetan Plateau Agaricus bisporus

Agaricus bisporus distributed in the Tibetan Plateau of China has high-stress resistance that is valuable for breeding improvements. However, its evolutionary history, specialization, and adaptation to the extreme Tibetan Plateau environment are largely unknown. Here, we performed de novo genome sequencing of a representative Tibetan Plateau wild strain ABM and comparative genomic analysis with the reported European strain H97 and H39. The assembled ABM genome was 30.4 Mb in size, and comprised 8,562 protein-coding genes. The ABM genome shared highly conserved syntenic blocks and a few inversions with H97 and H39. The phylogenetic tree constructed by 1,276 single-copy orthologous genes in nine fungal species showed that the Tibetan Plateau and European A. bisporus diverged ∼5.5 million years ago. Population genomic analysis using genome resequencing of 29 strains revealed that the Tibetan Plateau population underwent significant differentiation from the European and American populations and evolved independently, and the global climate changes critically shaped the demographic history of the Tibetan Plateau population. Moreover, we identified key genes that are related to the cell wall and membrane system, and the development and defense systems regulated A. bisporus adapting to the harsh Tibetan Plateau environment. These findings highlight the value of genomic data in assessing the evolution and adaptation of mushrooms and will enhance future genetic improvements of A. bisporus.

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.

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.

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.

Beyond ICOM8: perspectives on advances in mycorrhizal research from 2015 to 2017

This report reviews important advances in mycorrhizal research that occurred during the past 2 years. We highlight major advancements both within and across levels of biological organization and describe areas where greater integration has led to unique insights. Particularly active areas of research include exploration of the mechanisms underpinning the development of the mycorrhizal symbiosis, the mycorrhizal microbiome, comparisons among types of mycorrhizas from molecular to ecosystem scales, the extent and function of mycorrhizal networks and enhanced understanding of the role of mycorrhizas in carbon dynamics from local to global scales. The top-tier scientific journals have acknowledged mycorrhizas to be complex adaptive systems that play key roles in the development of communities and ecosystem processes. Understanding the mechanisms driving these large-scale effects requires integration of knowledge across scales of biological organization.

Accounting for local adaptation in ectomycorrhizas: a call to track geographical origin of plants, fungi, and soils in experiments

Local adaptation, the differential success of genotypes in their native versus foreign environments, can influence ecological and evolutionary processes, yet its importance is difficult to estimate because it has not been widely studied, particularly in the context of interspecific interactions. Interactions between ectomycorrhizal (EM) fungi and their host plants could serve as model system for investigations of local adaptation because they are widespread and affect plant responses to both biotic and abiotic selection pressures. Furthermore, because EM fungi cycle nutrients and mediate energy flow into food webs, their local adaptation may be critical in sustaining ecological function. Despite their ecological importance and an extensive literature on their relationships with plants, the vast majority of experiments on EM symbioses fail to report critical information needed to assess local adaptation: the geographic origin of the plant, fungal inocula, and soil substrate used in the experiment. These omissions limit the utility of such studies and restrict our understanding of EM ecology and evolution. Here, we illustrate the potential importance of local adaptation in EM relationships and call for consistent reporting of the geographic origin of plant, soil, and fungi as an important step towards a better understanding of the ecology and evolution of EM symbioses.