轉(zhuǎn)錄水平解析昆蟲(chóng)病原線蟲(chóng)共生菌菌株新型殺菌蛋白PPIA-L20基因表達(dá)調(diào)控——結(jié)果與分析
2 結(jié)果與分析
2.1 不同發(fā)酵時(shí)間ppia-l20基因的表達(dá)變化
分別取4～84 h不同發(fā)酵時(shí)間點(diǎn)的NN6菌株的菌體細(xì)胞，提取Total RNA進(jìn)行6次實(shí)時(shí)熒光定量PCR測(cè)定，結(jié)果(圖1)表明ppia-l20相對(duì)表達(dá)水平極不穩(wěn)定且趨勢(shì)多樣，即使是使用同一內(nèi)參基因16S rRNA。選取16S rRNA和ropB兩種內(nèi)參基因，共同定量得到相似基因表達(dá)趨勢(shì)的8和12 h 2個(gè)時(shí)間點(diǎn)的菌株進(jìn)行后續(xù)的轉(zhuǎn)錄組測(cè)序(圖1-F)。
圖1 不同發(fā)酵時(shí)間ppia-l20基因的表達(dá)水平
2.2 Xenorhabdus bovienii NN6發(fā)酵時(shí)間與菌體總數(shù)、活菌數(shù)、基因相對(duì)表達(dá)水平相關(guān)性分析
從圖2可知，隨著發(fā)酵時(shí)間的延長(zhǎng)，菌體量逐漸增加。從表3可見(jiàn)，對(duì)數(shù)期活菌數(shù)與菌體總數(shù)(D600表示)與發(fā)酵時(shí)間成正相關(guān)，說(shuō)明隨著發(fā)酵時(shí)間的延長(zhǎng)，對(duì)數(shù)期的菌體大量繁殖，但與ppia-l20基因表達(dá)水平無(wú)相關(guān)性。
圖2 發(fā)酵過(guò)程中菌體量變化
表3 對(duì)數(shù)期活菌數(shù)、菌體總數(shù)、ppia-l20基因相對(duì)表達(dá)情況相關(guān)性分析
2.3 RNA-Seq測(cè)序數(shù)據(jù)質(zhì)量分析
對(duì)8和12 h 2組不同發(fā)酵時(shí)間菌體細(xì)胞共6個(gè)樣本進(jìn)行高通量測(cè)序，并對(duì)測(cè)序質(zhì)量進(jìn)行統(tǒng)計(jì)分析。結(jié)果(表4)顯示:原始測(cè)序序列(raw reads)去除帶接頭序列、空序列及低質(zhì)量測(cè)序序列后分別得到15 352 102、14 498 440、18 087 468、15 466 120、15 821 766、16 257 026條純凈序列(clean reads)，堿基GC含量平均為48.86%，低質(zhì)量堿基測(cè)序序列(<Q20)占原始序列的比例很低，而高質(zhì)量堿基序列(≥Q20)平均比例達(dá)到98.18%，說(shuō)明測(cè)序數(shù)據(jù)的堿基組成情況和質(zhì)量良好。將過(guò)濾后的測(cè)序序列比對(duì)到參考基因組上，發(fā)現(xiàn)測(cè)序序列均勻覆蓋在參考基因組上，能定位到基因組上測(cè)序序列的平均比例達(dá)87.89%。這些結(jié)果表明此次測(cè)序序列的整體質(zhì)量良好，測(cè)序數(shù)據(jù)隨機(jī)性良好，可用于后續(xù)的分析。
表4 測(cè)序數(shù)據(jù)質(zhì)量統(tǒng)計(jì)
2.4 Xenorhabdus bovienii NN6兩組間差異表達(dá)基因(DEG)分析
以發(fā)酵8 h菌體中的基因表達(dá)水平為參照，利用DESeq2(V1.6.3)分析發(fā)酵12 h菌體中差異表達(dá)的基因，設(shè)置篩選閾值 P<0.05，log2(Fold change)>1，構(gòu)建火山圖(圖3-A)。RNA-seq共檢測(cè)到2 716個(gè)基因，其中466個(gè)DEG，基因組占比17.2%，其中193個(gè)DEG上調(diào)，273個(gè)DEG下調(diào)。DEG的熱圖聚類(lèi)分析表明重復(fù)組內(nèi)基因表達(dá)趨勢(shì)基本相同，組間差異顯著(圖3-B)。
圖3 Xenorhabdus bovienii NN6兩組間差異表達(dá)基因
2.5 RT-qPCR驗(yàn)證轉(zhuǎn)錄組數(shù)據(jù)
隨機(jī)選取8個(gè)DEG(包括4個(gè)上調(diào)基因和4個(gè)下調(diào)基因)和2個(gè)表達(dá)差異不顯著的基因，Real-time PCR驗(yàn)證結(jié)果(圖4)顯示基因表達(dá)趨勢(shì)與測(cè)序數(shù)據(jù)結(jié)果一致，說(shuō)明RNA-seq轉(zhuǎn)錄組測(cè)序結(jié)果數(shù)據(jù)可靠。
圖4 RNA-seq數(shù)據(jù)與RT-qPCR比對(duì)
2.6 Xenorhabdus bovienii NN6兩組間GO富集分析
進(jìn)一步對(duì)8和12 h 2組中的DEG進(jìn)行GO注釋分析，結(jié)果(圖5)顯示有442個(gè)DEG注釋到生物學(xué)過(guò)程、細(xì)胞組分和分子功能大類(lèi)的31個(gè)功能小類(lèi)中。分子功能分類(lèi)中含有最多的差異基因，分子功能中差異基因主要?dú)w屬于雜環(huán)化合物結(jié)合、轉(zhuǎn)錄調(diào)節(jié)因子活性、氮化合物代謝和轉(zhuǎn)移酶活性等。生物學(xué)功能和細(xì)胞組分分類(lèi)中的差異基因數(shù)量少于分子功能分類(lèi)，生物學(xué)過(guò)程中差異基因主要富集在代謝過(guò)程和應(yīng)激反應(yīng)。值得注意的是，細(xì)胞組分分類(lèi)中僅富集到2個(gè)功能小類(lèi)，主要?dú)w屬于細(xì)胞和膜部分，與細(xì)胞膜的功能相關(guān)。這表明Xenorhabdus bovienii NN6菌株可能在發(fā)酵過(guò)程中代謝活躍，ppia-l20表達(dá)量變化可能與轉(zhuǎn)移酶活性和膜蛋白轉(zhuǎn)運(yùn)過(guò)程有關(guān)。
為了探究DEG在ppia-l20基因表達(dá)變化中的潛在功能，對(duì)所有DEG進(jìn)行GO富集分析，分別統(tǒng)計(jì)生物過(guò)程(BP)、細(xì)胞成分(CC)和分子功能(MF)3個(gè)分類(lèi)富集程度前10的GO條目，圖6結(jié)果表明在這30個(gè)GO條目中，差異表達(dá)基因在生物過(guò)程中主要富集于細(xì)胞膜(membrane)、定位(localization)、水解酶(hydrolase activity)、膜整體成分(integral component of membrane)、膜內(nèi)成分(intrinsic component of membrane)、轉(zhuǎn)運(yùn)活性(transporter activity)。其中細(xì)胞膜(membrane)的富集率最高為77.5%，表明在發(fā)酵過(guò)程中Xenorhabdus bovienii NN6菌株與膜相關(guān)GO條目中的77.5%的基因差異表達(dá)，提示發(fā)酵過(guò)程中ppia-l20基因表達(dá)升高可能與細(xì)胞膜功能相關(guān)蛋白的表達(dá)密切相關(guān)。
圖5 Xenorhabdus bovienii NN6兩組間差異表達(dá)基因的GO注釋分類(lèi)圖
圖6 Xenorhabdus bovienii NN6組間差異表達(dá)基因各分類(lèi)前10的GO富集分析
2.7 Xenorhabdus bovienii NN6兩組間DEG的KEGG富集分析
KEGG注釋結(jié)果(圖7)顯示，有191個(gè)DEG注釋到了62個(gè)KEGG代謝通路。基于KEGG數(shù)據(jù)庫(kù)(https://www.kegg.jp/kegg/pathway.html)進(jìn)行通路分析，其中13條通路的可信度顯著(P<0.05)。上調(diào)基因主要富集在次級(jí)代謝產(chǎn)物合成、氨基酸合成、雙組分調(diào)節(jié)系統(tǒng)相關(guān)通路; 下調(diào)基因主要富集在次級(jí)代謝產(chǎn)物合成相關(guān)通路、鞭毛組裝與磷酸轉(zhuǎn)移酶系統(tǒng)相關(guān)通路。
圖7 Xenorhabdus bovienii NN6差異表達(dá)基因前20的KEGG富集分析
2.8 KEGG關(guān)鍵通路中代表性DEG分析
從表5可知，參與雙組分調(diào)節(jié)系統(tǒng)共有14個(gè)DEG，其中12個(gè)表達(dá)上調(diào)2個(gè)表達(dá)下調(diào)。ompR基因編碼的轉(zhuǎn)錄調(diào)控蛋白OmpR是雙組分調(diào)節(jié)系統(tǒng)的關(guān)鍵蛋白，qseC基因編碼雙組分傳感器激酶N端結(jié)構(gòu)。目的基因ppia-l20位于陽(yáng)離子抗菌肽抵制通路中，位于同一條通路上的3個(gè)基因中，arnF和eamA編碼的蛋白屬于DUF6轉(zhuǎn)錄因子家族，其表達(dá)上調(diào)，另一個(gè)基因nlpE負(fù)責(zé)編碼外膜脂蛋白，參與銅離子的攝取和黏附功能，其表達(dá)下調(diào)。
表5 關(guān)鍵通路中代表性差異基因表達(dá)情況
2.9 PPIA-L20相關(guān)蛋白互作網(wǎng)絡(luò)分析
通過(guò)STRING對(duì)PPIA-L20蛋白進(jìn)行蛋白互作網(wǎng)絡(luò)構(gòu)建，以獲取其潛在蛋白的互作關(guān)系，如圖8所示，PPIA-L20可能與AMPs結(jié)合蛋白TycC、GrsB、ATP酶HtpG、酮酰合成酶N端PpsE蛋白等蛋白相互作用，AMPs結(jié)合蛋白TycC、GrsB定位于細(xì)胞外膜，用于識(shí)別外界AMPs結(jié)合信號(hào)，ATP酶HtpG和酮酰合成酶N端PpsE蛋白參與細(xì)胞的能量代謝過(guò)程。
GO富集分析表明，ppia-l20基因變化與細(xì)胞膜的變化、轉(zhuǎn)移酶活性、代謝過(guò)程相關(guān)，KEGG通路分析預(yù)測(cè)，PPIA-L20蛋白位于陽(yáng)離子抗菌肽抵抗通路，與雙組分調(diào)節(jié)系統(tǒng)、磷酸轉(zhuǎn)移酶系統(tǒng)活動(dòng)相關(guān); 與PPIA-L20蛋白位于共同通路(陽(yáng)離子抗菌肽抵抗通路)的nlpE基因表達(dá)的外膜脂蛋白NlpE屬于DUF6轉(zhuǎn)錄因子家族，NlpE蛋白與雙組分調(diào)節(jié)系統(tǒng)中的cpxR基因、cpxA基因表達(dá)的蛋白具有互作關(guān)系。在腸桿菌科中，控制AMPs耐藥性的基因通常受雙組分信號(hào)通路以及磷酸轉(zhuǎn)移系統(tǒng)的控制，NlpE蛋白同時(shí)存在于PPIA-L20蛋白共同通路(陽(yáng)離子抗菌肽抵抗通路)的上游，可能起到調(diào)控ppia-l20基因表達(dá)的作用，該蛋白相互作用網(wǎng)絡(luò)的構(gòu)建為調(diào)控機(jī)制的推測(cè)提供線索。
圖8 蛋白互作的網(wǎng)絡(luò)預(yù)測(cè)
相關(guān)新聞推薦
1、巨大芽孢桿菌BM24生長(zhǎng)曲線、耐藥性、益生潛力評(píng)估及藥敏試驗(yàn)（一）
2、國(guó)內(nèi)知名品牌常溫酸奶中的乳酸菌菌落總數(shù)測(cè)定
3、李斯特菌ST87在25°C不同營(yíng)養(yǎng)條件下牛肉和蔬菜汁中的生長(zhǎng)曲線
4、一文了解細(xì)菌生長(zhǎng)曲線的定義、階段、影響因素以及應(yīng)用
5、優(yōu)化培養(yǎng)基配方減少無(wú)機(jī)硒對(duì)粗毛纖孔菌的抑制作用，同時(shí)提高富硒菌絲的產(chǎn)量（三）