汾酒酒醅發(fā)酵過程中真菌群落的變化規(guī)律

衛(wèi)春會，甄 攀，張?zhí)m蘭，任志強，黃治國，鄧 杰,

（1.四川輕化工大學(xué) 釀酒生物技術(shù)及應(yīng)用四川省重點實驗室，四川 宜賓 644000；2.山西杏花村汾酒廠股份有限公司，山西 汾陽 032205）

汾酒是清香型白酒的典型代表，汾酒入口綿、落口甜、尾味余香、回味悠長的特色與其采用傳統(tǒng)的地缸固態(tài)分離發(fā)酵和清蒸二次清工藝有關(guān)[1-2]。但是，傳統(tǒng)汾酒開放式的生產(chǎn)方式導(dǎo)致了汾酒發(fā)酵過程中微生物極其復(fù)雜的變化，微生物來源多種多樣[3]。甄攀等[4]研究了汾酒釀造微生物的來源，發(fā)現(xiàn)汾酒酒醅中微生物主要來自于釀造場地的空氣以及生產(chǎn)器具。韓莎等[5]通過可培養(yǎng)技術(shù)分析汾酒釀造過程中微生物群落結(jié)構(gòu)，并找出微生物種群與代謝物之間的關(guān)系。Xiao Ran等[6]基于核糖體高通量測序分析汾酒發(fā)酵過程中的微生物，認(rèn)為汾酒發(fā)酵過程中主要細(xì)菌主要是乳桿菌科，真菌主要是酵母科（Saccharomycetaceae）和復(fù)膜孢酵母科（Saccharomycopsidaceae）。雷振河[7]通過高通量測序手段分析清香型白酒地缸中心與邊緣位置酒醅微生物種群結(jié)構(gòu)差異。王雪山等[8]利用高通量測序技術(shù)分析不同層次清香型白酒酒醅中微生物種群結(jié)構(gòu)及演替規(guī)律，共檢測到161 個細(xì)菌屬和161 個真菌屬，其中乳桿菌屬（Lactobacillus）、畢赤酵母屬（Pichia）和假絲酵母屬（Candida）是酒醅中的絕對優(yōu)勢微生物。不同層次酒醅微生物多樣性變化規(guī)律不同，而且微生物種群結(jié)構(gòu)及演替規(guī)律存在明顯差異。目前，國內(nèi)對清香型白酒釀造微生態(tài)的研究多數(shù)是針對發(fā)酵過程中酒醅微生物種群演替。不同階段酒醅微生物存在很大的差異，從而造成了不同層次的白酒風(fēng)味特征，因此，研究汾酒發(fā)酵微生物群落種類和發(fā)酵規(guī)律，對提高汾酒品質(zhì)有重要理論支撐作用。

微生物傳統(tǒng)研究方法主要是菌種分離。隨著技術(shù)的進步，近幾年開始運用分子生物學(xué)方法對微生物進行研究，如變性梯度凝膠電泳技術(shù)、單鏈構(gòu)象多態(tài)性、溫度梯度凝膠電泳技術(shù)、基因文庫法和熒光原位雜交技術(shù)等[9-12]。傳統(tǒng)的菌種分離只能對少數(shù)能分離的菌種進行研究，對低頻突變檢出的敏感性不高，周轉(zhuǎn)時間長、工作量大。相較于此方法，高通量測序技術(shù)在微生物群落結(jié)構(gòu)的研究具有明顯的先進性，作用位點不斷增加，從單個基因分析到多基因分析，具有準(zhǔn)確定量、讀長長、實時檢測等優(yōu)勢[13-14]。

