基于宏轉錄組學的濃香型白酒酒醅活性微生物群落空間異質性研究

田瑞杰，張勇，馮大鴻，王康麗，遲雷，沈祥坤，胡曉龍，何培新

1.鄭州輕工業(yè)大學 食品與生物工程學院，河南 鄭州 450001;2.河南省食品工業(yè)科學研究所有限公司，河南 鄭州 450053

0 引言

濃香型白酒作為中國傳統(tǒng)蒸餾酒的典型代表，以其獨特的自然固態(tài)發(fā)酵工藝和典型的酒體風味特征而享有盛名[1]。酒醅主要是指處于整個固態(tài)發(fā)酵過程中的多種物料混合物，是白酒釀造微生物進行固態(tài)發(fā)酵的主要載體和白酒風味物質的直接來源。酒醅微生物(細菌、霉菌、酵母菌、產(chǎn)甲烷菌、乳酸菌等)主要來自酒曲、窖泥、操作工具、車間環(huán)境等[2-3]，它們在窖池內經(jīng)過一系列復雜的生化反應和能量代謝，形成多樣的風味物質(有機酸、醇類、酯類等)[3-4]。

近年來，基于DNA水平的免培養(yǎng)技術(DNA序列同源分析、PCR-DGGE(聚合酶鏈反應-變性梯度凝膠電泳)、克隆文庫、二代測序等)已廣泛應用于白酒酒醅微生物群落結構的解析[5-15]，濃香型白酒酒醅微生物群落的結構和多樣性特征逐漸清晰。例如，乳桿菌屬、芽孢桿菌屬、片球菌屬、地芽孢桿菌屬等是酒醅中的優(yōu)勢細菌屬[5,11,13-14]，Kazachstania、Thermoascus、Issatchenkia等是酒醅中的優(yōu)勢真菌屬[8,12,15]。但上述基于DNA水平方法檢測到的微生物大部分處于死亡或休眠狀態(tài)，不能準確反映濃香型白酒發(fā)酵過程中真正有貢獻的活性微生物[1,16-17]。與基于DNA水平的技術相比，基于RNA水平的宏轉錄組學技術能深入了解差異表達基因(DEG)，準確揭示目標環(huán)境樣品中的活性微生物群落，以及與風味物質合成相關的功能微生物及其代謝途徑，目前已廣泛應用于土壤、海洋、湖泊等生態(tài)環(huán)境[18]及豆醬、奶酪、泡菜、白酒等發(fā)酵食品的微生物研究中[19]。例如，Z.W.Song等[7]利用宏轉錄組學技術分析發(fā)現(xiàn)，酵母菌(釀酒酵母屬、畢赤酵母屬、裂殖酵母屬和接合酵母屬)和乳酸菌(乳桿菌屬)中的丙酮酸代謝活動在醬香型白酒風味成分的產(chǎn)生中分為兩個階段：第一階段為乙醇合成，裂殖酵母屬是核心功能微生物；第二階段為乳酸和乙酸合成，乳桿菌屬是核心功能微生物。因此，采用宏轉錄組學技術研究濃香型白酒酒醅中的活性微生物群落結構、代謝反應和功能特征，有助于進一步揭示其固態(tài)發(fā)酵機制和風味物質合成途徑的重要環(huán)節(jié)。

窖池中不同層次酒醅的物料配比、與黃水和窖泥的接觸程度及發(fā)酵后期理化性質的變化均會造成酒醅微生物群落的空間異質性[5,20-21]。例如，胡曉龍等[5]采用16S rRNA基因測序技術對發(fā)酵過程中不同空間位置的濃香型白酒酒醅細菌群落進行解析后發(fā)現(xiàn)，具有相同發(fā)酵節(jié)點的窖池中層，酒醅細菌群落的多樣性高于上層和下層。呂輝等[22]研究發(fā)現(xiàn)，不同空間位置的酒醅微生物生長情況不一樣，具體表現(xiàn)為上層酒醅細菌數(shù)量最多，下層酒醅酵母菌數(shù)量最多。而不同層次酒醅微生物群落的結構差異會導致不同層次酒醅的風味代謝物產(chǎn)生差異[23]。此外，25 d左右是白酒發(fā)酵過程中的重要時間節(jié)點，此時主發(fā)酵期結束，進入酯化階段，部分酒精在酯化反應中合成生香物質，使產(chǎn)酒下降、生香增多[24-25]。目前，對該發(fā)酵節(jié)點濃香型白酒酒醅中的活性微生物群落的結構、差異表達基因、代謝功能的空間異質性研究未見報道。

