基于高通量测序技术解析中高温制曲真菌群落的演替规律研究.pdf

332023No.10(Tol.352)2023年第1 0 期（总第352 期）LIQUOR-MAKINGSCIENCE&TECHNOLOGY酿酒科技D0I:10.13746/j.njkj.2023203基于高通量测序技术解析中高温制曲真菌群落的演替规律研究杨勇,葛向阳1,2*,张聪芝,李燕荣，贾亚伟（1.江苏洋河酒厂股份有限公司，江苏宿迁2 2 38 0 0；2.宿迁学院，江苏宿迁2 2 38 0 0)摘要：以洋河中高温大曲为研究对象，通过高通量测序技术分析不同培养阶段真菌微生物群落结构及其演替规律，并探索了大曲发酵、生香等功能与真菌种类之间的潜在关系。研究结果表明，子囊菌门和毛霉门占据绝对优势，子囊菌门包括半子囊菌纲的酵母目和散囊菌纲的散囊菌目，毛霉门即毛霉纲的毛霉目；属水平下，嗜热真菌属、根毛霉属、根霉属、曲霉属、生丝毕赤酵母属、横梗霉属、威克汉姆酵母属等真菌分别在不同培养阶段占据主导地位，前期主要是酵母属，后期主要是耐高温霉菌。其中，大曲发酵力与生丝毕赤酵母属相关性较高，酯化力与嗜热真菌属、红曲霉、威克汉姆酵母属相关性较高，水分则与孢圆酵母属的相关性较高。研究结果揭示了中高温大曲发酵过程中真菌群落的演替规律，解释了大曲发酵、生香等功能与真菌种类之间的潜在关系，有助于大曲微生物多样性形成机制的进一步解析，为优质大曲生产提供了理论依据。关键词：中高温大曲；高通量测序；真菌；演替；相关性中图分类号：TS262.3;TS261.1文献标识码：A文章编号：1 0 0 1-92 8 6(2 0 2 3)1 0-0 0 33-0 6High-Throughput Sequencing Revealed Succession Law of FungiCommunity in Medium/High-Temperature DaquYANG Yong,GE Xiangyang2,ZHANG Congzhi,LI Yanrong and JIA Yawei(1.Yanghe Distillery Co.Ltd.,Suqian,Jiangsu 223800;2.Suqian University,Suqian,Jiangsu 223800,China)Abstract:Taking Yanghe medium/high-temperature Daqu as the research object,high-throughput sequencing technology was used toanalyze the structure and succession law of fungi communities at different cultivation stages.The potential relationship between Daqufermentation,aroma production and fungal species was explored.The results showed that Ascomycota and Mucoromycota were domi-nant.At the genus level,Thermomyces,Rhizomucor,Rhizopus,Aspergillus,Hyphopichia,Lichtheimia,and Wickerhamomyces weredominant respectively at different cultivation stages.In the early stages,yeasts were dominant;in the later stages,high-temperature-re-sistant molds were dominant.The fermentation power of Daqu was highly correlated with Hyphopichia.The esterification power ofDaqu was highly correlated with Thermomyces,Monascus and Wickerhamomyces.The water content of Daqu was highly correlatedwith Torulaspora.The results revealed the succession law of fungi communities during the fermentation process of medium/high-tem-perature Daqu,explained the potential relationship between Daqu fermentation,aroma production and fungal species,contributed tofurther analysis of the formation mechanism of microbial diversity of Daqu,and provided a theoretical basis for the production ofhigh-quality Daqu.Key words:medium/high-temperature Daqu;high-throughput sequencing;fungi;succession;correlation大曲作为一种富含多酶多菌的微生态制品，在中国白酒的酿造过程中起着糖化、生香、发酵、提供基金项目：2 0 2 1 年宿迁市科技计划（基础研究计划自然科学资金）项目（k202135）。作者简介：杨勇（1 98 4-），男，江苏宿迁人，硕士，高级工程师，研究方向为微生物学与白酒酿酒工艺，E-mail:。通讯作者：葛向阳（1 97 4-），男，江苏新沂人，博士，正高级工程师，目前主要研究酿酒工程及生物工程，E-。342023 No.10(Tol.352)酿酒科技2023年第1 0 期（总第352 期）LIQUOR-MAKINGSCIENCE&TECHNOLOGY菌源的作用，是传统固态发酵蒸馏大曲酒的重要物质保障，被人们评价“曲为酒之骨。对大曲中真菌微生物的研究最早是通过传统的分离培养方法，葛向阳等2 经过多级富集、分离、纯化与鉴定，从洋河传统包包曲中分离得到毕赤酵母属、异常汉逊酵母、假丝酵母等酵母菌类群以及卷枝毛霉、寄生曲霉、黑曲霉和黄曲霉等霉菌类群；罗惠波等3 发现泸州老窖大曲内存在大量的曲霉属、根毛霉属、横梗霉属等霉菌类群以及酿酒酵母、德巴利氏酵母属、异常威克汉姆酵母、锁掷酵母属等酵母菌类群。对中高温制曲发酵过程中可培养微生物的消长及生化指标的变化规律也做了相关研究(4-,然而，传统的分离培养方法受到培养基的选择性、微生物的生长特性等因素影响,很难直接展现大曲中微生物的真实群落结构。随着分子生物学的发展，基于免培养的高通量测序技术已经越来越广泛地应用于微生物群落的研究。谭崇尧等7 利用高通量测序发现北方派浓香型白酒大曲优势真菌有热子囊菌属、嗜热真菌属、威克汉姆酵母属等，长江中游派优势真菌为曲霉属、热子囊菌属、嗜热真菌属，江淮派优势真菌有曲霉属、热子囊菌属、耐干霉菌属等，川派优势真菌有热子囊菌属、嗜热真菌属、曲霉属等；吴树坤等1 8 利用高通量测序发现四川不同地区浓香型白酒大曲中优势真菌为嗜热真菌属、曲霉属、嗜热子囊菌属、丝衣霉属、节担菌属；施思等9利用高通量测序技术发现浓香型白酒大曲在贮存过程中的真菌群落结构不断调整，毕赤酵母属、根霉属及横梗霉属成为最终的优势菌群；Lio对不同工艺的中、低温大曲的微生物群落变化进行研究，发现在不同环境作用下，微生物群落结构会发生改变。基于高通量测序技术，已经对细菌群落演替规律进行了相关分析叫。本文以中高温大曲为研究对象，采用高通量测序技术对培菌发酵几个主要阶段的真菌群落进行分析，研究各阶段优势真菌的分布、丰度及其在培菌发酵过程中的演替规律，为进一步剖析制曲发酵机理提供可靠的微生物信息。1材料与方法1.1材料、试剂及仪器大曲样品：中高温大曲培养过程中入房（编号：RF，发酵0 天）、并房（编号：BF，发酵5天）、上架（编号：SJ，发酵9天）、大火（编号：DH，发酵1 5天）、下架（编号：XJ，发酵2 1 天）、出房（编号：CF，发酵35天）各个阶段的混合样，分别粉碎、混匀、装袋，于-8 0 保存备用。试剂耗材:TaKaRa Premix Taq?Version 2.0、E.Z.N.A.?Gel Extraction Kit(Omega,USA)、NEB-Next?Ultra TM II DNA Library Prep Kit for Illumi-na?(New England Biolabs,USA)。仪器设备：Novaseq6000PE250平台、美国Bio-Rad Laboratory S1000 型PCR仪、Thermo Nano-Drop One。1.2实验方法按照标准流程对大曲样品进行DNA提取，以ITS5-1737F和ITS2-2043R为引物进行PCR扩增，每个样本进行3个重复，并将同一样本的PCR产物进行混合，然后电泳检测、切胶纯化、建库测序。2结果与分析2.1测序有效性分析中高温制曲各阶段真菌的richness曲线如图1所示。结果显示，随着测序深度的增加，各样本的richness指数曲线都是先增加后逐步趋于平缓，说明测序数据量足够大，可以反映中高温制曲各阶段样本中的绝大多数物种信息。80RFBF70-SJDH60-XJCFSJ50RF403020-BF10020406080100图1大曲培养各阶段真菌的richness曲线2.2Alpha多样性分析对中高温制曲各阶段的真菌群落进行Alpha多样性分析，结果见表1。所有样本的coverage均35杨勇，葛向阳，张聪芝贾亚伟基于高通量测序技术解析中高温制曲真菌群落的演替规律研究大于0.99，说明该结果能够真实反映真菌群落结构。刚入房时，真菌丰度与菌群多样性都较低，说明自然接种的微生物种类与数量都较低；主发酵期，曲心温度升高，真菌大多不耐热，导致并房后真菌的丰度与菌群多样性都降至最低；此后进入潮火期，随着曲心温度进一步升高，能够适应高温的真菌种类在大曲中不断富集，因此上架时的真菌丰度与菌群多样性都开始升高；随着曲心温度保持60左右进入占火，嗜热真菌进一步增殖，曲皮水分也丧失殆尽，大火期的真菌丰度与菌群多样性进一步升至最高值；经过大火期的高温驯化与淘汰，下架时的真菌丰度与大火期相当，但菌群多样性降低；此后进入养曲期，曲心水分进一步降低，干燥缺水的环境逐渐不适合真菌的生长，因此出房时的真菌丰度基本不变，菌群多样性也进一步降低。表1大曲培养各阶段真菌Alpha多样性指数样本序列数OTU数ChaolSimpsonAceCoverageRF29674749.00.2962.510.995BF577162023.10.8227.50.999SJ163855757.20.2464.320.999DH456227777.30.2082.020.999XJ857827676.20.3580.600.999CF750657878.30.4681.680.9992.3OTU分布Venn分析通过对中高温制曲各阶段的OTU进行统计分类，Venn分析结果如图2 所示。结果显示，真菌有11个OTUs是共有的；从入房到并房，真菌无变化；上架时缓慢增加；进入大火期，真菌快速增加至最高点；下架时，真菌缓慢下降；出房时，真菌有所回升。2.4不同分类水平的物种变化通过合并相同物种分类的OTUs，从门、纲、目、科、属不同分类水平，统计平均含量前1 5的物种，其余的物种归于others中，根据百分比绘制真菌的分布柱状图，结果见图3与图4。从图3中可以看出，不同培养阶段大曲中的真菌都是以子囊菌门和毛霉门为主，两者占比超过90%以上，其中子囊菌门包括半子囊菌纲的酵母目和散囊菌纲的散囊菌目，毛霉门即毛霉纲的毛霉目。毛霉目前期很少，上架后才大量繁殖，丰度占SJXJ4611BFRF1310DHCF图2大曲培养各阶段真菌OTUs交叠Venn图比高达42%，大火期间降至1 1%，下架后回升至32%，出房后又缓落至2 0%酵母目的丰度则先缓升后不断下降，入房时即高达8 8%，并房后升至99%以上，之后开始快速下降，上架、大火、下架时分别下降至52%、31%6%，出房时回升至8%；散囊菌目则先下降后上升，入房时占比1 0%，而后缓慢下降至3%，大火后迅速生长繁殖，丰度占比高达57%，而后维持较高水平，下架、出房时分别为53%、7 0%。从科水平看，以相对丰度1%为下限，共检测到1 2 个真菌科。入房时以酵母科、发菌科、曲霉科和法夫酵母科为主，占比分别为2 5%、2 2%、1 8%、16%，此外还有丰度不足5%的双足囊菌科、德巴列酵母科、根霉科、横梗霉。酵母科在并房后消失殆尽；发菌科先快速减少后逐渐增加，并房时消失殆尽，上架后开始逐渐增加，各个阶段的丰度分别为2%、2 6%、54%、7 0%；曲霉科先减少、后增加、而后又减少，仅大火期高达32%，其他阶段丰度都在1%2%；法夫酵母科呈减少趋势，并房时占比6%，之后降至1%以下；双足囊菌科从入房时的4%增加至并房时的1 1%，而后降至1%以下，出房时回升至5%；德巴列酵母科在大曲培养各阶段呈先增加后减少趋势，上架和大火期丰度较高，分别为32%、2 8%；根霉科则呈波浪式变化，出房时丰度达1 0%；属于毛霉目的横梗霉也大致呈先增加后减少趋势，上架时丰度最高达56%，出房时为11%。在属分类水平上，以相对丰度1%为下限，共362023 No.10(Tol.352)酿酒科技2023年第1 0 期（总第352 期）LIQUOR-MAKINGSCIENCE&TECHNOLOGY1.00-1.000.750.75classC-DothideomycetesC-EurotiomycetesPhylumC-LeotiomycetesP-AscomycotaC-Mucoromycetes0.50P-Basidiomycota0.50C-SaccharomycetesP-MucoromycotaC-SordariomycetesP-unidentifiedC-TremellomycetesC-unidentified0.250.25-C-Wallemiomycetes0.000.00-RFBFSJDHXJCFRFBFSJDHXJCF1.001.00-familyOrderF-AspergllaceaeF-Debaryomycetaceae0.75O-Eurotiales0.75O-HypocrealesF-DipodascaceaeO-MicroascalesF-LichtheimiaceaeO-MucoralesF-MetschnikowiaceaeO-O-TremellalesF-Others0.500.50O-PleosporalesF-PhaffomycetaceaeO-SaccharomycetalesF-RhizopodaceaeO-SordarialesF-SaccharomycetaceaeO-TrichosporonalesF-Saccharomycetales_fam_lncertae_sedis0.25O-unidentified0.25F-TrichocomaceaeO-WallemialesF-unidentifiedF-Wallemiaceae0.000.00SJXJ1RFBFDHCFRFBFSJ DHXJCF图3大曲培养各阶段真菌在门、纲、目与科水平上的分布变化genus1.00G-AspergillusG-CandidaG-HyphoPichiaG-Kodamaea0.75G-LichtheimiaG-MillerozymaG-Monascus0.50G-othersG-PenicilliumG-RasamsoniaG-RhizoG-MucorG-Rhizopus0.25G-ThermomycesG-TorulasporaG-unidentified0.00G-WallemiaG-WickerhamomycesRFBFSJDHXJCF图4大曲培养各阶段真菌在属水平上的分布变化检测到1 6 个真菌属，结果如图4所示。入房时以孢圆酵母属、嗜热真菌属、威克汉姆酵母属和曲霉属为主，占比分别为2 6%、2 1%、1 7%、1 4%；并房后，微生物群落发生显著变化，主体微生物调整为柯达酵母属、根毛霉属、横梗霉属、威克汉姆酵母属、生丝毕赤酵母属、假丝酵母属，丰度分别为35%、19%、1 0%、6%、5%、5%；上架后，根毛霉属、生丝毕赤酵母属、横梗霉属的丰度分别增长至45%、31%、1 1%；大火期，主体微生物转变为30%的曲霉属、2 9%的生丝毕赤酵母属、2 7%的嗜热真菌属和7%的根霉属；下架后，嗜热真菌属、根毛霉属大量增长至54%、31%；出房后，嗜热真菌属占据绝对优势，丰度高达7 3%，此外还包括1 1%的根毛霉属和1 0%的根霉属。2.5群落分布特征通过PCoA主坐标分析与丰度聚类热图研究中高温制曲各阶段真菌群落分布特征，结果显示，37杨勇，葛向阳，张聪高通量测序技术解析中高温制曲真菌群落的演替规律研究业伟图5中的第一轴能够解释数据6 4.6%的变量，第二轴能够解释数据1 8.5%的变量，两轴和为8 3.1%。结合样品在分布图中的距离及其丰度聚类热图，可将真菌变化趋势的类型分为四类：（1)在入房和并房阶段，距离接近、同处于第二象限且聚为一类，说明这两个阶段的真菌群落结构相似，可归为同一阶段，其中孢圆酵母属主要在入房阶段，柯达酵母属主要在并房阶段；（2）上架时，样本与入房和并房阶段距离较远，且单独位于第三象限，可单独归为一类，此时伴随温湿度等环境因子的变化，大曲中水分营养充足，微生物结构不断进行适应性调整，横梗霉属大量生长至最高丰度，假丝酵母属、粉状米勒氏酵母属等丰度也升高；（3)大火期，样本与入房、并房和上架阶段的距离都较远，且单独位于第四象限，结合丰度聚类热图，也单独归为一类，此时温度等环境因素达到极限，不耐热微生物逐渐被淘汰，生丝毕赤酵母属、曲霉属、节担菌属、青霉属大量生长至最高丰度；（4）下架虽然与大火期距离较近，但下架和出房样本同位于第一象限且距离不远，结合丰度聚类热图可将下架和出房样本聚为一类，两者可归为同一阶段，此时随着品温等环境因子渐趋温和，微生物利用曲心残留的水分和营养进行生长繁殖与演替，根毛霉属、粉状米勒氏酵母属、横梗霉属及未知属在下架时丰度占据绝对优势，柯达酵母属、假丝酵母属、威克汉姆酵母属、根霉属在出房时丰度占据绝对优势，嗜热真菌属、红曲霉属、罗萨氏菌丰度在下架和出房时都较高。2.6真菌群落与理化因子关联分析利用R语言中RDA分析和作图，分析中高温制曲过程中的真菌群落结构与理化指标之间的相关性。理化因子关联分析发现，水分、糖化力、酯化力与真菌群落结构的相关性较高，发酵力其次，酸度与真菌群落结构的相关性较低，说明真菌主要对水分和酶活有较大影响，对酸度的影响较小。就主要优势属看，发酵力与生丝毕赤酵母属相关性较高，与曲霉属、青霉属也有一定的相关性；酯化力与嗜热真菌属、红曲霉、威克汉姆酵母属相关性较高，与罗萨氏菌、柯达酵母属、根霉属也有一定的相关性；水分和糖化力则与孢圆酵母属的相关性较高。0.38OREBF+SJ0.2口ODH%s81lzvodAXJ+CF0-0.2-0.4-0.51-0.690.6-0.4-0.200.20.20.55PCoA164.6%图5大曲培养各阶段真菌主坐标分析图G-Kodamaea2G-CandidaG-WickerhamomycesG-RhizopusG-Thermomyces0G-MonascusG-RasamsoniaG-Lichtheimia-1G-MillerozymaG-RhizoG-Mucor2G-unidentifiedG-TorulasporaG-WallemiaG-HyphoG-PichiaG-AspergillusG-PenicilliumXJCFDHSJRFBF图6大曲培养各阶段真菌群落丰度聚类热图3总结与展望本研究采用高通量基因测序技术对中高温制曲培养各阶段真菌不同分类水平的组成与丰度进行了系统的研究，结果表明：在门水平上以子囊菌门和毛霉门为主，其中子囊菌门包括半子囊菌纲的酵母目和散囊菌纲的散囊菌目，毛霉门即毛霉纲的毛霉目；在科水平上，以酵母科、发菌科、曲霉科和法夫酵母科为主；属水平上，嗜热真菌属、根毛霉属、根霉属、曲霉属、生丝毕赤酵母属、横梗霉属、威克汉姆酵母属分别在制曲过程不同时期占据主导地位，前期主要是酵母属，为大曲的酿酒和生香提供菌源，后期主要是耐高温霉菌类微生物，为大曲提供糖化、液化以及生香相关的微生物。唐清兰等1 2)采用高通量测序技术跟踪浓香型白酒大曲制曲过程真菌群落结构变化，发现伊萨酵母属、嗜热子囊菌属、曲霉属分别在制曲过程不同时期占据主导地位。不难发现，两类大曲的真菌种类及其丰度存在明显差异，这是由于大曲采用开放式生产、自然接种，环境、原料、工艺等对大曲微生物区系均有直接影响。2023 No.10(Tol1.352)38酿酒科技2023年第1 0 期（总第352 期）LIQUOR-MAKINGSCIENCE&TECHNOLOGYOODHG9G8G2G14G7M4OCFG15G16G13M5G610BG12G11M2M34RFSJM1G3G4G1G5OXJ-1.01.0注：Ml：水分；M2：酸度；M3：糖化力；M4：发酵力；M5：酯化力；Gl:Rhizomucor;G2:Hyphopichia;G3:Lichtheimia;G4:unidentified;G5:Millerozyma;G6:Thermomyces;G7:Rhizopus;G8:Aspergillus;G9:Penicillium;G10:Wickerhamo-myces;Gll:Torulaspora;G12:Candida;G13:Monascus;G14:Wallemia;G15:Rasamsonia;G16:Kodamaea。图7理化因子及酶活指标与真菌优势属的RDA分析此外，本文通过RDA分析首次发现真菌与水分和酶活关系密切。其中发酵力与生丝毕赤酵母属相关性较高，生丝毕赤酵母属是白酒发酵中的重要功能微生物，可以代谢产生风味物质，如乙酸乙酯、4-羟基-2-丁酮等1 3-1 4。酯化力与嗜热真菌属、红曲霉、威克汉姆酵母属相关性较高，其中嗜热真菌属耐高温，最高生长温度为50 或以上，可产生嗜热酶，包括纤维素酶、蛋白酶、淀粉酶、脂肪酶等2 ，红曲霉能分泌大量酯酶，对碳酸十四乙酯和苯乙醇的形成有促进作用1 5；威克汉姆酵母属能够将原料中的前体物质转化成酯、酸、醇、醛等风味物质，对白酒风味形成有重要作用6 。水分与孢圆酵母属的相关性较高，孢圆酵母属是非酿酒酵母菌，但通过其独特的次生代谢物也可对酒的感官特性作出积极贡献1 5。通过研究，明确了中高温制曲过程中真菌群落结构的演替规律，并首次解析了真菌种类与大曲的发酵和生香之间的潜在关系，有助于中高温大曲真菌多样性形成机制的解析，从而提高对大曲微生物认识的深度和广度，为中高温大曲生产提供参考与借鉴。下一步，我们将结合微生物与非生物因素进行关联性分析，如温度、湿度、氧气、门窗管控等，以全面准确地研究中高温大曲微生态演变的内在驱动力。参考文献：1李李静心,王艳丽，何宏魁，等.基于高通量测序技术解析高温大曲和中高温大曲的真菌群落结构J.食品与发酵工业,2 0 1 8,44(1 2):52-59.2 葛向阳,徐岩，周新虎,等.应用现代生物技术实现大曲自动化生产的研究J.酿酒科技,2 0 1 4(3)：1-3.3 罗惠波，杨晓东,杨跃寰，等.浓香型大曲中可培养真菌的分离鉴定与系统发育学分析J.现代食品科技,2 0 1 3(9):2047-2052.4 李燕荣,杨勇,张龙云，等.中高温大曲发酵过程微生物消长规律探析J.酿酒科技,2 0 1 9(5)：8 9-93.5 杨勇，李燕荣，姜雷，等.中高温大曲曲块部位间生化指标的差异及变化规律J.食品与发酵工业,2 0 1 9,45(1 9)：73-78.6赵东，郑佳，彭志云，等.高通量测序技术解析五粮液窖泥原核微生物群落结构J.食品与发酵工业,2 0 1 7,43(9):1-8.7谭崇尧，徐军.高通量测序法对不同地域浓香型大曲微生物结构的分析.酿酒科技,2 0 1 8(9)：1 1 8-1 2 2.8 吴树坤，谢军,卫春会,等.四川不同地区浓香型大曲微生物群落结构比较J.食品科学,2 0 1 9,40(1 4)：1 44-1 52.9施思,彭智辅，乔宗伟，等.浓香型大曲贮藏过程中糖化力发酵力变化及真菌多样性分析J.食品与发酵工业2017,43(5):76-79.10LI P,LIN W,LIU X,et al.Environmental factoringaffecting microbiota dynamics during traditionalsolid-state fermentation of Chinese Daqu starter J.Frontier in microbiology,2016,7:1237.11 杨勇,葛向阳,张龙云,等.基于高通量测序技术解析中高温制曲细菌群落的演替规律.酿酒科技,2 0 2 2(8)：34-41.12唐清兰，周利，徐姿静.浓香型白酒大曲制曲过程真菌群落结构变化的研究J.酿酒科技,2 0 2 1(9):2 7-35.13LI X R,MA E B,YAN L Z,et al.Bacterial and fungaldiversity in the starter production process of Fen liquor,a traditional Chinese liquorJJ.J microbiol,2013,51(4):430-438.14LI H,HUANG J,LIU X P,et al.Characterization ofinterphase microbial community in Luzhou-flavoredliquor manufacturing pits of various ages by polyphasicdetection methodsJJ.J microbiol biotechnol,2016,27(1):130-140.15 赵驰，苏伟,母应春，等.不同原料对酱香大曲微生物群落结构及多样性的影响J.食品科学,2 0 2 2,43(1 0)：133-141.16简晓平，简章容，余丽辰，等.1 株Wickerhamomyces酵母菌对浓香型白酒发酵的影响J.酿酒科技,2 0 1 7(9)：33-38.