高通量測序技術(shù)已經(jīng)運用于多種分子生物學(xué)領(lǐng)域[15-16]，但在汾酒釀造微生物研究中鮮有報道[17]。本研究運用高通量測序技術(shù)，對汾酒發(fā)酵過程中酒醅樣品中的真菌菌落結(jié)構(gòu)進行分析，以求找到發(fā)酵過程中的優(yōu)勢微生物；再結(jié)合酒醅營養(yǎng)成分、乙醇和主要香味物質(zhì)的變化，分析優(yōu)勢微生物與營養(yǎng)消耗、物質(zhì)生成之間的關(guān)系，建立一套以高通量測序技術(shù)為基準(zhǔn)的探究汾酒微生物的方法，以期更加全面地解析汾酒發(fā)酵過程中真菌群落結(jié)構(gòu)變化規(guī)律。

1 材料與方法

1.1 材料與試劑

酒醅由山西杏花村汾酒廠提供：汾酒發(fā)酵過程酒醅材料，“大楂”指第1次發(fā)酵，“二楂”指第2次發(fā)酵。根據(jù)汾酒酒醅發(fā)酵過程中關(guān)鍵控制點的發(fā)酵時間取樣，大楂：入缸材料（D0d）、4 d材料（D4d）、7 d材料（D7d）、10 d材料（D10d）、15 d材料（D15d）、21 d材料（D21d）、出缸材料（D28d）；二楂：入缸材料（D0e）、4 d材料（D4e）、7 d材料（D7e）、10 d材料（D10e）、15 d材料（D15e）、21 d材料（D21e）、出缸材料（D28e），選取材料為同一生產(chǎn)班組的同一批次材料，取樣時選取地缸中心部位樣品。

Premix ExTaqTMHot Start Version、DL2000 DNA Marker寶生物工程（大連）有限公司；蛋白酶K 德國Merck公司；溶菌酶 美國Sigma公司；瓊脂糖、丙烯酰胺、N,N’-亞甲基雙丙烯酰胺、去離子甲酰胺 美國Solarbio公司；引物由上海英濰捷基生物技術(shù)公司合成。

1.2 儀器與設(shè)備

水平電泳儀、凝膠成像分析系統(tǒng)、Thermal Cycler S1000聚合酶鏈?zhǔn)椒磻?yīng)（polymerase chain reaction，PCR）儀 美國Bio-Rad公司；X1R高速冷凍離心機、酶標(biāo)儀 美國Thermo公司；MX-S型可調(diào)式混勻儀 美國賽洛捷克公司；L96G梯度型PCR儀 杭州朗基科學(xué)儀器有限公司；Illumina MiSeq PE300 美國Illumina公司；7890A-5975B氣相色譜-質(zhì)譜聯(lián)用儀 美國Agilent公司；DB-WAX色譜柱（60 m×0.25 mm，0.25 μm）、2410高效液相色譜儀 美國Waters公司。

1.3 方法

1.3.1 酒醅的理化分析

水分測定：采用熱干燥法，取10 g酒醅在120 ℃干燥至質(zhì)量恒定。

淀粉、乙醇、葡萄糖和麥芽糖測定：采用液相色譜法，準(zhǔn)確稱取50 g酒醅，用100 mL蒸餾水充分浸泡，過濾，取清液，然后色譜檢測，流動相：0.002 mol/L H2SO4溶液；色譜柱：離子交換柱；檢測器：示差檢測器。

代謝產(chǎn)物檢測：采用氣相色譜-質(zhì)譜（gas chromatography-mass spectrometry，GC-MS）聯(lián)用法，準(zhǔn)確稱取50 g酒醅，用100 mL 60%乙醇溶液充分浸泡，過濾，取清液檢測。色譜條件：進樣口溫度為230 ℃，載氣He（99.999%）；恒流1.0 mL/min，分流比50∶1，進樣量1 μL；采用程序升溫：初始溫度60 ℃，保持1 min，以5 ℃/min升至155 ℃，再以3 ℃/min升至180 ℃，保持10 min。質(zhì)譜條件：離子源溫度為230 ℃，電離方式為電子電離源，電子能量為70 eV，掃描質(zhì)量范圍為20～500 u，溶劑延遲3 min。

