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

    基于高通量測(cè)序的貯藏白牡丹表面細(xì)菌多樣性分析

    2020-08-27 07:28:30陳佳佳胡雲(yún)飛沈詩(shī)鈺王志華周喆唐琴孫威江1
    茶葉科學(xué) 2020年4期
    關(guān)鍵詞:白牡丹白茶菌門(mén)

    陳佳佳,胡雲(yún)飛,沈詩(shī)鈺,王志華,周喆,唐琴,孫威江1,*

    基于高通量測(cè)序的貯藏白牡丹表面細(xì)菌多樣性分析

    陳佳佳1,2,4,胡雲(yún)飛2,3,4,沈詩(shī)鈺1,2,4,王志華1,2,4,周喆1,2,4,唐琴1,2,4,孫威江1,2,3,4*

    1. 福建農(nóng)林大學(xué)園藝學(xué)院,福建 福州 352000;2. 福建省行業(yè)技術(shù)開(kāi)發(fā)基地,福建 福州 350002;3. 福建農(nóng)林大學(xué)安溪茶學(xué)院,福建 泉州 362400;4. 福建省茶產(chǎn)業(yè)工程技術(shù)研究中心,福建 福州 350002

    為分析貯藏白牡丹細(xì)菌多樣性,利用高通量測(cè)序?qū)?1個(gè)不同貯藏年份白牡丹細(xì)菌進(jìn)行16?S rDNA檢測(cè)。結(jié)果表明,貯藏茶樣的優(yōu)勢(shì)菌門(mén)為變形菌門(mén)()、擬桿菌門(mén)()、厚壁菌門(mén)(等,優(yōu)勢(shì)菌綱為-變形菌綱()、-變形菌綱()、鞘脂桿菌綱()等,優(yōu)勢(shì)菌屬為鞘氨醇單胞菌屬()、甲基桿菌屬()、土地桿菌屬()等。不同貯藏年份白牡丹細(xì)菌組成結(jié)構(gòu)相似,豐度存在差異;貯藏年份對(duì)白牡丹表面細(xì)菌組成未有顯著影響,而貯藏環(huán)境影響細(xì)菌多樣性。

    高通量測(cè)序;白牡丹;細(xì)菌;多樣性

    白茶是一種輕發(fā)酵茶類(lèi),制作簡(jiǎn)單,只有萎凋及干燥工序[1]。貯藏白茶是在適當(dāng)條件下倉(cāng)儲(chǔ)若干時(shí)間的陳年白茶。貯藏黑茶較為常見(jiàn),黑茶中已發(fā)現(xiàn)存在大量微生物,如冠突散囊菌()、曲霉屬()、假單胞菌屬()等[2-3]。微生物在黑茶加工過(guò)程中具有重要作用,關(guān)于黑茶中微生物的研究已經(jīng)大量開(kāi)展。研究證明,貯藏年份[4]、加工過(guò)程[5]、原料[3]、產(chǎn)地[6]均會(huì)對(duì)茶葉微生物群落組成產(chǎn)生影響。楊雅焯等[7]研究發(fā)現(xiàn),廣西六堡茶的優(yōu)勢(shì)菌為鞘脂單胞菌屬()、、乳酸菌屬()等,而重慶沱茶無(wú)明顯優(yōu)勢(shì)菌。在下關(guān)沱茶發(fā)酵過(guò)程,發(fā)酵前期霉菌為優(yōu)勢(shì)菌屬,發(fā)酵后期酵母則占主導(dǎo)[8]。對(duì)不同年份普洱茶微生物群落的研究發(fā)現(xiàn),芽孢桿菌屬()為其主要優(yōu)勢(shì)菌,其次為賴(lài)氏桿菌屬()[4]。茶葉中微生物來(lái)源豐富,黑茶中微生物可能來(lái)自于茶樹(shù)內(nèi)生菌、茶樹(shù)生長(zhǎng)環(huán)境和加工環(huán)境等。路偉堯[8]發(fā)現(xiàn)普洱茶渥堆過(guò)程中的部分優(yōu)勢(shì)菌與生長(zhǎng)環(huán)境中空氣微生物、茶樹(shù)內(nèi)生菌、車(chē)間環(huán)境(空氣、地面、發(fā)酵用水)微生物一致,推測(cè)普洱茶渥堆發(fā)酵過(guò)程中的微生物部分來(lái)自于茶樹(shù)的生長(zhǎng)環(huán)境、茶葉內(nèi)生菌以及發(fā)酵車(chē)間。陳麗瑩等[9]檢測(cè)到茶樹(shù)大量?jī)?nèi)生菌在土壤中存在,推測(cè)土壤是內(nèi)生菌的可能來(lái)源。因此,白茶在貯藏過(guò)程中的微生物也可能存在多種來(lái)源。

    在黑茶、煙草的貯藏過(guò)程中,均發(fā)現(xiàn)微生物多樣性變化對(duì)其品質(zhì)產(chǎn)生影響[10-12]。目前,國(guó)內(nèi)大部分的研究仍是基于黑茶,對(duì)于白茶微生物研究較少,且貯藏白茶細(xì)菌多樣性還尚未見(jiàn)報(bào)道,在白茶后期貯藏過(guò)程中,細(xì)菌多樣性的變化值得探討。因此,本試驗(yàn)選取了11個(gè)不同貯藏年份的白牡丹樣品,通過(guò)對(duì)其表面細(xì)菌多樣性進(jìn)行分析,揭示貯藏過(guò)程中細(xì)菌群落結(jié)構(gòu)變化,以期為探究白茶貯藏條件提供參考。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)樣品為福建品品香茶業(yè)有限公司、福建省裕榮香茶業(yè)有限公司、福建省天湖茶業(yè)有限公司、福鼎市張?jiān)洸铇I(yè)有限公司和福建鼎白茶業(yè)有限公司的不同貯藏年份白牡丹茶樣。貯藏地點(diǎn)均為福建省福鼎市,同一企業(yè)茶樣貯藏環(huán)境一致,樣品收集完畢后于實(shí)驗(yàn)室常溫密封保存,濕度小于50%。取樣過(guò)程的自封袋、錫箔紙、手套等均經(jīng)過(guò)無(wú)菌處理。具體樣品信息如表1所示,以2017年生產(chǎn)的白牡丹為對(duì)照樣品(CK)。

    表1 不同貯藏年份白牡丹樣品表

    1.2 試驗(yàn)試劑

    DNA Marker,日本Takara公司;E.Z.N.A.?Soil DNA Kitomega試劑盒,美國(guó)Omega公司;NaCl、Na2HPO4?12H2O、KH2PO4、50×TAE緩沖液,北京索萊寶科技有限公司;瓊脂糖,南京生興生物技術(shù)有限公司;rTaq DNA聚合酶試劑盒,北京全式金生物技術(shù)有限公司;AxyPrep DNA Gel Extraction Kit,美國(guó)Axygen公司。

    PBS緩沖液配制:稱(chēng)取NaCl 8?g,KCl 0.2?g,Na2HPO4?12H2O 3.62?g,KH2PO40.24?g,溶解于800?mL蒸餾水中,用鹽酸調(diào)節(jié)pH為7.4,然后加水定容至1?000?mL。

    1.3 試驗(yàn)方法

    1.3.1 不同年份白牡丹表面細(xì)菌富集及DNA提取

    稱(chēng)取5?g混合好的樣品于無(wú)菌容器中,加入50?mL PBS緩沖液,置于搖床上振蕩培養(yǎng)30?min(180?r·min-1,37℃),然后超聲提取15?min,取出搖晃15?s,再次超聲15?min,置于搖床上振蕩培養(yǎng)30?min(180?r·min-1,37℃),過(guò)濾后將濾液轉(zhuǎn)移至新的無(wú)菌管,常溫狀態(tài)12?000?r·min-1離心10?min,棄上清液,將下層液體轉(zhuǎn)移到2.0?mL無(wú)菌離心管中,12?000?r·min-1離心10?min,棄上清液,沉淀進(jìn)行DNA提取。根據(jù)E.Z.N.A.?Soil DNA Kit試劑盒說(shuō)明書(shū)進(jìn)行總DNA抽提。

    1.3.2 測(cè)序建庫(kù)

    提取樣品DNA后,根據(jù)V3和V4保守區(qū)設(shè)計(jì)得到引物(分別為338F和806R,338F:5'-ACTCCTACGGGAGGCAGCA-3',806R:5'-GGACTACHVGGGTWTCTAAT-3'),進(jìn)行PCR擴(kuò)增。PCR擴(kuò)增程序?yàn)椋?5℃預(yù)變性5?min;95℃變性30?s,50℃退火30?s,72℃延伸30?s,30個(gè)循環(huán);72℃穩(wěn)定延伸7?min;4℃保存。PCR產(chǎn)物使用2%凝膠回收,根據(jù)AxyPrep DNA Gel Extraction Kit說(shuō)明書(shū)純化產(chǎn)物,檢測(cè)定量,使用NEXTFLEX Rapid DNA-Seq Kit進(jìn)行建庫(kù),委托北京百邁克生物公司進(jìn)行測(cè)序。

    1.3.3 數(shù)據(jù)分析

    使用FLASH v1.2.7軟件,通過(guò)overlap對(duì)每個(gè)樣品的雙端reads(PE reads)進(jìn)行拼接,得到原始序列(Raw tags)。使用Trimmomatic v0.33軟件,對(duì)原始序列進(jìn)行過(guò)濾,得到優(yōu)化序列(Clean tags)。使用UCHIME v4.2軟件,鑒定并去除嵌合體序列,得到最終有效序列(Effective tags)。分別計(jì)算有效率(Effective,有效序列占原始序列的比例)、GC含量(G和C類(lèi)型的堿基占總堿基的百分比)、Q20(質(zhì)量值大于等于20的堿基占總堿基數(shù)的百分比)和Q30(質(zhì)量值大于等于30的堿基占總堿基數(shù)的百分比)。使用QIIME(version 1.8.0)軟件中的UCLUST對(duì)Tags在97%的相似度水平下進(jìn)行聚類(lèi)、獲得OTU(Operational taxonomic unit),并基于Silva(細(xì)菌)分類(lèi)學(xué)數(shù)據(jù)庫(kù)對(duì)OTU進(jìn)行分類(lèi)學(xué)注釋。使用Mothur v.1.30軟件進(jìn)行Alpha多樣性指數(shù)分析,作稀釋性曲線(xiàn)(Rarefaction Curve)[13],進(jìn)行Chao1、Ace、Shannon、Simpson指標(biāo)計(jì)算[14]。通過(guò)Graph Prism 8.0、R軟件作圖。

    2 結(jié)果與分析

    2.1 測(cè)序質(zhì)量分析

    由表2可知,11個(gè)樣品測(cè)序共獲得784?599對(duì)序列,703?966條優(yōu)化序列,每個(gè)樣品至少產(chǎn)生39?027條優(yōu)化序列;有效率均>85%,有效序列占原始序列比值大;樣品GC含量范圍為51.75%~55.20%,不存在明顯的GC偏好性;Q20值均≥96.24%,Q30值均≥92.71%。綜上所述,測(cè)序質(zhì)量較好,可以用于后續(xù)試驗(yàn)分析。

    2.2 OTU分析

    對(duì)OTU分析,共獲得1?366個(gè)OTUs,各樣品OTU數(shù)在147~891之間,樣品間OTU數(shù)存在差異(圖1)。對(duì)同一企業(yè)的白牡丹樣品進(jìn)行比較,來(lái)自裕榮香的白牡丹樣OTU數(shù)呈現(xiàn)為:YRX16

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

    圖1 11個(gè)樣品OTU個(gè)數(shù)分布圖

    2.3 樣品豐富度分析

    2.3.1 稀釋性曲線(xiàn)

    稀釋性曲線(xiàn)(Rarefaction curve)可驗(yàn)證測(cè)序數(shù)據(jù)量是否足以反映樣品中的物種多樣性,反映樣品中物種的豐富程度與持續(xù)抽樣下新物種出現(xiàn)的速率。如圖2所示,樣本測(cè)序數(shù)量小于1?000時(shí),隨測(cè)序數(shù)量增加,PPX07、YRX07、PPX08、PPX11、YRX11、DB13的OTU個(gè)數(shù)急速上升,表示群落中鑒定到大量物種,測(cè)序數(shù)量大于1?000時(shí),OTU增加趨于平緩,環(huán)境中物種不會(huì)隨測(cè)序數(shù)量增加而顯著增加,測(cè)序序列充分;樣本測(cè)序數(shù)量小于4?000時(shí),TH12、ZYJ14、PPX16、YRX16、CK的OTU個(gè)數(shù)急速上升,大于4?000后,OTU增加趨于平緩,該5個(gè)樣本測(cè)序數(shù)據(jù)量大于上述6個(gè)樣本。綜上所述,11個(gè)樣本序列充分,可以用于數(shù)據(jù)分析。

    2.3.2 樣品間Alpha多樣性分析

    Chao1和Ace指數(shù)用于衡量物種豐度,Shannon和Simpson指數(shù)用于衡量物種多樣性。如表3所示,11個(gè)樣品中DB13的Chao1和Ace指數(shù)最高,分別為917.538?5和904.756?9,表明DB13的物種豐度最高;其Shannon指數(shù)在11個(gè)樣品中最高,Simpson指數(shù)最低,分別為5.374?0和0.017?2,說(shuō)明DB13物種多樣性最高。PPX16的Chao1和Ace指數(shù)最低,分別為206.368?4和201.113?5,表明其物種豐度最低。物種多樣性最低的是ZYJ14,其Shannon指數(shù)最低,Simpson指數(shù)最高。整體上,各樣品的Chao1、Ace、Shannon指數(shù)均較高,Simpson指數(shù)較低,說(shuō)明物種豐富度及多樣性高,表明不同年份白牡丹中均存在著物種豐富的細(xì)菌群落。OTU覆蓋率反映樣本中微生物的真實(shí)情況,數(shù)值越高,則樣本中物種被測(cè)出的概率越高。各樣品的覆蓋率均大于0.998?7,說(shuō)明絕大部分物種可以被檢測(cè)出,證實(shí)測(cè)序質(zhì)量良好,可以滿(mǎn)足后續(xù)分析要求。

    圖2 11個(gè)樣品稀釋性曲線(xiàn)圖

    表3 貯藏白牡丹Alpha多樣性指數(shù)

    2.4 樣品間細(xì)菌群落結(jié)構(gòu)分析

    2.4.1 細(xì)菌群落分析

    11個(gè)白牡丹茶樣的高通量測(cè)序結(jié)果共鑒定出266種細(xì)菌,涉及2個(gè)界,24個(gè)門(mén),54個(gè)綱,102個(gè)目,194個(gè)科,373個(gè)屬。由圖3-A可知,在門(mén)水平上,優(yōu)勢(shì)菌門(mén)為變形菌門(mén)()、擬桿菌門(mén)()、厚壁菌門(mén)()等,其中變形菌門(mén)是最主要的優(yōu)勢(shì)菌門(mén),在所有樣品中占比均大于40%,在TH12、PPX07中的占比大于90%。由圖3-B可知,在綱水平上,貯藏茶樣中的優(yōu)勢(shì)菌綱為-變形菌綱()、-變形菌綱()、鞘脂桿菌綱()、-變形菌綱()、梭菌綱()、擬桿菌綱()、芽孢桿菌綱(),其余一些不可識(shí)別以及不可培養(yǎng)細(xì)菌的占比極小。對(duì)上述結(jié)果分析發(fā)現(xiàn),在門(mén)及綱水平上,相同貯藏年份的白牡丹樣品物種組成相似,豐度上存在顯著差異,表明物種豐度和年份之間無(wú)顯著相關(guān)性,推測(cè)白牡丹表面的細(xì)菌可能來(lái)自于環(huán)境,主要受環(huán)境因素的影響。

    注:A為細(xì)菌門(mén)分類(lèi);B為細(xì)菌綱分類(lèi)

    2.4.2 屬水平群落分析

    選取屬水平上豐度大于5%的物種進(jìn)行分析發(fā)現(xiàn)(圖4),不同貯藏年份的白牡丹茶樣細(xì)菌群落組成相似,豐度較高的菌屬均為鞘氨醇單胞菌屬()、甲基桿菌屬()、土地桿菌屬()、貪噬菌屬()、中慢生根瘤菌屬()、布赫納氏菌屬()、沙雷菌屬()、根瘤菌屬(),以及毛螺菌科()中未分類(lèi)的一些菌屬。但是在細(xì)菌豐度、優(yōu)勢(shì)菌屬上存在差異,CK中的優(yōu)勢(shì)菌屬為貪噬菌屬()、鞘氨醇單胞菌屬()、土地桿菌屬()、甲基桿菌屬(),而YRX07中的優(yōu)勢(shì)菌屬為中慢生根瘤菌屬()、鞘氨醇單胞菌屬()以及毛螺菌科下未分類(lèi)菌屬。屬水平下,不同貯藏年份中的細(xì)菌群落結(jié)構(gòu)組成相似,物種豐度上存在差異,無(wú)明顯規(guī)律,說(shuō)明貯藏年份對(duì)細(xì)菌菌落多樣性及結(jié)構(gòu)影響不顯著,其年份并非影響菌落多樣性的關(guān)鍵因素,與上述在門(mén)和綱水平上的分析結(jié)果一致。

    圖4 不同貯藏年份白牡丹細(xì)菌屬水平的相對(duì)豐度統(tǒng)計(jì)分析

    3 討論

    3.1 貯藏白牡丹細(xì)菌群落結(jié)構(gòu)

    貯藏白牡丹細(xì)菌群落多樣性高,豐度大;細(xì)菌群落結(jié)構(gòu)相似,物種豐度存在差異。高通量測(cè)序結(jié)果表明,貯藏茶樣中優(yōu)勢(shì)菌門(mén)主要為變形菌門(mén)()、擬桿菌門(mén)()、厚壁菌門(mén)()等;優(yōu)勢(shì)菌綱為-變形菌綱()、-變形菌綱()、鞘脂桿菌綱()、-變形菌綱();優(yōu)勢(shì)菌屬為鞘氨醇單胞菌屬()、甲基桿菌屬()。不同貯藏年份白牡丹樣品的細(xì)菌分析結(jié)果表明,群落結(jié)構(gòu)未有規(guī)律性變化,菌屬種類(lèi)與貯藏年份無(wú)顯著相關(guān)性,推測(cè)白茶表面的細(xì)菌主要來(lái)源于貯藏環(huán)境,而非加工過(guò)程,煙草的研究中也存在相似結(jié)果[12]。而黑茶后期貯藏中的優(yōu)勢(shì)菌與渥堆過(guò)程中微生物有直接相關(guān)性,渥堆促進(jìn)部分微生物大量繁殖成為優(yōu)勢(shì)菌屬,這些優(yōu)勢(shì)菌屬在后續(xù)干燥等加工過(guò)程中依然可以存活,隨后在后期貯藏過(guò)程中這些優(yōu)勢(shì)微生物再次繁殖,成為優(yōu)勢(shì)菌,使得后期貯藏與渥堆過(guò)程中黑茶表面呈現(xiàn)的微生物群落結(jié)構(gòu)相似。研究發(fā)現(xiàn)芽孢桿菌屬()是普洱茶渥堆發(fā)酵及后期貯藏過(guò)程中的優(yōu)勢(shì)菌屬,可證實(shí)黑茶后期貯藏過(guò)程中的微生物主要來(lái)源于加工過(guò)程[10,15]。已有研究顯示,這些微生物可促進(jìn)黑茶品質(zhì)的形成,改善后期加工及貯藏中黑茶的滋味[16]。因此,在白茶貯藏過(guò)程中,尤其是一些粗老的白茶,可利用微生物進(jìn)一步提升其品質(zhì),改善其苦澀的口感,增強(qiáng)甜醇的滋味。

    3.2 貯藏環(huán)境對(duì)白牡丹細(xì)菌群落的影響

    不同貯藏條件對(duì)微生物群落及動(dòng)態(tài)變化有顯著影響,不同貯藏環(huán)境下的煙草在細(xì)菌群落上存在明顯差異[17]。研究表明,適宜的溫度、水活度、氧氣、pH和充足的養(yǎng)分會(huì)促進(jìn)微生物的生長(zhǎng)[18-19]。大部分細(xì)菌最適培養(yǎng)溫度為25℃~37℃,水活度為0.94以上[20],本試驗(yàn)中所有白牡丹貯藏環(huán)境均為常溫密封保存,濕度小于50%(水活度小于0.5),該貯藏環(huán)境不利于微生物的生長(zhǎng);微生物生長(zhǎng)還需要豐富的營(yíng)養(yǎng)基質(zhì),茶葉中(碳水化合物、蛋白質(zhì)等)營(yíng)養(yǎng)基質(zhì)貧乏,同時(shí)白茶在加工過(guò)程中經(jīng)過(guò)干燥后,水分含量低,不適合微生物的繁殖,并且多酚等物質(zhì)具有顯著抑菌作用[21],抑制微生物的生長(zhǎng)。綜上所述,白牡丹貯藏時(shí)微生物不易大量繁殖,未形成優(yōu)勢(shì)菌屬。本研究發(fā)現(xiàn),貯藏白牡丹表面細(xì)菌中豐度較高的是一種廣泛存在于環(huán)境中的細(xì)菌——鞘氨醇單胞菌屬(),除此之外未檢測(cè)到其他在貯藏白茶中大量富集的細(xì)菌群落,證實(shí)在該白茶貯藏環(huán)境中細(xì)菌的群落相對(duì)穩(wěn)定。貯藏條件的改變均有可能引起細(xì)菌多樣性及豐度變化。近年來(lái),已有大量研究對(duì)貯藏條件進(jìn)行探究,旨在通過(guò)調(diào)控貯藏條件改善茶葉品質(zhì),研究結(jié)果顯示適當(dāng)溫、濕度可以加快貯藏白茶品質(zhì)的轉(zhuǎn)變[22],但是白茶貯藏轉(zhuǎn)化需要適度環(huán)境,要兼顧微生物安全。因?yàn)榄h(huán)境條件會(huì)影響細(xì)菌生長(zhǎng),因此在白茶貯藏過(guò)程中要嚴(yán)格把控環(huán)境因子,尤其是溫度、濕度與氧氣,加強(qiáng)安全監(jiān)管,減少有害細(xì)菌富集。

    [1] 劉誼健, 郭玉瓊, 詹梓金. 白茶制作過(guò)程主要化學(xué)成分轉(zhuǎn)化與品質(zhì)形成探討[J]. 福建茶葉, 2003(4): 13-14. Liu Y J, Guo Y Q, Zang Z J. Analysis of the transformation of main chemical components and quality formation of white tea [J]. Tea in Fujian, 2003(4): 13-14.

    [2] 皇甫瑞娟. 茯磚茶中冠突散囊菌的分離鑒定[J]. 現(xiàn)代食品, 2019(23): 183-185.Huangfu R J. Isolation and identification offrom Fuzhuan tea [J]. Modern Food, 2019(23): 183-185.

    [3] 曾橋, 呂生華, 李祥, 等. 不同原料茯磚茶活性成分及微生物多樣性分析[J/OL]. 食品科學(xué), 2020: 1-12. http://kns.cnki.net/kcms/detail/11.2206.TS.20191219.1855.018.html.Zeng Q, Lv S H, Li X, et al. Active ingredients and microbial diversity of Fu brick tea produced by different raw materials [J/OL]. Food Science, 2020: 1-12. http://kns.cnki.net/kcms/detail/11.2206.TS.20191219.1855.018.html.

    [4] 葛慈斌,劉波,陳梅春, 等. 不同年份陳年普洱茶中芽胞桿菌的多樣性[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 2017, 25(8): 1336-1348. Ge C B, Liu B, Chen M C, et al. Diversity of the bacillus-like species in the different year-old fermented Pu-erh tea [J]. Journal of Agricultural Biotechnology, 2017, 25(8): 1336-1348.

    [5] 方欣, 駱愛(ài)國(guó), 涂青, 等. 普洱茶(熟茶)發(fā)酵過(guò)程各層間真菌群落的動(dòng)態(tài)變化[J]. 食品科技, 2019, 44(5): 37-42.Fang X, Luo A G, Tu Q, et al. Fungal community dynamic change in different layers of solid-state fermentation of Pu-erh ripe tea [J]. Food Science and Technology, 2019, 44(5): 37-42.

    [6] 楊吉霞, 曾祥平, 蒲慧敏, 等. 陜西茯磚茶微生物多樣性和群落結(jié)構(gòu)的研究[J]. 食品與發(fā)酵工業(yè), 2020, 46(3): 50-57. Yang J X, Zeng X P, Pu H M, et al. An investigation of the microbial diversity and community structure in Shaanxi Fu brick tea [J]. Food and Fermentation Industries, 2020, 46(3): 50-57.

    [7] 楊雅焯, 汪迎, 李輝, 等. 廣西六堡茶和重慶沱茶的微生物多樣性分析[J]. 茶葉學(xué)報(bào), 2019(3): 93-98. Yang Y Z, Wang Y, Li H, et al. Microbial diversity of Guangxi Liubao and Chongqing Bowl teas [J]. Acta Tea Sinica, 2019(3): 93-98.

    [8] 路偉堯. 普洱茶發(fā)酵微生物的溯源分析[D]. 北京: 北京化工大學(xué), 2013. Lu W Y. Analysis of microbial tracking in Pu’er tea fermentation [D]. Beijing: Beijing University of Chemical Technology, 2013.

    [9] 陳麗瑩, 張玉滿(mǎn), 陳曉英, 等. 茶樹(shù)葉際選擇性富集的內(nèi)生細(xì)菌的鑒定[J]. 微生物學(xué)報(bào), 2018, 58(10): 1776-1785. Chen L Y, Zhang Y M, Chen X Y, et al. Identification of endophytic bacteria selectively enriched inleaf [J]. Microbiology China, 2018, 58(10): 1776-1785.

    [10] 林長(zhǎng)欣. 普洱茶中的風(fēng)味成分及微生物在貯藏過(guò)程中的變化[D]. 廣州: 暨南大學(xué), 2010. Lin C X. Flavor compounds and microorganisms in Pu'er tea at different storage period [D]. Guangzhou: Jinan University, 2010.

    [11] 龔俊. 烤后片煙儲(chǔ)存過(guò)程中微生物多樣性及變化動(dòng)態(tài)[D]. 上海: 華東師范大學(xué), 2015. Gong J. The diversity and dynamic of microorganism on fluecured tobacco leaves during different aged phases [D]. Shanghai: East China Normal University, 2015.

    [12] 周家喜. 不同倉(cāng)儲(chǔ)環(huán)境微生物群落對(duì)煙葉陳化品質(zhì)的影響[D]. 貴陽(yáng): 貴州大學(xué), 2017. Zhou J X. Effects of microbial communities on quality of tobacco leaves in different storage environments [D]. Guiyang: Guizhou University, 2017.

    [13] Wang Y, Sheng H F, He Y, et al. Comparison of the levels of bacterial diversity in freshwater, intertidal wetland, and marine sediments by using Millions of Illumina Tags [J]. Applied and Environmental Microbiology, 2012, 78(23): 8264-8271.

    [14] Kong H, Conlan S, Grice E A, et al. Topographical and temporal diversity of the human skin microbiome [J]. Science, 2009, 324(5931): 1190-1192.

    [15] 姚靜, 陳迪, 鄭曉燕, 等. 普洱茶渥堆發(fā)酵過(guò)程中細(xì)菌種群的分離與分子鑒定[J]. 安徽農(nóng)業(yè)科學(xué), 2013, 41(6): 2667-2671. Yao J, Chen D, Zheng X Y, et al. Isolation and molecular identification of bacterial colonization during the pile fermentation process of Pu-er tea [J]. Journal of Anhui Agricultural Sciences, 2013, 41(6): 2667-2671.

    [16] 徐正剛, 吳良, 劉石泉, 等. 黑茶發(fā)酵過(guò)程中微生物多樣性研究進(jìn)展[J]. 生物學(xué)雜志, 2019, 36(3): 92-95.Xu Z G, Wu L, Liu S Q, et al. Review for development of microbial diversity during dark tea fermentation period [J]. Journal of Biology, 2019, 36(3): 92-95.

    [17] Jessica C, Suhana C, Prachi K, et al. Temporal variations in cigarette tobacco bacterial community composition and tobacco-specific nitrosamine content are influenced by brand and storage conditions [J]. Frontiers in Microbiology, 2017, 8: 358. doi: 10.3389/fmicb.2017.00358.

    [18] Kinkel L. Microbial population dynamics on leaves [J]. Annual Review of Phytopathology, 1997, 35(1): 327-347.

    [19] Di Giacomo M, Paolino M, Silvestro D, et al. Microbial community structure and dynamics of dark fire-cured tobacco fermentation [J]. Applied and Environmental Microbiology, 2007, 73(3): 825-837.

    [20] 何國(guó)慶. 食品微生物學(xué)[M]. 北京: 中國(guó)農(nóng)業(yè)大學(xué)出版社, 2009. He G Q. Food microbiology [M]. Beijing: China Agricultural University Press, 2009.

    [21] Bansal S , Choudhary S , Sharma M , et al. Tea: A native source of antimicrobial agents [J]. Food Research International, 2013, 53(2): 568-584.

    [22] 祁丹丹, 陳維, 苗愛(ài)清, 等. 氣相色譜-串聯(lián)質(zhì)譜結(jié)合化學(xué)計(jì)量學(xué)分析快速陳化對(duì)白茶香氣的影響[J]. 浙江大學(xué)學(xué)報(bào)(農(nóng)業(yè)與生命科學(xué)版), 2018, 44(6): 704-710. Qi D D, Chen W, Miao A Q, et al. Effects of rapid ageing technology on the aroma quality of white tea using gas chromatography-mass spectrometry/mass spectrometry combined with chemometrics [J]. Journal of Zhejiang University (Agric. & Life Sci.), 2018, 44(6): 704-710.

    Analysis of Bacterial Diversity of Stored White Peony by High-Throughput Sequencing

    CHEN Jiajia1,2,4, HU Yunfei2,3,4, SHEN Shiyu1,2,4, WANG Zhihua1,2,4, ZHOU Zhe1,2,4, TANG Qin1,2,4, SUN Weijiang1,2,3,4*

    1. College of Horticulture, Fujian Agriculture and Forestry University, Fuzhou 352000, China; 2. Industrial Technology Development Base in Fujian Province, Fuzhou 350002, China; 3. Anxi College of Tea Science, Fujian Agriculture and Forestry University, Quanzhou 362400, China; 4. Tea Industry Engineering Technology Research Center of Fujian Province, Fuzhou 350002, China

    To investigate the diversity of the bacterial communities in stored white peony, 16?S rDNA of bacteria collected from 11 samples of different storage durations were analyzed by high-throughput sequencing. The results show that the species richness and diversity of DB13 were the highest. While in all stored white peony, the dominant phylum wereand. The dominant class wereand.And the dominant bacterial species at genus level include:.The composition structures of samples were similar, but there was difference in their abundance. Storage duration had no significant effect on the bacteria composition of white peony, but environments played important roles.

    high-throughput sequencing, white peony, bacteria, diversity

    S571.1

    A

    1000-369X(2020)04-519-09

    2019-12-23

    2020-03-02

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFF0201904-03)、國(guó)家自然科學(xué)基金項(xiàng)目(31770732)、福建農(nóng)林大學(xué)茶產(chǎn)業(yè)鏈科技創(chuàng)新與服務(wù)體系建設(shè)

    陳佳佳,女,碩士研究生,主要從事茶葉加工方面的研究,1160576736@qq.com。

    swj8013@126.com

    投稿平臺(tái):http://cykk.cbpt.cnki.net

    猜你喜歡
    白牡丹白茶菌門(mén)
    何冰歌
    收藏與投資(2022年9期)2022-10-13 06:55:32
    “白牡丹”
    繪本
    野生樹(shù)鼩與人工飼養(yǎng)樹(shù)鼩消化道不同部位微生物組成的比較研究
    《戲牡丹》中白牡丹的角色感想及情緒變化分析
    戲劇之家(2019年20期)2019-08-09 07:43:03
    饑餓與重?cái)z食對(duì)河蟹腸道菌群結(jié)構(gòu)的影響
    這只貓說(shuō)得好有道理
    意林(2019年10期)2019-06-03 16:37:12
    昆蟲(chóng)體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價(jià)值分析
    ANIMATING ‘IP’
    精品电影一区二区在线| 一个人看的www免费观看视频| 精品免费久久久久久久清纯| 午夜a级毛片| 全区人妻精品视频| 国产人伦9x9x在线观看| www.999成人在线观看| 国内精品美女久久久久久| 中文资源天堂在线| 级片在线观看| 国产单亲对白刺激| 99精品久久久久人妻精品| 一边摸一边抽搐一进一小说| 丰满人妻一区二区三区视频av | 亚洲av免费在线观看| 免费在线观看影片大全网站| 国内精品一区二区在线观看| av黄色大香蕉| 在线观看免费午夜福利视频| 欧美一区二区国产精品久久精品| av在线蜜桃| 精品久久蜜臀av无| 亚洲精品色激情综合| 啦啦啦观看免费观看视频高清| 久久久久久久精品吃奶| 国内精品一区二区在线观看| 国产亚洲精品久久久com| 神马国产精品三级电影在线观看| 久久久久性生活片| 美女黄网站色视频| 波多野结衣巨乳人妻| 嫩草影院精品99| 一个人看的www免费观看视频| 亚洲色图 男人天堂 中文字幕| 国产真人三级小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 黄色女人牲交| 一本综合久久免费| 午夜久久久久精精品| 免费在线观看亚洲国产| 成在线人永久免费视频| 久久精品国产亚洲av香蕉五月| 少妇人妻一区二区三区视频| 日韩人妻高清精品专区| 成年版毛片免费区| 久久久久免费精品人妻一区二区| 看免费av毛片| 精品99又大又爽又粗少妇毛片 | 巨乳人妻的诱惑在线观看| 最好的美女福利视频网| 午夜免费成人在线视频| 久久精品人妻少妇| 三级国产精品欧美在线观看 | 免费大片18禁| 亚洲中文av在线| 久久天堂一区二区三区四区| 黄色日韩在线| 国产1区2区3区精品| 久久国产精品人妻蜜桃| 日本在线视频免费播放| 久久午夜综合久久蜜桃| 日本三级黄在线观看| 成人三级做爰电影| 别揉我奶头~嗯~啊~动态视频| 久久这里只有精品19| 别揉我奶头~嗯~啊~动态视频| 女警被强在线播放| 国产日本99.免费观看| 操出白浆在线播放| av女优亚洲男人天堂 | 99国产精品99久久久久| 啦啦啦观看免费观看视频高清| 亚洲精品乱码久久久v下载方式 | 男插女下体视频免费在线播放| 搡老妇女老女人老熟妇| 观看美女的网站| 欧美性猛交╳xxx乱大交人| 男女床上黄色一级片免费看| 午夜a级毛片| 亚洲美女黄片视频| 日韩免费av在线播放| 午夜两性在线视频| 久久精品91蜜桃| 两个人的视频大全免费| 伊人久久大香线蕉亚洲五| 欧美一区二区精品小视频在线| 午夜福利在线观看免费完整高清在 | 精品福利观看| 三级男女做爰猛烈吃奶摸视频| h日本视频在线播放| 日日摸夜夜添夜夜添小说| 少妇的逼水好多| 一区二区三区国产精品乱码| 成人精品一区二区免费| 成人精品一区二区免费| 又黄又粗又硬又大视频| 少妇熟女aⅴ在线视频| 久久精品国产综合久久久| 亚洲无线观看免费| 亚洲欧美日韩东京热| 757午夜福利合集在线观看| 在线a可以看的网站| 欧美日韩精品网址| 国产亚洲av高清不卡| 亚洲成人久久性| 日韩欧美国产一区二区入口| 村上凉子中文字幕在线| 国产精品久久久av美女十八| 九九久久精品国产亚洲av麻豆 | 99久久精品热视频| 91麻豆精品激情在线观看国产| 欧美黑人欧美精品刺激| 一二三四社区在线视频社区8| 精品久久久久久久毛片微露脸| 香蕉丝袜av| 黄色视频,在线免费观看| a级毛片a级免费在线| 90打野战视频偷拍视频| 舔av片在线| 一夜夜www| 岛国在线免费视频观看| 色播亚洲综合网| 这个男人来自地球电影免费观看| 午夜福利在线观看吧| 色精品久久人妻99蜜桃| 久久性视频一级片| 国产欧美日韩一区二区精品| 99久久国产精品久久久| 亚洲欧美激情综合另类| 精品无人区乱码1区二区| 欧美丝袜亚洲另类 | 欧美黑人巨大hd| 久久久久久久精品吃奶| 在线观看午夜福利视频| av在线天堂中文字幕| 久久久久久国产a免费观看| 久久久久久大精品| 丁香欧美五月| 精品国产三级普通话版| 日韩欧美三级三区| 久久精品国产清高在天天线| 国产精品影院久久| 久久久久久人人人人人| 我的老师免费观看完整版| 每晚都被弄得嗷嗷叫到高潮| 久久国产精品影院| 九九在线视频观看精品| 亚洲熟女毛片儿| 亚洲男人的天堂狠狠| 精品久久久久久久久久久久久| 久久国产精品人妻蜜桃| 一区二区三区国产精品乱码| 国产亚洲精品一区二区www| 国产精品98久久久久久宅男小说| 夜夜看夜夜爽夜夜摸| 亚洲中文字幕日韩| 精品国产亚洲在线| 一区二区三区高清视频在线| 伊人久久大香线蕉亚洲五| 五月伊人婷婷丁香| 99久久精品热视频| 欧美黄色片欧美黄色片| 国产成人福利小说| 嫩草影院精品99| 深夜精品福利| www.999成人在线观看| 老汉色av国产亚洲站长工具| 视频区欧美日本亚洲| 俄罗斯特黄特色一大片| 国产成人精品久久二区二区91| 精品久久久久久久末码| or卡值多少钱| 久久久久久久午夜电影| 国内毛片毛片毛片毛片毛片| 欧美zozozo另类| 欧美成狂野欧美在线观看| 一级作爱视频免费观看| 在线免费观看不下载黄p国产 | 嫩草影院精品99| 琪琪午夜伦伦电影理论片6080| 男女之事视频高清在线观看| 男女视频在线观看网站免费| 午夜福利在线观看吧| 亚洲成人久久爱视频| 法律面前人人平等表现在哪些方面| 亚洲av成人精品一区久久| 超碰成人久久| 一二三四在线观看免费中文在| 99国产极品粉嫩在线观看| 日本五十路高清| 成人国产一区最新在线观看| ponron亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 国内揄拍国产精品人妻在线| 久久久水蜜桃国产精品网| 亚洲精品中文字幕一二三四区| 久久午夜亚洲精品久久| 免费看日本二区| 看免费av毛片| www日本在线高清视频| 12—13女人毛片做爰片一| 国产男靠女视频免费网站| 色吧在线观看| 视频区欧美日本亚洲| 中亚洲国语对白在线视频| 欧美3d第一页| 国产又黄又爽又无遮挡在线| 麻豆一二三区av精品| 成人国产综合亚洲| 国产极品精品免费视频能看的| 亚洲专区中文字幕在线| av欧美777| 久久久久久大精品| 国产成人影院久久av| 国产欧美日韩一区二区精品| 久久久久国产一级毛片高清牌| 中文字幕久久专区| 又黄又爽又免费观看的视频| 三级国产精品欧美在线观看 | 国产成人精品无人区| 久99久视频精品免费| 色综合欧美亚洲国产小说| 精品一区二区三区av网在线观看| 亚洲无线观看免费| 国产精品1区2区在线观看.| 欧美又色又爽又黄视频| 中文字幕人成人乱码亚洲影| АⅤ资源中文在线天堂| 亚洲欧美一区二区三区黑人| 麻豆成人午夜福利视频| 欧美日韩亚洲国产一区二区在线观看| 91老司机精品| 国产精品国产高清国产av| 国产一级毛片七仙女欲春2| 少妇人妻一区二区三区视频| 在线免费观看的www视频| 黄频高清免费视频| 1024手机看黄色片| 国产69精品久久久久777片 | 日韩欧美国产一区二区入口| 国产高清激情床上av| 欧美日韩中文字幕国产精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 操出白浆在线播放| 亚洲欧洲精品一区二区精品久久久| 91在线观看av| 国产精品久久久久久久电影 | 国产成人欧美在线观看| 亚洲 欧美一区二区三区| 小说图片视频综合网站| 嫩草影视91久久| 91麻豆精品激情在线观看国产| 中文字幕久久专区| 午夜精品一区二区三区免费看| 亚洲专区国产一区二区| 国产精品影院久久| 国产91精品成人一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美成狂野欧美在线观看| x7x7x7水蜜桃| 夜夜躁狠狠躁天天躁| 无限看片的www在线观看| 18禁观看日本| 色吧在线观看| 国产精品女同一区二区软件 | 亚洲一区二区三区不卡视频| 动漫黄色视频在线观看| 国产不卡一卡二| 国产精品女同一区二区软件 | 国产三级黄色录像| 日本一二三区视频观看| 亚洲人成伊人成综合网2020| 一本精品99久久精品77| 亚洲专区字幕在线| 久久久久久人人人人人| 宅男免费午夜| 国产精品久久久久久人妻精品电影| xxxwww97欧美| 亚洲在线自拍视频| 又黄又爽又免费观看的视频| 99国产极品粉嫩在线观看| 欧美性猛交黑人性爽| 男插女下体视频免费在线播放| 成熟少妇高潮喷水视频| 韩国av一区二区三区四区| 午夜视频精品福利| 他把我摸到了高潮在线观看| 欧美zozozo另类| 99久久精品国产亚洲精品| 成人一区二区视频在线观看| 欧美日韩综合久久久久久 | 精品一区二区三区四区五区乱码| 午夜久久久久精精品| 成人国产综合亚洲| 少妇人妻一区二区三区视频| 成在线人永久免费视频| 操出白浆在线播放| 天堂网av新在线| 精品一区二区三区视频在线观看免费| 国产精品香港三级国产av潘金莲| 国产精品自产拍在线观看55亚洲| 国产欧美日韩一区二区精品| 中文字幕最新亚洲高清| 国产在线精品亚洲第一网站| 老司机在亚洲福利影院| 男女下面进入的视频免费午夜| 免费高清视频大片| 亚洲欧美一区二区三区黑人| 狂野欧美白嫩少妇大欣赏| 久久热在线av| 全区人妻精品视频| ponron亚洲| 亚洲精品色激情综合| 国产亚洲精品久久久com| 精品久久久久久久久久久久久| 日本成人三级电影网站| 婷婷六月久久综合丁香| 夜夜爽天天搞| 国产精品永久免费网站| 老司机在亚洲福利影院| 一级a爱片免费观看的视频| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 国产成年人精品一区二区| 最近在线观看免费完整版| 国产aⅴ精品一区二区三区波| 熟女电影av网| 久久天堂一区二区三区四区| 欧美性猛交╳xxx乱大交人| 中文字幕精品亚洲无线码一区| 亚洲av中文字字幕乱码综合| 欧美日本亚洲视频在线播放| 少妇裸体淫交视频免费看高清| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美日韩东京热| 18禁国产床啪视频网站| 18禁黄网站禁片午夜丰满| 99热只有精品国产| 成人三级做爰电影| 欧美日韩一级在线毛片| 成人性生交大片免费视频hd| 听说在线观看完整版免费高清| 女生性感内裤真人,穿戴方法视频| 久久精品夜夜夜夜夜久久蜜豆| 久久久成人免费电影| 久久精品aⅴ一区二区三区四区| 美女免费视频网站| 无遮挡黄片免费观看| 国产精品,欧美在线| 欧美午夜高清在线| 美女午夜性视频免费| 国产男靠女视频免费网站| 2021天堂中文幕一二区在线观| av中文乱码字幕在线| 久久久国产欧美日韩av| 2021天堂中文幕一二区在线观| 男人舔女人下体高潮全视频| 美女高潮喷水抽搐中文字幕| or卡值多少钱| 亚洲国产欧美一区二区综合| 国产97色在线日韩免费| 久久午夜综合久久蜜桃| 一级a爱片免费观看的视频| 午夜视频精品福利| 一级a爱片免费观看的视频| 97碰自拍视频| 欧美日本亚洲视频在线播放| 亚洲午夜理论影院| 99久久99久久久精品蜜桃| 色老头精品视频在线观看| 国产精品一区二区精品视频观看| 白带黄色成豆腐渣| 久久久久免费精品人妻一区二区| 99久久99久久久精品蜜桃| 日本黄大片高清| 在线视频色国产色| 国产成年人精品一区二区| 午夜福利免费观看在线| 在线观看午夜福利视频| 国内精品久久久久久久电影| 成人亚洲精品av一区二区| 最近最新中文字幕大全电影3| 精品久久久久久久人妻蜜臀av| 亚洲中文字幕日韩| 人人妻人人澡欧美一区二区| 97超视频在线观看视频| 成人特级av手机在线观看| 亚洲国产精品sss在线观看| 亚洲国产精品sss在线观看| 一级a爱片免费观看的视频| 精品国产美女av久久久久小说| 日本精品一区二区三区蜜桃| 国产一区二区在线av高清观看| 亚洲国产精品成人综合色| 观看美女的网站| 久久久国产精品麻豆| 久久久水蜜桃国产精品网| 亚洲精品色激情综合| 欧美最黄视频在线播放免费| xxx96com| 老司机深夜福利视频在线观看| 精品福利观看| 黑人巨大精品欧美一区二区mp4| 亚洲成人中文字幕在线播放| 亚洲精品美女久久久久99蜜臀| a在线观看视频网站| 日本黄色片子视频| 男女那种视频在线观看| 小说图片视频综合网站| 黄色女人牲交| 国产成人精品无人区| 女人被狂操c到高潮| 久久精品国产综合久久久| 给我免费播放毛片高清在线观看| 欧美极品一区二区三区四区| 白带黄色成豆腐渣| 男女床上黄色一级片免费看| 69av精品久久久久久| 夜夜看夜夜爽夜夜摸| 美女高潮的动态| 老鸭窝网址在线观看| 亚洲av日韩精品久久久久久密| 中文字幕人妻丝袜一区二区| 久久久久久久精品吃奶| 欧美色视频一区免费| 一a级毛片在线观看| 亚洲国产精品成人综合色| 成熟少妇高潮喷水视频| 九色国产91popny在线| 亚洲成人免费电影在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品99久久99久久久不卡| 欧美成人性av电影在线观看| 757午夜福利合集在线观看| 搡老妇女老女人老熟妇| 成人国产综合亚洲| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 亚洲精品中文字幕一二三四区| 91在线观看av| 性色av乱码一区二区三区2| 欧美激情久久久久久爽电影| 亚洲精华国产精华精| 又紧又爽又黄一区二区| 青草久久国产| 此物有八面人人有两片| 欧美3d第一页| 婷婷六月久久综合丁香| 观看免费一级毛片| 国产午夜精品久久久久久| 在线国产一区二区在线| 中文字幕熟女人妻在线| 亚洲va日本ⅴa欧美va伊人久久| 一夜夜www| 国产精品国产高清国产av| 久久久久久久精品吃奶| 一级黄色大片毛片| 精品欧美国产一区二区三| 老汉色∧v一级毛片| 国产真实乱freesex| 91字幕亚洲| 麻豆成人午夜福利视频| 欧美又色又爽又黄视频| 久久精品国产99精品国产亚洲性色| 欧美在线黄色| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 久久久精品大字幕| 亚洲中文av在线| 少妇的逼水好多| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 成人国产综合亚洲| 久久九九热精品免费| 国产一级毛片七仙女欲春2| 国产精品美女特级片免费视频播放器 | av黄色大香蕉| 深夜精品福利| 国产激情偷乱视频一区二区| 午夜亚洲福利在线播放| 在线a可以看的网站| 亚洲五月婷婷丁香| 男人舔奶头视频| 最新在线观看一区二区三区| 婷婷六月久久综合丁香| 免费在线观看日本一区| 精品久久久久久久毛片微露脸| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 两性夫妻黄色片| 欧美又色又爽又黄视频| 校园春色视频在线观看| 久久中文看片网| 99国产精品99久久久久| 一夜夜www| 在线观看66精品国产| 国产淫片久久久久久久久 | 国产aⅴ精品一区二区三区波| 国产精品一区二区三区四区久久| 小蜜桃在线观看免费完整版高清| 国产成人精品久久二区二区91| 亚洲熟女毛片儿| 亚洲人成电影免费在线| 美女免费视频网站| 国产精品 欧美亚洲| 久久精品影院6| 黑人操中国人逼视频| 男女视频在线观看网站免费| 我的老师免费观看完整版| 美女 人体艺术 gogo| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 久久精品综合一区二区三区| 免费看光身美女| 欧美黑人巨大hd| 午夜福利在线观看免费完整高清在 | 亚洲,欧美精品.| 美女被艹到高潮喷水动态| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 99在线视频只有这里精品首页| 婷婷精品国产亚洲av在线| 中出人妻视频一区二区| 免费电影在线观看免费观看| 亚洲成av人片在线播放无| 欧美在线一区亚洲| 日韩欧美 国产精品| 亚洲第一欧美日韩一区二区三区| 国产私拍福利视频在线观看| 在线国产一区二区在线| 亚洲黑人精品在线| 欧美一区二区精品小视频在线| 狂野欧美激情性xxxx| 久久天堂一区二区三区四区| 色综合婷婷激情| 欧美在线一区亚洲| 久久久久久九九精品二区国产| 特大巨黑吊av在线直播| 欧美黑人欧美精品刺激| 欧美乱色亚洲激情| 99riav亚洲国产免费| 国产三级中文精品| 成人18禁在线播放| 成年女人永久免费观看视频| 精品久久久久久成人av| 一区二区三区激情视频| 国产精品av久久久久免费| 亚洲av电影在线进入| 桃红色精品国产亚洲av| 欧美乱码精品一区二区三区| 热99在线观看视频| 狂野欧美白嫩少妇大欣赏| 香蕉久久夜色| 亚洲精品色激情综合| 日韩有码中文字幕| 黄色丝袜av网址大全| av天堂中文字幕网| 欧美性猛交黑人性爽| 亚洲精品乱码久久久v下载方式 | 黄色丝袜av网址大全| 好男人电影高清在线观看| 欧美成人免费av一区二区三区| 一二三四社区在线视频社区8| 欧美乱色亚洲激情| 真人一进一出gif抽搐免费| 2021天堂中文幕一二区在线观| 女人被狂操c到高潮| 久久婷婷人人爽人人干人人爱| 国产1区2区3区精品| 日本 欧美在线| 久久精品人妻少妇| 精品久久久久久久久久久久久| 国产一区二区在线av高清观看| 免费看光身美女| 一级作爱视频免费观看| 日韩 欧美 亚洲 中文字幕| 99热精品在线国产| 搡老妇女老女人老熟妇| 免费在线观看影片大全网站| 99热这里只有精品一区 | 老司机福利观看| 久久精品影院6| 夜夜看夜夜爽夜夜摸| 美女午夜性视频免费| 亚洲精品色激情综合| 美女 人体艺术 gogo| 国产毛片a区久久久久| 久久天堂一区二区三区四区| 成年人黄色毛片网站| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 99国产精品一区二区蜜桃av| 可以在线观看毛片的网站| 国产又色又爽无遮挡免费看| 成人性生交大片免费视频hd| 亚洲午夜精品一区,二区,三区| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 淫秽高清视频在线观看| 日本成人三级电影网站| 国产69精品久久久久777片 | 亚洲精品乱码久久久v下载方式 | 99热6这里只有精品| 1024香蕉在线观看| 一本综合久久免费| 麻豆国产av国片精品| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久人妻蜜臀av| 18美女黄网站色大片免费观看| 伊人久久大香线蕉亚洲五| 国产激情久久老熟女|