歐揚(yáng)雯珊 劉松玲 劉治麟 陳巧燕 夏兵兵 趙 亮，*

（1 教育部功能乳品重點(diǎn)實(shí)驗(yàn)室 中國(guó)農(nóng)業(yè)大學(xué)食品科學(xué)與營(yíng)養(yǎng)工程學(xué)院 北京 100083

2 北京市高等學(xué)校畜產(chǎn)品工程技術(shù)研究中心 北京 100083

3 北京和益源生物技術(shù)有限公司 北京 100088

4 深圳市晨光乳業(yè)有限公司 廣東深圳 518107）

雙歧桿菌（Bifidobacterium）是一類安全有效的益生菌，然而其環(huán)境耐受與營(yíng)養(yǎng)供應(yīng)問題一直影響著它的有效利用。雙歧桿菌定植能力弱[1-3]，其豐度受腸道內(nèi)營(yíng)養(yǎng)物質(zhì)含量、種類的制約[4]。特定糖類則可促進(jìn)雙歧桿菌生長(zhǎng)[4]，幫助其定植。

寡糖（亦稱為低聚糖）是指由3～10個(gè)單糖組成的糖類[5]。根據(jù)生理功能，寡糖可分為消化性寡糖與無法被人體利用的非消化性寡糖[5-7]，后者可被雙歧桿菌代謝利用[8-11]，因此膳食會(huì)影響腸道微環(huán)境[12]。不同菌株對(duì)糖的利用范圍、種類、方式存在特異性[13]，如嬰兒雙歧桿菌主要利用母乳寡糖[14]、兩歧雙歧桿菌主要利用 mucin寡糖[15]。Vulevic等[9-10]發(fā)現(xiàn)低聚半乳糖可提高老年人與肥胖成人腸道中雙歧桿菌含量，分別提升IL-10水平與降低總膽固醇含量；Lewis等[11]發(fā)現(xiàn)喂養(yǎng)含巖藻糖基化寡糖母乳的嬰兒，腸道中雙歧桿菌豐度顯著高于對(duì)照組。

不同于多數(shù)嬰兒來源的雙歧桿菌，來自中國(guó)廣西巴馬長(zhǎng)壽老人腸道的BBMN68（Bifidobacterium longum subsp.longum BBMN68）對(duì)寡糖的利用具有特異性[15-17]。該菌株具有增強(qiáng)免疫[18]，緩解過敏[19]，改善便秘[20]的功能。生物信息學(xué)技術(shù)從基因?qū)用娼沂玖司陮?duì)環(huán)境的適應(yīng)特性，是研發(fā)菌株?duì)I養(yǎng)利用的有效手段[15-16]。本研究旨在通過生物信息學(xué)技術(shù)預(yù)測(cè)BBMN68糖苷水解酶及碳水化合物轉(zhuǎn)運(yùn)系統(tǒng)，獲得潛在可利用寡糖，對(duì)開發(fā)利用特有BBMN68益生元具有重要意義。

1 材料與試劑

1.1 試驗(yàn)菌株

B.longum subsp.longum BBMN68（CGMCC No.2265；GenBank assembly accession GCA_000166315.1）分離自廣西巴馬長(zhǎng)壽老人糞便，由教育部功能乳品重點(diǎn)實(shí)驗(yàn)室保藏。

1.2 主要儀器（表1）

表1 主要儀器與設(shè)備Table1 Main instruments and equipment

1.3 主要試劑及配制

1.3.1 菌株培養(yǎng) 碳源：低聚甘露糖、低聚木糖、低聚半乳糖、低聚異麥芽糖、低聚麥芽糖、低聚果糖、葡萄糖。

半合成培養(yǎng)基[21]：細(xì)菌蛋白胨10 g/L；酵母浸粉5 g/L；磷酸氫二鉀 2 g/L；乙酸鈉 5 g/L；檸檬酸氫二胺 2 g/L；MnSO4·4H2O 0.25 g/L；MgSO4·7H2O0.58 g/L；Tween-80 1 g/L；Cys-HCl 0.5 g/L。

1.3.2 定量PCR細(xì)菌RNAout試劑盒，天恩澤公司；cDNA第一鏈合成試劑盒、SYBR Green I Real-time PCR Master Mix，日本TAKARA公司。