理化分析利用Excel軟件對數(shù)據(jù)進行整理，然后利用Origin 8作圖。

1.3.2 真菌群落結(jié)構(gòu)分析

酒醅總DNA的提?。菏榛谆寤@法[5]，將提取的DNA加適量含10 ng/μL RNase的100 μL超純水溶解，37 ℃孵育1 h，分裝兩管，置于－20 ℃保存?zhèn)溆谩?/p>

真菌通用ITS1FI2（5’-GAA CCW GCG GAR GGA TCA-3’）和ITS2（5’-GCT GCG TTC TTC ATC GAT GC-3’）用于真菌ITS序列的擴增，序列長度205 bp。PCR擴增體系為：模板DNA 100 ng，引物1和引物2各2.5 μL，Premix ExTaqTMHot Start Version12.5 μL，加ddH2O補至25 μL。PCR程序：98 ℃預(yù)變性10 s后進行33 個循環(huán)，每個循環(huán)包括98 ℃變性10 s、52 ℃退火30 s、72 ℃延伸30 s；最后一次循環(huán)72 ℃延伸10 min，4 ℃保存[18-19]。

高通量測序：應(yīng)用Illumina MiSeq PE300分析平臺，得到的原始圖像數(shù)據(jù)文件，并經(jīng)堿基識別分析轉(zhuǎn)化為原始測序序列[20]。

數(shù)據(jù)分析：對原始測序序列進行過濾、雙端拼接，得到優(yōu)化序列（Tags），將優(yōu)化序列在97%的相似性水平下進行聚類，劃分可操作分類單元（operational taxonomic units，OTU），并根據(jù)OTU的序列組成得到其物種分類[21]。利用Qiime軟件和Mothur軟件[22]進行微生物群落測序數(shù)據(jù)分析，對測序數(shù)據(jù)進行總體分析后，再利用Silva[23]數(shù)據(jù)庫進行16S rRNA基因序列的比對，獲得OTU分類，確定序列對應(yīng)微生物的分類學(xué)地位和多樣性分析，同時利用R語言進行樣品間的主坐標(biāo)分析（principal coordinate analysis，PCoA）及作圖[21,24]，利用Canoco5進行微生物多樣性分析中物種與環(huán)境因子相關(guān)性分析（redundancy analysis/canonical correspondence analysis，RDA/CCA）及作圖[25-26]。

2 結(jié)果與分析

2.1 汾酒酒醅中物質(zhì)的消長規(guī)律

汾酒酒醅入缸發(fā)酵過程是一個在微生物作用下物質(zhì)消長的過程，在這個過程中主要的營養(yǎng)消耗物有淀粉、還原糖、葡萄糖和麥芽糖。酒醅中的淀粉、糊精和糖類是釀酒所需要的主要碳源，酒醅中的淀粉被轉(zhuǎn)化為葡萄糖、麥芽糖等還原糖類物質(zhì)，還原糖類是微生物直接利用的碳源，可以通過糖酵解途徑生成丙酮酸，然后進行發(fā)酵生成乙醇。水分是細(xì)胞的重要組成部分，也是汾酒固態(tài)發(fā)酵酒醅的主要控制指標(biāo)之一；根據(jù)真菌主要作用產(chǎn)生的酒類主要成分有乙醇、乙酸乙酯和高級醇等，乙醇是酒醅發(fā)酵的主要產(chǎn)物，生成量是鑒定酒醅質(zhì)量的主要指標(biāo)，也是判斷發(fā)酵是否正常的重要數(shù)據(jù)，乙酸乙酯是汾酒中的主要香味物質(zhì)之一，汾酒中的高級醇主要包括正丙醇、正戊醇、異丁醇和丁二醇等。檢測汾酒酒醅中的物質(zhì)變化反映發(fā)酵過程的進展，也反映微生物檢測時間點的物質(zhì)基礎(chǔ)。