基于此，本研究擬采用宏轉錄組學技術，以主發(fā)酵期結束后的第25 d作為時間節(jié)點，選取濃香型白酒生產(chǎn)窖池的上、中、下層不同空間位置的酒醅作為研究對象，深入揭示酒醅中具有生理活性的細菌和真菌群落在窖池不同空間位置的組成及差異表達基因，以期為進一步了解濃香型白酒整個發(fā)酵過程中的酒醅活性微生物群落的演替規(guī)律、代謝功能、發(fā)酵機理等提供理論依據(jù)。

1 材料和方法

1.1 材料與試劑

實驗樣品：酒醅樣品，取自河南省某濃香型白酒企業(yè)生產(chǎn)車間。

實驗試劑：十二烷基硫酸鈉(SDS)，上海百賽生物技術股份有限公司產(chǎn)；β-巰基乙醇、異戊醇、聚乙二醇(PEG)、焦碳酸二乙酯(DEPC)，上海麥克林生化科技有限公司產(chǎn)；水飽和酚(pH值為4.5)，上海源葉生物科技有限公司產(chǎn)；無水乙醇，天津市富宇精細化工有限公司產(chǎn)；三羥甲基氨基甲烷(Tris)、聚乙烯吡咯烷酮(PVP)、RNase-Free Water、2×Taq Master Mix、2000 Marker、6×loading buffer，北京索萊寶科技有限公司產(chǎn)；MgCl2、NaCl、CaCl2、草酸鈉、NaH2PO3、Na2HPO3、乙二胺四乙酸二鈉，北京北化精細化學品有限責任公司產(chǎn)。以上化學試劑均為分析純。

1.2 主要儀器與設備

MP200A型精密電子天平，上海良豐儀器儀表有限公司產(chǎn)；TGL-20M型高速離心機，上海盧湘儀離心機儀器有限公司產(chǎn)；Mixer4K微型渦旋混合儀，生工生物工程(上海)股份有限公司產(chǎn)；85-2型恒溫磁力攪拌器，江蘇金壇市中大儀器廠產(chǎn)；IMS-40型全自動雪花制冰機，常熟市雪科電器有限公產(chǎn)；LDZX-50KBS 型立式壓力蒸汽滅菌器，上海中安醫(yī)療器械廠產(chǎn)；SW-CJ-1型超凈工作臺，蘇州凈化設備有限公司產(chǎn)；DYY-8C型電泳儀、WD-9403C型紫外儀，北京六一生物科技有限公司產(chǎn)；NanoDrop2000型微量分光光度計，賽默飛世爾科技公司產(chǎn)；移液槍，上海大龍公司產(chǎn)。

1.3 實驗方法

1.3.1 酒醅樣品采集以發(fā)酵時間為25 d的濃香型白酒酒醅為研究對象，分別在3個不同窖池的上、中、下層位置各取樣50 g，將3個窖池的同一層酒醅樣品混合后得到該層的代表樣品，將其放入無菌袋并置于盛有干冰的泡沫箱中，立即帶回實驗室進行總RNA提取。上、中、下層的代表樣品分別編號為T25、M25、B25。

1.3.2 宏轉錄組測序首先，按照文獻[26]中的SDS-苯酚法對酒醅進行總RNA提??；然后，采用1%瓊脂糖凝膠電泳和微量分光光度計評估總RNA的質量濃度、純度和完整性；最后，選擇高質量的RNA樣本構建Illumina測序文庫。宏轉錄組測序在北京奧維森科技有限公司的 Illumina Hiseq 2500平臺上完成。

