Bioinformatics matters: The accuracy of plant and soil fungal community data is highly dependent on the metabarcoding pipeline

Mycorrhizal Fungi Influence on Mature Tree Growth: Stronger in High-Nitrogen Soils for an EMF-Associated Tree and in Low-Nitrogen Soils for Two AMF-Associated Trees

The plant-mycorrhizal fungi relationship can range from mutualistic to parasitic as a function of the fungal taxa involved, plant ontogeny, as well as the availability of resources. Despite the implications this relationship may have on forest carbon cycling and storage, we know little about how mature trees may be impacted by mycorrhizae and how this impact may vary across the landscape. We collected growth data of two arbuscular mycorrhizal fungi (AMF)-associated tree species, Acer rubrum and A. saccharum, and one ectomycorrhizal fungi (EMF)-associated tree species, Quercus rubra, to assess how the mycorrhizal fungi-plant association may vary along a gradient of nitrogen (N) availability. Individual assessments of fungal taxa relative abundances showed non-linear associations with tree growth; positive associations for the two AMF-associated trees were mostly under low N, whereas positive to neutral associations for the EMF-associated tree mainly took place at high N. Only A. rubrum exhibited greater tree growth with its tree soil-specific mycorrhizal community when compared with predictions under a random mycorrhizal soil community. Because mycorrhizal fungi are likely to mediate how plants respond to warming, increasing levels of N deposition and of atmospheric CO2, understanding these relationships is critical to accurately forecasting tree growth.

Marine Mycobiomes Colonize Mediterranean Sponge Hosts in a Random Fashion

Marine sponges are widespread, sessile, filter-feeding animals, known for living in association with complex prokaryotic communities structured by host species. Though marine fungi are ubiquitous across marine environments, little is known about sponge-associated fungal communities (mycobiome). Indeed, aside from a few studies based on the isolation of fungal strains for biotechnological purposes, little information is available to understand the diversity and structure of sponge mycobiome. Here, a metabarcoding approach based on the ITS1 marker was applied to examine the structure and composition of fungal communities associated with four Mediterranean sponges. The species: Petrosia ficiformis, Chondrosia reniformis, Crambe crambe, and Chondrilla nucula were analyzed along with the surrounding seawater, revealing Aspergillus (1-56%), Cladosporium (1-75%), Malassezia (1-38.5%), and Pennicillium (1.5-36%) as the most represented fungal genera. Our data showed high intra-specific variability and no clear core mycobiome within each of the sponge species host, suggesting stochastic and perhaps transient community membership. This study sheds light on one of the most abundant yet least understood components of the marine ecosystem. Unraveling the dynamics of fungal interactions within sponge holobionts is essential to advance our understanding of their ecological roles and functions. By addressing the enigmatic nature of sponge-associated fungi, this research opens new avenues for exploring their contributions to marine ecosystems and resolving the many unanswered questions in this field.

Protein Coding Low-Copy rpb2 and ef1-α Regions Are Viable Fungal Metabarcoding DNA Markers Which Can Supplement ITS for Better Accuracy

The nuclear ribosomal DNA Internal Transcribed Spacer (ITS) region is used as a universal fungal barcode marker, but often lacks a significant DNA barcoding gap between sister taxa. Here we tested the reliability of protein coding low-copy genes as alternative barcode markers. Mock communities of three unrelated agaric genera (Dermoloma, Hodophilus, and Russula) representing lineages of closely related species were sequenced by the Illumina platform targeting the ITS1, ITS2, the second largest subunit of RNA polymerase II gene (rpb2) and the transcription elongation factor 1-alpha gene (ef1-α) regions. Species representation and their relative abundances were similar across all tested barcode regions, despite a lower copy number in protein coding markers. ITS1 and ITS2 required more sophisticated sequence filtering because they produced a high number of chimeric sequences requiring reference-based chimera removal and had a higher number of sequence variants per species. Although clustering of filtered ITS sequences resulted in an average higher number of correctly clustered units at optimal similarity thresholds, these thresholds varied substantially among genera. Best-fitted thresholds of low-copy markers were more consistent across genera but frequently lacked species resolution due to low intraspecific variability. At some thresholds, we observed multiple species lumped together, and at the same time, species split into multiple partial clusters, which should be taken into consideration when assessing the best clustering thresholds and taxonomic identity of clusters. To achieve the best taxonomic resolution and improve species detection, we recommend combining different markers and applying additional reference-based sorting of clusters. The current availability of rpb2 and ef1-α reference sequences in public databases is far from being complete for all fungal groups, but a combined marker approach can be used for group-specific studies that can build reference data for their own purposes.

Fungi promote cross-domain interactions even in deep anoxic mangrove sediments

BACKGROUND

Microbial communities in mangrove sediments play vital ecological roles that underpin the functioning of the overall mangrove ecosystem. Fungal communities, in particular, are known to play crucial roles across sediment systems, yet their roles in mangrove sediments, especially in deeper layers, remain poorly understood without a comprehensive inter-domain characterization. To better understand fungal roles in sediment horizons, 10 sediment cores extending down to a depth of 1 m were taken in three mangrove sites to characterise the archaeal, bacterial, and fungal communities at 10 cm depth intervals.

RESULTS

We demonstrate that sediment depth has distinct effects on the three microbial communities. While fungal community compositions were similar across sediment depths, bacterial and archaeal community compositions were stratified into three distinct layers, surface (10-30 cm), subsurface (40-60 cm), and deep (70-100 cm). Co-occurrence networks were then constructed to investigate the roles of fungi in these sediment layers, where fungi were consistently identified as keystone taxa in maintaining the microbial network topology, with co-domain interactions constituting more than half of all interactions. Even in the deepest layer, fungal nodes still retained high betweenness centralities, acting as network hubs to potentially augment microbial interactions vital for the functioning of the overall ecosystem.

CONCLUSIONS

Overall, our results emphasise the important role of fungi in mediating microbial interactions across sediment depths even in deep, anoxic sediment layers, and highlight the importance of cross-domain interactions as integral to a more holistic understanding of the mangrove microbiome.

Watershed urbanization alters aquatic plant mycobiomes through the loss of rare taxa

Urban expansion, projected to triple globally from 2000 to 2030, significantly impacts biodiversity and ecosystem processes, including those of microbial communities. Microbes are key drivers of many ecosystem processes and affect the fitness and resilience of plants and animals, but research on the biotic effects of urbanization has focused primarily on macroorganisms. This study investigates host-associated fungal communities in the pollution-tolerant aquatic plant Ranunculus aquatilis along an urbanization gradient in the Provo River, Utah, USA, a rapidly urbanizing region. We collected plant and adjacent water samples from 10 locations along the river, spanning from rural to urbanized areas within a single watershed, and conducted DNA amplicon sequencing to characterize fungal community composition. Our results show a significant decline in fungal alpha diversity correlated with increased urbanization metrics such as impervious surface area and developed land cover. Specifically, fungal richness and Shannon diversity decreased as urbanization intensified, driven primarily by a reduction in rare taxa. Despite a stable core microbiome dominated by a few taxa, the overall community structure varied significantly along the urbanization gradient, with notable shifts in dominant fungal taxa. Contrary to expectations, no detectable levels of heavy metals were found in water samples at any location, suggesting that other urbanization-related factors, potentially including organic pollutants or plant stress responses, influence fungal endophyte communities. Our findings underscore the need for further investigation into the mechanisms driving these patterns, particularly the roles of organic pollution, nutrient loads, and plant stress. As global urbanized watershed area grows, the fate of aquatic plant life is tied to their fungal community. Understanding these interactions is crucial for predicting the impacts of continued urbanization on freshwater ecosystems.

Application of Amplicon Metagenomics to Identify Fungal Pathogens in Formalin-Fixed Paraffin-Embedded Samples: Proof of Concept in Animals with Fungal Pathologies

