生长相关的大黄鱼（Larimichthys crocea）肠道微生物群生态重组：微生物持续性（persistence）与生态位宽度（niche breadth）的再分配

《Biology》：Growth-Associated Ecological Reorganization of the Gut Microbiota in Larimichthys crocea: Shifts in Microbial Persistence and Niche Breadth

　　

肠道微生物群落（gut microbial communities）在鱼类发育过程中经历持续重组，但幼鱼后期生长阶段中微生物持续性（persistence）和生态位策略（ecological niche strategies）的变化仍认识不足。在此，研究人员表征了养殖大黄鱼（Larimichthys crocea）在三个连续体尺寸等级中的肠道细菌群落（gut bacterial communities）：EGC（5–8 g）、MGC（39–46 g）和LGC（160–162 g）。三个阶段依次采自同一养殖群体，并维持在同一个养殖池中。在每个阶段，采样15尾鱼，并将3个个体的肠道样品混合以生成1个混合肠道样品（pooled intestinal sample），每个阶段产生5个混合样品。鱼在整个养殖期间接受相同饲料配方，而投喂量根据生长进行调整。稀释分析（rarefaction analysis）表明测序覆盖度充分，并且ASV水平（amplicon sequence variants, ASVs）的群落结构在生长阶段间存在显著差异（ANOSIM，R = 0.4698，p = 0.001）。Pseudomonadota在整个生长过程中保持优势门，但在更精细的分类分辨率上出现了显著重构。EGC群落显示出相对一致的组成，以Methylobacterium为主，而MGC和尤其是LGC在优势细菌谱系之间表现出显著再分配。阶段特异性成员也随生长变化：特有属从EGC的36个增加到MGC的48个和LGC的50个，而阶段特异性ASV呈非线性模式并在MGC达到峰值（EGC、MGC和LGC中分别为118、173和124个特有ASV）；88个ASV在三个stage共享。持久性分类群（persistent taxa）在数量上仍重要，分别在EGC、MGC和LGC中占检出分类群的36.34%、31.59%和36.95%，但其对总群落丰度的贡献从90.30%下降到75.29%和68.06%。相反，间歇性分类群（intermittent taxa）的丰度贡献从6.80%增加到14.85%和29.17%。生态位宽度（niche-breadth）分析进一步显示，特化者（specialists）占已分类分类群的74.38%，广生者（generalists）占19.26%。尽管广生者在数量上仍占优势，其贡献从EGC到LGC下降，而特化者贡献增加。生长相关分类群进一步表现出相反轨迹，特别是Methylobacterium逐渐下降和Aquabacterium增加。总体而言，这些发现揭示了大黄鱼（L. crocea）在连续体尺寸和年龄相关采样阶段中肠道微生物群的显著生态重组。由于在单一养殖池中体尺寸和年龄与按时间顺序采样时间共变，这些模式应被解释为生长相关的时间性重组，而不是仅归因于宿主生长的效应。

一、研究背景与问题

鱼类胃肠道支持复杂微生物群落，这些微生物参与营养转化、肠道发育、免疫调节、定植抵抗及宿主对环境和饲料变化的响应。与陆生脊椎动物不同，鱼类持续与养殖水体、饲料、沉积物及其他水生环境组分来源的微生物相互作用。但肠道群落并非环境微生物池的被动反映，宿主生理和肠道生境条件可选择保留特定细菌种群。已有比较分析显示，宿主生态、饲料、发育阶段和环境条件共同影响肠道微生物群落结构。宿主生长是鱼类微生物组变异的重要维度，个体发育伴随肠道形态、消化能力、免疫发育、摄食、营养需求和肠道理化条件变化，从而改变微生物可利用的生态机会。Danio rerio、Ictalurus punctatus、Salmo salar和Gadus morhua等均报道了明显发育演替，但鱼类肠道微生物组变化不一定单向趋于稳定，生态组装过程及各阶段环境过滤重要性可能变化。核心微生物组研究表明，重复出现的细菌种群可跨宿主、环境、饲料或时间持续。生态位宽度提供互补生态维度：广生者可占据更广生态背景，特化者则与更窄的实际生态位相关。大黄鱼（L. crocea）是中国重要海水养殖物种，已有研究显示其肠道细菌组成存在明显年龄相关变化，近期整合转录组—微生物组分析揭示消化系统发育期间阶段特异性微生物演替，饲料干预也可改变其肠道微生物群。然而，在饲料配方保持不变时，后期生长阶段微生物占据和生态位策略如何变化仍知之甚少。仅凭组成周转不能揭示生长是否改变持久性群落成员的优势度，或改变微生物广生者与特化者之间的平衡。因此，研究人员开展本研究。

