王俊鋼，李宇輝*，劉成江，郭安民，岳建平

1(新疆農(nóng)墾科學(xué)院 農(nóng)產(chǎn)品加工研究所，新疆 石河子，832000) 2(農(nóng)產(chǎn)品加工重點(diǎn)實(shí)驗(yàn)室(新疆農(nóng)墾科學(xué)院)，新疆 石河子，832000) 3(額敏縣新大同創(chuàng)生物工程有限責(zé)任公司，新疆 額敏，834699)

風(fēng)干牛肉是新疆哈薩克族傳統(tǒng)的發(fā)酵肉制品，也是山區(qū)牧民過冬必備食品[1]。新疆北疆牧區(qū)的牧民通常在11月到來年1月制作風(fēng)干肉，牛肉經(jīng)過分割切條、腌制后自然風(fēng)干，制作時(shí)間約為2～8個(gè)月不等，其產(chǎn)品類型和我國的火腿[2]、土耳其的pastirma[3]和意大利的bresaola[4]相類似。近年來，隨著新疆旅游業(yè)的發(fā)展，風(fēng)干肉作為新疆的特色旅游食品，其消費(fèi)量逐年上升。傳統(tǒng)哈薩克族風(fēng)干肉的發(fā)酵主要依靠天然多菌種混合發(fā)酵，包括乳酸菌、微球菌、酵母菌、霉菌等各種微生物[5]，發(fā)酵體系中微生物菌群結(jié)構(gòu)復(fù)雜，其中真菌所產(chǎn)生的酶可以對(duì)肉中的糖類、脂肪和蛋白質(zhì)進(jìn)行分解，生成酮醛酯等風(fēng)味物質(zhì)提高產(chǎn)品品質(zhì)、賦予產(chǎn)品特殊的風(fēng)味和色澤[6-9]。FLORESM等[10]研究發(fā)酵香腸發(fā)現(xiàn)德巴利酵母菌(Debaryomycesspp.)對(duì)成熟過程中揮發(fā)性物質(zhì)的產(chǎn)生有重要影響，通過抑制脂質(zhì)氧化產(chǎn)物的產(chǎn)生、促進(jìn)乙酯的生成，可以產(chǎn)生香腸特有的香氣。SIRINE等[11]研究表明青霉菌具有抗氧化、防止酸敗和保持色澤的作用，還可以防止香腸表面形成黏性物質(zhì)。

傳統(tǒng)食品中的微生物菌落結(jié)構(gòu)更加復(fù)雜[12-14]。為了更加清楚地了解傳統(tǒng)發(fā)酵食品中微生態(tài)體系，可借鑒生態(tài)學(xué)中對(duì)微生物多樣性研究技術(shù)和方法[15]。高通量測(cè)序技術(shù)是微生態(tài)學(xué)研究方法的一種，能夠快速全面的分析微生物結(jié)構(gòu)復(fù)雜的樣品[16]，QUIJADA等[17]使用高通量測(cè)序技術(shù)對(duì)西班牙發(fā)酵香腸“Chorizo de León”中微生物菌群多樣性研究發(fā)現(xiàn)，不同地區(qū)發(fā)酵香腸中微生物在亞屬分類等級(jí)上有明顯差異；REBECCHI等[18]采用Illumina高通量測(cè)序研究豬肉、牛肉和羊肉臘腸中微生物的多樣性，結(jié)果顯示不同動(dòng)物來源的香腸微生物多樣性存在顯著差異。

由于地理環(huán)境的差異，使不同來源的風(fēng)干肉中微生物多樣性有所不同[19]。但目前對(duì)新疆風(fēng)干肉中微生物菌群多樣性還沒有系統(tǒng)的研究，尤其是真菌，其對(duì)環(huán)境的要求更低，適應(yīng)性更強(qiáng)，但在風(fēng)干肉中的群落結(jié)構(gòu)還不清楚，因此很有必要對(duì)新疆風(fēng)干肉中的真菌多樣性進(jìn)行研究，明確真菌結(jié)構(gòu)組成，為風(fēng)干肉安全性研究和下一步挖掘?qū)Ｓ冒l(fā)酵劑提供理論基礎(chǔ)，同時(shí)也可為后期產(chǎn)品的工業(yè)化生產(chǎn)提供依據(jù)。

1 材料與方法

1.1 樣品

