全局转录调控因子CodY在嗜热链球菌中的调控机制研究

【摘要】嗜热链球菌（Streptococcusthermophilus）是重要的工业菌种,与德氏乳杆菌保加利亚亚种（Lactobacillusdelbrueckiisubsp.bulgaricus）或者瑞士乳杆菌（Lactobacillushelveticus）混合发酵生产酸奶和硬奶酪。嗜热链球菌年消费量达1021菌体细胞,是公认的食品级安全性菌株（GenerallyRecognizedasSafe,GRAS）。长期生长于牛乳环境,嗜热链球菌形成了独特的乳糖转运和代谢系统,使其能够高效利用牛乳中的乳糖为菌体生长提供能量。由于嗜热链球菌的胞外蛋白酶活性较弱或者缺失,不能有效利用牛乳中的酪蛋白；同时牛乳中游离氨基酸和寡肽浓度低,限制了其生长。为了更好的适应牛乳环境,嗜热链球菌通过基因水平转移（HorizontalGeneTransfer,HGT）获得了多种氨基酸合成途径。此外,嗜热链球菌在生长和发酵过程中常常受到氧、温度和噬菌体等胁迫压力,影响菌体的生长和存活。为了能够充分利用环境中的营养物质和抵抗胁迫环境,细菌通过全局转录调控因子对胞内基因进行转录控制,形成一个庞大的代谢调控网络。全局转录调控因子CodY作为一个独特的转录调控蛋白普遍存在于革兰氏阳性菌,在细菌代谢、环境适应和致病性控制等方面都具有重要的调控作用。CodY通过与调控基因的启动子区域结合控制结构基因的转录,而能量信号分子鸟嘌呤三磷酸（Guanosine-5’-Triphosphate,GTP）和营养信号分子分支氨基酸（BranchedChainAminoAcid,BCAA）能够有效促进CodY和调控基因的作用,但是不同菌株CodY的作用机制不同。嗜热链球菌碳、氮代谢途径研究相对清楚,但是有关CodY在调节碳、氮代谢以及抗胁迫中的作用尚未见报道。通过比较基因组分析发现,嗜热链球菌基因组上存在一段保守的蛋白编码序列codY,推测其编码蛋白CodY具有重要的调控功能。本论文聚焦S.thermophlisST2017全局转录调控因子CodY的碳、氮代谢调控及抗胁迫机制,通过转录组分析,建立CodY对碳、氮代谢的调控网络,利用凝胶阻滞电泳（ElectrophoreticMobilityShiftAssay,EMSA）和等温滴定量热法(（IsothermalTitrationCalorimetry,ITC）分析其调控机制,阐明CodY在抗胁迫中的调控作用,开展工作及取得结果如下：1、S.thermophlisST2017全局转录调控因子CodY的调控网络分析通过比较基因组分析,确定S.thermophlisST2017基因组上的codY基因序列。利用温度敏感型质粒pG+host9,通过同源双交换敲除S.thermophlisST2017基因组上的codY基因,得到突变株S.thermophlisST2017△codY。通过转录组测序,分析比较野生株ST2017和突变株ST2017△codY基因差异表达谱,并通过荧光定量PCR验证。结果显示,codY基因敲除后共有87个基因发生了差异表达,其中62个基因的转录受CodY抑制,对这些差异表达基因进行分类,发现CodY调控以下代谢类型的基因转录：（1）控制氨基酸合成与转运：如分支氨基酸的合成（ilvA、LeuABCD和ilvBHC）及转运（livJHMGF）,色氨酸（TrpEGDCFBA）、半胱氨酸的合成（cysM2metB2cysE2）,天冬氨酸、丝氨酸和苏氨酸的合成,组氨酸分解代谢中的鸟苷水合酶基因（hutU）则是氨基酸代谢中唯一受CodY转录激活的基因；（2）碳代谢相关基因的转录,如EMP途径中的甘油醛3一磷酸脱氢酶和参与半乳糖代谢的磷酸变位酶；（3）CodY对谷氨酸脱氢酶基因gdhA具有强烈的转录抑制作用；（4）控制了多种ABC转运蛋白和二价离子转运蛋白的转录如STND0202-0203、STND1235-1232和STND0733;（5）激活次黄嘌呤（IMP）合成操纵子purCLFhNH的转录,并抑制尿嘧啶合成操纵子的转录；（6）参与环境胁迫相关基因的转录调控,如伴侣蛋白DnaK和GroESL操纵子在codY突变株中显著上调,与渗透压胁迫相关的STND0113、STND0135和与氧胁迫相关的STND1346则发生了转录下调,同时一些参与环境胁迫应答的转录调控因子也受CodY调控；（7）抑制二组分系统TCS01的转录；（8）其他基因的转录,在差异表达谱中,有大量基因的功能没有确切的注释。由此看出,CodY参与了细胞的多种代谢途径的调控,尤其对氮代谢的调节更为广泛,同时还参与调控环境胁迫相关基因的转录。2、S.thermophilusST2017全局转录调控因子CodY的作用机制研究在革兰氏阳性菌中,转录调控因子CodY通过和细胞内GTP和BCAA的协同作用来调控目的基因的转录,为了确定S.thermophilusST2017中CodY的作用机制,选取在转录组中明显上调的分支氨基酸转运蛋白操纵子（livJHMG）为研究对象,通过EMSA检测CodY和livJHMG启动子的相互作用。结果显示,CodY和livJHMG启动子具有明显的结合作用,当反应体系中添加2mM的GTP时,CodY和启动子结合能力没有变化,而15mM的BCAA能够有效促进CodY和启动子的结合,说明在S.thermophilusST2017中,GTP不能促进CodY和启动子的结合,BCAA是唯一的激活效应物,其作用模式类似于乳酸乳球菌（Lactococcuslactis）和肺炎链球菌（Streptococcuspneumoniae）。为了证明效应物BCAA的正调控作用,将livJHMG启动子和绿色荧光蛋白基因融合后在S.thermophilusST2017中表达。结果显示,BCAA能够有效降低菌体的荧光强度,从体内进一步证明了BCAA作为CodY效应物的正调控作用。CodY与调控基因的结合具有高度选择性,为了检测CodY与调控基因的结合位点,选取转录组测序中具有显著差异表达的24个基因,通过MEME预测其启动子区域的序列特征,结果发现存在一段保守的DNA序列A[AT]T[AT]TTC[TC]GA[ACT]A[AT]TT,该序列和L.lactis及枯草芽孢杆菌（Bacillussubtilis）的CodY结合域（AATTTTCNGAAAATT）具有高度的序列相似性,推测S.thermophilusST2017的CodY通过该保守序列和调控基因的启动子结合。为了验证这一假设,通过ITC实验检测CodY和该保守序列DNA的结合能力,结果证明嗜热链球菌CodY和保守序列具有显著的结合作用,并且以二聚体形式和保守序列按照化学计量数1：1结合；将保守序列的碱基突变后,这种结合作用随之消失,说明该保守序列在CodY的识别中必不可少,而BCAA能够有效提高CodY和保守序列的亲和力,使得解离常数Kd从780nM降低到8.33nM。3、S.thermophilusST2017转录调控因子CodY在碳、氮代谢转化中的作用不同于其他乳酸菌,嗜热链球菌具有完整的从丙酮酸到α-酮戊二酸的代谢途径,并且α-酮戊二酸在谷氨酸脱氢酶的作用下可以生成谷氨酸,实现碳代谢和氮代谢的相互转化。转录组学分析发现,S.thermophilusST2017△codY中谷氨酸脱氢酶基因（gdhA）的转录水平是野生菌株ST2017的3倍,说明转录调控因子CodY对gdhA具有明显的抑制作用。EMSA进一步证明CodY对gdhA启动子具有直接的转录调控作用。体外测定S.thermophilusST2017及ST2017△codY的谷氨酸脱氢酶酶活,野生菌株的酶活为190U,而突变株为330U,说明codY基因的缺失解除了CodY对gdhA的转录抑制；将codY基因回补到突变株ST2017△codY,并使之高效表达,谷氨酸脱氢酶的酶活降低到100U,进一步证明了CodY对gdhA表达的抑制作用。在嗜热链球菌中,谷氨酸脱氢酶催化α-酮戊二酸合成谷氨酸是一个可逆反应,但催化两种反应的Km值不同,α-酮戊二酸合成谷氨酸的Km值为0.131±0.030mM,远低于谷氨酸降解为α-酮戊二酸的Km值1.170±0.088mM,这使得在嗜热链球菌体内,具有从α-酮戊二酸合成谷氨酸的代谢优势。控制化学合成培养基（Chemical-definedMedium,CDM）中的氨基酸组分,只添加维持S.thermophilusST2017生长必须的组氨酸、半胱氨酸、脯氨酸和甲硫氨酸。当乳糖浓度从0.5%提高到1.0%时,菌株ST2017的生物量从0.768±0.011上升到1.048±0.061,而突变株ST2017△codY的生物量则从0.822±0.017增长到1.381±0.023,说明在氮源限制性条件下,菌株生物量的增加来自于乳糖向氨基酸的转化。为了验证谷氨酸脱氢酶在乳糖向氨基酸转化中的作用,分别敲除野生菌株ST2017和突变株ST2017△codY的gdhA基因,得到gdhA突变株ST2017△gdhA和ST2017△codY△gdhA.测定其在0.5%乳糖的CDM中的生物量,分别为0.305±0.034和0.531±0.041,相对菌株ST2017和ST2017△codY显著下降；而乳糖浓度增加到1.0%,生物量只略微提高到0.328±0.009和0.581±0.104。表明在氮源限制条件下,嗜热链球菌通过谷氨酸脱氢酶将充足的碳源转化为氨基酸供菌体生长,codY基因的缺失提高了这种转化能力,当切断碳、氮代谢偶联节点时（缺失gdhA基因）,菌体生长受到严重抑制。在碳源限制条件下,将胰蛋白胨从0.25%提高到0.5%时,菌株ST2017和ST2017△codY的生物量分别提高到1.340±0.067和1.403±0.053,而谷氨酸脱氢酶基因突变株ST2017△gdhA和ST2017△codY△gdhA的生物量只能提高到1.188±0.004和1.164±0.044。说明谷氨酸脱氢酶的失活,使得氮源转化为碳源的能力也降低,最终限制了菌体生长。综合以上结果,可以看出嗜热链球菌全局转录调控因子CodY通过控制谷氨酸脱氢酶来影响碳、氮代谢的相互转化,以此调控菌体的生长。4、S.thermophilusST2017转录调控因子CodY在菌体抗胁迫中的作用嗜热链球菌是重要的乳制品发酵菌株,在发酵剂生产、存放及乳品发酵过程中常常遭受氧、温度、噬菌体等多种胁迫压力,从而导致菌株生长缓慢、发酵能力减弱,直接影响到产品的结构和风味。转录组学分析发现,二组分系统TCS0l在突变株ST2017△codY中显著上调,而TCS01参与压力响应蛋白的转录调控。为了检测因TCS01上调而引起的相关压力响应蛋白的转录变化,构建了二组分系统TCS01过表达菌株ST2017-TCS01exp,比较分析野生菌株ST2017、突变株ST2017△cody及ST2017-TCSOlexp中压力响应蛋白的差异表达。结果显示,嗜热链球菌CodY不但自身可以调控压力响应蛋白DnaK、GroESLSTND0113和STND0135的转录,而且通过控制二组分系统TCS01间接调控STND1346的转录。嗜热链球菌是兼性厌氧菌株,可以通过超氧化物歧化酶和NADH氧化酶降低胞内超氧自由基和分子氧浓度,同时利用谷胱甘肽/谷氧还原蛋白维持胞内氧化还原平衡。而谷胱甘肽合成酶基因（STND1346）的转录水平直接受二组分系统TCS01的抑制,CodY通过控制TCS01来间接调控STND1346的转录水平。比较S.thermophilusST2017和ST2017△codY的抗氧胁迫能力。在CDM培养基中添加谷氨酸和甘氨酸,经过H202处理后,ST2017的存活数是突变株的3倍以上；菌株ST2017细胞内谷胱甘肽浓度达到3.40±0.28μ/1,而突变株ST2017△codY为2.47±0.04gM/1；当不添加谷氨酸和甘氨酸时,菌株ST2017和ST2017△codY的存活数都显著下降,没有明显差异,并且胞内基本测不到谷胱甘肽含量。由此说明突变株ST2017△codY抗氧胁迫能力的下降是因为体内谷胱甘肽含量的减少。进一步研究发现,突变株ST2017△codY中,DnaK和GroESL操纵子的转录上调提高了菌体的高温胁迫耐受性；而STND0113、STND0135基因的转录下调则降低了菌体对渗透压胁迫的抗性。由此看出,全局转录调控因子CodY在嗜热链球菌对抗氧、温度、渗透压等胁迫压力中具有重要的调控功能。噬菌体污染是乳制品发酵过程中遇到的不可避免的安全隐患,为此发酵剂菌株自身编码多种防御机制抵抗噬菌体的感染,其中,CRISPR/Cas是最近发现的一种新的对抗外源DNA入侵的免疫机制。转录学分析发现,S.thermophilusST2017转录调控因子CodY对CRISPR1/Cas系统的cas5基因具有强烈的转录抑制作用。研究显示,二组分系统TCS01也抑制cas5的转录。比较S.thermophilusST2017和ST2017△codY中pSEC质粒的丢失率,结果发现无选择压力下传代100代,ST2017△codY的丢失率为8/768,远高于ST2017的1/788。将Cas5基因过表达,外源质粒在S.thermophilusST2017中的丢失率大大升高；同时提高了菌体S.thermophilusST1368对噬菌体的抗性。测序发现,在得到的质粒丢失菌株和噬菌体抗性菌株中,CRISPR序列中并没有新的间隔序列（Spacer）插入,这与已报道的二型CRISPR/Cas系统的作用机制不同。体外测定发现Cas5蛋白具有不依赖于crRNA和tracrRNA的核酸酶活性,能够随意切割λ-DNA和质粒。推测正是Cas5蛋白的核酸酶活性,使得外源质粒在嗜热链球菌CodY缺失突变株和Cas5过表达菌株中的丢失率提高,同时更容易产生对噬菌体的抗性,这是首次发现嗜热链球菌CRISPR/Cas系统中Cas5的非特异性切割引起的对外源DNA的免疫性。5、基于丙氨酸消旋酶的食品级表达系统的构建及应用随着遗传操作工具的完善,乳酸菌作为活疫苗载体受到重视。传统表达系统因带有抗生素抗性选择标记而限制其在实际生产中的应用。通过食品级表达体系表达抗原蛋白,使得乳酸菌成为新一代的活菌疫苗。本文以丙氨酸消旋酶基因（alr）为选择标记,分别构建了诱导型和组成型食品级分泌表达体系。首先通过温度敏感型质粒pG+host9和自杀性整合载体pRV300分别敲除了LactococcuslactisNZ9000和LactobacilluscaseiBL23的丙氨酸消旋酶基因,使宿主细胞自身成为D-Ala营养缺陷型,然后以丙氨酸消旋酶基因作为互补选择标记,实现功能回补。以绿色荧光蛋白作为报告基因,与诱导型nisA启动子融合表达,在乳酸乳球菌中构建了严谨的食品级诱导分泌表达体系,实现了猪圆环病毒衣壳蛋白高效的分泌表达,其最大分泌量达到1,059±99μg/1；且nisin诱导2小时,重组菌就能持续分泌衣壳蛋白,到10小时其分泌量依然达到940±101μg/1。同时还将卷曲乳杆菌（Lactobacilluscrispatus）S-层蛋白的启动子和核酸酶基因（nucB）融合,构建了适用于乳酸乳球菌和干酪乳杆菌的组成型食品级分泌表达载体,实现了肠毒性大肠杆菌（EnterotoxigenicEscherichiacoli,ETEC）鞭毛蛋白（FaeG）在乳酸乳球菌和干酪乳杆菌中分泌表达。上述构建的食品级表达体系,为乳酸菌活疫苗的开发奠定了基础。
【作者】陆文伟；
【导师】孔健；
【作者基本信息】山东大学，微生物学，2014，博士
【关键词】嗜热链球菌；全局转录调控因子CodY；代谢调控；抗胁迫机制；食品级表达体系；

【参考文献】
[1]郝洁.新经济条件下跨国公司价值链研究[D].武汉理工大学,技术经济及管理,2004,硕士.
[2]张二盼.对外汉语初级精读教材里的“在、正、正在、着、呢”研究[D].广东外语外贸大学,汉语国际教育,2013,硕士.
[3]杜晓梅.公司法人格外部逆向否认研究[D].中国政法大学,民商法学,2014,硕士.
[4]廖斌.广州乡村旅游发展对策研究[D].复旦大学,旅游管理,2012,硕士.
[5]章芹敏.SCM及其对我国高科技企业管理绩效评价发展的启示[D].中国地质大学,科学技术史,2003,硕士.
[6]王柏茵.基于一家医院机械通气患儿综合管理分析[D].河北医科大学,儿科学,2014,硕士.
[7]刘艾蕊.二十二冶机电公司综合计划管理研究[D].燕山大学,公共管理,2012,硕士.
[8]李保胜.不同种植体表面形貌对神经再生影响的初步研究[D].吉林大学,口腔医学,2013,硕士.
[9]李智博.第三人利益合同中受益第三人的履行请求权[D].吉林大学,法律,2014,硕士.
[10]秦世耀,李少林,王瑞明,陈晨,孙勇.大容量永磁同步风电机组系统谐振与抑制策略[J].电力系统自动化,2014,22:11-16.
[11]吴强乐.慢病毒介导的HER2-siRNA对人卵巢癌的基因治疗及干扰效应显像研究[D].苏州大学,影像医学与核医学,2013,硕士.
[12]张辰飞.我国主权信用的影响因素研究[D].北京交通大学,2014.
[13]郭玲.威斯特《弗吉尼亚人》中的牛仔形象研究[D].四川外国语大学,英语语言文学,2014,硕士.
[14]姚海琴.大众参与非专业类旅游节庆的影响因素研究[D].浙江大学,2006.
[15]杨少兵,吴命利.基于广义柯西分布的系统侧谐波阻抗估计[J].中国电机工程学报,2014,07:1159-1166.
[16]王东辉,韩光宇.安徽天长地区构造圈闭封堵特性的分析[J].石油地球物理勘探,1989,03:350-358+366.
[17]李冬菊.沙门氏菌Ⅲ型分泌系统效应蛋白SipA导致肠炎的分子机制[D].复旦大学,病原生物学,2011,博士.
[18]王月琴.基于改进近邻传播聚类的无线传感器网络分簇协议及优化[D].华东理工大学,控制科学与工程,2013,硕士.
[19]刘平昌唐娟.国内首条石墨烯薄膜生产线投产[N].中国化工报,2013-05-21002.
[20]袁舒.山东HL公司托管经营模式及托管绩效研究[D].山东大学,会计（专业学位）,2013,硕士.
[21]陈磊.“行家”“利家”语义流变考[D].渤海大学,美术学,2013,硕士.
[22]田书栋.20世纪以来景德镇陶瓷艺术风格变迁研究[D].景德镇陶瓷学院,美术,2014,硕士.
[23]朱英兰.尝试教学法在初中数学课堂中的实践研究[D].湖南师范大学,学科教学（专业学位）,2014,硕士.
[24]朱爱明.英汉朗读语篇调群切分模式对比研究[D].江苏科技大学,外国语言学及应用语言学,2014,硕士.
[25]张冠宇.高锰酸盐指数水质自动分析仪检定规程的研究[D].河北大学,控制工程,2014,硕士.
[26]孔德凌.《诗经》宴饮诗与周代礼乐文化的变迁[D].曲阜师范大学,中国古代文学,2004,硕士.
[27]王晓红,陈浩.1999-2006年我国各省市高校科研效率的实证研究——基于科技成果指标变化的对比分析[J].科研管理,2011,04:94-101+135.
[28]毛小林.岩石圈有效弹性厚度及其各向异性的计算和构造意义[D].南京大学,固体地球物理学（专业学位）,2013,硕士.
[29]陈砚雄.事务间量化关联规则挖掘的研究及应用[D].华中师范大学,2003.
[30]罗雪姣.小鼠真皮间充质干细胞局部注射对硬皮病小鼠模型局部皮损和转化生长因子-β1影响的实验研究[D].河北医科大学,皮肤病与性病学,2014,硕士.
[31]袁松.聚丙烯/纳米颗粒复合材料流变及剪切诱导结晶行为研究[D].天津大学,材料学,2013,硕士.
[32]余燕钟.铁路装卸向现代物流发展的分析与对策[D].西南交通大学,交通运输工程,2013,硕士.
[33]程逸宇.应用cDNA-AFLP技术分析铜胁迫下水稻根系基因差异表达[D].南京农业大学,植物营养,2012,硕士.
[34]刘翠.妊娠期肝内胆汁淤积症围产儿不良结局的研究[D].遵义医学院,儿科学（专业学位）,2014,硕士.
[35]王建源,冯晓云,薛德强,刘欣,谢爱红.GIS在泰安市板栗农业气候区划中的应用[J].中国农业资源与区划.2003(05)
[36]吴晨昊.六年级学生数感的发展[D].华东师范大学,课程与教学论,2013,硕士.
[37]顾容,何益明,金丽燕.基于Kinect骨骼追踪技术的智能教学控制系统的设计与实现[J].计算机应用与软件,2013,06:276-280.
[38]吴海龙,林晓青,严密,陈彤,陆胜勇,李晓东,严建华.基于硫基循环抑制技术的危险废物焚烧炉二噁英排放的控制[J].化工学报,2014,11:4593-4598.
[39]艾爱.GN公司在新能源环境下的战略管理研究[D].华中科技大学,工商管理,2013,硕士.
[40]朱宏伟.石墨烯:单原子层二维碳晶体——2010年诺贝尔物理学奖简介[J].自然杂志,2010,06:326-331+308.
[41]程渤芮.高新技术企业员工关系管理影响因素的实证研究[D].成都理工大学,企业管理,2013,硕士.
[42]李连杰.基于柔性配体构筑的金属有机骨架的合成、结构和荧光性质研究[D].东北师范大学,物理化学,2012,硕士.
[43]子华明.大学生吸烟行为影响因素分析与干预对策研究[D].华东师范大学,应用心理学,2013,硕士.
[44]蓝利纺.左右Robertshaw支气管导管插管后咽喉损伤并发症的比较[D].广西医科大学,麻醉学（专业学位）,2013,硕士.
[45]耿冠世.基于地震台阵数据的场地效应研究[D].中国地震局地球物理研究所,固体地球物理学,2014,博士.
[46]王敏.从费唐报告看近代上海崛起之路[J].学术月刊,2014,06:127-137.
[47]刘少林,刘有山,汪文帅,李小凡,张美根.最优化广义离散Shannon奇异核褶积微分算子地震波场模拟[J].石油地球物理勘探,2013,03:410-416+506+328-329.
[48]马昌科.模拟人体三维环境下微观沟槽结构对细胞取向性生长的调控[D].苏州大学,骨外科（专业学位）,2014,硕士.
[49]王倩倩.反共振振动机虚拟样机及性能预测的研究[D].东北大学,机械电子工程,2010,硕士.
[50]朱琪.LC调谐磁悬浮轴承研究[D].南京航空航天大学,机械设计及理论,2013,硕士.