功能宏基因组学驱动环境样本活性纤维素酶高效挖掘：基于多重PCR与滚环扩增的毕赤酵母（K. phaffii）快速筛选体系

《Applied Microbiology and Biotechnology》：A fast workflow to explore active enzymes from environmental samples through functional metagenomics

针对传统宏基因组酶挖掘依赖大肠杆菌中间转化、流程繁琐且难以平衡效率与新颖性的痛点，本研究以印尼濒危鹿Axis kuhlii粪便为样本，开发“宏基因组测序-多重PCR构建基因库-RCA（滚环扩增）简化克隆-半高通量筛选” workflow，在Komagatella phaffii中实现功能表达，最终获得5种内切葡聚糖酶（endoglucanase）和3种β-葡萄糖苷酶（β-glucosidase），为环境微生物组功能酶发现提供了“计算预测-实验验证”的高效桥梁。

在印度尼西亚巴韦安岛的热带雨林里，生活着一种全球仅存约300只的极度濒危鹿——Axis kuhlii（巴韦安鹿）。这种小鹿的粪便里，藏着微生物世界的“宝藏”：能分解木质纤维素的酶。但长期以来，科学家想从环境样本（比如粪便、土壤）中挖出这些“宝藏酶”，却总被两道坎卡住：一是传统培养法只能搞定1%的可培养微生物，剩下99%的“不可培养微生物”基因像锁在黑箱里；二是就算用宏基因组学测出了基因序列，要把这些基因变成能实际用的活性酶，还得经过大肠杆菌转化、质粒扩增等一堆步骤，又慢又贵，还容易漏掉那些“长得不像已知酶”的新颖候选者。

这时候，功能宏基因组学站了出来——它不只看基因序列，更直接找“能干活”的酶。但怎么让这个过程更快、更省、更能挖到“新货”？来自国外的研究团队把目光投向了Axis kuhlii的粪便：这种鹿吃植物纤维，肠道微生物肯定有厉害的纤维素酶；而且它是濒危物种，样本里的微生物可能带着独特的酶特性。他们琢磨出一个“组合拳”：用宏基因组测序找候选基因，用多重PCR（multiplex PCR）一次性扩多个目标，用滚环扩增（RCA）跳过大肠杆菌步骤直接给毕赤酵母（Komagatella phaffii）准备DNA，再用96孔板半高通量筛选——这套流程能不能打破传统瓶颈？

研究结果没让人失望：他们真的从粪便里挖到了5种内切葡聚糖酶和3种β-葡萄糖苷酶，而且整个过程比传统方法快了不止一倍。这篇论文后来发表在《Applied Microbiology and Biotechnology》上，给环境酶挖掘指了条“高效路”。

关键技术方法

研究以印尼巴韦安岛野生动物保护区的9份Axis kuhlii粪便（3份自由放养、2份圈养群体）为样本，先通过纳米孔宏基因组测序组装注释纤维素酶基因；再用多重PCR扩增目标基因（内切葡聚糖酶选退火温度62-64°C的16个，β-葡萄糖苷酶选64-66°C的10个）；接着用RCA（滚环扩增）直接扩增重组质粒，跳过大肠杆菌转化；最后用96孔板深孔板培养毕赤酵母，通过DNS法（羧甲基纤维素底物）、4-MUC（纤维二糖水解酶活性）、4-MUG（β-葡萄糖苷酶活性）进行半高通量筛选，结合SDS-PAGE、Western blot验证蛋白表达，BLAST和MEGA12.1做序列分析与系统发育树。

研究结果

Amplification of cellulase genes by multiplex PCR and RCA

通过宏基因组测序从9份粪便中鉴定出31个内切葡聚糖酶和154个β-葡萄糖苷酶基因，选退火温度适配的多重PCR引物（内切16个、β-葡萄糖苷酶10个）扩增，产物经HiFi组装入pPICZαA载体后，用RCA扩增得到线性化DNA（内切~4400bp、β-葡萄糖苷酶~5600bp，对应载体+基因），电泳验证了片段大小正确性。