2 試驗(yàn)方法

2.1 BBMN68寡糖相關(guān)基因功能預(yù)測(cè)

2.1.1 糖苷水解酶分析 使用CAZy數(shù)據(jù)庫(kù)（http：//www.cazy.org）比對(duì)分析BBMN68基因組，獲取糖苷水解酶編碼基因。利用PSORTb數(shù)據(jù)庫(kù)（http：//www.psort.org/psortb/）定位上述基因編碼的糖苷水解酶在細(xì)胞中的位置。通過EC數(shù)據(jù)庫(kù)（http：//www.ebi.ac.uk/）對(duì) BBMN68 中糖苷水解酶進(jìn)行編碼歸類。

2.1.2 寡糖轉(zhuǎn)運(yùn)子分析 獲取BBMN68基因組中所有COG（G）蛋白，使用 TCDB 數(shù)據(jù)（http：//www.tcdb.org/）庫(kù)BLAST功能，對(duì)所有蛋白進(jìn)行檢索分析。根據(jù)TCDB數(shù)據(jù)庫(kù)分類信息：分類號(hào)3.A.1.1及3.A.1.2為碳水化合物ABC轉(zhuǎn)運(yùn)子；2.A.1.X為MFS轉(zhuǎn)運(yùn)子；8.A.7及8.A.8為PTS轉(zhuǎn)運(yùn)系統(tǒng)[22]。

2.2 功能驗(yàn)證

2.2.1 生長(zhǎng)試驗(yàn)驗(yàn)證 按照2%接種量，將BBMN68接種于含不同寡糖的半合成培養(yǎng)基中，37℃厭氧培養(yǎng)12 h，測(cè)定OD600條件下的吸光值。

2.2.2 定量PCR 總RNA的提取及cDNA合成：使用TIANDZ細(xì)菌RNAout試劑盒提取對(duì)數(shù)期長(zhǎng)雙歧桿菌BBMN68的總RNA。收集5～8mL新鮮菌液于離心管中，4℃，12 000×g離心1min，棄盡上清；加入1mL細(xì)菌RNAout并用槍充分吹打，確保細(xì)菌充分裂解，沒有塊狀物；加入0.2mL氯仿，充分振蕩混勻，以打斷DNA并去除蛋白質(zhì)；4℃，12 000×g離心3min，吸取無色上清液（注意不要吸取中間層的蛋白質(zhì)）于新的離心管中，并加入等體積異丙醇，輕輕顛倒混勻；4℃，12 000×g離心10min，吸棄上清；加入75%乙醇，充分振蕩混勻，4℃，12 000×g離心 3min，吸棄上清；加入 50μL的0.1%DEPC水沉淀溶解。將RNA用1.2%的瓊脂糖凝膠電泳檢驗(yàn)所提取RNA的完整度及純度。使用cDNA第一鏈合成試劑盒，作為后續(xù)實(shí)時(shí)熒光定量PCR的模板。

引物的設(shè)計(jì)、合成與實(shí)時(shí)熒光定量PCR：利用NCBI數(shù)據(jù)庫(kù)中PRIMER BLAST在線引物設(shè)計(jì)軟件設(shè)計(jì)引物，并對(duì)選自文獻(xiàn)中的引物進(jìn)行比對(duì)驗(yàn)證。引物由上海英駿生物技術(shù)有限公司合成（表3）。PCR 反應(yīng)體系 25μL：cDNA 模板 1μL，上下游引物各 1 μL，SYBR Green I Real-time PCR Master Mix 12.5 μL，ddH2O 9.5 μL。反應(yīng)條件：95℃預(yù)變性10 s，95℃變性5 s，61℃退火延伸31 s，45個(gè)循環(huán)。相對(duì)表達(dá)量計(jì)算公式：相對(duì)表達(dá)量=2[（CT目的基因-CT內(nèi)參基因）對(duì)照組-（CT目的基因-CT內(nèi)參基因）處理組]。在以葡萄糖為唯一碳源的半合成培養(yǎng)基中BBMN68基因表達(dá)量作為對(duì)照組，在以其它寡糖為唯一碳源的半合成培養(yǎng)基中其基因表達(dá)量為處理組。CT目的基因指BBMN68中糖苷水解酶基因拷貝數(shù)，CT內(nèi)參基因指 BBMN68 16S rRNA拷貝數(shù)。每個(gè)目標(biāo)基因進(jìn)行3個(gè)生物重復(fù)。

