高级检索+

稻蟹共作水体细菌群落的生态稳定性维持机制及其环境驱动解析

Environmental drivers and mechanisms maintaining the ecological stability of bacterial communities in rice-crab co-culture systems

  • 摘要: 稻蟹共作水体微生态稳定性直接关系到养分循环、水质维持和养殖安全。为阐明全养殖周期水体细菌群落的稳态维持机制及环境驱动,本研究于4—8月连续采集5个关键阶段水样,采用16S rRNA高通量测序,结合α、β多样性、中性群落模型(neutral community model,NCM)、改进随机性比率(modified stochasticity ratio,MST)、PER-SIMPER、扩散—生态位连续体指数(Dispersal–Niche Continuum Index,DNCI)、共现网络和结构方程模型(structural equation model,SEM),定量解析群落时序演替、组装过程及网络稳定性。结果表明,不同月份细菌α和β多样性均存在显著差异(P<0.05),变形菌门、拟杆菌门和放线菌门合计相对丰度超过60%;群落差异主要由物种周转驱动,各月周转贡献率均高于80%,嵌套贡献率仅为3.3%~12.5%。NCM拟合度由4月的0.481升至6月的0.678,全周期R2为0.651;MST由4月的44.65%升至5—8月的54.93%~67.89%,表明群落组装由养殖初期确定性过程主导转变为中后期随机过程主导,PER-SIMPER和DNCI结果与之相互印证。共现网络在6月连接最紧密,平均度、图密度和平均聚类系数分别达到29.817、0.229和0.642,综合复杂性指数和稳定性指数分别由4月的5.02和16.67升至73.95和83.33。SEM拟合良好(X2=9.586,P=0.388,RMSEA=0.000),环境因子与组装过程共同解释网络稳定性59.0%的变异;亚硝态氮(NO2-N)和总氮(TN)对稳定性呈正向作用(λ=0.559、0.475),正磷酸盐(PO4-P)显著抑制随机性组装(λ=-0.869),随机性组装则正向促进稳定性(λ=0.397)。研究表明,随机过程增强及氮磷因子的协同调控是维持稻蟹共作水体微生态稳态的关键,可为西北干旱区稻蟹系统控磷稳氮、水质调控及健康养殖提供理论依据。

     

    Abstract: The ecological stability of aquatic microecosystems in rice–crab co-culture systems is essential for regulating nutrient cycling, maintaining water quality, and sustaining the long-term functioning of integrated agricultural ecosystems. However, the mechanisms underlying bacterial community stability under seasonal environmental fluctuations and intensive aquaculture management remain poorly understood. This study investigated the temporal dynamics of bacterial community succession, assembly mechanisms, co-occurrence network characteristics, and environmental regulatory pathways throughout the rice–crab cultivation cycle. Water samples were collected at five key growth stages from April to August, and bacterial communities were characterized using 16S rRNA gene high-throughput sequencing. Multiple ecological approaches, including α and β diversity analyses, the neutral community model (NCM), modified stochasticity ratio (MST), PER-SIMPER, dispersal–niche continuum index (DNCI), co-occurrence network analysis, and structural equation modeling (SEM), were integrated to quantify bacterial community assembly and ecological stability. The results showed significant temporal variations in bacterial α and β diversity during the cultivation period (P < 0.05). Proteobacteria, Bacteroidota, and Actinobacteriota were the dominant phyla, collectively accounting for more than 60% of total bacterial abundance. Community dissimilarity was mainly driven by species turnover, which contributed more than 80% of β diversity in all sampling months, whereas nestedness accounted for only 3.3%–12.5%, indicating rapid bacterial species replacement under changing environmental conditions. NCM analysis revealed that bacterial communities generally followed neutral assembly patterns, with model explanatory power increasing from R2 = 0.481 in April to R2 = 0.678 in June, and the overall model explaining 65.1% of community variation (R2 = 0.651). Meanwhile, MST values increased from 44.65% in April to 54.93%, 67.89%, 57.64%, and 64.29% from May to August, respectively, demonstrating a transition from deterministic environmental filtering during the early stage to stochastic processes dominating during the middle and late stages. PER-SIMPER and DNCI analyses further confirmed the increasing contribution of stochastic processes, particularly dispersal and ecological drift, to bacterial community assembly. Co-occurrence network analysis revealed pronounced temporal restructuring of microbial interactions. Network complexity and stability were relatively low during the early cultivation stage but increased substantially during the middle and late stages. In June, the bacterial network exhibited the strongest connectivity, with average degree, graph density, and average clustering coefficient reaching 29.817, 0.229, and 0.642, respectively, while average path length decreased to 2.315. Correspondingly, the comprehensive network complexity and stability indices increased from 5.02 and 16.67 in April to maximum values of 73.95 and 83.33 in June, indicating enhanced microbial interactions and improved ecological resilience. SEM further revealed the relationships among environmental factors, community assembly processes, and ecological stability. The model showed excellent fitting performance (X2 = 9.586, P = 0.388, RMSEA = 0.000), with environmental variables and assembly processes jointly explaining 59.0% of bacterial network stability variation. Nitrite nitrogen and total nitrogen positively influenced network stability, with standardized path coefficients of 0.559 and 0.475, respectively, whereas orthophosphate phosphorus strongly inhibited stochastic assembly (λ = −0.869). In contrast, stochastic assembly positively promoted network stability (λ = 0.397), highlighting its role as an intrinsic ecological buffering mechanism. Overall, this study demonstrates that enhanced stochastic assembly processes, coordinated nitrogen–phosphorus regulation, and microbial network restructuring jointly maintain ecological stability in rice–crab co-culture aquatic ecosystems. These findings provide new insights into microbial resilience mechanisms and offer quantitative guidance for nutrient management, water-quality regulation, and sustainable aquaculture practices in arid and environmentally sensitive regions.

     

/

返回文章
返回