母應春，姜 麗，蘇 偉

（貴州大學釀酒與食品工程學院，貴州 貴陽 550025）

酒曲是在蒸煮的谷物及其輔料上接種或天然定型發(fā)酵而成。我國用曲釀酒已有幾千年的歷史，歷代勞動人民傳承下來的民間酒曲是我國寶貴的微生物資源，隨著科學技術與釀酒業(yè)的發(fā)展，對酒曲的研究也日益增多[1-2]，酒曲是中國酒釀造的糖化、發(fā)酵和生香劑，含多種微生物及其產(chǎn)生的酶，酒曲品質對酒的產(chǎn)率和品質影響極大，故有“曲是酒之骨”之稱[3]。

貴州省位置獨特，少數(shù)民族眾多，酒文化深厚。很多地區(qū)至今保持著傳統(tǒng)的釀酒工藝以及傳統(tǒng)的酒曲制作技術。傳統(tǒng)曲一般是自然接種微生物，因此微生物種類繁多。貴州省花溪、盤州、安順酒曲是采用傳統(tǒng)工藝制成，糖化力、液化力都相對較強，且出酒率高，品質香醇。近年來，國內外有許多的酒曲微生物多樣性的研究。Hui Wenyan等[4]使用單分子實時測序技術發(fā)現(xiàn)3 個日本酒曲中存在大量隸屬于Ochrobactrum的細菌。沈馨等[5]采用高通量測序技術對3 個孝感鳳窩酒曲的細菌多樣性分析發(fā)現(xiàn)，鳳窩酒曲中存在大量的核心細菌菌群。相飛[6]采用變性梯度凝膠電泳技術發(fā)現(xiàn)Weissella kimchi、Enterococcus faecium、Herbaspirilum sp.為辣蓼甜酒曲中的優(yōu)勢菌。劉孟華等[7]利用宏基因組學的方法對劍南春中細菌群落多樣性的研究得出，劍南春酒曲微生物主要分布在Proteobacteria、Firmicutes、Actinobacteria三大門中，占總量的95%以上。然而有關貴州酒曲的微生物多樣性的研究報道尚少。

高通量測序技術具有通量高、成本低、準確率高、可進行雙向測序等優(yōu)點[8-10]，可以方便快捷地分析樣品中復雜的微生物菌群結構，是分析復雜多樣微生物樣品的首選[11-12]。本研究利用高通量測序技術及相關性分析方法對貴州3 個地區(qū)生產(chǎn)的酒曲中微生物組成差異及多樣性進行比較，以期為不同酒曲中微生物資源開發(fā)利用提供參考。

1 材料與方法

1.1 材料與試劑

花溪酒曲，編號酒曲A；盤州酒曲，編號酒曲B；安順酒曲，編號酒曲C。

DNA試劑提取盒 德國Qiagen公司；2%瓊脂糖凝膠 西班牙Biowest公司；E.Z.N.A.?Soil試劑盒上海玉博生物科技有限公司；AxyPrep DNA凝膠提取試劑盒 AxyPrep DNA Gel Extraction Kit 美國Axygen Biosciences公司；QuantiFluor?-ST微型熒光劑 美國Promega公司；相關生理生化實驗所用試劑均為進口或國產(chǎn)分析純。

1.2 儀器與設備

Fresco21型高速冷凍離心機 美國Thermo公司；GeneAmp?9700型聚合酶鏈式反應（polymerase chain reaction，PCR）儀 美國ABI公司；Gel DOC XR凝膠成像系統(tǒng) 美國Bio-Rad公司。

1.3 方法

1.3.1 樣品采集及前處理

采集來自貴州省3 個地區(qū)（花溪、盤州、安順）不同酒廠生產(chǎn)的酒曲，采集時間2017年9月，所有樣品經(jīng)無菌密封好后放入-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

1.3.2 酒曲的理化指標檢測

水分測定：直接干燥法；酸度測定：酸堿中和法（以乳酸計）；糖化酶活力測定：婓林快速法；液化酶活力測定：分光光度法。以上方法參考文獻[13]。

1.3.3 DNA抽提和PCR擴增

