DNA中的潜力与真正转录的RNA
关键问题不是存在多少BGC，而是有多少被使用。BGC是一组在染色体上成簇排列的基因，它们共同合成一种次级代谢产物；在DNA中找到BGC只说明具备这种能力，并不说明它正在发挥作用。为区分两者，作者从宏基因组（群落DNA）中鉴定BGC，并用取自同一泥炭层的宏转录组（群落RNA）测量BGC转录。
研究地点是巴伐利亚北部的Schlöppnerbrunnen矿质泥炭沼泽，这是一处弱酸性泥炭地（上层30 cm内pH约4.7），由缺氧、富铁的地下水补给。2020年7月采集的两根平行土芯按5 cm间隔分样。分析区间为5至30 cm，覆盖从含氧表层到缺氧、饱水层的过渡。由于野外条件非常干燥，顶部0–5 cm未纳入；20–30 cm层虽多次尝试仍无法回收RNA，因此最深层只提供DNA数据，没有表达数据。
近万个基因簇，多数没有近缘参考
作者用antiSMASH在长度超过5 kb的组装contig上鉴定BGC，共得到9,694个簇，其中1,594个（16.4%）为完整簇，平均长度31,073 bp。由于短读长组装可能把大型模块化簇拆开、从而虚增计数，作者把生物合成核心基因全部落在contig末端5 kb以内的不完整NRPS、PKS及NRPS-PKS杂合BGC标记为可疑片段。这一保守过滤把9,694个BGC中的1,544个（15.9%）列为可能高估，其中包括41.5%的NRPS簇和39.2%的PKS簇。
片段化对多样性的影响大于对原始计数的影响。剔除这些弱片段后，BGC总数减少16%，但基因簇家族（GCF，即BiG-SCAPE对相关簇的分组）减少26%，从4,669降至3,452。即便经过这一过滤，七个宏基因组仍产生3,452个GCF；作者将其规模与采样量大得多的全球调查相比较，同时指出不同研究的采样设计和阈值并不严格等同。多数簇在参考数据库中缺乏近缘序列，论文将此解读为广泛的新颖性。
产物类别以核糖体合成并经翻译后修饰的肽（RiPP，n = 1,361）为主，其次是萜烯（620）、非核糖体肽合成酶（NRPS，510）、聚酮合酶（PKS，367）和NRPS-PKS杂合簇（107）。GCF在分类上分布广泛，以Acidobacteriota（1,404）和Pseudomonadota（754）最多，另有40个古菌GCF。
仅9–27%的编码潜力被转录
核心定量结果是DNA与RNA之间的落差：所回收BGC潜力中只有9–27%在原位表达。这是本数据集内已转录簇与已编码簇的比值，并非对泥炭地整体BGC的断言，且仅基于一处沼泽的七个宏基因组。作者把这一范围视为环境与生理约束、而非基因组容量主导次级代谢的证据。
正是这一落差让丰度与活性的比较具有意义。统计簇的数量只能说明群落可能合成哪些化学产物；统计转录本则说明采样时哪些簇正在被使用。论文自身的表述是：仅凭基因组容量很难预测实际发生的化学过程。
Acidobacteriota编码更少却转录更多
Acidobacteriota的基因簇储备量只是中等，却贡献了全部BGC转录的一半以上。与之对照的是BGC丰富的Pseudomonadota：论文称这类类群为“高产者”，每个基因组携带多达24个簇，但大多保持沉默。在基因组层面，BGC储备中被表达的比例与每个基因组的BGC数量呈负相关，用Spearman秩相关量化。
作者把这种解耦归因于生活史策略。Pseudomonadota的预测倍增时间更短（由gRodon2基于全基因组密码子使用偏性估计），复杂碳加工潜力也更低。Acidobacteriota则专门适应复杂碳加工，并表达了其储备中更大的比例。作者给出的解释是：快速生长类群的大储备代表潜在能力，只有在调控与能量条件允许时才会启用，而作者认为这类条件在泥炭地中并不常见或持续时间很短。这是作者对相关关系的解释，而非已证实的因果。
表达与碳水化合物活性酶活性同步
在基因组层面，BGC表达程度与碳水化合物活性酶（CAZyme）表达关系最紧密，尤其是糖苷水解酶（GH），即负责分解复杂碳水化合物的酶。在所有活跃MAG中，97%的BGC表达观测同时出现GH表达，99.4%出现CAZyme表达。
在控制基因组与深度的模型中，这一关联依然成立。GH表达与BGC表达在混合逻辑回归模型中相关（优势比57.9，95% CI 19.7–170）。在控制MAG总转录、深度和MAG身份后，GH表达越高，表达的BGC储备比例也越大（优势比4.78，95% CI 2.73–8.35）。以MAG身份和深度为随机截距的线性混合效应模型给出，归一化GH表达对BGC核心基因表达的群体水平斜率为β = 0.447（95% CI 0.404–0.489）。CAZyme总表达几乎无处不在，因此作为二分类预测变量的区分度较低。
相比之下，氧化还原背景更为异质。在BGC活跃的重点MAG观测中，表达与微氧或兼性末端氧化酶、好氧氧化酶、厌氧功能标志物以及低表达或未归类的氧化还原模块同时出现。15个重点MAG跨越Acidobacteriota、Pseudomonadota、Chloroflexota、Desulfobacterota、Desulfobacterota_B、Verrucomicrobiota和Actinomycetota，它们之间没有共享单一的氧化还原特征。
生物合成能力随深度下降，表达并未停止
来自更深泥炭层的MAG编码的BGC更少，但那里的基因组也更小，因此作者按基因组大小做了归一化。BGC密度（每Mbp的BGC数）仍随深度下降（Kruskal–Wallis H(4) = 16.2，P = 0.00275，η²H = 0.0216）。在Acidobacteriota内部也出现相同方向（H(4) = 12.8，P = 0.0125，η²H = 0.0695）。被预测为厌氧的基因组每个基因组携带的BGC显著少于预测为好氧的基因组（Welch t检验，t(290.17) = 7.62，P = 3.56×10⁻¹³，Hedges' g = 0.68）。除深层MAG完整度略高、GC含量略高外，各深度的MAG质量指标没有差异。
然而，在所有采样深度都检测到表达，包括更深的缺氧层，且伴随多种共表达的氧化还原背景。作者将此解读为：即使能力下降、能量受限，次级代谢仍然持续；并认为这证明Acidobacteriota能在含氧表层和更深缺氧区都表达BGC。
这处沼泽能支持什么、不能定论什么
证据基础是2020年7月采自一处沼泽两根土芯的七个宏基因组，最深20–30 cm区间缺少表达数据，表层0–5 cm被排除。9–27%的表达范围和Acidobacteriota的主导地位属于该地点与该采样窗口，并非对整个泥炭地的估计。作者也指出，由于采样设计和分析阈值不同，跨研究的GCF丰富度比较并不严格等同。
机制层面的论断是相关性。每个基因组BGC数与表达比例之间的负相关，以及BGC表达与GH表达的耦合，都是跨MAG和深度量化的关联。快速生长的Pseudomonadota持有潜在储备、等待罕见许可条件的说法，是与数据相符的解释，而非经过检验的因果路径。
论文确立的实际优先级是一个比较，而非一份目录：在这处沼泽内，转录BGC最多的类群并不是编码BGC最多的类群。Acidobacteriota储备中等、复杂碳加工活跃，却贡献了一半以上的BGC转录；而簇丰富的生产者大多沉默。任何仅按簇数量排序生物合成重要性的调查，在这里都会把顺序颠倒——这正是作者主张多样性编目不能替代原位表达测量的具体理由。
本文引用的文献
High diversity but limited expression of biosynthetic gene clusters in a peatland microbial community
作者从七个泥炭地宏基因组中回收了9,694个BGC，但仅9–27%在原位被转录；Acidobacteriota储备中等却贡献了一半以上的BGC转录，且BGC表达与糖苷水解酶表达关系最紧密。