二、研究开展、结论与意义

研究人员在相同养殖群体、同一养殖池中连续检查三个体尺寸阶段：EGC（5–8 g）、MGC（39–46 g）和LGC（160–162 g）。三阶段为连续的体尺寸和年龄相关采样阶段，而非独立实验处理。鱼持续接受相同配方饲料，投喂量按鱼体大小和摄食需求调整。16S rRNA基因扩增子测序显示，ASV水平群落结构在生长阶段间显著分化（ANOSIM，R = 0.4698，p = 0.001）。Pseudomonadota全程为优势门，但精细分类水平发生显著重构。持久性分类群数量上仍重要，但其对总群落丰度贡献从90.30%降至68.06%；间歇性分类群丰度贡献从6.80%升至29.17%。特化者在数量上占多数，并在后期贡献增加；广生者仍数量上占优势，但贡献下降。研究意义在于揭示大黄鱼肠道微生物群沿连续体尺寸和年龄相关采样阶段发生显著生态重组，而非简单微生物稳定化；微生物持续性（persistence）和生态位宽度（niche breadth）提供描述大黄鱼幼鱼后养殖期肠道微生物组变异的互补生态维度。由于体尺寸、宿主年龄、养殖时长、投喂量和按时间顺序采样时间在单一养殖池中共同变化，观察到的微生物组重组不能因果归因于体质量本身。

三、关键技术方法

研究队列来自同一养殖群体、单一养殖池，连续三个体尺寸阶段EGC、MGC、LGC；每阶段15尾鱼，3尾肠道组织混合为1个混合肠道样本，每阶段5个，共15个。采用16S rRNA基因V4区扩增子测序（Illumina MiSeq PE300），经QIIME 2和DADA2获得ASV；用Bray–Curtis距离、主坐标分析（PCoA）和相似性分析（ANOSIM）评估群落结构；Venn分析共享与阶段特异分类群；按跨样本出现率将ASV分为瞬时（transient）、间歇（intermittent）和持久（persistent）；用Levins生态位宽度指数分为特化者（specialists）、中性者（neutralists）和广生者（generalists）；用Kruskal–Wallis H检验探索属和ASV差异。样本队列来源为同一养殖群体单一养殖池，伦理批准号ND-010-2025。

四、研究结果

3.1 测序覆盖度与生长相关肠道细菌群落分化：15个混合肠道样本获得909,832条高质量序列，平均60,655 ± 20,006条；稀释曲线和Good’s coverage表明测序覆盖充分。ASV水平主坐标分析中PC1和PC2共解释54.13%变异，ANOSIM显示EGC、MGC、LGC显著差异（R = 0.4698，p = 0.001），为阶段相关群落分化。研究未进行正式alpha多样性推断。

3.2 保守门水平框架掩盖精细分类重构：Pseudomonadota全程优势；EGC门水平相似，MGC次优势门增加，LGC再分配最大，包括Thermodesulfobacteriota、Bacillota等。纲、目、科水平亦显示后期优势谱系配置多样，如Beijerinckiaceae、Comamonadaceae、Desulfovibrionaceae、Mycoplasmataceae等变化。

3.3 属水平组成揭示优势成员强再分配：EGC中Methylobacterium占约一半或更多，Pseudomonas等亦常见；MGC中Methylobacterium优势下降、Aquabacterium增加；LGC多谱系扩展，Aquabacterium、Desulfovibrionaceae相关类群、Malacoplasma、Lawsonia等在不同池样富集，非单一替代。

