沈 毅，程 偉，鄧小波，陳 波，王 西，甘浪飛，張亞?wèn)|

（四川省古藺郎酒廠有限公司，四川古藺646523）

醬香型白酒是我國(guó)的一種主要香型白酒，具有醬香突出、酒體醇厚、回味悠長(zhǎng)和空杯留香等特點(diǎn)[1]。郎酒作為經(jīng)典醬香型白酒的代表之一，其生產(chǎn)依賴傳統(tǒng)的高溫制曲、原料清蒸、入窖前堆積、多次發(fā)酵取酒等工藝。由于傳統(tǒng)的開(kāi)放式生產(chǎn)過(guò)程及固態(tài)發(fā)酵中的不可控因素，來(lái)自高溫大曲、酒醅和窖泥的多種復(fù)雜菌種參與混合發(fā)酵，郎酒微生物菌群組成一直不夠明確，其菌群多樣性有待闡明[2-3]。早期研究集中于白酒發(fā)酵過(guò)程中高含量的細(xì)菌菌群，對(duì)其種類和多樣性進(jìn)行研究，發(fā)現(xiàn)細(xì)菌菌群對(duì)白酒品質(zhì)和風(fēng)味物質(zhì)形成具有重要影響[4-5]。近幾年研究發(fā)現(xiàn)，真菌菌群作為釀造微生物群落的重要組成部分亦具有非常重要的作用（產(chǎn)酒精、產(chǎn)酶和產(chǎn)香），如真菌菌群在制曲前期和高溫堆積發(fā)酵中具有高于細(xì)菌菌群的數(shù)量，且釀酒酵母是酒精生成的主要功能菌，曲霉菌和根霉菌具有很強(qiáng)的酯化酶和糖化酶活力，產(chǎn)酯酵母可形成以乙酸乙酯為主體的復(fù)合香氣[6-10]。本研究致力于探究醬香型郎酒釀造過(guò)程中真菌菌群的組成特征及其多樣性分析，為明確郎酒酒香的來(lái)源，保持酒品的穩(wěn)定性奠定基礎(chǔ)。

目前，圍繞醬香型白酒發(fā)酵微生物的研究，多采用分離培養(yǎng)的方法通過(guò)表型和生理生化特征對(duì)微生物種類進(jìn)行鑒定，隨后通過(guò)變性梯度凝膠電泳(PCR-DGGE)技術(shù)、rRNA基因克隆文庫(kù)技術(shù)等解析微生物群落結(jié)構(gòu)。鑒于一些微生物菌群不可培養(yǎng)的特征以及解析深度的匱乏，釀酒微生物的豐富組成和變化仍有待深入研究[11-12]。近年來(lái)，分子生物學(xué)技術(shù)廣泛應(yīng)用于微生物多樣性的研究，尤其是高通量測(cè)序技術(shù)。高通量測(cè)序技術(shù)不僅具有數(shù)據(jù)量高、覆蓋面廣、準(zhǔn)確率高的特點(diǎn)，還可以全面地定性和定量揭示樣品微生物菌群的多樣性和組成信息[13-14]。該技術(shù)非常適合用于白酒釀造過(guò)程中復(fù)雜微生物體系的研究。因此，本研究應(yīng)用高通量測(cè)序技術(shù)，以醬香型高溫大曲、酒醅和窖泥為研究對(duì)象，解析其真菌菌落結(jié)構(gòu)和多樣性，以期能夠全面準(zhǔn)確地分析醬香型郎酒在高溫制曲、高溫堆積和多次發(fā)酵取酒的生產(chǎn)過(guò)程中不同材料中的真菌菌落結(jié)構(gòu)和差異，為醬香型郎酒高溫制曲、高溫堆積發(fā)酵工藝的改進(jìn)和完善提供理論基礎(chǔ)。

1 材料與方法

1.1 材料、試劑及儀器

曲樣：本實(shí)驗(yàn)所有樣品均采集自郎酒公司，其中高溫大曲取自貯存1個(gè)月的新制成品曲，酒醅取自第4輪次堆積的酒醅中間部分，窖泥為窖底泥，上述樣品均裝入無(wú)菌采樣袋后置于低溫保存。

