枯草芽孢杆菌(Bacillus subtilis)作为一种模式菌株，具有生理生化和遗传背景清晰、非致病性和分泌性良好的优良特性，因此被用作食品安全级别的微生物。随着Bacillus subtilis 168菌株全基因组测序的完成，越来越多的外源基因已经在枯草芽孢杆菌及其亲缘菌株中成功高效表达，这使其成为基因表达系统的研究热点之一[1]。

相比大肠杆菌表达系统，枯草芽孢杆菌的表达系统研究起步较晚，其本身具有一些独特的优势[2]：1)安全无毒，已被美国食品药品监督管理局(FDA)评为生物安全(Generally regarded as safe，GRAS)菌株；2)较强的蛋白分泌能力，方便下游表达产物的分离；3)发酵工艺成熟，容易进行高密度发酵。但是枯草芽孢杆菌表达系统也存在许多不足之处，主要表现在[3]：启动子转录水平低，某些目的蛋白不能大量表达；重组蛋白折叠不完全，部分没有活性；菌株感受态细胞转化率低，内源性质粒结构不稳定；表达产物在胞外容易被本身分泌的蛋白酶降解。这些因素限制了枯草芽孢杆菌表达系统不能完全满足工业化生产的需求，因此构建新的蛋白表达系统仍是研究枯草芽孢杆菌的重点。

稳定的载体骨架和高效可控的启动子是保证工业上枯草芽孢杆菌高效表达外源蛋白的关键因素，本文对枯草芽孢表达系统的特点做简单的介绍，着重介绍表达系统中的转录调控元件——启动子，详细叙述启动子的种类，新启动子的克隆方法，改造思路以及启动子的预测与分析。

克隆载体和表达载体在分子遗传操作领域中具有十分重要的作用，是基因重组技术的重要工具。根据类型的不同可以将枯草芽孢杆菌表达载体分为3类[4]：独立自主复制型载体、整合型载体和噬菌体载体。

枯草芽孢杆菌内源性质粒DNA含量很少，现有的枯草芽孢杆菌表达系统采用的质粒主要来源于葡萄球菌和链球菌[5]等其他微生物。根据不同的复制方式可以将这些质粒载体分为滚环复制型载体和θ复制型载体。只有一些大的质粒载体以θ复制方式进行复制，而一大部分质粒载体则通过滚环复制方式复制，这些质粒载体大小均在在12 Kb以下。来源于金黄色葡萄球菌的质粒pC194、pE194、pT181和pUB110均是以滚环复制方式进行复制的。

pC194大小为2 906 bp，抗氯霉素特性，质粒的拷贝数约15个染色体[6]，该质粒较小，容易克隆，当插入外源片段时容易产生大分子的畸形复制产物，从而高度不稳定，这就决定了该质粒载体在宿主中的应用性偏低；pE194大小为3 728 bp，温度敏感型质粒载体，拷贝数大约为10个染色体，该质粒载体在32℃正常复制，随着温度的升高拷贝数会逐渐降低，达到45℃时停止拷贝，所以在枯草芽孢杆菌中也不容易稳定地遗传；pT181具有抗四环素特性，每个细胞中拷贝数约为20个染色体[2]；pUB110大小为4 548 bp，抗卡那霉素和博来霉素，每个细胞约有30~50个拷贝数[7]，该质粒已被完全测序，但是因其容易产生一些不稳定的衍生物，所以应用价值也不大。

以上这些种类载体能在枯草芽孢杆菌中独立自主地复制表达其抗性基因，但是在传代培养中会出现结构不稳定性和分离不稳定性的问题[8]。载体的不稳定性严重限制了其在研究工作中的实际应用，由此造成在枯草芽孢杆菌中分子遗传操作的不便。

由于枯草芽孢杆菌提取质粒和构建质粒较为困难，而在大肠杆菌中操作则很简单方便，所以可通过构建大肠杆菌-枯草芽孢杆菌穿梭载体，先将外源目的基因在大肠杆菌中构建完成后再转入枯草芽孢杆菌中表达。目前常见的大肠杆菌-枯草芽孢杆菌穿梭载体有pEB10[9]、pEB20[10]、pEB60[10]、pUB18[11]、pUB19[12]和pWB980[13]等，详情见表 1。