1.3.3 宏轉錄組數(shù)據(jù)分析過濾測序獲得的原始序列，剔除帶有測序接頭、未知堿基含量大于1%和低質量堿基(Q≤20)含量大于50%的序列得到質控后序列，利用Diamond軟件將其與NR數(shù)據(jù)庫進行比對分析，并利用Megan軟件進行微生物分類學解析。在物種分類結果中，挑選屬于細菌界和真菌界的微生物，并對其門、綱、屬水平上微生物的豐度進行換算后，以屬水平數(shù)據(jù)為基礎，按照公式①計算所有酒醅樣品中細菌界和真菌界的活性微生物群落的多樣性指數(shù)(Shannon指數(shù))。柱狀圖和熱圖分別采用Origin和Heml軟件進行繪制。

①

式中，H為Shannon指數(shù)，s為微生物屬總數(shù)，Pi為第i個屬的數(shù)量占總數(shù)的比例。

采用 Trinity(v0.27)軟件對質控后序列進行拼接，以拼接后得到的轉錄組作為參考序列(ref)，采用RSEM軟件[26]將每個樣品的質控后序列對ref做映射，得到每個樣品比對每個基因的序列數(shù)目，并對其進行FPKM(每百萬片段中來自某一基因每千堿基長度的片段數(shù)目)轉換，進而分析基因的表達水平。使用統(tǒng)計軟件R中的edgeR包和DESeq2包，分析不同樣品間表達有差異的基因，并使用火山圖對其分布情況進行分析，閾值設定為:|log2(fold change)|>1且q<0.05。使用KOBAS(2.0)軟件，設置參數(shù)--fdr為BH(即使用BH校正)對差異表達基因進行 KEGG 代謝通路富集分析，結果以散點圖形式展現(xiàn)。

2 結果與討論

2.1 測序數(shù)據(jù)結果統(tǒng)計分析

經(jīng)宏轉錄組測序及質量控制后，3個樣品共獲得33.1 G的堿基，各個樣品雙端測序后的序列獲得率均在99%以上，平均測序錯誤率≤0.03%，且Phred質量值分別大于20和30的堿基數(shù)占總堿基數(shù)的比例(Q20、Q30)均大于93.59%(見表1)，表明測序數(shù)據(jù)質量可靠，可用于下一步分析。此外，由于普通的轉錄組測序文庫中Read1和Read2 會以相等的概率取自基因的正反鏈，所以Read1和Read2上的鳥嘌呤(G)和胞嘧啶(C)分布傾向于等比例。在本研究測序數(shù)據(jù)中，同一樣品雙端測序所得的Read1和Read2數(shù)據(jù)中G與C的數(shù)量和占總堿基數(shù)的比例相差較小，即測序數(shù)據(jù)沒有出現(xiàn)明顯的G和C分離，表明測序質量較好，有利于下一步的數(shù)據(jù)分析。

表1 測序數(shù)據(jù)統(tǒng)計Table 1 Sequencing data statistics

2.2 各層酒醅中的活性微生物群落的組成分析

