Elevational dynamics in soil microbial co-occurrence: Disentangling biotic and abiotic influences on bacterial and fungal networks on Mt. Seorak

Abstract

● Fungi outperformed bacterial in maintaining the microbial co-occurrence networks.

● Fungi showed different elevational network co-occurrence pattern from bacteria.

● Distinct biotic/abiotic factors influenced bacterial and fungal network dynamics.

The interplay between soil micro-organisms in mountain ecosystems critically influences soil biogeochemical cycles and ecosystem processes. However, factors affecting the co-occurrence patterns of soil microbial communities remain unclear. In an attempt to understand how these patterns shift with elevation and identify the key explanatory factors underpinning these changes, we studied soil bacterial and fungal co-occurrence networks on Mt. Seorak, Republic of Korea. Amplicon sequencing was used to target the 16S rRNA gene and ITS2 region for bacteria and fungi, respectively. In contrast to bacteria, we found that fungi were predominantly situated in the core positions of the network, with significantly weakened co-occurrence with increasing elevation. The different co-occurrence patterns of fungal and bacterial communities could be resulted from their distinct responses to various environments. Both abiotic and biotic factors contributed significantly to shaping co-occurrence networks of bacterial and fungal communities. Fungal richness, bacterial community composition (indicated by PCoA1), and soil pH were the predominant factors influencing the variation in the entire microbial co-occurrence network. Biotic factors, such as the composition and diversity of bacterial communities, significantly influenced bacterial co-occurrence networks. External biotic and abiotic factors, including climatic and vegetative conditions, had a significant influence on fungal co-occurrence networks. These findings enhance our understanding of soil microbiota co-occurrences and deepen our knowledge of soil microbiota responses to climatic changes across elevational gradients in mountain ecosystems.

Graphical abstract

Keywords

co-occurrence network / soil microbial community / elevational gradient / soil pH / plant vegetation / climate

References

[1] Adhikari, D., Barik, S.K., Upadhaya, K., 2012. Habitat distribution modelling for reintroduction of Ilex khasiana Purk., a critically endangered tree species of northeastern India. Ecological Engineering40, 37–43.

[2] Ashe, E.C., Comeau, A.M., Zejdlik, K., O’Connell, S.P., 2021. Characterization of bacterial community dynamics of the human mouth throughout decomposition via metagenomic, metatranscriptomic, and culturing techniques. Frontiers in Microbiology12, 689493.

[3] Banerjee, S., Kirkby, C.A., Schmutter, D., Bissett, A., Kirkegaard, J.A., Richardson, A.E., 2016. Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biology and Biochemistry97, 188–198.

[4] Banerjee, S., Walder, F., Büchi, L., Meyer, M., Held, A.Y., Gattinger, A., Keller, T., Charles, R., van der Heijden, M.G.A., 2019. Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. The ISME Journal13, 1722–1736.

[5] Barberán, A., Bates, S.T., Casamayor, E.O., Fierer, N., 2012. Using network analysis to explore co-occurrence patterns in soil microbial communities. The ISME Journal6, 343–351.

[6] Berry, D., Widder, S., 2014. Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Frontiers in Microbiology5, 219.

[7] Bertness, M.D., Callaway, R., 1994. Positive interactions in communities. Trends in Ecology & Evolution9, 191–193.

[8] Cai, Z.Q., Zhang, Y.H., Yang, C., Wang, S., 2018. Land-use type strongly shapes community composition, but not always diversity of soil microbes in tropical China. CATENA165, 369–380.

[9] Chen, B.B., Jiao, S., Luo, S.W., Ma, B.B., Qi, W., Cao, C.D., Zhao, Z.G., Du, G.Z., Ma, X.J., 2021. High soil pH enhances the network interactions among bacterial and archaeal microbiota in alpine grasslands of the Tibetan Plateau. Environmental Microbiology23, 464–477.