游离型质粒的不稳定性是制约枯草芽孢杆菌高效表达外源蛋白的主要障碍，一般是通过构建整合型载体来克服这一问题。构建枯草芽孢杆菌的整合型载体主要部分包括：大肠杆菌质粒pBR322或其衍生质粒的复制起点，供筛选的抗性标记基因，一段到两段与枯草芽孢杆菌基因组序列同源的基因序列。整合型载体也需要先构建大肠杆菌-枯草芽孢杆菌穿梭载体，再导入枯草芽孢中进行同源交换把表达单元整合到枯草芽孢杆菌染色体上，根据整合方式的不同可以将整合型载体分为同源单交换整合和同源双交换整合。

(1) 同源单交换整合

同源单交换整合是将外源基因某侧和枯草芽孢杆菌基因组同源的序列与目的序列进行同源重组，利用这种方式将外源基因整合到载体上面。单交换时同源臂不要低于25 bp，同源臂达到70 bp时整合效率较高[14]。

(2) 同源双交换整合

同源双交换整合是将外源基因两侧的同源序列与枯草芽孢杆菌基因组目标序列进行同源重组，为了发生有效的同源重组，同源臂必须要在400~500 bp[15]，其次还要考虑重组DNA的线性稳定性。一般在构建重组载体时在同源臂中加入Chi位点就可以解决该类问题[4]。

在基因转录过程中RNA聚合酶具有重要的作用，它能特异性识别启动子并且与之结合后转录才能开始。RNA聚合酶由5个亚基构成：α亚基、β亚基、β'亚基、ω亚基和一个Sigma(σ)因子[16]。全酶包含5个亚基，核心酶不包括σ因子，σ因子的作用是帮助核心酶识别和结合到启动子的转录起始位点上面。枯草芽孢杆菌表达基因时不同生长阶段的基因表达不同，σ因子出现的种类和数量也不同。目前已知的枯草芽孢杆菌σ因子一共有14种，这些σ因子功能不尽相同，有的与芽孢的形成有关，有的与非重叠、重叠启动子和转录级有关[17-18]，具体描述如下(表 2)。

大肠杆菌只含有7种σ因子(σ70、σ54、σ38、σ32、σ28、σ24、σ19)，与枯草芽孢杆菌有极大的差别。因此，枯草芽孢杆菌的启动子在大肠杆菌中可能很容易被识别并且表达，反之，大肠杆菌中表达原件大多数不能在枯草芽孢杆菌中使用，因为很难被识别，或者活性低甚至没有活性。

启动子作为表达系统的一个重要元件，对基因的高水平表达起着关键的调控作用。原核启动子的转录起始位点上游200 bp到下游100 bp范围为启动子的核心区域，该区域的显著特征就是存在小段的保守序列，可以用来预测和分析启动子。对于典型的原核启动子，其保守区域特征有转录起始位点(TSS)、Sextama-35区(TTGACA)、Pri-bmow-10区(TATAAT)[19]，Sextama-35区和Pri-bmow-10区序列之间的间隔距离(图 1)。

转录起始位点TSS一般超过90%是嘌呤碱基，CAT序列的中心。几乎所有的启动子TSS上游10 bp处都有一个6 bp的保守区域，这就是Pri-bmow-10区，详细情况为T80A95T45A60A50T96(下标为出现该碱基的最大频率百分数)。由此可知Pri-bmow-10区序列中开始的两个碱基TA和最后一个碱基T高度保守，为该区域最重要的碱基。另外一个6 bp的保守区域为Sextama-35，详细序列为T82T84G78A65C54A45，其碱基的保守性较Pri-bmow-10低[20]。Sextama-35区和Pri-bmow-10区的间隔序列一般在(17±1)bp，也有的小于15 bp或者大于20 bp，尽管间隔区域的序列构成并不重要，但是距离却十分重要，间隔太小或者太大均会不同程度地影响启动子的强度，该距离对影响着RNA聚合酶的几何构像。综上得出，启动子核心区域的保守序列按顺序排列，间隔一般在17 bp左右，这个距离被称为最佳距离。

除此之外，原核启动子的区域内还存在这些“非典型”元件：Pri-bmow-10区延伸序列[21]、UP元件[22]、TSS下游AT富集区[23]，与典型特征相比它们的保守性要差一些。实际上，各段特征序列会有灵活的变化。除了Pri-bmow-10区的保守性强一些以外，几乎没有其他哪段序列是必须的。因为丢失了某段特征性弱的序列而损失的启动子活性，可以通过加强其他位置进行一定程度的弥补。如在枯草芽孢杆菌中，丰富某些启动子Sextama-35区上游的某处的AT含量可以增加转录强度[24]。

