于學(xué)健，馮慧軍，翟磊，白秀彬，許玲，于盼盼，程池，姚粟*

1(中國(guó)食品發(fā)酵工業(yè)研究院有限公司，北京，100015) 2(天津科技大學(xué)，天津，300222)3(山東扳倒井股份有限公司，山東 高青，256300)

白酒獨(dú)特的生產(chǎn)工藝和風(fēng)味特征使其成為世界蒸餾酒的典型代表，是中國(guó)最具民族文化特色的產(chǎn)品之一，在我國(guó)擁有龐大的消費(fèi)群體和文化基礎(chǔ)[1]。芝麻香型白酒作為我國(guó)傳統(tǒng)白酒的創(chuàng)新香型，以獨(dú)特的清濃醬融合的風(fēng)味特點(diǎn)和優(yōu)雅的芝麻香而聞名，符合白酒淡雅、爽凈的消費(fèi)趨勢(shì)[2-4]。高溫大曲作為芝麻香型白酒釀造的微生物載體和糖化發(fā)酵劑，為白酒釀造提供了豐富的微生物菌種、生物酶類(lèi)以及獨(dú)特的芝麻香味成分和前體物質(zhì)，是芝麻香型白酒特殊風(fēng)味形成的關(guān)鍵因素之一[5-7]。微生物作為大曲制備的關(guān)鍵性因素，以功能基因?yàn)榛A(chǔ)產(chǎn)生多種生物酶類(lèi)和香味成分及其前體物質(zhì)，對(duì)白酒獨(dú)特風(fēng)味的形成具有重要作用[8-11]。

高溫大曲中的功能菌與白酒產(chǎn)香關(guān)系密切，采用16S rRNA基因克隆文庫(kù)技術(shù)分析芝麻香型白酒高溫大曲，結(jié)果表明，在大曲入房時(shí)以泛菌屬(Pantoea)為優(yōu)勢(shì)類(lèi)群，隨著制曲時(shí)間的延長(zhǎng)，溫度升高，高溫放線(xiàn)菌(Thermoactinomyces)等耐高溫菌屬成為優(yōu)勢(shì)菌系，其變化趨勢(shì)與品溫較一致，呈現(xiàn)先升高再降低的趨勢(shì)[12]；16S rRNA基因高通量測(cè)序結(jié)果也表明在高溫大曲中高溫放線(xiàn)菌為優(yōu)勢(shì)類(lèi)群[13]。通過(guò)ARDRA技術(shù)和分離培養(yǎng)方法可從高溫大曲中獲得普通高溫放線(xiàn)菌(T.vulgaris)、中間型高溫放線(xiàn)菌(T.intermedius)、大曲高溫放線(xiàn)菌(T.daqus)等多株高溫放線(xiàn)菌屬菌種[14]。該屬微生物常見(jiàn)于自然界高溫場(chǎng)所，具備耐高溫特性，適應(yīng)大曲發(fā)酵中的高溫環(huán)境，具備開(kāi)發(fā)高溫淀粉酶、蛋白酶、酯酶、脫氫酶和多酚氧化酶等潛力[15-18]，可促進(jìn)蛋白降解產(chǎn)生多種氨基酸，為酯化反應(yīng)及美拉德反應(yīng)提供前體物質(zhì)，促進(jìn)芝麻香風(fēng)味的形成，在芝麻香白酒釀造中起到重要作用[19-23]。這些微生物作用均與功能基因緊密相關(guān)，高溫大曲為白酒釀造提供豐富的微生物資源和功能基因，明析高溫放線(xiàn)菌中功能基因在大曲制備過(guò)程中的分布及變化規(guī)律，有助于深入解析該菌屬在大曲及白酒釀造過(guò)程的作用機(jī)制。

但目前對(duì)芝麻香型白酒高溫大曲的研究主要集中在菌種的分離篩選及微生物群落結(jié)構(gòu)分析上，對(duì)于芝麻香型大曲制備過(guò)程中特定菌種的功能基因組成及動(dòng)態(tài)變化規(guī)律研究還處于起步階段。為此，本研究針對(duì)芝麻香型白酒高溫大曲的制備過(guò)程，應(yīng)用宏基因組學(xué)分析技術(shù)，揭示白酒釀造高溫大曲中高溫放線(xiàn)菌及其功能基因的組成及動(dòng)態(tài)變化規(guī)律，為高溫大曲中特定功能酶系及菌株的篩選提供理論基礎(chǔ)，為解析高溫放線(xiàn)菌在白酒釀造中的具體功能作用及定向調(diào)控提供技術(shù)支撐，對(duì)探究芝麻香白酒風(fēng)味物質(zhì)形成機(jī)理、提升芝麻香型白酒生產(chǎn)品質(zhì)具有重要意義。

