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.

The sizes of life.

Recent research has revealed the diversity and biomass of life across ecosystems, but how that biomass is distributed across body sizes of all living things remains unclear. We compile the present-day global body size-biomass spectra for the terrestrial, marine, and subterranean realms. To achieve this compilation, we pair existing and updated biomass estimates with previously uncatalogued body size ranges across all free-living biological groups. These data show that many biological groups share similar ranges of body sizes, and no single group dominates size ranges where cumulative biomass is highest. We then propagate biomass and size uncertainties and provide statistical descriptions of body size-biomass spectra across and within major habitat realms. Power laws show exponentially decreasing abundance (exponent -0.9±0.02 S.D., R2 = 0.97) and nearly equal biomass (exponent 0.09±0.01, R2 = 0.56) across log size bins, which resemble previous aquatic size spectra results but with greater organismal inclusivity and global coverage. In contrast, a bimodal Gaussian mixture model describes the biomass pattern better (R2 = 0.86) and suggests small (~10-15 g) and large (~107 g) organisms outweigh other sizes by one order magnitude (15 and 65 Gt versus ~1 Gt per log size). The results suggest that the global body size-biomass relationships is bimodal, but substantial one-to-two orders-of-magnitude uncertainty mean that additional data will be needed to clarify whether global-scale universal constraints or local forces shape these patterns.

Fungal Community Shift Along Steep Environmental Gradients from Geothermal Soils in Yellowstone National Park.

Geothermal soils offer unique insight into the way extreme environmental factors shape communities of organisms. However, little is known about the fungi growing in these environments and in particular how localized steep abiotic gradients affect fungal diversity. We used metabarcoding to characterize soil fungi surrounding a hot spring-fed thermal creek with water up to 84 °C and pH 10 in Yellowstone National Park. We found a significant association between fungal communities and soil variable principal components, and we identify the key trends in co-varying soil variables that explain the variation in fungal community. Saprotrophic and ectomycorrhizal fungi community profiles followed, and were significantly associated with, different soil variable principal components, highlighting potential differences in the factors that structure these different fungal trophic guilds. In addition, in vitro growth experiments in four target fungal species revealed a wide range of tolerances to pH levels but not to heat. Overall, our results documenting turnover in fungal species within a few hundred meters suggest many co-varying environmental factors structure the diverse fungal communities found in the soils of Yellowstone National Park.

Phylogenetic analysis of the distribution of deadly amatoxins among the little brown mushrooms of the genus Galerina.

Some but not all of the species of 'little brown mushrooms' in the genus Galerina contain deadly amatoxins at concentrations equaling those in the death cap, Amanita phalloides. However, Galerina's ~300 species are notoriously difficult to identify by morphology, and the identity of toxin-containing specimens has not been verified with DNA barcode sequencing. This left open the question of which Galerina species contain toxins and which do not. We selected specimens for toxin analysis using a preliminary phylogeny of the fungal DNA barcode region, the ribosomal internal transcribed spacer (ITS) region. Using liquid chromatography/mass spectrometry, we analyzed amatoxins from 70 samples of Galerina and close relatives, collected in western British Columbia, Canada. To put the presence of toxins into a phylogenetic context, we included the 70 samples in maximum likelihood analyses of 438 taxa, using ITS, RNA polymerase II second largest subunit gene (RPB2), and nuclear large subunit ribosomal RNA (LSU) gene sequences. We sequenced barcode DNA from types where possible to aid with applications of names. We detected amatoxins only in the 24 samples of the G. marginata s.l. complex in the Naucoriopsis clade. We delimited 56 putative Galerina species using Automatic Barcode Gap Detection software. Phylogenetic analysis showed moderate to strong support for Galerina infrageneric clades Naucoriopsis, Galerina, Tubariopsis, and Sideroides. Mycenopsis appeared paraphyletic and included Gymnopilus. Amatoxins were not detected in 46 samples from Galerina clades outside of Naucoriopsis or from outgroups. Our data show significant quantities of toxin in all mushrooms tested from the G. marginata s.l. complex. DNA barcoding revealed consistent accuracy in morphology-based identification of specimens to G. marginata s.l. complex. Prompt and careful morphological identification of ingested G. marginata s.l. has the potential to improve patient outcomes by leading to fast and appropriate treatment.

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.

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

Human-altered environments can shape the evolution of organisms. Fungi are no exception, although little is known about how they withstand anthropogenic pollution. Here, we document adaptation in the mycorrhizal fungus Suillus luteus driven by soil heavy metal contamination. Genome scans across individuals from recently polluted and nearby unpolluted soils in Belgium revealed low divergence across isolates and no evidence of population structure based on soil type. However, we detected single nucleotide polymorphism divergence and gene copy-number variation, with different genetic combinations potentially conferring the ability to persist in contaminated soils. Variants were shared across the population but found to be under selection in isolates exposed to pollution and located across the genome, including in genes involved in metal exclusion, storage, immobilization and reactive oxygen species detoxification. Together, our results point to S. luteus undergoing the initial steps of adaptive divergence and contribute to understanding the processes underlying local adaptation under strong environmental selection.