枯草芽孢杆菌启动子筛选研究：为合成生物学与工业应用奠基

一、 研究团队与发表信息

本研究由来自中国天津工业生物技术研究所（Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences）、中国科学院系统微生物工程重点实验室（Key Laboratory of Systems Microbial Biotechnology）、国家非粮生物质能源工程技术研究中心（National Engineering Research Center for Non-food Biorefinery）以及广西科学院广西生物质产业化工程院（Guangxi Biomass Industrialization Engineering Institute, Guangxi Academy of Sciences）的研究人员共同完成。主要作者包括Yafeng Song、Jonas M. Nikoloff、Gang Fu、Jingqi Chen、Qinggang Li、Nengzhong Xie、Ping Zheng、Jibin Sun和Dawei Zhang（通讯作者）。该研究成果于2016年7月5日发表于学术期刊《PLOS ONE》上，文章标题为“Promoter screening from Bacillus subtilis in various conditions hunting for synthetic biology and industrial applications”。

二、 学术背景与研究目的

本研究属于合成生物学与代谢工程领域，具体聚焦于模式微生物枯草芽孢杆菌（Bacillus subtilis）的遗传工具开发。枯草芽孢杆菌因其高效的蛋白分泌能力和公认安全（GRAS）特性，在洗涤剂、纺织和制药等工业中具有广泛应用前景。然而，其作为合成生物学底盘细胞的潜力尚未被充分挖掘，主要原因之一在于缺乏经过系统表征和标准化的生物元件库，特别是启动子。启动子是调控基因表达强度和时序的核心元件，其性能的精确评估是实现代谢通路精细调控和异源蛋白高效表达的前提。

尽管已有研究对大肠杆菌和酵母等模式生物的启动子进行了大量筛选，但对枯草芽孢杆菌启动子的系统性、定量化评估仍相对匮乏。已有的研究多关注特定少数启动子，或依赖于转录组学（omics）数据的预测，缺乏在统一遗传背景下通过实验直接测量其驱动报告基因表达强度的全面分析。因此，本研究旨在填补这一空白。具体目标包括：1) 从枯草芽孢杆菌基因组中筛选并实验评估一批潜在的启动子序列；2) 在多种胁迫条件（高温、高渗、乙醇）下测试这些启动子的响应性；3) 验证所筛选启动子在驱动不同基因（胞内蛋白与分泌蛋白）表达时的通用性；4) 最终构建一个经过定量表征的启动子库，为利用枯草芽孢杆菌进行合成生物学设计和工业应用提供关键的标准化元件。

三、 详细研究流程

本研究流程严谨，主要包括以下几个核心步骤：

1. 启动子候选序列选择与载体构建：

   • 研究对象与样本量： 研究人员从枯草芽孢杆菌168菌株基因组中，选取了84个预测的启动子候选序列。这些序列位于不同类别基因的上游，包括热激蛋白、细胞膜蛋白、推测的金属抗性蛋白以及其他功能未知的蛋白编码基因。选择依据是它们在特定条件下（如胁迫）可能具有高转录活性。
   • 研究方法： 构建了一个启动子探针载体 pDL-gfp。该载体基于原有的 pDL 质粒（含有可整合到染色体 amyE 位点的同源臂），但将其原有的报告基因 bgab（编码热稳定性β-半乳糖苷酶）替换为 gfpmut1（绿色荧光蛋白突变体）。这种替换简化了高通量筛选过程。使用延伸重叠PCR（prolonged-overlap-PCR）技术进行基因替换，避免了寻找兼容性酶切位点的麻烦。随后，将84个约2 kb的启动子候选片段通过酶切连接的方式，分别克隆到 pDL-gfp 载体的多克隆位点，置于 gfp 基因上游。
   • 关键创新方法： 采用染色体单拷贝整合策略。将构建好的 pDL-gfp-pc（pc代表启动子候选）质粒线性化后，通过同源重组整合到枯草芽孢杆菌1A751菌株的染色体 amyE 位点。这种方法确保了每个细胞中报告基因的拷贝数恒定，避免了质粒拷贝数不均一可能带来的测量偏差，从而能更精确地反映启动子的本征强度。

2. 启动子强度定量评估：

   • 研究对象： 上述84个构建的枯草芽孢杆菌工程菌株，每个菌株携带一个特定的启动子-gfp 融合单元，整合在染色体相同位置。
   • 实验方法： 在标准LB培养基中培养工程菌株至对数中期（OD600 0.4-0.6），使用多功能酶标板读取器（microplate reader）测量全细胞荧光强度（激发波长488 nm，发射波长523 nm）。以不含启动子的菌株作为阴性对照，以含有已知强启动子 P43 的菌株作为参考基准（设定其相对强度为1）。所有实验均进行至少三次独立的生物学重复。
   • 数据分析： 计算每个启动子驱动的相对荧光单位（RFU），并归一化为相对于 P43 启动子的倍数变化。此外，还使用流式细胞仪（FACS）对部分高、中、低活性的启动子菌株进行了验证，以确保酶标板读数的可靠性。