核糖体结合位点(Ribosomebinding site，RBS)，位于起始密码子AUG上游的一段富含嘌呤的DNA序列，也称作Shine-Dalg-arno (SD)序列，该序列在1974年被Shine J和Dalgarno L一起发现[25]，所以以此命名。细菌RBS位于细菌mRNA起始密码子AUG上游10 bp的一段富含嘌呤的序列，该序列能识别核糖体16S rRNA并与之进行互补配对，促使mRNA与核糖体结合，从而促进翻译的起始。大肠杆菌的RBS序列为“AAGG”，枯草芽孢杆菌的RBS序列为“GGAGG”，其中以GGAG 4个碱基比较重要。Kozak[26]研究发现突变起始密码子AUG周围的碱基在一定程度上会影响转录和翻译的效率。有研究者通过随机突变RBS提高了目的蛋白的表达量；也有研究者通过对RBS的长度进行增加或者删减达到了高水平调控目的蛋白表达的目的。在线工具RBS Calculator 2.0可以对启动子自身的RBS进行预测和设计。目前利用该软件根据需求设计的RBS进行实验取得了不错的效果[27]。

枯草芽孢杆菌表达载体常用的启动子按照其诱导的方式可以分为组成型和诱导型启动子[28]。

组成型启动子是指不需诱导也可以持续表达的启动子，已经报道的枯草芽孢杆菌启动子有Papre、P43、PHpaⅡ[29]、Plapsd等。其中P43启动子为强启动子，在枯草芽孢杆菌中应用最为广泛，它来自枯草芽孢杆菌胞苷脱氨酶(cdd)基因的启动子，有两个重叠的Sextama-35区和Pri-bmow-10区，分别由σA和σB识别[30]。以W800为宿主，以pMA0911-PHpaⅡ为载体表达支链淀粉酶(PUL)，最高酶活力可达3.9 U/mL。通过替换P43启动子后pMA0911-P43-PUL同等条件下最高酶活力可达8.7 U/mL[31]。Yang等[32]筛选到一个来自地衣芽孢杆菌的启动子Pshuttle-09，强度为P43的8倍，又在Pshuttle-09基础上构建了一个双启动子Plaps，强度为P43的13倍，是目前在枯草芽孢杆菌表达中启动能力最强的组成型启动子。组成型启动子能高质量表达目的蛋白，成本较低，但同时也存在一些不足之处，如利用组成型启动子表达目的蛋白时不容易控制表达量，而且组成型启动子构建成的载体不稳定，当用一个组成型启动子表达多个外源蛋白的时候，容易造成共抑制和基因沉默等问题。

诱导型启动子不同于组成型启动子，它具有独特的优势，可以根据需求给予特定的诱导剂。根据诱导方式可以将诱导型启动子分为两类：1)环境诱导型启动子，该类启动子受低温、热休克、氧气匮乏、渗透压、pH值等因素的影响而做出相应的基因表达来适应环境的变化；2)化学诱导型启动子，可以通过添加诱导物来调节目的基因的表达量。

最常见的是受异丙基-β-D-硫代半乳糖苷(IPTG)诱导的启动子Pspac，该启动子由大肠杆菌lac操纵子和枯草芽孢杆菌噬菌体SPO1融合而成[33]。随后Phan等[34]又在Pspac的基础上优化其调控元件，得到一个启动能力更强的启动子Pgrac100。这些启动子在IPTG的诱导下能高效表达目的蛋白，而且可以控制，但仍存在一些问题：1)IPTG具有毒性，不能用于表达食品或者医疗用途的蛋白产品，而且其价格昂贵，成本偏高；2)Pspac启动子表达外源蛋白量不够理想；3)Pspac启动子不具有严谨可控性，在没有诱导剂添加的时候仍有少许的目的蛋白进行表达。以上这些因素限制了IPTG诱导型启动子在工业上的大规模应用。

以木糖为诱导剂的启动子Pxyl也是枯草芽孢杆菌中常用的诱导型启动子。该启动子来源于巨大芽孢杆菌木糖操纵元件，严格受xylR编码的阻遏蛋白控制。木糖启动子对大部分代谢产物不敏感，受葡萄糖的抑制[35]。木糖价格昂贵，工业应用中成本较高。

