Inter-phylum negative interactions affect soil bacterial community dynamics and functions during soybean development under long-term nitrogen fertilization

Abstract

Understanding interspecies interactions is essential to predict the response of microbial communities to exogenous perturbation. Herein, rhizospheric and bulk soils were collected from five developmental stages of soybean, which grew in soils receiving 16-year nitrogen inputs. Bacterial communities and functional profiles were examined using high-throughput sequencing and quantitative PCR, respectively. The objective of this study was to identify the key bacterial interactions that influenced community dynamics and functions. We found that the stages of soybean development outcompeted nitrogen fertilization management in shaping bacterial community structure, while fertilization treatments significantly shaped the abundance distribution of nitrogen functional genes. Temporal variations in bacterial abundances increased in bulk soils, especially at the stage of soybean branching, which helps to infer underlying negative interspecies interactions. Members of Cyanobacteria and Actinobacteria actively engaged in inter-phylum negative interactions in bulk soils and soybean rhizosphere, respectively. Furthermore, the negative interactions between nitrogen-fixing functional groups and the reduction of nifH gene abundance were coupled during soybean development, which may help to explain the linkages between population dynamics and functions. Overall, these findings highlight the importance of inter-phylum negative interactions in shaping the correlation patterns of bacterial communities and in determining soil functional potential.

Keywords

Nitrogen management / Soybean development / Community spatiotemporal variations / Bacterial interactions / Functional groups / Nitrogen-cycling genes

References

[1] AbregoN, RoslinT, HuotariT, TackAJM, LindahlBD, TikhonovG, SomervuoP, SchmidtNM, OvaskainenO. Accounting for environmental variation in co-occurrence modelling reveals the importance of positive interactions in root-associated fungal communities. Mol Ecol, 2020, 29(14):2736-2746

[2] AglerMT, RuheJ, KrollS, MorhennC, KimST, WeigelD, KemenEM. Microbial hub taxa link host and abiotic factors to plant microbiome variation. PLoS Biol, 2016, 14(1):e1002352

[3] AsilogluR, SevilirB, SamuelSO, AycanM, AkcaMO, SuzukiK, MuraseJ, TurgayOC, HaradaN. Effect of protists on rhizobacterial community composition and rice plant growth in a biochar amended soil. Biol Fertil Soils, 2021, 57(2):293-304

[4] BanerjeeS, SchlaeppiK, van der HeijdenMGA. Keystone taxa as drivers of microbiome structure and functioning. Nat Rev Microbiol, 2018, 16(9):567-576

[5] BastianM, HeymannS, JacomyM. Gephi: an open source software for exploring and manipulating networks. ICWSM, 2009, 8: 361-362

[6] BensonHP, BoncompagniE, GuerinotML. An iron uptake operon required for proper nodule development in the Bradyrhizobium japonicum-soybean symbiosis. Mol Plant Microbe In, 2005, 18(9):950-959

[7] BokulichNA, SubramanianS, FaithJJ, GeversD, GordonJI, KnightR, MillsDA, CaporasoJG. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods, 2013, 10(1):57-59

[8] BrownSP, GrilloMA, PodowskiJC, HeathKD. Soil origin and plant genotype structure distinct microbiome compartments in the model legume Medicago truncatula. Microbiome, 2020, 8(1):139

[9] ButlerJL, WilliamsMA, BottomleyPJ, MyroldDD. Microbial community dynamics associated with rhizosphere carbon flow. Appl Environ Microbiol, 2003, 69(11):6793-6800

[10] CamenzindT, HattenschwilerS, TresederKK, LehmannA, RilligMC. Nutrient limitation of soil microbial processes in tropical forests. Ecol Monogr, 2018, 88(1):4-21

[11] ChenSM, WaghmodeTR, SunRB, KuramaeEE, HuCS, LiuBB. Root-associated microbiomes of wheat under the combined effect of plant development and nitrogen fertilization. Microbiome, 2019, 7(1):136

[12] CostaR, GotzM, MrotzekN, LottmannJ, BergG, SmallaK. Effects of site and plant species on rhizosphere community structure as revealed by molecular analysis of microbial guilds. FEMS Microbiol Ecol, 2006, 56(2):236-249