將質控后序列與NR數(shù)據(jù)庫進行比對分析后，各層酒醅樣品中具有轉錄活性的細菌和真菌群落在不同分類水平上的微生物數(shù)目及物種多樣性結果如表2和表3所示。由表2和表3可知，所有酒醅樣品中，細菌群落共檢測到87個門、78個綱、165個目、396個科、1612個屬和7234個種，真菌群落共檢測到8個門、26個綱、58個目、127個科、223個屬和394個種，真菌數(shù)量明顯低于細菌數(shù)量。本研究是從RNA水平揭示濃香型白酒酒醅中的活性微生物群落的分類情況，因此細菌屬(1612個)和真菌屬(223個)的數(shù)量遠高于已報道酒醅中細菌屬(496個)和真菌屬(155個)的數(shù)量[12-13]。此外，不同層酒醅之間的活性微生物數(shù)目差異較小，下層酒醅的細菌數(shù)目在屬和種兩個水平上高于上層和中層，在其他水平上低于其余兩層;真菌數(shù)目則是下層略高于其余兩層。在微生物多樣性方面，上層酒醅的細菌物種多樣性最高，中層的最低；真菌物種多樣性則是中層最高，上層最低，且真菌群落的物種多樣性明顯高于細菌群落。胡曉龍等[5]通過16S rRNA基因測序技術(DNA水平)對濃香型白酒酒醅細菌群落進行研究發(fā)現(xiàn)，發(fā)酵25 d時，微生物Shannon指數(shù)表現(xiàn)為中層最高，上層最低，與本研究結果正好相反，這可能是由于各研究所使用的測序技術不同，方法的局限性和測序的步驟也不同，給研究結果帶來了不同的影響[17]。

表2 各層酒醅樣品中細菌群落在不同分類水平上的微生物數(shù)目及物種多樣性結果Table 2 The number of microorganisms and species diversity of the bacterial communities in each layer of fermented grains samples at different taxonomic levels

表3 各層酒醅樣品中真菌群落在不同分類水平上的微生物數(shù)目及物種多樣性結果Table 3 The number of microorganisms and species diversity of the fungal communities in each layer of fermented grains samples at different taxonomic levels

分別對各層酒醅樣品中的細菌群落和真菌群落在門水平和綱水平上的組成進行進一步解析，結果如圖1和圖2所示。由圖1可以看出，各層酒醅均以厚壁菌門(Firmicutes，95.52%)為主要優(yōu)勢細菌門(圖1a))，以芽孢桿菌綱(Bacilli，94.64%)為主要優(yōu)勢細菌綱(圖1b))，另外，平均相對豐度>1.00%的其他優(yōu)勢細菌門和綱分別只有變形菌門(Proteobacteria)和γ-變形菌綱(Gammaproteobacteria)。本文基于RNA水平的研究結果與之前基于DNA水平解析得到的酒醅優(yōu)勢菌門相比，具有一定的相似性，例如，厚壁菌門和變形菌門在不同地區(qū)濃香型白酒酒醅中均被鑒定為優(yōu)勢細菌門[3,5,11]，表明酒醅的共性特征(如高醇、高酸)可能是富集上述菌門的主要因素[5]。不同層次酒醅之間相比，厚壁菌門和芽孢桿菌綱含量在中層酒醅中較高，變形菌門和γ-變形菌綱含量在上層酒醅中較高，但整體來看，各層酒醅的細菌群落的組成基本一致，無顯著差異。在真菌群落中也是如此，各層酒醅中的優(yōu)勢微生物基本一致，如：上、中、下層酒醅均以子囊菌門(Ascomycota，65.48%)為主要優(yōu)勢真菌門，其他平均相對豐度>1.00%的優(yōu)勢真菌門有擔子菌門(Basidiomycota)、壺菌門(Chytridiomycota)和毛霉菌門(Mucoromycota)(圖2a))，另外，上、中、下層酒醅均以酵母綱(Saccharomycetes，40.27%)為主要優(yōu)勢真菌綱，其他平均相對豐度>1.00%的優(yōu)勢真菌綱有傘菌綱(Agaricomycetes)、散囊菌綱(Eurotiomycetes)、糞殼菌綱(Sordariomycetes)、座囊菌綱(Dothideomycetes)、壺菌綱(Chytridiomycetes)、錘舌菌綱(Leotiomycetes)和毛霉菌綱(Mucoromycetes，除下層外)(圖2b))。胡曉龍等[5]研究發(fā)現(xiàn)，相同發(fā)酵時間節(jié)點下，各層酒醅中的原核微生物群落的結構不存在顯著差異，這與本研究結果具有一定的相似性。雖然本研究選取的酒醅在窖池不同層次的物料(大曲、原料、糟醅、輔料等)配比不同，但各層次酒醅的活性微生物群落的組成基本一致，表明物料配比可能不是影響活性微生物菌群在窖池不同空間位置生長和分布的主要因素。