根據(jù)E.Z.N.A.?Soil試劑盒說明書進行總DNA抽提，DNA濃度和純度利用NanoDrop2000進行檢測，利用1%瓊脂糖凝膠電泳檢測DNA提取質量；細菌用338F（5’-ACTCCTACGGGAGGCAGCAG-3’）和806R（5’-GGACTACHVGGGTWTCTAAT-3’）引物對V3-V4可變區(qū)進行PCR擴增；真菌用ITS1F（5’-CTTGGTCATTTAGAGGAAGTAA-3’）和ITS2R（5’-GCTGCGTTCTTCATCGATGC-3’）引物對ITS1可變區(qū)進行PCR擴增；擴增程序為：95 ℃預變性3 min，27 個循環(huán)（95 ℃ 變性30 s，55 ℃退火30 s，72 ℃延伸30 s），最后72 ℃延伸10 min。擴增體系為20 μL，4 μL 5×FastPfu緩沖液，2 μL 2.5 mmol/L dNTPs，0.8 μL引物（5 μmol/L），0.4 μL FastPfu聚合酶；10 ng DNA模板。

1.3.4 Illumina MiSeq測序

使用2%瓊脂糖凝膠回收PCR產(chǎn)物，利用AxyPrep DNA凝膠提取試劑盒進行純化，Tris-HCl洗脫，2%瓊脂糖電泳檢測。利用QuantiFluor?-ST微型熒光劑進行檢測定量。根據(jù) Illumina MiSeq平臺標準操作規(guī)程將純化后的擴增片段構建PE 2×300的文庫。

構建文庫步驟：1）連接“Y”字形接頭；2）使用磁珠篩選去除接頭自連片段；3）利用PCR擴增進行文庫模板的富集；4）氫氧化鈉變性，產(chǎn)生單鏈DNA片段。利用Illumina公司的MiSeq PE300平臺進行測序（上海美吉生物醫(yī)藥科技有限公司）。原始數(shù)據(jù)上傳至NCBI數(shù)據(jù)庫中。

1.3.5 數(shù)據(jù)處理

原始測序序列使用Trimmomatic軟件質控，使用FLASH軟件進行拼接：1）設置50 bp的窗口，如果窗口內的平均質量值低于20，從窗口開始截去后端堿基，去除質控后長度低于50 bp的序列；2）barcode需精確匹配，引物允許2 個堿基的錯配，去除模糊堿基；3）根據(jù)重疊堿基overlap將兩端序列進行拼接，overlap需大于10 bp。去除無法拼接的序列。使用UPARSE軟件（version 7.1 http://drive5.com/uparse/），根據(jù)97%的相似度對序列進行操作分類單元（operational taxonomic unit，OTU）聚類；使用UCHIME 軟件剔除嵌合體。利用RDP classifier（http://rdp.cme.msu.edu/）對每條序列進行物種分類注釋，比對Silva數(shù)據(jù)庫（SSU123），設置比對閾值為70%。

1.4 數(shù)據(jù)分析

利用Excel 2010、Origin 9、SPSS Statistics 20.0等統(tǒng)計軟件進行數(shù)據(jù)分析。

2 結果與分析

2.1 酒曲中基本理化性質分析

表1 酒曲基本理化性質分析Table 1 Basic physicochemical properties of rice wine kojis

由表1可知，3 種酒曲的水分質量分數(shù)在10.04%～13.273%之間，酒曲A和酒曲C不顯著（P＞0.05），酒曲B和酒曲C之間差異不顯著（P＞0.05）；酸度在0.421%～1.501%之間，酒曲A、B、C之間差異顯著（P＜0.05）；液化酶活力在77.420～91.880 U/g之間，酒曲A、B、C之間差異顯著（P＜0.05）；糖化酶活力在346.282～673.540 U/g之間，其中酒曲A和酒曲B 之間差異不顯著（P＞0.05），酒曲B和酒曲C之間差異不顯著（P＞0.05）。3 個地區(qū)的酒曲在酸度和液化酶活力上差異顯著。

2.2 酒曲中DNA提取與定量

提取3 個地方酒曲的DNA，測定其質量與純度，并觀察OD260nm/OD280nm和OD260nm/OD230nm值（表2）可知，酒曲樣品的DNA的OD260nm/OD280nm值在1.8～2.3之間，OD260nm/OD230nm值在0.04～0.17之間，其質量濃度在8.5～42.8 ng/μL之間，說明提取的DNA濃度較純，以便用于酒曲的細菌和真菌多樣性分析。

表2 酒曲中DNA質量濃度分析Table 2 Microbial DNA concentrations from rice wine koji