[10] Chun, J.H., Lee, C.B., 2013. Assessing the effects of climate change on the geographic distribution of Pinus densiflora in Korea using ecological niche model. Korean Journal of Agricultural and Forest Meteorology15, 219–233.

[11] Chun, J.H., Lee, C.B., 2018. Partitioning the regional and local drivers of phylogenetic and functional diversity along temperate elevational gradients on an East Asian peninsula. Scientific Reports8, 2853.

[12] Correa, H., Haltli, B., Duque, C., Kerr, R., 2013. Bacterial communities of the gorgonian octocoral Pseudopterogorgia elisabethae. Microbial Ecology66, 972–985.

[13] Coyte, K.Z., Schluter, J., Foster, K.R., 2015. The ecology of the microbiome: networks, competition, and stability. Science350, 663–666.

[14] Csardi, G., Nepusz, T., 2006. The igraph software package for complex network research. InterJournal, Complex Systems1695, 1–9.

[15] de Vries, F.T., Griffiths, R.I., Bailey, M., Craig, H., Girlanda, M., Gweon, H.S., Hallin, S., Kaisermann, A., Keith, A.M., Kretzschmar, M., Lemanceau, P., Lumini, E., Mason, K.E., Oliver, A., Ostle, N., Prosser, J.I., Thion, C., Thomson, B., Bardgett, R.D., 2018. Soil bacterial networks are less stable under drought than fungal networks. Nature Communications9, 3033.

[16] Dong, K., Tripathi, B., Moroenyane, I., Kim, W., Li, N., Chu, H.Y., Adams, J., 2016. Soil fungal community development in a high Arctic glacier foreland follows a directional replacement model, with a mid-successional diversity maximum. Scientific Reports6, 26360.

[17] Dong, K., Yu, Z., Kerfahi, D., Lee, S.S., Li, N., Yang, T., Adams, J.M., 2022. Soil microbial co-occurrence networks become less connected with soil development in a high Arctic glacier foreland succession. Science of the Total Environment813, 152565.

[18] Erdös, P., 1960. On sets of distances of n points in Euclidean space. Publications of the Mathematical Institute of the Hungarian Academy of Sciences5, 165–169.

[19] Fan, K.K., Cardona, C., Li, Y.T., Shi, Y., Xiang, X.J., Shen, C.C., Wang, H.F., Gilbert, J.A., Chu, H.Y., 2017. Rhizosphere-associated bacterial network structure and spatial distribution differ significantly from bulk soil in wheat crop fields. Soil Biology and Biochemistry113, 275–284.

[20] Fan, K.K., Weisenhorn, P., Gilbert, J.A., Shi, Y., Bai, Y., Chu, H.Y., 2018. Soil pH correlates with the co-occurrence and assemblage process of diazotrophic communities in rhizosphere and bulk soils of wheat fields. Soil Biology and Biochemistry121, 185–192.

[21] Fan, X.X., Pan, H.Y., Ping, Y., Jin, G.Z., Song, F.Q., 2022. The underlying mechanism of soil aggregate stability by fungi and related multiple factor: a review. Eurasian Soil Science55, 242–250.

[22] Faust, K., Raes, J., 2012. Microbial interactions: from networks to models. Nature Reviews Microbiology10, 538–550.

[23] Faust, K., Sathirapongsasuti, J.F., Izard, J., Segata, N., Gevers, D., Raes, J., Huttenhower, C., 2012. Microbial co-occurrence relationships in the human microbiome. PLoS Computational Biology8, e1002606.

[24] Fierer, N., Jackson, R.B., 2006. The diversity and biogeography of soil bacterial communities. Proceedings of the National Academy of Sciences of the United States of America103, 626–631.

[25] Genre, A., Lanfranco, L., Perotto, S., Bonfante, P., 2020. Unique and common traits in mycorrhizal symbioses. Nature Reviews Microbiology18, 649–660.

[26] Goberna, M., Verdú, M., 2022. Cautionary notes on the use of co-occurrence networks in soil ecology. Soil Biology and Biochemistry166, 108534.