耗材及試劑：高質(zhì)量DNA提取試劑盒Power-Soil DNA Isolation Kit購(gòu)自美國(guó)MO Bio公司，微生物高通量測(cè)序建庫(kù)試劑盒購(gòu)自北京百邁克公司，Phusion HF MM高保真PCR酶購(gòu)自美國(guó)BioLabs公司，MinElute?PCR Purification Kit購(gòu)自德國(guó)Qiagen公司，Sample Preparation Kit購(gòu)自美國(guó)Illumina公司，Taq DNA Polymerase購(gòu)自北京全式金公司，電泳相關(guān)試劑購(gòu)自北京全式金公司，引物（ITS1-F1：5'-CTTGGTCATTTAGAGGAAGTAA-3'/ITS1-R1：5'-GCTGCGTTCTTCATCGATGC-3'）合成于北京華大公司等。

儀器設(shè)備：Illumina HiSeq 2500測(cè)序平臺(tái)，美國(guó)Bio-Rad公司；S1000 PCR儀，美國(guó)Bio-Rad公司；CFX96熒光定量PCR儀；美國(guó)伯樂(lè)BIO-RAD Gel Doc?XR+凝膠成像系統(tǒng)；北京六一儀器廠DYY-8C水平電泳槽；德國(guó)IKA公司VORTEX-3旋渦混勻器；德國(guó)Eppendorf公司Research Plus移液器；德國(guó)Beckman公司Allegra X-22R低溫冷凍離心機(jī)；日本三洋公司超低溫冰箱，準(zhǔn)微量分析天平BT25S；德國(guó)Eppendorf BioPhotometer Plus核酸蛋白測(cè)定儀；美國(guó)Thermo公司NanoDrop 2000超微量分光光度計(jì)；優(yōu)普UPT純水機(jī)等。

1.2 試驗(yàn)方法

1.2.1 高溫大曲、酒醅和窖泥的菌群基因組DNA提取

本實(shí)驗(yàn)基因組DNA提取按照PowerSoil DNA Isolation Kit說(shuō)明書(shū)進(jìn)行，其中裂解步驟如下：首先將樣品置于滅菌烘干的研缽中，加入液氮研磨制備大小均一的粉末，準(zhǔn)確稱取0.25 g各樣品至Power-Bead Tubes中，加入裂解液Solution C1后，在最大轉(zhuǎn)速3200 r/min下連續(xù)渦旋振蕩15 min，經(jīng)碳化硅研磨珠的機(jī)械力和裂解液的化學(xué)反應(yīng)共同作用下，充分裂解微生物，之后按說(shuō)明書(shū)對(duì)微生物宏基因組DNA進(jìn)行提取。

1.2.2 高溫大曲、酒醅和窖泥的菌群DNA文庫(kù)建立與測(cè)序分析

高溫大曲、酒醅和窖泥的菌群基因組DNA用0.8%的瓊脂糖電泳檢測(cè)，利用超微量分光光度計(jì)測(cè)定其OD260/280值及其含量。同時(shí)我們?cè)O(shè)計(jì)了特異性引物（ITS1-F1/ITS1-R1）進(jìn)行PCR擴(kuò)增，PCR擴(kuò)增程序?yàn)椋?5℃預(yù)熱5 min；95℃變性30 s，50℃退火30 s，72 ℃延伸40 s，共計(jì)30 cycles；72 ℃延伸10 min，PCR產(chǎn)物經(jīng)1.8%瓊脂糖凝膠電泳檢測(cè)擴(kuò)增目標(biāo)產(chǎn)物的大小和條帶信息。

各樣品菌群DNA經(jīng)質(zhì)量檢測(cè)合格后，通過(guò)針對(duì)ITS1區(qū)域的特異性引物分別進(jìn)行PCR擴(kuò)增。PCR擴(kuò)增使用Phusion HF MM高保真PCR酶，PCR條件如下：98℃預(yù)熱1 min；98℃變性10 s，50℃退火30 s，72 ℃延伸60 s，共計(jì)30 cycles；最后72 ℃延伸5 min。PCR結(jié)束后，利用Qiagen公司的Min-Elute?PCR Purification Kit，按照試劑盒說(shuō)明書(shū)進(jìn)行PCR產(chǎn)物純化。針對(duì)純化產(chǎn)物，我們利用Illumina公司的Tru-Seq?DNA PCR-Free Sample Preparation Kit制備測(cè)序文庫(kù)并添加接頭序列，文庫(kù)經(jīng)過(guò)Thermo NanoDrop 2000測(cè)定質(zhì)量后，利用Illumina HiSeq 2500進(jìn)行高通量測(cè)序。

