王香君,段 杉
(1.四川省農(nóng)業(yè)科學院蠶業(yè)研究所,四川南充 637000; 2.華南農(nóng)業(yè)大學食品學院,廣東廣州 510642)
魚露,又名魚醬油,是一種在高鹽條件下利用魚體自身酶系和多種微生物的協(xié)同作用,將海產(chǎn)低值魚等原料中蛋白質(zhì)、脂肪等成分進行數(shù)年發(fā)酵分解而成的營養(yǎng)豐富、風味獨特的液體調(diào)味品[1]。
目前,關于魚露中微生物的研究多采用平板計數(shù)、分離等傳統(tǒng)的純培養(yǎng)法,基于此方法在魚露微生物的研究中取得了一定的成果[2-3]。Kilinc等[4]和Harikedua等[5]利用平板計數(shù)法研究魚露發(fā)酵過程中細菌、真菌的數(shù)量變化情況,而一些文獻報道從發(fā)酵魚產(chǎn)品中篩選出一些耐鹽或嗜鹽的乳酸菌[3-4]。但由于自然環(huán)境樣品中僅有不到1%的微生物可培養(yǎng),因此用純培養(yǎng)法不能很好反映微生物群落存在的真實狀態(tài)。而基于高通量測序的16S rDNA分子生態(tài)分析技術,通過對細菌16S rDNA序列進行擴增、分析,可以準確分析樣品中細菌群落的組成,利于鑒定低豐度物種,現(xiàn)已成為研究細菌群落結構及多樣性的重要手段[6]。胡小喜等[7]利用此技術研究了不同發(fā)酵條件的蝦油細菌群落結構的組成情況,結果檢測出蝦油中細菌多達1285屬,全面的揭示蝦油中的細菌組成,相對于篩選和計數(shù)微生物的平板計數(shù)法具有明顯優(yōu)勢。
傳統(tǒng)魚露生產(chǎn)常采用露天發(fā)酵方式,工藝受氣候條件影響較大,工藝穩(wěn)定性差,衛(wèi)生條件不易控制,大量雜菌滋生,安全性低。胡小喜等[7]研究指出,35 ℃條件下,蝦油中蛋白酶失活較少,且存在大量產(chǎn)生良好風味的細菌。張豪等[8]、黃紫燕等[9]對(30±5) ℃發(fā)酵條件下的魚露進行研究。王香君等[10]用FQ-PCR法和平板計數(shù)法檢測到發(fā)酵魚露樣品的真菌總數(shù)遠低于細菌總數(shù)。Kilinc等[4]和Harikedua等[5]研究指出在魚露發(fā)酵過程中真菌數(shù)量很少,在魚露發(fā)酵過程中真菌作用很小。因此試驗采用16S rDNA高通量測序技術研究35 ℃發(fā)酵條件下魚露的細菌多樣性及群落結構,從而揭示魚露發(fā)酵過程中細菌群落的演替規(guī)律。
魚露原材料:5種新鮮魚(包括頭、內(nèi)臟等)混合而成,質(zhì)量百分比分別是雙斑瓦鰈(Poecilopsettaplinthus)2.3%、短尾大眼鯛(Priacanthusmacracanthus)5.2%、翼紅娘魚(Lepidotriglaalata)5.4%、鲬(Platycephalusindicus)6.3%、眶棘雙邊魚(Ambassisgymnocephalus)80.8% 陽江市海陵島,加冰塊方式及時運到實驗室快速處理;D6005真菌/細菌DNA提取試劑盒 ZYMO RESEARCH生物科技公司;GeneJET膠回收試劑盒 Thermo Scientific公司;NEB Next? UltraTMDNA Library Prep Kit for Illumina建庫試劑盒 New England Biolabs公司;食鹽 廣東省鹽業(yè)集團有限公司;硝酸銀、鉻酸鉀等均為分析純 國藥集團化學試劑有限公司。
DPX-9082B-2型電熱恒溫培養(yǎng)箱 上海?,攲嶒炘O備公司;LX-C50L立式壓力蒸汽滅菌器 合肥華泰醫(yī)療設備有限公司;JA2003電子天平 上海精密科學儀器有限公司;SW-CJ-1FDA 型超凈工作臺 蘇州安泰空氣技術有限公司;5415D高速冷凍離心機 德國Eppendorf公司;TL-MM-1型微量振蕩器 廣州云薈貿(mào)易有限公司;DYCP-31CN水平電泳儀 北京六一儀器廠;GEL-DOL 2000電泳凝膠成像分體系統(tǒng) BIO-RAD公司。
1.2.1 魚露樣品的制備 按照工廠生產(chǎn)魚露的工藝進行制備,將原料魚與食鹽按3∶1的質(zhì)量比直接混勻,然后再將魚露樣品放置于35 ℃的恒溫箱中恒溫發(fā)酵,發(fā)酵時間12個月,發(fā)酵過程中每隔5 d取樣測定含鹽量后,用無菌水補加蒸發(fā)的水分以保持鹽度恒定(25%,w/w)。并于發(fā)酵0、2、8、12個月時用滅菌的100 mL離心管取樣(取樣時,將魚露發(fā)酵液混勻),備用。
1.2.2 魚露樣品中細菌DNA的提取 將1.2.1取得的樣品在無菌條件下先用快速濾紙(滅菌處理)過濾掉粗顆粒后,用0.22 μm無菌濾膜過濾菌體,然后用無菌水將濾膜上的菌體洗下來。使用D6005真菌/細菌DNA提取試劑盒提取菌體總DNA。所有DNA樣品均于-20 ℃保存?zhèn)溆谩?/p>
1.2.3 細菌計數(shù) 根據(jù)Castillo等[11]報道,采用FQ-PCR法對魚露中的細菌總數(shù)進行計數(shù),并照王香君等[10]方法對樣品進行菌落計數(shù)。
1.2.4 16S rDNA 分子生態(tài)分析 以稀釋后的細菌總DNA為模板,用帶Barcode的16S rDNA V4區(qū)特異引物515F(5′-GTGCCAGCMGCCGCGGTAA-3′)、806R(5′-GGACTACHVGGGTWTCTAAT-3′)擴增16S rDNA 的V4區(qū),用GeneJET 膠回收試劑盒回收純化后的PCR產(chǎn)物,再使用NEB Next? UltraTMDNA Library Prep Kit for Illumina建庫試劑盒進行文庫構建,構建好的文庫經(jīng)過Qubit定量和文庫檢測,合格后,使用MiSeq進行上機由北京諾禾致源科技股份有限公司完成測序,測序得到原始數(shù)據(jù)(Raw Data)。
原始數(shù)據(jù)包含一定的接頭污染序列、低質(zhì)量數(shù)據(jù)等,會影響后續(xù)的分析,為了使信息分析的結果更加準確、可靠,需要對原始數(shù)據(jù)進行過濾。數(shù)據(jù)過濾得到clean data后,根據(jù)PE reads之間的重疊關系,將成對雙端reads拼接為一條序列,再以97%序列相似性為閾值,對clean data進行OTUs(Operational Taxonomic Units)聚類和物種分類分析,并將OTUs和物種注釋結合,從而得到每個樣品的OTUs和分類譜系的基本分析結果。通過對OTUs進行豐度、α-多樣性、β-多樣性分析,同時對物種注釋在各個分類水平上進行群落結構的統(tǒng)計分析,得到細菌群落組成及結構。
1.2.5 指標測定 含鹽量(以氯化鈉計):參照GB/T 5009.39-2003《醬油衛(wèi)生標準的分析方法》。
測定和分析結果采用Excel進行處理,結果采取均值±標準差形式。
圖1、表1顯示,發(fā)酵0、2、8、12個月的魚露中細菌分別測得78524、44031、29512、26744條序列,在科水平上,各樣品的序列占比在70%以上,在屬水平上,各樣品的序列占比在50%以下。其中分別有69869、33520、21531、19742條注釋到科,分別有34919、21182、12453、11777條注釋到屬。4個樣品測得的序列分別聚類了1197、1385、1286、1170個OTUs,圖2中有492個OTUs是所有發(fā)酵魚露樣品共同的,而8個月和12個月共同的OTU數(shù)量最多,達到818個,這說明發(fā)酵8個月和發(fā)酵12個月的魚露中細菌群落構成差異小。
表1 各分類水平上的序列數(shù)以及各自占比Table 1 The numbers of tags and their proportions at the different levels
圖1 不同樣品的序列和OTUs數(shù)量統(tǒng)計Fig.1 Statistics of the numbers of tags and OTUs in different samples
圖2 OTU交集韋恩圖Fig.2 Intersected Venn diagram analysis of OTU注:圖中數(shù)字代表樣品的OTUs數(shù)量。
熱圖展示了每個樣品中的物種組成及豐度信息,圖3顯示,不是所有的序列聚類到OTU后都能注釋到種,部分OTUs只能注釋到目、科、屬,有的只能注釋到界。結合表1可知,由于注釋到種的序列較少,并且種與種之間的差異性不明顯,會存在測序錯誤等情況,同時所有細菌OTUs中有50%左右的未能注釋到屬,因此本文選擇對注釋到目和科的物種進行研究。
圖3 不同樣品的OTU及分類注釋結果的熱圖Fig.3 Heatmap analysis of OTU and species annotation from different samples注:熱圖中的顏色表示一個OTU或分類所包括的序列數(shù)目的高低水平(顏色越深表示越高), 數(shù)字表示某個樣品中某個OTU或分類所包括的序列數(shù);k:界,p:門,c:綱,o:目,f:科,g:屬,s:種。
多樣性分析用于分析樣品內(nèi)的群落多樣性,主要包含指標:稀釋曲線、Chao1指數(shù)曲線、Shannon指數(shù)曲線,結果分別見圖4、圖5、圖6。
由稀釋曲線(見圖4)可知,當隨機抽取的測序條數(shù)逐漸增加時,4個樣品的稀釋曲線均逐漸趨于平緩,說明4個樣品的測序結果基本反映了魚露樣品中的細菌群落結構。分析Chao1指數(shù)曲線(見圖5)可知,當樣品測序量在5000~10000條之間時,Chao1指數(shù)曲線逐漸趨于平緩,當達到10000條后,Chao1指數(shù)曲線幾乎保持水平,進入一個穩(wěn)定的平臺期。韓亞飛等[12]指出Chao1指數(shù)值越高表明群落物種的豐富度越高,而除了發(fā)酵0個月的樣品外,各發(fā)酵魚露樣品的Chao1指數(shù)值是發(fā)酵2個月>發(fā)酵8個月>發(fā)酵12個月,這說明發(fā)酵過程中魚露中細菌群落的豐富度逐漸降低。樣品的Shannon指數(shù)曲線(見圖6)顯示,當樣品的測序條數(shù)在1000左右時,曲線逐漸趨向平坦,并進入平臺期,說明測序數(shù)據(jù)量足夠大,可以反映樣品中絕大多數(shù)的細菌群落信息。韓亞飛等[12]指出Shannon指數(shù)值越大說明群落多樣性越高,魚露發(fā)酵0個月時Shannon指數(shù)值為4.4,魚露發(fā)酵2個月、8個月、12個月時Shannon指數(shù)值為6.7,其中8個月時Shannon指數(shù)值最高,這說明2個月、8個月、12個月的發(fā)酵魚露中細菌群落多樣性明顯比0個月的高。
圖4 稀釋曲線Fig.4 Rarefaction curve
圖5 Chao1指數(shù)曲線Fig.5 Chao1 index curve
圖6 Shannon指數(shù)曲線Fig.6 Shannon index curve
β-多樣性分析是對不同樣品的微生物群落構成進行比較,反應不同樣品的微生物群落構成的主要指標有加權Unifrac(定量)/非加權Unifrac(定性)距離、主成分/主坐標分析、樣品聚類分析等,其中以(非)加權Unifrac距離繪制的β-多樣性指數(shù)熱圖的結果見圖7。
圖7 β-多樣性指數(shù)熱圖Fig.7 β-diversity heatmap analysis注:同一方格中,上下兩個值分別代表加權unifrac和非加權unifrac距離。
由于兩個樣品間的相異系數(shù)是由加權Unifrac距離來衡量,距離越小,表示這兩個樣品在物種多樣性方面存在的差異越小。圖7顯示,魚露發(fā)酵0個月和2個月之間、0個月和8個月之間、0個月和12個月之間的加權Unifrac和非加權Unifrac距離值分別為(0.249,0.432)、(0.332,0.423)、(0.318,0.426);發(fā)酵12個月和2個月之間、12個月和8個月之間的加權Unifrac和非加權Unifrac距離值分別為(0.134,0.357)、(0.042,0.328)??梢?魚露發(fā)酵前(0個月)的細菌群落結構與后期(8、12個月)的群落結構差異較大;而發(fā)酵后期,8個月與12個月的群落結構差異較小。這說明在魚露發(fā)酵初期,細菌群落結構的變化較大,越到發(fā)酵后期,細菌群落結構的變化越小,群落結構趨于穩(wěn)定。
2.4.1 發(fā)酵魚露中細菌群落在各分類水平上的組成 根據(jù)物種注釋結果顯示,發(fā)酵魚露中細菌多達37個門、99個綱、144個目、200個科,有的序列甚至無法鑒定,細菌群落構成復雜。由于部分分類相對豐度較低,推測這些細菌屬污染的細菌,因此我們只選取了相對豐度排名前10的分類進行統(tǒng)計,并根據(jù)FQ-PCR測得細菌總數(shù)算出它們各自的數(shù)量,結果分別見圖8、圖9、圖10、表2、表3、表4。
表2 魚露發(fā)酵各階段中細菌群落在門分類水平上的各自數(shù)量(lg(個/mL))Table 2 Comparison of the number of bacteria in fish sauce at the level of phylum(lg(個/mL))
表3 魚露發(fā)酵個階段中細菌群落在目分類水平上的各自數(shù)量(lg(個/mL))Table 3 Comparison of the number of bacteria in fish sauce at the level of order
表4 魚露發(fā)酵個階段中細菌群落在科分類水平上的各自數(shù)量(lg(個/mL))Table 4 Comparison of the number of bacteria in fish sauce at the level of family
圖8 魚露中細菌群落在門分類水平上的組成Fig.8 Composition of bacterial community in fish sauce at the level of phylum
圖9 魚露中細菌群落在目分類水平上的組成Fig.9 Composition of bacterial community in fish sauce at the level of order
圖10 魚露中細菌群落在科分類水平上的組成Fig.10 Composition of bacterial community in fish sauce at the level of family
圖8和表2顯示,通過16S rDNA高通量測序分析得出發(fā)酵魚露在各個發(fā)酵階段中的細菌主要分布在變形菌門和厚壁菌門,其他優(yōu)勢菌所屬的門還有酸桿菌門、放線菌門、擬桿菌門。變形菌門的相對豐度逐漸增加,由發(fā)酵0個月的43.69%逐漸增加到發(fā)酵12個月的68.79%,相對豐度波動相對較小,數(shù)量從1.07×105個/mL逐漸減少到2.82×104個/mL。而厚壁菌門的相對豐度波動相對較大,從發(fā)酵0個月的52.63%逐漸減少到發(fā)酵12個月的19.39%,數(shù)量從1.29×105個/mL逐漸減少到7.94×103個/mL,推測是由于魚露發(fā)酵是開放式的,不斷受到環(huán)境中各種微生物的污染,同時在發(fā)酵過程中要定時攪拌,促使魚鹽混勻,因此可能使得魚露中污染的細菌數(shù)量增加、厭氧微生物和兼性厭氧微生物的生長受到抑制。在好氧微生物和厭氧微生物的分類中,變形菌門中假單胞菌目、交替單胞菌目基本都是好氧微生物,而厚壁菌門中乳桿菌目、芽孢桿菌目的大部分是厭氧或兼性厭氧微生物[13]。
在目類水平上(見圖9、表3),魚露在各個發(fā)酵階段中的細菌主要分布在假單胞菌目、芽孢桿菌目、乳桿菌目。芽孢桿菌目的相對豐度在發(fā)酵初期降低明顯,從0個月的39.42%降低到2個月的11.84%,數(shù)量從9.55×104個/mL降到4.17×103個/mL;假單胞菌目的相對豐度變化不明顯,0個月和12個月相對豐度分別為31.31%和 30.12%,但數(shù)量呈降低趨勢,數(shù)量分別為7.59×104個/mL和1.23×104個/mL;乳桿菌目的相對豐度呈先增加后降低、再趨于平緩的變化趨勢,數(shù)量從0個月的2.82×104個/mL降低到12個月的3.24×103個/mL。
圖10和表4顯示,在科類水平上,主要分布在葡萄球菌科、莫拉菌科、鏈球菌科、乳桿菌科、假單胞菌科,其中芽孢桿菌目所包含的葡萄球菌科的相對豐度波動較大,從0個月的38.77%降低到12個月的9.12%,數(shù)量從9.33×104個/mL降到3.72×103個/mL,數(shù)量降低明顯,降低了1.4個數(shù)量級,推測可能是由于發(fā)酵初期高鹽抑制了大部分的葡萄球菌科細菌的生長。假單胞菌目所包含的莫拉菌科相對豐度和數(shù)量在發(fā)酵初期迅速降低,然后再緩慢增加,而假單胞菌科的相對豐度一直呈增長趨勢,數(shù)量變化不明顯,數(shù)量在3.0~5.9×103個/mL范圍內(nèi),但這些細菌部分是食品的優(yōu)勢腐敗菌,會產(chǎn)生腥臭味、硫化氫味等不良氣味,影響魚露的風味和安全性[14]。乳桿菌目所包含的鏈球菌科、乳桿菌科的相對豐度呈先降低、后增加、再降低、后再趨于平緩的曲線增減,而數(shù)量分別從0個月的1.32×104個/mL、7.59×103個/mL降低到12個月的1.48×103個/mL、1.38×103個/mL,數(shù)量降低了1個數(shù)量級左右,這說明鏈球菌科、乳桿菌科能在高鹽的環(huán)境中生存,并隨著發(fā)酵時間的延長逐漸成為優(yōu)勢菌。
2.4.2 發(fā)酵魚露中乳酸菌類群細菌的構成分析 16S rDNA測序結果顯示發(fā)酵魚露中細菌多達238個屬,同時在測序得到的所有細菌OTUs中有50%左右的未能注釋到屬,細菌群落結構復雜。而現(xiàn)有報道指出一些耐鹽或嗜鹽的乳酸菌是發(fā)酵魚產(chǎn)品的優(yōu)勢菌群,產(chǎn)生的蛋白酶能將蛋白質(zhì)分解成活性肽、游離氨基酸等風味物質(zhì),對魚露發(fā)酵過程中風味的形成具有重要作用[2-3]。參考《乳酸菌:生物學基礎及應用》[15]中Bergey氏細菌鑒定手冊、《乳酸細菌分類鑒定及實驗方法》[16]以及盛海圓等[17]對于乳酸菌類群細菌的分類,對魚露發(fā)酵過程中乳酸菌的菌屬構成進行分析,結果見表5。
表5 魚露中乳酸菌類群細菌的構成分析結果Table 5 Composition of lactic acid bacteria in fish sauce at different fermentation stages
由表5可知,檢測到魚露樣品中有10個為乳酸菌的屬或其中部分種為乳酸菌的屬,主要分布在鏈球菌屬、乳桿菌屬、葡萄球菌屬,他們的相對豐度分別由發(fā)酵0個月的5.34%、3.05%、2.11%逐漸變?yōu)?2個月的3.62%、3.23%、1.20%,鏈球菌屬、乳桿菌屬的相對豐度變化呈先增加、后降低、再趨于平緩的曲線增減趨勢,其中鏈球菌屬的相對豐度波動較大;葡萄球菌屬呈先降低、后增加的趨勢,但相對豐度波動較小。乳桿菌很少有致病菌,是食品發(fā)酵的主要功能菌,部分種能水解精氨酸產(chǎn)NH3;鏈球菌部分種能產(chǎn)生吡咯烷基芳酰胺酶以及水解精氨酸、馬尿酸;葡萄球菌屬能在高鹽的環(huán)境中生存,具有較強耐鹽性,部分種具有蛋白酶和脂肪酶活性,對于發(fā)酵肉制品風味成分的形成具有重要作用[18]。同時在魚露中還檢測到低豐度的芽孢桿菌屬、片球菌屬、氣球菌屬、營養(yǎng)缺陷菌屬、魏斯氏菌屬、奇異菌屬、雙歧桿菌屬等乳酸菌。其中芽孢桿菌屬中少數(shù)種屬于乳酸菌,如凝結芽孢桿菌(Bacilluscoagulans)、斯密氏芽孢桿菌(Bacillussmithii)、嗜熱脂肪芽孢桿菌(Bacillusstearothermophilus)等,部分種會產(chǎn)吲哚;氣球菌屬是鹽水腌肉的菌系,能耐高鹽環(huán)境,部分種會產(chǎn)亮氨酸氨肽酶、吡咯烷基芳酰胺酶;營養(yǎng)缺陷菌屬部分種會產(chǎn)亮氨酸氨肽酶;雙歧桿菌屬的部分種能產(chǎn)生丙氨酸、纈氨酸、天冬氨酸、蘇氨酸等多種氨基酸;片球菌屬、魏斯氏菌屬、奇異菌屬中部分種能水解精氨酸產(chǎn)NH3[16]。此外,由于50%左右的OTUs無法注釋到屬,在這些OTUs中是否是乳酸菌不得而知。
2.4.3 發(fā)酵12個月時魚露中細菌群落分析 對每個樣品進行物種分類,可以篩選出特別關注的物種。現(xiàn)有研究表明微生物對于傳統(tǒng)發(fā)酵魚產(chǎn)品風味和生物活性成分的形成具有重要作用[19]。選擇對發(fā)酵12個月的魚露中細菌進行物種分類樹統(tǒng)計(默認在屬水平上,選擇樣品中最大相對豐度前10的屬所對應的結果進行展示),見圖11。
發(fā)酵到12個月時,發(fā)酵魚露的細菌在變形菌門(68.79%)、厚壁菌門(19.39%)、放線菌門(3.19%)、酸桿菌門(2.19%)、擬桿菌門(1.17%)等10個門類中均有分布,其中主要是變形菌門和厚壁菌門(見2.4.1)。變形菌門中以乳桿菌屬、鏈球菌屬、葡萄球菌屬為主,而厚壁菌門中以不動菌屬、嗜冷桿菌屬、假單胞菌屬為主(見圖11)。多數(shù)研究報道指出乳桿菌屬中部分細菌即有產(chǎn)蛋白酶能力又能對發(fā)酵食品風味的形成發(fā)揮重要作用,所以國內(nèi)外學者都致力于從魚露里篩選乳酸菌,并將它們應用于魚露發(fā)酵中[20]。但在本試驗中,乳酸菌目(7.93%)相對豐度和數(shù)量并不算高,推測是由于魚露在發(fā)酵過程中要定時攪拌,促使魚鹽混勻,使得乳桿菌屬中部分厭氧和兼性厭氧細菌的生長受到抑制。
圖11 發(fā)酵12個月的發(fā)酵魚露中細菌群落的特定物種分類樹Fig.11 Classification tree of bacterial community in fish sauce during fermentation from 12th month注:不同顏色的圓圈表示不同的分類水平,對應左側圖例;圓圈的大小表示該分類的相對豐度大小;分類名下方的兩個數(shù)字均表示相對豐度百分率,前者表示該分類占所有分類中的百分率(相對豐度前10的屬占此分類的百分率);后者則表示該分類占所選取的分類中的百分率(以相對豐度前10的屬為百分制計)。
現(xiàn)有研究在魚露中發(fā)現(xiàn)的主要細菌還包括嗜鹽四鏈球菌屬、枝芽孢桿菌屬、芽孢桿菌屬、鹽單胞菌屬、片球菌屬等[3,21-22]。本文也發(fā)現(xiàn)了高豐度的葡萄球菌屬、乳桿菌屬、鏈球菌屬、嗜鹽單胞菌屬、片球菌屬,但未發(fā)現(xiàn)嗜鹽四鏈球菌、枝芽孢桿菌屬等細菌,這可能是由于原料、發(fā)酵環(huán)境等不同導致細菌菌相的差異。
研究表明,用16S rDNA技術檢測出在35 ℃發(fā)酵條件下制得的魚露中細菌多達37個門、238個屬,能夠比較真實地反映其細菌群落構成,較篩選和計數(shù)微生物的平板計數(shù)法具有明顯優(yōu)勢。發(fā)酵前期菌相結構差異性較大,到發(fā)酵后期,菌相結構逐漸趨于穩(wěn)定。在魚露發(fā)酵各階段都是優(yōu)勢菌的菌屬有鏈球菌屬、乳桿菌屬、不動菌屬、假單胞菌屬、希瓦氏菌屬。在發(fā)酵過程中,魚露中假單胞菌目中嗜冷桿菌屬、不動菌屬、假單胞菌屬的相對豐度波動較大,而耐鹽和嗜鹽的乳酸菌目中鏈球菌屬、乳桿菌屬的相對豐度波動較小。