3.4 阶段特异性成员在属水平增加但ASV水平非线性：门13共享；科45共享，独特科20/22/31；属60共享，独特属36/48/50；ASV 88共享，独特ASV 118/173/124，MGC峰值。说明阶段特异性不呈单调轨迹。

3.5 持久性分类群保持数量代表性但失去丰度优势：EGC transient/intermittent/persistent数量30.04%、33.61%、36.34%，丰度2.90%、6.80%、90.30%；MGC数量35.74%、32.67%、31.59%，丰度9.86%、14.85%、75.29%；LGC数量32.83%、30.22%、36.95%，丰度2.77%、29.17%、68.06%。关键变化是持久类群丰度优势下降，间歇类群丰度贡献上升。

3.6 特化者数量占优并增加后期贡献：全数据集特化者74.38%，广生者19.26%，其余为中性者。EGC广生者约95%丰度；MGC广生者下降、特化者增加；LGC广生者约80%，特化者接近五分之一。持久性下降与特化者增加平行，但二者独立分类，不能证明持久者即广生者或间歇者即特化者。

3.7 优势属呈现相反生长相关轨迹：Kruskal–Wallis识别十个名义差异属：Methylobacterium、Aquabacterium、Pseudomonas、Acinetobacter、Burkholderia-Caballeronia、Agrobacterium、Methyloversatilis、Delftia、Stenotrophomonas、Sphingobium。Methylobacterium下降，Aquabacterium上升；Pseudomonas等下降；Stenotrophomonas复杂。ASV水平十个差异ASV，ASV2最多且下降，ASV1上升，ASV8在MGC最高。说明多轨迹同时发生。

五、讨论总结

讨论认为，生长伴随显著但非线性肠道微生物演替，MGC可能是过渡性微生物配置，但未对应正式生理成熟阶段。门水平稳定不意味稳定，精细分类重构明显。持久性基于跨样本出现，非个体纵向持续；持久类群失去数量优势而非消失。特化者贡献增加表示窄生态位类群数量代表性增加，但未实验证明肠道生态位复杂性增加。优势类群相反轨迹显示精细尺度替换。16S rRNA基因丰度不能推断功能，差异类群仅为候选指示。设计局限包括：单池、体尺寸/宿主年龄/养殖时长/按时间顺序采样时间共变；每阶段5个混合样为池内生物学亚样，非独立池重复；统计功效有限；Kruskal–Wallis探索性且未多重校正；环境指标未测；投喂量随体尺寸变化；混合样无法评估个体间变异。未来需重复养殖单元、同步环境监测和宿主生理测量。

六、研究结论翻译

大黄鱼（L. crocea）肠道微生物群在三个连续体尺寸和年龄相关采样阶段中表现出显著生态重组。尽管Pseudomonadota仍为优势门，但在更精细分类分辨率上发生了显著周转，包括Methylobacterium优势逐渐降低以及若干替代细菌谱系代表性增加。微生物成员呈非线性变化，阶段特异性属向后期采样阶段增加，而特有ASV在MGC达到峰值，表明中间采样点存在尤其显著的精细尺度周转。然而，由于生长阶段根据体尺寸和养殖时长而非组织学验证的发育或生理成熟标准定义，MGC不应被解释为正式定义的生理过渡阶段。持久性ASV在整个数据集中数量代表性仍良好，但其对总群落丰度的贡献从EGC的90.30%降至LGC的68.06%，而间歇性ASV从6.80%增至29.17%。同时，广生者在数量上仍占优势，而操作分类的特化者贡献在后期阶段群落中增加。总之，这些结果证明养殖大黄鱼连续采样阶段中微生物优势度、出现模式和推断生态位策略发生显著再分配。由于在单一池中宿主年龄和体尺寸与按时间顺序采样时间共变，这些模式应解释为生长相关的时间性微生物组重组，而不是仅由宿主生长引起效应的证据。尽管如此，微生物出现率和生态位宽度为描述大黄鱼幼鱼后养殖期肠道微生物组变异提供了互补生态维度。由于在本单池设计中体尺寸、宿主年龄、养殖时长、投喂量和按时间顺序采样时间共变，观察到的微生物组重组不能因果归因于体质量本身。