作为中国人常喝的白酒，浓香型最受人们喜欢，有什么魔力吗？
阅读文章前，辛苦您点下右上角的“关注”，方便您讨论分享，持续关注每日优质内容~
作为中国市场占有率最高的白酒，浓香型白酒因其独特的酿造工艺和绵柔醇厚的口感受到了广大消费者的喜爱。
浓香型白酒在酿造过程中使用大曲多为中温大曲或中高温大曲，其中，中高温大曲以软质小麦为原料，蕴含了丰富的微生物群和酶类，为白酒的发酵生香提供了动力。
王清龙采用MiSeq高通量测序技术对河南南部、东部、西部和中部地区浓香型白酒大曲微生物群落结构进行了差异研究，发现南部和西部大曲的细菌属以芽孢杆菌属（Bacillus）为主，东部和中部大曲的细菌属以青枯菌属（Ralstonia）和乳杆菌属（Lactobacillus）为主。
MA采用MiSeq高通量测序技术对山东潍坊、江苏宿迁和四川绵阳地区浓香型大曲微生物群落结构进行了差异研究，发现山东潍坊大曲中细菌属以魏斯氏菌属（Weissella）和未分类的欧文菌科（unclassifiedErwiniaceae）为主，江苏宿迁大曲中细菌属以糖多孢菌属（Saccharopolyspora）、Bacillus和Weissella为主，四川绵阳大曲细菌属以变形杆菌属（Proteus）、Lactobacillus和Weissella为主。
吴树坤采用MiSeq高通量测序技术对四川宜宾、泸州和遂宁地区浓香型大曲的微生物群落结构进行了差异研究，发现宜宾和泸州大曲中细菌属以Weisslla和高温放线菌属（Thermoactinomyces）为主，遂宁大曲中细菌属以Bacillus和Staphylococcus为主。
由此可见，地域环境的不同会对大曲所含微生物的类群造成一定的影响。
同处于黄淮流域超级产区的山东临沂和河南周口地区，孕育了兰陵和宋河等多数知名浓香型白酒。
其中临沂市位于山东省东南部，地近黄海，该地区以沂、沭河为中心，西、北、东三面群山环抱，光照充足，雨量充沛，属温带季风气候。
周口市位于豫东平原，河南省东南部，该地区拥有沙颍河、涡河、西肥河和洪汝河四大水系，降水分布不均，日照、气温变幅大，属暖温带半湿润季风型气候。
由此可见，两地气候存在明显差异，且两地相差500公里左右，其环境亦存在一定差异。
因此，对临沂和周口地区的中高温大曲细菌类群进行比较分析，不仅可以加深研究人员对2个地区中高温大曲中微生物资源进行全面了解，亦可以进一步论证地域环境的因素是否会对大曲细菌类群结构产生影响。
近年来，随着第二代高通量测序技术的迅速发展，众多学者不断地利用该项技术来解析酒醅、酒曲和窖泥的微生物多样性，同时该方法也被广泛应用于解析不同香型白酒酿造过程中酒醅微生物群落结构的动态变化趋势。
与传统的微生物培养方法相比，第二代测序技术具有无需培养、高通量、检测速度快和结果准确等优点，此外，该技术还能准确识别难以培养、丰度较低的微生物，全面和客观地揭示微生态系统中微生物群落结构，为探索白酒大曲中微生物的多样性和功能提供了新的工具。
本研究采用IlluminaMiSeq高通量测序技术对临沂和周口地区共28份中高温大曲细菌群落结构进行比较分析，同时基于直系同源蛋白数据库对2个地区中高温大曲细菌潜在的功能差异进行解析，以期为黄淮流域超级产区酿酒微生物资源的挖掘提供数据支撑和理论指导。
中高温大曲分别采集自山东省临沂市兰陵县（E117°41′~118°18′，N34°37′~35°06′）和河南省周口市鹿邑县（E115°02′~115°37′，N33°43′~34°05′），样品随机选取，及时粉碎后备用。
其中，临沂地区大曲样品采集13份，周口地区大曲样品采集15份。
中高温大曲制作步骤基本如下：选取麦粒完整、饱满和无农药污染的小麦为生产原料，加入深井水对其进行润料再辊磨粉碎，粉碎后二次加入深井水进行拌料，拌料均匀后输送到压曲机制作曲坯，待曲坯成型后立即入房排列进行卧曲，再发酵培养28-30d，曲坯培养至曲心水分基本全部蒸发，品温降至接近室温时，曲已培养成熟，即可将曲块出房入库，贮存3～6个月。
QIAGENDNeasymericonFoodKitDNA基因组提取试剂盒德国QIAGEN公司；5×TransStartTM FastPfu Buffer、dNTPsMix和FastPfu FlyDNAPolymerase北京全式金生物技术有限公司；引物338F/806R（正向引物序列338F：5’-ACTCCTACGGGAGGCAGCA-3’；反向引物序列806R：5’-GGACTACHVGGGTWTCTAAT-3’）上海桑尼生物科技有限公司；MiSeq高通量测序配套试剂美国Illumina公司。
VeritiFASTPCR仪美国ABI公司；DYY-12电泳仪北京六一仪器厂；UVPCDS8000凝胶成像分析系统美国ProteinSimple公司；ND-2000C微量紫外分光光度计美国NanoDrop公司；PE300MiSeq高通量测序平台美国Illumina公司；R930机架式服务器美国DELL公司。
按照试剂盒的操作步骤对中高温大曲中细菌的基因组DNA进行提取，并以其为模板，使用已添加核苷酸标签（barcode）的引物338F/806R对总DNA的16SrRNAV3~V4区域进行扩增，扩增体系与方法参照GUO的进行，之后将检测合格的PCR产物寄至上海美吉生物医药科技有限公司完成2×300bp的双末端测序。
本研究利用QIIME（v1.9.1）平台进行生物信息学分析。
首先遵循重叠区的碱基数≥10bp、重叠区错配率≤0.2、barcode碱基错配或有效序列长度≥50bp的原则对测序返回的双端序列进行拼接和高质量序列筛选，使用PyNAST进行校准并排齐序列，按照100%和97%相似度进行序列划分构建分类操作单元（operationaltaxonomicunits，OTU），使用ChimeraSlayer去除含有嵌合体的OTU，最后选取代表性序列在RDP、Greengenes和SILVA数据库中进行同源性比对并确定其分类学地。
计算各样品中细菌类群的超1、发现物种数、香农和辛普森等α多样性指数；基于加权和非加权UniFrac距离进行主坐标分析2个地区大曲细菌类群的β多样性；采用PICRUSt（phylogeneticinvestigationofcommunitiesbyreconstructionofunobservedstates）软件对大曲中细菌基因功能进行预测。
使用两独立样本的非参数（Mann-Whitney）检验和多元方差分析（multivariateanalysisofvariance，MANOVA）解析2个地区细菌类群的差异。
使用R软件（v4.1.3）绘制箱型图和PCoA散点图；使用在线网站（http://jvenn.toulouse.inra.fr/app/example.html）绘制Venn图；使用Origin2017软件绘制堆积柱状图；使用在线网站（http://huttenhower.sph.harvard.edu/galaxy/）绘制LEfSe图；使用STAMP软件绘制菌群功能差异图。
基于高通量测序技术，28份中高温大曲共获得1477414条高质量16srRNA基因序列，每份样品平均52764条，采用100%和97%相似度进行序列划分共得到9330个OTU。
超1指数、发现物种数和香农指数、辛普森指数分别反映了微生物物种的丰富度和多样性，指数越高代表物种的丰富度和多样性越高。
本研究首先对2个地区中高温大曲的α多样性进行了比较分析，结果如图1所示。
2个地区中高温大曲的平均超1指数分别为2310和2625，平均发现物种数分别为1353和1611，平均香农指数分别为5.35和5.68，平均辛普森指数分别为0.85和0.86，经Mann-Whitney检验分析发现，2个地区大曲菌群在4个α多样性指数上差异均不显著（P>0.05）。
经序列比对发现，2个地区中高温大曲共鉴定到40个门和904个属，其中，临沂地区大曲鉴定到22个门和522个属，周口地区大曲鉴定到39个门和797个属。
本研究将平均相对丰度>1.0%的细菌门和属定义为优势细菌门和属。
基于门和属水平2个地区中高温大曲细菌群落结构组成。
2个地区中高温大曲中优势细菌门有4个，分别为硬壁菌门（Firmicutes）、变形菌门（Proteobacteria）、放线菌门（Actinobacteria）和拟杆菌门（Bacteroidetes）。
在临沂和周口地区大曲中，Actinobacteria的平均相对丰度分别为1.05%和4.38%，Bacteroidetes分别为0.81%和2.86%，经Mann-Whitney检验分析发现两者差异显著（P<0.05）。
2个地区中高温大曲中平均相对含量大于1.0%的属有6个，分别为芽孢杆菌属（Bacillus）、高温放线菌属（Thermoactinomyces）、不动杆菌属（Acinetobacter）、克罗彭斯特德菌属（Kroppenstedtia）、地杆菌属（Geobacter）和泛菌属（Pantoea）。
在临沂和周口地区大曲中，Acinetobacter的平均相对丰度分别为8.08%和4.28%，Geobacter分别为0.002%和3.44%，Pantoea分别为0.43%和1.83%，经Mann-Whitney检验分析发现差异均显著（P<0.05）。
此外，魏斯氏菌属（Weissella）和黏液乳杆菌属（Limosilactobacillus）在周口地区大曲中的平均相对丰度仅为0.34%和0.04%，而其作为优势细菌属在临沂地区大曲中平均相对丰度分别为1.74%和1.55%。
有研究报道，在大曲发酵过程中Weissella和Limosilactobacillus等乳酸菌的代谢活动为发酵提供了酸性环境，有助于微生物群落的有效自我驯化，其中Weissella还参与了2，3-丁二醇的代谢，促进了大曲风味物质的积累，除此之外，乳酸菌是大曲中的优势细菌，在白酒发酵过程中具有促进美拉德反应和保持酿酒微生态环境等作用，是白酒风味物质演替的推动者之一。
短波单胞菌属（Brevundimonas）、葡萄球菌属（Staphylococcus）和金黄杆菌属（Chryseobacterium）虽为周口地区大曲中的优势细菌属，平均相对丰度分别为1.50%、1.44%和1.07%，但其在临沂地区大曲中的平均相对丰度仅为0.26%、0.18%和0.25%。
ZHANG发现Staphylococcus参与了2，3-丁二醇的代谢，且其代谢产生的水解长链脂肪酸和甘油的脂肪酶可生成游离脂肪酸和甘油，因而构成了白酒中香味物质本身或前体物质。
姚亚林发现Brevundimonas具有较强水解淀粉和蛋白质的能力，在中温和高温大曲中与水分和淀粉的生成呈现正相关。
亦有研究指出，Thermoactinomyces、Pantoea、Staphylococcus和Chryseobacterium在大曲开始发酵阶段只存在于空气中，之后随着发酵温度的降低才逐渐进入大曲中得到积累。
由此可见，2个地区大曲细菌类群在属水平上存在明显差异，且在后续研究中积极探索周围环境因素对大曲发酵过程中细菌群落的累积与代谢是极为必要的。
本研究进一步通过LefSe对2个地区中高温大曲中的生物标志物进行了甄别。
柱长显示了线性判别值，值越大其富集程度越高。
在LDA阈值得分为4.0时，2个地区中高温大曲中共有9个细菌类群存在显著差异（P<0.05），其中各有2个、3个、1个、2个和1个分别体现在门、纲、目、科和属水平上。
在属水平，Geobacter在周口地区大曲中富集，由此可见，Geobacter可作为周口地区大曲中的生物标志物；在属水平，没有细菌类群在临沂地区大曲中富集，因而临沂地区大曲中的生物标志物在属水平未能体现。
通过上述分析发现，2个地区中高温大曲细菌类群在门属构成上存在一定的差异，因而为进一步直观的展现其群落结构异同，本研究在OTU水平上绘制了Venn图并基于UniFrac距离进行了主坐标分析。
本研究共划分了9，330个细菌OTU，其中2，912个OTU为临沂特有，1，212个OTU为周口特有，5，206个OTU为2个地区共有。
由此可见，2个地区大曲既具有各自独特的细菌类群，同时还含有共有的细菌群系。
在考虑各OTU序列相对丰度的情况下进行主坐标分析时，第一主成分和第二主成分的总贡献值为56.42%，2个地区样品间虽存在明显的重叠现象，但亦具有一定的分离趋势，经多元方差分析其方差值为5.81，结果可信度大于99.99%；由在仅考虑OTU中是否含有序列的情况下进行主坐标分析时，临沂地区大曲位于X轴的正方向，周口地区大曲位于X轴负方向，2个地区样品间存在明显的分离趋势，经多元方差分析其方差值为4.83，结果可信度大于99.99%。
由此可见，2个地区中高温大曲中细菌类群在整体上存在明显差异。
在比较了2个地区中高温大曲细菌菌群结构的差异后，本研究对其菌群的基因功能亦进行了差异分析，将序列上传至PICRUSt软件进行比对后得到具有差异的功能类群。
相较于周口地区大曲，临沂地区大曲菌群基因功能在未知功能、转录功能、氨基酸转运与代谢功能、一般功能预测和核苷酸转运与代谢功能上的表达均显著偏高（P<0.05），在细胞内运输功能、翻译后修饰，蛋白质周转，伴侣功能、信号转导机制功能和能量生产和转换功能上的表达均显著偏低（P<0.05）。
焦晶凯的研究发现，乳酸菌的蛋白水解系统可为其生长提供所需氨基酸，结合大曲中细菌属的分析结果来看，Weissella和Limosilactobacillus在临沂地区大曲中平均相对丰度均>1.0%，而其在周口地区大曲中平均相对丰度均不足0.4%，这可能是造成临沂地区大曲中细菌氨基酸转运与代谢功能较强的原因。
由此可见，2个地区中高温大曲的细菌菌群基因功能亦存在明显差异。
本研究通过MiSeq高通量测序技术对黄淮流域超级产区临沂和周口地区中高温大曲细菌群落结构和菌群基因功能进行了比较分析发现，结果发现，相较于周口地区大曲，临沂地区大曲中Weissella和Limosilactobacillus平均相对丰度明显偏高，Brevundimonas、Staphylococcus和Chryseobacterium平均相对丰度明显偏低，且其在转录、氨基酸和核苷转运与代谢功能表达上明显偏高；Geobacter为周口地区大曲中的生物标志物，而临沂地区大曲中的生物标志物在属水平未能体现。
由此可见，临沂和周口地区中高温大曲的细菌类群结构和菌群基因功能均存在明显差异，并进一步证实了地域因素会对大曲细菌类群结构产生影响。