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

    不同品種釀酒葡萄表皮微生物群落多樣性分析

    2017-04-06 05:46:24張世偉陳曦鐘其頂黃占斌孟鎮(zhèn)羅金學石玲白志輝
    生物技術通報 2017年3期
    關鍵詞:釀酒葡萄酒群落

    張世偉 陳曦 鐘其頂 黃占斌 孟鎮(zhèn) 羅金學 石玲 白志輝

    (1. 中國礦業(yè)大學(北京)化學與環(huán)境工程學院,北京 100083;2. 中國食品發(fā)酵工業(yè)研究院,北京 100015;3. 中國科學院生態(tài)環(huán)境研究中心,北京 100085;4. 中糧長城桑干酒莊(懷來)有限公司,張家口 075400)

    不同品種釀酒葡萄表皮微生物群落多樣性分析

    張世偉1,2,3陳曦3鐘其頂2黃占斌1孟鎮(zhèn)2羅金學3石玲4白志輝3

    (1. 中國礦業(yè)大學(北京)化學與環(huán)境工程學院,北京 100083;2. 中國食品發(fā)酵工業(yè)研究院,北京 100015;3. 中國科學院生態(tài)環(huán)境研究中心,北京 100085;4. 中糧長城桑干酒莊(懷來)有限公司,張家口 075400)

    葡萄酒釀造是多種微生物參與代謝的過程,其中的微生物在沒有外源接種的情況下通常認為源自釀酒葡萄本身。應用高通量測序技術,分析沙城地區(qū)不同品種釀酒葡萄表皮的微生物群落,旨在從源頭上了解葡萄酒釀造原料的品質,為研究釀酒葡萄微生物對釀酒品質的影響提供理論依據。結果顯示,無論細菌還是真菌豐富度最大的均為雷司令,細菌主要分布在變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)等8個門,優(yōu)勢細菌主要是歐文氏菌屬(Erwinia)、假單胞菌屬(Pseudomonas)、芽孢桿菌屬(Bacillus)等6個屬。真菌僅分布在子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)和接合菌門(Zygomycota)3個門,鏈格孢菌屬(Alternaria)、枝孢屬(Cladosporium)、莖點霉屬(Phoma)、鐮孢屬(Fusarium)等9個主要的屬。細菌中心OTU最相似菌種分別為節(jié)桿菌(Arthrobacter sp.)、假單胞菌(Pseudomonas sp.)等,真菌為枝孢菌(Cladosporium sp.)、莖點霉(Phoma sp.)、鏈格孢菌(Alternaria sp.)以及大量未知種屬的OTU。研究表明,葡萄品種是影響微生物群落的最重要的因素,推測釀酒葡萄上的一些微生物會對葡萄植株的健康、葡萄果實的品質以及葡萄酒釀造產生有益或有害的影響。

    釀酒葡萄;高通量測序;微生物群落;葡萄酒品質;細菌多樣性;真菌多樣性

    葡萄酒釀造過程的本質是從葡萄收集、破碎成汁、葡萄酒一次、二次發(fā)酵直到包裝、儲藏整個工藝中多種微生物參與代謝的過程。新鮮采摘釀酒葡萄上的微生物數量大約為107-108CFU/g(葡萄鮮重),這些微生物會對后續(xù)的葡萄酒發(fā)酵過程產生重要影響[1]。Leveau等[2]通過焦磷酸測序方法研究美國加利福尼亞Clarksburg地區(qū)的釀酒葡萄發(fā)現假單胞桿菌屬(Pseudomonas)是霞多麗葡萄主要的細菌類群。Bokulich等[3]通過Illumina Miseq高通量測序方法研究了美國加利福尼亞兩個不同地區(qū)霞多麗、赤霞珠、仙粉黛釀酒葡萄汁中的微生物群落,發(fā)現其中的葡糖桿菌屬(Gluconobacter)等在葡萄酒釀造過程中能夠促進酯類物質的生成,改善葡萄酒的感官,而乳桿菌屬(Lactorbacillus)則會導致葡萄酒產生異味。Martin等[4]通過T-RFLP技術研究了法國Lussac St. Emilion地區(qū)不同葡萄園中美樂葡萄上的附生細菌,發(fā)現其優(yōu)勢菌群為假單胞桿菌屬(Pseudomonas)和鞘脂單胞菌屬(Sphingomonas),其中的假單胞桿菌屬(Pseudomonas)是葡萄中普遍存在的細菌類群;同時該研究還發(fā)現釀酒葡萄園土壤中的微生物和葡萄上的微生物存在聯系,因此推測植物不同器官上的附生菌可能有著共同的來源。迄今為止,雖已有少量有關釀酒葡萄微生物的報道,但對其群落結構特點、對后續(xù)葡萄酒發(fā)酵的影響等仍然知之甚少。因此,探明不同品種釀酒葡萄微生物群落結構,從源頭上了解葡萄酒釀造原料的品質,對于保證葡萄酒產業(yè)健康發(fā)展具有重要意義。

    沙城地區(qū)多丘陵山地,土壤為褐土,光照充足,晝夜溫差大,干燥少雨,是釀酒葡萄種植最為理想的黃金地帶。目前,沙城地區(qū)種植著幾十個品種的釀酒葡萄,幾乎涵蓋了我國釀酒葡萄的全部品種,是我國優(yōu)質葡萄酒的主要產區(qū)之一。本研究應用Illumina平臺的Hiseq 2500高通量測序技術,分析了沙城地區(qū)不同品種釀酒葡萄表面的細菌、真菌群落組成、結構以及多樣性,旨在為研究不同品種釀酒葡萄微生物對釀酒品質的影響提供理論依據。

    1 材料與方法

    1.1 材料

    1.1.1 樣地概況 采樣樣地位于河北省懷來縣沙城鎮(zhèn)某釀酒葡萄園內,該地區(qū)位于北緯40°4'-40°35',東經115°16'-115°58',屬于溫帶大陸性季風氣候,雨熱同期,光照充足,日照率68%,年均氣溫9.1℃,白天太陽輻射強,夜間地面散熱快,晝夜溫差大,5-9月是農作物生長季節(jié),期間晝夜溫差可達11℃-12.4℃。樣地屬于規(guī)模化釀酒葡萄園,葡萄樹齡最短的為10年,最長的為37年。

    1.1.2 樣品采集 分別采集霞多麗(Chardonnay,Cha)、 雷 司 令(Riesling,Rie)、 黑 品 諾(Pinot nior,Pin)、寶石(Gem,Gem)、赤霞珠(Cabernet Sauvignon,Cab)、增芳德(Zinfandel,Zin)、龍眼(Longan,Lon)、西拉(Syrah,Syr)、美樂(Merlot,Mer)共9個品種的釀酒葡萄果實。釀酒葡萄樣品用經75%酒精消毒后的剪刀采集新鮮、飽滿、無病害的果實,置于無菌自封袋中,再將裝有葡萄樣品的自封袋置于低溫采樣箱中,帶回實驗室,-20℃冷凍備用。每個品種釀酒葡萄樣品采集均按照“五點取樣法”收集5個植株的果實,混合后用于后續(xù)實驗。部分品種釀酒葡萄發(fā)酵制成葡萄酒的理化指標如表1所示,均采用GB/15038-2006方法測定。其中龍眼是我國特有的葡萄品種,和其他釀酒葡萄品種相比,其含糖量較高,既可作為釀酒葡萄,釀造半甜型葡萄酒,也可作為鮮食葡萄。

    表1 不同品種釀酒葡萄發(fā)酵制成葡萄酒的理化指標

    1.2 方法

    1.2.1 總DNA提取 每個釀酒葡萄品種果實樣品稱取約30 g,加入100 mL PBS(0.1 mol/L,pH7.0)緩沖液,180 r/min渦旋震蕩30 min,超聲15 min,0.22 μm微孔濾膜(Millipore,USA)抽濾,濾膜剪碎放入滅菌離心管中,用DNA提取試劑盒(FastDNA SPIN Kit for Soil,MP)提取總DNA,用核酸蛋白測定儀(Nano Drop,Thermo fisher)檢測其濃度和純度。

    1.2.2 PCR擴增和測序 采用細菌16S rRNA基因V5-V7區(qū)引物799F(5'-AACGGATTAGATACCG-3')和1193R(5'-ACGTCATCCCCACCTTCC-3')[5], 以及真菌ITS1區(qū)引物(5'-CTTGGTCATTTAGAGGAAGTAA-3')和(5'-GCTGCGTTCTTCATCGAT GC-3')[4]進行PCR擴增。50 μL反應體系如下:DNA模板1 μL,正反向引物(10 μmol/L)各2 μL,dNTPs(2.5 μmol/L)4 μL,10×Pyrobest Buffer 5 μL,Pyrobest DNA Polymerase(2.5 U/μL)0.3 μL,ddH2O 35.7 μL。擴增程序如下:94℃預變性5 min;94℃變性30 s,55℃退火 30 s,72℃延伸1 min,28個循環(huán);最后72℃延伸10 min。每個樣品3個重復,將同一樣品3個重復的PCR產物混勻后用1.2%的瓊脂糖凝膠電泳進行檢測,再用AxyPrep DNA 凝膠回收試劑盒(AXYGEN,USA)切膠回收。將回收的PCR產物用Qubit(Life Technology,USA)進行定量化和均一化,用Illumina建庫試劑盒進行文庫構建,并用Agilent 2100 (Agilent,USA)檢測文庫片段分布,將構建好的合格文庫送至北京新科開源基因有限公司進行測序(Hiseq 2500,PE250)。

    1.2.3 數據分析

    1.2.3.1 數據質控及OTU劃分 將下機原始數據進行質控,刪除含有連續(xù)>5 bp的Q值<15的低質量堿基及含有接頭序列的測序數據。然后,根據barcode標簽提取每個樣本序列數據,用軟件Mothur(version 1.36.1)[6]進行數據清洗及操作分類單元(Operational taxonomic unit,OTU)劃分,OTU序列相似度設為97%,統(tǒng)計每個樣品含有OTU的數目以及每個OTU含有序列的條數,以每個OTU中重復條數最多的序列作為每個OTU的代表序列[7]。

    1.2.3.2 豐富度和多樣性分析 選擇在序列97%相似度閾值下繪制稀釋曲線(rarefaction curve),稀釋曲線反映了樣品的取樣深度,可以用來評價測序數量是否足以覆蓋所有類群。計算豐富度指數Chao 1和ACE,二者反映樣品的豐富度,數值越高說明群落物種的豐富度越高。計算多樣性指數 Shannon和Simpson,二者反映樣品的多樣性程度,Shannon數值越高說明群落物種的多樣性越高,Simpson則相反,數值越高說明群落物種的多樣性越低[8-10]。

    1.2.3.3 物種信息注釋 為得到每個OTU對應的物種分類信息,采用RDP算法對97%閾值的OTU代表序列進行物種注釋。用數據庫Silva(Release115 http://www.arb-silva.de) 和Unite(https://unite.ut.ee/)分別對細菌和真菌進行比對,在門(Phylum)和屬(Genus)的水平上統(tǒng)計每個樣品的群落組成,為了得到主要菌群信息,規(guī)定細菌門的水平上選取相對豐度大于0.5%的門,屬的水平上選取相對豐度大于1%的屬;真菌門的水平上菌群分類較少,不設定閾值,屬的水平上選取相對豐度大于1%的屬。提取序列數占序列總數1%以上的OTU為中心OTU,在NCBI上進行Blast序列比對。其他數據分析及統(tǒng)計均采用軟件R 3.0和Origin 9.0。

    2 結果

    2.1 有效序列及OTU統(tǒng)計

    對9個釀酒葡萄品種果實表面的細菌16S rRNA基因V5-V7區(qū)以及真菌ITS1區(qū)進行PCR擴增和測序,過濾掉原始數據中的低和霞多麗的真菌測序量不足,舍棄,其余 8個細菌樣品共得到154 612條有效序列,在質量序列,去除線粒體DNA、葉綠體DNA以及嵌合體等,其中葡萄品種西拉的細菌測序量97%相似度下劃分OTU,共得到1 897個OTU;8個真菌樣品共得到714 707條有效序列,在97%相似度下將其聚為不同的OTU,共得到1 313個OTU。

    2.2 豐富度和多樣性分析

    如圖1所示,相似度在97%條件下不同品種釀酒葡萄表面細菌16S rRNA基因V5-V7區(qū)和真菌ITS1區(qū)的稀釋曲線,當細菌OTU數量達到600個左右,真菌OTU數量達到300個左右時,曲線開始趨于平緩,說明本研究中樣品的取樣深度足夠,測序數量足以覆蓋所有類群。

    圖1 97%相似度條件下的稀釋曲線

    如表2所示,不同品種釀酒葡萄表面細菌16S rRNA基因V5-V7區(qū)樣本的測序深度指數coverage均較高(最低的為龍眼94.87%),說明此次測序能夠覆蓋樣本中絕大多數的微生物信息,可以代表樣本的真實情況。不同品種釀酒葡萄表面細菌的豐富度指數Chao1和ACE最大的為雷司令,最小的為赤霞珠,說明這8種釀酒葡萄表面雷司令細菌群落最豐富,赤霞珠則最稀少。細菌的多樣性指數Shannon最大的為霞多麗,最小的為赤霞珠,Simpson指數則基本相反,說明這8種釀酒葡萄表面霞多麗細菌群落多樣性最高,赤霞珠最低。

    如表3所示,不同品種釀酒葡萄表面真菌ITS1區(qū)測序深度指數coverage更高,即使是最低的龍眼也可達99.23%,表明此次測序數據量可以代表樣本的真實情況。不同品種釀酒葡萄表面真菌的豐富度指數Chao 1和ACE最大的均為雷司令,Chao1指數最小的為龍眼,ACE指數最小的為寶石,龍眼也較小,Chao1和ACE基本一致,說明這8種釀酒葡萄表面雷司令真菌群落最豐富,龍眼和寶石則最稀少。真菌的多樣性指數Shannon最大的為增芳德,最小的為美樂,Simpson指數則基本相反,說明這8種釀酒葡萄表面增芳德真菌群落多樣性最高,美樂最低。

    表2 不同品種釀酒葡萄表面細菌16S rRNA基因V5-V7區(qū)測序深度、豐富度及多樣性指數

    表3 不同品種釀酒葡萄表面真菌ITS1區(qū)測序深度、豐富度及多樣性指數

    2.3 細菌和真菌類群分析

    在門的分類水平上,不同品種釀酒葡萄表面細菌除未被分類(unclassfied)類群外,主要分布在8個已知的門(圖2-A):變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)、擬 桿 菌 門(Bacteroidetes)、 芽 單 胞 菌 門(Gemmatimonadetes)、醋桿菌門(Acidobacteria)、綠彎菌門(Chloroflexi)和Thermi。在這8個門的細菌中,占比最高的是變形菌門(Proteobacteria),可占霞多麗總細菌的75.30%,占比最低的黑品諾也可達31.30%。在屬的分類水平上,不同品種釀酒葡萄表面細菌主要分布在6個屬(圖2-B):歐文氏菌屬(Erwinia)、假單胞菌屬(Pseudomonas)、芽孢桿菌屬(Bacillus)、丙酸桿菌屬(Propionibacterium)、Kaistobacter和小雙孢菌屬(Microbispora)。在這6個屬的細菌中,占比最高的是歐文氏菌屬(Erwinia),可占赤霞珠總細菌的46.20%,占比最低的黑品諾則僅為4.20%。

    如圖2-C所示,不同品種釀酒葡萄表面真菌多樣性較少,僅分布在3個門:子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)和接合菌門(Zygomycota)。在這3個門的真菌中,占比最高的是子囊菌門(Ascomycota),可占雷司令總真菌的92.50%,占比最低的寶石也可達89.38%。占比最低的是接合菌門(Zygomycota),僅占龍眼總真菌的1.43%,增芳德則低達0.63%。如圖2-D所示,不同品種釀酒葡萄表面真菌主要分布在9個屬(排序不分先后):鏈格孢菌屬(Alternaria)、枝孢屬(Cladosporium)、莖點霉屬(Phoma)、鐮孢屬(Fusarium)、耐冷酵母屬(Guehomyces)、孢霉屬(Mortierella)、隱球菌屬(Cryptococcus)、黑附球菌屬(Epicoccum)和漆斑霉屬(Myrothecium)。在這9個屬的真菌中,占比最高的是鏈格孢菌屬(Alternaria),可占黑品諾總真菌的46.24%,占比最低的霞多麗也可達20.57%。

    圖2 不同品種釀酒葡萄表面微生物在門和屬水平上的比較

    2.4 中心OTU及最相似序列分析

    如表4所示,對于不同品種釀酒葡萄表面細菌提取序列數占序列總數1%以上的OTU為中心OTU,得到16個相對豐度較高的OTU。細菌中心OTU所含序列數占序列總數比例最高的是赤霞珠(69.02%),最低的是霞多麗(43.13%)。將這些OTU的代表序列在NCBI上進行Blast序列比對得到的最相似菌種分別為節(jié)桿菌(Arthrobacter sp., OTU4、OTU7和OTU9)、假單胞菌(Pseudomonas sp.,OTU1和OTU10)、芽孢桿菌(Bacillus sp.,OTU2)、動性球菌(Planococcus sp.,OTU5)、泛生菌(Pantoea sp.,OTU6)、短小桿菌(Curtobacterium sp.,OTU8)、纖維菌(Cellulomonas sp.,OTU11)、柄桿菌(Caulobacter sp.,OTU12)、迪茨氏菌(Dietzia sp.,OTU-21)、Uncultured malikia sp.(OTU22)、甲基桿菌(Methylobacterium sp.,OTU37)、根瘤菌(Mesorhizobium sp.,OTU47)和未知種屬的放線菌(OTU23)。

    表4 不同品種釀酒葡萄表面細菌中心OTU及最相似序列

    如表5所示,對于不同品種釀酒葡萄表面真菌提取序列數占序列總數1%以上的OTU為中心OTU,得到21個相對豐度較高的真菌OTU,真菌中心OTU所含序列數占序列總數的比例均在80%以上。將這些OTU的代表序列在NCBI上進行Blast序列比對得到最相似菌種分別為枝孢菌(Cladosporium sp.,OTU2、OTU5、OTU6和OTU9)、莖點霉(Phoma sp.,OTU3、OTU7和OTU37)、鏈格孢菌(Alternaria sp.,OTU1和OTU4)、Gibberia sp.(OTU8)、 鐮 孢菌(Fusarium sp.,OTU14)、黑附球菌(Epicoccum sp.,OTU15)、隱球菌(Cryptococcus sp.,OTU17)、小 穴 殼 菌(Dothioraceae sp.,OTU33)、Plesporales sp.(OTU55)、Kwoniella sp.(OTU82) 和 未 知 種屬 的OTU(OTU10、OTU11、OTU16、OTU18和OTU23)。

    表5 不同品種釀酒葡萄表面真菌中心OTU及最相似序列

    3 討論

    3.1 高通量測序技術是分析微生物群落的重要工具

    隨著分子生物學技術的發(fā)展,新一代高通量測序技術憑借其通量高、速度快、成本低等特點為研究微生物多樣性提供了新的技術平臺。本研究采用Illumina平臺的Hiseq 2500高通量測序技術對不同品種釀酒葡萄表面細菌16S rRNA基因V5-V7區(qū)和真菌ITS1區(qū)進行了擴增和測序,結果顯示8個細菌樣品共得到154 612條有效序列,最少的為龍眼,有5 191條,最多的為霞多麗,有24 995條;8個真菌樣品共得到 714 707條有效序列,最少的為龍眼,有17 466條,最多的為梅璐哲,有152 569條。采用細菌16S rRNA基因V5-V7區(qū)引物(799F和1193R)是為了在擴增過程中避免線粒體DNA和葉綠體DNA的干擾,這樣雖然避免了植物來源物種信息的干擾,但是也無法捕捉到藍細菌等的物種信息;而且,實際上,引物799F僅能和RDP數據庫中62%的序列匹配[11],因此,本研究中細菌的有效序列條數(154 612條)較真菌(714 707條)低。在97%相似度下將其聚為不同的OTU,細菌共得到1 897個OTU,真菌共得到1 313個OTU。各樣品的覆蓋度(coverage)結果顯示在相似度為97%的條件下,OTU涵蓋了釀酒葡萄表面94%以上的細菌和真菌。以上結果說明此次測序能夠覆蓋樣本中絕大多數的微生物信息,可以代表樣本的真實情況。因此,高通量測序技術是分析釀酒葡萄微生物群落的有力工具。

    高通量測序技術較傳統(tǒng)的分子生物學手段(如聚丙烯酰胺凝膠電泳(Denaturing gradient gel electrophoresis,DGGE)、克隆文庫、末端限制性片段長度多態(tài)性(Terminal restriction fragment length polymorphism,T-RFLP)和熒光原位雜交(Fluorescence in situ hybridization,FISH)等以及磷脂脂肪酸(Phospholipid fatty acid,PLFA)等技術,具有明顯的優(yōu)勢。例如,Redford等[12-14]通過DGGE技術研究棉白楊葉際微生物,結果發(fā)現每個樣品的可視條帶僅為幾十條,很難反映葉際微生物群落的全貌,而用高通量測序技術每個樣品能夠得到上萬條序列,其檢測靈敏度遠遠超過DGGE,可以得到更為全面的微生物群落信息。T-RFLP技術可通過和已有數據庫進行比對,獲得所擴增片段的系統(tǒng)進化信息,如Berg等[15]利用此技術研究了馬鈴薯抵抗植物病原真菌的附生細菌的群落結構,但是該技術對微生物種屬的定性判斷存在一定誤差。PLFA技術可通過不同微生物含有不同種類和數量的PLFAs對微生物進行分析,如Zhang等[16,17]利用此技術研究了氯氰菊酯殺蟲劑對辣椒和黃瓜葉際微生物群落的影響,但是該技術僅能檢測到葉際中活體微生物群落的動態(tài)變化??梢姡咄繙y序技術較傳統(tǒng)的非培養(yǎng)方法可以得到更為全面、準確、豐富的微生物群落信息。

    高通量測序技術因其突出的優(yōu)點已經在蔬菜、水果等植物表面微生物群落研究中得以應用。例如,2010年,Lopez-Velasco等[18]將焦磷酸測序技術應用于菠菜葉際細菌16S rRNA 基因多樣性的研究,在新鮮菠菜葉際發(fā)現了8 800個獨有的序列,其中75%是以前從未報道過的。2012年,Bokulich等[19]通過焦磷酸測序技術研究了貴腐葡萄酒發(fā)酵過程中細菌16S rRNA基因V4-V5區(qū)的動態(tài)變化,從發(fā)酵結束的酒泥中分離出α-變形菌(Alphaproteobacteria)兩個種的細菌,揭示了低豐度微生物群落的多樣性,這些結果表明貴腐葡萄酒較之前的認識更為復雜,同時為監(jiān)測葡萄酒發(fā)酵過程的動態(tài)變化奠定了基礎。2013年,Leff等[20]通過焦磷酸測序技術對不同品種新鮮水果和蔬菜表面細菌16S rRNA基因V5-V6區(qū)進行了分析,結果發(fā)現蘋果、桃、葡萄以及蘑菇表面主要的細菌類群為放線菌門(Actinobacteria)、芽孢桿菌門(Bacillus)、厚壁菌門(Firmicutes)和變形菌門(Proteobacteria)。如今,高通量測序技術仍然在不斷改進,必將推動著微生物生態(tài)學研究更加迅猛的發(fā)展。

    3.2 不同品種釀酒葡萄優(yōu)勢菌群存在差異

    本研究發(fā)現釀酒葡萄表面細菌在門的水平上最優(yōu)勢的類群是變形菌門(Proteobacteria),放線菌門(Actinobacteria)和厚壁菌門(Firmicutes)。Leveau等[2]通過焦磷酸測序技術發(fā)現霞多麗葡萄上53.5%的細菌為變形菌門(Proteobacteria),其他主要有厚壁菌門(Firmicutes,15.1%),擬桿菌門(Bacteroidetes,10.1%)和放線菌門(Actinobacteria,8.0%)。Martin等[4]通過T-RFLP技術發(fā)現美樂葡萄上的優(yōu)勢細菌也主要歸屬于變形菌門(Proteobacteria),本研究結果與上述結果基本一致。另外,本研究發(fā)現釀酒葡萄表面細菌在屬的水平上最優(yōu)勢的類群是歐文氏菌屬(Erwinia)、假單胞菌屬(Pseudomonas)和芽孢桿菌屬(Bacillus)。Leveau等[2]發(fā)現霞多麗葡萄上的優(yōu)勢菌屬為假單胞菌屬(Pseudomonas)。Martin等發(fā)現美樂葡萄上的優(yōu)勢菌屬為假單胞菌屬(Pseudomonas)和鞘脂單胞菌屬(Sphingomona)。本研究結果與這些結果具有相似之處,只是本研究發(fā)現鞘脂單胞菌屬(Sphingomona)并不是沙城地區(qū)美樂葡萄的優(yōu)勢菌屬。在真菌方面,有關葡萄品種和葡萄表面真菌之間關系的研究較少[20-23],Bokulish等[4,19]研究發(fā)現增芳德葡萄表面上的優(yōu)勢真菌菌屬為C. zemplinina,C. zemplinina能夠增強葡萄酒的感官品質,但是本研究發(fā)現增芳德主要的真菌菌屬為鏈格孢菌屬(Alternaria),這與前人研究結果差異較大。造成本研究結果與上述文獻報道結果不一致的原因可能是多方面的,如葡萄產地、環(huán)境以及采樣時間等,Martin等[24]的采樣時間是在秋季葡萄成熟期,本研究采樣的時間是在七月葡萄生長旺盛期,此時正處于葡萄漿果著色的轉色期,這個時期,漿果已不再膨大,果皮中的葉綠素大量分解,紅葡萄品種果皮開始積累色素,基本轉變?yōu)榧t色和深藍色,而白葡萄品種透明,漿果含糖量直線上升,這個時期基本是大部分葡萄品種含酸量下降的轉折點。

    本研究還發(fā)現釀酒葡萄表面微生物的優(yōu)勢類群因葡萄品種不同存在較大差異。無論是細菌還是真菌在門的水平上差異相對較小,尤其是真菌,3個優(yōu)勢門的真菌可占雷司令總真菌的92.50%,占比最低的寶石也可達89.38%,但在屬的水平上差異較大。在細菌屬的水平上,霞多麗上的歐文氏菌屬(Erwinia)可占細菌總數的46.20%,而雷司令僅為4.30%;雷司令上假單胞菌屬(Pseudomonas)為10.17%,黑品諾僅為3.03%;黑品諾上的芽孢桿菌屬(Bacillus)為19.96%,赤霞珠僅為1.43%;西拉上的丙酸桿菌屬(Propionibacterium)為16.39%,赤霞珠僅為0.03%。在真菌屬的水平上,龍眼上的鏈格孢菌屬(Alternaria)可占真菌總數的39.40%,最少的霞多麗也占到20.57%;龍眼上的枝孢屬(Cladosporium)為24.14%,而西拉僅為6.74%;赤霞珠上的莖點霉屬(Phoma)為30.62%,而雷司令僅為0.64%;增芳德中的隱球菌屬(Cryptococcus)為1.17%,霞多麗僅為0.29%。由此可見,釀酒葡萄表皮微生物群落因葡萄品種不同差異較大,葡萄品種是影響微生物群落的最重要的因素。

    3.3 微生物群落及其對葡萄酒釀造的影響

    研究還發(fā)現,細菌中心OTU結果顯示釀酒葡萄表面假單胞菌(Pseudomonas sp.)和芽孢桿菌(Bacillus sp.)含量較多,這兩種細菌是自然生態(tài)系統(tǒng)中廣泛存在的類群,其中的假單胞菌能夠產生胞外多糖,有利于微生物膜(biofilm)的形成,對于微生物在葡萄表面定殖有著重要的作用。細菌中心OTU結果還顯示葡萄表面還有較多短小桿菌(Curtobacterium sp.)和根瘤菌(Mesorhizobium sp.),有研究表明短小桿菌(Curtobacterium sp.)分泌物對葡萄酒品質無明顯影響[24],而根瘤菌(Mesorhizobium sp.)具有固氮的獨特作用,能夠促進葡萄樹的生長[25]。

    真菌中心OTU結果顯示釀酒葡萄表面含有較多鏈格孢菌(Alternaria sp.),有研究表明其可產生鏈格孢醇、鏈格孢醇單甲醚、鏈孢菌酮酸等70多種有毒的真菌代謝物從而引發(fā)葡萄黑斑病,人及動物攝入這些真菌毒素可導致急性或慢性中毒,有些還有致癌、致畸、致突變等的作用[26,27]。鏈格孢菌(Alternaria sp.)在其他水果如蘋果、西紅柿、藍莓以及谷物微生物研究中也有報道[28]。真菌中心OTU結果還顯示葡萄表面還有較多枝孢菌(Cladosporium sp.)、莖點霉(Phoma sp.)和隱球菌(Cryptococcus sp.),枝孢菌(Cladosporium sp.)是環(huán)境中廣泛存在的絲狀真菌,莖點霉(Phoma sp.)和隱球菌(Cryptococcus sp.)通常主要來源于葡萄園。枝孢菌(Cladosporium sp.)和莖點霉(Phoma sp.)可使葡萄腐爛而產生毒素,會造成葡萄中化學成分改變并且在葡萄酒發(fā)酵過程中可影響釀酒酵母的生長繁殖。隱球菌(Cryptococcus sp.)可產生果膠酶,提高葡萄酒釀造過程中的出汁率[29-31]。葡萄酒發(fā)酵過程中的微生物在沒有外源接種的情況下通常認為源自釀酒葡萄本身[32]。生長在葡萄上的霉菌產生的多種代謝產物(如真菌毒素赭曲霉素A)會干擾葡萄的微生物生態(tài)環(huán)境,從而影響葡萄酒發(fā)酵過程中釀酒酵母的生長,改變酒的風味[33,34]。真菌還能夠在葡萄表面創(chuàng)造一種有利于醋酸菌生長的環(huán)境,產生酸性物質(如醋酸等),從而改變葡萄汁和葡萄酒的pH,并且產生酸味,對葡萄酒的感官、生化和微生物的穩(wěn)定性有直接影響[35-37],甚至會導致酒的腐敗變質[38]。

    4 結論

    不同品種釀酒葡萄表皮微生物群落結構差異明顯,微生物豐度最大的為雷司令,細菌主要分布在變形菌門(Proteobacteria)、放線菌門(Actinobacteria)、厚壁菌門(Firmicutes)等8個門,優(yōu)勢細菌主要分布在歐文氏菌屬(Erwinia)、假單胞菌屬(Pseudomonas)、芽孢桿菌屬(Bacillus)等6個屬。真菌僅分布在子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)和接合菌門(Zygomycota)3個門,鏈格孢菌屬(Alternaria)、枝孢屬(Cladosporium)、莖點霉屬(Phoma)、鐮孢屬(Fusarium)等9個屬。

    [1]孫煒寧, 張巧格, 李興興, 等. 葡萄酒釀造過程中微生物多樣性的研究現狀[J]. 食品研究與開發(fā), 2014, 35(18):365-368.

    [2]Leveau JHJ, Tech JJ. Grapevine microbiomics:bacterial diversity on grape leaves and berries revealed by high-throughput sequence analysis of 16S rRNA amplicons[C]. International Symposiumon Biological Control of Postharvest Diseases:Challenges and Opportunities 905, 2010:31-42.

    [3]Bokulich NA, Thorngate JH, Richardson PM, et al. Microbial biogeography of wine grapes is conditioned by cultivar, vintage, and climate[J]. Proceedings of the National Academy of Sciences, 2014, 111(1):E139-E148.

    [4]Martins G, Lauga B, Miot-Sertier C, et al. Characterization of epiphytic bacterial communities from grapes, leaves, bark and soil of grapevine plants grown, and their relations[J]. PLoS One, 2013, 8(8):e73013.

    [5]Chelius MK, Triplett EW. The diversity of archaea and bacteria in association with the Roots of Zea mays L. [J]. Microbial Ecology, 2001, 41(3):252-263.

    [6]Schloss PD, Westcott SL, Ryabin T, et al. Introducing mothur:open-source, platform-independent, community-supported software for describing and comparing microbial communities[J]. Applied and Environmental Microbiology, 2009, 75(23):7537-7541.

    [7]Zarraonaindia I, Owens SM, Weisenhorn P, et al. The soil microbiome influences grapevine-associated microbiota[J]. MBio, 2015, 6(2):e02527-14.

    [8]Chao A. Nonparametric estimation of the number of classes in a population[J]. Scandinavian Journal of Statistics, 1984, 11(4):265-270.

    [9]郭正剛, 王根緒, 沈禹穎, 等. 青藏高原北部多年凍土區(qū)草地植物多樣性[J]. 生態(tài)學報, 2004, 24(1):149-155.

    [10]Gupta C, Dubey R, Maheshwari D. Plant growth enhancement and suppression of Macrophomina phaseolina causing charcoal rot of peanut by fluorescent Pseudomonas[J]. Biology and Fertility of Soils, 2002, 35(6):399-405.

    [11]Redford AJ, Bowers RM, Knight R, et al. The ecology of the phyllosphere:geographic and phylogenetic variability in the distribution of bacteria on tree leaves[J]. Environmental Microbiology, 2010, 12(11):2885-2893.

    [12]Bodenhausen N, Horton MW, Bergelson J. Bacterial communities associated with the leaves and the roots of Arabidopsis thaliana[J]. PLoS One, 2013, 8(2):e56329.

    [13]Redford AJ, Fierer N. Bacterial succession on the leaf surface:a novel system for studying successional dynamics[J]. Microbial Ecology, 2009, 58(1):189-198.

    [14]夏圍圍, 賈仲君. 高通量測序和DGGE分析土壤微生物群落的技術評價[J]. 微生物學報, 2014, 54(12):1489-1499.

    [15]Berg G, Krechel A, Ditz M, et al. Endophytic and ectophytic potato-associated bacterial communities differ in structure and antagonistic function against plant pathogenic fungi[J]. FEMS Microbiology Ecology, 2005, 51(2):215-229.

    [16]Zhang B, Bai Z, Hoefel D, et al. The impacts of cypermethrin pesticide application on the non-target microbial community of the pepper plant phyllosphere[J]. Science of the Total Environment, 2009, 407(6):1915-1922.

    [17]Zhang B, Zhang H, Bo Jin, et al. Effect of cypermethrin insecticide on the microbial community in cucumber phyllosphere[J]. Journal of Environmental Sciences, 2008, 20(11):1356-1362.

    [18]Lopez-Velasco G, Welbaum GE, Boyer RR, et al. Changes in spinach phylloepiphytic bacteria communities following minimal processing and refrigerated storage described using pyrosequencing of 16S rRNA amplicons[J]. Journal of Applied Microbiology, 2011, 110(5):1203-1214.

    [19]Bokulich NA, Joseph CML, Allen G, et al. Next-generation sequencing reveals significant bacterial diversity of botrytized wine[J]. PLoS One, 2012, 7(5):e36357.

    [20]Leff JW, Fierer N. Bacterial communities associated with the surfaces of fresh fruits and vegetables[J]. PLoS One, 2013, 8(3):e59310.

    [21]Sabate J, Cano J, Esteve-Zarzoso B, et al. Isolation and identification of yeasts associated with vineyard and winery by RFLP analysis of ribosomal genes and mitochondrial DNA[J]. Microbiological Research, 2002, 157(4):267-274.

    [22]Cordero-Bueso G, Arroyo T, Serrano A, et al. Influence of the farming system and vine variety on yeast communities associated with grape berries[J]. International Journal of Food Microbiology, 2011, 145(1):132-139.

    [23]Tofalo R, Schirone M, Torriani S, et al. Diversity of Candida zemplinina strains from grapes and Italian wines[J]. Food Microbiology, 2012, 29(1):18-26.

    [24]Martins G, Miot-Sertier C, Lauga B, et al. Grape berry bacterial microbiota:impact of the ripening process and the farming system[J]. International Journal of Food Microbiology, 2012, 158(2):93-100.

    [25]姜怡, 楊穎, 陳華紅, 等. 植物內生菌資源[J]. 微生物學通報, 2005, 32(6):146-147.

    [26]da Motta S, Valente Soares L M. Survey of Brazilian tomato products for alternariol, alternariol monomethyl ether, tenuazonic acid and cyclopiazonic acid[J]. Food Additives & Contaminants, 2001, 18(7):630-634.

    [27] 李鳳琴. 鏈格孢毒素及其食品衛(wèi)生問題[J]. 中國食品衛(wèi)生雜志, 2001, 13(6):45-49.

    [28]Chulze S, Torres A, Dalcero A, et al. Production of alternarioi and alternariol monomethyl ether in natural substrates in comparison with semisynthetic culture medium[J]. Mycotoxin Research, 1994, 10(2):79-84.

    [29]Pretorius IS, Van der Westhuizen TJ, Augustyn OPH. Yeast biodiversity in vineyards and wineries and its importance to the South African wine industry. A review[J]. South African Journal of Enology and Viticulture, 1999, 20(2):61-68.

    [30]Barata A, Malfeito-Ferreira M, Loureiro V. The microbial ecology of wine grape berries[J]. International Journal of Food Microbiology, 2012, 153(3):243-259.

    [31]李華. 現代葡萄酒工藝學[M]. 西安:陜西人民出版社, 2000.

    [32]Bokulich NA, Ohta M, Richardson PM, et al. Monitoring seasonal changes in winery-resident microbiota[J]. PLoS One, 2013, 8(6):e66437.

    [33]Stander MA, Steyn PS. Survey of ochratoxin A in South African wines[J]. S Afr J Enol. Vitic, 2002, 23(1):9-13.

    [34]Jackson RS. Wine science:principles and applications[M]. Academic Press, 2008.

    [35] Boulton RB, Singleton VL, Bisson LF, et al. Principles and practices of winemaking[M]. Springer Science & Business Media, 2013.

    [36] Swiegers JH, Bartowsky EJ, Henschke PA, et al. Yeast and bacterial modulation of wine aroma and flavour[J]. Australian Journal of Grape and Wine Research, 2005, 11(2):139-173.

    [37]Innerebner G, Knief C, Vorholt JA. Protection of Arabidopsis thaliana against leaf-pathogenic Pseudomonas syringae by Sphingomonas strains in a controlled model system[J]. Applied and Environmental Microbiology, 2011, 77(10):3202-3210.

    [38]Fowles GWA. Acids in grapes and wines:A review[J]. Journal of Wine Research, 1992, 3(1):25-41.

    (責任編輯 馬鑫)

    Microbial Communities on the Wine Grape Surfaces of Different Cultivars

    ZHANG Shi-wei1,2,3CHEN Xi3ZHONG Qi-ding2HUANG Zhan-bin1MENG Zhen2LUO Jin-xue3SHI Ling4BAI Zhi-hui3
    (1. School of Chemical & Environmental Engineering,China University of Mining &Technology(Beijing),Beijing 100083;2. China National Research Institute of Food & Fermentation Industries,Beijing 100015;3. Research Center for Eco-Environmental Sciences,Chinese Academy of Sciences,Beijing 100085;4. Chateau SunGod GreatWall,Zhangjiakou 075400)

    Brewing grape wine is a process involved by many microorganisms that generally come from the grape berries if without exogenous inoculation. The high-throughput sequencing technology was used to analyze the microbial communities on the surface of different cultivars of wine grapes in Shacheng,China,which aims to understand the quality of the wine grapes,and provide the theoretical evidences for investigating the influences of microorganisms on the quality of brewing. The results showed that the richest microbial communities among these wine grapes were Riesling no matter of bacteria or fungi. The bacteria mainly were distributed in 8 phyla of Proteobacteria,Actinobacteria,Firmicutes,etc.,and 6 dominant genera of Erwinia,Pseudomonas,Bacillus,and so on. The 3 phyla of the fungi were Ascomycota,Basidiomycota,and Zygomycota;and the 9 genera of fungi were Alternaria,Cladosporium,Phoma,Fusarium,etc. The most similar species of the bacterial core OTUs were Arthrobacter sp. and Pseudomonas sp.,as well as Cladosporium sp.,Phoma sp.,Alternaria sp. and a lot of unknown ones for the fungal core OTUs. This study illustrates that the cultivar is the most important factor affecting the microbial communities. And from the results,it is inferred that the microorganisms on the surfaces of the wine grapes may be beneficial or adverse to the health of the grape vine,the quality of the grape berries and the process of wine brewing.

    wine grape;high-throughput sequencing;microbial community;grape wine quality;bacterial diversity;fungal diversity

    10.13560/j.cnki.biotech.bull.1985.2017.03.019

    2016-06-20

    國家自然科學基金項目(31570494),國際科技合作專項項目(2015FA31720),北京市科技新星計劃(xx2016B079)

    張世偉,女,博士研究生,研究方向:微生物代謝流;E-mail:zhangdan1987114@163.com

    鐘其頂,男,高級工程師,研究方向:食品安全;E-mail:zhongqiding@163.com

    白志輝,男,研究員,研究方向:環(huán)境生物技術與微生物生態(tài)學;E-mail:zhbai@rcees.ac.cn

    猜你喜歡
    釀酒葡萄酒群落
    大學生牙齦炎齦上菌斑的微生物群落
    上半年釀酒產業(yè)產、銷、利均增長
    釀酒科技(2021年8期)2021-12-06 15:28:22
    自制的葡萄酒為啥愛“上頭”?
    合成微生物群落在發(fā)酵食品中的應用研究
    為什么酵母菌既能做面包也能釀酒?
    十款葡萄酒與十塊石頭
    收藏界(2018年3期)2018-10-10 05:34:08
    法國葡萄酒何以譽滿天下
    中國商界(2017年4期)2017-05-17 04:36:48
    釀酒忘米
    春季和夏季巢湖浮游生物群落組成及其動態(tài)分析
    調虧灌溉在釀酒葡萄上的應用
    久久热精品热| 天天一区二区日本电影三级| 天堂网av新在线| 亚洲精品乱码久久久v下载方式| 久久久久性生活片| 欧美日韩亚洲国产一区二区在线观看| 啦啦啦观看免费观看视频高清| 别揉我奶头 嗯啊视频| 久久草成人影院| 人妻久久中文字幕网| 国产精品久久电影中文字幕| 热99re8久久精品国产| 日本欧美国产在线视频| 国产一区二区激情短视频| 亚洲精品日韩av片在线观看| 亚洲国产精品合色在线| 免费观看在线日韩| 人妻夜夜爽99麻豆av| 亚洲av二区三区四区| 成人国产综合亚洲| 身体一侧抽搐| 成人永久免费在线观看视频| 欧美不卡视频在线免费观看| 国产精品国产高清国产av| 国产欧美日韩精品一区二区| 亚洲中文日韩欧美视频| 国产精品野战在线观看| 国产视频内射| 一级毛片久久久久久久久女| 99国产精品一区二区蜜桃av| 俄罗斯特黄特色一大片| 99热网站在线观看| 99热只有精品国产| 男女下面进入的视频免费午夜| 伊人久久精品亚洲午夜| 可以在线观看毛片的网站| 美女被艹到高潮喷水动态| 亚洲性夜色夜夜综合| 久久久久久九九精品二区国产| 一个人看的www免费观看视频| 欧美高清成人免费视频www| 淫妇啪啪啪对白视频| 真人做人爱边吃奶动态| 亚洲成人久久爱视频| 美女cb高潮喷水在线观看| 我要搜黄色片| 88av欧美| 露出奶头的视频| 琪琪午夜伦伦电影理论片6080| a级毛片a级免费在线| 亚洲熟妇熟女久久| 国产伦精品一区二区三区视频9| 美女免费视频网站| 日韩欧美 国产精品| 亚洲,欧美,日韩| 国产精品电影一区二区三区| 婷婷亚洲欧美| 看片在线看免费视频| 国产精品乱码一区二三区的特点| 我的老师免费观看完整版| 美女免费视频网站| 色哟哟·www| 日日干狠狠操夜夜爽| 毛片一级片免费看久久久久 | 久久久色成人| 久9热在线精品视频| 最近在线观看免费完整版| 男女视频在线观看网站免费| 亚洲av熟女| 亚洲国产欧美人成| 国产国拍精品亚洲av在线观看| 18禁黄网站禁片免费观看直播| 欧美+亚洲+日韩+国产| 婷婷六月久久综合丁香| 男女边吃奶边做爰视频| www.www免费av| 国模一区二区三区四区视频| 午夜精品在线福利| 97超级碰碰碰精品色视频在线观看| 日本在线视频免费播放| 亚洲精华国产精华液的使用体验 | 我的女老师完整版在线观看| 国国产精品蜜臀av免费| av在线天堂中文字幕| 最近中文字幕高清免费大全6 | 亚洲图色成人| 人妻丰满熟妇av一区二区三区| 欧美不卡视频在线免费观看| 日本-黄色视频高清免费观看| 中文字幕av成人在线电影| www日本黄色视频网| 色综合亚洲欧美另类图片| 少妇人妻一区二区三区视频| 亚洲av美国av| 毛片一级片免费看久久久久 | 精品午夜福利视频在线观看一区| 欧美一级a爱片免费观看看| 亚洲国产日韩欧美精品在线观看| 久久久久国产精品人妻aⅴ院| 麻豆久久精品国产亚洲av| 2021天堂中文幕一二区在线观| 国产高清激情床上av| 听说在线观看完整版免费高清| 亚洲欧美日韩卡通动漫| xxxwww97欧美| 成人精品一区二区免费| 亚洲中文字幕一区二区三区有码在线看| 久久精品国产亚洲av香蕉五月| 亚洲七黄色美女视频| 国产亚洲欧美98| 少妇的逼水好多| 麻豆久久精品国产亚洲av| 国产91精品成人一区二区三区| 久久久久久久精品吃奶| 国产人妻一区二区三区在| 99热6这里只有精品| 欧美极品一区二区三区四区| 亚洲欧美日韩高清在线视频| 国产精品亚洲美女久久久| 国产精品98久久久久久宅男小说| 免费在线观看影片大全网站| 午夜影院日韩av| 日韩精品有码人妻一区| 有码 亚洲区| 我要搜黄色片| 精品人妻熟女av久视频| 久久人妻av系列| av在线蜜桃| 国产一区二区三区在线臀色熟女| 亚洲美女搞黄在线观看 | 精品一区二区免费观看| 麻豆成人午夜福利视频| 日韩在线高清观看一区二区三区 | 三级毛片av免费| 成人永久免费在线观看视频| 深夜精品福利| 欧美成人免费av一区二区三区| 波野结衣二区三区在线| 天堂动漫精品| 桃色一区二区三区在线观看| 日韩中文字幕欧美一区二区| 欧美性感艳星| 免费av毛片视频| 国模一区二区三区四区视频| 免费人成在线观看视频色| 午夜免费成人在线视频| 成人无遮挡网站| 黄色欧美视频在线观看| 成人三级黄色视频| 麻豆国产av国片精品| 亚洲精品一区av在线观看| 18禁在线播放成人免费| 国产成人a区在线观看| 大型黄色视频在线免费观看| 丝袜美腿在线中文| 免费一级毛片在线播放高清视频| 精品人妻1区二区| 黄色视频,在线免费观看| 特大巨黑吊av在线直播| 中文字幕av成人在线电影| 国产高清有码在线观看视频| 成人综合一区亚洲| 日韩精品有码人妻一区| 热99re8久久精品国产| 又黄又爽又免费观看的视频| 亚洲国产欧美人成| 级片在线观看| 精品久久久久久,| 最近最新免费中文字幕在线| 亚洲午夜理论影院| www.色视频.com| 精品人妻熟女av久视频| 在线观看66精品国产| 变态另类丝袜制服| 99热这里只有精品一区| 婷婷精品国产亚洲av在线| 女人十人毛片免费观看3o分钟| 男女那种视频在线观看| 91av网一区二区| 桃红色精品国产亚洲av| 成人欧美大片| 久久精品国产亚洲av香蕉五月| 在线免费十八禁| 国产精品乱码一区二三区的特点| 国产一区二区激情短视频| 亚洲欧美精品综合久久99| 成人特级黄色片久久久久久久| 中文字幕熟女人妻在线| 国产在视频线在精品| 男女啪啪激烈高潮av片| 高清毛片免费观看视频网站| 国产精品99久久久久久久久| 有码 亚洲区| 校园春色视频在线观看| 大型黄色视频在线免费观看| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 亚洲18禁久久av| 久久精品国产99精品国产亚洲性色| 久久香蕉精品热| 欧美一区二区国产精品久久精品| 国产综合懂色| 亚洲精品成人久久久久久| 桃红色精品国产亚洲av| 亚洲性久久影院| 99久久中文字幕三级久久日本| 国产成人av教育| 亚洲精品一卡2卡三卡4卡5卡| 国产成人福利小说| 国产精品久久电影中文字幕| 日本免费一区二区三区高清不卡| 桃红色精品国产亚洲av| 免费无遮挡裸体视频| 琪琪午夜伦伦电影理论片6080| 蜜桃久久精品国产亚洲av| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 精品人妻一区二区三区麻豆 | 少妇熟女aⅴ在线视频| 不卡一级毛片| 99热网站在线观看| 啪啪无遮挡十八禁网站| 国产老妇女一区| 亚洲精品成人久久久久久| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 成人精品一区二区免费| 12—13女人毛片做爰片一| 欧美黑人欧美精品刺激| 久久久久九九精品影院| 亚洲美女搞黄在线观看 | 十八禁网站免费在线| 欧美丝袜亚洲另类 | 亚洲三级黄色毛片| 国国产精品蜜臀av免费| 日日夜夜操网爽| 免费在线观看影片大全网站| 黄色视频,在线免费观看| 国产精品久久久久久久久免| 成年版毛片免费区| 精品午夜福利在线看| 热99re8久久精品国产| 黄色一级大片看看| 99久久精品热视频| 亚洲精品久久国产高清桃花| 国产欧美日韩精品一区二区| 男人舔奶头视频| 亚洲欧美日韩卡通动漫| 午夜精品一区二区三区免费看| 精品人妻1区二区| 成人无遮挡网站| 麻豆一二三区av精品| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 最近最新中文字幕大全电影3| 九九在线视频观看精品| 午夜a级毛片| 成人午夜高清在线视频| 人妻久久中文字幕网| 韩国av一区二区三区四区| 国产三级在线视频| 91久久精品国产一区二区三区| 简卡轻食公司| 亚洲图色成人| 熟女电影av网| 天堂动漫精品| 在现免费观看毛片| 成人美女网站在线观看视频| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 亚洲av免费高清在线观看| 国产av麻豆久久久久久久| av.在线天堂| 国产精品一及| 能在线免费观看的黄片| 啦啦啦韩国在线观看视频| 日日啪夜夜撸| 一级毛片久久久久久久久女| 欧美极品一区二区三区四区| 永久网站在线| 搡老妇女老女人老熟妇| 久久久久精品国产欧美久久久| 亚洲一区二区三区色噜噜| 天堂√8在线中文| 国产单亲对白刺激| 国产大屁股一区二区在线视频| 国产淫片久久久久久久久| 女的被弄到高潮叫床怎么办 | 欧美成人a在线观看| 日韩欧美在线二视频| 国内精品美女久久久久久| АⅤ资源中文在线天堂| a在线观看视频网站| 真实男女啪啪啪动态图| 波野结衣二区三区在线| 97超视频在线观看视频| 在现免费观看毛片| 天堂影院成人在线观看| 永久网站在线| 日本免费a在线| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 午夜爱爱视频在线播放| 国产女主播在线喷水免费视频网站 | 国产精品三级大全| 免费人成在线观看视频色| 我的女老师完整版在线观看| 国产aⅴ精品一区二区三区波| 日日夜夜操网爽| 搡老熟女国产l中国老女人| 亚洲18禁久久av| 亚洲人成伊人成综合网2020| av天堂中文字幕网| 色精品久久人妻99蜜桃| 一级黄色大片毛片| 亚洲五月天丁香| 看免费成人av毛片| 岛国在线免费视频观看| 黄色配什么色好看| 岛国在线免费视频观看| 黄色一级大片看看| 校园人妻丝袜中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 九九热线精品视视频播放| 啦啦啦韩国在线观看视频| 五月玫瑰六月丁香| 日韩大尺度精品在线看网址| 哪里可以看免费的av片| 成人欧美大片| 欧美日韩中文字幕国产精品一区二区三区| 国产欧美日韩精品一区二区| 国产精品日韩av在线免费观看| 亚洲三级黄色毛片| 床上黄色一级片| 亚洲人成网站高清观看| 欧美精品国产亚洲| 最近最新中文字幕大全电影3| 精品一区二区免费观看| 日韩欧美国产在线观看| 久久婷婷人人爽人人干人人爱| 校园春色视频在线观看| 精品国产三级普通话版| 黄色丝袜av网址大全| 日韩精品有码人妻一区| 亚洲第一电影网av| 色吧在线观看| 一本精品99久久精品77| 丰满人妻一区二区三区视频av| 亚洲第一电影网av| 国产精品免费一区二区三区在线| 成熟少妇高潮喷水视频| 久久热精品热| 国产精品爽爽va在线观看网站| 99热网站在线观看| 春色校园在线视频观看| 日本黄大片高清| 岛国在线免费视频观看| 欧美最黄视频在线播放免费| 2021天堂中文幕一二区在线观| 久久久久九九精品影院| 亚洲午夜理论影院| 久久中文看片网| 亚洲欧美日韩高清专用| 欧美人与善性xxx| 麻豆国产97在线/欧美| 女的被弄到高潮叫床怎么办 | 深爱激情五月婷婷| 久久久久性生活片| 国产精品国产高清国产av| 精品午夜福利在线看| 啦啦啦啦在线视频资源| 麻豆成人av在线观看| 天堂√8在线中文| 午夜福利成人在线免费观看| 两个人的视频大全免费| 丝袜美腿在线中文| 中文字幕av在线有码专区| 99久久中文字幕三级久久日本| 国产v大片淫在线免费观看| 成人三级黄色视频| 女的被弄到高潮叫床怎么办 | 成人美女网站在线观看视频| 国产一区二区在线av高清观看| 日日夜夜操网爽| 久久久久久国产a免费观看| 波多野结衣巨乳人妻| 麻豆国产97在线/欧美| 看片在线看免费视频| 日韩强制内射视频| 国产精品亚洲美女久久久| 观看美女的网站| 人妻制服诱惑在线中文字幕| 两个人视频免费观看高清| 欧美不卡视频在线免费观看| 99riav亚洲国产免费| 国产精品久久久久久久电影| 午夜精品久久久久久毛片777| 国产真实伦视频高清在线观看 | 亚洲欧美日韩无卡精品| 精品99又大又爽又粗少妇毛片 | 国产精品电影一区二区三区| 露出奶头的视频| 日日摸夜夜添夜夜添av毛片 | 国产熟女欧美一区二区| 给我免费播放毛片高清在线观看| 12—13女人毛片做爰片一| 亚洲人与动物交配视频| 精品久久久久久,| 特大巨黑吊av在线直播| 一进一出好大好爽视频| 免费在线观看日本一区| or卡值多少钱| 欧美日韩瑟瑟在线播放| 18禁裸乳无遮挡免费网站照片| 亚洲中文字幕日韩| 可以在线观看的亚洲视频| 午夜福利在线在线| 如何舔出高潮| 日日撸夜夜添| 成人欧美大片| 我要搜黄色片| 免费av不卡在线播放| 男人的好看免费观看在线视频| 国产午夜福利久久久久久| 简卡轻食公司| 无遮挡黄片免费观看| 亚洲欧美精品综合久久99| 久久精品国产亚洲av香蕉五月| 国产又黄又爽又无遮挡在线| 国产国拍精品亚洲av在线观看| 中文字幕高清在线视频| 欧美成人性av电影在线观看| 国产精品亚洲一级av第二区| www.www免费av| 国产精品98久久久久久宅男小说| 91在线精品国自产拍蜜月| 九色成人免费人妻av| 国产美女午夜福利| 麻豆成人av在线观看| 亚洲无线在线观看| 真实男女啪啪啪动态图| 日韩精品有码人妻一区| 免费看av在线观看网站| 男女啪啪激烈高潮av片| 一本一本综合久久| 成年女人看的毛片在线观看| 免费搜索国产男女视频| 熟妇人妻久久中文字幕3abv| 99在线人妻在线中文字幕| 国产激情偷乱视频一区二区| 超碰av人人做人人爽久久| 老熟妇乱子伦视频在线观看| 一区福利在线观看| 精品人妻一区二区三区麻豆 | 久久精品国产亚洲av香蕉五月| 国产在线精品亚洲第一网站| 男女边吃奶边做爰视频| 搡老熟女国产l中国老女人| 欧美在线一区亚洲| 亚洲国产高清在线一区二区三| 亚洲成人久久爱视频| 日韩欧美精品v在线| 别揉我奶头 嗯啊视频| 波多野结衣高清作品| 中文字幕高清在线视频| 99国产精品一区二区蜜桃av| 国产一区二区三区视频了| 国产真实乱freesex| 国产午夜精品久久久久久一区二区三区 | 国产精华一区二区三区| 午夜福利在线在线| 亚洲国产精品久久男人天堂| 美女黄网站色视频| 亚洲成av人片在线播放无| 欧美zozozo另类| 成年女人看的毛片在线观看| 国产视频内射| 国产精品久久久久久精品电影| 亚洲欧美日韩高清在线视频| 偷拍熟女少妇极品色| 国产精品一区二区三区四区免费观看 | 99精品久久久久人妻精品| 午夜影院日韩av| 内地一区二区视频在线| 长腿黑丝高跟| 免费看光身美女| 成人无遮挡网站| 国产乱人视频| 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| 久久久久免费精品人妻一区二区| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区性色av| 最新在线观看一区二区三区| 男女做爰动态图高潮gif福利片| 淫妇啪啪啪对白视频| 乱人视频在线观看| 又紧又爽又黄一区二区| 国产美女午夜福利| 男女下面进入的视频免费午夜| 午夜精品一区二区三区免费看| 俄罗斯特黄特色一大片| 亚洲人成网站在线播放欧美日韩| 亚洲欧美激情综合另类| 99热网站在线观看| 免费av观看视频| 毛片一级片免费看久久久久 | 一区二区三区免费毛片| 日韩精品青青久久久久久| 婷婷丁香在线五月| 天堂√8在线中文| 久久精品国产亚洲av涩爱 | 狠狠狠狠99中文字幕| 成年女人看的毛片在线观看| 久久久国产成人精品二区| 国产91精品成人一区二区三区| 午夜免费男女啪啪视频观看 | 国内久久婷婷六月综合欲色啪| 亚洲av成人av| 日韩国内少妇激情av| 国产亚洲av嫩草精品影院| 国产高清视频在线观看网站| 成人永久免费在线观看视频| 国产亚洲av嫩草精品影院| 国产精品福利在线免费观看| 国产aⅴ精品一区二区三区波| 嫩草影院精品99| 欧美最新免费一区二区三区| x7x7x7水蜜桃| 久久精品国产亚洲av香蕉五月| 真人做人爱边吃奶动态| 亚洲五月天丁香| 日韩精品青青久久久久久| 别揉我奶头 嗯啊视频| 琪琪午夜伦伦电影理论片6080| 免费无遮挡裸体视频| 国产精品美女特级片免费视频播放器| 国产不卡一卡二| 免费人成在线观看视频色| 女人十人毛片免费观看3o分钟| 日韩欧美 国产精品| 亚洲avbb在线观看| 亚洲av日韩精品久久久久久密| 黄色日韩在线| x7x7x7水蜜桃| 美女黄网站色视频| 国产免费男女视频| 一本一本综合久久| 国产一区二区在线观看日韩| 免费一级毛片在线播放高清视频| 男女视频在线观看网站免费| 国产三级中文精品| 极品教师在线视频| 99在线人妻在线中文字幕| 99久久九九国产精品国产免费| 国产成人aa在线观看| 男女做爰动态图高潮gif福利片| 日韩欧美在线乱码| 男人的好看免费观看在线视频| 亚洲成人精品中文字幕电影| 黄色欧美视频在线观看| 91久久精品国产一区二区成人| 老熟妇仑乱视频hdxx| 国产久久久一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产高清不卡午夜福利| xxxwww97欧美| 天堂动漫精品| 色综合色国产| 亚洲成人久久爱视频| 色综合亚洲欧美另类图片| 国产午夜精品论理片| 一本一本综合久久| 日韩精品有码人妻一区| 亚洲成av人片在线播放无| 春色校园在线视频观看| www日本黄色视频网| 国产伦精品一区二区三区四那| 性插视频无遮挡在线免费观看| 97碰自拍视频| 国产久久久一区二区三区| 国产精品自产拍在线观看55亚洲| 国产午夜福利久久久久久| 午夜福利成人在线免费观看| 综合色av麻豆| 亚洲第一区二区三区不卡| 女的被弄到高潮叫床怎么办 | 成年版毛片免费区| 日本撒尿小便嘘嘘汇集6| 日本精品一区二区三区蜜桃| 欧美不卡视频在线免费观看| 久久精品国产亚洲网站| 亚洲自偷自拍三级| 成年版毛片免费区| 欧美色欧美亚洲另类二区| 亚洲精品色激情综合| 久久久久国产精品人妻aⅴ院| 少妇人妻精品综合一区二区 | 亚洲中文字幕一区二区三区有码在线看| 人妻久久中文字幕网| 淫秽高清视频在线观看| 欧美三级亚洲精品| 99热这里只有精品一区| av在线天堂中文字幕| 一区二区三区激情视频| 人妻久久中文字幕网| 亚洲狠狠婷婷综合久久图片| 亚洲第一区二区三区不卡| 国产精品久久电影中文字幕| 国产精品精品国产色婷婷| 两个人的视频大全免费| 精品久久久久久久久久久久久| 少妇裸体淫交视频免费看高清| 在线观看av片永久免费下载| 观看免费一级毛片|