表2 RT-PCR引物Table2 The primer of RT-PCR

3 結(jié)果與分析

3.1 BBMN68寡糖相關(guān)基因功能預(yù)測(cè)

3.1.1 糖苷水解酶預(yù)測(cè) CAZy數(shù)據(jù)庫(kù)顯示BBMN68基因組中含有58個(gè)編碼糖苷水解酶的基因。通過PSORTb數(shù)據(jù)庫(kù)進(jìn)一步對(duì)58個(gè)編碼蛋白進(jìn)行亞細(xì)胞定位，發(fā)現(xiàn)其中有16個(gè)蛋白位于胞外，占總數(shù)的28%。全部糖苷水解酶可歸納為24個(gè)蛋白家族，其中 GH13家族數(shù)量最多（13），GH43家族其次（9），且GH43家族糖苷水解酶均位于胞外或錨定于細(xì)胞壁。

結(jié)果顯示，BBMN68主要以細(xì)胞壁多糖及淀粉為底物。以植物細(xì)胞壁為底物的糖苷水解酶數(shù)目最多，包括26個(gè)編碼基因，其中12個(gè)位于胞外或細(xì)胞壁；其次是以淀粉為底物的糖苷水解酶，包括15個(gè)編碼基因，全部屬于胞內(nèi)酶。

EC數(shù)據(jù)庫(kù)顯示，33個(gè)糖苷水解酶可注釋EC編碼，共編碼14種酶。水解獲得的單糖產(chǎn)物包括半乳糖、葡萄糖、N-乙酰葡萄糖胺、N-乙酰半乳糖胺、果糖、甘露糖、阿拉伯糖和木糖。

表3 BBMN68糖苷水解酶Table3 Glycoside hydrolyses of BBMN68

（續(xù)表3）

3.1.2 寡糖轉(zhuǎn)運(yùn)系統(tǒng)預(yù)測(cè) 與糖苷水解酶不同，多數(shù)寡糖的轉(zhuǎn)運(yùn)蛋白為膜蛋白，限制了該類蛋白的深入研究。目前，學(xué)者主要通過蛋白質(zhì)組學(xué)[23]、轉(zhuǎn)錄組學(xué)[24]和糖組學(xué)[25]，開展對(duì)雙歧桿菌寡糖轉(zhuǎn)運(yùn)系統(tǒng)底物的研究。TCDB數(shù)據(jù)庫(kù)顯示BBMN68寡糖轉(zhuǎn)運(yùn)系統(tǒng)主要包括ABC（ATP binding cassette）轉(zhuǎn)運(yùn)系統(tǒng)和PEP-PTS（Phosphoenolpyruvate-phosphotransferase）轉(zhuǎn)運(yùn)系統(tǒng)，以及MFS（Major facilitator superfamily）系統(tǒng)及少量的MIP（Major intrinsic protein）[26]系統(tǒng)。

表4 BBMN68碳水化合物轉(zhuǎn)運(yùn)子Table4 Carbohydrate transporters in BBMN68

（續(xù)表4）

TCDB數(shù)據(jù)庫(kù)顯示，BBMN68基因組共編碼57個(gè)碳水化合物轉(zhuǎn)運(yùn)子基因（表5），包括ABC轉(zhuǎn)運(yùn)子41個(gè)，MFS轉(zhuǎn)運(yùn)基因11個(gè)，PTS轉(zhuǎn)運(yùn)系統(tǒng)基因3個(gè)和MIP基因2個(gè)，介于已報(bào)道的雙歧桿菌碳水化合物轉(zhuǎn)運(yùn)基因數(shù)目（21～68）[24]。通常雙歧桿菌MFS基因數(shù)目在1～6個(gè)之間，而BBMN68具有11個(gè)MFS基因。