2.3 樣品中細菌的分析

2.3.1 細菌豐度和多樣性分析

表3 酒曲樣品細菌測序數(shù)據(jù)統(tǒng)計分析Table 3 Statistical analysis of bacterial sequencing data of koji samples

除去不合格序列后，3 種酒曲樣品中共獲得129 510 條高質量細菌序列（表3），其中從酒曲A、B、C中分別獲得45 308、40 940、43 262 條有效序列，基于97%相似度其OTU數(shù)量分別為474、89和476。ACE指數(shù)在246.09～476.27之間，Shannon指數(shù)在1.64～4.42之間，Chao 1指數(shù)在139.65～477.67之間，覆蓋率范圍99.87%～99.97%。說明所測得的OTU基本能反映所研究的酒曲樣品細菌群落組成。每個樣品的Chao 1指數(shù)和Shannon指數(shù)用于評估物種豐度和多樣性[14]。通過計算樣品的多樣性指數(shù)發(fā)現(xiàn)，酒曲A和酒曲C的Chao 1指數(shù)和Shannon指數(shù)明顯高于酒曲B，這說明酒曲A和酒曲C中的物種豐度明顯高于酒曲B，且酒曲A和酒曲C的細菌群落結構比較復雜。由圖1可知，隨著測序深度的增加物種的數(shù)量也隨之增加，但與此同時每個樣品的Shannon曲線已進入平臺期[15]。由此可知，雖然隨著測序量的增加有可能會發(fā)現(xiàn)新的細菌種系型，但樣品細菌多樣性已經(jīng)不再隨測序量的增加而發(fā)生變化，說明測序數(shù)量足夠充分合理[16]。因此，本研究細菌的測序量滿足后續(xù)的生物信息學分析要求。

圖1 3 種酒曲中細菌的稀釋性曲線（a）和Shannon指數(shù)（b）圖Fig. 1 Rarefaction curve (a) and Shannon index (b) of bacteria in three rice wine kojis

2.3.2 細菌中微生物群落組成分析

在門水平上，3 個樣品中共檢測出8 個細菌門，分別為厚壁菌門（Firmicutes）、變形桿菌門（Proteobacteria）、放線菌門（Actinobacteria）、擬桿菌（Bacteroides）、綠彎菌門（Chloroflexi）、芽單胞菌門（Gemmatimonadetes）、異常球菌-棲熱菌門（Deinococcus-Thermus）、未分類菌群（others）。如圖2a所示，相對豐度高且3 個樣品均有分布的厚壁菌門、變形菌門、放線菌門。不同地區(qū)的酒曲樣品之間微生物菌群明顯不同，其中酒曲A和酒曲C中優(yōu)勢的細菌門為變形菌門，比例高達31.86%，39.98%；放線菌門為酒曲A中第2大優(yōu)勢菌群，占樣品序列的29.31%；酒曲B中厚壁菌門為優(yōu)勢細菌門，占樣品序列的98.08%，而酒曲C中厚壁菌門的相對豐度僅次于變形菌門，占樣品序列的29.28%，居第2位。不同的酒曲樣品之間微生物菌群不同，3 個樣品中酒曲A和酒曲C的細菌微生物群落組成相似，但相對豐度存在差異。

在屬水平上，3 個樣本中共檢測到19 個細菌屬，分別為芽孢桿菌屬（Bacillus）、片球菌屬（Pediococcus）、紅球菌屬（Rhodococcus）、羅爾斯頓菌屬（Ralstonia）、不動細菌屬（Acinetobacter）、腸球菌屬（Enterococcus）、乳酸桿菌屬（Lactobacillus）、Burkholderia-Paraburkholderia、雙歧桿菌屬（Bifidobacterium）、norank-f-Mitochondria、鞘氨醇單胞菌屬（sphingomonas）、Hydrogenispora、norankc-Gemmatimonadetes、Faecalibaculum、unclassified-f-Micromonospraceae、棲熱菌屬（Thermus）、norank-f-OM1-clade、Alistipesothers、未分類菌群（others）。如圖2b所示，相對豐度較高且3 個樣品中均分布的主要細菌屬有芽孢桿菌屬、片球菌屬、紅球菌屬、羅爾斯頓菌屬和不動細菌屬。其中酒曲A和酒曲C中細菌微生物組成較為復雜，且相似。但所占比例不同，主要以紅球菌屬、青枯菌屬、芽孢桿菌屬、不動細菌屬等為主。酒曲B中微生物多樣性最低，主要有芽孢桿菌屬、不動細菌屬和紅球菌屬等為主。紅球菌屬為酒曲A的優(yōu)勢菌屬，占樣本序列比例為18.50%；其次是羅爾斯頓菌屬，占樣本序列比例為8.46%；芽孢桿菌屬為酒曲B的優(yōu)勢菌屬，占樣本序列比例為13.09%；其次是片球菌屬和腸球菌屬，分別占樣本序列比例為11.15%，11.34%。而酒曲C中的優(yōu)勢菌屬也是芽孢桿菌屬，占樣本序列比例61.49%，其次是不動細菌屬，占樣本序列比例的23.74%。