由圖1可見，發(fā)酵過程中，大楂酒醅在4～10 d水分含量快速上升，后水分含量穩(wěn)定，二楂酒醅0～4 d水分含量增加明顯，后緩慢上升。大楂酒醅的初始淀粉質(zhì)量分?jǐn)?shù)為30%左右，在4～10 d之間淀粉含量迅速下降，10 d后緩慢下降，第28天大楂酒醅淀粉為13%，大楂階段利用淀粉17%左右；二楂酒醅的初始淀粉為13.5%左右，在0～4 d淀粉含量快速下降，后緩慢下降，第28天二楂酒醅淀粉為8%左右，二楂階段利用淀粉5.5%左右。汾酒大楂發(fā)酵過程中，0～4 d酒醅內(nèi)液化、糖化作用明顯，大量的葡萄糖、麥芽糖生成，4～7 d葡萄糖、麥芽糖含量快速下降，說明微生物在快速消耗還原糖，進行發(fā)酵，7 d后，葡萄糖含量緩慢上升，麥芽糖含量較為平穩(wěn)，汾酒二楂發(fā)酵過程中，0～7 d葡萄糖、麥芽糖含量迅速下降，7 d后含量較為穩(wěn)定，說明大楂發(fā)酵的殘余還原糖直接影響了二楂發(fā)酵變化規(guī)律。

圖1 酒醅中主要營養(yǎng)物質(zhì)變化規(guī)律Fig.1 The changes of main nutrients in fermented grains

由圖2可見，在汾酒大楂發(fā)酵過程中，4～10 d進入乙醇主發(fā)酵期，含量迅速增加，后緩慢上升，說明4～10 d為酒醅產(chǎn)生乙醇的最重要時期，汾酒二楂發(fā)酵過程，入缸就進入快速發(fā)酵階段，在0～4 d，乙醇含量迅速上升，后緩慢上升。乙酸乙酯在大楂發(fā)酵中，隨著時間的延長，含量從無到有逐步上升，7 d后上升速率有所增加，第28天時能達(dá)到75 mg/100 g左右，二楂開始階段乙酸乙酯產(chǎn)生速率較快，中期放緩，后期略有下降，第28天時能達(dá)到38 mg/100 g左右。高級醇的變化較為相似，都是在發(fā)酵前期產(chǎn)生速率較快，后期基本趨于平穩(wěn)，在發(fā)酵10 d之后基本趨于穩(wěn)定，二楂酒醅中高級醇含量達(dá)到35 mg/100 g左右，大楂酒醅中高級醇含量達(dá)到30 mg/100 g左右，二楂中的高級醇含量明顯高于大楂，這是二楂酒醅產(chǎn)出的酒質(zhì)不如大楂的重要原因之一。

圖2 酒醅中微生物主要代謝物質(zhì)變化規(guī)律Fig.2 Changes of main microbial metabolites in fermented grains

酒醅中的淀粉包括糊精和糖類，在大楂酒醅發(fā)酵的0～4 d中，淀粉相對穩(wěn)定，糖類（葡萄糖和麥芽糖）有明顯的增加，這可能是淀粉轉(zhuǎn)化為糖的過程，這也是酒醅形成一個營養(yǎng)穩(wěn)態(tài)的時期（有利于后期被利用），而7～10 d酒醅中淀粉、糖類含量迅速下降，這時乙醇含量明顯增加，這是汾酒大楂酒醅乙醇發(fā)酵的主要時期；對于二楂酒醅，由于僅在大楂酒醅蒸餾酒后加入大曲然后進行入缸發(fā)酵，酒醅中的揮發(fā)性物質(zhì)流出，而其余環(huán)境基本不變，酒醅入缸后迅速進入發(fā)酵期，但由于后期營養(yǎng)條件的限制，發(fā)酵放緩。綜合可以看出，大楂相較于二楂進入發(fā)酵期稍晚，造成這樣的原因可能是酒醅條件的差異。同時，大楂比二楂物質(zhì)利用更高，產(chǎn)生乙醇和乙酸乙酯也較高，說明在這些物質(zhì)產(chǎn)生階段主要作用微生物比例也較高。