[13] Cruz-Morales P, Ramos-Aboites HE, Licona-Cassani C, Selem-Mojica N, Mejia-Ponce PM, Souza-Saldivar V, Barona-Gomez F (2017) Actinobacteria phylogenomics, selective isolation from an iron oligotrophic environment and siderophore functional characterization, unveil new desferrioxamine traits. FEMS Microbiol Ecol 93:fix086

[14] CsardiG, NepuszT. The igraph software package for complex network research. InterJ, 2006, 1695: 1-9

[15] EvansSE, WallensteinMD. Climate change alters ecological strategies of soil bacteria. Ecol Lett, 2014, 17(2):155-164

[16] FanKK, CardonaC, LiYT, ShiY, XiangXJ, ShenCC, WangHF, GilbertJA, ChuHY. Rhizosphere-associated bacterial network structure and spatial distribution differ significantly from bulk soil in wheat crop fields. Soil Biol Biochem, 2017, 113: 275-284

[17] FanKK, Delgado-BaquerizoM, GuoXS, WangDZ, ZhuYG, ChuHY. Biodiversity of key-stone phylotypes determines crop production in a 4-decade fertilization experiment. ISME J, 2021, 15(2):550-561

[18] FanKK, WeisenhornP, GilbertJA, ShiY, BaiY, ChuHY. Soil pH correlates with the co-occurrence and assemblage process of diazotrophic communities in rhizosphere and bulk soils of wheat fields. Soil Biol Biochem, 2018, 121: 185-192

[19] FaustK, RaesJ. Microbial interactions: from networks to models. Nat Rev Microbiol, 2012, 10(8):538-550

[20] FredricksonA, StephanopoulosG. Microbial competition. Science, 1981, 213(4511):972-979

[21] GarlandJL. Patterns of potential C source utilization by rhizosphere communities. Soil Biol Biochem, 1996, 28(2):223-230

[22] GhaiR, MizunoCM, PicazoA, CamachoA, Rodriguez-ValeraF. Metagenomics uncovers a new group of low GC and ultra-small marine Actinobacteria. Sci Rep-Uk, 2013, 3(1):2471

[23] GhaiR, MizunoCM, PicazoA, CamachoA, Rodriguez-ValeraF. Key roles for freshwater Actinobacteria revealed by deep metagenomic sequencing. Mol Ecol, 2014, 23(24):6073-6090

[24] GriffithsBS, RitzK, WheatleyR, KuanHL, BoagB, ChristensenS, EkelundF, SorensenSJ, MullerS, et al.. An examination of the biodiversity-ecosystem function relationship in arable soil microbial communities. Soil Biol Biochem, 2001, 33(12-13):1713-1722

[25] GrosskopfT, SoyerOS. Synthetic microbial communities. Curr Opin Microbiol, 2014, 18: 72-77

[26] GuimeraR, AmaralLAN. Functional cartography of complex metabolic networks. Nature, 2005, 433(7028):895-900

[27] HeY, SenB, ZhouS, XieN, ZhangY, ZhangJ, WangG. Distinct seasonal patterns of bacterioplankton abundance and dominance of phyla α-Proteobacteria and Cyanobacteria in Qinhuangdao coastal waters off the Bohai Sea. Front Microbiol, 2017, 8: 1579

[28] HerrenCM, McMahonKD. Keystone taxa predict compositional change in microbial communities. Environ Microbiol, 2018, 20(6):2207-2217

[29] HinsingerP, GobranGR, GregoryPJ, WenzelWW. Rhizosphere geometry and heterogeneity arising from root-mediated physical and chemical processes. New Phytol, 2005, 168(2):293-303

[30] HuangLN, KuangJL, ShuWS. Microbial ecology and evolution in the acid mine drainage model system. Trends Microbiol, 2016, 24(7):581-593

[31] HuangS, JafféPR. Identification of an Actinobacteria species capable to oxidize ammonium under iron reducing conditions. Biogeoences Discussions, 2014, 11: 12295-12321

[32] KanaparthiD, PommerenkeB, CasperP, DumontMG. Chemolithotrophic nitrate-dependent Fe (II)-oxidizing nature of actinobacterial subdivision lineage TM3. ISME J, 2013, 7(8):1582-1594