3. 胁迫条件下的启动子活性分析：

   • 研究对象： 选择了一组被认为可能响应特定胁迫的启动子对应的菌株进行测试，例如热激蛋白基因（psigI, pgroES, pftsH）、细胞膜胁迫相关蛋白基因（pyuPA, pytpA, pmurF）以及推测的金属抗性蛋白基因（pyusI, pycbR, pywrK）的启动子。
   • 实验方法：
     • 转录水平分析： 使用实时荧光定量PCR（qRT-PCR）检测胁迫处理后 gfp mRNA的瞬时变化。处理包括：热激（43°C，15分钟）、乙醇胁迫（终浓度4% v/v，15分钟）、钴离子胁迫（CoCl2，终浓度8 mM，15分钟）。
     • 蛋白表达水平分析： 在更长时间的胁迫条件下（如43°C热激1小时，或含4%乙醇的培养基中培养6.5小时），测量GFP荧光强度的变化。同时，还在含有不同浓度NaCl（0.4 M， 0.6 M， 0.8 M）的高渗培养基中培养所有84个菌株，测量其对高渗条件的响应。

4. 启动子通用性验证：

   • 研究对象： 从84个启动子中挑选了10个活性范围不同的启动子（如 pradA, ptrnQ, psigX, psecDF 等）。
   • 实验方法： 将这些启动子克隆回原始的 pDL 载体（报告基因为 bgab），构建 bgab 报告系统。在相同条件下培养菌株，测定β-半乳糖苷酶活性。通过比较同一启动子驱动 gfp 和 bgab 表达的相对强度，评估启动子活性是否受下游报告基因类型的影响（即“上下文无关性”）。

5. 工业应用潜力验证：

   • 研究对象： 选择活性最强的启动子 ptrnQ 与常用的强启动子 P43 进行对比。
   • 实验方法： 将 ptrnQ 和 P43 分别插入分泌表达载体 pMA5 中，替换其原有的启动子，用于驱动分泌蛋白α-淀粉酶（α-amylase）的表达。在枯草芽孢杆菌1A751中进行摇瓶发酵（72小时），定期取样测定发酵上清液中的α-淀粉酶酶活，并通过SDS-PAGE分析蛋白分泌情况。

四、 主要研究结果

1. 启动子强度谱的建立： 测量结果显示，84个启动子候选序列的活性范围非常广，相对于 P43 启动子的倍数变化从0.0023到4.53不等，跨度接近2000倍。其中，ptrnQ（位于精氨酸tRNA基因上游）的活性最强，是 P43 的4.53倍。psigX 和 pgroES 也表现出较强的活性，分别为 P43 的3.03倍和1.55倍。同时，也发现了许多活性很弱的启动子，其中一些位于功能未知的假设蛋白基因上游。序列分析显示，这些强启动子的-10区与经典的“TATAAT”序列相似，但-35区与“TTGACA”一致性较差，且间隔区（spacer）长度并不总是标准的17 bp，这可能解释了其强度差异。

2. 胁迫响应的有限性： qRT-PCR结果显示，短时间（15分钟）的特定胁迫处理确实能引起部分启动子转录水平的显著变化（如 psigI 和 pgroES 在热激后mRNA水平升高，pftsH 和 pmurF 则下降）。然而，在更长时间（小时级）的胁迫条件下，GFP蛋白表达水平的变化远不如mRNA水平变化显著。例如，在43°C热激1小时后，大多数启动子驱动的GFP表达仅有微小变化。在含4%乙醇的培养基中培养，也未观察到显著的GFP表达变化。最显著的发现出现在高渗条件下：当在含0.8 M NaCl的培养基中培养时，84个启动子中有11个（如 pmreB, ppssA, poatA, pytpA 等）的活性出现了相对明显的上调。这表明，尽管许多预测的“胁迫响应”启动子对急性转录刺激有反应，但在持续胁迫下，其蛋白输出水平的调控可能更为复杂或受到全局性调控的掩盖。

3. 启动子活性的上下文无关性： 对10个不同强度启动子的验证实验表明，同一启动子驱动 gfp 和 bgab 表达的相对强度呈良好的线性相关（R²值高）。这一结果至关重要，它证明这些启动子作为标准化“生物砖块”（BioBrick）的潜力，即其强度特征在不同基因背景下是相对稳定和可预测的，这为它们在合成生物学电路中的模块化应用提供了实验依据。

