李 婷，林玉暖，林茂娟，牛曉磊，李春霞，陳銀華，陶 均

(海南大學(xué) 海南省熱帶生物資源可持續(xù)利用重點(diǎn)實(shí)驗(yàn)室/熱帶作物學(xué)院，海南 ?？?570228)

雙組分信號傳導(dǎo)系統(tǒng)(Two-component signal transduction systems，TCSTSs)是細(xì)菌感應(yīng)內(nèi)外界信號最主要的方式[1-2]，通常由感應(yīng)信號的受體組氨酸激酶(histidine kinase, HK)和執(zhí)行功能的響應(yīng)調(diào)節(jié)因子(response regulator, RR)組成[1-2].幾乎所有的細(xì)菌都編碼TCSTSs，只是數(shù)目上有差異，少則幾個(gè)，多則上百個(gè)，但都調(diào)控著細(xì)菌絕大多數(shù)的生理過程[3].此外，一個(gè)HK可能對應(yīng)多個(gè)RR；同樣一個(gè)RR可能被多個(gè)HK磷酸化，TCSTSs之間存在著復(fù)雜的信號交叉[4-5].調(diào)控相同或相似代謝途徑的TCSTSs在進(jìn)化過程中可能趨向于融合或產(chǎn)生級聯(lián)作用，進(jìn)而更加精確地調(diào)控細(xì)胞行為，以適應(yīng)不同的環(huán)境變化[3].

RavA/RavR是在水稻白葉枯病菌(Xanthomonasoryzaepv.oryzae,Xoo)和野油菜黃單胞菌(Xanthomonascampestrispv.campestris,Xcc)中發(fā)現(xiàn)的，參與細(xì)胞第二信使環(huán)雙鳥苷酸Cyclic di-GMP(c-di-GMP)代謝轉(zhuǎn)換，調(diào)控細(xì)菌的群體感應(yīng)、運(yùn)動(dòng)能力及胞外多糖，是影響病原菌致病力的重要雙組分信號傳導(dǎo)系統(tǒng)[6-7].RaxH/RaxR是Xoo中調(diào)控細(xì)菌效應(yīng)因子AvrXa21加工分泌的主要調(diào)節(jié)因子，其可能感應(yīng)群體信號，調(diào)控Ⅰ型分泌系統(tǒng)raxAB及酪氨酸硫化酶raxST的表達(dá)，影響病原菌與植物的互作[8-11].在Xcc中，RavA/RavR可以調(diào)控Ⅰ型分泌系統(tǒng)raxAB的表達(dá)[6]，暗示RaxH/RaxR和RavA/RavR可能共同調(diào)控Ⅰ型分泌系統(tǒng)，或處于信號傳導(dǎo)通路的不同位置.但RaxH/RaxR和RavA/RavR是否在進(jìn)化上有關(guān)聯(lián)性以及兩種之間在功能上是否相關(guān)還未有報(bào)道.本研究首先分析了RaxH與RaxR以及RavA與RavR的共進(jìn)化關(guān)系，然后分析了這兩個(gè)雙組分系統(tǒng)在進(jìn)化上的相關(guān)性以及它們在調(diào)控上的關(guān)聯(lián)性，發(fā)現(xiàn)了這兩個(gè)雙組分系統(tǒng)高度協(xié)同進(jìn)化，RaxH/RaxR在表達(dá)水平調(diào)控RavA/RavR，形成一個(gè)級聯(lián)的調(diào)控通路.

1 材料與方法

