细菌基因组中的菌株特异性位点：全基因组发现、基因组背景及在多菌株qPCR监测中的应用

《Microorganisms》：Strain-Specific Loci in Bacterial Genomes: Whole-Genome Discovery, Genomic Context, and Application for Multi-Strain qPCR Monitoring

　　

在复杂微生物群落中监测单个菌株仍是微生物生态学与生物技术领域的根本挑战。研究人员提出了一种整合式流程，用于在来自三个属（Stenotrophomonas、Bacillus 和 Pseudomonas）的四个具有生物技术相关性的植物促生菌株中鉴定并验证菌株特异性位点（SSL）。该流程采用两轮特异性过滤策略，结合全基因组比较与针对NCBI核苷酸数据库的高灵敏度BLASTn验证。SSL数量随分析所用菌株的平均核苷酸一致性（ANIb）升高而下降，从 Bacillus halotolerans（ANIb = 98.91%）的1个位点，到 Stenotrophomonas rhizophila（ANIb = 86.49%）的15个位点。全部25个SSL均普遍富含AT，主要关联于附属基因组，其侧翼区域富集功能未知基因（34.6%）与可移动遗传元件（19.2%）。靶向SSL的TaqMan qPCR检测体系表现出高特异性——在十份地理分布不同的土壤样本及一份原生根际宏基因组中均未检测到靶序列——且灵敏度高，基因组DNA检测限为0.01–0.1 pg。土壤加标实验的方法检测限（MDL）为850–15,000 CFU/g。田间试验施用菌群七天后，小麦根际检出全部四种菌株，验证了该流程在多菌株田间监测中的适用性。

研究背景与立项依据

在微生物生态学与生物技术应用中，对复杂群落中单个引入菌株的持久性、定殖与追踪能力是评估生物制剂功效的核心需求。传统物种级PCR或LAMP依赖保守区段，无法区分同物种内不同菌株；RAPD、AFLP、rep-PCR等指纹图谱法仅适用于纯培养；微卫星或ORF内搜索策略覆盖面有限。随着高通量测序成本下降，平均核苷酸一致性（ANIb）与四核苷酸相关性搜索（TCS）等工具使近缘基因组比对成为可能，但如何从全基因组连续序列中发掘真正“株级唯一”且适于qPCR的靶标，仍缺乏兼顾计算严谨性与田间验证的完整流程。研究人员以三属四株小麦根际促生菌为对象，建立并实验验证一套SSL发现—验证—监测管线，论文发表于《Microorganisms》。

主要关键技术方法

研究人员采用目标菌株全基因组与TCS筛选的两株最近缘基因组经KEC（K-mer Elimination by Cross-reference）算法初筛120–250 bp特异片段（SL），将BLASTn高频命中基因组补入参照组进行第二轮KEC过滤得候选preSSL，再以高灵敏度blastn-short参数对nt库终审得SSL；SSL侧翼±2 kb功能注释借助Pfam/InterProScan，可移动元件背景用PHASTEST预测；检测端基于SSL设计TaqMan探针，用纯化基因组做灵敏度曲线，用非灭菌田间土加标做方法检测限，用冬小麦田间施菌后第7天根际土做实证；转录活性以rsfS为内参行RT-qPCR；样本队列含四株模式菌纯培养、十份异地土壤、一份小麦根际宏基因组及田间施菌小区根际序列。

研究结果

3.1 全基因组序列中菌株特异性区域鉴定：四株菌ANIb介于86.49%–98.91%，SSL数随ANIb升高由15降至1，证实SSL主要源自附属基因组且与近缘度呈负相关。

3.2 SSL基因组环境分析：25个SSL均AT富集且GC低于宿主；12个完全位于CDS内，9个跨CDS边界，4个位于非编码区；侧翼基因34.6%功能未知，19.2%属可移动遗传元件，并反复出现限制修饰、RHS、C4反义RNA、HNH核酸酶等“冲突系统”组件，仅1个SSL落入前噬菌体区。

3.3 qPCR靶区与引物设计：每株选1个低E值、宜设探针对的SSL，偏移≤100 bp优化引物，Tm≥57℃，Sanger复核一致。

3.4 选定SSL转录活性：静止期RT-qPCR显示相对rsfS表达自0.45%（MGMM119，远距CDS）至266%（MGMM121，邻启动子5′UTR），证明非编码SSL仍可被邻近调控转录。

3.5 灵敏度评估：纯化DNA LOD为0.01 pg（三株）与0.1 pg（一株），合1–14基因组当量；土壤MDL为6.5×10³–1.5×10⁴ CFU/g，基因组拷贝/CFU系数0.48–2.12。

3.6 特异性检验：十份异地土壤总DNA无扩增，16S V3–V4 SYBR对照排除抑制。

3.7 宏基因组中SSL搜寻：小麦根际鸟枪数据未命中任一SSL，正交印证唯一性。

3.8 田间检测：施菌7天后根际四株全检出，定殖量2.5×10⁵–8.0×10⁶ CFU/g，超MDL 1–3个数量级。

讨论部分总结

研究人员指出，SSL数量不取决于数据库规模而取决于目标菌与最近测序近缘菌的生物学分化度；SSL普遍AT偏高契合自发突变偏倚与可移动元件整合热点，侧翼“冲突系统”提示其位于选择松弛的附属区；静止期低转录可由异源沉默蛋白（如H-NS）抑制AT丰富外源段解释，MGMM121因邻5′调控区而例外高表达；三段式特异性验证（nt库/BLAST、根际宏基因组、异地土壤）优于既往纯计算研究；土壤MDL与宏基因组阈值相当但更快更廉；基因组拷贝/CFU>1反映假单胞菌及芽孢杆菌寡倍体，偏差不超一个量级；SSL“可错唯一性”随数据库扩充可能失效，需周期重验，但短片段在数十至数百代内碱基置换极少，适于中长期监测。

结论部分翻译

本研究提出一套经实验验证的菌株特异性位点识别流程，基于全基因组KEC搜索、纳入BLAST筛选近缘菌的两轮过滤及高灵敏度BLASTn验证。以三属四株生物技术相关菌为案例，研究人员表明最终SSL数量主要由种内生物学分化决定，而该种内分化目前只能从近缘菌株已有基因组推断。本研究鉴定的全部SSL均为附属基因组中偏AT区段，转录活性低，提示其在松弛选择下向唯一性加速演化。基于SSL设计的TaqMan qPCR体系通过三级特异性验证，并在田间小麦根际可靠检出四株菌。就特异与灵敏而言，该方案可与宏基因组法比肩，但更快、更经济，有望成为农业环境研究与工业微生物应用中的常规菌株监测工具。该流程的普适性仍需在更广细菌类群中进一步验证。