選取新疆北疆4個(gè)哈薩克族聚居的牧區(qū)為采樣對(duì)象，全部為自然發(fā)酵。制作時(shí)間都為當(dāng)年11月份制作的新鮮風(fēng)干肉，且自然風(fēng)干時(shí)間都為40 d，采集的樣品單獨(dú)裝入采樣袋，放在車載冰箱(0～4 ℃)條件下24 h內(nèi)運(yùn)回實(shí)驗(yàn)室，存貯在-80 ℃冰箱，備用。詳細(xì)分組情況如表1所示。其中采樣地方相同的用相同的大寫字母，數(shù)字代表來自于不同的家庭。

表1 樣品來源及分組Table 1 Source and grouping of air-dried meat samples

1.2 試劑

強(qiáng)力土壤?DNA提取試劑盒,美國OMEGA公司；2×TaqMaster Mix (Dye Plus),南京諾唯贊生物科技有限公司；Magnetic Soil and Stool DNA Kit試劑盒,天根生化科技(北京)有限公司；Axygen凝膠回收試劑盒,Axygen美國；瓊脂糖,西班牙Biowest公司。

1.3 儀器與設(shè)備

高速冷凍離心機(jī)(Eppendorf 5810R)，德國Eppendorf AG公司；精密電子天平(BSM-220.4)，上海卓精；Bio-Rad電泳儀、PCR儀(T100)、凝膠成像系統(tǒng)(HP1020)，美國Bio-Rad伯樂公司；立式壓力蒸汽滅菌器(LDZX-30KBS)，上海申安醫(yī)療器械廠。

1.4 方法

1.4.1 真菌DNA提取與檢測(cè)

直接使用DNA提取試劑盒提取樣品中真菌的總DNA。提取方案按照試劑盒說明書進(jìn)行，每個(gè)樣品提取3次，提取的DNA用乙醇沉淀并提純。DNA濃度用分光光度法檢測(cè)，之后在1.2%瓊脂糖凝膠電泳上檢測(cè)DNA的完整性。

1.4.2 真多樣性分析

選取真菌的ITS1 區(qū)序列進(jìn)行高通量測(cè)序分析。使用PCR儀進(jìn)行PCR擴(kuò)增，引物選擇ITS1 區(qū)通用引物ITS1F：5′-CTTGGTCATTTAGAGGAAGTAA-3′；ITS2：5′-GCTGCGTTCTTCATCGATGC-3′，擴(kuò)增出的片段用1.2%瓊脂糖凝膠電泳檢測(cè)，使用凝膠成像系統(tǒng)進(jìn)行觀察，效果較好的樣本于2%瓊脂糖凝膠電泳切膠回收，以回收產(chǎn)物為模板進(jìn)行1次8 循環(huán)的PCR 擴(kuò)增，將Illumina平臺(tái)測(cè)序所需要的接頭、測(cè)序引物、標(biāo)簽序列添加到目的片段兩端。全部PCR 產(chǎn)物采用Axy Prep DNA 凝膠回收試劑盒進(jìn)行回收，并用FTC-3000TMReal-Time PCR 儀進(jìn)行熒光定量，均一化混勻后完成文庫構(gòu)建，在IlluminaMiSeq 2×300 bp平臺(tái)上完成測(cè)序。

PCR反應(yīng)體系為：4 μL 5×FastPfu Buffer，2 μL 2.5 mmol/L dNTPs，0.8 μL 5 μmol/L Forward Primer，0.8 μL 5 μmol/L Reverse Primer，0.4 μL FastPfu Polymerase，0.2 μL BSA，10 ng Template DNA補(bǔ)dd H2O至20 μL。程序如下：95 ℃預(yù)變性5 min，95 ℃變性60 s，50 ℃退火60 s，72 ℃延伸60 s，循環(huán)數(shù)為25，72 ℃最終延伸7 min。擴(kuò)增結(jié)束后，采用2%瓊脂糖凝膠電泳檢測(cè)擴(kuò)增產(chǎn)物的純度。使用凝膠回收試劑盒切膠回收PCR產(chǎn)物進(jìn)行文庫構(gòu)建后，在IlluminaMiSeq PE 300平臺(tái)上進(jìn)行高通量測(cè)序分析。

1.5 數(shù)據(jù)處理及分析