圖2 3 種酒曲在門水平（a）和屬水平（b）的細菌組成Fig. 2 Bacterial community composition of three rice wine kojis at the phylum level (a) and genus level (b)

2.3.3 細菌微生物群落相似性分析

通過非度量多維尺度（nonmetric multidimensional scaling，NMDS）分析細菌群落間存在的關系，根據(jù)樣品中包含的物種信息，以點的形式反映在多維空間上，而對不同樣品間的差異程度，則是最終通過點與點的距離表現(xiàn)在圖中，其距離遠近代表樣品中菌體群落的差異程度[17]。由圖3a可知，酒曲A和酒曲C距離較近，聚集在圖的左側，群落差異較?。欢魄鶥的細菌微生物聚集在右側部分，距酒曲A和酒曲C較遠，差異較大，構成一個獨立的組群。通過對樣品細菌Beta多樣性距離矩陣進行層級聚類分析，構建樣品層級聚類樹研究不同樣品的相似性和差異性[18]。如圖3b所示，Beta分析結果與NMDS相似，3 種酒曲共聚為兩大支，酒曲A和酒曲C的細菌屬聚為一類，差異較小，酒曲B的細菌屬另聚為一類。

圖3 細菌非度量多維尺度分析圖（a）和基于Beta多樣性距離的樣品層級聚類樹（b）Fig. 3 Bacterial non-metric multi-dimensional analysis (a) and hierarchical dendrogram based on Beta diversity distance (b)

2.4 樣品中真菌的分析

2.4.1 真菌的豐度和多樣性分析

酒曲A、B、C經(jīng)過測序后得出3 種酒曲樣品中共獲得60 945 條高質量真菌序列（表4），3 種酒曲中獲得有效序列相同，均為20 315 條，基于97%相似度，其OTU數(shù)量分別為2.4、1和15.8。ACE指數(shù)在0～16.09之間，Shannon指數(shù)在0.06～1.06之間，Chao 1指數(shù)在5.5～25之間，覆蓋率均為99.99%。說明所測得的OTU基本能反映所研究酒曲樣品真菌群落組成。通過計算樣品的多樣性指數(shù)發(fā)現(xiàn)，酒曲A的Chao 1指數(shù)25最大，說明酒曲A的真菌微生物豐度高于酒曲B和酒曲C。在多樣性方面，酒曲C的Shannon指數(shù)1.06最大，說明酒曲C的真菌群落結構較為復雜。由圖4可知，隨著測序深度的增加，發(fā)現(xiàn)物種的數(shù)量也隨之增加，與此同時，每個樣品的Shannon曲線己進入平臺期。由此可知，雖然隨著測序量的增加有可能會發(fā)現(xiàn)新的真菌種系型，但樣品真菌的多樣性己經(jīng)不再隨測序量的增加而發(fā)生變化。因此，本研究真菌的測序量滿足后續(xù)的生物信息學分析要求。

表4 酒曲樣品真菌測序數(shù)據(jù)統(tǒng)計分析Table 4 Statistical analysis of fungal sequencing data of koji samples

圖4 3 種酒曲中真菌的稀釋性曲線（a）和Shannon指數(shù)（b）圖Fig. 4 Dilution curve (a) and Shannon index (b) of fungi in three kinds of koji

2.4.2 真菌中微生物群落組成分析

圖5 3 種酒曲在門水平（a）和屬水平（b）上的真菌組成Fig. 5 Fungal community composition of three kojis at the phylum level (a) and genus level (b)