1.1 材 料RaxH/RaxR與RavA/RavR蛋白序列下載于NCBI(Http://www.ncbi.nlm.nih.gov).XooPXO99由本實(shí)驗(yàn)室保存.化學(xué)試劑購自廣州化學(xué)試劑公司(廣州)和Sigma-Aldrich(上海).PCR相關(guān)試劑購自全式金生物科技有限公司(北京).引物由華大基因生物科技有限公司(北京)合成，序列見表1.限制性內(nèi)切酶購自NEB China(北京).

表1 本研究所用的引物

1.2 生物信息學(xué)分析序列比對及進(jìn)化樹分析采用MEGA6[12]，進(jìn)化分析采用最大似然法.蛋白質(zhì)二級結(jié)構(gòu)分析采用SMART法(http://smart.embl-heidelberg.de/)[13-14].共進(jìn)化分析采用Paralog Matching法[15]，采用KEGG數(shù)據(jù)庫[16]和String在線分析法(https://string-db.org/)分析共存在[17].

1.3 基因敲除采用同源重組兩步交換法構(gòu)建缺失突變體[6,18,19].以菌株P(guān)XO99的基因組為模板，擴(kuò)增待敲除基因的左右兩端各400～500 bp的片段，構(gòu)建到自殺質(zhì)粒pK18mobSacB上.將構(gòu)建正確的質(zhì)粒載體用電擊轉(zhuǎn)化方法轉(zhuǎn)入PXO99中，涂布在含有卡那霉素的PGA(即PGK)平板上，28℃培養(yǎng)2-3天.將PGK平板上生長的菌落，選取4～8個(gè)單克隆同時(shí)接種到PGK平板和含有15%蔗糖的PSA(即PSSU)平板上.選取在PGK平板上生長但在PSSU平板上不生長的2個(gè)單克隆，在PGA液體培養(yǎng)基中，28℃振蕩培養(yǎng)3天.取適量菌液涂布在PSSU平板上，28℃培養(yǎng)3天.選取在PSSU平板上生長的20個(gè)黃色單菌落，同時(shí)接種于PGK和PSSU平板上.選取在PSSU平板上生長但在PGK上不生長的菌體，進(jìn)行菌落PCR驗(yàn)證，獲得基因缺失突變體.

1.4 RNA提取及qPCR分析PXO99，ΔravA，ΔravR及ΔraxR接種于PSA培養(yǎng)基(1%牛肉膏，1%酵母提取物，1%蔗糖)中，1天后按1∶100轉(zhuǎn)接到新鮮的PSA中，待OD600≈0.8時(shí)，收集菌體，液氮速凍后進(jìn)行RNA提取.RNA提取和qPCR分析參照已發(fā)表文獻(xiàn)[18-19].PCR引物見表1.

1.5 數(shù)據(jù)分析及作圖采用GraphPad Prism 5作圖，圖像后期處理采用Adobe Photoshop 8.0軟件.

2 結(jié)果與分析

2.1 RavA/RavR的共進(jìn)化分析以野油菜黃單胞菌(Xcc)ravA在基因組的位置為標(biāo)準(zhǔn)，分析其附近的基因在其它菌株中是否也與ravA共同存在.首先通過同源搜索尋找同源蛋白，并定位每個(gè)相應(yīng)基因在其基因組中的位置.然后將獲得的物種數(shù)目作比較，分析這些基因與ravA共同存在的比率.結(jié)果顯示，ravR與ravA共同存在于同一基因組的概率是97%，而Xcc基因組上ravA相鄰的其它基因在其它物種中同時(shí)存在的概率都小于70%(圖1A)，表明RavR與RavA在通常情況下組成一對雙組分系統(tǒng).由于RavS也被認(rèn)為可磷酸化RavR[7]，但其與RavR共同進(jìn)化率低，因此RavS可能只在某些特定細(xì)菌中與RavR組成雙組分系統(tǒng)，不具有廣泛性，可能是RavR在進(jìn)化過程中獲得的新調(diào)控模式.同源性分析發(fā)現(xiàn)RavA和RavR都可以分為兩大類(圖1B和1C)，黃單胞菌中的RavA和RavR都在一大類中，且序列非常相似.但是在大多數(shù)細(xì)菌中，RavA和RavR都是相當(dāng)保守的，預(yù)示其調(diào)控功能也是保守的.

2.2 RaxH/RaxR的共進(jìn)化分析按照上述RavA/RavR共進(jìn)化分析方法，以Xcc中raxH在基因組的位置為標(biāo)準(zhǔn)，分析其附近的基因在其它菌中是否也與raxH共同存在.結(jié)果顯示，raxR與raxH共同存在于同一基因組的概率是96%，而Xcc基因組上raxH相鄰的基因在其它物種中同時(shí)存在的概率均小于50%(圖2A)，同樣說明RaxH與RaxR在通常情況下組成一對雙組分系統(tǒng).raxR和raxH組成一個(gè)操縱子，其基因組上下游基因的共進(jìn)化概率低，說明raxR/raxH的功能可能與其相鄰的基因沒有多大的關(guān)系.同RavA和RavR一樣，RaxH/RaxR在同源性上也可以分為兩大類(圖2B和2C)，黃單胞菌編碼蛋白都在一大類中.RaxH/RaxR的序列在不同菌中也是高度保守的，其調(diào)控功能也應(yīng)該是保守的.

圖1 RavA和RavR共進(jìn)化(A)和聚類分析(B和C)

圖2 RaxH和RaxR共進(jìn)化(A)和聚類分析(B和C)

2.3 RaxH/RaxR與RavA/RavR之間的共進(jìn)化分析RaxH/RaxR與RavA/RavR在細(xì)菌中都高度保守，各自之間的共進(jìn)化關(guān)系明顯(圖1和圖2)，而且在水稻黃單胞菌中，RaxH/RaxR可能感應(yīng)和群體感應(yīng)信號，而RavA/RavR則參與細(xì)胞第二信使c-di-GMP的代謝.c-di-GMP是調(diào)控細(xì)菌群體行為的主要信使[6, 8, 10, 20].因此，RaxH/RaxR與RavA/RavR可能在功能或調(diào)控上存在關(guān)聯(lián)，在細(xì)菌進(jìn)化過程中共同進(jìn)化.為此，我們分析了這兩個(gè)系統(tǒng)間的共進(jìn)化關(guān)系.如圖3A所示，兩者的共進(jìn)化系數(shù)接近80%，說明兩者之間關(guān)聯(lián)性很強(qiáng).在選擇的67個(gè)菌株中，56個(gè)編碼RaxH/RaxR，50個(gè)編碼RavA/RavR，編碼RavA/RavR的菌株均編碼RaxH/RaxR(表2)，進(jìn)一步表明兩者之間關(guān)聯(lián)性很強(qiáng).另外，幾乎所有的黃單胞菌都編碼RaxH/RaxR和RavA/RavR(圖3B)，說明這兩個(gè)系統(tǒng)在黃單胞菌的生存過程中起到非常重要的作用.從這些基因在基因組中的位置及其轉(zhuǎn)錄方向看，它們的保守性也很強(qiáng)：ravA和ravR轉(zhuǎn)錄方向相對，擁有自身的啟動(dòng)子，而raxR和raxH組成一個(gè)操縱子，一起轉(zhuǎn)錄(表2，最后行).這些結(jié)果說明RaxH/RaxR和RavA/RavR存在著共進(jìn)化的關(guān)系，在功能或調(diào)控上有相關(guān)性.

表2 部分細(xì)菌同時(shí)編碼RaxHR和RavAR

2.4 RaxH/RaxR在轉(zhuǎn)錄水平調(diào)控RavA/RavR以上分析表明RaxH/RaxR與RavA/RavR的共進(jìn)化關(guān)系明顯(圖3).有研究表明，在水稻黃單胞菌中RaxH/RaxR可能感應(yīng)和群體感應(yīng)信號，在介導(dǎo)細(xì)菌與植物互作過程中起到非常重要的作用，RavR可分解c-di-GMP，調(diào)控細(xì)菌群體行為，且兩者突變都影響細(xì)菌的致病能力[6,8,10,20].因此，RaxH/RaxR與RavA/RavR可能在調(diào)控上存在相關(guān)性.為此，我們首先構(gòu)建了這些基因的突變體，然后比較了突變體和野生型菌株中這些基因的表達(dá)情況.結(jié)果顯示，RaxH/RaxR可負(fù)調(diào)控RavA/RavR的表達(dá)(圖4)，表明兩者之間確實(shí)存在著相關(guān)性.

3 討 論

單細(xì)胞特點(diǎn)決定細(xì)菌需要對外界環(huán)境做出快速響應(yīng).TCSTSs是原核生物最重要的感應(yīng)響應(yīng)系統(tǒng)，其數(shù)目的多少是其應(yīng)變能力的衡量標(biāo)準(zhǔn)[21].不同細(xì)菌編碼數(shù)目不同的TCSTSs，如Xanthomonas屬共編碼92-121個(gè)TCSTSs蛋白，調(diào)控細(xì)胞分裂、運(yùn)動(dòng)能力、脅迫響應(yīng)、胞外多糖、生物膜等多個(gè)生理生化過程，影響細(xì)菌的生存能力[22-23].有些TCSTS共同調(diào)控一個(gè)代謝過程，而有的HK可能調(diào)控多個(gè)RR的活性[21].這些TCSTSs可能在進(jìn)化過程中被一起保留了下來.RaxH同RaxR的高度共進(jìn)化特性(圖2)表明這對雙組分系統(tǒng)對多數(shù)細(xì)菌的生存是必要的.同樣RavA與RavR的共進(jìn)化特征也很明顯(圖1)，說明它們組成的TCSTS對細(xì)菌來說同樣非常重要.在Xoo中，RaxH/RaxR與RavA/RavR調(diào)控相似的信號途徑，影響細(xì)菌與宿主的互作[20]，暗示它們之間可能在進(jìn)化上有協(xié)同進(jìn)化的特性.通過對不同細(xì)菌中這兩個(gè)TCSTSs共同進(jìn)化特征分析發(fā)現(xiàn)兩者共同存在于一個(gè)細(xì)菌基因組的概率很高(圖3和表2)，證實(shí)了它們之間的關(guān)聯(lián)性.同時(shí)raxH/raxR突變會(huì)導(dǎo)致ravA/ravR的表達(dá)上升(圖4)，說明這兩個(gè)雙組分系統(tǒng)確實(shí)存在著調(diào)控上的關(guān)聯(lián)性.RaxH/RaxR可能位于整個(gè)信號傳導(dǎo)的上游，而RavA/RavR處于中游，其調(diào)控的c-di-GMP作為第二信使再引起下游基因的表達(dá)變化，從而影響諸如生物膜、運(yùn)動(dòng)性及AvrXa21的加工及分泌等生理過程，調(diào)控病原菌的致病能力.我們后面的研究將關(guān)注這兩個(gè)TCSTSs之間的調(diào)控網(wǎng)絡(luò)，并進(jìn)一步分析RaxH/RaxR調(diào)控RavA/RavR表達(dá)的分子機(jī)制，為最終解析兩者之間的關(guān)系打下基礎(chǔ).