1.2.3 生物信息學(xué)分析

利用FLASH[15]對(duì)測(cè)序獲得的原始數(shù)據(jù)進(jìn)行拼接，之后使用Trimmomatic[16]軟件對(duì)拼接得到的序列進(jìn)行質(zhì)量過(guò)濾，并去除嵌合體（UCHIME[17]），獲得高質(zhì)量的Tags序列。利用UCLUST[18]軟件在相似性97%的水平上對(duì)序列進(jìn)行聚類，以測(cè)序所有序列數(shù)的0.005%作為閾值過(guò)濾OTU[19]?；赨NITE[20]和RDP（Ribosomal Database Project）[21]對(duì) OTU 進(jìn)行分類學(xué)注釋分析，同時(shí)計(jì)算該OTU在各樣品中的相對(duì)含量、利用ClustalW2[22]軟件以鄰接法進(jìn)行系統(tǒng)進(jìn)化樹(shù)分析，利用Mothur[23]軟件對(duì)樣品進(jìn)行Alpha多樣性分析等。

2 結(jié)果與分析

2.1 菌群基因組DNA的提取及ITS擴(kuò)增

如表1所示，本實(shí)驗(yàn)成功從高溫大曲、酒醅和窖泥中制備高質(zhì)量的菌落基因組DNA模板，其OD260/280值在1.8～1.9之間，含量均高于50 ng/μL。由圖1A可見(jiàn)，各菌落基因組DNA在電泳圖中有非常清晰的單一條帶，表1和圖1結(jié)果說(shuō)明，成功提取獲得高質(zhì)量的DNA。如圖1B所示，從3個(gè)樣品中各取用20 ng基因組DNA作為模板，利用靶向ITS1序列的特異性引物，從高溫大曲、酒醅和窖泥的基因組DNA中擴(kuò)增得到大小約250 bp的條帶，符合預(yù)期片段大小，可以進(jìn)行后續(xù)建庫(kù)測(cè)序?qū)嶒?yàn)。

表1 菌群基因組DNA的OD260/280值及濃度

圖1 基因組DNA和ITS1基因片段擴(kuò)增的電泳結(jié)果圖

2.2 樣品測(cè)序數(shù)據(jù)處理結(jié)果

本研究收集的郎酒高溫大曲、酒醅和窖泥樣品共產(chǎn)出598,550對(duì)reads，雙端reads拼接、過(guò)濾后共產(chǎn)生236,402條clean tags，平均每個(gè)樣品產(chǎn)出78,801條clean tags（見(jiàn)表2）。去除引物和barcode后的序列平均長(zhǎng)度分別為304 bp、258 bp和280 bp，Q30(%)為質(zhì)量值大于或等于30的堿基占總堿基數(shù)的百分比，3種樣品的Q30值均高于96%，說(shuō)明數(shù)據(jù)質(zhì)量較高。

表2 樣品測(cè)序數(shù)據(jù)處理結(jié)果統(tǒng)計(jì)

2.3 序列豐度和多樣性分析

如圖2所示，在97%的相似度水平下，針對(duì)236,402條代表性序列進(jìn)行聚類得到的各樣品OTU個(gè)數(shù)分別為83、70和308，合計(jì)OTU數(shù)量為350個(gè)。同時(shí)對(duì)3種樣品真菌菌群的構(gòu)成進(jìn)行比較分析，繪制得到的OTU-Venn圖如圖3所示，其中高溫大曲、酒醅和窖泥共有的OTU個(gè)數(shù)為13，高溫大曲和酒醅共有的OTU個(gè)數(shù)為25，高溫大曲和窖泥共有的OTU個(gè)數(shù)為63，酒醅和窖泥共有的個(gè)數(shù)為46，而高溫大曲、酒醅和窖泥特有的OTU個(gè)數(shù)分別為8個(gè)、22個(gè)和222個(gè)。

圖2 高溫大曲、酒醅和窖泥樣品中OTU個(gè)數(shù)分布圖

圖3 高溫大曲、酒醅和窖泥樣品的OTU-Venn圖