如圖5a所示，在門水平上，檢測出3 個真菌門為接合菌門（Zygomycota）、子囊菌門（Ascomycota）和擔子菌門（Basidiomycota）。在門水平上相對豐度高且3個樣品均有分布的主要真菌門為子囊菌門。不同的酒曲樣品之間微生物菌群明顯不同。酒曲A和酒曲B中接合菌門為優(yōu)勢真菌門，比例達64.57%、98.91%。酒曲C中子囊門為優(yōu)勢真菌門，占樣品序列比例95.26%。

如圖5b所示，在屬水平上，3 個樣本中共檢測到4 個真菌屬，分別為米根霉屬（Rhizopus）、酵母菌屬（Saccharomyce）、曲霉菌屬（Aspergillus）和嗜熱菌屬（Thermomyces）。酒曲A和酒曲B中米根霉屬為優(yōu)勢真菌屬，占樣本序列比例的66.58%和95.25%。酒曲C中優(yōu)勢真菌屬為酵母屬，占樣本序列比例的95.25%。3 種酒曲樣品中盡管微生物群落組成相似，但相對豐度仍存在顯著性差異。

2.4.3 真菌微生物群落相似性分析

圖6 真菌非度量多維尺度分析圖（a）和基于Beta多樣性距離的樣品層級聚類樹（b）Fig. 6 Fungal non-metric multidimensional scaling analysis (a) and hierarchical clustering dendrogram based on beta diversity distance (b)

如圖6a可知，真菌NMDS聚類與細菌結果不同，酒曲A和酒曲B距離較近，聚集在圖的右側，群落差異較小。而酒曲C的細菌微生物聚集在左側部分，距酒曲A和酒曲B較遠，差異較大，構成一個獨立的組群。如圖6b所示，真菌Beta多樣性距離矩陣層級聚類結果與NMDS結果一樣，3 種酒曲同樣共聚為兩大支，酒曲A和酒曲B聚為一類，酒曲C另聚為一類。

2.5 典型關聯(lián)分析（canonical correlation analysis，CCA）

CCA揭示不同地區(qū)酒曲樣品微生物群落組成與不同環(huán)境因子之間的相關性[19]。選擇3 個地方酒曲理化指標作為環(huán)境因子與酒曲中屬水平上細菌豐度前10的微生物進行典范對應分析。如圖7所示，酸度與液化酶活力之間夾角為鈍角，表明兩個理化指標之間不具有協(xié)同效應；水分質量分數(shù)與糖化酶活力之間夾角為銳角，表明兩個指標之間具有協(xié)同效應。酒曲A和酒曲C與酸度相關性最高，酒曲B與液化酶活力、水分質量分數(shù)、糖化酶活力相關性最高，Acinetobacter、Ralstonia、Rhodococcus、norank-f-Mitochondria、Bifidobacterium、Lactobacillus與酸度變化呈正相關；Bacillus、Enterococcus、Pediococcus與液化酶活力、糖化酶活力變化呈正相關，與水分質量分數(shù)變化也呈正相關但影響不顯著。由此說明，酒曲中微生物的種群結構對酒曲理化指標具有顯著影響。

圖7 細菌群落與環(huán)境因子典范對應分析Fig. 7 Canonical correspondence analysis between bacterial community and environmental factors

3 討 論

通過高通量測序分析3 種不同酒曲的群落結構及物種多樣性之間的差異，從3 種酒曲中共檢出8 個細菌門和3 個真菌門；19 個細菌屬和4 個真菌屬。細菌門主要有厚壁菌門、變形菌門、放線菌門、擬桿菌門、綠彎菌門、芽單胞菌門等。細菌屬主要有紅球菌屬、羅爾斯頓菌屬、不動細菌屬、乳酸桿菌屬、雙歧桿菌屬、芽孢桿菌屬、伯克爾德屬、片球菌屬、腸球菌屬等。真菌門有接合菌門、子囊菌門、擔子菌門。真菌屬有米根霉屬、酵母菌屬、曲霉菌屬、嗜熱菌屬。據(jù)報道，酒曲中的細菌屬有芽孢桿菌屬、乳酸桿菌屬、乳球菌屬、鏈球菌屬、芽孢梭菌屬、葡萄球菌屬、微細菌屬、節(jié)細菌屬、歐文氏菌屬、醋桿菌屬、柯匝克氏菌屬、顆粒鏈菌屬和希瓦氏菌屬13 個菌屬[20-22]，由此可知，酒曲樣品中微生物群落組成豐富，不同地區(qū)酒曲微生物群落組成雖然相似，但是豐度仍存在差異。3 個地區(qū)中，酒曲A和酒曲C的細菌種群相似，酒曲A和酒曲B的真菌種群相似，這應該與環(huán)境因素有所關系，并且揭示了酒曲中微生物組成種類的復雜程度。

