The crayfish-rice coculture model contributes to regulating the soil fertility of rice fields and maintaining the stability of soil microbial community composition and function

Abstract

Rice-fish coculture represents a classic sustainable agricultural paradigm; however, the microecological mechanisms underlying its capacity to maintain soil fertility and microbial community stability remain poorly understood. We conducted a 13-month field experiment comparing three cultivation systems:crayfish-rice coculture (CRCE), crayfish-waterweed coculture (CWCE), and rice monoculture (RME)-by integrating physicochemical analysis, 16S rRNA sequencing, metagenomics, microbial network analysis, and null model simulations. Our results demonstrated that coculture systems, particularly CRCE, enhanced soil fertility through carbon sequestration (total carbon: 25.0–45.0 mg/g; total organic carbon: 15.0–35.0 mg/g) and sustained redox homeostasis (consistently low oxidation–reduction potential: − 150 to − 50 mV), in stark contrast to the extreme redox fluctuations observed in RME. These stable edaphic conditions imposed deterministic selection on microbial communities (homogeneous selection contribution: 30%–50% in CRCE vs. 10%–20% in RME), shifting community assembly from stochastic drift dominance toward predictable succession. This assembly shift enriched functionally coupled keystone taxa, including iron reducers (Geobacter), sulfur oxidizers (Sulfuricurvum), and nitrifiers (Nitrospira), which formed ecological networks characterized by 98.6% positive interactions and enhanced functional gene repertoires associated with carbon, nitrogen, and sulfur biogeochemical cycles. Metagenomic analysis corroborated these findings, revealing enrichment of functional genes involved in polymer degradation, nitrification, and sulfate reduction in CRCE, supporting enhanced nutrient cycling capacity. We establish a hierarchical causal pathway in which bioturbation-induced environmental stabilization drives deterministic community assembly, which in turn promotes keystone taxon enrichment and functional integration. This framework provides a mechanistic explanation for how crayfish-rice coculture regulates soil fertility and sustains microbial community compositional and functional stability in anthropogenically designed agricultural ecosystems.

Keywords

Rice-crayfish coculture / Soil fertility / Microbial community stability / Community assembly / Co‑occurrence network / Functional metagenomics

References

[1] Achenbach LA, Michaelidou U, Bruce RA, Fryman J, Coates JD. Dechloromonas agitata gen. nov., sp. nov. and Dechlorosoma suillum gen. nov., sp. nov., two novel environmentally dominant (per)chlorate-reducing bacteria and their phylogenetic position. Int J Syst Evol Microbiol, 2001, 51(Pt 2): 527-533

[2] Barbaresi S, Tricarico E, Gherardi F. Factors inducing the intense burrowing activity of the red-swamp crayfish, Procambarus clarkii, an invasive species. Naturwissenschaften. 2004;91(7):342–5. https://doi.org/10.1007/s00114-004-0533-9.

[3] Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Proceed Int AAAI Conference Web Soc Media, 2009, 3(1): 361-362

[4] Chang J, Costa OYA, Sun Y, Wang J, Tian L, Shi S, Wang E, Ji L, Wang C, Pang Y, Yao Z, Ye L, Zhang J, Chen H, Cai Y, Chen D, Song Z, Rong J, Raaijmakers JM, Tian C, Kuramae EE. Domesticated rice alters the rhizosphere microbiome, reducing nitrogen fixation and increasing nitrous oxide emissions. Nat Commun, 2025, 16(1): 2038

[5] Chen S, Guo Y, Yuan P, Jiang Y, Cao C. Carbon sequestration, emission reduction, and technical strategies of rice-crayfish farming in central China. Sci Rep, 2025, 15(1): 29352

[6] Cron B, Henri P, Chan CS, Macalady JL, Cosmidis J. Elemental sulfur formation by Sulfuricurvum kujiense is mediated by extracellular organic compounds. Front Microbiol, 2019, 10: 2710

[7] Cui J, Liu H, Wang H, Wu S, Bashir MA, Reis S, Sun Q, Xu J, Gu B. Rice-animal co-culture systems benefit global sustainable intensification. Earths Future, 2023, 11(2): e2022EF002984

[8] Dorn NJ, Wojdak JM. The role of omnivorous crayfish in littoral communities. Oecologia. 2004;140(1):150–9. https://doi.org/10.1007/s00442-004-1548-9.

[9] Duffner C, Holzapfel S, Wunderlich A, Einsiedl F, Schloter M, Schulz S. Dechloromonas and close relatives prevail during hydrogenotrophic denitrification in stimulated microcosms with oxic aquifer material. FEMS Microbiol Ecol, 2021, 97(3): fiab004

