張守梅 郭玉秋 張娜娜 李克文 孫桂蓮 龔魁杰 王興亞 劉開昌

摘要：為了更加全面了解玉米浸泡副產(chǎn)物中微生物多樣性，比較其中細菌群落結構，利用高通量測序技術，對玉米浸泡液、玉米漿及玉米噴漿粉中細菌16SrRNAV4區(qū)進行測序，進而分析細菌多樣性和群落組成結構。結果表明：共獲得190069條有效數(shù)據(jù)，通過聚類獲得352個OTUs，隸屬于12門、68綱、183屬。α多樣性分析結果表明，玉米噴漿粉中Shannon指數(shù)明顯高于玉米浸泡液和玉米漿。微生物群落組成分析發(fā)現(xiàn)，在門水平上，厚壁菌門（Firmicutes）是玉米浸泡液、玉米漿及玉米噴漿粉的共同優(yōu)勢菌門;在屬水平上，乳桿菌屬（Lactobacillus）為玉米浸泡液、玉米漿及玉米噴漿粉的優(yōu)勢屬;在種水平上，玉米浸泡液中優(yōu)勢菌為耐酸乳桿菌（L.acetotolerans），玉米漿及玉米噴漿粉中優(yōu)勢菌為解淀粉乳桿菌（L.amylolyticus）。

關鍵詞：玉米浸泡副產(chǎn)物;高通量測序;細菌;群落結構;多樣性

中圖分類號：S513：Q939.1 文獻標識號：A 文章編號：1001-4942（2021）12-0143-06

玉米是我國主要糧食和飼料作物，隨著玉米工業(yè)化用途的增加，出現(xiàn)了更多的加工副產(chǎn)品，副產(chǎn)品以廢水的形式直接排放，降低了附加值，并且加劇了環(huán)境的污染。近年來，玉米深加工副產(chǎn)品的綜合利用引起了人們廣泛關注[1]。玉米浸泡液是濕法生產(chǎn)淀粉的副產(chǎn)物之一，玉米經(jīng)亞硫酸浸泡后，顏色為黃色透明，pH值為3.0左右，其中含有大量的蛋白質(zhì)、糖類、維生素等物質(zhì)。浸泡液濃縮即制成玉米漿[2-4]。玉米漿噴到玉米皮上，干燥后制成玉米噴漿粉，可以作為飼料原料出售[1]，或者用于微生物發(fā)酵行業(yè)[5-10]。另外，可以直接從玉米浸泡液中提取蛋白肽、黃色素等物質(zhì)[11-13]。

玉米浸泡液會造成某些特定微生物群體的大量繁殖，主要包括乳桿菌、短粗桿菌、芽孢桿菌、酵母菌等[14]，這些微生物與人和動物的健康密切相關。在玉米浸泡過程中若有乳酸菌和酵母菌大量生長，則可以提高玉米漿的質(zhì)量;若有腐敗性細菌則降低玉米漿的質(zhì)量。過去對微生物的研究主要采用直接鏡檢和傳統(tǒng)分離培養(yǎng)技術，直接鏡檢優(yōu)點是快速得到結果，缺點是無法確定具體的種類，同時對檢驗者知識儲備要求較高。利用傳統(tǒng)的分離培養(yǎng)方法操作簡單且直接，還可以獲得具有利用價值的菌株，是分析微生物群落結構多樣性最基本的方法。但是，該方法工作量大，且僅能培養(yǎng)出極少量微生物，具有很大的局限性，許多微生物在普通培養(yǎng)基上不能生長，不能準確反映環(huán)境中總體微生物群落結構特征。隨著現(xiàn)代分子技術的發(fā)展，人們開始利用非培養(yǎng)技術研究復雜環(huán)境中的微生物群落。高通量測序技術作為新一代全新的測序技術，能夠快速、全面分析不同樣本的群落結構，已被廣泛應用到傳統(tǒng)發(fā)酵食品、釀酒、土壤、糞便、腸道、水體中等微生物多樣性的研究[15-23]。

本研究從玉米浸泡液、玉米漿及玉米噴漿粉中提取細菌基因組DNA，對16SrRNAV4區(qū)進行高通量測序，獲得不同處理階段微生物組成和多樣性，為玉米副產(chǎn)物進一步開發(fā)利用奠定理論基礎。

1 材料與方法

1.1 材料與試劑

玉米浸泡液副產(chǎn)物（玉米浸泡液、玉米漿及玉米噴漿粉）來源于西王集團。

