Deadwood decomposition and other environmental processes mediated by microbial communities are generally studied with composite sampling strategies, where deadwood is collected from multiple locations in a large volume, that produce an average microbial community. In this study, we used amplicon sequencing to compare fungal and bacterial communities sampled with either traditional, composite samples, or small, 1 cm3 cylinders from a discrete location within decomposing European beech (Fagus sylvatica L.) tree trunks. We found that bacterial richness and evenness is lower in small samples when compared to composite samples. There was no significant difference in fungal alpha diversity between different sampling scales, suggesting that visually defined fungal domains are not restricted to a single species. Additionally, we found that composite sampling may obscure variation in community composition and this affects the understanding of microbial associations that are detected. For future experiments in environmental microbiology, we recommend that scale is explicitly considered as a factor and properly selected to correspond with the questions asked. Studies of microbial functions or associations may require samples to be collected at a finer scale than is currently practised.

Laboratory of Environmental Microbiology Institute of Microbiology of the Czech Academy of Sciences Vídeňská 1083 142 20 Prague Czechia

See more in PubMed

Ahlmann-Eltze C. ggsignif: significance brackets for “ggplot2”. 2019. https://CRAN.R-project.org/package=ggsignif. (last accessed on September 01, 2022).

Alteio LV, Séneca J, Canarini Aet al. . A critical perspective on interpreting amplicon sequencing data in soil ecological research. Soil Biol Biochem. 2021;160:108357.

Auguie B. gridExtra: miscellaneous functions for “grid” graphics. 2017. https://CRAN.R-project.org/package=gridExtra. (last accessed on September 01, 2022).

Avis PG, Branco S, Tang Yet al. . Pooled samples bias fungal community descriptions. Mol Ecol Resour. 2010;10:135–41. PubMed

Baldrian P, Zrůstová P, Tláskal Vet al. . Fungi associated with decomposing deadwood in a natural beech-dominated forest. Fung Ecol. 2016;23:109–22.

Baldrian P. Forest microbiome: diversity, complexity and dynamics. Banin E (ed.), FEMS Microbiol Rev. 2016;41:fuw040. doi: 10.1093/femsre/fuw040. PubMed

Baldrian P. Microbial activity and the dynamics of ecosystem processes in forest soils. Curr Opin Microbiol. 2017;37:128–34. PubMed

Boddy L. Interspecific combative interactions between wood-decaying basidiomycetes. FEMS Microbiol Ecol. 2000;31:185–94. PubMed

de Boer W, Folman LB, Summerbell RCet al. . Living in a fungal world: impact of fungi on soil bacterial niche development. FEMS Microbiol Rev. 2005;29:795–811. PubMed

Bolyen E, Rideout JR, Dillon MRet al. . Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7. PubMed PMC

Brabcová V, Nováková M, Davidová Aet al. . Dead fungal mycelium in forest soil represents a decomposition hotspot and a habitat for a specific microbial community. New Phytol. 2016;210:1369–81. PubMed

Cai Y-J, Liu Z-A, Zhang Set al. . Microbial community structure is stratified at the millimeter-scale across the soil–water interface. ISME Commun. 2022;2:53. PubMed PMC

Callahan BJ, McMurdie PJ, Rosen MJet al. . DADA2: high-resolution sample inference from Illumina amplicon data. Olsson IAS (ed.),Nat Methods. 2016;13:581–3. PubMed PMC

Caporaso JG, Lauber CL, Walters WAet al. . Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4. PubMed PMC

Carmona-Hernandez S, Reyes-Pérez J, Chiquito-Contreras Ret al. . Biocontrol of postharvest fruit fungal diseases by bacterial antagonists: a review. Agronomy. 2019;9:121.

Charrad M, Ghazzali N, Boiteau Vet al. . NbClust: an R package for determining the relevant number of clusters in a data set. J Stat Soft. 2014;61:1–36.

Christofides SR, Bettridge A, Farewell Det al. . The influence of migratory Paraburkholderia on growth and competition of wood-decay fungi. Fung Ecol. 2020;45:100937.

Cordero OX, Datta MS. Microbial interactions and community assembly at microscales. Curr Opin Microbiol. 2016;31:227–34. PubMed PMC

Crowther TW, Boddy L, Jones TH. Outcomes of fungal interactions are determined by soil invertebrate grazers: grazers alter fungal community. Ecol Lett. 2011;14:1134–42. PubMed

Deveau A, Bonito G, Uehling Jet al. . Bacterial–fungal interactions: ecology, mechanisms and challenges. FEMS Microbiol Rev. 2018;42:335–52. PubMed

Dunnington D. ggspatial: spatial data framework for ggplot2. 2021. https://CRAN.R-project.org/package=ggspatial. (last accessed on September 01, 2022).

