企业官网建设流程全解析

摘要

生物合成基因簇（biosynthetic gene cluster，BGC）的原位生态学意义仍不清楚。作者采用宏基因组与宏转录组联用策略，沿泥炭地氧化还原梯度表征BGC多样性与活性。基于7个宏基因组共回收9,694个BGC，分类涵盖多类微生物，其中大多数在参考数据库中缺乏近缘亲属，表明其具有高度新颖性。原位条件下仅有9%–27%的潜在BGC被转录；酸杆菌门（Acidobacteriota）尽管整体BGC储量中等，却贡献超过一半的BGC转录量。携带最多可达24个BGC的“高产区”基因组多数处于沉默状态，且BGC表达比例与每基因组BGC数量呈负相关。酸杆菌特化于复杂碳加工潜力，其BGC曲目表达比例高于BGC丰富但倍增时间更短、复杂碳潜力更弱的假单胞菌门（Pseudomonadota）。在基因组水平，BGC表达程度与碳水化合物活性酶（carbohydrate-active enzyme，CAZyme）表达，尤其是糖苷水解酶（glycoside hydrolase，GH）表达耦合最强，提示次生代谢与泥炭地微生物主动碳加工相整合。结果表明，BGC虽广泛存在，但其原位表达更多体现各分类群个体碳循环活性，而非统一的生理状态。

keyword：生物合成基因簇（BGC）；泥炭地；宏基因组；宏转录组；酸杆菌门（Acidobacteriota）；假单胞菌门（Pseudomonadota）；碳水化合物活性酶（CAZyme）；糖苷水解酶（GH）；氧化还原梯度；次级代谢

文献信息

• Human Z R, Hoover R, Chan C, Küsel K, Wegner C E（2026）. High diversity but limited expression of biosynthetic gene clusters in a peatland microbial community. Nature Communications 17, 10529. https://doi.org/10.1038/s41467-026-78121-8
• 期刊：Nature Communications（Nature Portfolio，ISSN 2041-1723，完全开放获取）
• 影响因子：2025 Journal Impact Factor 18.1，五年影响因子18.9，JCR多学科Q1（2026年JCR发布）
• 接收时间：2026年1月16日投稿，2026年9月15日接受，2026年10月6日在线发表
• 研究地点：德国Schloppnerbrunnen贫营养潜育泥炭沼（fen）
• 数据：NCBI SRA BioProject PRJNA1374247；MAG与BGC数据Zenodo https://doi.org/10.5281/zenodo.17812908

研究总结

背景介绍

次级代谢产物由BGC编码，包括金属螯合剂、抗真菌物、核糖体肽与抗生素等，具有生态竞争、胁迫防御和资源获取功能。培养研究多集中于放线菌门（Actinomycetota）与假单胞菌门，认为次生代谢常在营养限制、氧化应激或 stationary phase 激活；但该类结论来自可培养菌株，对环境微生物覆盖面有限。泥炭地是全球重要碳储，水文形成从表层好氧到深层永久缺氧的陡峭氧化还原梯度，并伴随电子受体竞争与营养变化。酸杆菌门、疣微菌门（Verrucomicrobiota）、浮霉菌门（Planctomycetota）等类群在培养中代表性不足，但已有证据显示其具有丰富BGC潜力。作者假设：泥炭地不仅是BGC新颖性储库，其氧化还原异质性会对广泛分布的BGC曲目进行选择性过滤，使分布与表达均受到环境驱动。研究以德国Schloppnerbrunnen泥炭沼为对象，联合宏基因组组装基因组（metagenome-assembled genome，MAG）与转录组，系统评估BGC多样性、原位表达及与碳代谢的关系。

重要结果

结果1 泥炭地群落组成与沿深度的氧化还原分层

作者基于k-mer覆盖估计确认测序深度可捕获群落基因组多样性；以SingleM单拷贝看家基因、CoverM比对MAG，获得573个合格MAG（完整度中位94.4%，污染2.2%）。表层以酸杆菌门、放线菌门、疣微菌门为主，随深度增加假单胞菌门、脱硫菌门（Desulfobacterota）、古菌占比上升。基于Aerobot氨基酸组成模型与Oxyphen含氧酶计数，厌氧类群比例由表层6%–34%增至25–30厘米层29%–50%。该分层为后续BGC分布与表达分析提供氧化还原背景。

图1 泥炭地不同深度微生物组成与多样性。A 基于原始读段单拷贝标记基因获得的分类剖面；B 重构MAG相对丰度。两图共同显示酸杆菌门在表层占优、深层厌氧类群增加，印证从好氧到厌氧的氧化还原过渡。

结果2 BGC总量、分类多样性与随深度下降的生物合成密度