1 材料與方法

1.1 實(shí)驗(yàn)材料

1.1.1 試劑

E.Z.N.A.TMSoil DNA Kit DNA提取試劑盒，Omega Bio-tek, Inc., USA；MGIEasy DNA文庫(kù)制備試劑盒，深圳華大智造科技有限公司；Qubit?dsDNA HS Assay Kit試劑盒，賽默飛世爾科技(中國(guó))有限公司；DNA Marker，北京全式金生物技術(shù)有限公司；Gold View，北京塞百盛基因技術(shù)有限公司；實(shí)驗(yàn)中所使用的各種其他無(wú)特別說(shuō)明的化學(xué)試劑均為進(jìn)口或國(guó)產(chǎn)分析純產(chǎn)品。

1.1.2 儀器與設(shè)備

BGISEQ-500基因測(cè)序儀，深圳華大智造科技有限公司； Centrifuge 5424 R臺(tái)式高速冷凍離心機(jī)，德國(guó)Eppendorf公司；Eppendorf微量移液器，德國(guó)Eppendorf公司；HH-3A電熱恒溫水浴鍋，常州國(guó)華電器；DG-Ⅲ型雙穩(wěn)數(shù)顯電泳儀，北京鼎國(guó)生物技術(shù)發(fā)展公司；DYY-Ⅲ型水平電泳槽，北京六一儀器廠(chǎng)；紫外/可見(jiàn)凝膠成像系統(tǒng)，Bio-Rad公司。

1.2 實(shí)驗(yàn)方法

1.2.1 樣品采集

2017年7月開(kāi)展山東扳倒井股份有限公司18#曲房大曲的跟蹤采樣工作，根據(jù)不同發(fā)酵時(shí)期的發(fā)酵特點(diǎn)結(jié)合大曲發(fā)酵過(guò)程的操作工藝，對(duì)大曲制備過(guò)程中各個(gè)工藝階段樣品進(jìn)行采集，工藝階段包括入房(A，32 ℃)、入房第2天(B，34.5 ℃)、一翻前(C，入房4 d，41.3 ℃)、一翻(D，入房5 d，42.7 ℃)、一翻后(E，入房6 d，40.7 ℃)、二翻前(F，入房7 d，41.0 ℃)、二翻(G，入房8 d，41.7 ℃)、二翻后(H，入房9 d，45.5 ℃)、頂溫(I，入房29 d，60.0 ℃)、出房(J，入房36 d，36.3 ℃)；采用五點(diǎn)法在曲房的四周和中間各取一塊曲樣，共5塊大曲，粉碎曲樣，采用四分法進(jìn)行隨機(jī)取樣，樣品均分為兩部分，一部分用于微生物菌群分析，另一部分用于理化指標(biāo)分析[17-18]。對(duì)取樣時(shí)發(fā)酵時(shí)間、曲樣品溫、曲房環(huán)境溫度及相對(duì)濕度變化進(jìn)行記錄。采用山東扳倒井股份有限公司企業(yè)標(biāo)準(zhǔn)《大曲評(píng)分標(biāo)準(zhǔn)》(BDJ/JW—WSW-06)及《釀酒大曲通用分析方法》(QB/T 4257—2011)對(duì)大曲樣品水分、酸度、糖化力、發(fā)酵力、液化力和酯化力等理化性質(zhì)進(jìn)行測(cè)定。

1.2.2 大曲總DNA提取

采用E.Z.N.A.TMSoil DNA Kit提取第2批取樣的18#曲房大曲總DNA。基因組提取操作流程按照試劑盒說(shuō)明書(shū)進(jìn)行。

1.2.3 宏基因組測(cè)序

1.2.3.1 DNA文庫(kù)構(gòu)建及測(cè)序

采用Qubit?dsDNA HS Assay Kit試劑盒檢測(cè)大曲總DNA樣品的濃度，進(jìn)而采用MGIEasy DNA文庫(kù)制備試劑盒構(gòu)建測(cè)序文庫(kù)，操作流程嚴(yán)格按照使用說(shuō)明書(shū)進(jìn)行。