The identification of fungal pathogens in formalin-fixed paraffin-embedded (FFPE) tissues is an unmet need in human and animal medicine, and sequence-agnostic approaches are needed to identify emerging pathogens. Eleven FFPE biopsy specimens with etiologic diagnoses of fungal disease based on standard testing of paired fresh tissue samples were utilized here to evaluate metabarcoding approaches. The cases included tissues from three dogs, three cats, one box turtle, one goat, one common loon, and one gray tree frog. The diagnoses from the fresh tissues in these cases were Microsporum canis, Penicillium sp., Exophiala sp. (likely E. jeanselmei), Verticillium sp., Rhizopus sp., atypical Cryptococcus neoformans, Conidiobolus spp., Aspergillus fumigatus, Cryptococcus neoformans var grubii, Batrachochytrium dendrobatidis, Fusarium solani, Blastomyces dermatitidis, Coccidiodes immitis, and Histoplasma capsulatum. We compared the ITS1 and 28S D1 rRNA gene genetic markers in combination with several bioinformatic strategies to identify fungal pathogens in the FFPE tissue samples, with a success rate of 9/11. These methods could allow diagnosticians who receive only FFPE tissues and see fungal pathogens to speciate the pathogens and could be of value in retrospective studies wherein FFPE tissue is the only archived tissue. Furthermore, these techniques could be of use to researchers investigating polymicrobial communities where DNA preservation is suboptimal.

The Establishment of a Terrestrial Macroalga Canopy Impacts Microbial Soil Communities in Antarctica

Prasiola is a genus of foliose green algae that forms extensive cryptogamic canopies that contribute to the greening of ice-free areas in the Antarctic tundra. To better understand the impact of Prasiola canopy establishment on colonization in these areas, this study compared the taxonomic and functional structures of bacterial and fungal communities in adjacent soils with and without extensive Prasiola colonization. DNA metabarcoding was employed to analyze the microbial community structure in these soils and in the canopy. Additionally, a phylogenetic study of Prasiola samples was conducted to characterize the taxonomic composition of the analyzed canopies, revealing the presence of Prasiola crispa (Lightfoot) Kützing and P. antarctica Kützing. Key soil attributes were assessed to examine the canopy's influence. Higher pH and carbon, nitrogen, and organic matter contents were found in Prasiola-covered soils than in bare soils. Furthermore, Prasiola canopy establishment not only influenced abiotic soil properties but also shaped soil microbial community structure and its functions. For instance, while Actinobacteriota predominated in bacterial communities both within the Prasiola canopy and beneath it, Bacteroidota dominated in the bare soil. Despite significant variability across soil types, fungal communities showed a trend of higher abundances in certain Ascomycetes, such as Helotiales, Hypocreales, or Xylariales, in soils beneath Prasiola compared to bare soils. Regarding functional diversity, covered soils exhibited a statistically significant lower potential for bacterial methanogenesis and autotrophic CO2 fixation compared to bare soils. Finally, lichenized fungi, plant pathogens, and fungal wood saprotrophs tended to be more abundant in covered soils.

Cover Cropping Attenuates Population Growth of Macrophomina phaseolina by Limiting Weed Biomass, Despite Asymptomatic Colonization of Cover Crops