同時(shí)，經(jīng)Alpha多樣性分析得到高溫大曲、酒醅和窖泥微生物的豐度指數(shù)（Chao1）和多樣性指數(shù)（Shannon），其中三者的Chao1指數(shù)分別為123、91和 310，而 Shannon指數(shù)分別為 1.122、1.764和4.239。由此可見(jiàn)，郎酒窖泥中真核微生物具有最高的豐富度和多樣性，豐富度其次是高溫大曲和酒醅，但是酒醅中真核微生物的多樣性高于高溫大曲。

不同樣品用不同顏色表示，不同顏色圖形之間交疊部分?jǐn)?shù)字為兩個(gè)樣品之間共有的OTU個(gè)數(shù)，多個(gè)顏色圖形之間交疊部分?jǐn)?shù)字為多個(gè)樣品之間共有OTU個(gè)數(shù)，非交疊部分為各樣品特有OTU個(gè)數(shù)。

2.4 基于各分類學(xué)地位對(duì)高溫大曲、酒醅和窖泥的分析

在OTU劃分的基礎(chǔ)上，通過(guò)將OTU代表序列與微生物參考數(shù)據(jù)庫(kù)進(jìn)行比對(duì)分析，本研究對(duì)高溫大曲、酒醅和窖泥進(jìn)行分類學(xué)鑒定（表3），將序列鑒定為5個(gè)門，16個(gè)綱，40個(gè)目，77個(gè)科，109個(gè)屬和148個(gè)種，其中仍有部分序列經(jīng)過(guò)比對(duì)無(wú)法鑒定到已知的屬和種。

表3 郎酒高溫大曲、酒醅和窖泥樣品的各分類學(xué)地位數(shù)量統(tǒng)計(jì)

針對(duì)樣品中平均含量大于1%的真菌在屬和種水平上進(jìn)行相對(duì)含量比較分析，在350個(gè)OTU中，平均相對(duì)含量＞1%的OTU有20個(gè)，其中包含4個(gè)無(wú)法比對(duì)到已鑒定真菌菌屬的未知OTU。在已知的屬水平上，高溫大曲中的優(yōu)勢(shì)菌屬主要隸屬于嗜熱絲孢菌屬（Thermomyces）、曲霉屬（Aspergillus）和嗜熱子囊菌屬（Thermoascus），其平均相對(duì)含量分別為54.9%、24.2%和18.5%；酒醅中的優(yōu)勢(shì)菌屬主要為假絲酵母屬（Candida）、擬青霉屬（Paecilomyces）和嗜熱絲孢菌屬（Thermomyces），其相對(duì)含量分別為55.1%、16.5%和15.1%；窖泥中的主要菌屬為馬拉色菌屬（Malassezia）、被孢霉屬（Mortierella）、枝孢屬（Cladosporium）和曲霉屬（Aspergillus），其相對(duì)含量分別為17.7%、10.1%、6.5%和5.3%。

在已知的種水平（圖4），高溫大曲的主要菌種包括疏棉狀嗜熱絲孢菌Thermomyces lanuginosus（54.9%）、曲霉屬的Aspergillus cibarius（24.2%）和熱子囊菌屬的Thermoascus sp.（18.5%）；酒醅中的主要菌種包括假絲酵母屬的Candida xylopsoci（52.3%）、擬青霉屬的Paecilomyces formosus（16.5%）和疏棉狀嗜熱絲孢菌Thermomyces lanuginosus（15.1%）；而窖泥中的優(yōu)勢(shì)菌種包含限制性馬拉色菌Malassezia restricta（17.3%）、枝孢屬的Cladosporium halotolerans（4.4%）、被孢霉屬的Mortierella antarctica（3.9%）和雜色曲霉Aspergillus versicolor（3.5%）。

圖4 在種水平上高溫大曲、酒醅和窖泥中優(yōu)勢(shì)真菌相對(duì)含量的比較分析

3 結(jié)論

本研究采用HiSeq 2500高通量測(cè)序技術(shù)，對(duì)郎酒高溫大曲、酒醅和窖泥樣品中的真菌微生物進(jìn)行解析，研究顯示三者之間的真菌菌群組成和多樣性存在明顯差異。