2.2 酒醅中真菌群落結(jié)構(gòu)分析

14 個酒醅樣品，經(jīng)過過濾、拼接后共產(chǎn)生1 131 292 條Tags，平均每個樣品產(chǎn)生80 852 條Tags，最多129 661 條Tags，最少36 942 條Tags。在97%的相似度水平下對Tags進行OTU劃分后，基于Silva數(shù)據(jù)庫對OTU進行分類學(xué)物種注釋。酒醅樣品Tags聚類共獲得98 個真菌OTU，其中大楂98 個OTU，二楂93 個OTU。采用對測序序列進行隨機抽樣的方法，以抽到的序列數(shù)與它們代表OTU數(shù)目構(gòu)建曲線，即稀釋性曲線[27-28]。當(dāng)曲線趨向平坦時，說明測序數(shù)量合理，更多的數(shù)據(jù)量對于發(fā)現(xiàn)新的OTU貢獻很小，反之則表明繼續(xù)測序還可能產(chǎn)生較多新的OTU[29]。因此，通過作稀釋性曲線（圖3），對測序深度評估。

圖3 稀釋曲線Fig.3 Rarefaction curves

為了在物種水平中更少占比出現(xiàn)在others中，在科分類學(xué)水平上對酒醅樣品真菌物種組成及相對豐度進行統(tǒng)計分析，得到各樣品中真菌科水平的分布情況（圖4），在屬水平對最優(yōu)勢的各屬進行了統(tǒng)計分析（圖5）。酒醅樣本總共檢測到29 個科的真菌，其中，大楂酒醅樣品中Incertae_sedis、Saccharomycetaceae、Mucoraceae與Lichtheimiaceae是優(yōu)勢菌科，共占已知菌科的90.0%。二楂酒醅樣品中Incertae_sedis、Pichiaceae、Saccharomycetaceae、Rhizopodaceae與Lichtheimiaceae是優(yōu)勢菌科，共占已知菌科的94.6%。而相對豐度最高的Incertae_sedis在大楂發(fā)酵時期，從0 d占比11.0%，至7 d上升至88.8%，后下降至50%左右，在二楂發(fā)酵時期，從0 d占比29.7%，至7 d上升至90.7%，后穩(wěn)定在90%以上（除D21e為60%左右），說明該科真菌在二楂發(fā)酵時的生長繁殖速度明顯高于大楂發(fā)酵時期。從二楂的樣品中可以看出Candida隨著發(fā)酵的推進，占比逐漸增加，到第15天樣品的時期達(dá)到最高（94%），在后期發(fā)酵中略有下降；Pichia在發(fā)酵開始的樣品中比例達(dá)到60%以上，隨著發(fā)酵的推進，占比下降。從大楂樣品中可以看出，Candida和酵母屬（Saccharomycessp.）的占比呈現(xiàn)先上升后下降最后趨于穩(wěn)定的態(tài)勢，Candida在7 d占比最高，達(dá)到87%，Saccharomycessp.在4 d占比最高，達(dá)到41%，Pichia在大楂樣品中隨著發(fā)酵的推進占比逐漸下降，釀酒酵母（S.cerevisiae）在大楂樣品中隨著發(fā)酵的推進占比呈現(xiàn)先上升后下降的趨勢。說明大楂酒醅發(fā)酵過程的主要優(yōu)勢菌為Candida、Saccharomycessp.、Pichia和S.cerevisiae，二楂酒醅發(fā)酵過程的主要優(yōu)勢菌為Candida和Pichia。

圖4 大楂（A）和二楂（B）真菌結(jié)構(gòu)及分布規(guī)律Fig.4 Fungal community structure and distribution in primary (A) and secondary (B) fermentation