BBMN68中ABC轉(zhuǎn)運(yùn)子基因多位于基因簇內(nèi)，其中，36個(gè)基因分布于11個(gè)基因簇中，僅5個(gè)基因單獨(dú)存在。通常，每個(gè)基因簇包括1～3個(gè)糖特異性連接蛋白，1～2個(gè)通透酶和1～2個(gè)ATP酶。結(jié)果顯示，BBMN68_925-BBMN68_927編碼3個(gè)糖連接蛋白，BBMN68_933-BBMN68_934編碼2個(gè)連續(xù)的通透酶。不同于已有文獻(xiàn)，11個(gè)ABC轉(zhuǎn)運(yùn)子基因簇中僅有2個(gè)基因簇（BBMN68_73-BBMN68_78 與 BBMN68_1723-BBMN68_1730）包含ATP連接蛋白基因。而BBMN68基因組中發(fā)現(xiàn)一個(gè)單獨(dú)存在的ATP連接蛋白（BBMN68_1403），該蛋白氨基酸序列與已報(bào)道的BL0673相似度大于60%。該蛋白在B.longum NCC2705與Streptomycescoelicolor A3中的注釋信息為ABC轉(zhuǎn)運(yùn)子通用ATP酶，該酶可被其它缺少ATP酶的ABC轉(zhuǎn)運(yùn)子共用[27]，因此 BBMN68_1403可能作為BBMN68 ABC轉(zhuǎn)運(yùn)子的通用ATP酶。該現(xiàn)象可幫助有限的基因發(fā)揮最大的作用，體現(xiàn)了BBMN68的遺傳優(yōu)勢(shì)。此外，ABC轉(zhuǎn)運(yùn)子基因簇常與糖苷水解酶基因及轉(zhuǎn)錄調(diào)節(jié)子基因毗鄰，表明轉(zhuǎn)錄子的特異性可能與糖苷水解酶的特異性相關(guān)，基因簇的表達(dá)可能受底物水平的調(diào)控。

圖1 碳水化合物轉(zhuǎn)運(yùn)子基因簇Fig.1 Gene clusters of carbohydrate transporters

3.2 BBMN68寡糖相關(guān)基因功能驗(yàn)證

3.2.1 生長(zhǎng)試驗(yàn)驗(yàn)證 生長(zhǎng)試驗(yàn)結(jié)果顯示BBMN68對(duì)不同來源的寡糖具備利用能力。BBMN68在植物細(xì)胞壁來源的低聚半乳糖、低聚木糖培養(yǎng)基中生長(zhǎng)良好。菌株在淀粉不完全降解寡糖產(chǎn)物——低聚異麥芽糖和低聚麥芽糖中亦可良好生長(zhǎng)，而在未降解的淀粉中不生長(zhǎng)（數(shù)據(jù)未顯示）。結(jié)果顯示BBMN68可在低聚果糖中生長(zhǎng)，無法在動(dòng)物來源的低聚甘露糖中生長(zhǎng)。

3.2.2 碳水化合物轉(zhuǎn)運(yùn)基因的驗(yàn)證 RT-PCR結(jié)果顯示，在寡糖半合成培養(yǎng)基中，碳水化合物轉(zhuǎn)運(yùn)子基因mRNA水平表達(dá)量發(fā)生了不同程度的變化，存在以下特點(diǎn)：（1）不同基因?qū)烟堑姆磻?yīng)不同，例如BBMN68_925，BBMN68_1022，BBMN68_1724基因?qū)υ囼?yàn)底物的反應(yīng)較低，表明試驗(yàn)所用寡糖可能并非其真正底物，或者上述基因?qū)儆诮M成型表達(dá)。（2）單糖組成相似，碳鏈長(zhǎng)度不同的底物誘導(dǎo)表達(dá)的基因相似，如低聚果糖、果糖和蔗糖（數(shù)據(jù)未顯示）均可誘導(dǎo)BBMN68_74，BBMN68_75表達(dá)。（3）同種底物可誘導(dǎo)多個(gè)轉(zhuǎn)運(yùn)子基因表達(dá)，例如：低聚麥芽糖可誘導(dǎo)BBMN68_1170，BBMN68_1422，BBMN68_1024，BBMN68_1169，BBMN68_1170，BBMN68_1422，BBMN68_1427，BBMN68_1597，BBMN68_1671基因表達(dá)。