1.2.3.2 DNA文庫(kù)測(cè)序

采用BGISEQ-500基因測(cè)序儀對(duì)構(gòu)建的DNA文庫(kù)進(jìn)行測(cè)序，操作流程按照使用說(shuō)明書(shū)進(jìn)行。

1.2.4 宏基因組數(shù)據(jù)分析

1.2.4.1 數(shù)據(jù)質(zhì)控過(guò)濾

利用metaWRAP軟件(v1.0.2) 結(jié)合SOAPnuke(v1.5.6)軟件完成下機(jī)數(shù)據(jù)的質(zhì)控過(guò)濾，由于大曲中主要原料為小麥，在提取DNA過(guò)程中會(huì)引入小麥DNA，從而造成宿主污染，通過(guò)比對(duì)小麥基因組序列去除大曲DNA樣品中宿主序列，免除小麥宿主序列對(duì)后續(xù)分析的干擾。

1.2.4.2 宏基因組De Novo拼接

metaWRAP軟件(v1.0.2)整合了metaSPAdes和megaHIT兩個(gè)用于宏基因組組裝的軟件，使用metaSPAdes完成宏基因組拼接。

1.2.4.3 基因預(yù)測(cè)及基因集構(gòu)建

首先使用metaGeneMark(v2.10)在組裝結(jié)果中預(yù)測(cè)開(kāi)放閱讀框(ORFs)，然后將不用樣品預(yù)測(cè)出的基因合并在一起，使用CD-Hit進(jìn)行聚類(lèi)(序列相似度在95%以上)，以去除冗余序列。以大曲宏基因組注釋基因?yàn)榛A(chǔ)，將高溫放線(xiàn)菌屬中典型種的基因組序列在大曲宏基因組注釋基因數(shù)據(jù)中進(jìn)行比對(duì)，得到高溫放線(xiàn)菌屬功能基因分布情況。

1.2.4.4 物種組成分析

利用metaWRAP軟件(v1.0.2)結(jié)合Kraken(V2.0.7-beta)軟件完成物種組成分析。

1.2.4.5 基因注釋分析

利用GhostKOALA將數(shù)據(jù)在線(xiàn)提交至KEGG數(shù)據(jù)庫(kù)完成KEGG pathway分析；預(yù)測(cè)的基因使用blast與CAZy數(shù)據(jù)庫(kù)比對(duì)得到碳水化合物活性酶基因信息。

1.2.5 數(shù)據(jù)統(tǒng)計(jì)處理

相關(guān)數(shù)據(jù)的統(tǒng)計(jì)及繪圖利用Excel和R(v3.5.0)等完成。利用R語(yǔ)言中的Hmisc包開(kāi)展數(shù)據(jù)的相關(guān)性分析，利用R語(yǔ)言中的performance analytics包開(kāi)展相關(guān)性分析結(jié)果的可視化。

2 結(jié)果與分析

2.1 大曲中高溫放線(xiàn)菌屬的分布及變化規(guī)律

采用BGISEQ-500基因測(cè)序儀對(duì)大曲宏基因組樣品進(jìn)行測(cè)序，原始數(shù)據(jù)經(jīng)質(zhì)控過(guò)濾、去宿主序列、拼接組裝后，所得contig數(shù)量維持在5 000以上，平均contig長(zhǎng)度大于2 000 bp，滿(mǎn)足后續(xù)物種及基因注釋要求。

高溫放線(xiàn)菌屬(Thermoactinomyces)菌種在各樣品中隨入房時(shí)間的變化趨勢(shì)與品溫變化有較高的一致性，入房時(shí)基本無(wú)高溫放線(xiàn)菌屬菌種分布，隨著入房時(shí)間延長(zhǎng)，大曲品溫升高，高溫放線(xiàn)菌屬也隨之富集，在大曲頂溫過(guò)后溫度下降時(shí)，其相對(duì)豐度也呈下降趨勢(shì)(圖1)。這可能與高溫環(huán)境對(duì)不耐熱微生物的淘汰作用有關(guān)，也與高溫放線(xiàn)菌屬菌種生長(zhǎng)特性一致，在品溫40 ℃以上時(shí)，高溫放線(xiàn)菌屬菌種生長(zhǎng)良好，有利于其在大曲中的富集。

