Deconvoluting the relative contributions made by specific biotic and abiotic drivers to soil fungal community compositions facilitates predictions about the functional responses of ecosystems to environmental changes, such as losses of plant diversity, but it is hindered by the complex interactions involved. Experimental assembly of tree species allows separation of the respective effects of plant community composition (biotic components) and soil properties (abiotic components), enabling much greater statistical power than can be achieved in observational studies. We therefore analyzed these contributions by assessing, via pyrotag sequencing of the internal transcribed spacer (ITS2) rDNA region, fungal communities in young subtropical forest plots included in a large experiment on the effects of tree species richness. Spatial variables and soil properties were the main drivers of soil fungal alpha and beta-diversity, implying strong early-stage environmental filtering and dispersal limitation. Tree related variables, such as tree community composition, significantly affected arbuscular mycorrhizal and pathogen fungal community structure, while differences in tree host species and host abundance affected ectomycorrhizal fungal community composition. At this early stage of the experiment, only a limited amount of carbon inputs (rhizodeposits and leaf litter) was being provided to the ecosystem due to the size of the tree saplings, and persisting legacy effects were observed. We thus expect to find increasing tree related effects on fungal community composition as forest development proceeds.

Background Amplicon sequencing of phylogenetic marker genes, e.g., 16S, 18S, or ITS ribosomal RNA sequences, is still the most commonly used method to determine the composition of microbial communities. Microbial ecologists often have expert knowledge on their biological question and data analysis in general, and most research institutes have computational infrastructures to use the bioinformatics command line tools and workflows for amplicon sequencing analysis, but requirements of bioinformatics skills often limit the efficient and up-to-date use of computational resources. Results We present dadasnake, a user-friendly, 1-command Snakemake pipeline that wraps the preprocessing of sequencing reads and the delineation of exact sequence variants by using the favorably benchmarked and widely used DADA2 algorithm with a taxonomic classification and the post-processing of the resultant tables, including hand-off in standard formats. The suitability of the provided default configurations is demonstrated using mock community data from bacteria and archaea, as well as fungi. Conclusions By use of Snakemake, dadasnake makes efficient use of high-performance computing infrastructures. Easy user configuration guarantees flexibility of all steps, including the processing of data from multiple sequencing platforms. It is easy to install dadasnake via conda environments. dadasnake is available at https://github.com/a-h-b/dadasnake.

No abstract