Fungal heavy metal adaptation through single nucleotide polymorphisms and copy‐number variation

Inter- and intra-specific metal tolerance variation in ectomycorrhizal fungal Suillus species

Soil metal contamination can affect growth, metabolism, and reproduction of organisms, and can lead to death. However, some fungi have evolved metal tolerance and are able to live in contaminated soils. Species in the ectomycorrhizal genus Suillus from Europe and Asia display variation in metal tolerance, yet it is unknown whether this is a widespread trait in the genus and whether it occurs in North America. Here we investigate cadmium (Cd) and zinc (Zn) tolerance in S. brevipes and S. tomentosus isolates collected from sites in the Rocky Mountains of Colorado displaying different metal content. In line with previous findings for other Suillus species, we hypothesized (1) S. brevipes and S. tomentosus to display intra-specific metal tolerance variation, (2) Zn and Cd tolerance to be correlated to soil metal content, and (3) tolerant isolates to show lower metal tissue content compared to sensitive isolates (due to increased metal exclusion). We found ample intra- and inter-specific Zn and Cd tolerance variation in both S. brevipes and S. tomentosus, but no correlation between soil metal content and tolerance. There was a negative correlation between tolerance level and Zn uptake, indicating an exclusion-based Zn tolerance strategy. Sensitive and tolerant isolates showed no difference in Cd accumulation, indicating that Cd tolerance in these species is likely not dependent on exclusion. Our study sets the groundwork for further investigation into the genetic basis of Suillus metal tolerance and whether and how it impacts pine mycorrhizal 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.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

That's no moon, it's a Starship: Giant transposons driving fungal horizontal gene transfer

To date, most reports of horizontal gene transfer (HGT) in fungi rely on genome sequence data and are therefore an indirect measure of HGT after the event has occurred. However, a novel group of class II-like transposons known as Starships may soon alter this status quo. Starships are giant transposable elements that carry dozens of genes, some of which are host-beneficial, and are linked to many recent HGT events in the fungal kingdom. These transposons remain active and mobile in many fungal genomes and their transposition has recently been shown to be driven by a conserved tyrosine-recombinase called 'Captain'. This perspective explores some of the remaining unanswered questions about how these Starship transposons move, both within a genome and between different species. We seek to outline several experimental approaches that can be used to identify the genes essential for Starship-mediated HGT and draw links to other recently discovered giant transposons outside of the fungal kingdom.

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.

Re-thinking the environment in landscape genomics

Detecting the extrinsic selective pressures shaping genomic variation is critical for a better understanding of adaptation and for forecasting evolutionary responses of natural populations to changing environmental conditions. With increasing availability of geo-referenced environmental data, landscape genomics provides unprecedented insights into how genomic variation and underlying gene functions affect traits potentially under selection. Yet, the robustness of genotype-environment associations used in landscape genomics remains tempered due to various limitations, including the characteristics of environmental data used, sampling designs employed, and statistical frameworks applied. Here, we argue that using complementary or new environmental data sources and well-informed sampling designs may help improve the detection of selective pressures underlying patterns of local adaptation in various organisms and environments.

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.

Water Quality, Heavy Metals, and Antifungal Susceptibility to Fluconazole of Yeasts from Water Systems

