基于高通量测序技术解析高温大曲和中高温大曲的真菌群落结构

李静心*,王艳丽,何宏魁,刘国英,李晓欢,陈斌,吴翠芳,王录

(安徽古井贡酒股份有限公司,安徽 亳州,236820)

采用高通量测序技术对白酒高温大曲和中高温大曲进行分析,比较2种大曲真菌群落结构的差异。结果表明,霉菌是构成高温大曲的主要真菌,其中疏绵状嗜热丝孢菌(Thermomyces lanuginosus)和Wallemia sp.为高温大曲所特有,分别占总水平的80.49%和4.25%;酵母菌是构成中高温大曲的主要真菌,中高温大曲中霉菌含量仅有1/4左右,东方伊萨酵母(Issatchenkia orientalis)、扣囊复膜酵母(Saccharomycopsis fibuligera)、米曲霉(Aspergillus oryzae)、热带假丝酵母(Candida tropicalis)、阿氏丝孢酵母(Trichosporon asahii)、异常威克汉姆酵母(Wickerhamomyces anomalus)、微小根毛霉(Rhizomucor pusillus)为中高温大曲所特有,分别占总水平的38.22%、25.66%、21.40%、3.82%、2.75%、1.5%和1.29%。此外,Aspergillus intermedius在高温大曲和中高温大曲中的含量分别为10.95%和3.45%。该研究初步揭示了高温大曲和中高温大曲真菌群落结构的差异,为大曲在酿造生产中的应用提供参考。

关键词 白酒;大曲;高通量测序;真菌群落结构

DOI:10.13995/j.cnki.11-1802/ts.018043

引用格式李静心,王艳丽,何宏魁,等.基于高通量测序技术解析高温大曲和中高温大曲的真菌群落结构[J].食品与发酵工业,2018,44(12):52-59. LI Jing-xin, WANG Yan-li, HE Hong-kui,et al.High-throughput sequencing revealed fungal community structures at high temperature Daqu and medium temperature Daqu[J].Food and Fermentation Industries,2018,44(12):52-59.

第一作者:硕士研究生,工程师(本文通讯作者,E-mail:xx1516217@163.com)。

收稿日期:2018-06-15, 改回日期:2018-07-04

中国白酒是世界六大蒸馏酒之一,具有悠久的生产历史,并通过独特的发酵工艺生产,其中浓香型白酒具有“绵柔甘洌、芳香浓郁、香味协调、尾净余长、入口甜、落口绵”的酒体风格,在我国白酒行业中占据主导地位[1]。大曲作为酿酒过程中的糖化剂和发酵剂在浓香型白酒的酿造过程中起着糖化、生香、发酵、提供菌源的作用,被人们评价“曲为酒之骨”。大曲中含有丰富的微生物种群,其中霉菌可以产生淀粉酶、糖化酶、蛋白酶等丰富的酶系,酵母则具有酒精发酵能力及产酯能力,对白酒中酒精及香味物质的产生起着十分重要的作用[2]

对大曲中真菌微生物的研究最早是通过传统的分离培养方法,罗惠波等[3]对泸州老窖大曲中的真菌进行分离培养发现,大曲内存在大量的曲霉属、根毛霉属、横梗霉属等霉菌类群以及酿酒酵母、德巴利氏酵母属、异常威克汉姆酵母、锁掷酵母属等酵母菌类群。刘婷婷等[4]从白云边大曲及小麦原料中分离得到壮观拟青霉(Paecilomyces spectabilis),伞枝犁头霉(Absidia corymbifera)和微小根毛霉(Rhizomucor pusillus)。然而,传统的分离培养方法受到培养基的选择性、微生物的生长特性等因素影响[5],很难直接展现大曲中微生物的真实群落结构。随着分子生物学技术的发展,针对真菌内转录间隔区(internal transcribed spacerITS)研究技术的出现,大大提升了大曲微生物的研究效率。WANG[6]借助PCR-DGGE技术分析了不同大曲中微生物组成的差异,结果表明Saccharomycopsis fibuligeraPichia anomal是大曲中主要的酵母类群,Aspergillus oryzaeAbsidia blakesleeana是大曲中主要的霉菌类群。乔晓梅等[7]利用高通量测序技术对清香型大曲真菌群落结构进行分析,结果表明伞枝梨头霉(Lichtheimia corymbifera)、米根霉(Rhizopus oryzae)、库德毕赤酵母(Pichia kudriavzevii)和扣囊复膜酵母(Saccharomycopsis fibuligera)是清香型大曲中的优势菌。胡佳音等[8]对清香、酱香、浓香3种大曲真菌进行多样性分析,结果表明扣囊复膜酵母(Saccharomycopsis fibuligera)和伞枝横梗霉(Lichtheimia corymbifera)是清香型大曲的主要真菌,谢瓦散囊菌(Eurotium chevalieri)、米曲霉(Aspergillus oryzae)、毕赤酵母(Pichia kudriavzevii)、嗜热子囊菌(Thermomyces lanugin)是浓香型大曲主要真菌,谢瓦散囊菌(Eurotium chevalieri)和嗜热子囊菌(Thermomyces lanuginosus)是酱香型大曲的主要真菌。