圖1 高溫放線(xiàn)菌屬在各樣品中的分布及變化情況Fig.1 Distribution and variation of Thermoactinomyces in Daqu samples

以高溫放線(xiàn)菌屬的分布情況與大曲理化性相關(guān)性分析表明，水分、酸度、液化力、品溫、入房時(shí)間均與大曲高溫放線(xiàn)菌屬顯著相關(guān)，其與大曲理化性質(zhì)等環(huán)境因子的相關(guān)性處于較強(qiáng)的水平，與水分、酸度呈顯著負(fù)相關(guān)關(guān)系，與液化力、品溫及入房時(shí)間呈顯著正相關(guān)關(guān)系(表1)。

表1 大曲核心微生物與理化性質(zhì)相關(guān)性分析Table 1 Correlation analysis between core microbiome and physicochemical properties of Daqu

注：只標(biāo)注顯著相關(guān)的因素，非顯著相關(guān)的因素為空白。

由表1可見(jiàn)，大曲中高溫放線(xiàn)菌屬與大曲理化性質(zhì)具備相互影響的可能性，也表明大曲的理化因素如發(fā)酵力、液化力、品溫等對(duì)于大曲中微生物群落具備潛在的調(diào)節(jié)作用，詳細(xì)了解兩者的相關(guān)關(guān)系，有助于解析微生物對(duì)大曲品質(zhì)的作用，可為定向調(diào)控大曲中微生物組成進(jìn)而提升大曲品質(zhì)提供理論基礎(chǔ)。

2.2 大曲中高溫放線(xiàn)菌功能基因的分布及變化規(guī)律

以大曲宏基因組注釋基因?yàn)榛A(chǔ)，將高溫放線(xiàn)菌屬中典型種的基因組序列在大曲宏基因組注釋基因數(shù)據(jù)中進(jìn)行比對(duì)，得到高溫放線(xiàn)菌屬的功能基因分布，結(jié)果表明在大曲宏基因數(shù)據(jù)中共比對(duì)到1 449個(gè)KEGG注釋基因(圖2)。KEGG注釋的基因主要分為genetic information processing(210)，environmental information processing(160)，metabolism(666)(global and overview maps、carbohydrate metabolism、amino acid metabolism、nucleotide metabolism、engry metabolism、lipid metabolism、metabolism of cofactors and vitamins)，cellular process(78)、organismal system(41)、human diseases(81)6個(gè)主要功能組別。其中與metabolism相關(guān)的基因數(shù)量約占被注釋全部基因的45.96%。metabolism中具有較多基因數(shù)目的代謝途徑有carbohydrate metabolism、amino acid metabolism、nucleotide metabolism等途徑，與細(xì)胞基本代謝相關(guān)。

圖2 高溫放線(xiàn)菌在大曲宏基因組中比對(duì)到的KEGG注釋基因數(shù)量Fig.2 Annotation of Daqu metagenomic genes in KEGG database

對(duì)高溫放線(xiàn)菌基因在各工藝階段大曲樣品中的分布和變化情況統(tǒng)計(jì)如圖3，入房至二翻的工藝階段，其占宏基因組的相對(duì)豐度總量相似，在二翻后、頂溫及出房大曲中功能基因相對(duì)豐度升高，這與該菌屬的相對(duì)豐度分布較為一致。入房大曲中，轉(zhuǎn)運(yùn)ATP酶(K02112)、精氨基琥珀酸裂解酶(K01755)、丙酮酸激酶(K00873)為優(yōu)勢(shì)基因，這些基因主要參與微生物的基礎(chǔ)代謝，涉及微生物的氧化磷酸化反應(yīng)、氨基酸代謝反應(yīng)以及糖酵解反應(yīng)。

圖3 各工藝階段大曲中高溫放線(xiàn)菌KEGG注釋基因分布Fig.3 Distribution of Thermoactinomyces genes annotated in KEGG database in Daqu with different fermentation stages (10 genes with the highest relative abundance in each sample were selected)