Aquatic environments could be reservoirs of pathogenic yeasts with acquired antifungal resistance. The susceptibility to antifungal agents of yeasts present in the wastewater and natural waters of the city of Cali was evaluated. Samples were taken from two types of water: drinking water (Meléndez River, drinking water treatment plant "Puerto Mallarino" in the Cauca River) and wastewater (South Channel of the Cauca River, "Cañaveralejo-PTAR" wastewater treatment plant). Physico-chemical parameters, heavy metal concentration, and yeast levels were determined using standard procedures. Yeasts were identified using API 20 C AUX (BioMérieux) and sequence analysis of the ITS1-5.8S-ITS2 and D1/D2 regions of the large subunit of the ribosome. Susceptibility assays against fluconazole and amphotericin B using the minimum inhibitory concentration (MIC) test were determined using the microdilution method. The influence of physico-chemical parameters and heavy metals was established using principal component analysis (PCA). Yeast counts were higher at WWTP "PTAR" and lower at Melendez River, as expected. A total of 14 genera and 21 yeast species was identified, and the genus Candida was present at all locations. Susceptibility tests showed a 32.7% resistance profile to fluconazole in the order DWTP "Puerto Mallarino = WWTP "PTAR" > South Channel "Navarro". There were significant differences (p < 0.05) in the physico-chemical parameters/concentration of heavy metals and yeast levels between the aquatic systems under study. A positive association was observed between yeast levels and total dissolved solids, nitrate levels, and Cr at the "PTAR" WWTP; conductivity, Zn, and Cu in the South Channel; and the presence of Pb in the "Puerto Mallarino" DWTP. Rhodotorula mucilaginosa, Candida albicans, and Candida sp. 1 were influenced by Cr and Cd, and Diutina catelunata was influenced by Fe (p < 0.05). The water systems explored in this study showed different yeast levels and susceptibility profiles, and, therefore, possible genetic differences among populations of the same species, and different physico-chemical and heavy metals concentrations, which were probably modulating the antifungal-resistant yeasts. All these aquatic systems discharge their content into the Cauca River. We highlight the importance to further investigate if these resistant communities continue to other locations in the second largest river of Colombia and to determine the risk posed to humans and animals.

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.

A large transposable element mediates metal resistance in the fungus Paecilomyces variotii

The horizontal transfer of large gene clusters by mobile elements is a key driver of prokaryotic adaptation in response to environmental stresses. Eukaryotic microbes face similar stresses; however, a parallel role for mobile elements has not been established. A stress faced by many microorganisms is toxic metal ions in their environment. In fungi, identified mechanisms for protection against metals generally rely on genes that are dispersed within an organism's genome. Here, we discover a large (∼85 kb) region that confers tolerance to five metal/metalloid ions (arsenate, cadmium, copper, lead, and zinc) in the genomes of some, but not all, strains of a fungus, Paecilomyces variotii. We name this region HEPHAESTUS (Hφ) and present evidence that it is mobile within the P. variotii genome with features characteristic of a transposable element. HEPHAESTUS contains the greatest complement of host-beneficial genes carried by a transposable element in eukaryotes, suggesting that eukaryotic transposable elements might play a role analogous to bacteria in the horizontal transfer of large regions of host-beneficial DNA. Genes within HEPHAESTUS responsible for individual metal tolerances include those encoding a P-type ATPase transporter-PcaA-required for cadmium and lead tolerance, a transporter-ZrcA-providing tolerance to zinc, and a multicopper oxidase-McoA-conferring tolerance to copper. In addition, a subregion of Hφ confers tolerance to arsenate. The genome sequences of other fungi in the Eurotiales contain further examples of HEPHAESTUS, suggesting that it is responsible for independently assembling tolerance to a diverse array of ions, including chromium, mercury, and sodium.

Cryptic genetic structure and copy-number variation in the ubiquitous forest symbiotic fungus Cenococcum geophilum