[10] Edgar R. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10. https://doi.org/10.1038/nmeth.2604

[11] Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics, 2004, 5: 113

[12] Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics, 2011, 27(16): 2194-2200

[13] FAO. Twenty years of Globally Important Agricultural Heritage Systems – Success stories of dynamic conservation for sustainable rural development. Rome. 2022. https://doi.org/10.4060/cc2385en.

[14] Faust K, Raes J. Microbial interactions: from networks to models. Nat Rev Microbiol, 2012, 10(8): 538-550

[15] Gemayel K, Lomsadze A, Borodovsky M. MetaGeneMark-2: Improved Gene Prediction in Metagenomes. bioRxiv 2022:2022.2007.2025.500264. https://doi.org/10.1101/2022.07.25.500264

[16] Guo Y, Hu Q, Liu T, Du Y, Li C, Zhang X, Liu J, Cao C. Long-term rice–crayfish coculture increases plant lignin but not microbial necromass contribution to soil organic carbon. Soil Tillage Res, 2025, 248: 106424

[17] Hartmann M, Frey B, Mayer J, Mäder P, Widmer F. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J, 2015, 9(5): 1177-1194

[18] Hoffman JIE. Chapter 25 - Analysis of Variance. I. One-Way, in: Hoffman, J.I.E. (Ed.), Basic Biostatistics for Medical and Biomedical Practitioners (Second Edition). Academic Press, 2019; pp. 391–417.

[19] Hou Y, Jia R, Sun W, Li B, Zhu J. Influences of the integrated rice-crayfish farming system with different stocking densities on the paddy soil microbiomes. Int J Mol Sci, 2024, 25(7): 3786

[20] Hu L, Zhang J, Ren W, Guo L, Cheng Y, Li J, Li K, Zhu Z, Zhang J, Luo S, Cheng L, Tang J, Chen X. Can the co-cultivation of rice and fish help sustain rice production?. Sci Rep, 2016, 6(1): 28728

[21] Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res, 2012, 40: D109-114

[22] Lehmann J, Hansel CM, Kaiser C, Kleber M, Maher K, Manzoni S, Nunan N, Reichstein M, Schimel JP, Torn MS, Wieder WR, Kögel-Knabner I. Persistence of soil organic carbon caused by functional complexity. Nat Geosci, 2020, 13(8): 529-534

[23] Leps J, Smilauer P. Multivariate analysis of ecological data. 87, 2006, 3: 193-193

[24] Li D, Luo R, Liu C-M, Leung C-M, Ting H-F, Sadakane K, Yamashita H, Lam T-W. MEGAHIT v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods, 2016, 102: 3-11

[25] Li P, Wu G, Li Y, Hu C, Ge L, Zheng X, Zhang J, Chen J, Zhang H, Bai N, Zhang H, Song L, Sun Y, Jiang W, Jia J, Chen Y, Wang C, Lv B, Wu X, Pan A, Li S, Lv W. Long-term rice-crayfish-turtle co-culture maintains high crop yields by improving soil health and increasing soil microbial community stability. Geoderma, 2022, 413: 115745

[26] Li Q, Xie L, Lin S, Cheng X, Liu Q, Li Y. Effects of rice–fish coculture on greenhouse gas emissions: a case study in terraced paddy fields of Qingtian, China. Agronomy, 2025, 15(6): 1480

[27] Liu Y, Ren W, Zhao L, Guo L, Liu Y, Hu L, Chen X. Community structure and functional features of soil microbes in the rice-fish co-culture system. J South China Agric Univ, 2024, 45(6): 865-877

[28] Loreau M, Naeem S, Inchausti P, Bengtsson J, Grime JP, Hector A, Hooper DU, Huston MA, Raffaelli D, Schmid B, Tilman D, Wardle DA. Biodiversity and ecosystem functioning: current knowledge and future challenges. Science, 2001, 294(5543): 804-808

[29] Lovley DR, Holmes DE, Nevin KP. Dissimilatory Fe(III) and Mn(IV) reduction. Adv Microb Physiol, 2004, 49: 219-286

[30] Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics, 2011, 27(21): 2957-2963

[31] McAllister SM, Vandzura R, Keffer JL, Polson SW, Chan CS. Aerobic and anaerobic iron oxidizers together drive denitrification and carbon cycling at marine iron-rich hydrothermal vents. ISME J, 2021, 15(5): 1271-1286

[32] Myers, L., Sirois, M.J., 2006. Spearman Correlation Coefficients, Differences between, in: Encyclopedia of Statistical Sciences.