Engel M, Behnke A, Bauerfeld Set al. . Sample pooling obscures diversity patterns in intertidal ciliate community composition and structure. FEMS Microbiol Ecol. 2012;79:741–50. PubMed

Eskalen A, Stouthamer R, Lynch SCet al. . Host range of Fusarium dieback and its ambrosia beetle (Coleoptera: scolytinae) vector in Southern California. Plant Dis. 2013;97:938–51. PubMed

Faust K, Sathirapongsasuti JF, Izard Jet al. . Microbial co-occurrence relationships in the human microbiome. Ouzounis CA (ed.), PLoS Comput Biol. 2012;8:e1002606. PubMed PMC

Fleishman SM, Eissenstat DM, Bell THet al. . Functionally-explicit sampling can answer key questions about the specificity of plant–microbe interactions. Environ Microbiome. 2022;17:51. PubMed PMC

FOREST EUROPE . State of Europe's forests 2020. 2020. https://foresteurope.org/state-of-europes-forests/. (last accessed August 10, 2022).

Hagge J, Bässler C, Gruppe Aet al. . Bark coverage shifts assembly processes of microbial decomposer communities in dead wood. Proc R Soc B. 2019;286:20191744. PubMed PMC

Heilmann-Clausen J, Boddy L. Inhibition and stimulation effects in communities of wood decay fungi: exudates from colonized wood influence growth by other species. Microb Ecol. 2005;49:399–406. PubMed

Hiscox J, Boddy L. Armed and dangerous – chemical warfare in wood decay communities. Fung Biol Rev. 2017;31:169–84.

Hoppe B, Kahl T, Karasch Pet al. . Network analysis reveals ecological links between N-fixing bacteria and wood-decaying fungi. Desvaux M (ed.), PLoS ONE. 2014;9:e88141. PubMed PMC

Huber N. ggdendroplot: create dendrograms for ggplot2. 2021. https://github.com/NicolasH2/ggdendroplot. (last accessed September 1, 2022).

Ihrmark K, Bödeker ITM, Cruz-Martinez Ket al. . New primers to amplify the fungal ITS2 region - evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol Ecol. 2012;82:666–77. PubMed

Kleine Bardenhorst S, Vital M, Karch Aet al. . Richness estimation in microbiome data obtained from denoising pipelines. Comput Struct Biotechnol J. 2022;20:508–20. PubMed PMC

Kõljalg U, Nilsson HR, Schigel Det al. . The taxon hypothesis paradigm—on the unambiguous detection and communication of taxa. Microorganisms. 2020;8:1910. PubMed PMC

Ladau J, Eloe-Fadrosh EA. Spatial, temporal, and phylogenetic scales of microbial ecology. Trends Microbiol. 2019;27:662–9. PubMed

McKnight DT, Huerlimann R, Bower DSet al. . Methods for normalizing microbiome data: an ecological perspective. Jarman S (ed.), Methods Ecol Evol. 2019;10:389–400.

McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. Watson M (ed.), PLoS ONE. 2013;8:e61217–. PubMed PMC

McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. McHardy AC (ed.), PLoS Comput Biol. 2014;10:e1003531. PubMed PMC

Maechler M, Rousseeuw P, Struyf Aet al. . Cluster: cluster analysis basics and extensions. 2021. https://CRAN.R-project.org/package=cluster. ( last accessed September 1, 2022).

Manter DK, Weir TL, Vivanco JM. Negative effects of sample pooling on PCR-based estimates of soil microbial richness and community structure. Appl Environ Microbiol. 2010;76:2086–90. PubMed PMC

Moll J, Kellner H, Leonhardt Set al. . Bacteria inhabiting deadwood of 13 tree species are heterogeneously distributed between sapwood and heartwood: bacteria in deadwood of 13 different tree species. Environ Microbiol. 2018;20:3744–56. PubMed

Neuwirth E. RColorBrewer: colorBrewer palettes. 2014. https://CRAN.R-project.org/package=RColorBrewer. (last accessed September 1, 2022).

Nilsson RH, Larsson K-H, Taylor AFSet al. . The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2019;47:D259–64. PubMed PMC

O'Brien SL, Gibbons SM, Owens SMet al. . Spatial scale drives patterns in soil bacterial diversity: spatial scale drives soil diversity. Environ Microbiol. 2016;18:2039–51. PubMed PMC

O'Leary J, Eastwood D, Müller Cet al. . Emergent properties arising from spatial heterogeneity influence fungal community dynamics. Fung Ecol. 2018;33:32–9.

Oksanen J, Blanchet FG, Friendly Met al. . vegan: community ecology package. 2019. https://CRAN.R-project.org/package=vegan. (last accessed September 1, 2022).

Paap T, de Beer ZW, Migliorini Det al. . The polyphagous shot hole borer (PSHB) and its fungal symbiont fusarium euwallaceae: a new invasion in South Africa. Australasian Plant Pathol. 2018;47:231–7.