DNA提取試劑盒、瓊脂糖凝膠DNA回收試劑盒（德國Qiagen公司）;TruSeq DNAPCRFreeSamplePreparationKit建庫試劑盒（美國Illumina公司）;PhusionDNA聚合酶（美國NewEnglandBiolabs公司）

1.2 儀器設備

PCR儀（力康生物醫(yī)療科技有限公司）;臺式高速冷凍離心機（美國ThermoFisherScientific公司）;恒溫金屬?。ê贾菝讱W儀器有限公司）;NUGenius凝膠成像系統(tǒng)（GeneCompanyLimited）;XW-80A旋渦混合儀（寧波新芝生物科技有限公司）;水平電泳儀（Bio-Rad公司）;分析天平（德國賽多利斯公司）;IlluminaNovaSeq6000（美國Illumina公司）。

1.3 試驗方法

1.3.1 樣品采集 玉米浸泡液、玉米漿：將采集樣品的玻璃瓶經(jīng)121℃高壓滅菌，裝入無菌自封袋。取樣時，佩戴一次性無菌手套，靠近出樣口時打開瓶蓋接黃漿水，之后擰緊瓶蓋，裝入自封袋中帶回實驗室，一部分于4℃冰箱暫存，一部分于-20℃冰箱保存。

玉米噴漿粉：利用無菌鏟取出，放于無菌塑料袋中，帶回實驗室后一部分于4℃冰箱暫存，一部分于-20℃冰箱保存。

1.3.2 樣品總DNA的提取及PCR擴增 按照試劑盒說明書對樣品基因組DNA進行提取，利用1%瓊脂糖凝膠電泳檢測DNA純度，取適量樣品DNA于離心管中，使用無菌水稀釋至100ng/μL。

以稀釋后的基因組DNA 為模板，選擇細菌16srRNAV4區(qū)進行PCR擴增。擴增產(chǎn)物經(jīng)2%瓊脂糖凝膠電泳檢測，并利用凝膠回收試劑盒回收。

PCR 擴增引物：515F（5′GTGYCAGCMGCCGCGGTAA3′）和806R（5′GGACTACNVGGGTWTCTAA3′）;PCR體系：10×Buffer2μL，dNTP2μL，正反引物各0.8μL，PhusionDNA聚合酶0.2μL，模板DNA1μL，ddH2O補足至20μL。PCR擴增條件：95℃預變性5min;95℃變性30s，50℃退火30s，72℃延伸1min，30個循環(huán);72℃延伸10min;4℃保存。

1.3.3 高通量測序 使用TruSeq DNAPCRFreeSamplePreparationKit建庫試劑盒進行文庫構建，構建好的文庫經(jīng)Qubit和qPCR定量，文庫合格后，使用NovaSeq6000進行上機測序（北京諾禾致源科技股份有限公司）。

1.3.4 數(shù)據(jù)處理 根據(jù)Barcode序列和PCR擴增引物序列從下機數(shù)據(jù)中拆分出各樣本數(shù)據(jù)，截去Barcode和引物序列后使用FLASH（V1.2.7，http：//ccb.jhu.edu/software/FLASH/）對每個樣本的Reads進行拼接，得到的拼接序列為原始Tags數(shù)據(jù)（RawTags）;拼接得到的RawTags經(jīng)過嚴格的過濾處理[19]得到高質(zhì)量的Tags數(shù)據(jù)（CleanTags）。參照Qiime（V1.9.1，http：//qiime.org/scripts/split_libraries_fastq.html）Tags質(zhì)量控制流程，得到最終的有效數(shù)據(jù)（EffectiveTags）。利用Uparse軟件對所有樣本的有效數(shù)據(jù)進行聚類，默認以97%的一致性將序列聚類成為OTUs（OperationalTaxonomicUnit）的代表序列。

對OTUs序列進行物種注釋，用Mothur方法比對數(shù)據(jù)庫進行物種注釋分析，同時計算該OTUs在各樣本中的相對含量。使用Qiime軟件對樣品的多樣性指數(shù)進行分析。

2 結果與分析

2.1 細菌多樣性分析

