2.2鼠李糖乳桿菌Probio-M9連續(xù)傳代過程中比較基因組學(xué)分析
2.2.1基因組信息本研究中不同代時(shí)菌株組裝的基因組信息如表3所示，鼠李糖乳桿菌Probio-M9原始菌株基因組大小3.00 Mb，GC含量46.76%。傳代過程中菌株的平均基因組大小為（3.00±0.002）Mb，GC含量（46.76±0.007）%。
表3 5個(gè)不同代時(shí)15株菌基因組信息
2.2.2系統(tǒng)發(fā)育分析
本研究中利用系統(tǒng)發(fā)育樹來描述傳代過程中菌株的進(jìn)化關(guān)系，進(jìn)而判斷菌株在連續(xù)培養(yǎng)過程中的穩(wěn)定性。以鼠李糖乳桿菌CNCM-I-3698作為外群，使無根樹轉(zhuǎn)為有根樹，鼠李糖乳桿菌Probio-M9原始基因組為參考基因組，使用Soapsnp方法得到的15個(gè)菌株的SNP位點(diǎn)連成序列，通過鄰接法（NJ法）構(gòu)建系統(tǒng)發(fā)育樹。
如圖4所示，系統(tǒng)發(fā)育樹共分為兩個(gè)分支，5個(gè)代時(shí)的15株菌聚在同一系統(tǒng)發(fā)育分支，親緣關(guān)系極近，且明顯與外群基因組處于不同分支。
圖4基于SNP位點(diǎn)構(gòu)建系統(tǒng)發(fā)育樹
2.2.3平均核苷酸一致性分析
平均核苷酸一致性（ANI）是在核苷酸的水平上，用以比較兩個(gè)基因組之間親緣關(guān)系的常用指標(biāo)，其特點(diǎn)是對(duì)于近緣種群的區(qū)分度更高。在比較基因組學(xué)分析中，常用ANI評(píng)估基因組間的多態(tài)性和相似性，當(dāng)ANI大于95%時(shí)，即表示為同一個(gè)種。通過計(jì)算鼠李糖乳桿菌Probio-M9原始基因組和15株不同代時(shí)菌株基因組間的ANI值，并繪制成熱圖。
由圖5可知，鼠李糖乳桿菌Probio-M9原始基因組和15株不同代時(shí)菌株基因組，兩兩菌株間的ANI值均在99.00%，菌株之間表現(xiàn)出極高的相似性，且高于種間鑒定標(biāo)準(zhǔn)，可準(zhǔn)確鑒定為同一亞種，表明連續(xù)培養(yǎng)100代過程中菌株均屬同一亞種即鼠李糖乳桿菌Probio-M9。
圖5基于5個(gè)不同代時(shí)基因組序列繪制ANI熱圖
2.2.4共線性分析
共線性可用來描述同一染色體上基因的位置關(guān)系，對(duì)不同的基因組間通過基因排列順序?qū)ζ涔餐嫦冗M(jìn)行探索，基因組相關(guān)性越高則其中同源序列一致性越高。本研究選取鼠李糖乳桿菌Probio-M9原始基因組序列作為參考，利用Mauve軟件對(duì)15株不同代時(shí)基因組進(jìn)行共線性分析。
結(jié)果如圖6所示，15株不同代時(shí)培養(yǎng)的菌株與鼠李糖乳桿菌Probio-M9原始基因組的共線性較好，具有極高的相似性，每個(gè)基因組中均未檢測(cè)到基因片段的倒位。與鼠李糖乳桿菌Probio-M9比較，A50-3和A75-2有大片段缺失，約為15 kb，主要包含整合酶家族蛋白、組氨酸激酶、核糖核酸酶、假設(shè)蛋白、tuf、cobB_1和ATP結(jié)合酶蛋白等基因。其余菌株均有小片段缺失，且基因缺失位置隨機(jī)。
圖6基于15株不同代時(shí)序列組裝結(jié)果與鼠李糖乳桿菌Probio-M9原始基因組的共線性分析
2017年，Zhang等發(fā)現(xiàn)一種由asp23基因編碼的堿性休克蛋白參與了干酪乳桿菌Zhang對(duì)慶大霉素的適應(yīng)性，從蛋白質(zhì)組學(xué)對(duì)其適應(yīng)性機(jī)制進(jìn)行探究，通過比較含有慶大霉素的培養(yǎng)基中生長(zhǎng)的asp23基因缺失突變體及其親本菌株的蛋白質(zhì)組，發(fā)現(xiàn)與asp23突變株相比，一些膜相關(guān)蛋白在親本菌株中顯著上調(diào)，表明asp23編碼的堿性休克蛋白能夠通過調(diào)節(jié)細(xì)胞膜蛋白來適應(yīng)慶大霉素的環(huán)境脅迫。本研究中15株不同代時(shí)培養(yǎng)的菌株均有一定的基因片段缺失現(xiàn)象發(fā)生，因此對(duì)于本研究中缺失片段所表達(dá)的基因功能還有待進(jìn)一步研究。
2.2.5 SNP突變位點(diǎn)以鼠李糖乳桿菌Probio-M9為參考菌株，利用Soapsnp軟件對(duì)15株不同代時(shí)鼠李糖乳桿菌Probio-M9菌株進(jìn)行SNP位點(diǎn)識(shí)別，同時(shí)使用Samtools軟件對(duì)測(cè)序數(shù)據(jù)檢測(cè)到的突變位點(diǎn)進(jìn)行過濾，剔除假陽性結(jié)果，結(jié)果見表4。
表4 15株不同代時(shí)菌株提取的SNP突變位點(diǎn)及注釋信息
以鼠李糖乳桿菌Probio-M9原始菌株基因組為參照，鼠李糖乳桿菌Probio-M9在連續(xù)培養(yǎng)100代時(shí)的不同代時(shí)菌株中共發(fā)現(xiàn)16個(gè)SNP位點(diǎn)，SNP位點(diǎn)分布及功能如表4所示，包括6個(gè)同義突變和4個(gè)非同義突變，6個(gè)位于非編碼區(qū)。非同義突變主要位于編碼酪氨酸受體激酶YwqD（Tyrosine-protein kinase）和假設(shè)蛋白（Hypothetical protein）的基因上。根據(jù)突變位點(diǎn)對(duì)其穩(wěn)定性進(jìn)行分析發(fā)現(xiàn)，連續(xù)培養(yǎng)100代過程中高頻突變數(shù)量較少，檢測(cè)到2個(gè)SNP，且首次發(fā)生在第0代基因組與參考基因組之間，一旦建立將穩(wěn)定遺傳。菌株低頻突變檢測(cè)到的個(gè)數(shù)均小于21。研究表明，SNP個(gè)數(shù)小于21時(shí)，可將其鑒定為同一物種。不同代時(shí)15株菌所識(shí)別到的SNP位點(diǎn)均小于21個(gè)，表明15株不同代時(shí)菌株可鑒定為同一物種。
5個(gè)代時(shí)菌株SNP位點(diǎn)隨機(jī)產(chǎn)生，從0代到25代、25代到50代、50代到75代、75代到100代的菌株中，均未發(fā)現(xiàn)穩(wěn)定存在的突變位點(diǎn)。表明鼠李糖乳桿菌Probio-M9在連續(xù)傳代100代時(shí)過程中，遺傳穩(wěn)定性良好。
2.2.6碳水化合物活性酶注釋
碳水化合物活性酶數(shù)據(jù)庫(kù)（CAZy）是物質(zhì)對(duì)碳水化合物及其衍生物的合成及分解的酶類數(shù)據(jù)資源庫(kù)。本研究通過該數(shù)據(jù)庫(kù)對(duì)鼠李糖乳桿菌Probio-M9碳水化合物活性酶功能基因進(jìn)行分析，在基因組水平觀察鼠李糖乳桿菌Probio-M9菌株在連續(xù)培養(yǎng)100代時(shí)對(duì)碳水化合物利用能力的穩(wěn)定程度。
如圖7所示，從15株不同代時(shí)基因組中共注釋到6大類碳水化合物活性酶的68個(gè)小類基因家族，主要以GHs（糖苷水解酶）和GTs（糖苷轉(zhuǎn)移酶）酶類的含量最為豐富。大部分基因家族在不同代時(shí)無顯著差異，雖然極個(gè)別菌株注釋到的酶類存在不一致情況，但不會(huì)穩(wěn)定遺傳，且總體碳水化合物利用情況相同，遺傳穩(wěn)定性較好，與表型中不同代時(shí)碳水化合物利用能力相同的結(jié)果相吻合。
圖7不同代時(shí)菌株碳水化合物活性酶注釋
3、結(jié)論
遺傳穩(wěn)定性即子代產(chǎn)生的性狀與親本擁有相同的基因型。在生物進(jìn)化過程中，遺傳變異會(huì)經(jīng)常出現(xiàn)，如不嚴(yán)格進(jìn)行人工選擇，就會(huì)導(dǎo)致菌種的衰退，從而導(dǎo)致菌種在生產(chǎn)應(yīng)用過程中出現(xiàn)低產(chǎn)、不穩(wěn)定現(xiàn)象。研究表明，菌株在培養(yǎng)及使用過程中，因保藏方法不同、培養(yǎng)基選擇、傳代代次和人工操作等原因會(huì)導(dǎo)致菌種出現(xiàn)衰退現(xiàn)象。早在1983年，Clements等發(fā)現(xiàn)不同批次用相同方法制備的乳桿菌在治療腹瀉的效果上存在顯著的差異。引發(fā)了人們對(duì)于遺傳穩(wěn)定性的關(guān)注。
近年來，益生菌愈加的受大眾歡迎，而在選擇益生菌時(shí)，不應(yīng)該僅僅考慮菌株的功能，更要參考菌株的遺傳穩(wěn)定性和生存能力等，這些性能關(guān)系到益生菌在腸道中的生存和定殖以及在生產(chǎn)應(yīng)用中其益生特性的發(fā)揮。鼠李糖乳桿菌Probio-M9是2017年從中國(guó)多地的健康母乳中分離出的一株益生菌，多項(xiàng)研究證實(shí)，鼠李糖乳桿菌Probio-M9具有耐酸、耐膽鹽特性，其顯著的益生特性體現(xiàn)在能夠活著進(jìn)入人體腸道，對(duì)腸道菌群起到調(diào)節(jié)作用并增強(qiáng)機(jī)體免疫力，其在食品工業(yè)中廣泛應(yīng)用，而目前所做研究更多的停留在生理生化水平，對(duì)菌株的遺傳穩(wěn)定性研究少之又少。
本研究中，對(duì)鼠李糖乳桿菌Probio-M9在MRS培養(yǎng)基（液體培養(yǎng)基）中連續(xù)培養(yǎng)5個(gè)代時(shí)（0，25，50，75，100代）過程中的表型特征及基因組穩(wěn)定性進(jìn)行分析。采用全基因組及重測(cè)序方法對(duì)不同代時(shí)15株菌進(jìn)行測(cè)序，基因組大小、GC含量表明15株連續(xù)培養(yǎng)的菌株與鼠李糖乳桿菌Probio-M9原始菌株基因組均無明顯差異，系統(tǒng)發(fā)育關(guān)系表明15株菌與原始基因組具有較近的親緣性。在識(shí)別到的SNP位點(diǎn)來看，每?jī)蓚€(gè)傳代代時(shí)間均未發(fā)現(xiàn)穩(wěn)定遺傳的突變位點(diǎn)，且菌株共線性良好。同時(shí)，從基因組水平編碼基因的突變水平證實(shí)了鼠李糖乳桿菌Probio-M9在連續(xù)培養(yǎng)100代過程中具有穩(wěn)定表型特征及代謝特性。
綜上所述，本試驗(yàn)結(jié)果表明鼠李糖乳桿菌Probio-M9在連續(xù)傳代過程中具有良好的遺傳穩(wěn)定性，可為其后期開發(fā)應(yīng)用提供理論參考。
相關(guān)新聞推薦
1、江米酒曲中高產(chǎn)凝乳酶菌篩選、生長(zhǎng)曲線及凝乳活力及蛋白水解活力測(cè)定（三）
2、表達(dá)紅色熒光蛋白的重組鴨腸炎病毒一步生長(zhǎng)曲線、穩(wěn)定性研究（一）
3、草魚腦組織的制備、培養(yǎng)條件及生長(zhǎng)曲線繪制
4、全自動(dòng)生長(zhǎng)曲線分析儀揭示結(jié)核分枝桿菌通過代謝途徑抑制宿主的適應(yīng)性免疫反應(yīng)