Logo: to the web site of Uppsala University

Identifying genuine polymorphic variants is a significant challenge in sequence data analysis, although detecting low-frequency variants in sequence data is essential for estimating demographic parameters and investigating genetic processes, such as selection, within populations. Arbuscular mycorrhizal (AM) fungi are multinucleate organisms, in which individual nuclei collectively operate as a population, and the extent of genetic variation across nuclei has long been an area of scientific interest. In this study, we investigated the patterns of polymorphism discovery and the alternate allele frequency distribution by comparing polymorphism discovery in 2 distinct genomic sequence datasets of the AM fungus model species, Rhizophagus irregularis strain DAOM197198. The 2 datasets used in this study are publicly available and were generated either from pooled spores and hyphae or amplified single nuclei from a single spore. We also estimated the intraorganismal variation within the DAOM197198 strain. Our results showed that the 2 datasets exhibited different frequency patterns for discovered variants. The whole-organism dataset showed a distribution spanning low-, intermediate-, and high-frequency variants, whereas the single-nucleus dataset predominantly featured low-frequency variants with smaller proportions in intermediate and high frequencies. Furthermore, single nucleotide polymorphism density estimates within both the whole organism and individual nuclei confirmed the low intraorganismal variation of the DAOM197198 strain and that most variants are rare. Our study highlights the methodological challenges associated with detecting low-frequency variants in AM fungal whole-genome sequence data and demonstrates that alternate alleles can be reliably identified in single nuclei of AM fungi.

Background: The colonization of land and the diversification of terrestrial plants is intimately linked to the evolutionary history of their symbiotic fungal partners. Extant representatives of these fungal lineages include mutualistic plant symbionts, the arbuscular mycorrhizal (AM) fungi in Glomeromycota and fine root endophytes in Endogonales (Mucoromycota), as well as fungi with saprotrophic, pathogenic and endophytic lifestyles. These fungal groups separate into three monophyletic lineages but their evolutionary relationships remain enigmatic confounding ancestral reconstructions. Their taxonomic ranks are currently fluid.

Results: In this study, we recognize these three monophyletic linages as phyla, and use a balanced taxon sampling and broad taxonomic representation for phylogenomic analysis that rejects a hard polytomy and resolves Glomeromycota as sister to a clade composed of Mucoromycota and Mortierellomycota. Low copy numbers of genes associated with plant cell wall degradation could not be assigned to the transition to a plant symbiotic lifestyle but appears to be an ancestral phylogenetic signal. Both plant symbiotic lineages, Glomeromycota and Endogonales, lack numerous thiamine metabolism genes but the lack of fatty acid synthesis genes is specific to AM fungi. Many genes previously thought to be missing specifically in Glomeromycota are either missing in all analyzed phyla, or in some cases, are actually present in some of the analyzed AM fungal lineages, e.g. the high affinity phosphorus transporter Pho89.

Conclusion: Based on a broad taxon sampling of fungal genomes we present a well-supported phylogeny for AM fungi and their sister lineages. We show that among these lineages, two independent evolutionary transitions to mutualistic plant symbiosis happened in a genomic background profoundly different from that known from the emergence of ectomycorrhizal fungi in Dikarya. These results call for further reevaluation of genomic signatures associated with plant symbiosis.

The development of carbon dioxide removal methods, coupled with decreased CO2 emissions, is fundamental to achieving the targets outlined in the Paris Agreement limiting global warming to 1.5 degrees C. Here we are investigating the importance of the organic carbon feedstock to support silicate mineral weathering in small-scale flow through bioreactors and subsequent CO2 sequestration. Here, we combine two bacteria and two fungi, widely reported for their weathering potential, in simple flow through bioreactors (columns) consisting of forsterite and widely available, cheap organic carbon sources (wheat straw, bio-waste digestate of pig manure and biowaste, and manure compost), over six weeks. Compared to their corresponding abiotic controls, the inoculated straw and digestate columns release more total alkalinity (similar to 2 times more) and produce greater dissolved and solid inorganic carbon (29% for straw and 13% for digestate), suggesting an increase in CO2 sequestration because of bio-enhanced silicate weathering. Microbial biomass is higher in the straw columns compared to the digestate and manure compost columns, with a phospholipid fatty acid derived total microbial biomass 10 x greater than the other biotic columns. Scanning Electron Microscopy imaging shows the most extensive colonisation and biofilm formation on the mineral surfaces in the straw columns. The biotic straw and digestate columns sequester 50 and 14 mg C more than their abiotic controls respectively, while there is no difference in the manure columns. The selection of organic carbon sources to support microbial communities in the flow through bioreactors controlls the silicate weathering rates and CO2 sequestration.