[33] KonopkaA, LindemannS, FredricksonJ. Dynamics in microbial communities: unraveling mechanisms to identify principles. ISME J, 2015, 9(7):1488-1495

[34] KrupkeA, MusatN, LaRocheJ, MohrW, FuchsBM, AmannRI, KuypersMMM, FosterRA. In situ identification and N₂ and C fixation rates of uncultivated cyanobacteria populations. Syst Appl Microbiol, 2013, 36(4):259-271

[35] LeggettJE, FrereMH. Growth and nutrient uptake by soybean plants in nutrient solutions of graded concentrations. Plant Physiol, 1971, 48(4):457-460

[36] LiJ, CooperJM, LinZA, LiYT, YangXD, ZhaoBQ. Soil microbial community structure and function are significantly affected by long-term organic and mineral fertilization regimes in the North China plain. Appl Soil Ecol, 2015, 96: 75-87

[37] LiY, WangM, ChenS. Application of N₂-fixing Paenibacillus triticisoli BJ-18 changes the compositions and functions of the bacterial, diazotrophic, and fungal microbiomes in the rhizosphere and root/shoot endosphere of wheat under field conditions. Biol Fertil Soils, 2021, 57(3):347-362

[38] LiangYT, JiangYJ, WangF, WenCQ, DengY, XueK, QinYJ, YangYF, WuLY, ZhouJ, SunB. Long-term soil transplant simulating climate change with latitude significantly alters microbial temporal turnover. ISME J, 2015, 9(12):2561-2572

[39] LiawA, WienerM. Classification and regression by randomForest. R News, 2002, 2: 18-22

[40] LiebmanM, GraefRL, NettletonD, CambardellaCA. Use of legume green manures as nitrogen sources for corn production. Renew Agric Food Syst, 2012, 27(3):180-191

[41] LingN, ZhangWW, TanSY, HuangQW, ShenQR. Effect of the nursery application of bioorganic fertilizer on spatial distribution of fusarium oxysporum f. sp niveum and its antagonistic bacterium in the rhizosphere of watermelon. Appl Soil Ecol, 2012, 59: 13-19

[42] LiuJH, ZhangML, ZhangRY, ZhuWY, MaoSY. Comparative studies of the composition of bacterial microbiota associated with the ruminal content, ruminal epithelium and in the faeces of lactating dairy cows. Microb Biotechnol, 2016, 9(2):257-268

[43] LiuJW, ZhuSQ, LiuXY, YaoP, GeTT, ZhangXH. Spatiotemporal dynamics of the archaeal community in coastal sediments: assembly process and co-occurrence relationship. ISME J, 2020, 14(6):1463-1478

[44] LiuX, HuGQ, HeHB, LiangC, ZhangW, BaiZ, WuYY, LinGF, ZhangXD. Linking microbial immobilization of fertilizer nitrogen to in situ turnover of soil microbial residues in an agro-ecosystem. Agric Ecosyst Environ, 2016, 229: 40-47

[45] LoucaS, PolzMF, MazelF, AlbrightMBN, HuberJA, O'ConnorMI, AckermannM, HahnAS, SrivastavaDS, et al.. Function and functional redundancy in microbial systems. Nat Ecol Evol, 2018, 2(6):936-943

[46] MagerDM, ThomasAD. Extracellular polysaccharides from cyanobacterial soil crusts a review of their role in dryland soil processes. J Arid Environ, 2011, 75(2):91-97

[47] MaidakBL, ColeJR, LilburnTG, ParkerCT, SaxmanPR, StredwickJM, GarrityGM, LiB, OlsenGJ, et al.. The RDP (ribosomal database project) continues. Nucleic Acids Res, 2000, 28(1):173-174

[48] MartinFM. NehlsU, DasA, NebD. Carbohydrate metabolism in ectomycorrhizal Symbiosis. Molecular Mycorrhizal Symbiosis, 2016 New York John Wiley & Sons, Inc. 159-177

[49] MaslovS, SneppenK. Population cycles and species diversity in dynamic kill-the-winner model of microbial ecosystems. Sci Rep-Uk, 2017, 7(1):39642