Macrophomina phaseolina is a fungus that causes charcoal rot in strawberry and a wide variety of crop species. Little is known about its potential to asymptomatically colonize crop plants or grow saprophytically on their tissues, both of which would create a potential for alternate, asymptomatic hosts to lead to increases in inoculum. To test the impact of cover cropping on M. phaseolina abundance, we conducted randomized-block field experiments in soils infested by M. phaseolina. None of the 15 cover crop varieties showed symptoms of charcoal rot. All Fabaceae and Brassicaceae varieties were asymptomatically colonized at varying rates, but among Poaceae M. phaseolina was recovered from only one individual oat plant. Soil samples collected at the time of planting, tillage, and 8 weeks after tillage showed that cover cropping attenuated the growth of M. phaseolina relative to fallow plots harboring the weedy legume Medicago polymorpha. This weed species was abundantly colonized by this pathogen in both living root samples and plant residue collected 8 weeks after tillage. Cover cropping also influenced the diversity and composition of bulk soil bacterial and fungal communities, but these effects were not associated with M. phaseolina population density. Although M. phaseolina was not detected in living wheat tissues, it was recovered from wheat residue, suggesting that it may be facultatively saprophytic. These results suggest that cover cropping does not pose a risk for increasing disease caused by M. phaseolina and could be beneficial as conducive weed species, such as M. polymorpha, are suppressed.[Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Symbiotic fungi alter plant resource allocation independent of water availability

PREMISE

The ability of plants to adapt or acclimate to climate change is inherently linked to their interactions with symbiotic microbes, notably fungi. However, it is unclear whether fungal symbionts from different climates have different impacts on the outcome of plant-fungal interactions, especially under environmental stress.

METHODS

We tested three provenances of fungal inoculum (originating from dry, moderate or wet environments) with one host plant genotype exposed to three soil moisture regimes (low, moderate and high). Inoculated and uninoculated plants were grown in controlled conditions for 151 days, then shoot and root biomass were weighed and fungal diversity and community composition determined via amplicon sequencing.

RESULTS

The source of inoculum and water regime elicited significant changes in plant resource allocation to shoots versus roots, but only specific inocula affected total plant biomass. Shoot biomass increased in the high water treatment but was negatively impacted by all inoculum treatments relative to the controls. The opposite was true for roots, where the low water treatment led to greater proportional root biomass, and plants inoculated with wet site fungi allocated significantly more resources to root growth than dry- or moderate-site inoculated plants and the controls. Fungal communities of shoots and roots partitioned by inoculum source, water treatment, and the interaction of the two.

CONCLUSIONS

The provenance of fungi can significantly affect total plant biomass and resource allocation above- and belowground, with fungi derived from more extreme environments eliciting the strongest plant responses.

Fungal degradation of complex organic carbon supports denitrification in saturated woodchip bioreactors

Woodchip bioreactor (WBR) is a promising technology for the removal of nitrate from agricultural drainage, although the performance of WBRs is dependent on the decomposition of lignocellulosic biomass and the carbon availability for microbial denitrification. Fungal species are more efficient than bacterial counterparts in driving wood decomposition; however, little is known about the fungal community structure and functions in saturated WBRs. In this study, we investigated the dynamics of the mycobiome in field-scale, constantly saturated WBRs located in Willmar, Minnesota, USA. Fungal community analysis suggested that wood-rotting fungi were abundant in WBRs, especially near their inlet locations where microbial denitrification was most active. Complex network structures of fungal hyphae associated with a decayed cavity on the woodchip surface was further evidenced by confocal and scanning electron microscopy. These results suggest that fungi play a major role in wood degradation in WBRs, thereby promoting denitrification activity.

Untrimmed ITS2 metabarcode sequences cause artificially reduced abundances of specific fungal taxa

Advances in DNA metabarcoding have greatly expanded our knowledge of microbial communities in recent years. Pipelines and parameters have been tested extensively for bacterial metabarcoding using the 16S rRNA gene and best practices are largely established. For fungal metabarcoding using the internal transcribed spacer (ITS) gene, however, only a few studies have considered how such pipelines and parameters can affect community prediction. Here, we report a novel bias uncovered during ITS region 2 (ITS2) sequencing of Trichoderma-infected ant fungus gardens and confirmed this bias using mock communities. Abnormally low forward read quality caused Trichoderma ITS2 reads to be computationally filtered before and during read pair merging, thus almost entirely eliminating Trichoderma amplicon sequence variants from the resulting fungal community profiles. Sliding window quality trimming before filtering allowed most of these reads to pass filtering and merge successfully, producing community profiles that now correlated with visual signs of Trichoderma infection and matched the composition of the mock communities. Applying such sliding window trimming to a previously generated environmental ITS2 data set increased the detected fungal diversity and again overcame read quality biases against Trichoderma to detect it in nearly every sample instead and often at high relative abundances. This analysis additionally identified a similar, but distinct, bias against a second fungal genus Meyerozyma. The prevalence of such quality biases against other fungal ITS sequences is unknown but may be widespread. We, therefore, advocate for the routine use of sliding window quality trimming as a best practice in ITS2 metabarcoding analysis.

IMPORTANCE

Metabarcode sequencing produces DNA abundance profiles that are presumed to reflect the actual microbial composition of their corresponding input samples. However, this assumption is not always tested, and taxon-specific biases are often not apparent, especially for low-abundance taxa in complex communities. Here, we identified internal transcribed spacer region 2 (ITS2) read quality aberrations that caused dramatic reductions in the relative abundances of specific taxa in multiple data sets characterizing ant fungus gardens. Such taxon-specific biases in read quality may be widespread in other environments and for other fungal taxa, thereby causing incorrect descriptions of these mycobiomes.

Microbial changing patterns across lateral and vertical horizons in recently formed permafrost after the outburst of Zonag Lake, Tibetan Plateau

In polar and alpine regions, global warming and landform changes are draining lakes, transforming them into permafrost with altered microbial communities and element cycling. In this study, we investigated bacterial and archaeal (prokaryotic) community changes in the newly exposed sediment of Zonag Lake (Tibetan Plateau), focusing on prokaryotic diversity, community structure, and genes involved in carbon fixation and nitrogen cycling across lateral (up to 800 m) and vertical (up to 80 cm) horizons. The results showed that prokaryotic richness decreased across the lateral horizons, coinciding with reductions in carbon concentrations. Dramatic changes in community structure were also observed, primarily influenced by the distance from the lake and then by sediment depth, with environmental filtering and dispersal limitations shaping the lateral and vertical distributions, respectively. Based on PICRUSt2 results, the relative abundance of genes related to carbon fixation increased along the lateral horizon, suggesting that microbial carbon fixers are counteracting the carbon loss during permafrost formation. In contrast, the genes related to denitrification also increased, which may lead to nitrogen loss and contribute to global warming by releasing nitric oxide gas. This study highlights the resilience of prokaryotic communities in drained lake basins and their ecological implications under global warming.

Airborne DNA: State of the art - Established methods and missing pieces in the molecular genetic detection of airborne microorganisms, viruses and plant particles

Bioaerosol is composed of different particles, originating from organisms, or their fragments with different origin, shape, and size. Sampling, analysing, identification and describing this airborne diversity has been carried out for over 100 years, and more recently the use of molecular genetic tools has been implemented. However, up to now there are no established protocols or standards for detecting airborne diversity of bacteria, fungi, viruses, pollen, and plant particles. In this review we evaluated commonalities of methods used in molecular genetic based studies in the last 23 years, to give an overview of applicable methods as well as knowledge gaps in diversity assessment. Various sampling techniques show different levels of effectiveness in detecting airborne particles based on their DNA. The storage and processing of samples, as well as DNA processing, influences the outcome of sampling campaigns. Moreover, the decisions on barcode selection, method of analysis, reference database as well as negative and positive controls may severely impact the results obtained. To date, the chain of decisions, methodological biases and error propagation have hindered DNA based molecular sequencing from offering a holistic picture of the airborne biodiversity. Reviewing the available studies, revealed a great diversity in used methodology and many publications didn't state all used methods in detail, making comparisons with other studies difficult or impossible. To overcome these limitations and ensure genuine comparability across studies, it is crucial to standardize protocols. Publications need to include all necessary information to enable comparison among different studies and to evaluate how methodological choices can impacts the results. Besides standardization, implementing of automatic tools and combining of different analytical techniques, such as real-time evaluation combined with sampling and molecular genetic analysis, could assist in achieving the goal of accurately assessing the actual airborne biodiversity.

Differential response of subterranean microbiome to exogenous organic matter input in a cave ecosystem

As a recurrent climatic phenomenon in the context of climate change, extreme rainstorms induce vertical translocation of organic matter and increase moisture content in terrestrial ecosystems. However, it remains unclear whether heavy rainstorms can impact microbial communities in the deep biosphere by modulating organic matter input. In this study, we present findings on the different responses of bacterial and fungal communities in a subsurface cave to rainstorms and moisture variations through field surveys and microcosm experiments. During periods of rainstorms, the influx of dissolved organic matter (DOM) from soil overlying the cave into cave sediments significantly enhanced the correlation between core bacteria and environmental factors, particularly fluorescence spectral indices. The resource utilization of core bacteria was diminished, while the functional diversity of core fungi remained relatively unaltered. We also performed simulated experiments with restricted external DOM inputs, in which DOM content was observed to decrease and microbial diversity increase in response to artificially increased moisture content (MC). The niche breadth of core bacteria decreased and became more closely associated with DOM as the MC increased, while the niche breadth of core fungi remained predominantly unchanged. Compared to fungi, cave bacteria exhibited higher sensitivity towards variations in DOM. The core microbiome can efficiently utilize the available organic matter and participate in nitrogen- and sulfur-related metabolic processes. The study systematically revealed distinct microbial responses to rainstorm events, thereby providing valuable insights for future investigations into energy utilization within deep biospheres.

Comparison of commonly used software pipelines for analyzing fungal metabarcoding data

BACKGROUND

Metabarcoding targeting the internal transcribed spacer (ITS) region is commonly used to characterize fungal communities of various environments. Given their size and complexity, raw ITS sequences are necessarily processed and quality-filtered with bioinformatic pipelines. However, such pipelines are not yet standardized, especially for fungal communities, and those available may produce contrasting results. While some pipelines cluster sequences based on a specified percentage of base pair similarity into operational taxonomic units (OTUs), others utilize denoising techniques to infer amplicon sequencing variants (ASVs). While ASVs are now considered a more accurate representation of taxonomic diversity for prokaryote communities based on 16S rRNA amplicon sequencing, the applicability of this method for fungal ITS sequences is still debated.

RESULTS

Here we compared the performance of two commonly used pipelines DADA2 (inferring ASVs) and mothur (clustering OTUs) on fungal metabarcoding sequences originating from two different environmental sample types (fresh bovine feces and pasture soil). At a 99% OTU similarity threshold, mothur consistently identified a higher fungal richness compared to DADA2. In addition, mothur generated homogenous relative abundances across multiple technical replicates (n = 18), while DADA2 results for the same replicates were highly heterogeneous.

CONCLUSIONS

Our study highlights a potential pipeline-associated bias in fungal metabarcoding data analysis of environmental samples. Based on the homogeneity of relative abundances across replicates and the capacity to detect OTUs/ASVs, we suggest using OTU clustering with a similarity of 97% as the most appropriate option for processing fungal metabarcoding data.

Microbiomes of urban trees: unveiling contributions to atmospheric pollution mitigation

Urban trees are crucial in delivering essential ecosystem services, including air pollution mitigation. This service is influenced by plant associated microbiomes, which can degrade hydrocarbons, support tree health, and influence ecological processes. Yet, our understanding of tree microbiomes remains limited, thus affecting our ability to assess and quantify the ecosystem services provided by trees as complex systems. The main hypothesis of this work was that tree microbiomes concur to hydrocarbon biodegradation, and was tested through three case studies, which collectively investigated two tree micro-habitats (phyllosphere and tree cavity organic soil-TCOS) under various conditions representing diverse ecological scenarios, by applying different culture-based and molecular techniques and at different scales. The integration of all results provided a more comprehensive understanding of the role of microbiomes in urban trees. Firstly, bacterial strains isolated from the phyllosphere of Quercus ilex were characterized, indicating the presence of Plant-Growth Promoting bacteria and strains able to catabolize PAHs, particularly naphthalene and phenanthrene. Secondly, naphthalene biodegradation on artificially spiked Hedera helix leaves was quantified in greenhouse experiments on inoculated and untreated plants. The persistence of the inoculated strain and community structure of epiphytic bacteria were assessed by Illumina sequencing of V5-V6 hypervariable regions of 16S rRNA gene. Results showed that naphthalene degradation was initially faster on inoculated plants but later the degradation rates became similar, probably because bacterial populations with hydrocarbon-degrading abilities gradually developed also on non-inoculated plants. Finally, we explored bacterial and fungal biodiversity hosted by TCOS samples, collected from six large trees located in an urban park and belonging to different species. Microbial communities were characterized by Illumina sequencing of V5-V6 hypervariable regions of bacterial gene 16S rRNA and of fungal ITS1. Results indicated TCOS as a distinct substrate, whose microbiome is determined both by the host tree and by canopy environmental conditions and has a pronounced aerobic hydrocarbon degradation potential. Overall, a better assessment of biodiversity associated with trees and the subsequent provision of ecosystem services constitute a first step toward developing future new microbe-driven sustainable solutions, especially in terms of support for urban green planning and management policy.

New insights into biologic interpretation of bioinformatic pipelines for fish eDNA metabarcoding: A case study in Pearl River estuary

Environmental DNA (eDNA) metabarcoding is an emerging tool for monitoring biological communities in aquatic ecosystems. The selection of bioinformatic pipelines significantly impacts the results of biodiversity assessments. However, there is currently no consensus on the appropriate bioinformatic pipelines for fish community analysis in eDNA metabarcoding. In this study, we compared three bioinformatic pipelines (Uparse, DADA2, and UNOISE3) using real and mock (constructed with 15/30 known fish) communities to investigate the differences in biological interpretation during the data analysis process in eDNA metabarcoding. Performance evaluation and diversity analyses revealed that the choice of bioinformatic pipeline could impact the biological results of metabarcoding experiments. Among the three pipelines, the operational taxonomic units (OTU)-based pipeline (Uparse) showed the best performance (sensitivity: 0.6250 ± 0.0166; compositional similarity: 0.4000 ± 0.0571), the highest richness (25-102) and minimal inter-group differences in alpha diversity. It suggested the OTU-based pipeline possessed superior capability in fish diversity monitoring compared to ASV/ZOTU-based pipeline. Additionally, the Bray-Curtis distance matrix achieved the highest discriminative effect in the PCoA (43.3%-53.89%) and inter-group analysis (P < 0.01), indicating it was better at distinguishing compositional differences or specific genera of fish community at different sampling sites than other distance matrices. These findings provide new insights into fish community monitoring through eDNA metabarcoding in estuarine environments.

Shifts in mycorrhizal types of fungi and plants in response to fertilisation, warming and herbivory in a tundra grassland

Climate warming is severely affecting high-latitude regions. In the Arctic tundra, it may lead to enhanced soil nutrient availability and interact with simultaneous changes in grazing pressure. It is presently unknown how these concurrently occurring global change drivers affect the root-associated fungal communities, particularly mycorrhizal fungi, and whether changes coincide with shifts in plant mycorrhizal types. We investigated changes in root-associated fungal communities and mycorrhizal types of the plant community in a 10-yr factorial experiment with warming, fertilisation and grazing exclusion in a Finnish tundra grassland. The strongest determinant of the root-associated fungal community was fertilisation, which consistently increased potential plant pathogen abundance and had contrasting effects on the different mycorrhizal fungal types, contingent on other treatments. Plant mycorrhizal types went through pronounced shifts, with warming favouring ecto- and ericoid mycorrhiza but not under fertilisation and grazing exclusion. Combination of all treatments resulted in dominance by arbuscular mycorrhizal plants. However, shifts in plant mycorrhizal types vs fungi were mostly but not always aligned in their magnitude and direction. Our results show that our ability to predict shifts in symbiotic and antagonistic fungal communities depend on simultaneous consideration of multiple global change factors that jointly alter plant and fungal communities.

Arbuscular mycorrhizal diversity increases across a plant productivity gradient driven by soil nitrogen availability

Arbuscular mycorrhizal fungi (AMF) are widespread obligate symbionts of plants. This dynamic symbiosis plays a large role in successful plant performance, given that AMF help to ameliorate plant responses to abiotic and biotic stressors. Although the importance of this symbiosis is clear, less is known about what may be driving this symbiosis, the plant's need for nutrients or the excess of plant photosynthate being transferred to the AMF, information critical to assess the functionality of this relationship. Characterizing the AMF community along a natural plant productivity gradient is a first step in understanding how this symbiosis may vary across the landscape. We surveyed the AMF community diversity at 12 sites along a plant productivity gradient driven by soil nitrogen availability. We found that AMF diversity in soil environmental DNA significantly increased along with the growth of the host plants Acer rubrum and A. saccharum., a widespread tree genus. These increases also coincided with a natural soil inorganic N availability gradient. We hypothesize photosynthate from the increased tree growth is being allocated to the belowground AMF community, leading to an increase in diversity. These findings contribute to understanding this complex symbiosis through the lens of AMF turnover and suggest that a more diverse AMF community is associated with increased host-plant performance.

Measurements of soil protist richness and community composition are influenced by primer pair, annealing temperature, and bioinformatics choices

Protists are a diverse and understudied group of microbial eukaryotic organisms especially in terrestrial environments. Advances in molecular methods are increasing our understanding of the distribution and functions of these creatures; however, there is a vast array of choices researchers make including barcoding genes, primer pairs, PCR settings, and bioinformatic options that can impact the outcome of protist community surveys. Here, we tested four commonly used primer pairs targeting the V4 and V9 regions of the 18S rRNA gene using different PCR annealing temperatures and processed the sequences with different bioinformatic parameters in 10 diverse soils to evaluate how primer pair, amplification parameters, and bioinformatic choices influence the composition and richness of protist and non-protist taxa using Illumina sequencing. Our results showed that annealing temperature influenced sequencing depth and protist taxon richness for most primer pairs, and that merging forward and reverse sequencing reads for the V4 primer pairs dramatically reduced the number of sequences and taxon richness of protists. The data sets of primers that targeted the same 18S rRNA gene region (e.g., V4 or V9) had similar protist community compositions; however, data sets from primers targeting the V4 18S rRNA gene region detected a greater number of protist taxa compared to those prepared with primers targeting the V9 18S rRNA region. There was limited overlap of protist taxa between data sets targeting the two different gene regions (80/549 taxa). Together, we show that laboratory and bioinformatic choices can substantially affect the results and conclusions about protist diversity and community composition using metabarcoding.IMPORTANCEEcosystem functioning is driven by the activity and interactions of the microbial community, in both aquatic and terrestrial environments. Protists are a group of highly diverse, mostly unicellular microbes whose identity and roles in terrestrial ecosystem ecology have been largely ignored until recently. This study highlights the importance of choices researchers make, such as primer pair, on the results and conclusions about protist diversity and community composition in soils. In order to better understand the roles protist taxa play in terrestrial ecosystems, biases in methodological and analytical choices should be understood and acknowledged.

A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

Arthropods as Vectors of Grapevine Trunk Disease Pathogens: Quantification of Phaeomoniella chlamydospora on Arthropods and Mycobiome Analysis of Earwig Exoskeletons

Viticulture worldwide is challenged by grapevine trunk diseases (GTDs). Involvement of arthropods in the dissemination process of GTD pathogens, notably esca pathogens, is indicated after detection of associated pathogens on arthropod exoskeletons, and demonstration of transmission under artificial conditions. The present study is the first to quantify spore loads via qPCR of the esca-relevant pathogen Phaeomoniella chlamydospora on arthropods collected in German vineyards, i.e., European earwigs (Forficula auricularia), ants (Formicidae), and two species of jumping spiders (Marpissa muscosa and Synageles venator). Quantification of spore loads showed acquisition on exoskeletons, but most arthropods carried only low amounts. The mycobiome on earwig exoskeletons was described for the first time to reveal involvement of earwigs in the dispersal of GTDs in general. Metabarcoding data support the potential risk of earwigs as vectors for predominantly Pa. chlamydospora and possibly Eutypa lata (causative agent of Eutypa dieback), as respective operational taxonomical unit (OTU) assigned genera had relative abundances of 6.6% and 2.8% in total reads, even though with great variation between samples. Seven further GTD-related genera were present at a very low level. As various factors influence the successful transmission of GTD pathogens, we hypothesize that arthropods might irregularly act as direct vectors. Our results highlight the importance of minimizing and protecting pruning wounds in the field.

The Effects of Salinity and Genotype on the Rhizospheric Mycobiomes in Date Palm Seedlings

Salinity severely affects the health and productivity of plants, with root-associated microbes, including fungi, potentially playing a crucial role in mitigating this effect and promoting plant health. This study employed metagenomics to investigate differences in the structures of the epiphyte mycobiomes in the rhizospheres of seedlings of two distinct date palm cultivars with contrasting salinity tolerances, the susceptible cultivar, 'Zabad', and the tolerant cultivar, 'Umsila'. Next-generation sequencing (NGS) of the internal transcribed spacer (ITS) rRNA was utilized as a DNA barcoding tool. The sequencing of 12 mycobiome libraries yielded 905,198 raw sequences of 268,829 high-quality reads that coded for 135 unique and annotatable operational taxonomic units (OTUs). An OTU analysis revealed differences in the rhizofungal community structures between the treatments regardless of genotype, and non-metric dimensional scaling (N-MDS) analyses demonstrated distinct separations between the cultivars under saline stress. However, these differences were not detected under the control environmental conditions, i.e., no salinity. The rhizospheric fungal community included four phyla (Ascomycota, Basidiomycota, Chytridiomycota, and Mucoromycota), with differences in the abundances of Aspergillus, Clonostachys, and Fusarium genera in response to salinity, regardless of the genotype. Differential pairwise comparisons showed that Fusarium falciforme-solani and Aspergillus sydowii-versicolor increased in abundance under saline conditions, providing potential future in vitro isolation guidelines for plant growth-promoting fungi. This study highlights the intricate dynamics of the rhizosphere microbial communities in date palms and their responses to salt stress. Additionally, we found no support for the hypothesis that indigenous epiphytic fungal communities are significantly involved in salinity tolerance in date palms.

Primer pairs, PCR conditions, and peptide nucleic acid clamps affect fungal diversity assessment from plant root tissues

High-throughput sequencing has become a prominent tool to assess plant-associated microbial diversity. Still, some technical challenges remain in characterising these communities, notably due to plant and fungal DNA co-amplification. Fungal-specific primers, Peptide Nucleic Acid (PNA) clamps, or adjusting PCR conditions are approaches to limit plant DNA contamination. However, a systematic comparison of these factors and their interactions, which could limit plant DNA contamination in the study of plant mycobiota, is still lacking. Here, three primers targeting the ITS2 region were evaluated alone or in combination with PNA clamps both on nettle (Urtica dioica) root DNA and a mock community. PNA clamps did not improve the richness or diversity of the fungal communities but increased the number of fungal reads. Among the tested factors, the most significant was the primer pair. Specifically, the 5.8S-Fun/ITS4-Fun pair exhibited a higher OTU richness but fewer fungal reads. Our study demonstrates that the choice of primers is critical for limiting plant and fungal DNA co-amplification. PNA clamps increase the number of fungal reads when ITS2 is targeted but do not result in higher fungal diversity recovery at high sequencing depth. At lower read depths, PNA clamps might enhance microbial diversity quantification for primer pairs lacking fungal specificity.

Microplastic fibres affect soil fungal communities depending on drought conditions with consequences for ecosystem functions

Microplastics affect soil functions depending on drought conditions. However, how their combined effect influences soil fungi and their linkages with ecosystem functions is still unknown. To address this, we used rhizosphere soil from a previous experiment in which we employed microplastic fibres addition and drought in a factorial design, and evaluated their effects on soil fungal communities. Microplastics decreased soil fungal richness under well-watered conditions, likely linked to microplastics leaching toxic substances into the soil, and microplastic effects on root fineness. Under drought, by contrast, microplastics increased pathogen and total fungal richness, likely related to microplastic positive effects on soil properties, such as water holding capacity, porosity or aggregation. Soil fungal richness was the attribute most affected by microplastics and drought. Microplastics altered the relationships between soil fungi and ecosystem functions to the point that many of them flipped from positive to negative or disappeared. The combined effect of microplastics and drought on fungal richness mitigated their individual negative effect (antagonism), suggesting that changes in soil water conditions may alter the action mode of microplastics in soil. Microplastic leaching of harmful substances can be mitigated under drought, while the improvement of soil properties by microplastics may alleviate such drought conditions.

Mock microbial community meta-analysis using different trimming of amplicon read lengths

Trimming of sequencing reads is a pre-processing step that aims to discard sequence segments such as primers, adapters and low quality nucleotides that will interfere with clustering and classification steps. We evaluated the impact of trimming length of paired-end 16S and 18S rRNA amplicon reads on the ability to reconstruct the taxonomic composition and relative abundances of communities with a known composition in both even and uneven proportions. We found that maximizing read retention maximizes recall but reduces precision by increasing false positives. The presence of expected taxa was accurately predicted across broad trim length ranges but recovering original relative proportions remains a difficult challenge. We show that parameters that maximize taxonomic recovery do not simultaneously maximize relative abundance accuracy. Trim length represents one of several experimental parameters that have non-uniform impact across microbial clades, making it a difficult parameter to optimize. This study offers insights, guidelines, and helps researchers assess the significance of their decisions when trimming raw reads in a microbiome analysis based on overlapping or non-overlapping paired-end amplicons.

Forest swamp succession alters organic carbon composition and survival strategies of soil microbial communities

Forest swamp ecosystems plays crucial role in the global carbon cycle. However, the effects of forest swamp succession on soil organic matter (SOM) and microbial community structure remain unclear. To determine the drivers of SOM change and soil microbial communities in forest swamp succession, a 'space instead of time' approach was used. Soil samples from 0 to 40 cm were collected along forest swamp (early stage), dried-up forest swamp (middle stage), and forest (late stage) ecosystems. Our findings reveal that as succession progresses, the relative content of aromatics decreases and SOM undergoes a transition towards a more readily degradable form. These changes affect soil carbon sequestration and nutrient availability. Bacterial diversity was significantly influenced by succession and changes in soil depth, with fungi exhibiting higher resilience. Soil properties and environmental conditions exert influence over the structure and function of microorganisms. As succession occurred, microbial interactions shifted from cooperation to competition, with bacteria displaying a deterministic distribution pattern and fungi exhibiting a random distribution pattern. SOM quality plays a key role in shaping microbial communities and influencing their growth strategies. Microorganisms are the major drivers of soil respiration, with K-strategist dominated communities in early succession exhibiting slower degradation rates, whereas r-strategists dominated in later stages, leading to faster decomposition.

Peat loss collocates with a threshold in plant-mycorrhizal associations in drained peatlands encroached by trees

Drainage-induced encroachment by trees may have major effects on the carbon balance of northern peatlands, and responses of microbial communities are likely to play a central mechanistic role. We profiled the soil fungal community and estimated its genetic potential for the decay of lignin and phenolics (class II peroxidase potential) along peatland drainage gradients stretching from interior locations (undrained, open) to ditched locations (drained, forested). Mycorrhizal fungi dominated the community across the gradients. When moving towards ditches, the dominant type of mycorrhizal association abruptly shifted from ericoid mycorrhiza to ectomycorrhiza at c. 120 m from the ditches. This distance corresponded with increased peat loss, from which more than half may be attributed to oxidation. The ectomycorrhizal genus Cortinarius dominated at the drained end of the gradients and its relatively higher genetic potential to produce class II peroxidases (together with Mycena) was positively associated with peat humification and negatively with carbon-to-nitrogen ratio. Our study is consistent with a plant-soil feedback mechanism, driven by a shift in the mycorrhizal type of vegetation, that potentially mediates changes in aerobic decomposition during postdrainage succession. Such feedback may have long-term legacy effects upon postdrainage restoration efforts and implication for tree encroachment onto carbon-rich soils globally.

The gut microbiota composition is linked to subsequent occurrence of ventilator-associated pneumonia in critically ill patients

Ventilator-associated pneumonia (VAP) is the most frequent nosocomial infection in critically ill-ventilated patients. Oropharyngeal and lung microbiota have been demonstrated to be associated with VAP occurrence, but the involvement of gut microbiota has not been investigated so far. Therefore, the aim of this study is to compare the composition of the gut microbiota between patients who subsequently develop VAP and those who do not. A rectal swab was performed at admission of every consecutive patient into the intensive care unit (ICU) from October 2019 to March 2020. After DNA extraction, V3-V4 and internal transcribed spacer 2 regions deep-sequencing was performed on MiSeq sequencer (Illumina) and data were analyzed using Divisive Amplicon Denoising Algorithm 2 (DADA2) pipeline. Among 255 patients screened, 42 (16%) patients with invasive mechanical ventilation for more than 48 h were included, 18 (43%) with definite VAP and 24 without (57%). Patients who later developed VAP had similar gut bacteriobiota and mycobiota α-diversities compared to those who did not develop VAP. However, gut mycobiota was dissimilar (β-diversity) between these two groups. The presence of Megasphaera massiliensis was associated with the absence of VAP occurrence, whereas the presence of the fungal genus Alternaria sp. was associated with the occurrence of VAP. The composition of the gut microbiota, but not α-diversity, differs between critically ill patients who subsequently develop VAP and those who do not. This study encourages large multicenter cohort studies investigating the role of gut-lung axis and oropharyngeal colonization in the development of VAP in ICU patients. Trial registration number: NCT04131569, date of registration: 18 October 2019. IMPORTANCE The composition of the gut microbiota, but not α-diversity, differs between critically ill patients who subsequently develop ventilator-associated pneumonia (VAP) and those who do not. Investigating gut microbiota composition could help to tailor probiotics to provide protection against VAP.

Profiling Walnut Fungal Pathobiome Associated with Walnut Dieback Using Community-Targeted DNA Metabarcoding

Walnut dieback can be caused by several fungal pathogenic species, which are associated with symptoms ranging from branch dieback to fruit necrosis and blight, challenging the one pathogen-one disease concept. Therefore, an accurate and extensive description of the walnut fungal pathobiome is crucial. To this end, DNA metabarcoding represents a powerful approach provided that bioinformatic pipelines are evaluated to avoid misinterpretation. In this context, this study aimed to determine (i) the performance of five primer pairs targeting the ITS region in amplifying genera of interest and estimating their relative abundance based on mock communities and (ii) the degree of taxonomic resolution using phylogenetic trees. Furthermore, our pipelines were also applied to DNA sequences from symptomatic walnut husks and twigs. Overall, our results showed that the ITS2 region was a better barcode than ITS1 and ITS, resulting in significantly higher sensitivity and/or similarity of composition values. The ITS3/ITS4_KYO1 primer set allowed to cover a wider range of fungal diversity, compared to the other primer sets also targeting the ITS2 region, namely, GTAA and GTAAm. Adding an extraction step to the ITS2 sequence influenced both positively and negatively the taxonomic resolution at the genus and species level, depending on the primer pair considered. Taken together, these results suggested that Kyo set without ITS2 extraction was the best pipeline to assess the broadest fungal diversity, with a more accurate taxonomic assignment, in walnut organs with dieback symptoms.

Salt altered rhizosphere fungal community and induced soybean recruit specific species to ameliorate salt stress

Different crop genotypes showed different adaptability to salt stress, which is partly attributable to the microorganisms in the rhizosphere. Yet, knowledge about how fungal communities of different genotypes in soybean respond to salt stress is limited. Here, qPCR and ITS sequencing were used to assess the response of rhizobial fungal communities of resistant and susceptible soybean to salt stress. Moreover, we isolated two fungal species recruited by resistant soybeans for validation. The assembly of fungal community structure might be strongly linked to alterations in fungal abundance and soil physicochemical properties. Salt stress derived structural differences in fungal communities of resistant and susceptible genotypes. The salt-resistant genotype appeared to recruit some fungal taxa to the rhizosphere to help mitigating salt stress. An increase of fungal taxa with predicted saprotrophic lifestyles might help promoting plant growth by increasing nutrient availability to the plants. Compared with the susceptible genotypes, the resistant genotypes had more stronger network structure of fungi. Lastly, we verified that recruited fungi, such as Penicillium and Aspergillus, can soybean adapt to salt stress. This study provided a promising approach for rhizospheric fungal community to enhance salt tolerance of soybean from the perspective of microbiology and ecology.

Biodiversity of Basidiomycetous Yeasts Associated with Cladonia rei Lichen in Japan, with a Description of Microsporomyces cladoniophilus sp. nov

For more than a century, lichens have been used as an example of dual-partner symbiosis. Recently, this has been challenged by the discovery of various basidiomycetous yeasts that coexist in multiple lichen species, among which Cladonia lichens from Europe and the United States were discovered to be highly specifically associated with the basidiomycetous yeast of the family Microsporomycetaceae. To verify this highly specific relationship, we investigated the diversity of basidiomycetous yeasts associated with Cladonia rei, a widely distributed lichen in Japan, by applying two approaches: yeast isolation from the lichen thalli and meta-barcoding analysis. We obtained 42 cultures of Cystobasidiomycetous yeast which were grouped into six lineages within the family Microsporomycetaceae. Unexpectedly, although the cystobasidiomycetes-specific primer was used, not only the cystobasidiomycetous yeasts but species from other classes were also detected via the meta-barcoding dataset; in particular, pucciniomycetous yeasts were found at a high frequency in some samples. Further, Halobasidium xiangyangense, which was detected in every sample with high abundance, is highly likely a generalist epiphytic fungus that has the ability to associate with C. rei. In the pucciniomycetous group, most of the detected species belong to the scale insect-associated yeast Septobasidium genus. In conclusion, even though Microsporomyces species are not the only yeast group associated with Cladonia lichen, our study demonstrated that the thalli of Cladonia rei lichen could be a suitable habit for them.

Bridging gut microbiota composition with extended-spectrum beta-lactamase Enterobacteriales faecal carriage in critically ill patients (microbe cohort study)

BACKGROUND

The worldwide dissemination of extended spectrum beta-lactamase producing Enterobacteriales (ESBL-E) is of major concern. Microbiota may play a role in the host resistance to colonization with ESBL-E, but the underlying mechanisms remain unknown. We aimed to compare the gut microbiota composition between ESBL-producing E. coli or K. pneumoniae carriers and ESBL-E non-carriers according to the bacterial species.

RESULTS

Among 255 patients included, 11 (4,3%) were colonized with ESBL-producing E. coli and 6 (2,4%) with ESBL-producing K. pneumoniae, which were compared with age- and sex-matched ESBL-E non carriers. While no significant differences were found between ESBL-producing E. coli carriers and non-carriers, gut bacteriobiota α-diversity was decreased in ESBL-K. pneumoniae faecal carriers compared both with non-carriers (p = 0.05), and with ESBL-producing E. coli carriers. The presence of Sellimonas intestinalis was associated with the absence of ESBL-producing E. coli fecal carriage. Campylobacter ureolyticus, Campylobacter hominis, bacteria belonging to Clostridium cluster XI and Saccharomyces sp. were associated with the absence of ESBL-producing K. pneumoniae faecal carriage.

CONCLUSIONS

The composition of the gut microbiota differs between ESBL-producing E. coli and K. pneumoniae faecal carriers suggesting that microbial species should be taken into account when investigating the role of gut microbiota in resistance to gut colonization with ESBL-E.

TRIAL REGISTRATION NUMBER

NCT04131569, date of registration: October 18, 2019.

Effects of Lumacaftor-Ivacaftor on Airway Microbiota-Mycobiota and Inflammation in Patients with Cystic Fibrosis Appear To Be Linked to Pseudomonas aeruginosa Chronic Colonization

The management of cystic fibrosis has been transformed recently by the advent of CFTR modulators, including lumacaftor-ivacaftor. However, the effects of such therapies on the airway ecosystem, particularly on the microbiota-mycobiota and local inflammation, which are involved in the evolution of pulmonary damage, are unclear.

Metabarcoding Reveals Impact of Different Land Uses on Fungal Diversity in the South-Eastern Region of Antioquia, Colombia

Changes in soil fungal communities caused by land use have not been sufficiently studied in South American Andosols, which are considered key food production areas. Since fungal communities play an important role in soil functionality, this study analysed 26 soil samples of Andosols collected from locations devoted to conservation, agriculture and mining activities in Antioquia, Colombia, to establish differences between fungal communities as indicators of soil biodiversity loss using Illumina MiSeq metabarcoding on nuclear ribosomal ITS2 region. A non-metric multidimensional scaling allowed to explore driver factors of changes in fungal communities, while the significance of these variations was assessed by PERMANOVA. Furthermore, the effect size of land use over relevant taxa was quantified. Our results suggest a good coverage of fungal diversity with a detection of 353,312 high-quality ITS2 sequences. We found strong correlations of Shannon and Fisher indexes with dissimilarities on fungal communities (r = 0.94). These correlations allow grouping soil samples according to land use. Variations in temperature, air humidity and organic matter content lead to changes in abundances of relevant orders (Wallemiales and Trichosporonales). The study highlights specific sensitivities of fungal biodiversity features in tropical Andosols, which may serve as a basis for robust assessments of soil quality in the region.

A 17-year time-series of fungal environmental DNA from a coastal marine ecosystem reveals long-term seasonal-scale and inter-annual diversity patterns

Changing patterns in diversity are a feature of many habitats, with seasonality a major driver of ecosystem structure and function. In coastal marine plankton-based ecosystems, seasonality has been established through long-term time-series of bacterioplankton and protists. Alongside these groups, fungi also inhabit coastal marine ecosystems. If and how marine fungi show long-term intra- and inter-annual diversity patterns is unknown, preventing a comprehensive understanding of marine fungal ecology. Here, we use a 17-year environmental DNA time-series from the English Channel to determine long-term marine fungal diversity patterns. We show that fungal community structure progresses at seasonal and monthly scales and is only weakly related to environmental parameters. Communities restructured every 52-weeks suggesting long-term stability in diversity patterns. Some major marine fungal genera have clear inter-annual recurrence patterns, re-appearing in the annual cycle at the same period. Low relative abundance taxa that are likely non-marine show seasonal input to the coastal marine ecosystem suggesting land-sea exchange regularly takes place. Our results demonstrate long-term intra- and inter-annual marine fungal diversity patterns. We anticipate this study could form the basis for better understanding the ecology of marine fungi and how they fit in the structure and function of the wider coastal marine ecosystem.

VTAM: A robust pipeline for validating metabarcoding data using controls

To obtain accurate estimates for biodiversity and ecological studies, metabarcoding studies should be carefully designed to minimize both false positive (FP) and false negative (FN) occurrences. Internal controls (mock samples and negative controls), replicates, and overlapping markers allow controlling metabarcoding errors but current metabarcoding software packages do not explicitly integrate these additional experimental data to optimize filtering. We have developed the metabarcoding analysis software VTAM, which uses explicitly these elements of the experimental design to find optimal parameter settings that minimize FP and FN occurrences. VTAM showed similar sensitivity, but a higher precision compared to two other pipelines using three datasets and two different markers (COI, 16S). The stringent filtering procedure implemented in VTAM aims to produce robust metabarcoding data to obtain accurate ecological estimates and represents an important step towards a non-arbitrary and standardized validation of metabarcoding data for conducting ecological studies. VTAM is implemented in Python and available from: https://github.com/aitgon/vtam. The VTAM benchmark code is available from: https://github.com/aitgon/vtam_benchmark.

Extracting abundance information from DNA-based data

The accurate extraction of species-abundance information from DNA-based data (metabarcoding, metagenomics) could contribute usefully to diet analysis and food-web reconstruction, the inference of species interactions, the modelling of population dynamics and species distributions, the biomonitoring of environmental state and change, and the inference of false positives and negatives. However, multiple sources of bias and noise in sampling and processing combine to inject error into DNA-based data sets. To understand how to extract abundance information, it is useful to distinguish two concepts. (i) Within-sample across-species quantification describes relative species abundances in one sample. (ii) Across-sample within-species quantification describes how the abundance of each individual species varies from sample to sample, such as over a time series, an environmental gradient or different experimental treatments. First, we review the literature on methods to recover across-species abundance information (by removing what we call "species pipeline biases") and within-species abundance information (by removing what we call "pipeline noise"). We argue that many ecological questions can be answered with just within-species quantification, and we therefore demonstrate how to use a "DNA spike-in" to correct for pipeline noise and recover within-species abundance information. We also introduce a model-based estimator that can be used on data sets without a physical spike-in to approximate and correct for pipeline noise.

Characterization of Seed Mycobiota Using Culture-Dependent and Culture-Independent Approaches

Seed fungi are potentially important for their roles in seedling microbiome assembly and seedling health, but surveys of full seed fungal communities remain limited. While culture-dependent methods have been used to characterize some members of the seed mycobiota, recent culture-independent studies have improved the ease in identifying and characterizing full seed fungal communities. In this chapter, we describe how to survey seed fungi using both traditional culture-based methods and culture-free metabarcoding. We first describe protocols for the isolation and long-term preservation of fungal strains from individual seeds and for the extraction and amplification of DNA from such fungal isolates for identification with Sanger sequencing. We also detail how to extract, amplify, and sequence fungal DNA directly from individual seeds. Finally, we provide suggestions for troubleshooting media choices, PCR inhibition by isolates and plant tissue, and PCR limitation by low fungal DNA.

Microbiota Modulation in Blueberry Rhizosphere by Biocontrol Bacteria

Microbial interactions in agricultural soils can play important roles in the control of soil-borne phytopathogenic diseases. Yields from blueberry plantations from southern Spain have been impacted by the pathogenic fungus, Macrophomina phaseolina. The use of chemical fungicides has been the common method for preventing fungal infections, but due to their high environmental impact, legislation is increasingly restricting its use. Biocontrol alternatives based on the use of microorganisms is becoming increasingly important. Using the metabarcoding technique, fungi and bacteria were characterized (via 16S and ITS regions, respectively) from rhizosphere soils of healthy and dead blueberry plants infected by M. phaseolina, and which had undergone three different treatments: two biocontrol strategies—one of them a mix of Pseudomonas aeruginosa and Bacillus velezensis and the other one with Bacillus amyloliquefaciens—and a third treatment consisting of the application of a nutrient solution. The treatments produced changes in the bacterial microbiota and, to a lesser extent, in the fungi. The abundance of Fusarium was correlated with dead plants, likely favoring the infection by M. phaseolina. The presence of other microorganisms in the soil, such as the fungi Archaeorhizomyces or the bacteria Actinospica, were correlated with healthy plants and could promote their survival. The different genera detected between dead and healthy plants opens the possibility of studying new targets that can act against infection and identify potential microorganisms that can be used in biocontrol strategies.

Effects of short-term soil exposure of different doses of ZnO nanoparticles on the soil environment and the growth and nitrogen fixation of alfalfa

The extensive application of nanomaterials has increased their levels in soil environments. Therefore, clarifying the process of environmental migration is important for environmental safety and human health. In this study, alfalfa was used to determine the effects of different doses of ZnO nanoparticles (NPs) on the growth of alfalfa and the soil environment. Results showed that the alfalfa biomass was inversely proportional to the exposure concentration of ZnO NPs. The Zn concentration in the alfalfa tissue and the exposure dose presented a significant positive correlation. A high concentration of ZnO NPs decreased the nitrogen-fixing area of root nodules while the number of bacteroids and root nodules, which in turn affected the nitrogen-fixing ability of alfalfa. At the same time, it caused different degrees of damage to the root nodules and root tip cells of alfalfa. A high dose of ZnO NPs decreased the relative abundance and diversity of the soil microorganisms. Therefore, short-term and high-dose exposure of ZnO NPs causes multiple toxicities in plants and soil environments.

Advancing biological invasion hypothesis testing using functional diversity indices

Pioneering investigations on the effects of introduced populations on community structure, ecosystem functioning and services have focused on the effects of invaders on taxonomic diversity. However, taxonomic-based diversity metrics overlook the heterogeneity of species roles within and among communities. As the homogenizing effects of biological invasions on community and ecosystem processes can be subtle, they may require the use of functional diversity indices to be properly evidenced. Starting from the listing of major functional diversity indices, alongside the presentation of their strengths and limitations, we focus on studies pertaining to the effects of invasive species on native communities and recipient ecosystems using functional diversity indices. By doing so, we reveal that functional diversity of the recipient community may strongly vary at the onset of the invasion process, while it stabilizes at intermediate and high levels of invasion. As functional changes occurring during the lag phase of an invasion have been poorly investigated, we show that it is still unknown whether there are consistent changes in functional diversity metrics that could indicate the end of the lag phase. Thus, we recommend providing information on the invasion stage under consideration when computing functional diversity metrics. For the existing literature, it is also surprising that very few studies explored the functional difference between organisms from the recipient communities and invaders of the same trophic levels, or assessed the effects of non-native organism establishment into a non-analogue versus an analogue community. By providing valuable tools for obtaining in-depth diagnostics of community structure and functioning, functional diversity indices can be applied for timely implementation of restoration plans and improved conservation strategies. To conclude, our work provides a first synthetic guide for their use in hypothesis testing in invasion biology.

Ecological Observations Based on Functional Gene Sequencing Are Sensitive to the Amplicon Processing Method

Until recently, the de facto method for short-read-based amplicon reconstruction was a sequence similarity threshold approach (operational taxonomic units [OTUs]). This has changed with the amplicon sequence variant (ASV) method where distributions are fitted to abundance profiles of individual genes using a noise-error model. While OTU-based approaches are still useful for 16S rRNA/18S rRNA genes, where thresholds of 97% to 99% are used, their use for functional genes is still debatable as there is no consensus on clustering thresholds. Here, we compare OTU- and ASV-based reconstruction approaches and taxonomy assignment methods, the naive Bayesian classifier (NBC) and Bayesian lowest common ancestor (BLCA) algorithm, using a functional gene data set from the microbial nitrogen-cycling community in the Brouage mudflat (France). A range of OTU similarity thresholds and ASVs were used to compare amoA (ammonia-oxidizing archaea [AOA] and ammonia-oxidizing bacteria [AOB]), nxrB, nirS, nirK, and nrfA communities between differing sedimentary structures. Significant effects of the sedimentary structure on weighted UniFrac (WUniFrac) distances were observed for AOA amoA when using ASVs, an OTU at a threshold of 97% sequence identity (OTU-97%), and OTU-85%; AOB amoA when using OTU-85%; and nirS when using ASV, OTU-90%, and OTU-85%. For AOB amoA, significant effects of the sedimentary structures on UniFrac distances were observed when using OTU-97% but not ASVs, and the inverse was found for nrfA. Interestingly, conclusions drawn for nirK and nxrB were consistent between amplicon reconstruction methods. We also show that when the sequences in the reference database are related to the environment in question, the BLCA algorithm leads to more phylogenetically relevant classifications. However, when the reference database contains sequences more dissimilar to the ones retrieved, the NBC obtains more information. IMPORTANCE Several analysis pipelines are available to microbial ecologists to process amplicon sequencing data, yet to date, there is no consensus as to the most appropriate method, and it becomes more difficult for genes that encode a specific function (functional genes). Standardized approaches need to be adopted to increase the reliability and reproducibility of environmental amplicon-sequencing-based data sets. In this paper, we argue that the recently developed ASV approach offers a better opportunity to achieve such standardization than OTUs for functional genes. We also propose a comprehensive framework for quality filtering of the sequencing reads based on protein sequence verification.

Relationships between Sphaerulina musiva Infection and the Populus Microbiome and Metabolome

Pathogenic fungal infections in plants may, in some cases, lead to downstream systematic impacts on the plant metabolome and microbiome that may either alleviate or exacerbate the effects of the fungal pathogen. While Sphaerulina musiva is a well-characterized fungal pathogen which infects Populus tree species, an important wood fiber and biofuel feedstock, little is known about its systematic effects on the metabolome and microbiome of Populus. Here, we investigated the metabolome of Populus trichocarpa and Populus deltoides leaves and roots and the microbiome of the leaf and root endospheres, phylloplane, and rhizosphere to understand the systematic impacts of S. musiva abundance and infection on Populus species in a common garden field setting. We found that S. musiva is indeed present in both P. deltoides and P. trichocarpa, but S. musiva abundance was not statistically related to stem canker onset. We also found that the leaf and root metabolomes significantly differ between the two Populus species and that certain leaf metabolites, particularly the phenolic glycosides salirepin and salireposide, are diminished in canker-infected P. trichocarpa trees compared to their uninfected counterparts. Furthermore, we found significant associations between the metabolome, S. musiva abundance, and microbiome composition and α-diversity, particularly in P. trichocarpa leaves. Our results show that S. musiva colonizes both resistant and susceptible hosts and that the effects of S. musiva on susceptible trees are not confined to the site of canker infection. IMPORTANCE Poplar (Populus spp.) trees are ecologically and economically important trees throughout North America. However, many western North American poplar plantations are at risk due to the introduction of the nonnative fungal pathogen Sphaerulina musiva, which causes leaf spot and cankers, limiting their production. To better understand the interactions among the pathogen S. musiva, the poplar metabolome, and the poplar microbiome, we collected leaf, root, and rhizosphere samples from poplar trees consisting of 10 genotypes and two species with differential resistance to S. musiva in a common garden experiment. Here, we outline the nuanced relationships between the poplar metabolome, microbiome, and S. musiva, showing that S. musiva may affect poplar trees in tissues distal to the site of infection (i.e., stem). Our research contributes to improving the fundamental understanding of S. musiva and Populus sp. ecology and the utility of a holobiont approach in understanding plant disease.

New-Generation Sequencing Technology in Diagnosis of Fungal Plant Pathogens: A Dream Comes True?

The fast and continued progress of high-throughput sequencing (HTS) and the drastic reduction of its costs have boosted new and unpredictable developments in the field of plant pathology. The cost of whole-genome sequencing, which, until few years ago, was prohibitive for many projects, is now so affordable that a new branch, phylogenomics, is being developed. Fungal taxonomy is being deeply influenced by genome comparison, too. It is now easier to discover new genes as potential targets for an accurate diagnosis of new or emerging pathogens, notably those of quarantine concern. Similarly, with the development of metabarcoding and metagenomics techniques, it is now possible to unravel complex diseases or answer crucial questions, such as "What's in my soil?", to a good approximation, including fungi, bacteria, nematodes, etc. The new technologies allow to redraw the approach for disease control strategies considering the pathogens within their environment and deciphering the complex interactions between microorganisms and the cultivated crops. This kind of analysis usually generates big data that need sophisticated bioinformatic tools (machine learning, artificial intelligence) for their management. Herein, examples of the use of new technologies for research in fungal diversity and diagnosis of some fungal pathogens are reported.

Potential of Meta-Omics to Provide Modern Microbial Indicators for Monitoring Soil Quality and Securing Food Production

Soils are fundamental resources for agricultural production and play an essential role in food security. They represent the keystone of the food value chain because they harbor a large fraction of biodiversity-the backbone of the regulation of ecosystem services and "soil health" maintenance. In the face of the numerous causes of soil degradation such as unsustainable soil management practices, pollution, waste disposal, or the increasing number of extreme weather events, it has become clear that (i) preserving the soil biodiversity is key to food security, and (ii) biodiversity-based solutions for environmental monitoring have to be developed. Within the soil biodiversity reservoir, microbial diversity including Archaea, Bacteria, Fungi and protists is essential for ecosystem functioning and resilience. Microbial communities are also sensitive to various environmental drivers and to management practices; as a result, they are ideal candidates for monitoring soil quality assessment. The emergence of meta-omics approaches based on recent advances in high-throughput sequencing and bioinformatics has remarkably improved our ability to characterize microbial diversity and its potential functions. This revolution has substantially filled the knowledge gap about soil microbial diversity regulation and ecology, but also provided new and robust indicators of agricultural soil quality. We reviewed how meta-omics approaches replaced traditional methods and allowed developing modern microbial indicators of the soil biological quality. Each meta-omics approach is described in its general principles, methodologies, specificities, strengths and drawbacks, and illustrated with concrete applications for soil monitoring. The development of metabarcoding approaches in the last 20 years has led to a collection of microbial indicators that are now operational and available for the farming sector. Our review shows that despite the recent huge advances, some meta-omics approaches (e.g., metatranscriptomics or meta-proteomics) still need developments to be operational for environmental bio-monitoring. As regards prospects, we outline the importance of building up repositories of soil quality indicators. These are essential for objective and robust diagnosis, to help actors and stakeholders improve soil management, with a view to or to contribute to combining the food and environmental quality of next-generation farming systems in the context of the agroecological transition.

Changes in precipitation patterns can destabilize plant species coexistence via changes in plant-soil feedback

Climate change can alter species coexistence through changes in biotic interactions. By describing reciprocal interactions between plants and soil microbes, plant-soil feedback (PSF) has emerged as a powerful framework for predicting plant species coexistence and community dynamics, but little is known about how PSF will respond to changing climate conditions. Hence, the context dependency of PSF has recently gained attention. Water availability is a major driver of all biotic interactions, and it is expected that precipitation patterns will change with ongoing climate change. We tested how soil water content affects PSF by conducting a full factorial pairwise PSF experiment using eight plant species common to southeastern United States coastal prairies under three watering treatments. We found coexistence-stabilizing negative PSF at drier-than-average conditions shifted to coexistence-destabilizing positive PSF under wetter-than-average conditions. A simulation model parameterized with the experimental results supports the prediction that more positive PSF accelerates the erosion of diversity within communities while decreasing the predictability in plant community composition. Our results underline the importance of considering environmental context dependency of PSF in light of a rapidly changing climate.

Seed microbiota revealed by a large-scale meta-analysis including 50 plant species

Seed microbiota constitutes a primary inoculum for plants that is gaining attention owing to its role for plant health and productivity. Here, we performed a meta-analysis on 63 seed microbiota studies covering 50 plant species to synthesize knowledge on the diversity of this habitat. Seed microbiota are diverse and extremely variable, with taxa richness varying from one to thousands of taxa. Hence, seed microbiota presents a variable (i.e. flexible) microbial fraction but we also identified a stable (i.e. core) fraction across samples. Around 30 bacterial and fungal taxa are present in most plant species and in samples from all over the world. Core taxa, such as Pantoea agglomerans, Pseudomonas viridiflava, P. fluorescens, Cladosporium perangustum and Alternaria sp., are dominant seed taxa. The characterization of the core and flexible seed microbiota provided here will help uncover seed microbiota roles for plant health and design effective microbiome engineering.