樣品高通量測序后共獲得序列312240條，經(jīng)過對原始數(shù)據(jù)拼接、過濾，共得到190069條有效數(shù)據(jù)，通過聚類獲得352個OTUs。為了直觀表現(xiàn)出玉米浸泡液（N1）、玉米漿（N2）、玉米噴漿粉（N3）共有OTUs和特有OTUs，對獲得的OTUs進行了統(tǒng)計和比較分析。由圖1可知，玉米浸泡液、玉米漿、玉米噴漿粉共有OTUs數(shù)量為50個，玉米浸泡液中特有OTUs為52個，玉米漿中特有OTUs為33個，玉米噴漿粉中特有OTUs為128個，表明玉米噴漿粉中微生物群落較玉米浸泡液和玉米漿豐富。

α多樣性分析反映了單個樣品內(nèi)物種的豐富度和多樣性。樣本的覆蓋率表明測序的數(shù)據(jù)能夠覆蓋目前狀態(tài)下樣品中細菌的種類，由表1可知，樣品中覆蓋率均為1，說明樣本中被測出的概率較高，能真實反映樣本中細菌物種豐度及多樣性[24]。OTUs數(shù)可以代表樣品物種的豐度[16]，ACE指數(shù)是用來估計群落中OTUs數(shù)目的指數(shù)，玉米噴漿粉的OTUs數(shù)量最多，ACE指數(shù)最高，其次是玉米浸泡液，因此玉米浸泡液副產(chǎn)物由濃縮到制粉過程中，微生物豐度發(fā)生了先降低后增加的變化。Chao1指數(shù)可估計樣品群落中包含的物種總數(shù)，值越大表明物種總數(shù)越多，由浸泡液到制成噴漿粉的過程物種總數(shù)發(fā)生了先降低后增加的變化。Shannon指數(shù)表示樣品中群落多樣性，指數(shù)越大群落多樣性越高，多樣性由高到低依次為玉米噴漿粉、玉米漿、玉米浸泡液。Simpson指數(shù)反映了樣品群落中物種均勻度，玉米噴漿粉中的均勻度最高。

2.2 細菌群落結構差異分析

2.2.1 門水平下玉米浸泡液副產(chǎn)物中細菌群落結構分析 在門水平下共鑒定出12個門類細菌，圖2顯示了玉米浸泡液（N1）、玉米漿（N2）、玉米噴漿粉（N3）中豐富度前10的物種，其余的合并為others。優(yōu)勢細菌門是厚壁菌門（Firmicutes），從濃縮到制粉的過程中，厚壁菌門的相對豐度逐漸減少，變形菌門（Proteobacteria）和藍菌門（Cyanobacteria）的相對豐度增加。變形菌門在玉米浸泡液、玉米漿、玉米噴漿粉中相對豐度分別為1.1%、14.2%和28.0%，呈逐漸上升趨勢，原因可能為變形菌門中的細菌為好氧型，不能適應玉米浸泡液高溫、高酸、缺氧的極端環(huán)境;玉米噴漿粉中營養(yǎng)物質(zhì)豐富，有氧的環(huán)境促進了變形菌門的生長。

2.2.2 屬水平下玉米浸泡液副產(chǎn)物中的細菌群落結構分析 在屬水平下共檢測出183個屬的細菌，圖3顯示屬級豐富度前30的細菌，其他的合并為others，圖4顯示豐富度前100屬的細菌系統(tǒng)發(fā)育分析。由圖3、圖4可知，玉米浸泡液、玉米漿、玉米噴漿粉的微生物主要為厚壁菌門（Firmicutes）下的乳桿菌屬（Lactobacillus），其相對豐度分別為97%、83%、23%。由浸泡液到噴漿粉的過程中許多好氧型細菌增加，如假單胞菌屬（Pseudomonas）、沙雷氏菌屬（Serratia）等。玉米浸泡液中微生物群落結構是一個動態(tài)變化的過程，初期乳桿菌大量繁殖，產(chǎn)生乳酸，抑制了異種微生物生長，后期乳酸對自身起到了限制作用;在玉米漿中，除乳桿菌外的其他細菌生長，相對豐度增加;由玉米漿到噴漿粉，許多好氧菌迅速增加，微生物豐富度增加。