圖1 各層酒醅樣品中細菌群落在門水平和綱水平上的組成Fig.1 The composition of bacterial communities in each layer of fermented grains samples at phylum level and class level

圖2 各層酒醅樣品中真菌群落在門水平和綱水平上的組成Fig.2 The composition of fungal communities in each layer of fermented grains samples at phylum level and class level

對所有酒醅樣品中平均相對豐度排序前20名的活性細菌屬(94.59%～95.27%)和真菌屬(90.59%～93.22%)分別進行熱圖分析(圖3)，通過樣品聚類，兩者均呈現(xiàn)下層與中層酒醅更為相似的特征。在細菌屬中(圖3a))，3層酒醅均以乳桿菌屬(Lactobacillus，83.34%)為主要優(yōu)勢細菌屬，其他平均相對豐度>1.00%的優(yōu)勢細菌屬為片球菌屬(Pediococcus)和酒球菌屬(Oenococcus)，兩者的平均相對豐度在不同層次酒醅中的差別均很小，如片球菌屬的平均相對豐度在上層為3.13%，在中層為3.24%，在下層為3.02%。與胡曉龍等[5]研究結果相比，濃香型白酒酒醅發(fā)酵后期，以乳桿菌屬為主的乳酸菌不僅在DNA水平上是優(yōu)勢細菌屬，在RNA水平上也處于絕對優(yōu)勢。乳酸菌的主要產(chǎn)物乳酸，能增加白酒的厚重感，降低其刺激感，對白酒風味的形成具有重要影響[14]。另外，第一次在濃香型白酒酒醅中發(fā)現(xiàn)酒球菌屬也為優(yōu)勢菌屬，而針對該菌作用的研究在葡萄酒領域較多，其中的一類菌種(Oenococcusoeni)可觸發(fā)葡萄酒中蘋果酸-乳酸的發(fā)酵，該發(fā)酵過程可降低葡萄酒自身的酸度，提高葡萄酒的品質和穩(wěn)定性[27]。酒球菌屬在濃香型白酒發(fā)酵過程中可能也起到了一定作用，具體影響有待進一步研究。

在真菌屬中，各層酒醅中的優(yōu)勢微生物組成也基本一致，但平均相對豐度略有差異(圖3b))。酵母菌屬(Saccharomyces)在上層(18.87%)和下層(17.92%)酒醅中均為主要優(yōu)勢真菌屬，中層酒醅中則以Scheffersomyces(17.99%)為主要優(yōu)勢真菌屬，而Scheffersomyces在上層(18.86%)和下層(17.45%)酒醅中的平均相對豐度只略低于酵母菌屬。在3層酒醅中共檢測到11個平均相對豐度>1.00%的其他優(yōu)勢細菌屬，包括鵝膏菌屬(Amanita)、籃狀菌屬(Talaromyces)、盤二孢屬(Marssonina)、畢赤酵母屬(Pichia)、滑銹傘屬(Hebeloma)、曲霉屬(Aspergillus)、木霉屬(Trichoderma)、Gelatoporia、Batrachochytrium、Baudoinia和Phlebiopsis，而在下層和中層酒醅中各細菌屬的平均相對豐度更為接近。之前的報道[12-13,15]顯示，發(fā)酵25 d左右時，酒醅中的優(yōu)勢真菌屬在DNA水平上主要為Kazachstania和Thermoascus，而本研究酒醅中的優(yōu)勢真菌屬在RNA水平上主要為酵母菌屬或Scheffersomyces，且Thermoascus的平均相對豐度很低(<0.50%)，Kazachstania未檢測到。造成此差異的原因可能包括以下三方面：1)在DNA水平上，某些含量較高的酒醅微生物可能處于死亡或休眠狀態(tài)，其轉錄活性沒有或較低，而一些含量較低的酒醅微生物則可能具有較高的轉錄活性[1,17]；2)兩種水平上，方法的局限性和測序的步驟不同，可能會給結果帶來偏差[17,28]；3)白酒釀造地區(qū)和釀造工藝存在差異。