[50] MeiX, LiuJ, GuoZW, LiPP, BiSQ, WangY, YangY, ShenWT, WangYH, XiaoY, YangX, ZhouB, LiuH, WuS. Simultaneous p-nitrophenol and nitrogen removal in PNP wastewater treatment: comparison of two integrated membrane-aerated bioreactor systems. J Hazard Mater, 2019, 363: 99-108

[51] MendesLW, KuramaeEE, NavarreteAA, van VeenJA, TsaiSM. Taxonomical and functional microbial community selection in soybean rhizosphere. ISME J, 2014, 8(8):1577-1587

[52] MoYY, ZhangWJ, YangJ, LinYS, YuZ, LinSJ. Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes. ISME J, 2018, 12(9):2198-2210

[53] MurphyMA, EvansJS, StorferA. Quantifying Bufo boreas connectivity in Yellowstone National Park with landscape genetics. Ecology, 2010, 91(1):252-261

[54] NazariesL, SinghBP, SarkerJR, FangYY, KleinM, SinghBK. The response of soil multi-functionality to agricultural management practices can be predicted by key soil abiotic and biotic properties. Agric Ecosyst Environ, 2021, 307: 107206

[55] NelsonMB, MartinyAC, MartinyJBH. Global biogeography of microbial nitrogen-cycling traits in soil. P Natl Acad Sci USA, 2016, 113(29):8033-8040

[56] OrtezOA, SalvagiottiF, EnricoJM, PrasadPVV, ArmstrongP, CiampittiIA. Exploring nitrogen limitation for historical and modern soybean genotypes. Agron J, 2018, 110(5):2080-2090

[57] Perez-GarciaO, LearG, SinghalN. Metabolic network modeling of microbial interactions in natural and engineered environmental systems. Front Microbiol, 2016, 7: 673

[58] PhilippotL, RaaijmakersJM, LemanceauP, van der PuttenWH. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol, 2013, 11(11):789-799

[59] RenLJ, HeD, ChenZ, JeppesenE, LauridsenTL, SondergaardM, LiuZW, WuQL. Warming and nutrient enrichment in combination increase stochasticity and beta diversity of bacterioplankton assemblages across freshwater mesocosms. ISME J, 2017, 11(3):613-625

[60] RogersEE, WuXL, StaceyG, NguyenHT. Two MATE proteins play a role in iron efficiency in soybean. J Plant Physiol, 2009, 166(13):1453-1459

[61] RognesT, FlouriT, NicholsB, QuinceC, MaheF. VSEARCH: a versatile open source tool for metagenomics. Peerj, 2016, 4: e2584

[62] SankariS, O'BrianMR. The Bradyrhizobium japonicum ferrous Iron transporter FeoAB is required for ferric Iron utilization in free living aerobic cells and for Symbiosis. J Biol Chem, 2016, 291(30):15653-15662

[63] SchmidtJE, KentAD, BrissonVL, GaudinACM. Agricultural management and plant selection interactively affect rhizosphere microbial community structure and nitrogen cycling. Microbiome, 2019, 7(1):146

[64] ShuDT, GuoJH, ZhangBG, HeYL, WeiGH. rDNA- and rRNA-derived communities present divergent assemblage patterns and functional traits throughout full-scale landfill leachate treatment process trains. Sci Total Environ, 2019, 646: 1069-1079

[65] SinghBP, GuptaVK, PassariAK. YadavAN, VermaP, KumarS, KumarV, KumarM, SugithaTCK, SinghBP, SaxenaAK, DhaliwalHS. Actinobacteria from rhizosphere: molecular diversity, distributions, and potential biotechnological applications. New and future developments in microbial biotechnology and bioengineering, 2018 Amsterdam Elsevier B. V 13-41

[66] Soil-Survey-Staff. Keys to Soil Taxonomy, 2010 Washington, DC USDA-Natural Resources Conservation Service

[67] SomanC, LiD, WanderMM, KentAD. Long-term fertilizer and crop-rotation treatments differentially affect soil bacterial community structure. Plant Soil, 2017, 413(1-2):145-159

[68] StricklandMS, LauberC, FiererN, BradfordMA. Testing the functional significance of microbial community composition. Ecology, 2009, 90(2):441-451