Forest fire is known to positively affect bark beetle populations by providing fire-damaged trees with impaired defenses for infestation. Tomicus piniperda, the common pine shoot beetle, breeds and lays eggs under the bark of stressed pine trees and is considered a serious forest pest within its native range. Wood-colonizing fungi have been hypothesized to improve substrate quality and detoxify tree defensive chemistry to indirectly facilitate tree colonization by beetles. While some bark beetle species form symbiotic associations with fungi and actively vector their partners when colonizing new trees, T. piniperda does not have mycangia or body hairs for specific vectoring of fungi. To explore the T. piniperda-associated fungal community for signs of specific association, we used ITS metabarcoding to separately characterize fungal communities associated with surface and gut of male and female beetles. We also characterized the temporal changes in fungal community and nutrient status of pine phloem with and without beetle galleries. Sampling was performed 2 years after a natural forest fire and included both burnt and unburnt sites. In our study system, we find that forest fire significantly impacts the fungal community composition associated with T. piniperda and that fire may also indirectly change nutrient availability in phloem to beetle galleries. We conclude that T. piniperda can vector fungi to newly colonized trees but the absence of positive effects on substrate quality and minimal effects of sex indicate that vectoring of associated fungal communities is not a strategy associated with the T. piniperda life cycle.

BackgroundArbuscular mycorrhizal (AM) fungi are arguably the most important symbionts of plants, offering a range of benefits to their hosts. However, the provisioning of these benefits does not appear to be uniform among AM fungal individuals, with genetic variation between fungal symbionts having a substantial impact on plant performance. Interestingly, genetic variation has also been reported within fungal individuals, which contain millions of haploid nuclei sharing a common cytoplasm. In the model AM fungus, Rhizophagus irregularis, several isolates have been reported to be dikaryotes, containing two genetically distinct types of nuclei recognized based on their mating-type (MAT) locus identity. However, their extremely coenocytic nature and lack of a known single nucleus stage has raised questions on the origin, distribution and dynamics of this genetic variation.ResultsHere we performed DNA and RNA sequencing at the mycelial individual, single spore and single nucleus levels to gain insight into the dynamic genetic make-up of the dikaryote-like R. irregularis C3 isolate and the effect of different host plants on its genetic variation. Our analyses reveal that parallel spore and root culture batches can have widely variable ratios of two main genotypes in C3. Additionally, numerous polymorphisms were found with frequencies that deviated significantly from the general genotype ratio, indicating a diverse population of slightly different nucleotypes. Changing host plants did not show consistent host effects on nucleotype ratio's after multiple rounds of subculturing. Instead, we found a major effect of host plant-identity on allele-specific expression in C3.ConclusionOur analyses indicate a highly dynamic/variable genetic organization in different isolates of R. irregularis. Seemingly random fluctuations in nucleotype ratio's upon spore formation, recombination events, high variability of non-tandemly repeated rDNA sequences and host-dependent allele expression all add levels of variation that may contribute to the evolutionary success of these widespread symbionts.

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large-scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus's time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic-dry part that has a species rich grass-dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%–40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.

Arbuscular mycorrhizal (AM) fungi are ubiquitous mutualistic symbionts of most terrestrial plants and many complete their lifecycles underground. Whole genome analysis of AM fungi has long been restricted to species and strains that can be maintained under controlled conditions that facilitate collection of biological samples. There is some evidence suggesting that AM fungi can adapt to culture resulting in phenotypic and possibly also genotypic changes in the fungi. In this study, we used field isolated spores of AM fungi and identified them as Funneliformis geosporum based on morphology and phylogenetic analyses. We separately assembled the genomes of two representative spores using DNA sequences of 19 and 22 individually amplified nuclei. The genomes were compared with previously published data from other members of Glomeraceae including two strains of F. mosseae. No significant differences were observed among the species in terms of gene content, while the single nucleotide polymorphism density was higher in the strains of F. geosporum than in the strains of F. mosseae. In this study, we demonstrate that it is possible to sequence and assemble genomes from AM fungal spores sampled in the field, which opens up the possibility to include uncultured AM fungi in phylogenomic and comparative genomic analysis and to study genomic variation in natural populations of these important plant symbionts.