蔗糖启动子PsacB是枯草芽孢杆菌第3个常用的启动子系统。该启动子是枯草芽孢杆菌基因组上SACB基因编码的蔗糖果聚糖酶的启动子。PsacB启动子的转录不受蔗糖的严格调控，在没有诱导剂时也能进行转录，但较有蔗糖诱导时强度低100倍[36]。与前两个诱导型启动子相比PsacB启动子强度较弱，不符合工业上大量表达目的蛋白的要求。

除此之外，麦芽糖诱导型启动子Pglv也是枯草芽孢杆菌中应用广泛的诱导型启动子，包括以下3个部分[37]：1)6-磷酸-α-葡萄糖苷酶(glvA)基因，含有449个氨基酸，在NAD(H)(还原型辅酶Ⅰ或辅酶Ⅱ)和二价金属离子存在的情况下才具有活性。2)调控蛋白(glvR)基因，含有254个氨基酸，与RpiR/YebK/YfhH家族相似性很高。3)磷酸烯醇丙酮酸酯-糖磷酸转移系统透性酶(glvC)基因，是一个有12个跨膜区域和527个氨基酸组成的转运酶。Pglv启动子上存在典型的代谢反应元件cre序列，该序列位于启动子的转录起始位点和核糖体的结合位点，它的存在使启动子受葡萄糖特异性抑制，如图 2。

Pglv启动子以麦芽糖为诱导物，价格便宜，工业成本低，安全无毒，使用价值比较大。但是发酵中产生的葡萄糖会对启动子有强烈的反馈抑制，这是工业应用上面急需解决的一个问题。

枯草芽孢杆菌启动子的克隆方法大体可以分为启动子探针筛选法和PCR技术克隆法。

(1) 启动子探针法筛选启动子

启动子探针又称启动子诱捕系统，核心部分的探针质粒载体由一个筛选标记基因和不含启动子成分的报告基因构成，是一种简单有效而且精准的筛选启动子的方法[39]。筛选启动子的时候用合适的限制性内切酶来消化基因组DNA，启动子探针载体也同时用相同的酶切处理，然后将酶切处理好的基因组片段和相同方式处理的启动子探针连接，根据报告基因的表达与否、表达的强弱来确定启动子的筛选情况[40]。在构建探针中报告基因的选取也至关重要，一般要求报告基因序列已知，目的蛋白特性清晰，检测方便快捷，且在宿主中表达时不存在背景干扰。目前常用的报告基因有[41]：1)β-半乳糖苷酶基因，能够水解底物X-gal，可以通过蓝白斑来筛选具有启动子的阳性克隆；2)氯霉素转乙酰酶(CAT)基因，该酶可以使宿主失去抗性，通过检测放射性可确定启动子插入情况；3)绿色荧光蛋白(GFP)基因，无需底物和辅助因子，可以在生物体内高效表达并且发出荧光，对机体无毒性、稳定、检测直观便捷，被广泛用于原核生物和真核生物中[16]；4)抗性基因，不需要底物，通过设置抗性梯度可以筛选出强弱不同的启动子。

利用启动子探针系统筛可随机筛选到大量的启动子，简单轻松，但是该方法也存在着工作量大，需要构建启动子文库，后期还要对启动子进行分析等问题。

(2) PCR技术克隆法筛选启动子

该方法结合了现代生物学技术和PCR技术，先用计算机相关软件对基因序列进行分析和预测，然后利用PCR技术获取已知序列中的启动子，该方法逐渐成为获取启动子的研究热点。

目前国内外有很多利用这种方法获取启动子的报道，Wang等[42]利用生物信息学预测到了一些粗略的RNA的转录起始位点，预测分析了该转录起始位点上游500 bp处的序列，并用pri-miRNA PCR(pPCR)技术克隆出该片段，得到了所要的启动子，并在此基础上筛选获得苏云金芽孢杆菌20多个不同时期表达的启动子；黄君健等[43]利用Panhandle PCR技术扩增得到人体外周血单核细胞端粒酶纯化亚基hTERT上游2 090 bp处的DNA序列；王新国等[44]通过将胡萝卜的基因组酶切后，将酶切片段作为模板，利用衔接头PCR方法扩增，成功得到一个新的sⅡ基因启动子序列。

除了以上介绍的几种PCR扩增技术外，常用于扩增启动子序列的PCR技术还有定向反向PCR、交错式热不对称(PCR TAIL-PCR)、Universal Fast Walking(UFW)等。利用PCR技术获得启动子的速度快、特异性强，避免了DNA本身的连接和环化的问题。其缺点是对启动子的预测和分析要求比较精准。