该研究拟通过高通量测序技术,对白酒酿造所需高温大曲和中高温大曲进行研究,分析大曲中的主要真菌微生物构成,为优化浓香型白酒大曲生产工艺和提升原酒品质提供理论指导。

1 材料与方法

1.1 样品采集

大曲样品均采集自安徽省北部某知名白酒生产企业,样品为出房高温大曲和中高温大曲,每种大曲采集5个样品,置于-80℃保存。样品编号及分组如下:CZ1、CZ2、CZ3、CZ4、CZ5(中高温大曲);CG1、CG2、CG3、CG4、CG5(高温大曲)。

1.2 试剂与仪器

高级化学发光成像工作站,美国UVP;G-26C高速冷冻离心机,德国Sartorius;DYY-4C电泳仪,北京六一生物科技有限公司;MaestroNano微量分光光度计,美国MaestroGen。

DNeasy PowerSoil Kit土壤DNA提取试剂盒,美国QIAGEN;50×TAE溶液,生工生物工程(上海)股份有限公司;Agarose M,生工生物工程(上海)股份有限公司。

1.3 样品处理方法

1.3.1 DNA提取

分别称取0.25 g高温大曲和中高温大曲样品置于2.0 mL离心管中,参照DNeasy PowerSoil Kit提取大曲微生物的基因组DNA。琼脂糖凝胶电泳检测所提DNA的完整性,MaestroNano微量分光光度计测定DNA浓度及纯度。将符合测序要求的DNA -80℃保存。

1.3.2 高通量测序

利用深圳华大基因科技有限公司Illumina HiSeq 2500 2×250 bp测序平台对真菌ITS1进行高通量测序。

1.3.3 数据分析

