来源：CRISPR干扰菌群丰度变化定量校正（DNA Lab Space 实验教程）| 编译：DNA Lab Space
菌群实验里，你用CRISPRi把目标菌株的某个基因干扰掉了，预期它的丰度该掉下来，结果Biological replicate三组数据横着飞，要么趋势翻转，要么标准偏差大到连你自己都心虚。这不是CRISPRi没用，大概率是你的定量校正路径出了岔子。
这篇教程直接聚焦一个场景：混合菌群中，CRISPRi干扰株的丰度变化到底怎么测才可信。覆盖从菌株活化、诱导条件、DNA提取到qPCR定量的排查流程，加上两个进阶疑难场景。别急着重新设计sgRNA，先从数据源头走一遍。
第一步：快速排查低级错误
先别打开分析软件。这几个坑，新手十个有九个踩过，检查只要五分钟。
- 诱导剂忘了加或浓度加错：aTc工作液终浓度应为100 ng/mL。你如果用的是储液1 mg/mL，稀释1000倍；如果加成了10 ng/mL，相当于没诱导，沉默效率骤降，丰度当然不变。直接重新补加aTc，37°C再孵育6小时后取样。
- 取样时间点不一致：OD600到0.4-0.6再加aTc，诱导6小时后收菌。有人3小时收，有人8小时收，丰度差值能到2倍以上。把Biological replicate的取样时间精确到±15分钟内。
- 菌群样本没破壁完全：革兰氏阳性菌占比高时（如产丁酸菌），只用试剂盒裂解液处理等于白做。加溶菌酶（20 mg/mL，37°C孵育30分钟）再走试剂盒流程，DNA得率立刻不一样。
- 质粒拷贝数混入定量口径：CRISPRi系统如果以质粒形式存在，你qPCR定的是基因组DNA还是总DNA？质粒在菌体内拷贝数不稳定，会直接污染dCas9/sgRNA序列的定量信号。务必用DNase消化去除游离DNA后再提取基因组，或改用针对染色体单拷贝基因的引物做校正基准。
第二步：分维度系统排查
低级错误排除了，数据还是有你解释不了的趋势别慌。按实验流程拆分，逐个模块筛查。
1. 样本制备模块：DNA得率不均
现象：三组平行样本，同一个内参基因的Cq值差了3个循环以上。
原理：裂解不完全或DNA纯化回收率不同，会直接放大丰度差异。菌群中不同菌的细胞壁结构差异很大，统一用同一条件裂解必然偏颇。
解决步骤：
- 每个样本取1 mL菌液，OD600统一调节到1.0。
- 加入20 mg/mL溶菌酶，37°C水浴30分钟，每10分钟颠倒一次。
- 后续提取按试剂盒说明书执行，但洗脱体积固定为50 μL，提前预热洗脱缓冲液至60°C。
- 使用Nanodrop检测A260/A280和浓度，DNA浓度统一稀释到10 ng/μL再进行qPCR。
验证方法：取同一样本，独立提取两次DNA。两次得到的基因组DNA各做管家基因（如16S rRNA通用引物）qPCR，看Cq差值是否≤0.5。如果大于0.5，说明提取流程本身波动太大，不能用于后续定量。
2. 扩增与诱导模块：基因表达时序不清
现象：诱导前后靶标菌丰度没有显著变化，或者抑制效果特别弱。
原理：CRISPRi的转录抑制效果取决于sgRNA的表达量以及dCas9-sgRNA复合物与靶DNA结合的有效浓度。若在菌群中诱导时间不够长，已经表达的靶蛋白还有很长半衰期，代谢表型不会立刻转变成生长速率的可测差异。
解决步骤：
- 诱导条件先用纯培养验证：将携带CRISPRi质粒的靶标菌单菌落接种至5 mL LB培养基（含相应抗生素），37°C振摇过夜。
- 次日按1:100转接至10 mL新鲜培养基，37°C振摇至OD600为0.4-0.6。
- 加入aTc使其终浓度为100 ng/mL，37°C继续培养6小时。
- 以未诱导培养物作为对照，提取RNA做逆转录定量PCR（RT-qPCR），确认靶基因的mRNA水平降至对照的30%以下。达不到这个阈值，不要进入混合菌群实验。
验证方法：做一个时间梯度——诱导后0小时、2小时、4小时、6小时、8小时分别取样，算靶标菌丰度占比。你会看到，多数CRISPRi系统在4小时后才出现可测的生长抑制。如果到8小时都没变化，换sgRNA结合位点设计，不纠结。
3. 定量方法模块：qPCR引物专一性差
现象：溶解曲线出现两个峰或主峰前有明显杂峰，扩增效率超过105%。
原理：菌群样本背景复杂，非靶标菌的基因组DNA可能与引物部分配对，带来非特异扩增，从而干扰靶标菌基因组DNA的绝对定量折算。
解决步骤：
- 在NCBI Primer BLAST里，把你的引物序列与样品中可能存在的菌种基因组序列比对，剔除3'端有2个及以上碱基与非靶菌种完全匹配的候选引物。
- 梯度稀释靶标菌基因组DNA，做标准曲线。斜率应在-3.6到-3.1之间，对应扩增效率90%-110%。如果斜率偏平或偏陡，调整退火温度（60°C附近做2°C的梯度），或增加镁离子浓度至2.5 mM（现成PCR预混液除外）。
- 用1.5%琼脂糖凝胶电泳跑qPCR产物验证产物大小，单一条带才算通过。
验证方法：把靶标菌DNA与等量的另一种非靶菌DNA混合，做加标回收实验。正常情况下，靶标菌qPCR拷贝数应回收80%-120%。如果回收率低于80%，你的引物大概率被杂菌干扰了。
4. 数据处理模块：内参基因缺乏独立校准
现象：靶标菌的绝对拷贝数波动大，但相对占比稳定，然而你报告的是绝对拷贝数变化。
原理：直接比较各组靶标菌基因拷贝数，有一个前提——各样本DNA提取总量一致。而菌群实验往往做不到，因为整体菌量可能因为CRISPRi的影响而发生了变化。
解决步骤：
- 引入外源内标：在取样后提取DNA前，向每管样本中加入等量的异源微生物（如乳酸乳球菌纯基因组DNA，或质粒载体作为spike-in），建议加到每管样本体积的1/100（如1 µL spike-in至100 µL样本），浓度不超过10 ng/mL，保证在qPCR线性范围内。
- 以spike-in对应基因为参考，用ΔΔCq法计算靶标菌DNA相对含量变化。计算公式：ΔΔCq =（靶标菌Cq⁻靶标菌内参基因spike-in Cq）处理组 −（靶标菌Cq⁻spike-in Cq）对照组，结果以2^（−ΔΔCq）表示。
- 同时测定样本中总菌的16S rRNA基因拷贝数，将靶标菌DNA相对变化换算为占总量比例变化。
验证方法：把每组样本各取2份，其中1份加spike-in后立即提取DNA，另1份不加，平行跑qPCR对比。加spike-in的组，Cq在板间和三组间的波动应≤0.5个循环。
第三步：进阶疑难问题
基础系统排查都能通过，但数据结果依然不符合预期？考虑两个深水区问题。
sgRNA脱靶导致非靶标菌群体显着下降
现象：CRISPRi干扰组中，非靶菌种（理论上不受影响的背景菌）丰度也同步下降。你甚至能看到整个菌群的Shannon指数在下降。
原理：菌群是网络结构，你干扰了靶标菌的某个必需基因（比如负责合成某种氨基酸或短链脂肪酸），靶标菌自身的丰度掉了，它分泌的代谢物也随之消失，依赖它的其他菌群成员也跟着受到波及。这个叫生态连锁效应。
解决步骤：
- 设置一个“空载对照”——携带不含靶向序列的sgRNA支架的菌株，走完全相同的实验流程。这是目前最直接的控制。
- 把受影响的非靶标菌单独挑出来，在纯培养条件下与CRISPRi菌株共培养，排除物理接触或代谢物依赖的间接效应。
- 如果你确认sgRNA序列与这些非靶菌的基因组有16 bp以上连续同一序列，那就要重新设计。使用BLAST（短链参数：word size 7，E-value 100）筛掉潜在匹配靶点。
验证方法：全基因组测序验证dCas9结合位点，或做全转录组测序看这些非靶菌基因的下调幅度是否只集中在某个代谢通路模块——后者提示生态级联，前者提示脱靶，处理逻辑完全不同。
质粒在混菌传代中丢失，你测的丰度其实是没带系统的背景菌生长
现象：qPCR数据显示靶标菌丰度回来了，第4天定量高于第1天，阳性对照却正常。
原理：混菌培养存在传代压力。没有抗生素维持的体系，CRISPRi质粒常常在8-12代内丢失，不带质粒的“逃逸菌”迅速抢占生态位。CRISPRi系统会拖慢携带者生长速度，如果没有筛选压力，丢失质粒的细胞就是最大的赢家。
解决步骤：
- 取样前转移培养物前，每2天补充一次抗生素。比如携带氯霉素抗性，浓度保持在5 µg/mL，不超过其最低抑菌浓度的一半，避免彻底杀掉非靶标背景菌。
- 每次传代时，涂板检查质粒保留率。取100 µL适当稀释的培养物，涂布于含抗生素和不含抗生素的LB平板。在含抗生素平板上长出来的菌落数 ÷ 不含抗生素平板的菌落数即为保留率，间期应≥80%。若低于这个阈值，加大抗生素筛选浓度或缩短传代时间。
- 尝试用染色体整合型CRISPRi系统（如转座子介导）替代自由复制质粒，把dCas9或sgRNA表达框整合进靶标菌基因组，消除分离稳定性问题。
验证方法：传代第0天、第4天、第8天，菌液划线至含抗生素平板并挑取单菌落做菌落PCR验证sgRNA序列存在。同时测定培养物密度对应的保留率变化曲线，方便定位丢失节点。
避坑总结 + 实操建议
1. 定量校正必须先于生物学解释：无论CRISPRi效果多惊艳，没有校准过的qPCR或测序定量数据，一切结论都可被审稿人一击致命。外源内标DNA加进去，才叫真的生信校正。
2. 诱导效果检查排在所有实验前：一个aTc终浓度错误，时间加错，或者sgRNA设计本身效率低，会直接毁掉丰度差异，后面刷多少组重复都救不回来，纯培养先验证CRISPRi系统的干扰效率。
3. 引物专一性是丰度定量唯一的地基：琼脂糖切胶跑一跑，加标回收做一做，这两项花不了半小时，结果可信度是呈指数上升的。别把时间用来反复重做生物重复，先搞定引物。
日常维持一个基本习惯：每次CRISPRi干扰实验都写“诱导时间-取样量-提取试剂盒批次”记录表，培养条件的改变直接影响丰度定量结果，细节写下来能帮你快速回溯异常实验。另外一个习惯是准备一支备用靶菌株，冻存在-80°C甘油管中，每批实验都用同一株复苏扩培，避免不同批次的菌株生理状态不同的干扰。
实验方案里的引物设计或sgRNA靶点选择吃不准时，用站内的引物设计工具自动跑BLAST和二级结构筛选，大批量设计直接省下半天时间。菌群样本中含质粒高拷贝的干扰系统，做定量校准遇到多峰溶解曲线的，效率低不要硬凑，直接换靶向区域再做一轮，会比你想象中快得多。