圖5 酵母的分布規(guī)律Fig.5 Yeast community structure and distribution

酵母是所有樣品中的主要微生物，最少的70%，也是兩個發(fā)酵過程差異性最大的微生物群體。發(fā)酵初期大楂樣品和二楂樣品的酵母結(jié)構(gòu)相似，都是以Pichia為主，Candida其次，S.cerevisiae較少（Saccharomycetaceae包括S.cerevisiae和漢遜酵母（Hanseniaspora），而S.cerevisiae占比更高）；隨著發(fā)酵的推進，大楂樣品和二楂樣品的酵母結(jié)構(gòu)出現(xiàn)了較大的差異性，隨著發(fā)酵的推進大楂樣品中的Candida和S.cerevisiae均有所上升，在大楂發(fā)酵后期S.cerevisiae比例不低于20%，而二楂樣品中Candida是主要優(yōu)勢菌，占比均超過60%（除21 d外均達(dá)到80%以上），S.cerevisiae比例最高的28 d也僅6%左右。綜上可以看出大楂樣品和二楂樣品的酵母結(jié)構(gòu)差異主要是Candida和S.cerevisiae的比例，這也可能是兩楂生產(chǎn)產(chǎn)酒率差異的主要原因。

2.3 酒醅樣品中的物質(zhì)成分與微生物差異性和關(guān)聯(lián)性分析

從酒醅樣品中真菌群落結(jié)構(gòu)出發(fā)，利用PCoA探究樣品間存在的差異性；酒醅中的淀粉和糖為微生物生長和代謝提供了營養(yǎng)物質(zhì)，微生物的代謝產(chǎn)物又影響其生存環(huán)境[30]，采用RDA/CCA分析探究酒醅中物質(zhì)成分和微生物的關(guān)系。

如圖6所示，大楂和二楂酒醅樣品的真菌群落結(jié)構(gòu)在發(fā)酵開始時（D0d和D0e）相似性高，隨著發(fā)酵的推進到7 d以后大楂和二楂酒醅樣品中的真菌群落結(jié)構(gòu)在發(fā)酵開始出現(xiàn)較大的差異性。發(fā)酵過程中二楂酒醅樣品的真菌群落結(jié)構(gòu)與大楂7 d樣品（D7d）真菌群落結(jié)構(gòu)較為相似，大楂酒醅從發(fā)酵第10天（D10d）以后的樣品真菌群落較為相似。如圖7所示，大楂酒醅樣品的真菌群落結(jié)構(gòu)，4、10、21 d材料相似性高，0、7、15 d材料相似性低。而二楂酒醅樣品的真菌群落結(jié)構(gòu)，7、10、15、28 d材料相似性高，0、4、21 d材料相似性低。也就是說，汾酒發(fā)酵過程中，大楂和二楂真菌群落在發(fā)酵中后期有較大的差異性，其各自真菌群落在發(fā)酵中后期相對較為穩(wěn)定。

圖6 酒醅樣品的PCoAFig.6 PCoA plot for similarity in fungal community structure in fermented grains

將群落結(jié)構(gòu)信息作決策曲線分析，軸中梯度最長為2.43，大于1，小于3，選擇對酒醅真菌菌落和酒醅中物質(zhì)成分作RDA（圖7），其分析結(jié)果蒙特卡洛檢驗值P值為0.10，PC1為56.7%，PC2為36.3%，二維成分占比達(dá)到93%。8 個物質(zhì)因子（水分、葡萄糖、淀粉、麥芽糖、乙醇、高級醇、乙酸乙酯和乳酸乙酯）全部用于進行關(guān)聯(lián)性分析。

圖7 酒醅中物質(zhì)成分與真菌群落的關(guān)聯(lián)性Fig.7 RDA analysis of correlation between fungal community structure and nutrients in fermented grains