使用QIIME 中的PyNAST[20]對(duì)代表性序列進(jìn)行比對(duì)并去除低質(zhì)量、接頭污染和低復(fù)雜的序列，用FLASH[21]對(duì)有效序列進(jìn)行拼接，根據(jù)序列首尾兩端的barcode 和引物序列區(qū)分樣品。再用UPARSE[22]軟件將序列聚類到97%的操作分類單位(operational taxonomic unit,OTU)[23]。篩選重復(fù)序列，使用Uchime軟件去除嵌合體[24]。所有數(shù)據(jù)分析使用R語言[25]和Python(https：//www.python.org/)進(jìn)行。結(jié)合熱圖或聚類分析計(jì)算物種豐度與環(huán)境因子的相關(guān)性及其顯著性[26]。

2 結(jié)果與分析

2.1 樣品ITS1基因序列質(zhì)量評(píng)估及Alpha多樣性分析

26個(gè)樣品測(cè)序共獲得 2 059 240對(duì)有效的基因序列條帶，平均每個(gè)樣品產(chǎn)生35 455條有效序列。為了得到更高質(zhì)量及更準(zhǔn)確的生物分析結(jié)果，對(duì)雙端 Reads拼接、過濾后共產(chǎn)生1 936 392條優(yōu)化序列，每個(gè)樣品平均產(chǎn)生74 477條優(yōu)化序列。

根據(jù)97%相似性水平下的OTU信息，采用Alpha多樣性指標(biāo)的Chao1、Ace、Shannon、Simpson 指數(shù)以及群落覆蓋度指數(shù)Coverage等[27]對(duì)樣品中物種豐度和多樣性進(jìn)行評(píng)估，結(jié)果如表2所示。所有樣品的OTU覆蓋率為99.98%～100.00%，說明樣品中真菌種類幾乎都能被檢測(cè)到，測(cè)序結(jié)果可以很好的反應(yīng)樣品中真菌的真實(shí)情況。OTU數(shù)量越多說明樣品中的真菌豐度越高[28]，由表2可知，樣品G1、G2、J1、J2、L1和L2的OTU數(shù)量在100以上，其他樣品中的OTU數(shù)量基本都在100以下。這說明G1、G2、J1、J2、L1和L2這6個(gè)樣品中真菌種類較多，其中G1、G2、J1、J2源自于富蘊(yùn)縣，L1和L2屬于巴里坤地區(qū)，這可能是由于當(dāng)?shù)丨h(huán)境所引起的。來源于富蘊(yùn)縣的4個(gè)樣品均采自水草豐富、環(huán)境相對(duì)較好的牧區(qū)，微生物更加容易生長(zhǎng)，因此真菌的種類相對(duì)較高。來源于巴里坤山區(qū)的2個(gè)風(fēng)干牛肉樣品采自牧民家，這可能是和牧民的制作方式有關(guān)，衛(wèi)生條件直接關(guān)系到產(chǎn)品中的微生物情況。圖1-a表示真菌群落OTU 數(shù)目的Venn
圖，該圖可用于統(tǒng)計(jì)多個(gè)樣本中共有和獨(dú)有的OTU 數(shù)目，4組共有的OTU 數(shù)目84個(gè)，26個(gè)樣品種共有的OUT數(shù)目?jī)H有13個(gè)(圖1-b)，其中塔城(Group1)、額敏(Group2)、富蘊(yùn)(Group3)和巴里坤(Group4)地區(qū)樣品中OTU數(shù)目分別為125、133、261和316， Group1和 Group2之間差異不顯著(P>0.05)，其他組之間差異顯著(P<0.05)。

表2 Alpha多樣性指數(shù)統(tǒng)計(jì)Table 2 Alpha diversity index statistics