Cultivation and screening for active enzymes in K. phaffii using 96-well plate format

转化毕赤酵母后得到超150个转化子/组，筛选到32个内切葡聚糖酶阳性菌落、24个β-葡萄糖苷酶阳性菌落；Sanger测序确认31个内切菌落含目标基因、11个β-葡萄糖苷酶菌落含目标基因，且转化子数量与宏基因组模板基因拷贝数正相关（如内切EG Q4有35拷贝→12菌落，β-葡萄糖苷酶BGL1有196拷贝→8菌落）。

SDS PAGE of five endoglucanases and three β-glucosidases

内切葡聚糖酶的SDS-PAGE显示~34kDa主导条带（EG11偶有~24kDa条带），表达良好；β-葡萄糖苷酶表达较弱，SDS-PAGE条带不清，但Western blot（抗His标签抗体）确认~83kDa目标蛋白存在。

Rapid screening of colonies for other properties using a 96-well plate

对32个内切转化子用4-MUC筛外切酶活性，发现13个含EG4/13/25/28的克隆有活性；1个内切克隆能水解木聚糖（具混杂木聚糖酶活性）；还筛到在酸性（pH5.0）、碱性（pH8.0）及60°C下仍保持高活性的克隆。

BLAST sequence comparisons

BLAST结果显示：内切葡聚糖酶与UniProtKB/SwissProt相似性24%-45%、与PDB相似性41%-51%；β-葡萄糖苷酶与UniProtKB/SwissProt相似性32%-51%、与PDB相似性33%-52%。系统发育树显示这些酶与已知序列差异大，具新颖性。

研究结论与讨论

研究用“宏基因组测序-多重PCR-RCA-半高通量筛选”workflow，从Axis kuhlii粪便中快速挖到5种内切葡聚糖酶和3种β-葡萄糖苷酶，全程仅需基因注释后2周，比传统“逐个克隆+大肠杆菌扩增”的方法省了大量时间和 labor。

为什么这套方法厉害？首先是快：RCA跳过了大肠杆菌转化步骤，直接给毕赤酵母准备DNA，不用再养大肠杆菌、提质粒；其次是省：多重PCR一次扩多个基因，比定制合成每个基因便宜多了；最重要的是能挖新货：多重PCR可以同时试高相似（低风险）和低相似（高风险）的候选基因，不会只盯着已知酶——比如研究中得到的酶与UniProt/PDB相似性大多低于52%，都是“新面孔”。

再看宿主选择：Komagatella phaffii是真核表达系统，能做翻译后修饰（比如糖基化），分泌能力强，本来适合真核基因，但这研究里连原核来源的纤维素酶都能表达成功——说明它也能搞定原核宏基因组基因，打破了“真核宿主只表达真核基因”的刻板印象。

当然，方法也有局限：比如β-葡萄糖苷酶的10个扩增基因里，只有3个产生了活性克隆，原因可能是扩增偏好（CG含量高的基因难扩）、密码子不适应（没做优化导致表达量低），或者蛋白折叠失败——这也给后续优化指了方向：比如调整多重PCR引物浓度比，或者给基因做密码子优化。

但整体来看，这套方法的潜力很大：它能处理不可培养微生物的基因，用96孔板同时筛多个属性（不同底物、pH、温度），要是遇到符合要求的酶，还能快速转到工艺优化和放大——比如文中的酶已经能在酸性、高温下工作，说不定能用在生物燃料或洗涤剂里。

说到底，这项研究的意义不只是找到8个纤维素酶，更是给环境微生物组的“功能基因挖掘”提供了一套“高效模板”：不用再被大肠杆菌绑架，不用再一个个试基因，而是用多重PCR+RCA把“选基因-扩基因-表达”串成流水线，让藏在环境样本里的“酶宝藏”更快见光。对于濒危物种的样本来说，这更是一举两得——既保护了物种，又利用了它们的微生物资源，说不定未来还能靠这些酶开发出新的工业应用，让“保护”和“利用”不再矛盾。