[33] Ning D, Yuan M, Wu L, Zhang Y, Guo X, Zhou X, Yang Y, Arkin AP, Firestone MK, Zhou J. A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming. Nat Commun, 2020, 11(1): 4717

[34] Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin P, O’Hara RB, Simpson G, Solymos P, Stevens MHH, Wagner H. Vegan: Community Ecology Package. R Package Version. 2.0–10. CRAN. 2013.

[35] Pang Z, Xu P. Probiotic model for studying rhizosphere interactions of root exudates and the functional microbiome. ISME J, 2024, 18(1): wrae223

[36] Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res, 2015, 25(7): 1043-1055

[37] Ren W, Hu L, Zhang J, Sun C, Tang J, Yuan Y, Chen X. Can positive interactions between cultivated species help to sustain modern agriculture? Front Ecol Environ. 2014;12(9):507–14. https://doi.org/10.1890/130162.

[38] Rühlemann MC, Wacker EM, Ellinghaus D, Franke A. MAGScoT: a fast, lightweight and accurate bin-refinement tool. Bioinformatics, 2022, 38(24): 5430-5433

[39] Six J, Conant RT, Paul E, Paustian K. Stabilization Mechanisms of Soil Organic Matter: implications for C-saturation of soils. Plant Soil, 2002, 241: 155-176

[40] Smith CA, O'Maille G, Want EJ, Qin C, Trauger SA, Brandon TR, Custodio DE, Abagyan R, Siuzdak G. METLIN: a metabolite mass spectral database. Ther Drug Monit, 2005, 27(6): 747-751

[41] Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol, 2007, 73(16): 5261-5267

[42] Wang A, Hao X, Chen W, Luo X, Huang Q. Rice-crayfish co-culture increases microbial necromass' contribution to the soil nitrogen pool. Environ Res, 2023, 216: 114708

[43] Warton DI, Wright ST, Wang Y. Distance-based multivariate analyses confound location and dispersion effects. Methods Ecol Evol, 2012, 3(1): 89-101

[44] Wei D, Xing C, Hou D, Zeng S, Zhou R, Yu L, Wang H, Deng Z, Weng S, He J, Huang Z. Distinct bacterial communities in the environmental water, sediment and intestine between two crayfish-plant coculture ecosystems. Appl Microbiol Biotechnol, 2021

[45] Xie J, Hu L, Tang J, Wu X, Li N, Yuan Y, Yang H, Zhang J, Luo S, Chen X. Ecological mechanisms underlying the sustainability of the agricultural heritage rice-fish coculture system. Proc Natl Acad Sci U S A, 2011, 108(50): E1381-1387

[46] Yang P, Lai DYF, Jin B, Bastviken D, Tan L, Tong C. Dynamics of dissolved nutrients in the aquaculture shrimp ponds of the Min River estuary, China: Concentrations, fluxes and environmental loads. Sci Total Environ. 2017;603–604:256–67. https://doi.org/10.1016/j.scitotenv.2017.06.074.

[47] Zhang J, Zhang N, Liu YX, Zhang X, Hu B, Qin Y, Xu H, Wang H, Guo X, Qian J, Wang W, Zhang P, Jin T, Chu C, Bai Y. Root microbiota shift in rice correlates with resident time in the field and developmental stage. Sci China Life Sci, 2018, 61(6): 613-621

[48] Zhang Z, Du L, Xiao Z, Li C, Wang Z, Zhou P, Sun G, Ye Y, Hu T, Wang H. Rice-crayfish farming increases soil organic carbon. Agric Ecosyst Environ, 2022, 329: 107857

[49] Zhang W, Song Y, Ma S, Lu J, Zhu J, Wang J, Li X. Rice-crayfish farming system promote subsoil microbial residual carbon accumulation and stabilization by mediating microbial metabolism process. Sci Total Environ, 2024, 946: 174188

[50] Zhou J, Ning D. Stochastic Community Assembly: Does It Matter in Microbial Ecology? Microbiol Mol Biol Rev: MMBR. 2017;81(4). https://doi.org/10.1128/mmbr.00002-17

[51] Zhu X, Yang P, Xiong G, Wei H, Zhang L, Wang Z, Ning K. Microbial biogeochemical cycling reveals the sustainability of the rice-crayfish co-culture model. iScience, 2023, 26(5): 106769

Funding

National Key Research and Development Program of China(2023YFD2401705)

National Key Research and Development Program of China(2024YFD2401202)

Earmarked Fund for CARS-48-20; Guangxi Science and Technology Major Special Project(AA23062047)

Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai)(SML2021SP203)

RIGHTS & PERMISSIONS

The Author(s)