根据华大基因提供的高通量测序下机数据(raw data),利用mothur软件进行拼接、过滤、嵌合体去除、OTU划分和物种注释[9],注释采用UNITE ITS数据库(Version No. 7.2, https://unite.ut.ee),置信度阈值设置为0.8。同时使用mothur软件对各样品的α多样性指数等指标进行计算。

2 结果与分析

2.1 高通量测序结果及真菌多样性分析

大曲样品提取基因组DNA后,对ITS1区进行高通量测序,10个样品共得到1 665 390条原始序列,质控后得到了1 549 585条优质序列,平均长度304 bp。将所有的序列按97% 的一致性将序列聚类成OTU,共得到200个OTU。使用UNITE真菌ITS数据库对每个OTU的代表序列进行物种注释,共注释得到3个门、16个纲、31个目、51个科、73个属和96个种水平的真核微生物。

对测序序列进行随机二次抽样,预测每个样品预期的OTU总数和相对丰度绘制稀释曲线[10](图1)。结果表明,随着测序序列的增加,10个样品的稀释曲线均趋于平滑,观测到的OTU数量趋于饱和,说明本试验的测序深度覆盖了绝大多数真菌微生物,本次测序结果能够真实反映中高温大曲和高温大曲真菌的群落结构和多样性。

图1 大曲真菌OTU数稀释曲线
Fig.1 Rarefaction curve of OTUs in Daqu

利用mothur软件计算各样品的α多样性指数(表1)。

表1 大曲的真菌多样性指数
Table 1 Fungus diversity index of Daqu

样品序列数目CoverageOTUChaoAceShannonSimpsonCG1157 0940.999 86298.11163.620.653 50.725 1CG2136 6800.999 93641.0044.290.337 40.868 7CG3156 2520.999 94159.2068.470.519 80.723 6CG4154 2530.999 96584.0099.190.582 70.765 6CG5154 2970.999 9104110.18112.161.373 50.364 1CZ1159 0840.999 98899.40106.161.681 70.237 1CZ2160 9520.999 97484.0089.251.393 10.308 5CZ3159 1940.999 96978.4376.991.047 90.537 5CZ4157 1130.999 998109.33112.041.949 80.189 2CZ5154 5500.999 98084.8892.061.613 40.278 3

注:Coverage-样品的序列覆盖度;OTU-样品的OTU数目;Chao指数和Ace指数-反映样品的群落丰富度;Shannon指数和Simpson指数-反映样品的群落多样性。

结果表明,所有样品的序列覆盖度(coverage)均大于0.999,再次证明本次测序的结果能够真实的反映所有样本真菌微生物的菌群结构。经分析,中高温大曲的OTU数量(81.80±5.14)多于高温大曲的OTU数量(61.60±12.02),但差异不显著(p≥0.05);中高温大曲的Chao指数(91.21±5.70)高于高温大曲的Chao指数(78.50±12.64),但差异不显著(p≥0.05),说明中高温大曲和高温大曲都具有较多的真菌种类,中高温大曲中的真菌种类稍多于高温大曲;中高温大曲的Shannon指数(1.537±0.151)高于高温大曲的Shannon指数(0.693±0.178),且差异极显著(p≤0.01),说明中高温大曲的真菌群落多样性高于高温大曲。综上所述,中高温大曲和高温大曲在真菌群落丰富度上差异较小,但中高温大曲的真菌群落多样性远高于高温大曲。

2.2 大曲真菌群落组成分析

2.2.1 属水平分析

使用mothur软件按照97%的相似度对测序结果进行OTU划分,采用UNITE ITS数据库(Version No. 7.2)对每个OTU进行classify,中高温大曲和高温大曲共得到73个属的真菌,其中核心真菌微生物(含量>1%)共有9个属(表2),分别是Aspergillus(曲霉属)、Thermomyces(嗜热真菌属)、Wickerhamomyces(威克汉姆酵母属)、Issatchenkia(伊萨酵母属)、Candida(假丝酵母属)、Saccharomycopsis(复膜孢酵母属)、Trichosporon(毛孢子菌属)、Wallemia属和Rhizomucor(根毛霉属),根据其在每个样品所占比例绘制柱状分布图(图2)。

表2 大曲真菌属水平百分比组成 单位:%

Table 2 Percentage of microbial community at genus levels in Daqu

真菌属名CG1CG2CG3CG4CG5CZ1CZ2CZ3CZ4CZ5Aspergillus11.483.2515.038.6425.0930.7132.2711.4427.5122.32Thermomyces84.6993.1183.7787.2153.700.080.050.010.000.01Wickerhamomyces0.000.000.000.010.301.290.220.392.103.51Issatchenkia0.622.810.200.711.6926.3341.5971.5627.0324.57Candida0.010.010.010.020.063.134.031.568.581.79Saccharomycopsis0.270.020.070.311.5929.9219.0211.6626.1341.59Trichosporon0.010.000.000.040.216.321.340.863.881.36Wallemia2.620.310.422.7215.190.010.000.000.020.03Rhizomucor0.020.020.020.140.080.810.721.011.812.11Others0.270.470.470.202.101.410.741.522.942.71

结果表明,高温大曲的优势真菌为Aspergillus((12.70±3.65)%)和Thermomyces((80.50±6.89)%),约占高温大曲真菌的78.78%~98.80%;中高温大曲的优势真菌为Aspergillus((24.85±3.76)%)、Issatchenkia((38.22±8.88)%)、Saccharomycopsis((25.66±5.06)%),约占中高温大曲真菌的80.67%~94.65%。

图2 大曲真菌属水平微生物菌群
Fig.2 Microbial community at genus levels in Daqu

2.2.2 种水平分析

在种水平对大曲真菌进行分类分析,共得到了96个种,其中优势真菌(>1%)有10个种(表3),分别是Aspergillus oryzae(米曲霉)、Aspergillus intermediusThermomyces lanuginosus(疏绵状嗜热丝孢菌)、Wickerhamomyces anomalus(异常威克汉姆酵母)、Issatchenkia orientalis(东方伊萨酵母)、Candida tropicalis(热带假丝酵母)、Saccharomycopsis fibuligera(扣囊复膜酵母)、Trichosporon asahii(阿氏丝孢酵母)、Wallemia sp.和Rhizomucor pusillus(微小根毛霉),根据其在每个样品所占比例绘制柱状分布图(图3)。

结果表明,高温大曲的优势真菌为2种霉菌Thermomyces lanuginosus((80.49±6.90)%)和Aspergillus intermedius((10.95±3.53)%),约占高温大曲真菌的74.46%~98.58%;中高温大曲的优势真菌为2种酵母菌Issatchenkia orientalis((38.22±8.88)%)和Saccharomycopsis fibuligera((25.66±5.06)%),以及1种霉菌Aspergillus oryzae((21.40±3.59)%),约占大曲真菌的70.39%~93.95%。

对高温大曲和中高温大曲优势真菌中霉菌与酵母菌的相对丰度进行分析,并绘制相对丰度柱形图(图4)。结果表明,高温大曲的霉菌含量较高,约占真菌的78.65%~98.73%,酵母菌含量较低,仅有0.29%~3.84%;中高温大曲的酵母菌所占比例约为66.21%~86.02%,霉菌所占比例稍低,约为12.46%~33.04%(表4)。霉菌是构成高温大曲的主要真菌,而酵母菌是构成中高温大曲的主要真菌。

表3 大曲真菌种水平百分比组成 单位:%

Table 3 Percentage of microbial community at species levels in Daqu

CG1CG2CG3CG4CG5CZ1CZ2CZ3CZ4CZ5Aspergillus oryzae3.370.110.142.482.1126.8931.2510.7417.2320.88Aspergillus intermedius7.963.1114.816.0822.773.821.020.7010.271.44Thermomyces lanugino-sus84.6993.1183.7787.2153.690.070.050.010.000.01Wickerhamomyces anom-alus0.000.000.000.010.301.290.220.392.103.51Issatchenkia orientalis0.622.810.200.711.6926.3341.5971.5627.0324.57Candida tropicalis0.010.010.010.020.063.134.031.568.581.79Saccharomycopsis fibu-ligera0.270.020.070.311.5929.9219.0211.6626.1341.59Trichosporon asahii0.010.000.000.040.216.321.340.863.881.35Wallemia sp.2.620.310.422.7215.190.010.000.000.020.03Rhizomucor pusillus0.020.020.020.130.080.810.721.011.812.11Others0.420.500.560.292.311.410.751.522.952.72

图3 大曲真菌种水平微生物菌群
Fig.3 Microbial community at species levels in Daqu

图4 霉菌与酵母菌含量柱形图
Fig.4 Contents of moulds and yeasts in Daqu

表4 大曲霉菌与酵母菌百分比组成单位:%

Table 4 Percentage of moulds and yeasts in Daqu

CG1CG2CG3CG4CG5CZ1CZ2CZ3CZ4CZ5霉菌96.0596.3598.7395.9078.6531.5933.0412.4629.3224.43酵母0.912.840.291.093.8466.9966.2186.0267.7272.82

2.3 物种差异分析

根据样品中不同种水平真菌含量的相对丰度,采用Bray-Curtis算法计算样本间距离矩阵,通过PCoA(principal coordinates analysis)分析高温大曲和中高温大曲群体间的差异。两组大曲真菌群落结构PCoA分析结果如图5所示,第一主成分PC1的贡献率为79.55%,权重最大,第二主成分PC2的贡献率为7.25%,累计贡献率为86.80%,说明两组间差异具有统计学意义。高温大曲和中高温大曲可明显的划分为两组,中高温大曲样品CZ2和CZ3虽然偏离较远,但根据主成分一(79.55%)可与其余样品聚为一类,说明高温大曲和中高温大曲真菌群落结构具有明显的不同。

图5 大曲PCoA分析图
Fig.5 PCoA analysis of Daqu

为了找出不同大曲中呈现显著性差异的物种,对大曲中的优势真菌进行LEfSe(LDA effect size)分析,通过LDA(linear discriminative analysis)值来估算每个物种的丰度对组间差异效果影响的大小[11-12],进一步明确高温大曲和中高温大曲独特的微生物类群(图6)。

图6 大曲微生物组LEfSe分析柱状图
Fig.6 Histogram based on LEfSe analyses of microbiome in Daqu

LEfSe分析柱状图展现了高温大曲和中高温大曲中具有显著性差异的真菌微生物,柱的长度代表差异物种的影响值大小。结果表明,高温大曲的Thermomyces lanuginosusWallemia sp.与中高温大曲具有显著性差异,中高温大曲的Issatchenkia orientalisSaccharomycopsis fibuligeraAspergillus oryzaeCandida tropicalisTrichosporon asahiiWickerhamomyces anomalusRhizomucor pusillus与高温大曲具有显著性差异。

疏绵状嗜热丝孢菌(Thermomyces lanuginosus)是一种分布广泛的、能够在高温条件下生存的特殊真菌,其最低生长温度高于20 ℃,在低于19 ℃的环境中无法生长,而最高生长温度可达50 ℃以上[13-14]。它能产生具有高热稳定性和活力的纤维素酶、蛋白酶、几丁质酶、β-葡聚糖酶等,有利于酿造过程中产酒生香[15-19]

Wallemia属是一种广泛分布的耐高渗透压真菌,包含耐旱、嗜热和嗜盐的物种,在各种环境中均有发现,如干燥、腌渍或高糖食物、干饲料、晒盐场、盐晶体、室内外空气中[20]。目前该属共有8个物种,国内最早于2006年报道从土壤中分离出1株Wallemia sebi,其在马铃薯葡萄糖琼脂培养基上对植物病原菌尖孢镰刀菌(Fusarium oxysporum)以及大丽花轮枝菌(Verticillium dahliae)具有很强的抑制作用[21]。在高温大曲中发现的Wallemia sp.属于未知的Wallemia物种,且含量不高((4.25±2.78)%),推测其是1种嗜热的Wallemia真菌,其在大曲中的作用需要进一步研究。

东方伊萨酵母(Issatchenkia orientalis)在汾酒[22]、白云边酒[23-24]和葡萄酒[25-26]的酿造过程中均有发现,是一种广泛存在于酿酒生产中的能够产生乙醇的酵母菌[27]。刘婷婷[23-24]在白云边酒醅中分离得到了1株耐高温、耐高酒精度且产乙酸乙酯能力强的东方伊萨酵母菌株,具有较强的产乙醇能力,是白云边酒酿造过程中的优势酵母菌株。LI等[22]对汾酒发酵生产过程中真菌多样性进行分析发现,东方伊萨酵母是汾酒发酵过程中含量最高的酵母。该研究发现,中高温大曲中东方伊萨酵母的含量最高(38.22%),表明中高温大曲可能是发酵过程中东方伊萨酵母的主要来源。

扣囊复膜酵母菌(Saccharomycopsis fibuligera)广泛分布于自然界中,能够分泌α-淀粉酶[28-29]、蛋白酶[30]在内的多种水解酶,在发酵行业中具有广泛的应用[31]。目前已从清香型大曲[32]、浓香型大曲[33]、酱香型大曲[34]、糯米酒曲[35]和黄酒曲[36]中分离得到该种酵母菌,本研究也在中高温大曲中发现了扣囊复膜酵母且含量较高(25.66%)。扣囊复膜酵母不但能够在窖池中发酵产生乙醇[29],还可以提高大曲酒中多种醇和酯的含量[37-38],使原酒品质有所提高。除东方伊萨酵母和扣囊复膜酵母菌外,中高温大曲中的热带假丝酵母(Candida tropicalis)、异常威克汉姆酵母(Wickerhamomyces anomalus)和阿氏丝孢酵母(Trichosporon asahii)也具有提高乙酸乙酯产量[39-40]、分泌β-葡萄糖苷酶[41]的能力,从而使得原酒的风味更加协调、丰满。

米曲霉(Aspergillus oryzae)是食品酿造中的常用菌株,尤其在酱油的酿造过程中的应用最为广泛[42-43],其安全性能已得到了广泛认可[44-46]。米曲霉具有酶系丰富,酶活高,适应能力强,不产毒素等特点[47],产淀粉酶[48-49]、蛋白酶[50-51]、脂肪酶[52]、纤维素酶[53]和果胶酶[54]的能力非常强。本研究发现米曲霉在中高温大曲中的含量较高(21.40%),是中高温大曲的主要霉菌,其较强的产淀粉酶能力有助于提高大曲的糖化力和液化力,而产生的脂肪酶则对酿造过程中乙酸乙酯、丁酸乙酯和己酸乙酯等骨架风味物质的形成有一定的促进作用。此外,虽然微小根毛霉(Rhizomucor pusillus)在中高温大曲中的含量较低,仅有(1.29±0.28)%,但其产生的α-淀粉酶则具有较好的热稳定性和pH耐受性[55]

3 结论与讨论

高温大曲和中高温大曲是以制曲过程中能达到的最高温度来区分的,高温大曲主要用于酱香型白酒的生产,其生产过程中最高温度约为60~70 ℃,中高温大曲主要用于浓香型白酒的生产,其生产过程中最高温度约为50~60 ℃,此外生产过程中最高温度低于50℃的低温大曲则主要用于清香型白酒的生产[56]。由于大曲生产过程中原料、环境和生产工艺的差异,不同大曲拥有不同的微生物菌系和酶系,使得糖化力、液化力和酯化力等理化指标有所差异。张春林[57]对泸州老窖生产的中、高温大曲进行研究发现,中高温大曲液化酶、糖化酶、蛋白酶和纤维素酶的活性以及酯化力、产酒能力均高于高温大曲。邢钢等[58]对3种不同温度大曲制曲过程的理化指标进行分析发现,高温大曲的水分、淀粉、酸度和中温大曲相差不大,但糖化力、液化力、酯化力和发酵力均低于中温大曲。邢爽等[59]针对3种大曲的酯化力进行分析,结果表明低温大曲合成乙酸乙酯和丁酸乙酯的酶活最高,中高温大曲合成乳酸乙酯和己酸乙酯的酶活最高,而高温大曲的酯化酶活力最小。

为了进一步分析高温大曲和中高温大曲菌系与酶系的不同,本研究借助高通量测序技术对2种大曲进行分析,找出了两者真菌微生物群落结构的差异。结果表明高温大曲和中高温大曲所含OTU数量接近,差异不显著,但中高温大曲真菌的群落多样性显著高于高温大曲;霉菌是构成高温大曲的主要真菌,其中疏绵状嗜热丝孢菌(Thermomyces lanuginosus)和Wallemia sp.为高温大曲所特有,Aspergillus intermedius也是高温大曲中的优势菌,含量在10.95%左右;酵母菌则是构成中高温大曲的主要真菌,霉菌仅占真菌含量的1/4左右,东方伊萨酵母(Issatchenkia orientalis)、扣囊复膜酵母(Saccharomycopsis fibuligera)、米曲霉(Aspergillus oryzae)、热带假丝酵母(Candida tropicalis)、阿氏丝孢酵母(Trichosporon asahii)、异常威克汉姆酵母(Wickerhamomyces anomalus)、微小根毛霉(Rhizomucor pusillus)为中高温大曲所特有。大量的酵母菌保证了中高温大曲具有较高的发酵力、液化力、酯化力及产酒能力,而米曲霉则保证了中高温大曲液化酶、糖化酶、蛋白酶和纤维素酶的活性;高温大曲在生产过程中最高温度超过60 ℃,酵母菌、米曲霉等不耐高温的微生物难以生长,使得高温大曲发酵力、液化力、酯化力及产酒能力较低,而耐高温疏绵状嗜热丝孢菌分泌的高热稳定性酶使得高温大曲具有独特的生香作用。

胡佳音等[8]采用高通量测序对清浓酱3种大曲真菌多样性进行分析发现,谢瓦散囊菌(Eurotium chevalieri)和疏绵状嗜热丝孢菌(Thermomyces lanuginosus)为高温大曲中的主要真菌,分别占43%和34.8%;谢瓦散囊菌(Eurotium chevalieri)、米曲霉(Aspergillus oryzae)、毕赤酵母(Pichia kudriavzevii)和嗜热子囊菌(Thermomyces lanugin)为中高温大曲的主要真菌,分别占52%、11.7%、8.9%和6.2%。徐占成等[60]对川派浓香型白酒中高温大曲进行高通量测序,结果表明RhizopusSaccharomycopsisAspergillusThermoascus等属为剑南春大曲的优势真菌。李申奥[61]对湖北白云边高温大曲真菌群落结构进行解析发现,Thermomyces属、Candida属和Aspergillus属为出房高温大曲的优势真菌,其占比分别为77.89%、12.60%和8.53%。本文所分析的大曲均采集于皖北某知名酒企,该企业是皖派白酒的代表企业,对其生产的高温大曲和中高温大曲进行分析能够反映皖派白酒大曲的真菌群落结构。与其他产区相比,该地区高温大曲的真菌群落结构与白云边高温大曲类似,嗜热真菌属(Thermomyces)均占据主要地位,与胡佳音[8]的研究结果相比并未发现谢瓦散囊菌(Eurotium chevalieri)。中高温大曲中曲霉属(Aspergillus)和伊萨酵母属(Issatchenkia)含量远高于剑南春中高温大曲,而根霉属(Rhizopus)含量远低于剑南春中高温大曲,也未发现胡佳音[8]研究结果中占主要地位的谢瓦散囊菌(Eurotium chevalieri)。

大曲中含有丰富的菌系和酶系,通过高通量测序技术对真菌群落结构进行分析明确了大曲中的微生物种类和相对丰度,了解大曲中主要优势菌的演变规律,更有效地帮助白酒企业认知大曲在酿酒过程中的作用。

参考文献

[1] 邓杰,黄治国,卫春会,等.基于高通量测序的浓香型白酒窖池细菌群落结构分析[J]. 现代食品科技, 2015(7):50-55.

[2] 罗惠波,侯海波,黄治国,等. 大曲真核微生物群落PCR-DGGE电泳条件优化[J]. 四川理工学院学报(自科版), 2011, 24(5):515-518.

[3] 罗惠波,杨晓东,杨跃寰,等. 浓香型大曲中可培养真菌的分离鉴定与系统发育学分析[J]. 现代食品科技, 2013(9):2 047-2 052.

[4] 刘婷婷,张明春,曾驰,等. 白云边酒大曲及小麦原料中主要微生物的分析[J]. 中国酿造, 2010, 29(11):32-35.

[5] 赵东,郑佳,彭志云,等. 高通量测序技术解析五粮液窖泥原核微生物群落结构[J]. 食品与发酵工业, 2017, 43(9):1-8.

[6] WANG H Y, GAO Y B, FAN Q W, et al. Characterization and comparison of microbial community of different typical Chinese liquor Daqus by PCR-DGGE.[J]. Letters in Applied Microbiology, 2011, 53(2):134-140.

[7] 乔晓梅,赵景龙,杜小威,等. 高通量测序法对清香大曲真菌群落结构的分析[J]. 酿酒科技, 2015,250(4):28-31.

[8] 胡佳音,周森,赵卫鹏,等. 清、浓、酱三种大曲真菌多样性初步分析[J]. 酿酒科技, 2016,266(8):87-90.

[9] SCHLOSS P D, WESTCOTT S L, RYABIN T, et al. Introducing mother: Open-source, platform-independent, community-supported software for describing and comparing microbial communities[J]. Applied & Environmental Microbiology, 2009, 75(23):7 537-7541.

[10] KEMP P F, ALLER J Y. Bacterial diversity in aquatic and other environments: what 16S rDNA libraries can tell us[J]. Fems Microbiology Ecology, 2004, 47(2):161-177.

[11] ZHANG C H, LI S F, YANG L, et al. Structural modulation of gut microbiota in life-long calorie-restricted mice[J]. Nature Communications, 2013, 4(1):2 163-2 173.

[12] SEGATA N, IZARD J, WALDRON L, et al. Metagenomic biomarker discovery and explanation[J]. Genome Biology, 2011, 12(6):R60-R78.

[13] 蔡杰华,刘斌,石冬梅,等. 四种嗜热真菌的分离与鉴定[J]. 现代生物医学进展, 2012, 12(6):1 059-1 064.

[14] 董瑞丽,罗惠波,叶光斌. 浓香型大曲中的嗜热真菌[J]. 中国酿造, 2011, 30(4):75-77.

[15] SINGH S, MSDLALA A M, PRIOR B A. Thermomyces lanuginosus: Properties of strains and their hemicellulases[J]. Fems Microbiology Reviews, 2003, 27(1):3-16.

[16] SINGH S, DU-PREEZ J C, PILLAY B, et al. The production of hemicellulases by Thermomyces lanuginosus strain SSBP: Influence of agitation and dissolved oxygen tension[J]. Applied Microbiology and Biotechnology, 2000, 54(5): 698-704.

[17] 王瑞杰, NOKUTHULA P M,牛丹丹,等. 疏绵状嗜热丝孢菌外切β-葡聚糖酶的基因克隆与酶学特征[J]. 微生物学通报, 2016, 43(2):285-291.

[18] 郭润芳,李多川,王荣. 疏绵状嗜热丝孢菌热稳定几丁质酶的纯化及其性质研究[J]. 微生物学报, 2005, 45(2):270-274.

[19] 夏玙,罗惠波,周平,等.不同处理方式的大曲真菌群落差异分析[J].食品科学,2018,39(22):166-172.

[20] EAJC J, GUNDECIMERMAN N. The genus Wallemia—from contamination of food to health threat[J]. Microorganisms, 2018; 6(2):46-56.

[21] 郭永霞,赵爽,刘庆洪,等. Wallemia sebi的RNase活性和抗植物病原真菌活性的初步研究[J]. 微生物学通报, 2006, 33(2):30-33.

[22] LI X R, MA E B, YAN L Z, et al. Bacterial and fungal diversity in the traditional Chinese liquor fermentation process[J]. International Journal of Food Microbiology, 2011, 146(1):31-37.

[23] 刘婷婷,曾驰,杨团元,等. 白云边酒优势酵母菌的分离鉴定及其发酵特征分析[J]. 中国酿造, 2013, 32(12):116-120.

[24] 刘婷婷. 白云边酒酿造微生物分析及东方伊萨酵母发酵特性研究[D]. 武汉:武汉工业学院, 2011.

[25] CLAVIJO A, CALDERON I L, PANEQUE P. Diversity of Saccharomyces and non-Saccharomyces yeasts in three red grape varieties cultured in the Serranía de Ronda, (Spain) vine-growing region[J]. International Journal of Food Microbiology, 2010, 143(3):241-245.

[26] LI S S, CHENG C, LI Z, et al. Yeast species associated with wine grapes in China[J]. International Journal of Food Microbiology, 2010, 138(1):85-90.

[27] GALLARDO J C M, SOUZA C S, CICARELLI R M B, et al. Enrichment of a continuous culture of Saccharomyces cerevisiae, with the yeast Issatchenkia orientalis, in the production of ethanol at increasing temperatures[J]. Journal of Industrial Microbiology & Biotechnology, 2011, 38(3):405-414.

[28] 郭玉婉. 酵母α-淀粉酶基因克隆、表达及其应用于麦芽三糖的制备[D]. 无锡:江南大学, 2013.

[29] 陈蕾. 扣囊复膜酵母A11菌株淀粉酶的生产和酒精发酵的研究[D]. 青岛:中国海洋大学, 2010.

[30] 于欣君. 扣囊复膜酵母菌酸性蛋白酶基因在解脂耶罗维亚酵母菌中表达和重组酸性蛋白酶的应用[D]. 青岛:中国海洋大学, 2011.

[31] 苏畅,马莹莹,杨建刚. 扣囊复膜酵母在酿酒中的应用研究进展[J]. 食品研究与开发, 2018(1):205-209.

[32] 马凯,崔哲男,郑晓卫,等. 汾酒大曲可培养真菌多样性的初步分析[J].中国酿造,2011,233(8):19-21.

[33] 明红梅,郭志,周健,等.浓香型大曲中产香微生物的筛选及鉴定[J].现代食品科技,2015,31(4):186-191.

[34] 吴徐建. 酱香型白酒固态发酵过程中酵母与细菌群落结构变化规律的研究[D].无锡:江南大学,2013.

[35] 应玲云,伍时华,赵东玲,等. 烧酒曲中扣囊复膜酵母的分离及鉴定[J].食品与发酵工业,2013,39(1):146-150.

[36] 臧威,谢广发,孙剑秋,等. 绍兴黄酒酒药中酵母菌的物种资源[J].菌物学报,2015,34(6):1 078-1 084.

[37] 罗小叶,邱树毅,陆安谋,等. 酱香大曲产香酵母的分离及鉴定[J].食品与发酵工业,2016,42(12):26-31.

[38] 王晓丹,陈美竹,班世栋,等. 茅台大曲中酵母的分离、鉴定及其功能初探[J]. 食品科学, 2017, 38(4):51-57.

[39] 王鹏昊,关统伟,张习超,等. 小曲中优质产酯酵母分离鉴定及其产酯条件优化[J]. 食品与发酵工业, 2018,44(2):62-67.

[40] 侯建光,郭富祥,樊建辉,等. 响应面法优化陶香型窖池中异常威克汉姆酵母发酵产乙酸乙酯条件[J]. 酿酒科技, 2016,268(10):77-81.

[41] 王玉霞. 阿氏丝孢酵母(Trichosporon asahii)β-葡萄糖苷酶及葡萄糖苷类风味物质水解机制的研究[D]. 无锡:江南大学, 2012.

[42] 李幼筠,周逦. 科学利用微生物推动中国酱油工艺大变革[J]. 中国酿造, 2011, 30(1):1-5.

[43] 包启安. 酱油科学与酿造技术[M]. 北京:中国轻工业出版社, 2011.

[44] TAYLOR M J, RICHARDSON T. Applications of microbial enzymes in food systems and in biotechnology [J]. Advances in Applied Microbiology, 1979, 25:7-35.

[45] ABE K, GOMI K, HASEGAWA F, et al. Impact of Aspergillus oryzae, genomics on industrial production of metabolites[J]. Mycopathologia, 2006, 162(3):143-155.

[46] MACHIDA M, YAMADA O, GOMI K. Genomics of Aspergillus oryzae: Learning from the history of Koji mold and exploration of its future[J]. Dna Research, 2008, 15(4):173-183.

[47] 赵中开,龙可,马莹莹,等. 米曲霉菌在酿酒工业中的研究进展[J]. 现代食品科技, 2013,29(4):932-935.

[48] 王璐,曹钰,陆健,等. 绍兴黄酒麦曲中α-淀粉酶的初步研究[J]. 中国酿造, 2007, 26(5):38-40.

[49] 刘晓蓉,连晓蔚,朱美娟. 米曲霉1228产α-淀粉酶特性的研究[J]. 中国调味品, 2009, 34(4):36-38.

[50] 刘志伟,谭兴和,周红丽,等. 米曲霉产中性蛋白酶条件的优化[J]. 中国酿造, 2011, 30(11):103-107.

[51] 邓静,徐静,吴华昌,等. 米曲霉高产酸性蛋白酶菌株的选育[J]. 中国调味品, 2010, 35(1):53-55.

[52] 王小花,洪枫,朱利民,等. 米曲霉产胞外脂肪酶培养条件的优化[J]. 食品与发酵工业, 2005, 31(4):25-28.

[53] 黄晓东,叶生梅,高新义,等. 米曲霉固态发酵豆渣及其提取物抗氧化性的研究[J]. 食品科技, 2005,82(5):80-82.

[54] 汤鸣强,卞矛,李琼华. 米曲霉固体发酵生产果胶酶的研究[J]. 福建师大福清分校学报, 2006(2):31-35.

[55] 何正贵. Rhizomucor pusillus耐热淀粉酶基因的克隆及其在毕赤酵母中的表达[D]. 上海:华东理工大学, 2013.

[56] 敖宗华,陕小虎,沈才洪,等. 国内主要大曲相关标准及研究进展[J]. 酿酒科技, 2010,188(2):104-108.

[57] 张春林. 泸州老窖大曲的质量、微生物与香气成分关系[D]. 无锡:江南大学, 2012.

[58] 邢钢,敖宗华,王松涛,等. 不同温度大曲制曲过程理化指标变化分析研究[J]. 酿酒科技, 2014,240(6):20-23.

[59] 邢爽,魏志阳,郭学武,等.不同香型大曲酯化酶的特性分析[J]. 现代食品科技,2018(05): 1-5.[2018-06-06]. http://kns.cnki.net/kcms/detail/44.1620.TS.20180411.1243.050.html.

[60] 徐占成,徐姿静,刘孟华,等. 高通量测序法对剑南春大曲真菌群落结构的分析[J]. 酿酒科技, 2018,286(4): 22-25;32.

[61] 李申奥. 兼香型白酒高温大曲微生物群落演替规律的研究[D]. 武汉:华中农业大学, 2016.

High-throughput sequencing revealed fungal community structures at high temperature Daqu and medium temperature Daqu

LI Jing-xin*, WANG Yan-li, HE Hong-kui, LIU Guo-ying, LI Xiao-huan, CHEN Bin, WU Cui-fang, WANG Lu

(Anhui Gujing Gongjiu Co. Ltd., Bozhou 236820, China)

Abstract The high-throughput sequencing was used to compare the difference of fungal community structures between high temperature Daqu and medium temperature Daqu in Baijiu brewing. The results showed that the moulds were the main fungus of high temperature Daqu, Thermomyces lanuginosus and Wallemia sp. were unique species in high temperature Daqu, which accounted for 80.49% and 4.25%, respectively. In medium temperature Daqu, the moulds accounts for 1/4 of fungi and the yeast were the main fungi. Issatchenkia orientalis, Saccharomycopsis fibuligera, Aspergillus oryzae, Candida tropicalis, Trichosporon asahii, Wickerhamomyces anomalus and Rhizomucor pusillus were unique species in medium temperature Daqu, which accounted for 38.22%, 25.66%, 21.40%, 3.82%, 2.75%, 1.5%, and 1.29%, respectively. Besides, Aspergillus intermedius accounted for 10.95% in high temperature Daqu and 3.45% in medium temperature Daqu, respectively. This results revealed the difference of fungal community structures between high temperature Daqu and medium temperature Daqu and provided a reference for the application of Daqu in the brewing process.

Key words Baijiu; Daqu; high-throughput sequencing; fungal community structure