通過酒醅中物質(zhì)成分與真菌群落主要微生物的RDA結(jié)果可以看出，RDA樣品真菌菌落的聚類與PCoA一致。3 個影響最大的微生物是S.cerevisiae、Pichia、Candida，影響大楂酒醅發(fā)酵后期最主要微生物是S.cerevisiae，影響二楂酒醅發(fā)酵最主要微生物是Candida，Pichia是發(fā)酵初期酒醅最主要微生物。發(fā)酵初期和大楂發(fā)酵后期酒醅中的營養(yǎng)物質(zhì)（淀粉、葡萄糖和麥芽糖）和真菌群落結(jié)構(gòu)呈正相關(guān)，說明營養(yǎng)越充分有利于該階段真菌群落生長。而二楂酒醅中真菌群落與營養(yǎng)物質(zhì)呈負(fù)相關(guān)，與高級醇和水分等物質(zhì)呈正相關(guān)，這種相關(guān)性說明二楂酒醅中的物質(zhì)成分（營養(yǎng)的缺乏）決定了其真菌群落結(jié)構(gòu)與大楂的差異性。高級醇和乳酸乙酯含量與Candida有較強的正相關(guān)，乙酸乙酯和乙醇含量與S.cerevisiae呈正相關(guān)。

3 結(jié)論與討論

本研究采用Illumina MiSeq測序方法，直接從汾酒發(fā)酵酒醅樣品中擴增高變區(qū)進行測序，避免了采用培養(yǎng)方式研究微生物的多種局限性，客觀還原了汾酒發(fā)酵過程中菌群結(jié)構(gòu)及豐度的變化規(guī)律，為汾酒微生物應(yīng)用研究提供了理論支撐。結(jié)果表明：在97%的相似度水平下，一共獲得98 個真菌OTU；在科的分類水平下，大楂酒醅樣品中Incertae_sedis、Saccharomycetaceae、Mucoraceae與Lichtheimiaceae是優(yōu)勢菌科。二楂酒醅樣品中Incertae_sedis、Pichiaceae、Saccharomycetaceae、Rhizopodaceae與Lichtheimiaceae是優(yōu)勢菌科，三大酵母菌屬（S.cerevisiae、Pichia、Candida）為發(fā)酵過程真菌的優(yōu)勢菌屬；大楂酒醅樣品的真菌多樣性指數(shù)呈現(xiàn)先上升后下降的規(guī)律，而二楂酒醅樣品的真菌多樣性指數(shù)在入缸發(fā)酵時最高，隨著發(fā)酵的進行，多樣性指數(shù)逐步降低；大、二楂樣品真菌群落結(jié)構(gòu)發(fā)酵后期存在顯著差異；從物質(zhì)的消長角度可以看出，大楂酒醅產(chǎn)生乙醇的時期在4～10 d之間，產(chǎn)生酯類從4 d左右到發(fā)酵結(jié)束，而二楂從發(fā)酵開始就進入乙醇和酯類產(chǎn)生時期，在發(fā)酵10 d左右物質(zhì)產(chǎn)生放緩。汾酒的兩楂酒醅發(fā)酵過程真菌群落呈現(xiàn)較大差異性，結(jié)合酒醅中物質(zhì)的變化和RDA可以推測這種差異性主要有兩方面的原因：