[27] Goslee, S.C., Urban, D.L., 2007. The ecodist package for dissimilarity-based analysis of ecological data. Journal of Statistical Software22, 1–19.

[28] Greenblum, S., Turnbaugh, P.J., Borenstein, E., 2012. Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proceedings of the National Academy of Sciences of the United States of America109, 594–599.

[29] Harrell, F.E.Jr., Dupont, C., 2012. Hmisc: Harrell Miscellaneous [Online]. R Package Version 3.9-3.

[30] He, D., Shen, W.J., Eberwein, J., Zhao, Q., Ren, L.J., Wu, Q.L., 2017. Diversity and co-occurrence network of soil fungi are more responsive than those of bacteria to shifts in precipitation seasonality in a subtropical forest. Soil Biology and Biochemistry115, 499–510.

[31] Hernandez, D.J., David, A.S., Menges, E.S., Searcy, C.A., Afkhami, M.E., 2021. Environmental stress destabilizes microbial networks. The ISME Journal15, 1722–1734.

[32] Herren, C.M., McMahon, K.D., 2018. Keystone taxa predict compositional change in microbial communities. Environmental Microbiology20, 2207–2217.

[33] Hirano, H., Takemoto, K., 2019. Difficulty in inferring microbial community structure based on co-occurrence network approaches. BMC Bioinformatics20, 329.

[34] Hooper, D.U., Bignell, D.E., Brown, V.K., Brussard, L., Dangerfield, J.M., Wall, D.H., Wardle, D.A., Coleman, D.C., Giller, K.E., Lavelle, P., Van Der Putten, W.H., De Ruiter, P.C., Rusek, J., Silver, W.L., Tiedje, J.M., Wolters, V., 2000. Interactions between Aboveground and Belowground Biodiversity in Terrestrial Ecosystems: patterns, Mechanisms, and Feedbacks: we assess the evidence for correlation between aboveground and belowground diversity and conclude that a variety of mechanisms could lead to positive, negative, or no relationship—depending on the strength and type of interactions among species. BioScience50, 1049–1061.

[35] Jarvis, S.G., Woodward, S., Taylor, A.F.S., 2015. Strong altitudinal partitioning in the distributions of ectomycorrhizal fungi along a short (300 m) elevation gradient. New Phytologist206, 1145–1155.

[36] Ji, L., Shen, F.Y., Liu, Y., Yang, Y.C., Wang, J., Purahong, W., Yang, L.X., 2022. Contrasting altitudinal patterns and co-occurrence networks of soil bacterial and fungal communities along soil depths in the cold-temperate montane forests of China. CATENA209, 105844.

[37] Jiang, Y.L., Lei, Y.B., Yang, Y., Korpelainen, H., Niinemets, Ü., Li, C.Y., 2018. Divergent assemblage patterns and driving forces for bacterial and fungal communities along a glacier forefield chronosequence. Soil Biology and Biochemistry118, 207–216.

[38] Jiao, S., Xu, Y.Q., Zhang, J., Hao, X., Lu, Y.H., 2019. Core microbiota in agricultural soils and their potential associations with nutrient cycling. mSystems4, e00313–18.

[39] Kim, H., Chun, J.H., Lee, C.B., 2020. Plant diversity and phylogenetic community structure along environmental gradients in a temperate forest, South Korea. The Journal of Animal & Plant Sciences,30( 4), 958–969.

[40] Langfelder, P., Horvath, S., 2008. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics9, 559.

[41] Lee, C.B., Chun, J.H., Ahn, H.H., 2014. Elevational patterns of plant richness and their drivers on an Asian mountain. Nordic Journal of Botany32, 347–357.

[42] Li, J., Wang, X., Wu, J.H., Sun, Y.X., Zhang, Y.Y., Zhao, Y.F., Huang, Z., Duan, W.H., 2023. Climate and geochemistry at different altitudes influence soil fungal community aggregation patterns in alpine grasslands. Science of the Total Environment881, 163375.