2.2.3 種水平下玉米浸泡液副產(chǎn)物中細菌群落結構分析 鑒定到種的細菌有101個，隸屬于乳桿菌屬的細菌有耐酸乳桿菌（L.acetotolerans）、解淀粉乳桿菌（L.amylolyticus），還包括少量的嗜淀粉乳桿菌（L.amylovorus）、開菲爾基質(zhì)乳桿菌（L.kefiranofacie，圖5）。玉米浸泡液的優(yōu)勢菌為耐酸乳桿菌，相對豐度為80.6%。玉米漿和玉米噴漿粉優(yōu)勢菌為解淀粉乳桿菌，相對豐度分別為82.5%、17.0%。玉米噴漿粉中還存在粘質(zhì)沙雷氏桿菌（Serratiamarcescens）、莓實假單胞菌（P.fragi）、耐熱芽孢桿菌（Bacillussporothermodurans），這些有害細菌不僅具有食物致腐能力，還是條件致病菌，在機體免疫力下降的情況下引起肺部、尿路感染甚至敗血癥。芽孢桿菌革蘭氏染色陽性，需氧或兼性厭氧，具有抵御不良環(huán)境的能力。許多芽孢桿菌具有解磷、固氮、廣譜性抗菌、高產(chǎn)淀粉酶、蛋白酶等功能，在農(nóng)業(yè)生產(chǎn)、食品加工、工業(yè)生產(chǎn)、醫(yī)療等方面起到了重要作用。但是也有一些芽孢桿菌具有致病性，產(chǎn)生耐熱毒素，引起食物中毒;因其獨特的耐高溫特性，普通的熱處理很難將其消滅，給食品質(zhì)量安全造成隱患。芽孢桿菌在玉米浸泡液和玉米漿中主要以芽孢的形式存在，常規(guī)的提取方法很難獲得高質(zhì)量的DNA用于PCR反應，因此，高通量測序顯示玉米浸泡液和玉米漿中芽孢桿菌數(shù)量極少。

在玉米浸泡液和玉米漿中的細菌主要為耐酸乳桿菌和解淀粉乳桿菌，其他微生物相對豐度較低。分析原因：其一是因為相對缺氧的、酸性較高的液體環(huán)境適合耐酸乳桿菌的生長;其二，許多研究發(fā)現(xiàn)乳桿菌可以產(chǎn)生細菌素等拮抗物質(zhì)，還可以通過與其它微生物競爭底物等途徑影響其它微生物的生長[25-28]。

3 討論與結論

玉米浸泡液是微生物生長、繁殖代謝的場所，本研究通過高通量測序技術對玉米浸泡液、玉米漿和玉米噴漿粉的細菌群落結構和多樣性進行分析，揭示了玉米副產(chǎn)物在不同加工過程中微生物群落結構變化。由于浸泡液含豐富的營養(yǎng)物質(zhì)，酸性環(huán)境較強，因此乳桿菌豐富，生長旺盛。濃縮成玉米漿后，隨著時間的增加，乳桿菌的生長受到抑制，其他種類的細菌相對增多。α多樣性分析表明，玉米噴漿粉中OUTs數(shù)量最多，群落結構中豐富度、多樣性、均勻度最高。該研究結果可為玉米副產(chǎn)物的保存、開發(fā)利用提供借鑒和參考。

本研究表明：玉米浸泡液、玉米漿和玉米噴漿粉的優(yōu)勢菌門是厚壁菌門，且厚壁菌門的相對豐度依次減少。變形菌門在玉米浸泡液、玉米漿、玉米噴漿粉中逐漸增多。在屬水平上，乳桿菌屬為玉米浸泡液、玉米漿、玉米噴漿粉的優(yōu)勢菌屬，且其相對豐度依次減少。劉慶艾等[29]研究了H2O2玉米浸泡工藝及其浸泡液中菌群結構的變化，結果表明，浸泡初期細菌結構組成較為豐富，一段時間后，Lactobacillaceae生長繁殖成為優(yōu)勢菌種并伴隨OTUs水平的降低，浸泡末期Lactobacillaceae進入衰亡期，其他細菌數(shù)量明顯增加。金渭武等[30]分析了逆流浸泡過程中微生物的變化，發(fā)現(xiàn)整個浸泡過程中，Lactobacillaceae占絕對優(yōu)勢并在浸泡前期發(fā)酵產(chǎn)生大量乳酸。乳桿菌作為環(huán)境中的優(yōu)勢菌種，能夠消耗有限的營養(yǎng)物質(zhì)，與病原菌競爭生存空間從而抑制病原菌的生長繁殖[31];另外，在生長過程中分泌的抗菌素物質(zhì)可抑制其他病原菌生長，產(chǎn)生的乳酸降低了環(huán)境中的pH值也能抑制其他病原菌的生長繁殖[32]。因此玉米浸泡液和玉米漿更有利于玉米浸泡液副產(chǎn)物的儲存和質(zhì)量保證。