《Environmental Microbiology》：Functional Metagenomics Reveals Gut Microbial Biosynthetic Potential for Essential Nutrients in Juvenile Chum Salmon
编辑推荐：
鱼类肠道微生物群有助于宿主营养、免疫和稳态，然而其代谢能力与宿主基因组缺口之间的关系仍不清楚。本研究采用鸟枪法宏基因组测序（shotgun metagenomic sequencing）和KEGG模块完成度分析（KEGG module completion a
鱼类肠道微生物群有助于宿主营养、免疫和稳态，然而其代谢能力与宿主基因组缺口之间的关系仍不清楚。本研究采用鸟枪法宏基因组测序（shotgun metagenomic sequencing）和KEGG模块完成度分析（KEGG module completion analysis），比较了幼年大麻哈鱼（Oncorhynchus keta）在孵化后第35天（即开始外源性摄食时）肠道微生物组与宿主基因组的生物合成潜力。由于幼鱼肠道体积小，微生物生物量低且宿主DNA比例高，研究人员首先优化了一种结合酶解裂解与基于NEBNext的宿主DNA去除方案。该方法将细菌序列占比从11.3%提高至41.2%，并将映射到KEGG模块的KEGG直系同源（KO）标识符数量从1128增加到1532。模块完成率（Module Completion Ratio, MCR）分析揭示了功能潜力上的互补性：大麻哈鱼基因组在必需氨基酸（苏氨酸、亮氨酸、缬氨酸/异亮氨酸、赖氨酸）和B族维生素（硫胺素、核黄素、生物素）的从头生物合成方面不完整，而肠道宏基因组编码了这些途径的完整或近乎完整的模块。归属于Bartonella样、Acinetobacter和Enterococcus谱系的序列贡献最大。这些观察表明肠道微生物组编码了宿主基因组所缺乏的生物合成潜力，但证明其对宿主的营养益处仍需实验操作。
肠道微生物群被誉为鱼类的“额外器官”，在营养供给、免疫成熟和代谢稳态中发挥广泛作用。然而，微生物组的代谢能力如何弥补宿主基因组自身无法合成的必需营养素，尤其是在幼鱼发育关键期，此前缺乏直接的功能证据。大麻哈鱼（Oncorhynchus keta）是研究这一问题的理想模型：它们比大多数太平洋鲑鱼更早、更小地进入海洋，孵化场释放的幼鱼质量直接关系到渔业资源补充。在孵化后第35天，仔鱼卵黄囊耗尽并转为外源性摄食，肠道首次接触外部底物与微生物，此时生长需求极高，而宿主基因组又无法从头合成若干必需氨基酸和B族维生素，因此肠道共生微生物可能承担弥补这些代谢缺口的任务。
为验证这一假说，研究人员对12尾35日龄幼年大麻哈鱼的全肠道匀浆进行混合取样，并比较四种DNA预处理策略：原始DNA、NEBNext宿主DNA去除、NEBNext去除后17轮PCR扩增、以及NEBNext去除后全基因组扩增（WGA）。优化后的方法（酶解裂解结合NEBNext甲基化依赖性宿主DNA去除）将细菌序列比例从11.3%提升至41.2%，并将可映射到KEGG模块的KO编号从1128增至1532。基于此，研究人员利用基因中心的KEGG模块完成率（MCR）分析方法，将大麻哈鱼宿主基因组（117,599条蛋白质编码序列）与肠道宏基因组的功能潜力在同一标准下进行比较。结果显示：宿主基因组在苏氨酸、亮氨酸、缬氨酸/异亮氨酸、赖氨酸以及硫胺素（B1）、核黄素（B2）、生物素（B7）的从头生物合成途径上不完整，而肠道宏基因组编码了这些途径的完整或近乎完整模块。序列归属分析表明，贡献主要来自Bartonella样、Acinetobacter和Enterococcus等谱系。这些发现首次从功能模块层面定量证明，幼年大麻哈鱼肠道微生物组编码了宿主基因组所缺乏的生物合成潜力，为理解微生物对早期发育的营养支持提供了直接证据，也为孵化场饲料优化和益生菌策略提供了微生物学框架。
研究人员为开展本研究主要采用了以下关键技术方法：基于酶解裂解与NEBNext甲基化结合蛋白的宿主DNA去除流程，以富集低生物量样本中的微生物序列；采用基因中心的宏基因组处理流程（Genomaple Submission Data Maker, GSDM）进行序列修剪、基因预测和氨基酸序列翻译；使用Genomaple version 2.3.2进行KEGG直系同源（KO）注释、模块映射和MCR计算；基于保守核糖体蛋白基因进行微生物群落结构谱分析。样本来源于日本岩手县水产技术中心提供的同一人工受精事件的大麻哈鱼胚胎，孵化后35天采集12尾个体全肠道（含内容物）混合为一个匀浆，再分为四份进行不同预处理。
3.1 测序结果与预处理方法比较：研究人员比较了四种DNA预处理策略。未处理原始DNA库产生653,543条序列，模块映射KO数为1128；NEBNext单独处理产生998,957条序列，模块KO数为1499；NEBNext联合17轮PCR产生999,458条序列，模块KO数最高（1532）；而WGA处理仅产生338,389条序列和666个模块KO。稀疏散点图表明宿主DNA去除显著提高了功能注释深度。
3.2 宿主DNA去除使细菌序列比例提升近四倍：通过Genomaple注释将各文库序列按来源域划分，未处理文库中仅11.3%为细菌序列，88.7%为真核（宿主）序列；NEBNext处理后细菌比例升至41.2%，富集3.6倍；两种扩增文库分别为35.4%和35.7%，未进一步改善。去除后仍残留约60%真核序列，说明甲基化法不能完全消除宿主DNA。同时，去除过程使DNA浓度从650.8降至28.1 ng/μL，约4%起始质量留存。
3.3 大麻哈鱼肠道微生物群落结构：基于核糖体蛋白基因的组成分析（宿主序列去除后）显示，α-变形菌纲和γ-变形菌纲持续占优。NEBNext处理文库中α-变形菌纲占80.5%（NEBNext+PCR为77.1%），而未处理文库为55.6%。物种水平上，检测到Mycoplasma亲和谱系（最接近参考基因组为Mycoplasma penetrans），但该鉴定仅为同源搜索最佳匹配而非ANI验证。低丰度类群在不同处理间不稳定，WGA文库偏差明显，因此低丰度类群仅作定性解读。基于功能深度，后续代谢分析选用NEBNext+17轮PCR文库。
3.4 大麻哈鱼基因组与肠道微生物组的功能组比较：MCR比较显示，糖酵解、三羧酸循环和核心核苷酸代谢在两者中均完整；而宿主基因组在若干必需氨基酸和B族维生素生物合成模块上MCR低或缺失，肠道宏基因组则高。此外还发现半胱氨酸生物合成和植物/细菌型血红素生物合成途径（与宿主自身动物/真菌型途径共存）存在类似互补。
3.5 必需氨基酸生物合成：苏氨酸模块（M00018）在宿主基因组中不完整，在宏基因组中完整，主要贡献来自Bartonella apis（20%）、Aureimonas sp. AU20（15%）、Enterococcus faecalis V583（10%）和Pseudovibrio sp. FO-BEG1（10%）。亮氨酸模块（M00432）由Erwinia amylovora CFBP1430（33%）、Enterococcus faecium Aus0004（17%）等贡献。缬氨酸/异亮氨酸模块（M00019）由Beijerinckia indica（17%）、B. apis（9%）等贡献。赖氨酸模块（M00016）经琥珀酰-DAP途径，由B. apis（33%）、Aureimonas sp. AU20（14%）等完成。所有物种名称均指最接近参考基因组而非ANI验证的物种鉴定。
3.6 必需维生素生物合成：硫胺素模块（M00127）由Acinetobacter calcoaceticus样序列（约28%）、Bartonella sp. JB63（14%）等贡献。核黄素模块（M00125）由A. calcoaceticus（11%）、Bartonella australis（11%）等贡献。生物素模块（M00123）在宏基因组中的MCR高于宿主基因组，但主要匹配为胞内内共生菌（Wolbachia、Candidatus Fokinia solitaria、Cardinium），研究人员对此持保留态度，仅报告模块存在而不指定具体类群。
在讨论部分，研究人员指出，酶解裂解与甲基化宿主DNA去除的结合使幼鱼肠道宏基因组测序变得可行，该方法可推广至其他低生物量鱼类样本。功能互补结果与脊椎动物普遍依赖共生微生物补充必需营养素的观点一致，但模块完整仅代表遗传能力，不证明转录或宿主吸收，需通过无菌模型或同位素示踪进一步验证。不同发育阶段鲑科鱼类优势菌群不同，本研究的Bartonella样、Acinetobacter和Enterococcus为主导类群，而成年鲑鱼常见Mycoplasma优势，提示功能收敛而非特定类群依赖。生物素模块的内共生菌匹配提示参考基因组数据库偏差，需通过基因组分箱和ANI验证纠正。研究局限性包括：混合样本无生物学重复、缺乏统计推断、物种分配依赖数据库覆盖、MCR仅代表潜力。结论方面，研究人员开发了适用于幼年大麻哈鱼极小肠道的基因中心宏基因组工作流，并比较了35日龄宿主基因组与肠道微生物组的代谢能力，发现必需氨基酸和B族维生素合成途径在宿主中不完整而在微生物组中高度编码。这些观察为孵化场微生物组意识策略（如调整饲料或使用益生菌）提供了依据，并强调了在海洋环境变化下理解野生鲑鱼微生物营养作用的重要性。