[43] Li, J.B., Li, C.N., Kou, Y.P., Yao, M.J., He, Z.L., Li, X.Z., 2020. Distinct mechanisms shape soil bacterial and fungal co-occurrence networks in a mountain ecosystem. FEMS Microbiology Ecology96, fiaa030.

[44] Li, J.R., Chen, L., Wang, H., Ouyang, S., Liu, X.H., Lu, J., 2022. Pattern and drivers of soil fungal community along elevation gradient in the Abies georgei forests of Segila mountains, Southeast Tibet. Global Ecology and Conservation39, e02291.

[45] Lodge, D.J., 1997. Factors related to diversity of decomposer fungi in tropical forests. Biodiversity & Conservation6, 681–688.

[46] Lupatini, M., Suleiman, A.K.A., Jacques, R.J.S., Antoniolli, Z.I., de Siqueira Ferreira, A., Kuramae, E.E., Roesch, L.F.W., 2014. Network topology reveals high connectance levels and few key microbial genera within soils. Frontiers in Environmental Science2, 10.

[47] Ma, B., Wang, H.Z., Dsouza, M., Lou, J., He, Y., Dai, Z.M., Brookes, P.C., Xu, J.M., Gilbert, J.A., 2016. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. The ISME Journal10, 1891–1901.

[48] Marizzoni, M., Gurry, T., Provasi, S., Greub, G., Lopizzo, N., Ribaldi, F., Festari, C., Mazzelli, M., Mombelli, E., Salvatore, M., Mirabelli, P., Franzese, M., Soricelli, A., Frisoni, G.B., Cattaneo, A., 2020. Comparison of bioinformatics pipelines and operating systems for the analyses of 16S rRNA gene amplicon sequences in human fecal samples. Frontiers in Microbiology11, 1262.

[49] Naranjo-Ortiz, M.A., Gabaldón, T., 2019. Fungal evolution: major ecological adaptations and evolutionary transitions. Biological Reviews94, 1443–1476.

[50] Neutel, A.M., Heesterbeek, J.A.P., van de Koppel, J., Hoenderboom, G., Vos, A., Kaldeway, C., Berendse, F., de Ruiter, P.C., 2007. Reconciling complexity with stability in naturally assembling food webs. Nature449, 599–602.

[51] Ni, Y.Y., Yang, T., Ma, Y.Y., Zhang, K.P., Soltis, P.S., Soltis, D.E., Gilbert, J.A., Zhao, Y.P., Fu, C.X., Chu, H.Y., 2021. Soil pH determines bacterial distribution and assembly processes in natural mountain forests of eastern China. Global Ecology and Biogeography30, 2164–2177.

[52] Nilsson, R.H., Larsson, K.H., Taylor, A.F.S., Bengtsson-Palme, J., Jeppesen, T.S., Schigel, D., Kennedy, P., Picard, K., Glöckner, F.O., Tedersoo, L., Saar, I., Kõljalg, U., Abarenkov, K., 2019. The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Research,47, D259–D264.

[53] Núñez, A., Moreno, D.A., 2020. The differential vertical distribution of the airborne biological particles reveals an atmospheric reservoir of microbial pathogens and aeroallergens. Microbial Ecology80, 322–333.

[54] Oksanen, J., 2015. Vegan: An Introduction to Ordination [Online]. available at the website of The Comprehensive R Archive Network

[55] Pollard, K.S., Dudoit, S., van der Laan, M.J., 2005. Multiple testing procedures: the multtest package and applications to genomics. In: Gentleman, R., Carey, V.J., Huber, W., Irizarry, R.A., Dudoit, S., eds. Bioinformatics and Computational Biology Solutions Using R and Bioconductor. New York: Springer, 249–271

[56] Romaní, A.M., Fischer, H., Mille-Lindblom, C., Tranvik, L.J., 2006. Interactions of bacteria and fungi on decomposing litter: differential extracellular enzyme activities. Ecology87, 2559–2569.