启动子自身构造决定着启动子的转录能力强弱，目前用来表达外源基因的启动子一般是天然的启动子。研究实验表明，改造启动子的主要结构区域，可以不同程度提高启动子的活性[45-46]。目前主要从以下两个思路对启动子进行改造：

(1) 将不同启动子的保守序列进行替换、杂交成为新的启动子

启动子典型的核心保守区域有：①Sextama-35区(Sextama盒)，它是RNA聚合酶的识别位点或松弛型结合位点，该区域的结构决定了它和RNA酶的亲和性，因为RNA聚合酶更容易识别强启动子；②Pri-bnow-10区(Pribnow盒)是启动子与RNA聚合酶的紧密结合部位，决定着转录的方向，RNA聚合酶在此区域与DNA序列形成稳定的结合物；③Sextama-35区和Pri-bnow-10的间隔区，该区的序列类型不重要，但是该区的间距长度决定了RNA聚合酶的构型，间距趋近17 bp时启动子启动能力表现为上调，远离17 bp时表现为下调趋势。随着现代生物信息学的发展，我们可以对得到的已知启动子进行预测分析，得出每个启动子核心区域中保守区域的具体情况，然后可以将多个启动子进行优势互补的杂交，融合每个启动子最有优势的核心的保守区域于一个启动子中，从而使启动子的能力得到提高，获得一加一大于二的效果。

(2) 将启动子的特征区域突变为典型启动子的保守区域

不同的启动子都含有保守的区域，与典型的原核启动子相比保守区域可能有所出入。Pri-bnow-10区详细情况为T80A95T45A60A50T96(下标出现该碱基的的最大频率百分数)，最开始的两个碱基TA和第6位的碱基T最为保守。Sextama-35区详细序列为T82T84G78A65C54A45, 其中前3个碱基TTG保守性较好。有研究表明通过将启动子的核心区域突变为保守区域，启动子的活性相应得到提高，重组蛋白的表达量也相应地提高了一倍[46]。

启动子分析和预测方法[47]有很多，如生物信息学预测法、酵母单杂交实验法、染色质免疫共沉淀法、瞬时转染法和凝胶阻滞分析实验法，以下主要介绍生物信息学预测法。

目前，启动子的预测软件有Cpgplot、CpG Prediction、EMBOSS、NNPP、Promoter 2.0、Promoter Scan等。此外，还需要参考许多数据库，如：1)真核生物启动子数据库(EPD)，主要收集和RNA聚合酶Ⅱ相关的启动子数据，用户只需输入核苷酸序列就可以查找到启动子的位置，数据库会提供该数据的相关数据链接、转录起始位点以及原始文献[48]；2)植物顺式作用调控元件(Plant-CARE)数据库[49]，收录的数据有植物顺式作用元件、植物的抑制子和增强子等；3)转录因子(TRANS-FAC)数据库，数据来自人类和酵母的转录因子，是有关转录因子与基因组结合位点的数据库[50]；4)Bacterial Promoter, Operon and Gene Finding(BPROM)(bacterial)，该数据库能够预测原核生物启动子的结构，此外还具有Single Nucleotide Polymorphisms(SNP)筛选、剪切位点寻找、表达谱分析等功能。现在很多在线工具可以用来分析和预测基因组的启动子、Untranslated Region(UTR)、内含子等结构，如表 3。

生物信息学提供了大量的参考信息，为在线预测和分析启动子奠定了坚实的基础[51]。通过生物信息学的相关软件分析和预测启动子的相关调控元件及其序列，为深入研究启动子结构和功能提供了便利[52]。

启动子是枯草芽孢杆菌高效表达目的蛋白的关键元件之一，近年来，国内外开展了大量的有关启动子的研究，并取得了较大的进展，获得大批可用于枯草芽孢杆菌的启动子，但是就现有的启动子而言，可以完美用于工业上的启动子数量有限，启动能力和调控方式存在着诸多问题。要获得更多转录强度高、诱导机制方便的启动子，还需要不断深入的研究。

启动子的研究是分子遗传操作十分重要的工作，对于研究枯草芽孢杆菌基因组的功能也有着重要的意义。未来启动子的研究主要集中在两个方面：1)新型高效可控启动子的的发现，这是解决工业和市场生产高要求的最佳的措施；2)对现有启动子的改造，启动子是控制目的基因表达的开关，如果能对启动子进行精确的定位，弄清其各部分元件的工作原理及意义，然后加以改造就能使启动子利用效果最大化。