[69] Sugiyama A, Yazaki K (2012) Root exudates of legume plants and their involvement in interactions with soil microbes. Secretions and exudates in biological systems. Springer, pp 27–48. https://doi.org/10.1007/978-3-642-23047-9_2

[70] SunRB, ZhangXX, GuoXS, WangDZ, ChuHY. Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw. Soil Biol Biochem, 2015, 88: 9-18

[71] TanakaN, KubotaF, ArimaS, TaguchiM. Quantitative analysis of root system in soybean plant (Glycine max (L.)) (Merrill). JPN J Crop Sci, 2008, 59(2):270-276

[72] TaoR, LiangYC, WakelinSA, ChuGX. Supplementing chemical fertilizer with an organic component increases soil biological function and quality. Appl Soil Ecol, 2015, 96: 42-51

[73] ThebaultE, FontaineC. Stability of ecological communities and the architecture of mutualistic and trophic networks. Science, 2010, 329(5993):853-856

[74] van der HeijdenMGA, BardgettRD, van StraalenNM. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett, 2008, 11(3):296-310

[75] WarshanD, EspinozaJL, StuartRK, RichterRA, KimSY, ShapiroN, WoykeT, KyrpidesNC, BarryK, et al.. Feathermoss and epiphytic Nostoc cooperate differently: expanding the spectrum of plant-cyanobacteria symbiosis. ISME J, 2017, 11(12):2821-2833

[76] WickhamH. The split-apply-combine strategy for data analysis. J Stat Softw, 2011, 40(1):1-29

[77] WuXJ, PanJ, LiM, LiY, BartlamM, WangYY. Selective enrichment of bacterial pathogens by microplastic biofilm. Water Res, 2019, 165: 114979

[78] XueYY, ChenHH, YangJR, LiuM, HuangBQ, YangJ. Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom. ISME J, 2018, 12(9):2263-2277

[79] YaoMY, EllingFJ, JonesC, NomosatryoS, LongCP, CroweSA, AntoniewiczMR, HinrichsKU, MarescaJA. Heterotrophic bacteria from an extremely phosphate-poor lake have conditionally reduced phosphorus demand and utilize diverse sources of phosphorus. Environ Microbiol, 2016, 18(2):656-667

[80] YaoT, ChenR, ZhangJ, FengY, HuangM, LinX. Divergent patterns of microbial community composition shift under two fertilization regimes revealed by responding species. Appl Soil Ecol, 2020, 154: 103590

[81] YeGP, BanerjeeS, HeJZ, FanJB, WangZH, WeiXY, HuHW, ZhengY, DuanCJ, WanS, ChenJ, LinY. Manure application increases microbiome complexity in soil aggregate fractions: results of an 18-year field experiment. Agric Ecosyst Environ, 2021, 307: 107249

[82] ZenglerK, PalssonBO. A road map for the development of community systems (CoSy) biology. Nat Rev Microbiol, 2012, 10(5):366-372

[83] ZhangX, WangZ, JiangC, XuJ-R. Regulation of biotic interactions and responses to abiotic stresses by MAP kinase pathways in plant pathogenic fungi. Stress Biol, 2021, 1(1):5

[84] ZhaoMX, XueK, WangF, LiuSS, BaiSJ, SunB, ZhouJZ, YangYF. Microbial mediation of biogeochemical cycles revealed by simulation of global changes with soil transplant and cropping. ISME J, 2014, 8(10):2045-2055

[85] ZhaoYN, MaoXX, ZhangMS, YangW, DiHJ, MaL, LiuWJ, LiBW. The application of bacillus Megaterium alters soil microbial community composition, bioavailability of soil phosphorus and potassium, and cucumber growth in the plastic shed system of North China. Agric Ecosyst Environ, 2021, 307: 107236

[86] ZhongWH, GuT, WangW, ZhangB, LinXG, HuangQR, ShenWS. The effects of mineral fertilizer and organic manure on soil microbial community and diversity. Plant Soil, 2010, 326(1-2):511-522

[87] ZhongYQW, LiuJ, JiaXY, ShangguanZP, WangRW, YanWM. Microbial community assembly and metabolic function during wheat straw decomposition under different nitrogen fertilization treatments. Biol Fertil Soils, 2020, 56(5):697-710

Funding

National Natural Science Foundation of China(42077049)