圖3 酒醅樣品中平均相對豐度前20名的活性細菌屬和真菌屬的熱圖Fig.3 Heat map of the top 20 active bacterial genera and fungal genera in the average relative abundance of fermented grains samples

2.3 各層酒醅的差異表達基因分析

2.3.1 差異表達基因篩選結果分析圖4為不同層次酒醅之間表達豐度具有差異的基因(DEGs)火山圖，其中上調基因是指表達量上升的基因，下調基因則與之相反。由圖4可以看出，酒醅樣品B25與M25，B25與T25，M25與T25之間的差異表達基因數(shù)量分別為90個、23個、67個，表明中層酒醅分別與下層和上層酒醅之間的差異表達基因數(shù)量較多。B25與M25之間(圖4a))，以中層酒醅作為對照組，下調基因數(shù)量(49個)略高于上調基因數(shù)量(41個)，即下層酒醅的基因表達量較中層酒醅低。M25與T25之間(圖4c))，以上層酒醅作為對照組，下調基因數(shù)量(59個)高于上調基因數(shù)量(8個)，即中層酒醅的基因表達量較上層酒醅低。綜上可知，上層酒醅的基因表達量較中層和下層酒醅高。

圖4 不同層次酒醅之間表達豐度具有差異的基因(DEGs)火山圖Fig.4 Volcano map of genes with different expression abundances (DEGs)between different levels of fermented grain

2.3.2 差異表達基因的KEGG富集分析對不同層次酒醅之間的差異上調基因和下調基因分別進行KEGG代謝通路富集分析，結果如圖5—7所示。其中，Rich factor指差異表達的基因中位于該代謝通路的基因數(shù)目與所有注釋基因中位于該代謝通路的基因總數(shù)的比值，q是經(jīng)多重假設檢驗校正之后的p，其取值范圍為[0,1]，越接近0，表示富集越顯著；圖5挑選了富集最顯著的20條代謝通路進行展示，不足20條的則全部展示。相較于中層酒醅樣品，下層酒醅樣品的上調基因富集到了25條代謝通路上(圖5a))。在一級代謝通路上，有22條通路屬于代謝(Metabolism)、2條通路屬于基因信息加工(Genetic information processing)、1條通路屬于環(huán)境信息加工(Environmental information processing)，其中富集最顯著的通路為RNA降解(RNA degradation)。RNA降解在一級代謝通路上屬于基因信息加工，下層酒醅中該通路的基因表達量較高，RNA降解現(xiàn)象較嚴重，這可能是因為下層酒醅離窖池底部較近，能接觸到較多窖泥，且長時間發(fā)酵后，底部的黃水水位會上升，而較多的水分和雜質均會導致RNA降解。此外，下調基因富集到了16條代謝通路上(圖5b))，在一級代謝通路上，有15條通路屬于代謝、1條通路屬于細胞過程(Cellular processes)，其中富集最顯著的通路為群體效應(Quorum sensing)。群體效應是一種調節(jié)系統(tǒng)，在一級代謝通路上屬于細胞過程，可使細菌共享有關細胞密度的信息并相應地調節(jié)基因表達[29]。中層酒醅中該通路的基因表達量較高，表明中層酒醅中的活性微生物菌群中的基因表達調節(jié)能力強于下層酒醅。

圖5 酒醅樣品B25與M25之間差異上調基因和下調基因的KEGG代謝通路富集圖Fig.5 KEGG pathway enrichment of up-regulated and down-regulated genes between B25 and M25 in fermented grains samples

