唐佳代,劉力萍,龍亞飛,趙益梅,趙騰飛,王芙蓉,李 喆
(1.茅臺(tái)學(xué)院 釀酒工程系,貴州 仁懷 564500;2.貴州茅臺(tái)酒股份有限公司,貴州 仁懷 564500)
醬香大曲是以小麥、水和母曲作為生產(chǎn)原料,經(jīng)過(guò)高溫發(fā)酵和貯存而成,用于生產(chǎn)大曲醬香型白酒的糖化發(fā)酵劑[1]。其生產(chǎn)工藝有“制曲黑白黃”一說(shuō),即成品曲按照色澤不同分為黃曲、黑曲和白曲,其外觀上具有黃曲皮厚,黑曲、白曲皮薄的特點(diǎn)[2]。黃曲醬香濃郁,糖化力較低;黑曲發(fā)酵品溫高常帶有糊味;白曲醬香弱,但其糖化率高[3]。在醬香大曲生產(chǎn)中,黃曲率常作為成品曲塊質(zhì)量分析重要指標(biāo)之一。通過(guò)控制制曲工藝提高黃曲率的研究較普遍,而黃曲與黑、白曲間微生物差異并不明晰,不利于應(yīng)用微生物技術(shù)對(duì)酒曲進(jìn)行改良。
大曲微生物組成結(jié)構(gòu)及代謝決定了大曲的品質(zhì),醬香大曲中微生物有細(xì)菌、酵母、霉菌以及放線菌。胡寶東等[3]利用傳統(tǒng)分離培養(yǎng)方法對(duì)黑、白、黃醬香大曲中微生物計(jì)數(shù),細(xì)菌總數(shù)較高,達(dá)到107數(shù)量級(jí)[3]。細(xì)菌在醬香大曲發(fā)酵初期提供釀造所需的酸性環(huán)境,產(chǎn)生蛋白酶、脂肪酶、糖化酶、纖維素酶等豐富多樣的酶系,且具有酯化、促美拉德反應(yīng)和產(chǎn)生吡嗪類物質(zhì)等能力,從而促進(jìn)酒中多種風(fēng)味物質(zhì)的形成[4-8]。高通量測(cè)序技術(shù)具備全面的表現(xiàn)樣品微生物豐度等優(yōu)勢(shì),近年來(lái)普遍被應(yīng)用于釀酒領(lǐng)域中不同香型大曲、酒醅、酒糟的微生物群落結(jié)構(gòu)分析,基于該技術(shù)分析醬香大曲中黑、白、黃曲微生物群落結(jié)構(gòu)的報(bào)道較少[9-12]。本研究應(yīng)用高通量測(cè)序分析不同感官特性醬香大曲中細(xì)菌群落結(jié)構(gòu),結(jié)合3種大曲理化指標(biāo),分析其微生物群落結(jié)構(gòu)與理化、發(fā)酵特性的關(guān)聯(lián)性。同時(shí),掌握更多醬香型白酒釀造過(guò)程中的細(xì)菌資源,為探索利用曲源微生物提升大曲的醬香和糖化力等性能奠定研究基礎(chǔ)。
黃曲(YQ)、黑曲(BQ)、白曲(WQ):貴州省仁懷市某酒廠出倉(cāng)曲。樣品采集時(shí)間為2020年11月,取樣后分別裝入密封袋,運(yùn)回實(shí)驗(yàn)室保藏,備用。
氫氧化鈉(分析純):天津大茂化學(xué)試劑廠;乙酸(分析純):川東化工集團(tuán);乙酸鈉(分析純):上海阿拉丁生化科技股份有限公司;脫氧核糖核酸(deoxyribonucleic acid,DNA)提取試劑盒:美國(guó)Zymo Research公司;KOD DNA聚合酶(KOD-401B)(1 U/μL):東洋紡(上海)生物科技有限公司;脫氧核糖核苷三磷酸(deoxy-ribonucleoside triphosphate,dNTP)混合液:東洋紡(上海)生物科技有限公司。
FA1204N電子天平:上海菁海儀器有限公司;DHG-9140A電熱恒溫鼓風(fēng)干燥箱、HH-8J數(shù)顯恒溫水浴鍋:常州潤(rùn)滑電器有限公司;DK-98-Ⅱ電爐:天津市泰斯特儀器有限公司;YXQ-LS-100SII立式壓力蒸汽滅菌鍋:上海博訊醫(yī)療股份有限公司;BCD-640WAGM冷藏柜:青島海爾股份有限公司。
1.3.1 不同感官特性醬香大曲理化特性分析
按照輕工行業(yè)標(biāo)準(zhǔn)QB/T 4257—2011《釀酒大曲通用分析方法》[13]對(duì)三種不同感官特性醬香大曲(黃曲、黑曲和白曲)的水分含量、酸度和液化力進(jìn)行測(cè)定。
1.3.2 醬香型大曲中DNA提取、PCR擴(kuò)增和測(cè)序
將醬香大曲樣品送至成都羅寧生物科技有限公司,使用BIOMICS DNA提取試劑盒提取黃曲、黑曲和白曲樣品中的DNA。在細(xì)菌16S rRNA V4區(qū)域進(jìn)行擴(kuò)增,擴(kuò)增引物F(5'-GTGYCAGCMGCCGCGGTAA-3')和R(5'-GGACTACHVGGGTWTCTAAT-3')。聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction,PCR)擴(kuò)增用酶為KOD PLUS Neo(KOD-401B),PCR擴(kuò)增體系(50 μL):5 μL 10×PCR Buffer for KODPlus-Neo,5 μL 2 mmol/L dNTP,5 μL 25 mmol/L MgSO4,1.5 μL上游引物U515F,1.5 μL下游引物U806R,5 μL DNA模板,31 μL超純水,1 μL KOD-Plus-Neo。PCR擴(kuò)增條件為:94 ℃預(yù)變性1 min;94 ℃變性2 s;54 ℃退火30 s;72 ℃延伸30 s,25~30個(gè)循環(huán),72 ℃延伸5 min;4 ℃保溫。擴(kuò)增產(chǎn)物純化后用2%的瓊脂糖凝膠,每個(gè)樣本進(jìn)行3次重復(fù),取線性期PCR產(chǎn)物等量混合后建庫(kù)。進(jìn)一步在Illumina HiSeq 2500測(cè)序平臺(tái)開(kāi)展高通量測(cè)序。
1.3.3 高通量測(cè)序數(shù)據(jù)處理與分析
高通量測(cè)序所得序列通過(guò)拼接雙端序列、過(guò)濾去除低質(zhì)量序列、使用Uchime算法和Gold數(shù)據(jù)庫(kù)去除嵌合體,進(jìn)而得到有效序列(Effiective Tag)?;赨search軟件,使用UPARSE算法在97%的一致性水平上進(jìn)行操作分類單元(operational taxonomic unit,OTU)聚類。并基于UCLUST分類法與Silva分類學(xué)數(shù)據(jù)庫(kù)進(jìn)行16S rRNA基因序列比對(duì)和注釋分析。使用QIIME(V 1.9.0)軟件對(duì)樣品Alpha多樣性指數(shù)(Alpha diversity)進(jìn)行分析,用于衡量Alpha多樣性指標(biāo)主要有Chao1、ACE、Shannon和Simpson指數(shù)。其中Chao1和ACE指數(shù)用于衡量樣本中物種的豐度,兩者數(shù)值與樣本物種總數(shù)成正比[14]:Chao1主要依據(jù)有單條或雙條序列OTU數(shù)目來(lái)判斷物種豐度;ACE指數(shù)依據(jù)有單條序列OTU數(shù)目與出現(xiàn)頻次低于10次的物種來(lái)衡量樣本中物種總數(shù)。Shannon和Simpson指數(shù)用于綜合衡量樣本物種多樣性[14]:Shannon指數(shù)表示樣本中各物種和物種間在樣本中分布的均勻性,Shannon數(shù)值與物種均勻性成正比;Simpson指數(shù)衡量樣本中物種多樣性,Simpson指數(shù)與樣本多樣性成正比。
不同感官特性醬香大曲的感官特性差異見(jiàn)表1。由表1可知,黃曲呈金黃色或棕黃色,曲香濃郁;黑曲呈黑色或棕褐色,帶有糊味;白曲呈麥粉色,有生麥味。制曲過(guò)程中,小麥磨碎度、拌料水分適中,入倉(cāng)堆積處于發(fā)酵倉(cāng)中層,溶氧較好易形成黃曲。若入倉(cāng)堆積的曲坯處于發(fā)酵倉(cāng)中下層,溶氧少,濕氣重,水分大,曲坯發(fā)酵過(guò)程中熱曲時(shí)間長(zhǎng)則易形成黑曲。當(dāng)入倉(cāng)堆積時(shí),曲坯處于發(fā)酵倉(cāng)上層,門窗密封性不好,水分揮發(fā)快,熱曲時(shí)間短則易形成白曲。
表1 不同感官特性的醬香大曲Table 1 Sauce-flavor Daqu with different sensory characteristics
醬香大曲的理化指標(biāo)是衡量大曲成熟與否的重要標(biāo)志,在一定程度上反映大曲品質(zhì)。3種不同感官特性醬香大曲的理化指標(biāo)結(jié)果見(jiàn)表2。
表2 不同感官特性醬香大曲的理化指標(biāo)Table 2 Physicochemical indexes of sauce-flavor Daqu with different sensory characteristics
水分含量過(guò)高會(huì)導(dǎo)致大曲生霉,從而導(dǎo)致大曲的質(zhì)量降低。由表2可知,所測(cè)不同感官特性醬香大曲樣品中水分含量由高到低依次為白曲>黑曲>黃曲,水分含量分別為白曲12.89%、黑曲12.44%、黃曲11.63%,均低于13%,符合地方標(biāo)準(zhǔn)DB52/T 871—2014《醬香型白酒釀酒用大曲》要求[15]。
酸度是醬香大曲的關(guān)鍵性指標(biāo)之一,是微生物綜合作用的結(jié)果,酸度適當(dāng)可一定程度上抑制雜菌,還可為有益微生物提供適宜的生長(zhǎng)繁殖環(huán)境。不同感官特性大曲樣品的酸度由高到低依次為黑曲>白曲>黃曲,分別為黑曲0.31 mmol/g、白曲0.262 mmol/g、黃曲0.252 mmol/g。醬香大曲酸度范圍應(yīng)在0.10~0.35 mmol/g,3種醬香大曲樣品均符合地方標(biāo)準(zhǔn)DB52/T 871—2014《醬香型白酒釀酒用大曲》要求[15]。
液化力是衡量大曲發(fā)酵性能高低的重要指標(biāo),液化力高將為產(chǎn)酒、產(chǎn)香奠定物質(zhì)基礎(chǔ)。不同感官特性大曲樣品液化力由高到低依次為黃曲>白曲>黑曲,分別為黃曲1.23 g/(g·h)、白曲0.98 g/(g·h)、黑曲0.84 g/(g·h)。
醬香型大曲樣品提取基因組DNA后,對(duì)V4區(qū)進(jìn)行測(cè)序,進(jìn)行雙端拼接獲得到優(yōu)質(zhì)序列,優(yōu)質(zhì)序列再經(jīng)過(guò)濾嵌合體后共獲得308 645條有效序列,平均長(zhǎng)度為297 bp。
韋恩圖直觀展現(xiàn)了不同感官特性醬香大曲樣品中的共有OTU和特有OTU[12]。不同感官特性醬香大曲細(xì)菌韋恩圖見(jiàn)圖1。由圖1可知,三種樣品共同的OTU有521個(gè),黃曲樣品特有的OTU為370個(gè),白曲樣品特有OTU為146個(gè),黑曲樣品特有OTU為970個(gè)。結(jié)果表明黃曲YQ、黑曲BQ和白曲WQ樣品中的OTU組成既有相似性又有特異性。
圖1 不同感官特性醬香大曲細(xì)菌韋恩圖Fig.1 Venn diagram of bacteria of sauce-flavor Daqu with different sensory characteristics
OTU數(shù)稀釋曲線是以樣本中隨機(jī)抽取一定數(shù)量序列數(shù)與序列所代表的物種數(shù)繪制的曲線。該曲線反映了樣本物種豐富度隨測(cè)序深度的變化趨勢(shì),以及持續(xù)抽樣下新物種出現(xiàn)的速率。不同感官特性醬香大曲OTU數(shù)稀釋曲線見(jiàn)圖2。由圖2可知,YQ、WQ和BQ樣品的稀釋曲線前期急劇上升,隨著抽取序列的增加后趨于平緩。表明在抽取序列前期隨著抽取序列數(shù)增加,有大量細(xì)菌物種被檢出;當(dāng)抽取檢測(cè)大于20 000個(gè)序列后,細(xì)菌物種不再隨測(cè)序數(shù)量的增加而顯著增多,曲線趨于平緩,表明樣本測(cè)序量充分,測(cè)序結(jié)果能真實(shí)反映不同感官特性醬香大曲的細(xì)菌多樣性。
圖2 不同感官特性醬香大曲細(xì)菌OTU數(shù)稀釋曲線Fig.2 OTU numbers rarefaction curves of bacteria of sauce-flavor Daqu with different sensory characteristics
Coverage數(shù)值表示樣本序列覆蓋率:Coverage數(shù)值數(shù)值與樣本中物種被測(cè)出的概率成正比,各樣本Coverage數(shù)值均大于0.99,認(rèn)為測(cè)序結(jié)果能夠真實(shí)的反映樣品物種豐度及物種多樣性[16]。在97%相似度水平下,大曲樣品Alpha多樣性指數(shù)值統(tǒng)計(jì)見(jiàn)表3。由表3可知,三種大曲樣品中黑曲BQ樣品的Chao1、Shannon和Ace指數(shù)均最高,說(shuō)明黑曲BQ中物種分布較為均勻且檢測(cè)出物種總數(shù)最高;Simpson指數(shù)與黃曲YQ接近,說(shuō)明兩者細(xì)菌的多樣性差異較小。細(xì)菌群落多樣性由高至低依次為:黑曲BQ>黃曲YQ>白曲WQ,細(xì)菌總數(shù)由高至低依次為:黑曲BQ>黃曲YQ>白曲WQ。
表3 不同感官特性醬香大曲細(xì)菌Alpha多樣性指數(shù)Table 3 Bacteria Alpha diversity indexes of sauce-flavor Daqu with different sensory characteristics
黃、黑、白曲屬水平上的細(xì)菌群落結(jié)構(gòu)分析結(jié)果見(jiàn)圖3。由圖3可知,醬香大曲樣品中相對(duì)豐度前十的物種有巖石芽孢桿菌屬(Scopulibacillus)、糖多孢菌屬(Saccharopolyspora)、Desmospora、Kroppenstedtia、乳酸桿菌屬(Lactobacillus)、海芽孢桿菌屬(Pontibacillus)、羅爾斯頓菌屬(Ralstonia)、枝芽孢桿菌屬(Virgibacillus)、大洋芽孢桿菌屬(Oceanobacillus)、芽孢桿菌屬(Bacillus)。
圖3 基于屬水平不同感官特性醬香大曲細(xì)菌群落結(jié)構(gòu)Fig.3 Bacterial community structure of sauce-flavor Daqu with different sensory characteristics based on genus level
白曲WQ中相對(duì)豐度>1%的包括大洋芽孢桿菌屬(Oceanobacillus)(36.16%)、枝芽孢桿菌屬(Virgibacillus)(27.67%)、芽孢桿菌屬(Bacillus)(21.47%)、Desmospora(5.45%)、乳酸桿菌屬(Lactobacillus)(1.64%)、羅爾斯頓菌屬(Ralstonia)(1.56%)。
黃曲YQ中相對(duì)豐度>1%有芽孢桿菌屬(Bacillus)(19.48%)、海芽孢桿菌屬(Pontibacillus)(14.40%)、乳酸桿菌屬(Lactobacillus)(14.02%)、羅爾斯頓菌屬(Ralstonia)(11.43%)、Kroppenstedtia(8.90%)、枝芽孢桿菌屬(Virgibacillus)(7.40%)、大洋芽孢桿菌屬(Oceanobacillus)(4.32%)、Desmospora(3.12%)、巖石芽孢桿菌屬(Scopulibacillus)(1.56%)、醋菌屬(Acetobacter)(1.15%)、糖多孢菌屬(Saccharopolyspora)(1.11%)。
黑曲BQ中相對(duì)豐度>1%有羅爾斯頓菌屬(Ralstonia)(22.21%)、枝芽孢桿菌屬(Virgibacillus)(11.91%)、芽孢桿菌屬(Bacillus)(10.47%)、乳酸桿菌屬(Lactobacillus)(6.39%)、糖多孢菌屬(Saccharopolyspora)(5.39%)、Desmospora(3.03%)、海芽孢桿菌屬(Pontibacillus)(2.34%)、巖石芽孢桿菌屬(Scopulibacillus)(1.68%)、根瘤桿菌屬(Rhizobacter)(1.37%)。
主成分分析(principal component analysis,PCA)通過(guò)將測(cè)序多維數(shù)據(jù)進(jìn)行降維,直觀反映所獲得主成分得分及變量得分,從而實(shí)現(xiàn)樣品細(xì)菌多樣性相似性可視化[9]。不同感官特性醬香大曲組間PCA分析結(jié)果見(jiàn)圖4。
圖4 不同感官特性醬香大曲樣品細(xì)菌屬水平群落主成分分析結(jié)果Fig.4 Principal component analysis results of the bacterial genus level community of sauce-flavor Daqu samples with different sensory characteristics
由圖4可知,兩個(gè)主要成分PC1、PC2分別占49.8%,26.4%,共能解釋76.2%的物種差異。黃曲YQ與白曲WQ、黑曲BQ樣品在橫縱坐標(biāo)上均有明顯距離,說(shuō)明黃曲YQ與白曲WQ、黑曲BQ樣品間細(xì)菌群落結(jié)構(gòu)在屬水平上有較大差異;白曲WQ與黑曲BQ樣品在第一主成分上距離較小,說(shuō)明白曲WQ與黑曲BQ樣品間細(xì)菌群落結(jié)構(gòu)較為相似。
線性判別分析效應(yīng)量(linear discriminant analysis effect size,LEfse)可以進(jìn)行組間類群的差異檢驗(yàn),找到在組間差異顯著的物種,也就是對(duì)組間差異貢獻(xiàn)大的類群。LEfSe的統(tǒng)計(jì)結(jié)果的顯著差異物種線性判別分析(linear discriminant analysis,LDA)值分布柱狀圖,用以展示每個(gè)組內(nèi)顯著富集的物種及其重要性程度。不同感官特性醬香大曲樣品細(xì)菌顯著差異物種線性判別分析見(jiàn)圖5。
圖5 不同感官特性醬香大曲樣品細(xì)菌顯著差異物種線性判別分析Fig.5 Linear discriminant analysis of bacteria significantly different species of sauce-flavor Daqu samples with different sensory characteristics
由圖5可知,從不同感官特性醬香大曲中檢出28個(gè)細(xì)菌類群,其中屬于高相對(duì)豐度菌屬的有:黃曲YQ中的海芽孢桿菌屬(Pontibacillus),黑曲BQ中的羅爾斯頓菌屬(Ralstonia)、糖多孢菌屬(Saccharopolyspora),白曲WQ中的大洋芽孢桿菌屬(Oceanobacillus)、枝芽孢桿菌屬(Virgibacillus)。造成不同感官特性醬香大曲細(xì)菌群落結(jié)構(gòu)差異的主要是海芽孢桿菌屬(Pontibacillus)、羅爾斯頓菌屬(Ralstonia)、糖多孢菌屬(Saccharopolyspora)、大洋芽胞桿菌屬(Oceanobacillus)、枝芽孢桿菌屬(Virgibacillus)和一些相對(duì)豐度較低的菌屬。
高通量測(cè)序結(jié)果顯示大曲樣品在Silva分類學(xué)數(shù)據(jù)庫(kù)比對(duì),鑒別出不同感官特性醬香大曲中主導(dǎo)細(xì)菌屬分別是:白曲WQ中為海洋芽孢桿菌屬(Oceanobacillus)、枝芽孢桿菌屬(Virgibacillus)和芽孢桿菌屬(Bacillus);黃曲YQ中為芽孢桿菌屬(Bacillus)、海芽胞桿菌屬(Pontibacillus)和乳酸桿菌屬(Lactobacillus);黑曲BQ中為羅爾斯頓菌屬(Ralstonia)、枝芽孢桿菌屬(Virgibacillus)和芽孢桿菌屬(Bacillus)。黃曲YQ主導(dǎo)細(xì)菌屬與白曲WQ、黑曲BQ在種類和豐度有明顯差異,白曲WQ和黑曲BQ主導(dǎo)細(xì)菌屬均包括枝芽孢桿菌屬(Virgibacillus)和芽孢桿菌屬(Bacillus)。結(jié)合主成分分析結(jié)果,黃曲YQ與白曲WQ、黑曲BQ中細(xì)菌多樣性有較大差異性,而白曲WQ和黑曲BQ的細(xì)菌菌屬多樣性較為相似。芽孢桿菌屬(Bacillus)是大曲中常見(jiàn)的功能菌,可產(chǎn)生中性、酸性蛋白酶將原料中蛋白質(zhì)降解為氨基酸和肽類,既為微生物生長(zhǎng)提供營(yíng)養(yǎng)物質(zhì),也為美拉德反應(yīng)提供原料[8]。芽孢桿菌屬(Bacillus)被證實(shí)具有產(chǎn)吡嗪類風(fēng)味物質(zhì)能力,特別是甲基類吡嗪化合物,被認(rèn)為是醬香型白酒中的關(guān)鍵風(fēng)味物質(zhì)之一[19-21]。黃曲YQ中主導(dǎo)細(xì)菌屬為芽孢桿菌屬(Bacillus),其可能與黃曲YQ液化力較高,醬香濃郁相關(guān)。有很多研究將大曲中的芽孢桿菌應(yīng)用于制備醬香強(qiáng)化曲[22-24]。此外,芽孢桿菌屬(Bacillus)可產(chǎn)生有機(jī)酸、纖維素酶等,適量的有機(jī)酸可使酒曲pH達(dá)到微酸性為酵母和霉菌創(chuàng)造生存環(huán)境,纖維素酶可提高出酒率、增加酯含量[25-27]。
為探究三種不同感官特性醬香大曲中細(xì)菌群落結(jié)構(gòu)差異,應(yīng)用LEfse分析進(jìn)行組間類群的差異檢驗(yàn),找到組間差異顯著的物種。分析結(jié)果顯示,造成細(xì)菌群落結(jié)構(gòu)差異的主要是Pontibacillus、羅爾斯頓菌屬(Ralstonia)、糖多孢菌屬(Saccharopolyspora)、海洋芽胞桿菌屬(Oceanobacillus)、枝芽孢桿菌屬(Virgibacillus)和一些豐度較低的菌屬。黃曲YQ中造成細(xì)菌群落結(jié)構(gòu)差異的菌屬是Pontibacillus,該菌屬曾被證明具有產(chǎn)生α-淀粉酶能力[28],α-淀粉酶可將淀粉降解為糊精,使淀粉液化,故推測(cè)其與黃曲有較好液化力有一定程度相關(guān)性。
醬香出倉(cāng)大曲表現(xiàn)為不同感特性,該研究通過(guò)理化分析表明,不同感官特性醬香大曲樣品中水分含量、酸度均符合醬香成品的相關(guān)標(biāo)準(zhǔn)。黃曲YQ醬香濃郁且液化力較高,達(dá)到1.23 g/(g·h)。醬香濃郁可能與黃曲YQ中芽孢桿菌屬(Bacillus)為主導(dǎo)細(xì)菌屬相關(guān),而導(dǎo)致黃曲與黑曲、白曲間的液化力差別可能與其差異貢獻(xiàn)較大海芽孢桿菌屬(Pontibacillus)相關(guān)。
基于高通量測(cè)序技術(shù)結(jié)合理化指標(biāo)分析,初步探究了不同感官特性醬香大曲理化特性差異與細(xì)菌多樣性與群落結(jié)構(gòu)差異間的關(guān)聯(lián)性。進(jìn)一步將完善不同感官特性醬香大曲中微生物組學(xué)與理化特性相關(guān)性研究,獲取醬香大曲中功能微生物,為提高醬香大曲發(fā)酵性能奠定理論與物質(zhì)基礎(chǔ)。