高溫大曲作為醬香型白酒發(fā)酵中微生物的主要來(lái)源，掌握其中的主要風(fēng)味微生物的種類和數(shù)量對(duì)于提高白酒的品質(zhì)具有重要意義。鑒于醬香型郎酒高溫制曲的特殊工藝，本研究發(fā)現(xiàn),在高溫大曲真菌微生物中，嗜熱型真菌菌屬占優(yōu)勢(shì)地位，其中嗜熱絲孢菌屬（Thermomyces）和嗜熱子囊菌屬（Thermoascus）真菌的比例合計(jì)高達(dá)73.4%，上述結(jié)果與蔣英麗等[24]的最新研究報(bào)道相一致，通過(guò)追蹤不同階段的高溫大曲微生物組成，研究發(fā)現(xiàn)，高溫對(duì)醬香功能菌群的篩選和富集效應(yīng)明顯，在新制成品曲和儲(chǔ)存曲塊中，嗜熱菌屬具有非常高的相對(duì)含量。本研究中疏棉狀嗜熱絲孢菌Thermomyces lanuginosus不僅在高溫大曲中占據(jù)了高達(dá)55%的比例，也在酒醅和窖泥菌群中得到鑒定（15.1%和0.4%），而該菌種同樣廣泛存在于醬香型茅臺(tái)酒發(fā)酵過(guò)程中[25]，以及濃香型大曲[26]和瀘州老窖大曲中[27]。Thermomyces lanuginosus的生長(zhǎng)溫度據(jù)報(bào)道可以達(dá)到真菌生存的最高溫度60℃，該類群真菌可以分泌β-木聚糖酶，將植物細(xì)胞壁的主要組成成分木聚糖轉(zhuǎn)化成木糖等[28-29]。此外，在高溫大曲中還存在較高豐度的曲霉屬（Aspergillus）真菌微生物，該菌屬微生物廣泛應(yīng)用于醬油、日本清酒、中國(guó)黃酒等傳統(tǒng)釀造行業(yè)，具有一定的耐酸能力，可分泌糖化酶，具有較高糖化力，可提高出酒率[30-31]，而在郎酒醬味相關(guān)研究中亦推測(cè)高溫曲霉在醬香獨(dú)特風(fēng)味產(chǎn)生上具有重要的地位[2]。上述高溫大曲中的研究顯示，嗜熱菌屬和曲霉屬在高溫大曲中占據(jù)很高的豐度，對(duì)高溫大曲的獨(dú)特風(fēng)味物質(zhì)的組成和含量具有重要的影響。

酒醅在堆積過(guò)程中不僅利用來(lái)自高溫大曲中已存在嗜熱微生物菌屬，還會(huì)吸附來(lái)自于空氣和生產(chǎn)環(huán)境中的微生物，同時(shí)在堆積過(guò)程中的高溫可再次對(duì)優(yōu)勢(shì)菌群進(jìn)行富集，增加酒醅中風(fēng)味物質(zhì)的種類和含量，表明了酒醅中可形成獨(dú)特的產(chǎn)菌群成分。在酒醅的真菌微生物分析中，研究發(fā)現(xiàn)，優(yōu)勢(shì)菌種包括假絲酵母屬的Candida xylopsoci（52.3%）、擬青霉屬的Paecilomyces formosus（16.5%）和疏棉狀嗜熱絲孢菌Thermomyces lanuginosus（15.1%），其中，假絲酵母屬的Candida xylopsoci占據(jù)絕對(duì)主導(dǎo)地位，該菌種在釀造過(guò)程中作為主要的產(chǎn)酒功能菌發(fā)揮重要作用，同時(shí)，張春林等[32]在瀘州高溫大曲中篩選到酵母菌GY-3（即Candida xylopsoci）可作為產(chǎn)香的功能菌。Candida xylopsoci在新疆塔城牧區(qū)奶酪樣品中及乳清和凝乳中作為主要菌種存在[33]，同時(shí)作為優(yōu)勢(shì)菌種存在于寧紹地區(qū)的特色腌制蔬菜腌制冬瓜中[34]。因此，我們推測(cè)假絲酵母屬的Candida xylopsoci在郎酒堆積發(fā)酵過(guò)程可能參與香味物質(zhì)的生成等，對(duì)于醬酒特殊性風(fēng)味的形成具有重要的作用。對(duì)比醬香型郎酒高溫大曲和酒醅的分析結(jié)果，揭示二者之間的優(yōu)勢(shì)菌群存在一定的異同。其中，嗜熱絲孢菌屬（Thermomyces）在大曲和酒醅均可作為優(yōu)勢(shì)菌屬，而酒醅中的其他優(yōu)勢(shì)菌屬則由假絲酵母屬和擬青霉屬組成。我們推測(cè)，盡管酒醅中含有少量成品曲粉導(dǎo)致其與高溫大曲微生物菌群的初始組成非常相似，且堆積過(guò)程中酒醅核心溫度接近高溫制曲溫度，但是制曲過(guò)程采用生料，而酒醅發(fā)酵采用蒸煮后富含水分的熟料，即高溫大曲和酒醅中微生物處于兩種差異顯著的培養(yǎng)發(fā)酵環(huán)境，因此可導(dǎo)致二者微生物組成存在較大差異。