稀釋曲線反映了持續(xù)抽樣下新OTU(新物種)出現(xiàn)的速率；在一定范圍內(nèi)，隨著測(cè)序條數(shù)的加大，若曲線表現(xiàn)為急劇上升則表示群落中有大量物種被發(fā)現(xiàn)；當(dāng)曲線趨于平緩，則表示此環(huán)境中的物種并不會(huì)隨測(cè)序數(shù)量的增加而顯著增多(圖2-a)。Shannon指數(shù)越大，則表示OTU種類越多，物種越豐富，表明樣品中已涵蓋絕大多數(shù)的微生物信息，當(dāng)曲線趨于平坦時(shí)，說明測(cè)序數(shù)據(jù)量足夠大，OTU種類不會(huì)再隨著測(cè)序量的增加而增加。由圖2-b 可知，對(duì)26種風(fēng)干肉制品中的真菌進(jìn)行的高通量測(cè)序分析表明，當(dāng)測(cè)序量<40 000時(shí)，OTU數(shù)呈明顯上升趨勢(shì)，表明此時(shí)仍未檢測(cè)到更多種類的菌種。當(dāng)測(cè)序數(shù)量逐漸增加時(shí)，樣品的OTU數(shù)量增加，但趨勢(shì)趨于平穩(wěn)。當(dāng)測(cè)序深度達(dá)到60 000時(shí)，所有樣品的OTU變化均趨于平穩(wěn)。隨著測(cè)序量的增加，真菌的多樣性幾乎不變，這表明現(xiàn)有的測(cè)序量已經(jīng)可以反映樣品中真菌的豐度信息。采用隨機(jī)抽樣的方法對(duì)序列進(jìn)行測(cè)序，以繪制出繪制的序列數(shù)和可代表的OTU數(shù)的曲線。即可知多數(shù)樣本在序列數(shù)接近60 000時(shí)，OTU已經(jīng)接近飽和，說明現(xiàn)有測(cè)序量已經(jīng)可以反映出樣品中真菌豐度信息。

圖1 樣品中真菌群落OTU的Venn圖Fig.1 Venn diagram of OTU of fungal community in samples

圖2 樣品稀釋性、Shannon指數(shù)曲線Fig.2 Sample dilution and Shannon index curve

2.2 不同來源風(fēng)干肉中真菌群落組成分析

根據(jù)分類學(xué)分析結(jié)果，可以得到不同樣品在不同分類水平上的群落組成情況。圖3-a表示的是在門水平下，樣品中真菌豐度在0.1%以上的真菌用群落組成表示，豐度<0.1%的真菌用Other表示。由圖3-a可知，樣品中的優(yōu)勢(shì)真菌門有子囊菌門(Ascomycota)73.13%、擔(dān)子菌門(Basidiomycota)24.77%和被孢霉菌門(Mortierellomycota)2.76%。這3個(gè)菌門真菌含量占到測(cè)序序列總數(shù)的99%。不同地區(qū)風(fēng)干肉中真菌主要優(yōu)勢(shì)菌群略有差異，塔城地區(qū)優(yōu)勢(shì)真菌門為Basidiomycota，占該組真菌菌群的50.65%，而其他3個(gè)地區(qū)(額敏、富蘊(yùn)、巴里坤)主要優(yōu)勢(shì)真菌門為Ascomycota，分別占該組真菌菌群的70.33%、85.65%和80.47%。高慶超等[29]研究發(fā)現(xiàn)，同屬于自然發(fā)酵的黑果枸杞酵素，通過高通量測(cè)序在門分類水平下，子囊菌門(Ascomycota)、擔(dān)子菌門(Basidiomycota)占整個(gè)真菌門水平的99%。陳新[30]在對(duì)新疆新疆天山雪嶺土壤中微生物種群的研究中也發(fā)現(xiàn)，土壤真菌中被孢霉菌門(Mortierellomycota)、擔(dān)子菌門(Basidiomycota)和子囊菌門(Ascomycota)的相對(duì)豐度約占總和的70%，屬于優(yōu)勢(shì)種群。以上結(jié)果可以看出，雖然研究對(duì)象不同，但真菌群落結(jié)構(gòu)和優(yōu)勢(shì)菌群的差異不大，這是因?yàn)樽匀话l(fā)酵產(chǎn)品中的微生物主要來自于當(dāng)?shù)氐耐寥?、空氣和水，在門的分類水平上來看，自然發(fā)酵食品的真菌類型和土壤當(dāng)中的真菌基本一致。

圖3 各樣品在門、屬水平下真菌群落結(jié)構(gòu)分布Fig.3 Relative abundance of the major fungus phylum and genus in samples

