• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    傳統(tǒng)清香型白酒發(fā)酵過程中細菌群落演替及自組裝機理

    2023-12-14 12:45:22喬美靈任宇婷張桂蓮楊玉珍劉夢宇孫子羽陳忠軍滿都拉
    食品科學(xué) 2023年22期
    關(guān)鍵詞:物種差異

    喬美靈,任宇婷,張桂蓮,楊玉珍,劉夢宇,孫子羽,陳忠軍,滿都拉,*

    (1.內(nèi)蒙古農(nóng)業(yè)大學(xué)食品科學(xué)與工程學(xué)院,內(nèi)蒙古 呼和浩特 010018;2.內(nèi)蒙古馳原酒業(yè)有限公司,內(nèi)蒙古 呼和浩特 011500;3.內(nèi)蒙古河套酒業(yè)集團股份有限公司,內(nèi)蒙古 巴彥淖爾 015400)

    白酒是我國釀造歷史悠久的傳統(tǒng)酒類,開放的生產(chǎn)方式使其在發(fā)酵過程中涉及多種微生物共同參與。研究顯示,白酒中含有上千種風味物質(zhì)[1],這與其由多種微生物參與發(fā)酵密不可分[2]。近年來,基于純培養(yǎng)[3]、磷脂脂肪酸分析[4]、變性梯度凝膠電泳[5]和高通量測序[6-9]等技術(shù)對白酒釀造微生物物種組成、群落結(jié)構(gòu)及群落演替開展大量研究,為后續(xù)探究微生物與白酒釀造之間的關(guān)系提供理論參考。由于參與白酒釀造的微生物受地理、氣候及工藝影響較大,一些特定地區(qū)白酒發(fā)酵過程中微生物群落自組裝,及其朝著有序、有規(guī)律方向演替的推動因素還需進一步了解。

    鑒別微生物種類、準確了解白酒發(fā)酵過程中微生物群落結(jié)構(gòu)及其隨時間變化的規(guī)律是研究白酒發(fā)酵微生物的主要內(nèi)容。高分辨的菌種注釋為揭示白酒微生物群落結(jié)構(gòu)、演變及自組裝提供了基礎(chǔ)。與二代短片段測序相比,以PacBio SMRT為代表的三代測序技術(shù)具有長讀長、高通量的優(yōu)點,通過其特有的CCS(circular consensus sequencing)模式,可將同一條序列多次重復(fù)測序后的結(jié)果相互校正,從而獲得測序質(zhì)量非常高的一致性序列[10]。借助該技術(shù),研究人員可以進一步獲得種水平的物種組成,得到更好的物種注釋分辨率,進一步深入解析群體樣本的物種組成情況[11]。

    環(huán)境因素的變化可以驅(qū)動微生物群落組成變化[12]。以往有關(guān)白酒釀造微生物自組裝機制的研究通?;凇罢麄€群落”進行評估[13]。而有學(xué)者認為微生物群落自組裝過程的真實情況更可能是在基因型和種群水平上發(fā)生,而不是整個群落,Ning Daliang等[14]提出基于系統(tǒng)發(fā)育的零模型分析方法——iCAMP,這為量化生態(tài)過程在控制微生物群落多樣性和演替中的相對重要性提供了可靠的工具。研究揭示白酒發(fā)酵過程中微生物群落結(jié)構(gòu)組成、演變及自組裝規(guī)律是深入了解控制群落結(jié)構(gòu)潛在因素以及白酒風味物質(zhì)形成機理的關(guān)鍵,也是優(yōu)化多種微生物混合發(fā)酵過程,提高產(chǎn)品品質(zhì)的重要理論基礎(chǔ)。因此,本研究以中國傳統(tǒng)清香型白酒為研究對象,利用PacBio SMRT測序技術(shù)結(jié)合統(tǒng)計學(xué)方法深入解析清香型白酒發(fā)酵過程中細菌物種組成及群落結(jié)構(gòu),進一步利用iCAMP模型開展清香型白酒發(fā)酵過程中群落自組裝機理的研究。

    1 材料與方法

    1.1 材料與試劑

    酒醅樣品采集自內(nèi)蒙古一清香型白酒廠(地缸發(fā)酵),采集時間為2020年11、12月。分別采集發(fā)酵第0、3、5、7、10、12、15、18、21、24、28天的酒醅,每個時間點采集3 個地缸。分別采集地缸的上、中以及底層酒醅(共500 g),混合為一個樣品后立即用無菌密封袋包裝保存于-80 ℃冰箱。

    1.2 儀器與設(shè)備

    YXQ-LS-100SII立式壓力蒸汽滅菌 上海博訊實業(yè)有限公司醫(yī)療設(shè)備廠;DHG-9053A電熱恒溫鼓風干燥箱上海一恒科學(xué)儀器有限公司;SIGMA3-18KS高速冷凍離心機 德國Sigma公司;Veriti梯度聚合酶鏈式反應(yīng)儀美國應(yīng)用生物系統(tǒng)公司;SE602F電子天平 奧豪斯儀器(上海)有限公司;DYY-6D電泳儀 北京六一生物科技有限公司;XW-80A旋渦儀 上海馳唐電子有限公司。

    1.3 方法

    1.3.1 理化指標測定

    溫度測定:采用測溫計直接測量取樣點附近的溫度,30 s不變溫,記錄數(shù)值;水分含量、pH值和酸度的測定參考賈麗艷等[15]的方法;淀粉含量的測定參考李大和等[16]的方法;蛋白質(zhì)含量測定參考GB 5009.5—2010《食品中蛋白質(zhì)的測定》,采用凱氏定氮法;游離氨基酸態(tài)氮含量測定參考GB 5009.239—2016《食品酸度的測定》,采用比色法;乙醇含量的測定參考劉曉等[17]的方法。

    1.3.2 酒醅微生物基因組的提取及測序

    酒醅微生物基因組采用試劑盒直接提取。稱取0.5 g酒醅樣品,按照TIANamp Soil DNA Kit說明書進行。利用含有Barcode 引物27 F(5’-AGAGTTTGATCCTGGCTCAG-3’)和1492R(5’-GGTTACCTTGTTACGACTT-3’)擴增細菌16S rRNA V1~V7區(qū)域[18]。產(chǎn)物經(jīng)純化、定量和均一化形成測序文庫(SMRT Bell),質(zhì)檢合格的文庫用PacBio Sequel進行測序。

    1.4 數(shù)據(jù)處理與分析

    1.4.1 下機數(shù)據(jù)處理

    利用Smrtlink version8.0軟件從下機數(shù)據(jù)中導(dǎo)出CCS文件,根據(jù)Barcode序列識別不同樣品的數(shù)據(jù)并轉(zhuǎn)化為fastq格式數(shù)據(jù)。利用UCLUST軟件劃分序列(相似度100%)為無重復(fù)的單一16S rRNA基因全長序列集。使用USEARCH(version10.0)以相似性97%對序列進行操作分類單元(operational taxonomic units,OTU)劃分,用Chimera Slayer去掉含嵌合體的OTU序列,并抽取代表性序列與SILVA(Release132,http://www.arb-silva.de)數(shù)據(jù)庫進行比對注釋,確定各OTU的分類學(xué)信息,進而計算各樣品在屬和種水平的相對含量。使用eggNOG數(shù)據(jù)庫(http://eggnog.Embl.de/)對細菌群進行功能預(yù)測。

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

    通過Shannon指數(shù)、Simpson指數(shù)[19]、Chao1指數(shù)、覆蓋率以及Shannon-Wiener曲線[20]評估各樣本α多樣性和測序質(zhì)量。基于OTU水平采用加權(quán)UniFrac算法和Bray-Curtis算法計算樣本間距離矩陣,通過構(gòu)建UPGMA聚類樹[21]和主坐標分析(principal co-ordinates analysis,PCoA)研究不同發(fā)酵時間樣本間差異。進一步通過LEfSe分析評估各發(fā)酵階段的差異物種[22]。通過計算Spearman相關(guān)系數(shù)(r),選擇|r|>0.60、P<0.05作為網(wǎng)絡(luò)強相關(guān)節(jié)點的條件,構(gòu)建相關(guān)性網(wǎng)絡(luò)圖并繪制功能及理化因子與物種間的相關(guān)性熱圖。基于典范對應(yīng)分析(canonical correspondence analysis,CCA)將微生物群落與理化因子進行關(guān)聯(lián)分析。微生物群落變化使用β多樣性分解進行評估[23]。計算歸一化隨機比(normalized stochasticity ratio,NST)和iCAMP分析細菌群落自組裝機制[24]。使用R語言(v4.1.2)進行統(tǒng)計分析和繪圖,采用方差分析進行組間差異顯著性檢驗。

    2 結(jié)果與分析

    2.1 發(fā)酵過程中細菌群落α多樣性

    經(jīng)測序共獲得413 198 條細菌CCS序列(均長度為1 500 bp),以97%相似性劃分共得到1 761 個OTU。Shannon指數(shù)和Simpson指數(shù)表征樣品中群落的多樣性信息,Shannon指數(shù)越大、Simpson指數(shù)越小,表示群落多樣性越高[25]。如圖1A、B所示,細菌群落的Shannon指數(shù)0 d最高隨后逐漸降低,Simpson指數(shù)0 d最低隨后逐漸升高,說明隨發(fā)酵進行細菌群落多樣性逐漸降低,發(fā)酵起始(0 d)時上述兩個指數(shù)均與發(fā)酵結(jié)束(28 d)時存在顯著差異(P<0.05)??傮w來看,發(fā)酵12 d后是發(fā)酵過程中細菌多樣性變化的關(guān)鍵時間段。Chao1指數(shù)表征樣品中群落的豐度信息,其值越大代表物種總數(shù)越多。本研究中,隨發(fā)酵的進行,細菌菌群的Chao1指數(shù)總體呈下降趨勢,這表明隨著發(fā)酵進行物種數(shù)逐漸降低,除0 d與10 d的Chao1指數(shù)差異顯著(P<0.05)外,其余無顯著差異(圖1C)。如圖1D所示,覆蓋率均大于0.99,說明樣本文庫的覆蓋率較高。Shannon-Wiener曲線隨著測序量的擴大,最終均達到飽和狀態(tài)(圖1E),測序深度滿足要求。從α多樣性(Simpson、Shannon指數(shù)和Chao1指數(shù))變化可以看出,由于入缸時帶入不同來源的微生物[26],細菌群落的微生物多樣性在0 d達到最大,之后由于發(fā)酵過程的進行,微生物需要適應(yīng)環(huán)境的變化以及菌群之間的相互作用導(dǎo)致細菌多樣性發(fā)生動態(tài)變化。

    圖1 發(fā)酵過程中細菌菌群Shannon指數(shù)(A)、Simpson指數(shù)(B)、Chao1指數(shù)(C)、覆蓋率(D)和Shannon-Wiener曲線(E)Fig.1 Shannon index (A),Simpson index (B),Chao1 index (C),coverages (D) and Shannon-Wiener curve (E) of bacterial flora during fermentation process

    2.2 發(fā)酵過程中細菌群落組成分析

    如圖2A所示,利用所測得的OTU經(jīng)物種注釋共得到87 個菌屬,其中乳桿菌屬(Lactobacillus)、醋酸菌屬(Acetobacter)、片球菌屬(Pediococcus)和魏斯氏菌屬(Weissella)為優(yōu)勢屬(相對豐度大于1%)。發(fā)酵0~3 d,Acetobacter的相對豐度從35%下降至0.9%,Weissella從15%下降至0.5%,Pediococcus從10%下降至8%,而Lactobacillus的相對豐度從10%迅速上升至90%。發(fā)酵3~10 d,Acetobacter和Weissella的相對豐度一直維持在較低水平(相對豐度小于0.1%),Pediococcus的相對豐度從8%下降至1%。發(fā)酵10 d后,Lactobacillus成為絕對優(yōu)勢屬(相對豐度大于98%)。發(fā)酵12 d時,Acetobacter和其他屬的相對豐度呈現(xiàn)小幅度增加。明串珠菌屬(Leuconostoc)、葡糖桿菌屬(Gluconobacter)、芽孢桿菌屬(Bacillus)、駒形桿菌屬(Komagataelibacter)和擬桿菌屬(Bacteroides)在發(fā)酵0 d時相對豐度均較?。ㄐ∮?.5%)。

    圖2 發(fā)酵過程中細菌群落屬水平(A)和種水平(B)組成Fig.2 Composition of bacterial communities at the genus level (A) and species level (B) during the fermentation process

    在種水平上共得到124 個菌種,其中橋乳桿菌(L.pontis)、類消化乳桿菌(L.paralimentarius)、短乳桿菌(L.brevis)、耐酸乳桿菌(L.acetotolerans)、啤酒乳桿菌(L.cerevisiae)、蘋果醋桿菌(A.malorum)、戊糖片球菌(P.pentosaceus)、植物乳桿菌(L.plantarum)、希氏乳桿菌(L.hilgardii)和食竇魏斯氏菌(W.cibaria)為優(yōu)勢菌種(相對豐度大于1%)。在發(fā)酵0~3 d時,L.pontis、L.paralimentarius、L.brevis、L.cerevisiae和L.plantarum的相對豐度逐漸升高,而A.malorum、P.pentosaceus和W.cibaria的相對豐度逐漸下降至1%。發(fā)酵3~10 d時,L.cerevisiae和P.pentosaceus的相對豐度逐漸下降至1%,而L.pontis的相對豐度逐漸升高。發(fā)酵12~28 d時,L.acetotolerans相對豐度從1%上升至40%。發(fā)酵28 d時,L.acetotolerans和L.pontis成為優(yōu)勢菌。其余細菌在發(fā)酵0 d相對豐度較高,但隨著發(fā)酵過程的進行,相對豐度逐漸降低并維持在較低水平(相對豐度小于1%)(圖2B)。

    2.3 發(fā)酵過程中細菌群落β多樣性

    如圖3A所示,發(fā)酵0 d樣品與其他時間樣品存在顯著差異,成為單獨一個分支。其他時間樣品又分為3~10 d和12~28 d兩個分支,說明整個發(fā)酵可分為3 個階段:0(I)、3~10 d(II)和12~28 d(III)。PCoA結(jié)果(圖3B)顯示,PCo1的貢獻率為55.3%,權(quán)重最大,說明細菌群落隨發(fā)酵過程的推進主要沿PCo1軸方向發(fā)生變化,并且沿PCo1軸正方向相鄰時間點樣本的間距越來越小,說明群落結(jié)構(gòu)逐漸趨于一致,這與聚類分析結(jié)果一致。綜上,發(fā)酵第0、3、12天可能是反映整個發(fā)酵過程微生物群落動態(tài)變化的主要時間節(jié)點。進一步利用方差分析(圖3C)研究3 個發(fā)酵階段的差異性發(fā)現(xiàn),3 組組間差異大于組內(nèi)差異,且差異顯著(P<0.001),表明將發(fā)酵過程分為3 個發(fā)酵階段是合理的。

    圖3 發(fā)酵過程中細菌群落的β多樣性Fig.3 β-Diversity of bacterial flora during fermentation

    為進一步明確不同發(fā)酵階段中具有顯著差異物種信息,通過LEfSe分析研究各階段差異菌種。如圖3D所示,各發(fā)酵階段細菌群落的差異物種主要來自8 個屬,大部分為優(yōu)勢屬,但階段II中沒有差異菌屬。其中階段I有11 個顯著性差異物種,分別為A.malorum、W.cibaria、P.pentosaceus、巴氏醋桿菌(A.pasteurianus)、腸膜明串珠菌(L.mesenteroides)、氧化葡萄糖桿菌(G.oxydans)、白面粉乳桿菌(L.siliginis)、清酒乳酸桿菌(L.sakei)、中間駒形氏桿菌(K.intermedius)、地衣芽孢桿菌(B.licheniformis)和羅氏菌屬中未分類的種(uncultured_g_Roseburia)。階段II有4 個顯著性差異菌種,為L.paralimentarius、L.brevis、L.cerevisiae和L.plantarum。階段III有3 個顯著性差異菌種,分別為L.pontis、L.acetotolerans和L.hilgardii。階段I、II和III中A.malorum、L.paralimentarius和L.pontis的線性判別分析(linear discriminant analysis,LDA)評分最高,分別為其主要的差異菌種(圖3E)。

    2.4 差異菌種間相關(guān)性網(wǎng)絡(luò)分析

    通過計算Spearman相關(guān)系數(shù)研究發(fā)酵過程中差異細菌間的共存關(guān)系。選擇|r|>0.60、P<0.05作為網(wǎng)絡(luò)強相關(guān)節(jié)點的條件篩選出10 個菌種,進行繪圖,如圖4所示。細菌群落中L.acetotolerans與各發(fā)酵階段的差異細菌具有顯著相關(guān)性(P<0.05)。階段I的差異菌種之間均呈正相關(guān),其中A.malorum、W.cibaria、A.pasteurianus、K.intermedius、P.pentosaceus、G.oxydans和L.sakei,這些菌種具有較高的有機酸和乙酸生產(chǎn)能力,可以導(dǎo)致酒醅酸化,維持發(fā)酵的正常進行,也為乳酸乙酯的合成提供了前體物質(zhì);L.siliginis多分離于酸面團中;而L.mesenteroides和B.licheniformis可以釋放各種酶水解酒醅中的蛋白質(zhì)和淀粉,為發(fā)酵過程中微生物提供物質(zhì)及能量[27]。階段II的差異菌種之間呈正相關(guān),其中L.paralimentarius利用葡萄糖發(fā)酵產(chǎn)酸,L.plantarum可產(chǎn)丁酸[28],L.brevis具有較高的γ-氨基丁酸生產(chǎn)能力[29],而L.cerevisiae在白酒中的研究較少,其在白酒釀造過程中的具體作用還需進一步研究。階段III的3 個差異菌種之間也呈正相關(guān),L.pontis為異型發(fā)酵乳酸菌,可改善面包的風味和味道[30-31],L.acetotolerans可以促進具有水果味、花香和甜香酯類物質(zhì)風味的產(chǎn)生[32],L.hilgardii參與葡聚糖的產(chǎn)生[33],促使某些細菌的生長。發(fā)酵過程中微生物群落結(jié)構(gòu)的差異影響微生物之間的相互作用。不同發(fā)酵階段的差異菌種之間呈負相關(guān),同一發(fā)酵階段的差異菌種之間呈正相關(guān),這是可能由于各階段的“分工”不同導(dǎo)致微生物之間存在競爭和共生等關(guān)系[34]。

    2.5 發(fā)酵過程中環(huán)境因子變化及其與細菌群落間的關(guān)系

    如圖5所示,隨著發(fā)酵的進行,理化因子呈明顯動態(tài)變化。水分含量、酸度和蛋白含量呈現(xiàn)上升趨勢,淀粉含量和pH值呈現(xiàn)下降趨勢。而溫度呈現(xiàn)先上升后下降趨勢,氨基酸態(tài)氮和乙醇含量呈現(xiàn)先升高后降低再升高的變化規(guī)律。

    圖5 發(fā)酵過程中理化指標動態(tài)變化Fig.5 Dynamic changes in physicochemical indicators during fermentation

    基于CCA將微生物群落與理化因子進行關(guān)聯(lián)分析,排序軸CC1貢獻率為54.19%,排序軸CC2貢獻率為28.44%(圖6A)。pH值和淀粉含量主要影響階段I的細菌群落,其中淀粉含量的影響最大,在這一階段中淀粉分解成還原糖,充足的糖供應(yīng)有利于微生物的生長和代謝[35]。溫度和氨基酸態(tài)氮含量主要影響發(fā)酵階段II的細菌群落,其中溫度對細菌群落的影響最大,溫度的變化反應(yīng)酒醅發(fā)酵是否正常,這是微生物的生長代謝產(chǎn)生生物熱與酒醅發(fā)酵環(huán)境熱能傳遞平衡的結(jié)果[36]。酸度、水分和蛋白質(zhì)含量主要影響發(fā)酵階段III的細菌群落,其中水分含量影響較大,水分含量是釀酒過程中關(guān)鍵指標之一,酒醅所含水分含量影響發(fā)酵空間的氣密性,水分含量偏高導(dǎo)致發(fā)酵的氣密性差,抑制好氧菌的生長,增加厭氧菌的生長,容易發(fā)生酸??;如果酒醅水分含量偏低,對出酒率和酒品質(zhì)都有影響[37]。差異菌種與環(huán)境因子相關(guān)性(圖6B)表明,階段I的差異菌種主要與pH值呈正相關(guān),與酸度呈負相關(guān),其中L.mesenteroides和P.pentosaceus呈高度正相關(guān)(P<0.001)。階段II的差異菌種主要與溫度呈正相關(guān),其中L.plantarum和L.cerevisiae與溫度高度相關(guān)(P<0.001)。在階段III的差異菌種中L.acetotolerans與水分和酸度呈高度正相關(guān)(P<0.001),與pH 值、溫度呈高度負相關(guān)(P<0.001),這表明不同種類的細菌受不同環(huán)境影響較大,具體情況還需對菌株進一步分離驗證。

    圖6 細菌群落與理化因子CCA(A)和差異菌種與理化指標相關(guān)性熱圖(B)Fig.6 CCA of bacterial community and physicochemical factors (A) and heatmap of correlation between differential species and physicochemical indicators (B)

    2.6 發(fā)酵過程中細菌群落自組裝機理

    微生物群落變化通常歸于兩種現(xiàn)象:物種替換(不同樣本間物種交換和替代)和物種流失(物種的增加或減少)[38]。如圖7A、B所示,本研究的發(fā)酵過程中物種替換(圖7B)幾乎貢獻了大部分的細菌群落演替。隨著發(fā)酵過程的進行,物種替換呈現(xiàn)先上升后下降再上升的變化趨勢,與物種流失(圖7A)呈現(xiàn)相反的變化。階段I細菌群落的替換最低,可能是由于環(huán)境細菌的進入導(dǎo)致細菌多樣性增加,而階段II、III,由于在發(fā)酵過程中微生物的代謝活動導(dǎo)致釀造環(huán)境發(fā)生變化,從而使替換迅速升高??傊鎿Q是群落演替的主要機制,同時流失對于細菌群落結(jié)構(gòu)的演替也至關(guān)重要。不同發(fā)酵階段的替換發(fā)生率呈現(xiàn)上升趨勢,表明發(fā)酵階段中的細菌群落中的物種在不斷發(fā)生交換和替代,細菌群落因發(fā)酵過程的進行而發(fā)生演替。

    圖7 β多樣性分解(A、B)、NST評價(C)和自組裝動態(tài)變化Sankey圖(D)Fig.7 Decomposition of beta-diversity (A and B),NST evaluation (C) and Sankey diagram showing dynamic changes in self-assembly (D)

    為探索確定性和隨機性對微生物群落組裝的相對重要性,基于零模型的NST進行評價,如圖7C所示,隨機過程主導(dǎo)了整個發(fā)酵過程中細菌群落的組裝結(jié)果。發(fā)酵過程中細菌群落的絕大部分NST>0.5,表明整個過程主要受到隨機性過程的影響。而在階段II和III,少數(shù)細菌群落也會受到確定性(NST<0.5)的影響,表明隨著發(fā)酵過程的進行,微生物對原料的分解和其自身代謝活動使細菌群落受到隨機性和確定性過程的影響。此外,iCAMP方法將確定性和隨機性過程分為5 個特定的生態(tài)過程(圖7D)。確定性過程包括均勻選擇(homogeneous selection,HoS)和異質(zhì)選擇(heterogeneous selection,HeS),而隨機過程包括擴散限制、均勻擴散和漂移(drift,DR)[14]。其中DR是發(fā)酵過程中細菌群落的主要過程,其次是HoS和HeS。然而,不同發(fā)酵階段呈現(xiàn)不同的過程,階段I小部分為HoS,可能因為發(fā)酵0 d引入大量的環(huán)境微生物導(dǎo)致酒醅中某些微生物發(fā)生變化。階段II大部分為HoS,這種均勻選擇是由于一定的環(huán)境選擇壓力導(dǎo)致,階段III大部分為DR,這很大可能是因為發(fā)酵環(huán)境的隨機性擾動。綜上,隨著發(fā)酵進行,細菌群落自組裝主要受到隨機過程中DR的影響。

    結(jié)合理化因子與細菌群落之間的相關(guān)性,因此推測淀粉含量是影響階段I的細菌群落演替的主要因素,溫度和水分含量分別是影響階段II和階段III的細菌群落演替的重要因素。

    3 討論

    微生物多樣性賦予傳統(tǒng)發(fā)酵食品多樣的菌群結(jié)構(gòu),其中核心微生物既為風味賦予特征性與多樣性,也為產(chǎn)品的質(zhì)量提供穩(wěn)定性[39]。精準定位參與傳統(tǒng)發(fā)酵食品中微生物物種、揭示發(fā)酵過程中微生物群落結(jié)構(gòu)變化及群落自組裝機理對進一步提升產(chǎn)品質(zhì)量及安全性意義重大。本實驗利用PacBio SMRT測序技術(shù)開展傳統(tǒng)清香型白酒發(fā)酵過程中細菌群落的多樣性及結(jié)構(gòu)演替規(guī)律的研究。

    在發(fā)酵過程中,細菌的多樣性和豐度表現(xiàn)出明顯的特征和動態(tài)變化。高豐度優(yōu)勢微生物通常被認為是食品發(fā)酵的重要組成部分[40]。經(jīng)注釋分析,在本研究中有4 個廣泛分布的細菌屬(Lactobacillus、Acetobacter、Pediococcus和Weissella),這與其他研究報道結(jié)果[41]有所差異,這可能與地區(qū)原料環(huán)境不同有關(guān)。擴增子全長測序能夠更精確地定位物種信息,也為進一步分析其功能奠定基礎(chǔ)。本研究中優(yōu)勢物種主要有L.pontis和L.acetotolerans,這與山西汾酒[42]、牛欄山二鍋頭[43]和青稞酒[44]等清香型白酒釀造系統(tǒng)中報道的優(yōu)勢乳酸菌種有明顯差異。此外,發(fā)酵0、3 d和12 d是細菌α多樣性、β多樣性和物種組成動態(tài)變化的主要時間節(jié)點。

    進一步根據(jù)細菌菌群的差異將整個發(fā)酵過程可分為3 個階段:I階段(0 d)、II階段(3~10 d)和III階段(12~28 d)。不同發(fā)酵階段具有不同的標志性差異菌種,這些差異菌種一部分來自優(yōu)勢菌種,另一部分來自稀有菌種。其中A.malorum(階段I)、L.paralimentarius(階段II)和L.pontis(階段III)是各發(fā)酵階段主要的差異菌種。微生物群落是在資源競爭、營養(yǎng)共生和群體感應(yīng)等相互作用下形成的復(fù)雜生態(tài)系統(tǒng)[45]。L.acetotolerans、L.pontis與其他菌種均相關(guān),不同發(fā)酵階段中差異菌種之間呈負相關(guān),屬于同一發(fā)酵階段的差異菌種之間呈正相關(guān),說明其中各物種在發(fā)酵過程中發(fā)揮的作用不同,此外還存在物種之間的相互競爭與共生等關(guān)系[43]。

    清香型白酒的釀造屬于固態(tài)發(fā)酵法,是邊糖化邊發(fā)酵工藝[46]。糖化主要是淀粉等多糖酶解轉(zhuǎn)化形成葡萄糖等可發(fā)酵性糖,發(fā)酵主要完成葡萄糖到乙醇及微量風味組分的轉(zhuǎn)變[47]。階段I中淀粉含量高,微生物繁殖快,短時間內(nèi)生成大量發(fā)酵熱,導(dǎo)致酒醅溫度升高,微生物利用氧氣進行生長繁殖及各類代謝活動產(chǎn)生CO2,引起pH值降低,酸度升高。階段II中隨著淀粉被消耗,微生物繁殖速度降低,酒醅溫度下降,厭氧及兼性厭氧細菌微生物的大量增殖生成多種有機酸,導(dǎo)致酒醅酸度增加,適量的酸能增加白酒的后味,使酒體更加豐滿。在階段III中一些具有酯化作用的微生物代謝產(chǎn)生較多的酶,對階段I生成的酸和醇產(chǎn)生催化酯化反應(yīng)的作用,酸度的增加有利于酯類物質(zhì)的形成,酯類物質(zhì)的生成會消耗酸類物質(zhì),導(dǎo)致pH值逐步降低。酯類的生成會消耗乙醇以及厭氧微生物利用乙醇進行生長繁殖[48],造成乙醇含量呈現(xiàn)動態(tài)變化。酒醅中的水分一部分來自于釀造用水,另一部分來自于釀酒微生物發(fā)酵代謝[49],水分為整個發(fā)酵順利進行奠定基礎(chǔ)。差異菌種與理化因子關(guān)聯(lián)分析發(fā)現(xiàn),階段I和II的微生物與淀粉含量、pH值和溫度正相關(guān),階段III的差異菌種主要與酸度、水分含量和氧化還原電位呈正相關(guān)。因此,發(fā)酵環(huán)境對發(fā)酵過程影響顯著[50]。

    發(fā)酵過程中,物種替換和物種流失兩種現(xiàn)象導(dǎo)致微生物群落多樣性的形成,而影響其形成和自組裝的生態(tài)進程包括基于環(huán)境選擇和生物競爭生態(tài)位理論的確定性過程,以及基于擴散和隨機出生死亡中性理論的隨機性過程[51]。清香型白酒發(fā)酵過程中大部分細菌群落幾乎都是替換現(xiàn)象[13],基于iCAMP群落組裝進一步挖掘發(fā)現(xiàn)主要為隨機性過程中的生態(tài)DR[52],即入缸時(階段I)帶入的微生物定植并占據(jù)酒曲微生物的生態(tài)位造成細菌菌群的隨機出生、死亡和繁殖,這種現(xiàn)象可能來源于淀粉含量的影響,白酒釀造中多數(shù)細菌依賴真菌降解淀粉產(chǎn)生的還原糖作為碳源,進行生長繁殖和代謝。釀造過程中因為原料和酒曲的不同導(dǎo)致所產(chǎn)生的還原糖不同,進而導(dǎo)致細菌所利用的碳源不同,最終影響細菌群落的結(jié)構(gòu)差異[53]。因此淀粉含量是影響階段I細菌群落自組裝的主要因素。溫度可以在短時間內(nèi)影響白酒釀造過程中功能微生物群落結(jié)構(gòu)(階段II),導(dǎo)致發(fā)酵環(huán)境發(fā)生擾動,細菌受到均勻選擇促使一些細菌群落發(fā)生系統(tǒng)發(fā)育性聚集,使得更適應(yīng)這種環(huán)境微生物的相對豐度提高。水分含量的控制是代謝活動和生命的關(guān)鍵組成部分,階段III作為細菌呼吸作用的最后階段,各種芳香化合物(主要是酸和酯)通過物理變化(范德華相互作用等)與化學(xué)反應(yīng)(如還原、氧化和酯化等)實現(xiàn)平衡[54]。所以水分含量是導(dǎo)致階段III的細菌群落自組裝的主要因素。而且不同生物體對環(huán)境變化的反應(yīng)也存在很大的差異。在不同微生物種群中,分散能力、多樣化和對DR的敏感性也有本質(zhì)上的不同。這些因素最終造成發(fā)酵過程由選擇的確定性過程逐漸向DR主導(dǎo)的隨機性過程轉(zhuǎn)變。

    綜上所述,由于開放的發(fā)酵環(huán)境和動態(tài)的發(fā)酵過程,清香型白酒發(fā)酵過程中所涉及的微生物群落表現(xiàn)出明顯的多樣性和動態(tài)性。本研究通過對清香型白酒發(fā)酵過程中細菌菌群物種進行精準定位,揭示細菌群落結(jié)構(gòu)的演替及自組裝機制,為進一步提升清香型白酒釀造工藝,優(yōu)化產(chǎn)品穩(wěn)定性和安全性提供了基礎(chǔ)理論。然而,本研究只關(guān)注了酒醅發(fā)酵過程中優(yōu)勢細菌菌種的演替,除此之外還存在大量的低豐度物種,這些物種隨著發(fā)酵過程的進行還可能具有一定的功能和對群落自組裝的貢獻。因此,未來應(yīng)結(jié)合測序技術(shù)的理論知識,將PacBio SMRT測序和傳統(tǒng)微生物培養(yǎng)方法相結(jié)合,根據(jù)酒醅發(fā)酵基質(zhì)環(huán)境和不同微生物生物學(xué)特征,進一步對中、低豐度的微生物物種進行分離和功能驗證。

    猜你喜歡
    物種差異
    物種大偵探
    物種大偵探
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    相似與差異
    音樂探索(2022年2期)2022-05-30 21:01:37
    找句子差異
    回首2018,這些新物種值得關(guān)注
    DL/T 868—2014與NB/T 47014—2011主要差異比較與分析
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    生物為什么會有差異?
    瘋狂的外來入侵物種
    亚洲,一卡二卡三卡| 999精品在线视频| 18禁观看日本| 欧美日本中文国产一区发布| av免费观看日本| 欧美少妇被猛烈插入视频| av电影中文网址| 天天躁夜夜躁狠狠躁躁| 亚洲精品美女久久av网站| freevideosex欧美| 丰满乱子伦码专区| 国产片特级美女逼逼视频| 精品人妻一区二区三区麻豆| 欧美日韩综合久久久久久| 日本av免费视频播放| 一级a做视频免费观看| 又黄又粗又硬又大视频| 51国产日韩欧美| 少妇精品久久久久久久| 看免费成人av毛片| 免费在线观看黄色视频的| 在线观看美女被高潮喷水网站| 欧美xxxx性猛交bbbb| 欧美97在线视频| 大陆偷拍与自拍| 日本av免费视频播放| 少妇被粗大猛烈的视频| 看免费av毛片| av免费在线看不卡| 欧美日韩视频精品一区| av.在线天堂| 男女无遮挡免费网站观看| 美女中出高潮动态图| 亚洲精品久久成人aⅴ小说| 久久久国产精品麻豆| 免费av不卡在线播放| 久久毛片免费看一区二区三区| 一区二区日韩欧美中文字幕 | 十八禁高潮呻吟视频| 国产色爽女视频免费观看| 精品一区二区三卡| 亚洲精华国产精华液的使用体验| 亚洲人成网站在线观看播放| 亚洲性久久影院| 午夜视频国产福利| 黄色怎么调成土黄色| av不卡在线播放| 一本久久精品| 色视频在线一区二区三区| 99热国产这里只有精品6| 男女啪啪激烈高潮av片| 国产精品偷伦视频观看了| 亚洲天堂av无毛| 欧美人与性动交α欧美精品济南到 | 一二三四在线观看免费中文在 | 欧美日本中文国产一区发布| 亚洲人成77777在线视频| 大陆偷拍与自拍| 亚洲国产精品专区欧美| 亚洲精品日本国产第一区| 国产淫语在线视频| 制服诱惑二区| 久久韩国三级中文字幕| av在线老鸭窝| 天天躁夜夜躁狠狠久久av| 九色成人免费人妻av| 亚洲国产色片| 伦精品一区二区三区| 街头女战士在线观看网站| 久久热在线av| 精品99又大又爽又粗少妇毛片| 女人久久www免费人成看片| 欧美bdsm另类| 99九九在线精品视频| 伊人久久国产一区二区| 国产av码专区亚洲av| 桃花免费在线播放| 一二三四在线观看免费中文在 | 啦啦啦在线观看免费高清www| 在现免费观看毛片| 成人综合一区亚洲| 亚洲av国产av综合av卡| 看免费av毛片| 欧美精品高潮呻吟av久久| 成年美女黄网站色视频大全免费| 久久久精品94久久精品| 美女视频免费永久观看网站| 黄色配什么色好看| 精品一区二区免费观看| 美女国产高潮福利片在线看| 国产不卡av网站在线观看| 最新的欧美精品一区二区| av天堂久久9| 激情视频va一区二区三区| 欧美人与善性xxx| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 亚洲美女搞黄在线观看| av.在线天堂| 少妇的逼好多水| 国产黄色视频一区二区在线观看| 捣出白浆h1v1| 亚洲综合色惰| 99久久精品国产国产毛片| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 亚洲av综合色区一区| 精品少妇黑人巨大在线播放| 亚洲国产av新网站| 国产精品99久久99久久久不卡 | 午夜视频国产福利| 国产精品久久久久久久电影| xxx大片免费视频| 我要看黄色一级片免费的| 免费黄频网站在线观看国产| 国产极品天堂在线| 99久久综合免费| 丁香六月天网| 亚洲精品日韩在线中文字幕| 美女福利国产在线| 久久婷婷青草| 欧美日韩精品成人综合77777| 午夜福利视频精品| 国产淫语在线视频| 免费在线观看黄色视频的| 一本久久精品| 人人妻人人爽人人添夜夜欢视频| 波多野结衣一区麻豆| 一二三四中文在线观看免费高清| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 一区二区三区四区激情视频| 国产在线一区二区三区精| 永久网站在线| 视频中文字幕在线观看| av在线老鸭窝| 亚洲熟女精品中文字幕| 亚洲精品日本国产第一区| 国产av一区二区精品久久| 国产精品偷伦视频观看了| 天堂8中文在线网| 在线亚洲精品国产二区图片欧美| 熟女电影av网| 久久久久视频综合| 亚洲精品第二区| 精品少妇黑人巨大在线播放| 大片电影免费在线观看免费| 久久精品久久久久久久性| 久久久久国产网址| 制服诱惑二区| 女性被躁到高潮视频| 国产午夜精品一二区理论片| 亚洲国产av新网站| 久久久精品区二区三区| 免费观看在线日韩| 人人妻人人添人人爽欧美一区卜| 亚洲国产看品久久| 看非洲黑人一级黄片| 国语对白做爰xxxⅹ性视频网站| 日韩欧美一区视频在线观看| 亚洲av成人精品一二三区| 日韩欧美精品免费久久| 国产精品国产av在线观看| 免费观看a级毛片全部| 精品少妇内射三级| 久久狼人影院| 午夜av观看不卡| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o| 久热这里只有精品99| 亚洲av日韩在线播放| 国产亚洲欧美精品永久| 蜜桃国产av成人99| 亚洲精品aⅴ在线观看| 国产片内射在线| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 午夜免费鲁丝| 成年女人在线观看亚洲视频| 搡女人真爽免费视频火全软件| 老司机影院成人| 亚洲精华国产精华液的使用体验| 亚洲一级一片aⅴ在线观看| 51国产日韩欧美| 性高湖久久久久久久久免费观看| 男女高潮啪啪啪动态图| 亚洲伊人色综图| 男女边摸边吃奶| 啦啦啦视频在线资源免费观看| 丝袜喷水一区| 五月伊人婷婷丁香| 狂野欧美激情性bbbbbb| 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 免费黄频网站在线观看国产| 日韩中字成人| 久久精品夜色国产| 国产熟女欧美一区二区| xxxhd国产人妻xxx| 国产成人精品一,二区| 国产精品偷伦视频观看了| 日韩欧美精品免费久久| 在线天堂中文资源库| 久久久国产一区二区| 国产无遮挡羞羞视频在线观看| 色94色欧美一区二区| 亚洲内射少妇av| 亚洲精品国产av蜜桃| 在线看a的网站| 观看美女的网站| 国产又爽黄色视频| 精品国产一区二区久久| 热re99久久精品国产66热6| 亚洲第一区二区三区不卡| 制服人妻中文乱码| 精品一区二区三卡| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 欧美日韩成人在线一区二区| 中文乱码字字幕精品一区二区三区| 亚洲在久久综合| 精品午夜福利在线看| 久久99热6这里只有精品| 亚洲成人av在线免费| 激情五月婷婷亚洲| 亚洲第一av免费看| 国产一区二区在线观看av| 免费黄色在线免费观看| 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 宅男免费午夜| 日韩av免费高清视频| 中文字幕制服av| 人成视频在线观看免费观看| 一本久久精品| 亚洲中文av在线| 欧美国产精品一级二级三级| 丝袜喷水一区| 精品国产一区二区久久| 18禁动态无遮挡网站| 伊人亚洲综合成人网| 毛片一级片免费看久久久久| 欧美日本中文国产一区发布| 两个人看的免费小视频| 午夜免费观看性视频| 国产成人91sexporn| 99热这里只有是精品在线观看| 久久女婷五月综合色啪小说| xxxhd国产人妻xxx| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男人操女人黄网站| 亚洲欧美一区二区三区国产| 久久久国产精品麻豆| 欧美性感艳星| 97超碰精品成人国产| 亚洲av.av天堂| 韩国高清视频一区二区三区| 久久99一区二区三区| xxxhd国产人妻xxx| 欧美日韩一区二区视频在线观看视频在线| 久久综合国产亚洲精品| 一区二区日韩欧美中文字幕 | 三级国产精品片| 免费女性裸体啪啪无遮挡网站| 啦啦啦中文免费视频观看日本| 考比视频在线观看| 亚洲欧洲日产国产| 日本欧美视频一区| 亚洲欧美成人综合另类久久久| 欧美日韩国产mv在线观看视频| av线在线观看网站| 免费大片黄手机在线观看| 久久久国产一区二区| 天天躁夜夜躁狠狠躁躁| 哪个播放器可以免费观看大片| 9热在线视频观看99| 亚洲成人av在线免费| 在线观看美女被高潮喷水网站| 国产色婷婷99| 中国三级夫妇交换| 乱码一卡2卡4卡精品| 观看av在线不卡| 久久人人爽人人片av| 午夜福利在线观看免费完整高清在| 美女福利国产在线| 性高湖久久久久久久久免费观看| 两个人看的免费小视频| kizo精华| 大片电影免费在线观看免费| 免费少妇av软件| 亚洲国产精品专区欧美| 亚洲国产精品国产精品| 精品国产乱码久久久久久小说| 秋霞在线观看毛片| 国产成人精品无人区| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 男女国产视频网站| 热99久久久久精品小说推荐| 男女国产视频网站| 国产又色又爽无遮挡免| 欧美国产精品一级二级三级| 人妻 亚洲 视频| 18禁裸乳无遮挡动漫免费视频| 欧美成人精品欧美一级黄| 婷婷色麻豆天堂久久| www日本在线高清视频| 亚洲av电影在线进入| xxxhd国产人妻xxx| 国产伦理片在线播放av一区| 狂野欧美激情性xxxx在线观看| 日韩精品免费视频一区二区三区 | 麻豆乱淫一区二区| 亚洲精品乱码久久久久久按摩| 一区二区三区乱码不卡18| 欧美激情极品国产一区二区三区 | 色网站视频免费| 久久97久久精品| 另类亚洲欧美激情| 一级毛片电影观看| 国产av一区二区精品久久| 亚洲综合色网址| av黄色大香蕉| 日日撸夜夜添| 午夜福利视频在线观看免费| 色吧在线观看| 老女人水多毛片| 国产精品秋霞免费鲁丝片| 亚洲国产成人一精品久久久| 日本午夜av视频| 麻豆精品久久久久久蜜桃| 男的添女的下面高潮视频| 女性被躁到高潮视频| 中文字幕人妻丝袜制服| 日韩免费高清中文字幕av| 成人亚洲欧美一区二区av| 91精品国产国语对白视频| 亚洲精品一二三| 亚洲精品国产色婷婷电影| 日韩电影二区| 久久久精品94久久精品| 97超碰精品成人国产| 国产视频首页在线观看| 国产又色又爽无遮挡免| 欧美性感艳星| 国产在线视频一区二区| 91精品国产国语对白视频| 亚洲国产日韩一区二区| 久久97久久精品| 在线天堂中文资源库| 在线观看美女被高潮喷水网站| 日韩成人伦理影院| 精品少妇久久久久久888优播| 久久久久久久久久成人| 亚洲欧美日韩卡通动漫| 大香蕉久久网| 午夜91福利影院| 久久精品人人爽人人爽视色| 一区二区三区乱码不卡18| 色网站视频免费| 精品久久久精品久久久| 寂寞人妻少妇视频99o| 春色校园在线视频观看| 午夜福利在线观看免费完整高清在| 中文字幕精品免费在线观看视频 | 亚洲欧美一区二区三区国产| 国产精品国产三级专区第一集| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 丰满饥渴人妻一区二区三| 最新的欧美精品一区二区| 欧美成人午夜精品| 日本色播在线视频| 18禁观看日本| 我的女老师完整版在线观看| 国产xxxxx性猛交| 精品午夜福利在线看| 亚洲三级黄色毛片| 伦理电影免费视频| 国产伦理片在线播放av一区| 午夜免费鲁丝| 久久久久久久久久久免费av| 一本久久精品| 成人毛片a级毛片在线播放| 下体分泌物呈黄色| 亚洲av中文av极速乱| 97精品久久久久久久久久精品| 18在线观看网站| 建设人人有责人人尽责人人享有的| 2022亚洲国产成人精品| 国产精品一二三区在线看| 欧美日韩成人在线一区二区| 欧美日韩亚洲高清精品| 少妇人妻 视频| 久久精品国产亚洲av天美| 极品人妻少妇av视频| 乱人伦中国视频| 中国三级夫妇交换| 午夜免费观看性视频| 亚洲欧洲精品一区二区精品久久久 | 日韩成人av中文字幕在线观看| 久久鲁丝午夜福利片| 91精品伊人久久大香线蕉| 成人亚洲欧美一区二区av| 日韩免费高清中文字幕av| 国产福利在线免费观看视频| 欧美xxxx性猛交bbbb| 一区二区三区精品91| 黄色 视频免费看| 少妇被粗大猛烈的视频| 国产亚洲av片在线观看秒播厂| 一本—道久久a久久精品蜜桃钙片| 国产色婷婷99| 青春草亚洲视频在线观看| 亚洲精品美女久久久久99蜜臀 | 精品福利永久在线观看| 99热6这里只有精品| 亚洲美女黄色视频免费看| 国产成人欧美| 精品久久久久久电影网| 亚洲国产看品久久| 两个人看的免费小视频| 久久精品人人爽人人爽视色| 亚洲av免费高清在线观看| 纵有疾风起免费观看全集完整版| 亚洲精品成人av观看孕妇| 午夜av观看不卡| 亚洲美女搞黄在线观看| av又黄又爽大尺度在线免费看| 亚洲精品久久成人aⅴ小说| 中国国产av一级| 久久久久国产网址| 高清欧美精品videossex| 国产成人精品无人区| 久久人人爽人人片av| av免费观看日本| 亚洲国产精品国产精品| 有码 亚洲区| kizo精华| 男男h啪啪无遮挡| 九九在线视频观看精品| a级毛片黄视频| 国产极品天堂在线| 男人舔女人的私密视频| 一本色道久久久久久精品综合| 欧美97在线视频| av天堂久久9| 岛国毛片在线播放| 午夜福利影视在线免费观看| 亚洲精品乱码久久久久久按摩| 天天影视国产精品| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| 亚洲成人手机| 交换朋友夫妻互换小说| 成人二区视频| 色吧在线观看| 免费看不卡的av| 日韩成人av中文字幕在线观看| 国产精品一区二区在线观看99| 一级毛片我不卡| 久久久久久久久久久久大奶| 国产精品久久久久久久久免| 一本—道久久a久久精品蜜桃钙片| 两个人免费观看高清视频| 中文字幕最新亚洲高清| videosex国产| 色婷婷久久久亚洲欧美| 国产精品无大码| 国产成人精品婷婷| 国产又色又爽无遮挡免| 久久韩国三级中文字幕| 边亲边吃奶的免费视频| 久热久热在线精品观看| 免费不卡的大黄色大毛片视频在线观看| 午夜福利乱码中文字幕| 一区二区三区四区激情视频| 寂寞人妻少妇视频99o| 黄网站色视频无遮挡免费观看| 成人国产麻豆网| 精品国产一区二区三区四区第35| 丰满迷人的少妇在线观看| 国产成人精品在线电影| 亚洲国产色片| av福利片在线| 久久狼人影院| 亚洲 欧美一区二区三区| 日本免费在线观看一区| 久久99热6这里只有精品| 久久这里只有精品19| 在线观看免费视频网站a站| 亚洲精品久久午夜乱码| 国产成人精品久久久久久| 韩国av在线不卡| 精品一区二区三区四区五区乱码 | 中国三级夫妇交换| 中国美白少妇内射xxxbb| 美女内射精品一级片tv| 纵有疾风起免费观看全集完整版| 欧美 日韩 精品 国产| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 春色校园在线视频观看| 中文欧美无线码| 高清在线视频一区二区三区| 欧美国产精品一级二级三级| 看十八女毛片水多多多| 日韩 亚洲 欧美在线| 观看av在线不卡| 久久国产精品男人的天堂亚洲 | 中文字幕人妻丝袜制服| 天堂俺去俺来也www色官网| 亚洲综合色惰| 尾随美女入室| 高清欧美精品videossex| av播播在线观看一区| 亚洲伊人色综图| 国产亚洲最大av| 伦理电影大哥的女人| 丝袜脚勾引网站| 日韩伦理黄色片| 成年动漫av网址| 免费观看无遮挡的男女| 精品人妻偷拍中文字幕| 久久久久久久国产电影| 日韩av免费高清视频| www.av在线官网国产| 国产精品熟女久久久久浪| 波多野结衣一区麻豆| 国产xxxxx性猛交| 久久精品久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 飞空精品影院首页| 2018国产大陆天天弄谢| 久久鲁丝午夜福利片| 国产精品久久久久久精品电影小说| 国产一区二区在线观看日韩| 亚洲国产看品久久| 国产成人免费无遮挡视频| 一级片免费观看大全| 国产欧美日韩一区二区三区在线| 大片电影免费在线观看免费| 国产精品.久久久| 久久久久精品性色| 亚洲美女视频黄频| 97超碰精品成人国产| av电影中文网址| 两性夫妻黄色片 | av片东京热男人的天堂| 成人综合一区亚洲| 卡戴珊不雅视频在线播放| 中文天堂在线官网| 熟女人妻精品中文字幕| 亚洲精品aⅴ在线观看| 国产乱来视频区| 久久精品久久久久久噜噜老黄| 人妻 亚洲 视频| 丰满少妇做爰视频| 91aial.com中文字幕在线观看| 国产精品麻豆人妻色哟哟久久| 久久青草综合色| 满18在线观看网站| 九九爱精品视频在线观看| 夫妻午夜视频| 国产在线免费精品| 性色avwww在线观看| 成年美女黄网站色视频大全免费| 欧美3d第一页| 亚洲内射少妇av| 久久久久精品人妻al黑| 热re99久久国产66热| 寂寞人妻少妇视频99o| 亚洲国产欧美日韩在线播放| 大片电影免费在线观看免费| 波野结衣二区三区在线| 麻豆乱淫一区二区| 日本爱情动作片www.在线观看| 免费高清在线观看日韩| 久久人人爽人人片av| 色吧在线观看| 99国产精品免费福利视频| 一区二区三区乱码不卡18| 成人手机av| 国产亚洲av片在线观看秒播厂| 青青草视频在线视频观看| 永久网站在线| 纯流量卡能插随身wifi吗| 欧美少妇被猛烈插入视频| 亚洲经典国产精华液单| 日韩一区二区三区影片| 免费不卡的大黄色大毛片视频在线观看| 91精品三级在线观看| 国产老妇伦熟女老妇高清| 波多野结衣一区麻豆| 美女中出高潮动态图| 美女福利国产在线| 国产极品天堂在线| 老女人水多毛片| 人人妻人人添人人爽欧美一区卜| 国产免费福利视频在线观看| 久久久久久人人人人人| 多毛熟女@视频| 欧美最新免费一区二区三区| 热re99久久精品国产66热6| 国产精品三级大全| 香蕉丝袜av| 久久久久久人妻| 久久人人爽av亚洲精品天堂| 国产福利在线免费观看视频| 伊人久久国产一区二区| 看免费av毛片| 亚洲精品久久久久久婷婷小说| 国产色婷婷99|