Pebesma E. Simple features for R: standardized support for spatial vector data. R J. 2018;10:439.

Peschel S, Müller CL, von Mutius Eet al. . NetCoMi: network construction and comparison for microbiome data in R. Briefings Bioinf. 2020;22:bbaa290. PubMed PMC

Quast C, Pruesse E, Yilmaz Pet al. . The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6. PubMed PMC

R Core Team . R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing, 2022. https://www.R-project.org/

Rinta-Kanto JM, Sinkko H, Rajala Tet al. . Natural decay process affects the abundance and community structure of bacteria and archaea in Picea abies logs. de Boer W (ed.), FEMS Microbiol Ecol. 2016;92:fiw087. PubMed

Rolling T, Zhai B, Frame Jet al. . Customization of a DADA2-based pipeline for fungal internal transcribed spacer 1 (ITS1) amplicon data sets. JCI Insight. 2022;7. doi: 10.1172/jci.insight.151663. PubMed PMC

Seibold S, Müller J, Baldrian Pet al. . Fungi associated with beetles dispersing from dead wood – let’s take the beetle bus!. Fung Ecol. 2019;39:100–8.

Shipunov A, Murrell P, D'Orazio Met al. . Shipunov: miscellaneous functions from alexey Shipunov. 2022. https://CRAN.R-project.org/package=shipunov

Simonsohn U, Gruson H. Groundhog: the simplest solution to version-control for CRAN packages. 2021. https://CRAN.R-project.org/package=groundhog. ( last accessed September 1, 2022).

Slowikowski K. ggrepel: automatically position non-overlapping text labels with “ggplot2”. 2021. https://CRAN.R-project.org/package=ggrepel. (last accessed September 1, 2022).

South A. Rnaturalearth: world map data from natural Earth. 2017a. https://CRAN.R-project.org/package=rnaturalearth. (last accessed September 1, 2022).

South A. Rnaturalearthdata: world vector map data from natural Earth used in “Rnaturalearth”. 2017b. https://CRAN.R-project.org/package=rnaturalearthdata. ( last accessed September 1, 2022).

Steenwyk JL, Rokas A. ggpubfigs: colorblind-friendly color palettes and ggplot2 graphic system extensions for publication-quality scientific figures. Newton ILG (ed.), Microbiol Resour Announc. 2021;10:e00871–21. PubMed PMC

Štursová M, Bárta J, Šantrůčková Het al. . Small-scale spatial heterogeneity of ecosystem properties, microbial community composition and microbial activities in a temperate mountain forest soil. de Boer W (ed.), FEMS Microbiol Ecol. 2016;92:fiw185. PubMed

Tláskal V, Brabcová V, Větrovský Tet al. . Complementary roles of wood-inhabiting fungi and bacteria facilitate deadwood decomposition. Faust K (ed.),Msystems. 2021a;6:e01078–20. PubMed PMC

Tláskal V, Brabcová V, Větrovský Tet al. . Metagenomes, metatranscriptomes and microbiomes of naturally decomposing deadwood. Sci Data. 2021b;8:198. PubMed PMC

de Vries A, Ripley BD. Ggdendro: create dendrograms and tree diagrams using “Ggplot2”. 2020. https://CRAN.R-project.org/package=ggdendro. (last accessed September 1, 2022).

West JR, Whitman T. Disturbance by soil mixing decreases microbial richness and supports homogenizing community assembly processes. FEMS Microbiol Ecol. 2022;98:fiac089. PubMed PMC

Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer, 2016. https://ggplot2.tidyverse.org. ( last accessed September 1, 2022).

Wickham H. stringr: simple, consistent wrappers for common string operations. 2019. https://CRAN.R-project.org/package=stringr. (last accessed September 1, 2022).

Willis AD. Rarefaction, alpha diversity, and statistics. Front Microbiol. 2019;10. doi: 10.3389/fmicb.2019.02407. PubMed PMC

Xiong B-J, Kleinsteuber S, Sträuber Het al. . Impact of fungal hyphae on growth and dispersal of obligate anaerobic bacteria in aerated habitats. Bailey MJ (ed.), Mbio. 2022;13:e00769–22. PubMed PMC

Yan L. Ggvenn: draw Venn diagram by “Ggplot2”. 2021. https://CRAN.R-project.org/package=ggvenn. (last accessed September 1, 2022).

Yilmaz P, Parfrey LW, Yarza Pet al. . The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucl Acids Res. 2014;42:D643–8. PubMed PMC

Zhang L, Xu M, Liu Yet al. . Carbon and phosphorus exchange may enable cooperation between an arbuscular mycorrhizal fungus and a phosphate-solubilizing bacterium. New Phytol. 2016;210: 1022–32. PubMed

Editorial: Theme issue on the ecology of soil microorganisms