4. 核酸提取工程
宏基因组测序对核酸提取的要求与常规单菌基因组测序有本质区别：样本中微生物种类繁多、丰度跨度可达数个数量级，同时混杂大量宿主细胞、腐殖酸、多糖、多酚等干扰物。提取环节一旦出现偏差，后续建库和生信分析都无法补救。本节围绕方法选择、裂解优化、抑制剂去除和质量评估四个环节展开。
一、DNA提取方法的选择依据
宏基因组DNA提取方法大致分为三类：机械裂解法（珠磨法）、化学裂解法（SDS、CTAB、酶解）和商品化柱式/磁珠试剂盒法。选择时需权衡三个维度：细胞壁裂解强度、片段完整性与得率、抑制剂残留水平。
珠磨法通过高速震荡使玻璃珠或锆珠对细胞产生剪切与碰撞，对革兰氏阳性菌、芽孢、真菌孢子等难裂解细胞效果显著，但剪切力会打断DNA，片段长度通常集中在10–30 kb，对需要长片段（如Nanopore、PacBio）的项目不利。化学裂解法温和，片段可达50 kb以上，但对厚壁菌和孢子裂解不彻底，容易造成革兰氏阳性菌丰度被系统性低估。商品化试剂盒通常将酶解、化学裂解与柱纯化整合，重复性好，但不同品牌对同一土壤样本的群落结构还原度差异可达20%以上。
实际选择建议：土壤、沉积物、粪便等复杂基质优先考虑珠磨+柱纯化组合；水体滤膜、培养物可用温和化学裂解；长读长项目应避免高强度珠磨，改用酶解+低剪切纯化。若研究目标是革兰氏阳性菌或古菌，必须验证裂解覆盖度，可通过加入已知量的难裂解菌（如Bacillus subtilis孢子）作为内标评估。
二、裂解效率优化
裂解效率直接决定群落结构的真实性。优化目标不是“裂解得越狠越好”，而是在覆盖度与片段完整性之间找到平衡点。
珠磨法的关键参数包括：珠材质与粒径、震荡速度、震荡时间、循环次数。以土壤样本为例，常用方案为0.1 mm锆珠+0.5 mm玻璃珠混合，6.0 m/s震荡40 s，间歇1 min，重复2–3次。震荡时间过短，革兰氏阳性菌裂解不足；过长则DNA片段严重碎片化，且产生过多气溶胶增加交叉污染风险。建议首次使用某方案时，设置时间梯度（20 s、40 s、60 s、90 s）并配合内标回收率评估。
化学裂解优化主要围绕裂解液配方：SDS浓度（通常0.5%–2%）、蛋白酶K用量（20–100 μg/mL）、CTAB浓度（针对多糖多酚样本）、以及孵育温度与时间（55–65°C，30–120 min）。对于腐殖酸含量高的样本，可在裂解液中加入PVPP（聚乙烯聚吡咯烷酮）或CaCl₂，通过竞争结合减少腐殖酸与DNA的共沉淀。
酶解预处理对真菌和孢子尤为重要。溶壁酶（lyticase）、几丁质酶、溶菌酶的组合可在机械裂解前先破坏细胞壁，降低后续珠磨强度，从而保护DNA完整性。典型流程为：样本重悬于酶解缓冲液 → 加入溶菌酶（1 mg/mL）和溶壁酶（0.5 U/μL）→ 37°C孵育30–60 min → 再进行化学裂解或温和珠磨。
三、抑制剂去除
宏基因组样本中最常见的抑制剂包括腐殖酸、富里酸、多糖、多酚、重金属离子和残留乙醇。这些物质会抑制PCR聚合酶活性、影响文库接头连接效率，甚至导致测序信号异常。
去除策略需根据抑制剂类型选择：
- 腐殖酸/富里酸：柱纯化时使用含PVPP的裂解液；或在纯化前用CaCl₂沉淀腐殖酸；也可采用琼脂糖凝胶电泳回收DNA，但得率损失较大。
- 多糖：CTAB裂解液可有效沉淀多糖；纯化时使用高盐缓冲液促进多糖与DNA分离；必要时用PEG沉淀替代乙醇沉淀。
- 多酚：裂解液中加入β-巯基乙醇或DTT防止氧化；PVPP可吸附多酚。
- 残留乙醇/盐：柱纯化后以70%乙醇洗涤两次，离心干燥至无乙醇气味，但避免过度干燥导致DNA难溶。
抑制剂去除效果不能仅凭A260/A280判断，必须通过PCR扩增测试验证。常用方法：取1 μL DNA模板进行16S rRNA基因V3–V4区PCR（25–30循环），同时设置阳性对照（已知纯DNA）和阴性对照（无模板）。若样本PCR无产物或条带极弱，而阳性对照正常，则提示抑制剂残留。此时可将DNA稀释5–10倍后重新PCR，若扩增恢复，说明抑制为可稀释型，通常为腐殖酸或盐类残留。
四、提取质量评估
提取质量评估应覆盖四个维度：浓度、纯度、完整性和功能可用性。
| 评估维度 | 检测方法 | 合格参考范围（示例） | 异常提示 |
|---|---|---|---|
| 浓度 | Qubit荧光定量 | ≥1 ng/μL（土壤）；≥0.5 ng/μL（水体） | 过低：裂解不足或纯化损失 |
| 纯度 | NanoDrop A260/A280、A260/A230 | A260/A280 1.8–2.0；A260/A230 ≥1.8 | A260/A230偏低：腐殖酸/多糖残留 |
| 完整性 | 琼脂糖凝胶电泳或片段分析仪 | 主带≥10 kb，无明显弥散 | 严重弥散：过度剪切或核酸酶污染 |
| 功能可用性 | 16S PCR扩增 | 25–30循环可见清晰条带 | 无条带：抑制剂残留或模板不足 |
需要强调的是，NanoDrop测得的浓度常因腐殖酸和RNA干扰而偏高，应以Qubit荧光定量为准。A260/A230比值低于1.5时，通常提示腐殖酸或碳水化合物残留，建议重新纯化或稀释后使用。完整性评估中，若主带模糊且低分子量弥散严重，说明DNA已发生明显降解，对于需要长片段建库的项目应重新提取。
功能可用性测试是最终判据。即使浓度和纯度指标合格，若PCR扩增失败，仍不能用于建库。建议每批提取设置一个已知群落结构的标准品（如模拟微生物群落），通过16S测序或qPCR评估提取效率的批次间重复性。对于关键项目，可进一步用qPCR定量特定菌属（如Escherichia coli、Bacillus subtilis）的回收率，判断是否存在选择性裂解偏差。
综上，宏基因组DNA提取工程的核心在于：根据样本类型和下游需求选择合适方法，通过参数优化平衡裂解覆盖度与片段完整性，针对性地去除抑制剂，并以功能可用性为最终判据进行质量评估。只有这四个环节协同控制，才能获得真实反映微生物群落结构的合格DNA模板。