一是營養(yǎng)物質(zhì)的差異性。大楂在發(fā)酵過程消耗淀粉大約為17%，而二楂酒醅發(fā)酵過程僅消耗約5.5%，從糖類的消耗也可以看出大楂酒醅中的葡萄糖和麥芽糖含量在各個時期均高于二楂酒醅，同時大楂酒醅中的葡萄糖含量在整個發(fā)酵過程均高于1%，而二楂酒醅在發(fā)酵15 d左右葡萄糖就處于最低的水平（0.25%），后期可利用的葡糖糖很少，這些都說明二楂酒醅中的營養(yǎng)條件比酒醅匱乏，二楂中最主要的微生物是Candida（最少為D21e號樣，達(dá)到61%），在營養(yǎng)匱乏的情況下以芽生孢子出芽繁殖，孢子伸長形成芽管，不與母體分離，形成較長的假菌絲，能夠克服酒醅營養(yǎng)匱乏，這可能是該類酵母在二楂中高于其他酵母含量的重要原因。Saccharomycetaceae包括S.cerevisiae、Hanseniaspora、有孢圓酵母（Torulaspora）等，該類菌種在營養(yǎng)條件充足情況下有較強的產(chǎn)乙醇和產(chǎn)酯能力，所以該類菌在大楂酒醅中有較高的比例，這有利于酒醅產(chǎn)乙醇和產(chǎn)酯，據(jù)報道Hanseniaspora是清香型白酒中乙酸乙酯能力最強的菌種[31]，這可能是大楂乙酸乙酯明顯高于二楂的重要原因。所以酒醅中營養(yǎng)物質(zhì)的動態(tài)變化不僅間接反映乙醇的生成情況和發(fā)酵狀況，而且可以特征反映發(fā)酵過程中微生物的生命活動狀況，也就是酒醅的環(huán)境影響了微生物的群落結(jié)構(gòu)。

二是酒醅環(huán)境條件的差異性。在二楂酒醅中酸度高于大楂酒醅，Candida有較強的耐受性，能在極端環(huán)境進行生長（微厭氧和高酸度），而Pichia的耐酸性和耐溫能力較低，隨著發(fā)酵的進行，酸度和溫度的升高導(dǎo)致Pichia的比例下降，同時由于霉菌的厭氧的耐受性較差，所以在所有的酒醅中霉菌的比例較低。這種營養(yǎng)條件和環(huán)境條件的差異性可以下一步進行深入研究，結(jié)合微生物變化和工藝條件，可能為汾酒的發(fā)酵生產(chǎn)提供更有價值的信息。

大楂酒醅和二楂酒醅入缸后進入乙醇主發(fā)酵期有一定的時間差異性，大楂是在發(fā)酵4 d后進入乙醇的主發(fā)酵期，二楂是在入缸后進入乙醇的主發(fā)酵期，這種差異性可能主要是微生物和酒醅環(huán)境的影響，大楂酒醅入缸后由于是糊化后的糧食和曲混合，酒醅比較疏松，氧氣含量充足，糖化速度快，底物濃度高，水分、溫度、酸度比較適宜，酵母細(xì)胞出芽、裂殖成倍增加，酵母菌迅速達(dá)到峰值，通過該時期后迅速進入乙醇的主發(fā)酵期，而二楂酒醅是在大楂酒醅蒸餾酒后伴曲直接發(fā)酵，酒醅的黏稠度、酸度、水分和大楂初始發(fā)酵酒醅完全不同，此環(huán)境是對接入的大曲微生物進行篩選，同時環(huán)境的還原糖含量較高，有利于微生物利用，酒醅迅速進入乙醇的主發(fā)酵期，但是由于營養(yǎng)物質(zhì)的缺乏，發(fā)酵10 d基本趨于平穩(wěn)，生成的乙醇量僅為大楂酒醅的一半。在乙醇和高級醇的產(chǎn)生時期（大楂4～15 d，二楂0～10 d），Candida含量較高，該菌屬鮮見有高產(chǎn)醇類的報道，關(guān)于Candida在汾酒酒醅主發(fā)酵期的作用是后期的一個研究點。

本研究解析了汾酒兩楂酒醅發(fā)酵過程的真菌群落結(jié)構(gòu)，找到了其差異性和差異的優(yōu)勢真菌屬，同時推測了真菌群落結(jié)構(gòu)差異原因和酒醅乙醇主發(fā)酵期差異性的原因，為汾酒的生產(chǎn)工藝優(yōu)化和提高生產(chǎn)質(zhì)量提供基礎(chǔ)數(shù)據(jù)支撐。