Morphological characters and nuclear ribosomal DNA (rDNA) phylogenies have so far been the basis of the current classifications of arbuscular mycorrhizal (AM) fungi. Improved understanding of the evolutionary history of AM fungi requires extensive ortholog sampling and analyses of genome and transcriptome data from a wide range of taxa. To circumvent the need for axenic culturing of AM fungi we gathered and combined genomic data from single nuclei to generate de novo genome assemblies covering seven families of AM fungi. We successfully sequenced the genomes of 15 AM fungal species for which genome data was not previously available. Comparative analysis of the previously published Rhizophagus irregularis DAOM197198 assembly confirm that our novel workflow generates genome assemblies suitable for phylogenomic analysis. Predicted genes of our assemblies, together with published protein sequences of AM fungi and their sister clades, were used for phylogenomic analyses. We evaluated the phylogenetic placement of Glomeromycota in relation to its sister phyla (Mucoromycota and Mortierellomycota), and found no support to reject a polytomy. Finally, we explored the phylogenetic relationships within Glomeromycota. Our results support family level classification from previous phylogenetic studies, and the polyphyly of the order Glomerales with Claroideoglomeraceae as the sister group to Glomeraceae and Diversisporales.

Fungi form diverse communities and play essential roles in many terrestrial ecosystems, yet there are methodological challenges in taxonomic and phylogenetic placement of fungi from environmental sequences. To address such challenges, we investigated spatiotemporal structure of a fungal community using soil metabarcoding with four different sequencing strategies: short-amplicon sequencing of the ITS2 region (300-400 bp) with Illumina MiSeq, Ion Torrent Ion S5 and PacBio RS II, all from the same PCR library, as well as long-amplicon sequencing of the full ITS and partial LSU regions (1200-1600 bp) with PacBio RS II. Resulting community structure and diversity depended more on statistical method than sequencing technology. The use of long-amplicon sequencing enables construction of a phylogenetic tree from metabarcoding reads, which facilitates taxonomic identification of sequences. However, long reads present issues for denoising algorithms in diverse communities. We present a solution that splits the reads into shorter homologous regions prior to denoising, and then reconstructs the full denoised reads. In the choice between short and long amplicons, we suggest a hybrid approach using short amplicons for sampling breadth and depth, and long amplicons to characterize the local species pool for improved identification and phylogenetic analyses.

Forests and woodlands in the West African Guineo-Sudanian transition zone contain many tree species that form symbiotic interactions with ectomycorrhizal (ECM) fungi. These fungi facilitate plant growth by increasing nutrient and water uptake and include many fruiting body-forming fungi, including some edible mushrooms. Despite their importance for ecosystem functioning and anthropogenic use, diversity and distribution of ECM fungi is severely under-documented in West Africa. We conducted a broad regional sampling across five West African countries using soil eDNA to characterize the ECM as well as the total soil fungal community in gallery forests and savanna woodlands dominated by ECM host tree species. We subsequently sequenced the entire ITS region and much of the LSU region to infer a phylogeny for all detected soil fungal species. Utilizing a long read sequencing approach allows for higher taxonomic resolution by using the full ITS region, while the highly conserved LSU gene allows for a more accurate higher-level assignment of species hypotheses, including species without ITS-based taxonomy assignments. We detect no overall difference in species richness between gallery forests and woodlands. However, additional gallery forest plots and more samples per plot would have been needed to firmly conclude this pattern. Based on both abundance and richness, species from the families Russulaceae and Inocybaceae dominate the ECM fungal soil communities across both vegetation types. The community structure of both total soil fungi and ECM fungi was significantly influenced by vegetation types and showed strong correlation within plots. However, we found no significant difference in fungal community structure between samples collected adjacent to different host tree species within each plot. We conclude that within plots, the fungal community is structured more by the overall ECM host plant community than by the species of the individual host tree that each sample was collected from.