表5 BBMN68預(yù)測(cè)寡糖底物中的生長(zhǎng)情況Table5 Growth profile of BBMN68 in predicted oligosaccharide substrates

圖2 轉(zhuǎn)運(yùn)子基因在預(yù)測(cè)寡糖底物中的表達(dá)變化Fig.2 Relative transcription levels of carbohydrate transporter-encoding genes from BBMN68 upon cultivation in MRS medium supplemented with various oligosaccharides

4 討論

本文預(yù)測(cè)并驗(yàn)證了BBMN68對(duì)多種寡糖的影響，發(fā)現(xiàn)BBMN68具有轉(zhuǎn)運(yùn)和降解低聚木糖、低聚半乳糖、低聚異麥芽糖、低聚麥芽糖、低聚果糖的酶系統(tǒng)，可有效利用上述寡糖。

研究顯示，BBMN68不能利用動(dòng)物來源的低聚甘露糖，推測(cè)原因：（1）該菌株可能缺少特異性的甘露糖轉(zhuǎn)運(yùn)子；（2）該菌株可能缺少利用甘露糖的基因。低聚甘露糖糖苷水解酶將低聚甘露糖水解成甘露糖單體，提供給其它腸道微生物，實(shí)現(xiàn)共生。

BBMN68糖類酶系明顯異于其它雙歧桿菌。BBMN68中26%（16/58）的糖苷水解酶位于胞外，比例顯著高于多數(shù)雙歧桿菌。其中，有12個(gè)胞外糖苷水解酶的預(yù)測(cè)底物為植物細(xì)胞壁多糖。結(jié)果顯示BBMN68可在常見植物細(xì)胞壁多糖降解寡糖產(chǎn)物中生長(zhǎng)，卻無法在未降解的多糖中生長(zhǎng)（數(shù)據(jù)未顯示）。推測(cè)原因：植物細(xì)胞壁多糖水解酶僅有外切活性，而植物細(xì)胞壁多糖結(jié)構(gòu)復(fù)雜，需多種糖苷水解酶共同降解。如B.longum B667中GH51家族的糖苷水解酶AbfB，僅可內(nèi)切阿拉伯聚糖、阿拉伯木聚糖等多糖，推測(cè)AbfB與其它腸道共生菌群分泌的糖苷水解酶，共同實(shí)現(xiàn)對(duì)多糖的降解[28]。研究發(fā)現(xiàn)，將長(zhǎng)雙歧桿菌與擬桿菌共同定植于無菌小鼠腸道后，擬桿菌細(xì)胞壁多糖降解基因表達(dá)量顯著高于擬桿菌單獨(dú)定植時(shí)，表明雙歧桿菌對(duì)擬桿菌的多糖代謝過程有促進(jìn)作用[29]。胞外植物細(xì)胞壁多糖降解基因幫助BBMN68實(shí)現(xiàn)與其它腸道微生物的和諧共生。

本文發(fā)現(xiàn)，雙歧桿菌優(yōu)勢(shì)底物——淀粉，無法被BBMN68利用。BBMN68不具備胞外淀粉酶，無法降解淀粉，僅能利用糊精等淀粉部分降解產(chǎn)物，該特點(diǎn)與BBMN68的特殊來源相關(guān)。廣西巴馬長(zhǎng)壽人群的飲食以未經(jīng)深加工的玉米及蔬菜為主，玉米是抗性淀粉的主要來源。攝入玉米后，小腸及糞便中可檢測(cè)到較高含量的糊精[30]。宿主為BBMN68提供足量的糊精，以保證其在腸道中穩(wěn)定的生態(tài)位。長(zhǎng)壽老人腸道中含有大量擬桿菌，可降解膳食纖維，為BBMN68提供除淀粉以外的纖維寡糖，保證BBMN68在腸道中的定植生長(zhǎng)。

受限于寡糖組成結(jié)構(gòu)的多樣性和糖苷水解酶的特異性兩點(diǎn)因素，本研究尚不能對(duì)全部寡糖進(jìn)行精準(zhǔn)預(yù)測(cè)。本文僅選取部分代表性寡糖，證明了寡糖代謝對(duì)BBMN68生態(tài)適應(yīng)的重要性，對(duì)開發(fā)菌株特異性寡糖有重要意義。