以contig大于5 kb运行antiSMASH，共识别9,694个BGC，平均31,073 bp；类型为RiPP（1,361）、NRPS（510）、PKS（367）、萜类（620）、NRPS-PKS杂合（107）及多种RiPP子类。针对短读长组装碎片，作者仅对缺核心基因且距contig端大于5 kb的NRPS/PKS/杂合簇记为弱碎片（1,544个，15.9%），其余保留。BiG-SCAPE聚类获得大量基因簇家族（GCF），分类以酸杆菌1,404、假单胞菌754、放线菌278、疣微菌175、黏球菌（Myxococcota）164为主，古菌40个。BGC密度（每Mbp BGC数）随深度显著下降（全体Kruskal-Wallis H=16.2，P=0.00275；酸杆菌H=12.8，P=0.0125），说明深层低生物合成能力不是组装假象，而与厌氧类群增加有关。

图2 BGC多样性。A 各BGC类别及其归入GCF的数量；B 各分类与深度对应的GCF计数；C 全部MAG按5厘米间隔的BGC密度箱线图；D 酸杆菌MAG同口径BGC密度。箱线图显示中位数、四分位与1.5倍四分位距须线；深度组间检验采用Kruskal-Wallis，事后Dunn-BH校正。

结果3 高BGC储量MAG与分类贡献谱系

在全部MAG中，52个基因组含至少5个BGC；其中假单胞菌门两个α-变形菌MAG各含24个BGC（含5个PKS、5个NRPS、5个RiPP），另有13个假单胞菌MAG含10–15个BGC，作者定义为“高产区（talented producer）”。蛭弧菌门（Bdellovibrionota）两个UBA1018 MAG分别含17、18个BGC，其中RiPP占10、12个；酸杆菌14个MAG含10–15个BGC，多归Terriglobia与Bryobacteriaceae。古菌中Halobacteriota一株含4个BGC（3个NRPS），Thermoproteota与Thermoplasmatota各3个，并检获氢氰酸与darobactin样RiPP潜力。该结果说明高BGC储量不限于经典放线菌/假单胞菌，而广泛分布于酸杆菌、蛭弧菌、古菌等类群。

图3 高BGC储量MAG系统发育与BGC组成。按分类着色系统树展示入选MAG；内圈条形为基因组大小（Mbp），外圈条形为BiG-SCAPE各类BGC数量（RiPP、NRPS、PKS、萜类等）。该图综合呈现高储量基因组的系统分布与产物类型偏好。

结果4 群落水平BGC转录与酸杆菌主导

contig水平以至少1个生物合成核心基因TPM大于0.25为表达判定。全体1,120个BGC（13.7%）有表达；每样以萜类表达数最多（均值129）、其次RiPP（103）、PKS（25）、NRPS（29）。酸杆菌每样平均表达153个BGC，占全部表达BGC的68.2%；经按宏基因组丰度归一后，其BGC转录份额log2富集中位约0.77–0.80，即相对其群落占比约1.7–1.75倍。假单胞菌原始贡献次之但归一后低于丰度预期，脱硫菌_B（Desulfobacterota_B）呈正向但绝对量较低。结果表明酸杆菌主导并非仅由高丰度造成，而是单位丰度转录活性更高。

图4 群落BGC转录贡献。A 各深度BGC核心基因总TPM按分类堆叠；B 各分类表达BGC数量均值；C 按表达数量份额减宏基因组丰度份额的log2富集；D 按核心基因TPM份额的同口径富集。箱线图展示12个样品三深度剖面，仅纳入全部样本均存在的门，未做推断统计。

结果5 基因组水平激活率、储量负相关与酸杆菌/假单胞菌性状分化

MAG水平仅统计活跃基因组（7个单拷贝看家基因均值TPM大于0.5），以核心基因达到家基因阈值、超过家基因均值或TPM大于0.8判定表达。全体280个BGC在MAG中表达，每样41–92个，占活跃群落潜力8.8%–27%。按产物分，ranthipeptide表达比例最高（均33%），其后RiPP-like 24.4%、I型PKS 20.1%、lassopeptide 18.9%、III型PKS 17.3%、NRPS 11.9%。酸杆菌MAG表达163个BGC，每样占其储量8.5%–31%；假单胞菌仅表达40个，每样2.4%–19%。两个各含24个BGC的假单胞菌MAG及两个18个BGC的蛭弧菌MAG均无核心基因转录。全部MAG中BGC储量与表达比例呈显著负相关（Spearman ρ=-0.522，P=1.26×10^-13）。性状上，假单胞菌预测倍增时间更短（均5.77 vs酸杆菌8.54小时，Welch t=7.35，P=1.45×10^-11），但GH/AA/PL类CAZyme注释数更低（均69.9 vs 152.6，t=7.84，P=6.22×10^-12）；酸杆菌以更小储量、更强复杂碳降解酶表达实现更高BGC激活。