不同樣品中真菌在分類屬水平上的分布情況如圖3-b所示。在26份風(fēng)干肉樣品中共鑒定出32個(gè)屬，相對(duì)含量>1%的10個(gè)屬分別為假絲酵母屬(Candida) 30.09%、Mrakia16.26%、Unclassified15.82%、德巴利酵母屬(Debaryomyces)11.14%、球腔菌屬(Mycosphaere)6.37%、青霉屬(Penicillium)4.23%、曲霉屬(Aspergillus)2.39%、(Mrakiella)2.88%、Naganishia1.71%和維希尼克氏酵母(Vishniacozyma)1.55%。其中，塔城牛肉(A1、A2)相較于其他幾個(gè)樣品所含真菌屬比較單一，優(yōu)勢(shì)真菌屬為Mrakia，占到85%左右。額敏牛肉(F1、F2)和托里牛肉(K1、K2)在真菌屬組成較為相似。不同地區(qū)風(fēng)干肉中真菌主要優(yōu)勢(shì)菌群不同，樣品塔城牛肉(C1、C2)中優(yōu)勢(shì)真菌屬為青霉屬(Penicillium),占40%左右;樣品富蘊(yùn)牛肉(G1、G2)中優(yōu)勢(shì)真菌屬為德巴利酵母屬(Debaryomyces),占60%左右;樣品巴里坤馬肉(M1、M2)中優(yōu)勢(shì)真菌屬為球腔菌屬(Mycosphaere),占65%左右。不同地區(qū)真菌在屬水平上有著一定差異(圖4-b)，塔城地區(qū)和額敏地區(qū)菌群結(jié)構(gòu)及種類比較類似，這可能是由于塔城和額敏都屬于塔額盆地，地理環(huán)境和氣候都比較類似。其他2個(gè)地區(qū)樣品中真菌多樣性則更加豐富，可能是由于這2個(gè)地區(qū)都屬于山區(qū)牧區(qū)，降水量多，物種相對(duì)豐富，因此，菌群多樣性也更加豐富。

a-門水平;b-屬水平圖4 四組樣品在門、屬水平下真菌群落結(jié)構(gòu)分布Fig.4 Relative abundance of the major fungus phylum and genus in four groups of samples

有研究表明，漢遜德巴利酵母(Dabaryomyces)和假絲酵母(Candida)是發(fā)酵肉制品的生產(chǎn)中最常用的酵母菌[31]，因?yàn)槠淇纱龠M(jìn)發(fā)酵香腸的呈色，改善發(fā)酵香腸的風(fēng)味，另外，有些酵母菌還能抑制腐敗菌的滋生[32]，趙永強(qiáng)等[33]用米曲霉和魯氏酵母對(duì)合浦珠母貝肉進(jìn)行發(fā)酵發(fā)現(xiàn)，其可以有效改善產(chǎn)品風(fēng)味及游離氨基酸含量。另外，文開勇等[34]對(duì)四川傳統(tǒng)臘肉研究發(fā)現(xiàn)，曲霉屬(Aspergillus)為四川臘肉中優(yōu)勢(shì)真菌屬。發(fā)酵肉在長(zhǎng)時(shí)間的發(fā)酵過程中，真菌的生長(zhǎng)繁殖可以產(chǎn)生多種酶類化合物，將肉中大分子的蛋白質(zhì)、脂肪分解成小分子的醇酮醛等化合物，并能增加氨基酸的種類和游離脂肪酸的形成，促使亞油酸轉(zhuǎn)化成共軛亞油酸，有利于風(fēng)干肉中特殊風(fēng)味的形成，顯著提升發(fā)酵肉的風(fēng)味品質(zhì)和營養(yǎng)價(jià)值[35]。另外，酵母菌和霉菌還可以促進(jìn)發(fā)酵肉中乳酸菌的生長(zhǎng)，加速肉制品中水分含量的降低，并且能夠降低風(fēng)干肉的酸度，抑制發(fā)酵肉中脂質(zhì)氧化、過氧化物的生成及腐敗微生物生長(zhǎng)繁殖，從而提高風(fēng)干肉的安全性能[36-37]。本實(shí)驗(yàn)中優(yōu)勢(shì)的真菌屬有假絲酵母屬(Candida,30.09%)、德巴利酵母屬(Debaryomyces,11.14%)和曲霉屬(Aspergillus,2.39%)，說明這3個(gè)屬的真菌在風(fēng)干肉的發(fā)酵過程中屬于優(yōu)勢(shì)菌群，對(duì)風(fēng)干肉品質(zhì)形成有一定的影響。