[57] Rooney, N., McCann, K., Gellner, G., Moore, J.C., 2006. Structural asymmetry and the stability of diverse food webs. Nature442, 265–269.

[58] Schloss, P.D., Westcott, S.L., Ryabin, T., Hall, J.R., Hartmann, M., Hollister, E.B., Lesniewski, R.A., Oakley, B.B., Parks, D.H., Robinson, C.J., Sahl, J.W., Stres, B., Thallinger, G.G., Van Horn, D.J., Weber, C.F., 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology75, 7537–7541.

[59] Shen, C.C., Gunina, A., Luo, Y., Wang, J.J., He, J.Z., Kuzyakov, Y., Hemp, A., Classen, A.T., Ge, Y., 2020. Contrasting patterns and drivers of soil bacterial and fungal diversity across a mountain gradient. Environmental Microbiology22, 3287–3301.

[60] Shen, C.C., Wang, J., Jing, Z.W., Qiao, N.H., Xiong, C., Ge, Y., 2022. Plant diversity enhances soil fungal network stability indirectly through the increase of soil carbon and fungal keystone taxa richness. Science of the Total Environment818, 151737.

[61] Shi, S.J., Nuccio, E.E., Shi, Z.J., He, Z.L., Zhou, J.Z., Firestone, M.K., 2016. The interconnected rhizosphere: high network complexity dominates rhizosphere assemblages. Ecology Letters19, 926–936.

[62] Song, L.Y., Yang, T., Xia, S.G., Yin, Z., Liu, X., Li, S.P., Sun, R.B., Gao, H.J., Chu, H.Y., Ma, C., 2022. Soil depth exerts stronger impact on bacterial community than elevation in subtropical forests of Huangshan Mountain. Science of the Total Environment852, 158438.

[63] Stark, S., Männistö, M.K., Eskelinen, A., 2014. Nutrient availability and pH jointly constrain microbial extracellular enzyme activities in nutrient-poor tundra soils. Plant and Soil,383, 373–385.

[64] Tripathi, B.M., Stegen, J.C., Kim, M., Dong, K., Adams, J.M., Lee, Y.K., 2018. Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. The ISME Journal12, 1072–1083.

[65] Tu, Q.C., Yan, Q.Y., Deng, Y., Michaletz, S.T., Buzzard, V., Weiser, M.D., Waide, R., Ning, D.L., Wu, L.Y., He, Z.L., Zhou, J.Z., 2020. Biogeographic patterns of microbial co-occurrence ecological networks in six American forests. Soil Biology and Biochemistry148, 107897.

[66] Vieira, L.C., da Silva, D.K.A., de Melo, M.A.C., Escobar, I.E.C., Oehl, F., da Silva, G.A., 2019. Edaphic factors influence the distribution of arbuscular mycorrhizal fungi along an altitudinal gradient of a tropical mountain. Microbial Ecology78, 904–913.

[67] Wagg, C., Bender, S.F., Widmer, F., van der Heijden, M.G.A., 2014. Soil biodiversity and soil community composition determine ecosystem multifunctionality. Proceedings of the National Academy of Sciences of the United States of America111, 5266–5270.

[68] Wang, S., Wang, X.B., Han, X.G., Deng, Y., 2018. Higher precipitation strengthens the microbial interactions in semi-arid grassland soils. Global Ecology and Biogeography27, 570–580.

[69] Wu, L.W., Yang, Y.F., Chen, S., Zhao, M.X., Zhu, Z.W., Yang, S.H., Qu, Y.Y., Ma, Q., He, Z.L., Zhou, J.Z., He, Q., 2016. Long-term successional dynamics of microbial association networks in anaerobic digestion processes. Water Research104, 1–10.

[70] Wu, X.F., Yang, J.J., Ruan, H., Wang, S.N., Yang, Y.R., Naeem, I., Wang, L., Liu, L.E., Wang, D.L., 2021. The diversity and co-occurrence network of soil bacterial and fungal communities and their implications for a new indicator of grassland degradation. Ecological Indicators129, 107989.