图5 基因组水平BGC表达。A 各样品表达BGC数量按产物类型分布；B 各门表达BGC绝对数；C 各门表达占比箱线；D 含至少20个编码BGC的门按“表达BGC数/编码BGC数”计算激活率。图中未做深度间推断检验，仅汇总12个生物学样本变化。

结果6 MAG水平代谢背景：GH表达最强耦合，氧化还原背景异质

作者选取15个高活跃焦点MAG（跨越酸杆菌、假单胞菌、绿弯菌、脱硫菌、疣微菌、放线菌），以45个MAG×深度观测建立线性混合模型（MAG与深度为随机效应）。BGC核心基因TPM与GH表达正相关最强（β=0.528，t=12.37，BH-P=8.98×10^{-22）；资源获取、厌氧能量代谢、好氧呼吸、胁迫保护亦正相关的但强度较低。扩展至127个活跃MAG、566观测，GH标准化TPM对BGC核心TPM斜率β=0.447（t=20.44，P=1.93×10}-69），加入看家基因归一后仍稳定（β=0.525）。逻辑模型显示BGC表达观测中97%伴GH表达、99.4%伴CAZyme表达；GH是表达比例的正向预测因子（z=7.38，OR=57.9）。氧化还原方面，焦点MAG同时出现微氧终端氧化酶、好氧氧化酶、厌氧标记、硫/铁/氮/发酵不同组合，无单一共有 redox 状态。因此BGC原位表达最稳定地耦合碳加工特别是GH，而不是某一种胁迫或氧态。

图6 MAG水平BGC表达与代谢背景。A 焦点MAG各功能模块标记基因汇总表达热图，红环为BGC核心基因均值；B 同批MAG氧化还原相关模块热图，红环同口径；C GH标准化表达对BGC核心基因标准化表达散点，实心为重点MAG、空心为其他活跃MAG，黑线为混合模型固定效应拟合。所有TPM均以单拷贝看家基因均值归一，功能类别低于0.5 TPM留空。

方法学参考

上述流程中，antiSMASH加核心基因位置过滤可降低短读长对NRPS/PKS边界的碎片高估；BiG-SCAPE聚类适合大样本去冗余；以98%ANI去重MAG后再做转录映射可减少多映射；看家基因TPM归一的混合模型适用于深度重复但非技术重复的生态转录数据。

总结

研究仅用7个宏基因组即回收近9,700个BGC，规模可与更大尺度海洋、土壤、反刍动物调查相比，但97.1%在BGC-Atlas、99.9%在MIBiG中缺乏近缘簇，证实泥炭地为次级代谢新颖性热点。原位转录仅激活小部分潜力：群落水平酸杆菌以中等储量贡献过半转录，基因组水平高储量假单胞菌与蛭弧菌多沉默，BGC比例随每基因组簇数增加而下降。机制上，酸杆菌具有更长预测倍增时间但更多GH/AA/PL类CAZyme，其BGC激活与复杂碳降解耦合；假单胞菌倍增快、营养摄取与群体感应调节丰富，但高储能BGC曲目在泥炭地营养与氧化还原条件下多不部署。MAG分辨率显示BGC表达可发生于好氧、微氧与厌氧背景，最一致信号为GH与总CAZyme而非单一redox或胁迫状态。该框架说明泥炭地BGC调查需同时结合基因组潜力、转录激活与碳加工表型，单靠BGC目录会高估生态功能。

参考文献

• Human Z R, Hoover R, Chan C, Küsel K, Wegner C E（2026）. High diversity but limited expression of biosynthetic gene clusters in a peatland microbial community. Nature Communications 17, 10529. https://doi.org/10.1038/s41467-026-78121-8
• NCBI SRA BioProject PRJNA1374247：https://www.ncbi.nlm.nih.gov/sra
• Zenodo MAG与BGC数据集：https://doi.org/10.5281/zenodo.17812908
• antiSMASH 7.0文档与BGC预测：https://antismash.secondarymetabolites.org/
• BiG-SCAPE与BiG-SLiCE聚类：https://github.com/NBChub/biG-SCAPE ；BGC-Atlas https://bgcat.rug.nl/
• MIBiG 4.0实验验证BGC库：https://mibig.secondarymetabolites.org/
• dbCAN2 CAZyme注释：http://bcb.unl.edu/dbCAN2/ ；KEGG KOfam：https://www.kegg.jp/
• GTDB与GTDB-Tk分类：https://gtdb.ecogenomic.org/ ；CheckM2分箱质量：https://github.com/chklovski/CheckM2
• Nature Communications期刊指标：ISSN 2041-1723，2025 JCR IF 18.1、五年18.9、Q1多学科（2026 JCR发布）