為了能夠更加直觀的看出不同來源不同樣品之間真菌屬水平的相似程度，通過物種豐度聚類熱圖(Heatmap)表示(圖5)。Heatmap是一顏色梯度來代表數(shù)據(jù)矩陣中數(shù)值的大小，并根據(jù)物種豐度相似性進(jìn)行聚類，將高豐度和低豐度的物種分類聚集，通過顏色梯度及相似程度反映多個(gè)樣品群落組成的相似性和差異性。由圖5可知，在真菌屬水平，Group1、Group2、Group3和Group4中樣品的聚類分析結(jié)果相互交叉，并有差異，其中L2、J2和D1樣品和其他樣品的群落組成差異較大，剩下的其他樣品之間菌群聚類結(jié)果比較相似。

圖5 樣品在屬水平物種豐度聚類熱圖Fig.5 Cluster heatmap of species abundance at genus level

26個(gè)樣品都屬于哈薩克族傳統(tǒng)手工制作，其中4個(gè)是馬肉，22個(gè)是牛肉，全部沒有添加發(fā)酵劑，都屬于自然發(fā)酵，但由于制作地區(qū)環(huán)境、風(fēng)干溫度以及制作人員的操作不同等原因?qū)е铝瞬煌瑯悠返恼婢郝浣M成差異。實(shí)驗(yàn)結(jié)果表明，風(fēng)干馬肉和風(fēng)干牛肉菌群多樣性差異不大。從PCoA圖可以看出(圖6)，26 個(gè)樣品的真菌群落組成大致可以分為2 個(gè)類型，其中Group1和Group2樣品分布相對(duì)較集中，說明這2組樣品真菌群落結(jié)構(gòu)比較相近，這與圖5的結(jié)果相互印證。Group3和Group4組中樣品分布在2個(gè)不同的位置，且相互之間距離較遠(yuǎn)，其中Group3組中J1、J2和G1、G2分別來自于富蘊(yùn)的牧區(qū)，而其他樣品則是山區(qū)，地理環(huán)境造成了其真菌菌群的差異。Group4組中K1、K2和其他樣品真菌群落差異較大則是因?yàn)長(zhǎng)1、L2和M1、M2采自于臨近村落的2種不同的原料肉，而K1、K2則是采自于另外較遠(yuǎn)的村莊。

3 結(jié)論

本研究采用IlluminaMiSeq第二代測(cè)序技術(shù)研究了新疆北疆4個(gè)地區(qū)26個(gè)傳統(tǒng)風(fēng)干肉中真菌多樣性，并結(jié)合高通量測(cè)序結(jié)果揭示了風(fēng)干肉中真菌的多樣性，不僅可以檢測(cè)到可培養(yǎng)的真菌，而且還檢測(cè)出了大量的不可培養(yǎng)真菌，尤其是對(duì)自然發(fā)酵肉制品，可以更加全面的反應(yīng)樣品中的菌群多樣性。

高通量結(jié)果表明，4組26 個(gè)樣品測(cè)序共獲得真菌原始序列 2 059 240條，過濾后得到的有效序列為 1 936 392 條。26份風(fēng)干肉樣品中真菌在真菌分類門水平上共鑒定出5個(gè)門，平均豐度值>1%的3個(gè)真菌門分別為：子囊菌門(Ascomycota，73.13%)、擔(dān)子菌門(Basidiomycota，24.77%)和被孢霉門(Mortierellomycota，2.76%)。這3個(gè)菌門真菌含量占到測(cè)序序列總數(shù)的99%。不同地區(qū)風(fēng)干肉中真菌主要優(yōu)勢(shì)菌群不同，塔城地區(qū)優(yōu)勢(shì)真菌門為Basidiomycota，而其他3個(gè)地區(qū)(額敏、富蘊(yùn)和巴里坤)主要優(yōu)勢(shì)真菌門為Ascomycota。另外，在26份風(fēng)干肉樣品中共鑒定出32個(gè)屬，相對(duì)含量>1%的屬分別為假絲酵母屬(Candida，30.09%)、Mrakia16.26%、德巴利酵母屬(Debaryomyces，11.14%)、球腔菌屬(Mycosphaere，6.37%)、青霉屬(Penicillium，4.23%)、曲霉屬(Aspergillus，2.39%)、(Mrakiella) 2.88%、Naganishia1.71%和維希尼克氏酵母(Vishniacozyma,1.55%)。豐度聚類熱圖和Beta多樣性分析表明，L2、J2和D1樣品和其他樣品的群落組成差異較大。