張雙燕等[23]研究發(fā)現(xiàn)，大曲中優(yōu)勢細菌屬有乳酸桿菌屬，片球菌屬，也檢測出醋酸桿菌屬和芽孢桿菌屬等。其中乳桿菌屬，片球菌屬和芽孢桿菌屬在本實驗3 種酒曲中均被檢測出，芽孢桿菌屬為酒曲B和酒曲C的第一優(yōu)勢菌屬。芽孢桿菌屬在白酒釀造過程中發(fā)揮著重要的作用，能產(chǎn)生淀粉酶和蛋白酶，分解淀粉和蛋白質等大分子物質。此外，芽孢桿菌能影響白酒的風味物質如酯、醇的形成[24-25]。李艷等[26]研究發(fā)現(xiàn)，羊羔美酒大曲中優(yōu)勢菌屬有片球菌屬和乳酸桿菌屬，甜酒中酸味較重可歸功于含量較高的片球菌屬和乳酸桿菌屬，表明了甜酒最終質量和風味與發(fā)酵微生物分類組成息息相關。本研究也發(fā)現(xiàn)，在酒曲B中片球菌屬的豐度較高，但乳酸桿菌屬豐度略低。Bora等[27]采用基于Illumina平臺的全基因組鳥槍測序法對印度的米酒曲xaj-pitha進行研究發(fā)現(xiàn)，米酒曲含有大量以乳酸菌為其主要細菌類群的細菌屬。3 種酒曲檢測的結果均有乳酸桿菌屬，但相對豐度各有不同，乳酸桿菌屬與乳酸乙醋的形成具有密切關系[28]，其對于有機酸具有一定的耐受能力[29]，能減弱酵母菌的好氧代謝速度，延長前發(fā)酵期，有利于發(fā)酵有益菌的生長，對于白酒香味物質增加有利[30-31]。另外3 種酒曲中也檢測到雙歧桿菌屬，雙歧桿菌是有益菌的代表，對于維持腸道微生態(tài)平衡有重要作用，具有調整腸道功能、凈化腸道環(huán)境、促進機體健康的效果，提高抗病能力等[32]。陳玲等[33]采用16S rDNA克隆文庫法和高通量測序法分析大曲中細菌微生物組成，檢測出在門水平上的細菌門有厚壁菌門、變形菌門、放線菌門、擬桿菌門等。王丹丹等[34]采用高通量測序技術對孝感鳳窩酒曲的真菌微生物多樣性進行評價，得出在門水平上真菌主要有子囊菌門。本研究也在3 種酒曲樣品中檢測到了厚壁菌門、變形菌門、放線菌門、擬桿菌門、子囊菌門的存在，其在酒曲A和酒曲C中子囊菌門的豐度最高，而在酒曲C中相對豐度較低。接合菌門為酒曲C中的優(yōu)勢菌門。喬曉梅等[35]研究發(fā)現(xiàn)，清香大曲真菌菌群有米根霉和酵母菌。羅惠波等[36]采用ITS基因文庫法研究清香型大曲真菌群落結構，得出清茬曲的優(yōu)勢真菌屬為畢赤酵母屬和米根霉屬。而本研究也檢測到酒曲A和酒曲B中的優(yōu)勢真菌屬為米根霉屬，酒曲C中的優(yōu)勢真菌屬為酵母屬。

本實驗對3 種酒曲微生物群落結構與組成及多樣性進行研究，采用CCA酒曲的基本理化指標與酒曲細菌種群結構關系，發(fā)現(xiàn)不同地區(qū)的酒曲微生物群落構成及多樣性存在明顯差異，且酒曲的基本理化指標與微生物菌群結構有一定的相關性，這可能與環(huán)境、氣候等因素有關[37]，通過本研究對貴州不同地區(qū)酒曲微生物種群結構有了初步的了解，也為酒曲微生物多樣性探究進一步提供參考。