隨入房時(shí)間延長(zhǎng)大曲品溫升高，堿性磷酸酶蛋白(Clp, K03695)、groEL分子伴侶(K04077)、磷酸轉(zhuǎn)運(yùn)系統(tǒng)ATP結(jié)合蛋白(pstB, K02036)等成為優(yōu)勢(shì)基因(圖4)，Clp屬于熱休克蛋白HSP100家族，在高溫脅迫下可與其他結(jié)構(gòu)不穩(wěn)定的蛋白結(jié)合促進(jìn)蛋白質(zhì)的穩(wěn)定性，防止蛋白質(zhì)凝聚變性，并可協(xié)同其他分子伴侶修復(fù)變形的蛋白，在抵抗高溫脅迫過(guò)程中起到重要作用[24-25]。groEL分子伴侶屬于熱休克蛋白HSP60家族，可提高蛋白質(zhì)的正確折疊與組裝效率，在熱脅迫的變性蛋白恢復(fù)過(guò)程中發(fā)揮重要作用[26-27]。pstB屬于ABC轉(zhuǎn)運(yùn)子家族，可以結(jié)合并水解ATP，可為細(xì)胞的磷酸轉(zhuǎn)運(yùn)系統(tǒng)提供能量，參與磷的轉(zhuǎn)運(yùn)功能，與細(xì)胞生物活性密切相關(guān)[28]。這些與微生物耐熱相關(guān)的基因菌存在隨溫度升高而上升的趨勢(shì)，說(shuō)明隨大曲品溫升高，高溫放線(xiàn)菌啟動(dòng)了相應(yīng)的耐熱應(yīng)答機(jī)制，相關(guān)耐熱基因含量提升，為其在高溫環(huán)境中的生存提供良好的物質(zhì)基礎(chǔ)。

圖4 各工藝階段大曲中高溫放線(xiàn)菌耐熱基因分布Fig.4 Distribution of Clp, groEL and pstB genes of Thermoactinomyces in Daqu samples

相關(guān)性分析顯示，Clp,groEL,pstB基因和高溫放線(xiàn)菌豐度及大曲品溫之間，均呈顯著正相關(guān)關(guān)系(圖5)。

圖5 Clp, groEL, pstB基因、高溫放線(xiàn)菌豐度和大曲品溫的相關(guān)性分析Fig.5 Correlation analysis of Clp, groEL, pstB genes,Thermoactinomyces and Daqu temperature

這也說(shuō)明了這些潛在耐熱基因的豐度與高溫放線(xiàn)菌保持較高的一致性，都在大曲品溫提升時(shí)，高溫放線(xiàn)菌相對(duì)豐度升高，相應(yīng)的耐熱基因數(shù)量也呈現(xiàn)增高趨勢(shì)，從而提升高溫放線(xiàn)菌在高溫下的生存能力進(jìn)而發(fā)揮其功能作用。

大曲作為白酒釀造的糖化發(fā)酵劑，含豐富的碳水化合物活性酶，解析其中相關(guān)酶的基因分布及變化規(guī)律，有助于明析其對(duì)原料碳水化合物的代謝能力及作用機(jī)制，為進(jìn)一步研究大曲中微生物的作用機(jī)制及調(diào)控技術(shù)奠定理論基礎(chǔ)。以大曲宏基因組的碳水化合物活性酶數(shù)據(jù)庫(kù)(CAZy)注釋基因?yàn)榛A(chǔ)，高溫放線(xiàn)菌在其中共比對(duì)到17個(gè)基因，主要分為糖基水解酶(GH，8個(gè))、糖酯酶(GE，4個(gè))、糖基轉(zhuǎn)移酶(GT，4個(gè))、碳水化合物結(jié)合結(jié)構(gòu)域(CBM，1個(gè))(圖6)。

圖6 高溫放線(xiàn)菌在大曲宏基因組中比對(duì)到的CAZy注釋基因數(shù)量Fig.6 Annotation of Daqu metagenomic genes in CAZy database

高溫放線(xiàn)菌CAZy注釋基因在入房至二翻前的工藝階段占宏基因組的相對(duì)豐度總量相似，在二翻之后大曲中功能基因相對(duì)豐度逐漸升高，出房曲中碳水化合物活性酶基因的種類(lèi)和豐度達(dá)到峰值(圖7)，預(yù)示高溫放線(xiàn)菌相關(guān)碳水化合物代謝能力在出房曲中大幅提升，也有助于大曲發(fā)揮相應(yīng)功能作用。入房大曲中，高溫放線(xiàn)菌只有α-淀粉酶基因分布，隨制曲時(shí)間延長(zhǎng)，大曲高溫放線(xiàn)菌中碳水化合物活性酶基因數(shù)量與種類(lèi)增多，出房大曲中相關(guān)基因達(dá)17種，主要由α-淀粉酶、乙?；揪厶酋ッ负推咸烟擒彰傅冉M成，這些基因的富集可促進(jìn)大曲原料中的碳水化合物的降解，為微生物生長(zhǎng)提供必要的能量及代謝底物。