相較于上層酒醅樣品，下層酒醅樣品的下調基因未富集到任何代謝通路上，上調基因富集到了14條代謝通路上(圖6)。在一級代謝通路上，所有通路均屬于代謝，其中富集最顯著的通路為糖酵解/糖異生(Glycolysis/Gluconeogensis)。糖酵解是將葡萄糖轉化為丙酮酸并產(chǎn)生少量ATP(能量)和NADH(還原力)的過程，也是產(chǎn)生重要前體代謝產(chǎn)物(葡萄糖-6P、果糖-6P、乙酰-CoA等)的主要途徑[30]；糖異生是非碳水化合物前體合成葡萄糖的途徑，本質上是糖酵解的逆轉[31]。兩者在代謝物質生成方面均可發(fā)揮重要作用，在發(fā)酵25 d時，相較于上層酒醅，下層酒醅的活性微生物菌群在代謝活動上貢獻更大。

圖6 酒醅樣品B25與T25之間差異上調基因的KEGG代謝通路富集圖Fig.6 KEGG pathway enrichment of the up-regulated genes between B25 and T25 in fermented grains samples

相較于上層酒醅樣品，中層酒醅樣品的上調基因富集到了3條代謝通路上(圖7a))。在一級代謝通路上，均屬于代謝，其中富集最顯著的通路為硒化合物代謝(Selenocompound metabolism)。此外，下調基因富集到了34條代謝通路上(圖7b))。在一級代謝通路上，有28條通路屬于代謝、3條通路屬于基因信息加工、2條通路屬于環(huán)境信息加工、1條通路屬于細胞過程，其中富集最顯著的通路為RNA降解，其次是糖酵解/糖異生。有研究表明[32]，許多微生物可以耐受高質量濃度的SeO42-和SeO32-，有些微生物還可以將低價態(tài)的硒氧化為高價態(tài)，因此，微生物對硒進行形態(tài)轉化是有效解決硒引起的環(huán)境和健康問題的手段之一。另外，微生物對硒的代謝主要包括硒的轉運、還原、氧化、同化、甲基化等。相較于上層酒醅，中層酒醅中硒化合物的代謝水平較高，這可能是因為其中存在一些可影響硒化合物代謝的活性微生物，具體是哪些微生物尚有待進一步研究確認。此外，上層酒醅中的RNA降解通路的基因表達較高，且活性微生物菌群在代謝活動上的貢獻較大。

圖7 酒醅樣品M25與T25之間差異上調基因和下調基因的KEGG 代謝通路富集圖Fig.7 KEGG pathway enrichment of up-regulated and down-regulated genes between M25 and T25 in fermented grains samples

綜上所述，中層酒醅中的活性微生物菌群對RNA降解的影響最低，下層酒醅中的活性微生物菌群在代謝活動上的貢獻最大。

3 結論

本研究基于宏轉錄組學技術分析了位于窖池不同空間位置的濃香型白酒酒醅(主發(fā)酵期結束后的第25 d)活性微生物群落的組成及差異表達基因。不同酒醅層之間的活性微生物數(shù)目差異較小，上層酒醅的細菌物種多樣性最高，中層酒醅的真菌物種多樣性最高。各層酒醅中的細菌和真菌群落在門、綱和屬水平上的物種組成基本一致，均以厚壁菌門和子囊菌門為主要優(yōu)勢細菌門和真菌門，以芽孢桿菌綱和酵母綱為主要優(yōu)勢細菌綱和真菌綱，以乳桿菌屬為主要優(yōu)勢細菌屬，以酵母菌屬(上層和下層)或Scheffersomyces(中層)為主要優(yōu)勢真菌屬。在差異表達基因表達水平及功能方面，上層酒醅的基因表達水平較中層和下層酒醅更高，各層之間的差異表達基因功能大部分為代謝，且中層酒醅中的活性微生物菌群對RNA降解的影響最低，下層酒醅中的活性微生物菌群在代謝活動上的貢獻最大。本研究可豐富人們對濃香型白酒酒醅菌群的認識，為揭示微生物生態(tài)系統(tǒng)的特定貢獻、開發(fā)生物強化等技術以改善濃香型白酒的品質提供參考。