4. 强启动子在工业蛋白生产中的应用验证： 应用实验证实，使用最强启动子 ptrnQ 驱动α-淀粉酶分泌表达时，其最终酶活产量是使用 P43 启动子时的约2.1倍。SDS-PAGE结果也直观显示了 ptrnQ 菌株分泌的α-淀粉酶条带更浓。值得注意的是，这个提升倍数（2.1倍）低于在胞内报告基因（gfp 和 bgab）中观察到的差异（约4倍）。作者分析认为，这可能是由于分泌蛋白的表达受到翻译后加工、转运、折叠等多重瓶颈的限制，导致启动子强度与最终分泌产量之间并非简单的线性关系。尽管如此，该结果仍明确证明了 ptrnQ 是一个优于常用强启动子 P43 的候选元件，可用于提高异源蛋白的产量。

五、 研究结论与价值

本研究的核心结论是成功构建并定量表征了一个来自枯草芽孢杆菌的、包含84个成员的基础启动子库，揭示了其广泛的强度范围和部分启动子对高渗条件的响应特性。研究证明了所筛选启动子在不同遗传背景下的功能一致性，并将最强启动子 ptrnQ 成功应用于提高分泌蛋白的产量。

科学价值： 1. 数据资源： 提供了枯草芽孢杆菌大量内源启动子的首次系统性实验定量数据，为理解该菌的转录调控网络提供了宝贵的实证资料。 2. 工具开发： 建立了一个经过验证的启动子“工具箱”，其强度跨度大（近2000倍），为合成生物学中实现基因表达的精细、梯度调控提供了关键元件。 3. 机理洞察： 揭示了转录水平响应与最终蛋白表达水平响应之间的不一致性，强调了在评估条件诱导型启动子时进行多层级（转录组、蛋白质组）验证的重要性。 4. 方法学贡献： 展示了基于染色体单拷贝整合和荧光报告基因的高通量启动子筛选流程的可靠性与高效性。

应用价值： 1. 代谢工程与合成生物学： 该启动子库可直接用于枯草芽孢杆菌的代谢通路优化，通过匹配不同强度的启动子与通路中各个基因的需求，实现代谢流的平衡，从而提高目标化合物（如维生素、氨基酸、生物燃料）的产量。 2. 工业酶生产： 鉴定出的超强启动子（如 ptrnQ）和潜在胁迫诱导型启动子，为工业酶（特别是分泌型酶）的高效、可控生产提供了新的遗传调控元件。 3. 标准化生物元件库建设： 这项工作是对iGEM等机构推动的生物元件标准化运动的重要补充，将枯草芽孢杆菌这一重要工业微生物纳入了标准化“生物砖块”的版图。

六、 研究亮点

1. 系统性与定量化： 首次对枯草芽孢杆菌如此大规模（84个）的内源启动子候选序列进行了统一的、基于实验的定量强度评估，并建立了完整的强度谱。
2. 严谨的评估体系： 采用染色体单拷贝整合策略，消除了质粒拷贝数变异的影响，确保了数据准确反映启动子的本征强度。同时使用两种报告基因（gfp 和 bgab）进行交叉验证，证明了启动子性能的上下文无关性。
3. 多维条件测试： 不仅评估了启动子在标准条件下的强度，还系统测试了其在多种胁迫条件（热、乙醇、高渗、金属离子）下的响应，提供了更全面的性能图谱。
4. 从元件表征到应用验证的闭环研究： 研究并未止步于元件的体外表征，而是进一步将筛选得到的最强启动子应用于实际的工业酶分泌表达中，验证了其应用潜力，完成了从基础研究到应用探索的完整链条。
5. 发现强效新元件： 鉴定出 ptrnQ 这一活性显著高于常用强启动子 P43 的新型强启动子，为枯草芽孢杆菌的高水平表达系统提供了新的有力工具。

七、 其他有价值内容

研究还通过生物信息学工具（如SoftBerry， WebLogo）对最强启动子的序列特征进行了分析，探讨了其-10区、-35区及间隔区的特点，虽然未发现与经典序列完全一致的保守模式，但为后续通过理性设计或定向进化改造启动子提供了初步的序列基础。此外，文中对结果讨论部分深入分析了胁迫响应不显著的可能原因（如全局调控掩盖特异性响应、mRNA与蛋白水平调控脱节等），体现了研究者对微生物生理的深刻理解，对后续相关研究具有启发意义。