圖7 各工藝階段大曲中高溫放線(xiàn)菌CAZy注釋基因分布Fig.7 Distribution of Thermoactinomyces genes annotated in CAZy database in Daqu with different fermentation stages

相關(guān)性分析顯示，高溫放線(xiàn)菌中糖基水解酶(GH)、糖基轉(zhuǎn)移酶(GT)基因與大曲糖化力顯著正相關(guān)，糖基水解酶(GH)、糖基轉(zhuǎn)移酶(GT)、糖酯酶(GE)基因與大曲酯化力顯著正相關(guān)(圖8)，可見(jiàn)高溫放線(xiàn)菌中碳水化合物活性酶相關(guān)基因?qū)Υ笄腔εc酯化力的提升有潛在的促進(jìn)作用，也預(yù)示了高溫放線(xiàn)菌具備潛在的碳水化合物降解能力，促進(jìn)大曲中原料的利用，有助于提升大曲發(fā)酵品質(zhì)。入房時(shí)間與糖基水解酶(GH)基因、糖基轉(zhuǎn)移酶(GT)基因、CZAy注釋基因、高溫放線(xiàn)菌豐度、糖化力及酯化力均顯著正相關(guān)，與水分、酸度顯著負(fù)相關(guān)，說(shuō)明隨大曲發(fā)酵時(shí)間延長(zhǎng)，大曲本身水分降低，酸度下降，其中高溫放線(xiàn)菌豐度上升，碳水化合物降解相關(guān)基因逐漸富集，糖化力與酯化力提升，促使大曲發(fā)酵成熟。

圖8 高溫放線(xiàn)菌CAZy注釋基因與理化性質(zhì)相關(guān)性分析Fig.8 Correlation analysis between Thermoactinomyces genes annotated in CAZy database and physicochemical properties of Daqu

3 結(jié)論

本研究以芝麻香型白酒高溫大曲為研究對(duì)象，通過(guò)宏基因組測(cè)序技術(shù)解析了其制備過(guò)程中高溫放線(xiàn)菌屬及其功能基因的動(dòng)態(tài)變化規(guī)律，結(jié)果表明高溫放線(xiàn)菌屬在大曲入房時(shí)基本無(wú)分布，隨著大曲品溫升高，高溫放線(xiàn)菌屬也隨之富集，整體呈現(xiàn)先升高再降低的變化過(guò)程，其變化趨勢(shì)與品溫呈顯著正相關(guān)關(guān)系；高溫放線(xiàn)菌功能基因在入房大曲中以轉(zhuǎn)運(yùn)ATP酶、精氨基琥珀酸裂解酶、丙酮酸激酶為優(yōu)勢(shì)基因，隨入房時(shí)間延長(zhǎng)大曲品溫升高，Clp堿性磷酸酶蛋白、groEL分子伴侶、磷酸轉(zhuǎn)運(yùn)系統(tǒng)ATP結(jié)合蛋白等潛在耐熱基因成為優(yōu)勢(shì)基因，從而有效提升高溫放線(xiàn)菌在高溫下的生存能力，為其在大曲中發(fā)揮相關(guān)功能奠定基礎(chǔ)。入房時(shí)間與CZAy注釋基因、高溫放線(xiàn)菌豐度、糖化力及酯化力均顯著正相關(guān)，與水分、酸度顯著負(fù)相關(guān)，糖基水解酶、糖基轉(zhuǎn)移酶基因與大曲糖化力與酯化力顯著正相關(guān)，說(shuō)明大曲發(fā)酵過(guò)程中水分與酸度降低，其中高溫放線(xiàn)菌豐度上升，碳水化合物降解相關(guān)基因逐漸富集，促進(jìn)大曲糖化力與酯化力提升，初步揭示了高溫放線(xiàn)菌及其功能基因與大曲理化性質(zhì)的關(guān)系，為微生物及高溫大曲在白酒釀造中的功能作用解析及定向調(diào)控提供技術(shù)支撐，對(duì)提升芝麻香型白酒生產(chǎn)品質(zhì)具有重要意義。