窖泥作為醬香型白酒發(fā)酵過(guò)程中的重要組成部分，其所在的窖池窖齡越長(zhǎng)，其富集的功能微生物越多，風(fēng)味物質(zhì)越豐富，而目前的應(yīng)用集中于窖泥微生物的復(fù)合應(yīng)用上，可有效提高原酒優(yōu)質(zhì)品率等[11]。窖池作為酒醅發(fā)酵的主要場(chǎng)所，提供獨(dú)特的厭氧環(huán)境，并通過(guò)窖泥提供豐富的釀造微生物，產(chǎn)生多種對(duì)酒體有貢獻(xiàn)的風(fēng)味成分[35]。本研究發(fā)現(xiàn)，在相同含量的基因組DNA模板中，窖泥樣品中通過(guò)建庫(kù)測(cè)序得到的有效reads數(shù)最少，ITS1片段的擴(kuò)增條帶含量最少，表明窖泥中真核微生物占總微生物含量的比例相對(duì)于高溫大曲和酒醅較低。有研究亦表明，在濃香型白酒發(fā)酵窖泥微生物中，原核微生物總類和豐度高于真核微生物，占據(jù)主導(dǎo)地位[36]。盡管窖泥真核菌群在總微生物中未占據(jù)優(yōu)勢(shì)地位，但本研究揭示窖泥中真核微生物的豐富性和多樣性遠(yuǎn)高于高溫大曲和酒醅，在所有鑒定到的305個(gè)OTU中，屬水平上相對(duì)含量超過(guò)1%的OTU有16個(gè)，呈現(xiàn)出豐富的真核微生物多樣性，與高溫大曲和酒醅中少數(shù)幾個(gè)優(yōu)勢(shì)菌屬占據(jù)超過(guò)近80%相對(duì)含量的現(xiàn)象存在顯著差異。同時(shí)，本研究發(fā)現(xiàn)，窖泥中的主要菌屬為馬拉色菌屬（Malassezia）、被孢霉屬（Mortierella）、枝孢屬（Cladosporium）和曲霉屬（Aspergillus），上述菌屬亦同樣存在于其他白酒真核微生物組成中，參與白酒釀造過(guò)程中復(fù)雜的物質(zhì)能量代謝和物質(zhì)形態(tài)變化過(guò)程[36-38]。

綜上所述，本研究采用高通量測(cè)序技術(shù)，對(duì)醬香型郎酒生產(chǎn)中的高溫大曲、酒醅和窖泥的真菌微生物進(jìn)行多樣性解析和比較分析，結(jié)果顯示，三者在真菌微生物多樣性和組成上存在較大差異，其中窖泥中微生物豐富度和多樣性最高，在窖泥中可以鑒定到大量共同存在于高溫大曲和酒醅中的真菌菌種；高溫大曲中主要由優(yōu)勢(shì)菌屬嗜熱絲孢菌屬（Thermomyces）和嗜熱子囊菌屬（Thermoascus）組成；酒醅中優(yōu)勢(shì)菌屬包括假絲酵母屬（Candida）、擬青霉屬（Paecilomyces）和嗜熱絲孢菌屬（Thermomyces），與高溫大曲和窖泥的微生物組成存在顯著差異。同時(shí)，本研究在郎酒高溫大曲、酒醅和窖泥中鑒定到部分無(wú)法比對(duì)到已知序列的菌屬和菌種，表明郎酒釀造過(guò)程中可能存在與風(fēng)味和品質(zhì)相關(guān)的獨(dú)特未知菌種，后續(xù)還將開(kāi)展更深入的研究進(jìn)行解析。