[71] Xiao, X., Liang, Y.T., Zhou, S., Zhuang, S.Y., Sun, B., 2018. Fungal community reveals less dispersal limitation and potentially more connected network than that of bacteria in bamboo forest soils. Molecular Ecology27, 550–563.

[72] Yang, T., Adams, J.M., Shi, Y., He, J.S., Jing, X., Chen, L.T., Tedersoo, L., Chu, H.Y., 2017. Soil fungal diversity in natural grasslands of the Tibetan Plateau: associations with plant diversity and productivity. New Phytologist215, 756–765.

[73] Yang, T., Tedersoo, L., Lin, X.W., Fitzpatrick, M.C., Jia, Y.S., Liu, X., Ni, Y.Y., Shi, Y., Lu, P.P., Zhu, J.G., Chu, H.Y., 2020. Distinct fungal successional trajectories following wildfire between soil horizons in a cold-temperate forest. New Phytologist227, 572–587.

[74] Yang, T., Tedersoo, L., Liu, X., Gao, G.F., Dong, K., Adams, J.M., Chu, H.Y., 2022. Fungi stabilize multi-kingdom community in a high elevation timberline ecosystem. iMeta1, e49.

[75] Yang, T., Tedersoo, L., Soltis, P.S., Soltis, D.E., Gilbert, J.A., Sun, M., Shi, Y., Wang, H.F., Li, Y.T., Zhang, J., Chen, Z.D., Lin, H.Y., Zhao, Y.P., Fu, C.X., Chu, H.Y., 2019. Phylogenetic imprint of woody plants on the soil mycobiome in natural mountain forests of eastern China. The ISME Journal13, 686–697.

[76] Yang, Y., Shi, Y., Kerfahi, D., Ogwu, M.C., Wang, J.J., Dong, K., Takahashi, K., Moroenyane, I., Adams, J.M., 2021. Elevation-related climate trends dominate fungal co-occurrence network structure and the abundance of keystone taxa on Mt. Norikura, Japan. Science of the Total Environment799, 149368.

[77] Yoon, S.H., Ha, S.M., Kwon, S., Lim, J., Kim, Y., Seo, H., Chun, J., 2017. Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. International Journal of Systematic and Evolutionary Microbiology67, 1613–1617.

[78] Yu, Z., Zou, S.Q., Li, N., Kerfahi, D., Lee, C., Adams, J., Kwak, H.J., Kim, J., Lee, S.S., Dong, K., 2021. Elevation-related climatic factors dominate soil free-living nematode communities and their co-occurrence patterns on Mt. Halla, South Korea. Ecology and Evolution11, 18540–18551.

[79] Yun, J.I., 2010. Agroclimatic maps augmented by a GIS technology. Korean Journal of Agricultural and Forest Meteorology12, 63–73.

[80] Zhang, B.G., Zhang, J., Liu, Y., Shi, P., Wei, G.H., 2018. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale. Soil Biology and Biochemistry118, 178–186.

[81] Zhou, J.Z., Deng, Y., Shen, L.N., Wen, C.Q., Yan, Q.Y., Ning, D.L., Qin, Y.J., Xue, K., Wu, L.Y., He, Z.L., Voordeckers, J.W., Van Nostrand, J.D., Buzzard, V., Michaletz, S.T., Enquist, B.J., Weiser, M.D., Kaspari, M., Waide, R., Yang, Y.F., Brown, J.H., 2016. Temperature mediates continental-scale diversity of microbes in forest soils. Nature Communications7, 12083.

[82] Zhu, B.J., Li, C.N., Wang, J.M., Li, J.B., Li, X.Z., 2020. Elevation rather than season determines the assembly and co-occurrence patterns of soil bacterial communities in forest ecosystems of Mount Gongga. Applied Microbiology and Biotechnology104, 7589–7602.

RIGHTS & PERMISSIONS

Higher Education Press