Ectomycorrhizal (ECM) fungi associated with plants constitute one of the most successful symbiotic interactions in forest ecosystems. ECM support trophic exchanges with host plants and are important factors for the survival and stress resilience of trees. However, ECM clades often harbour morpho-species and cryptic lineages, with weak morphological differentiation. How this relates to intraspecific genome variability and ecological functioning is poorly known. Here, we analysed 16 European isolates of the ascomycete Cenococcum geophilum, an extremely ubiquitous forest symbiotic fungus with no known sexual or asexual spore-forming structures but with a massively enlarged genome. We carried out whole-genome sequencing to identify single-nucleotide polymorphisms. We found no geographic structure at the European scale but divergent lineages within sampling sites. Evidence for recombination was restricted to specific cryptic lineages. Lineage differentiation was supported by extensive copy-number variation. Finally, we confirmed heterothallism with a single MAT1 idiomorph per genome. Synteny analyses of the MAT1 locus revealed substantial rearrangements and a pseudogene of the opposite MAT1 idiomorph. Our study provides the first evidence for substantial genome-wide structural variation, lineage-specific recombination and low continent-wide genetic differentiation in C. geophilum. Our study provides a foundation for targeted analyses of intra-specific functional variation in this major symbiosis.

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.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

Fungal-Metal Interactions: A Review of Toxicity and Homeostasis

Metal nanoparticles used as antifungals have increased the occurrence of fungal-metal interactions. However, there is a lack of knowledge about how these interactions cause genomic and physiological changes, which can produce fungal superbugs. Despite interest in these interactions, there is limited understanding of resistance mechanisms in most fungi studied until now. We highlight the current knowledge of fungal homeostasis of zinc, copper, iron, manganese, and silver to comprehensively examine associated mechanisms of resistance. Such mechanisms have been widely studied in Saccharomyces cerevisiae, but limited reports exist in filamentous fungi, though they are frequently the subject of nanoparticle biosynthesis and targets of antifungal metals. In most cases, microarray analyses uncovered resistance mechanisms as a response to metal exposure. In yeast, metal resistance is mainly due to the down-regulation of metal ion importers, utilization of metallothionein and metallothionein-like structures, and ion sequestration to the vacuole. In contrast, metal resistance in filamentous fungi heavily relies upon cellular ion export. However, there are instances of resistance that utilized vacuole sequestration, ion metallothionein, and chelator binding, deleting a metal ion importer, and ion storage in hyphal cell walls. In general, resistance to zinc, copper, iron, and manganese is extensively reported in yeast and partially known in filamentous fungi; and silver resistance lacks comprehensive understanding in both.

Tolerance of three fungal species to lithium and cobalt: Implications for bioleaching of spent rechargeable Li-ion batteries

AIMS

This paper aims to quantify the growth and organic acid production of Aspergillus niger, Penicillium chrysogenum and Penicillium simplicissimum when these fungi are exposed to varying levels of lithium (Li) and cobalt (Co). The study also tests whether pre-exposing the fungi to these metals enables the fungi to develop tolerance to Li or Co.

METHODS AND RESULTS

When cultures of A. niger, P. chrysogenum or P. simplicissimum were exposed to 250 mg l^-1 of Li or Co, biomass production and excretion of organic acids were significantly inhibited after 5 days of growth compared to cultures grown in the absence of these metals. Pre-exposing cultures of A. niger to 250 mg l^-1 of Li or Co for 20 days significantly increased biomass production when the fungus was subsequently sub-cultured into 250 or 500 mg l^-1 of Li or Co. However, pre-exposure of P. chrysogenum or P. simplicissimum to 250 mg l^-1 of Li or Co for 20 days did not increase biomass production.

CONCLUSIONS

Aspergillus niger, but not the Penicillium species, developed tolerance to Li and to Co during the 20-day pre-exposure period. Therefore, processes that utilize fungal bioleaching with A. niger to mobilize and recover valuable metals such as Li or Co should consider a pre-exposure step for fungi to improve their tolerance to metal toxicity.

SIGNIFICANCE AND IMPACT OF THE STUDY

Fungi may have the ability to extract valuable metals such as Li and Co from spent rechargeable batteries. However, the toxicity of the extracted metals can inhibit fungal growth and organic acid production. Pre-exposure to metals may alleviate toxicity for some fungal species. This knowledge can be used to improve the design of bioleaching protocols, increasing the potential for fungal bioleaching to become an economical and environmentally friendly method of recovering Li and Co from spent batteries.

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.