• <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’
    啦啦啦韩国在线观看视频| 日韩av在线大香蕉| 五月伊人婷婷丁香| 午夜福利在线观看免费完整高清在 | 日韩欧美一区二区三区在线观看| 日本-黄色视频高清免费观看| 成人欧美大片| 久久久精品大字幕| 精品人妻熟女av久视频| 成人av在线播放网站| 麻豆成人午夜福利视频| 最近中文字幕高清免费大全6 | 免费av不卡在线播放| av在线老鸭窝| 国产在线精品亚洲第一网站| 国产精品伦人一区二区| 性欧美人与动物交配| 一本精品99久久精品77| 午夜精品一区二区三区免费看| 中文字幕av在线有码专区| 伦精品一区二区三区| 免费av不卡在线播放| 三级国产精品欧美在线观看| 成年女人看的毛片在线观看| 国产一区二区三区视频了| 欧美激情久久久久久爽电影| 伊人久久精品亚洲午夜| 熟女电影av网| 日本黄色视频三级网站网址| 亚洲一区二区三区色噜噜| 久久欧美精品欧美久久欧美| 日韩一本色道免费dvd| 国产伦精品一区二区三区视频9| 中文在线观看免费www的网站| 国产精品久久久久久久久免| 老司机福利观看| 欧美成人免费av一区二区三区| 亚洲图色成人| 美女被艹到高潮喷水动态| 亚洲欧美日韩东京热| 国产视频一区二区在线看| 亚洲精品在线观看二区| 国产主播在线观看一区二区| 亚洲熟妇熟女久久| 亚洲内射少妇av| 国产伦精品一区二区三区视频9| 亚洲av.av天堂| 免费电影在线观看免费观看| 亚洲av.av天堂| 麻豆国产97在线/欧美| 精品免费久久久久久久清纯| 欧美一区二区亚洲| 免费观看人在逋| 国产精品爽爽va在线观看网站| 十八禁网站免费在线| 国产男人的电影天堂91| 在线观看美女被高潮喷水网站| 一本久久中文字幕| 很黄的视频免费| 精品久久久久久久久久久久久| 精品久久久噜噜| 国产伦在线观看视频一区| 婷婷亚洲欧美| 小说图片视频综合网站| 欧美最新免费一区二区三区| 我的老师免费观看完整版| 国产午夜福利久久久久久| 国产乱人视频| 午夜福利欧美成人| 国产高清视频在线观看网站| 亚洲七黄色美女视频| 国产欧美日韩一区二区精品| 男人狂女人下面高潮的视频| 久久草成人影院| 成熟少妇高潮喷水视频| 免费黄网站久久成人精品| 国产在线精品亚洲第一网站| 国产日本99.免费观看| 欧美成人免费av一区二区三区| 淫妇啪啪啪对白视频| 搡老熟女国产l中国老女人| 国产精品乱码一区二三区的特点| 黄色配什么色好看| 国产探花在线观看一区二区| 国产精品野战在线观看| 亚洲最大成人中文| 老熟妇乱子伦视频在线观看| 少妇高潮的动态图| 日韩精品中文字幕看吧| 波野结衣二区三区在线| 亚洲国产精品合色在线| 内射极品少妇av片p| 日韩精品中文字幕看吧| 国产精品综合久久久久久久免费| 成人三级黄色视频| 国产不卡一卡二| 成人特级av手机在线观看| 国产精品爽爽va在线观看网站| 天美传媒精品一区二区| 校园人妻丝袜中文字幕| 最近中文字幕高清免费大全6 | 午夜激情福利司机影院| 久久久国产成人精品二区| 99热这里只有是精品50| 老熟妇乱子伦视频在线观看| 午夜免费成人在线视频| 97超视频在线观看视频| 少妇裸体淫交视频免费看高清| 性插视频无遮挡在线免费观看| 别揉我奶头~嗯~啊~动态视频| 午夜激情欧美在线| 久99久视频精品免费| 亚洲精品456在线播放app | 国产一区二区在线av高清观看| 国产精品久久视频播放| 日韩亚洲欧美综合| 久久精品国产99精品国产亚洲性色| 久久久久久久亚洲中文字幕| 国产探花极品一区二区| 日韩欧美国产一区二区入口| 日韩中字成人| 欧美一区二区国产精品久久精品| 床上黄色一级片| 69人妻影院| 日本欧美国产在线视频| 免费在线观看成人毛片| 国产精品伦人一区二区| 亚洲天堂国产精品一区在线| 国内精品久久久久精免费| 黄色日韩在线| 精品久久久久久,| 精品久久国产蜜桃| 男插女下体视频免费在线播放| 成年女人毛片免费观看观看9| 久久久久久久亚洲中文字幕| 女同久久另类99精品国产91| 亚洲自拍偷在线| 淫秽高清视频在线观看| 国产成人影院久久av| 国产日本99.免费观看| 如何舔出高潮| 亚洲在线自拍视频| 亚洲午夜理论影院| 成人毛片a级毛片在线播放| 97人妻精品一区二区三区麻豆| 久久精品影院6| 美女cb高潮喷水在线观看| 久久久久精品国产欧美久久久| 成人高潮视频无遮挡免费网站| 男人舔奶头视频| 99热只有精品国产| 男女边吃奶边做爰视频| 嫩草影院入口| 在线免费观看不下载黄p国产 | 精品一区二区三区人妻视频| 精华霜和精华液先用哪个| 淫秽高清视频在线观看| 久久6这里有精品| 极品教师在线免费播放| 日韩高清综合在线| 三级毛片av免费| 日韩在线高清观看一区二区三区 | 亚洲成人免费电影在线观看| 久久国产乱子免费精品| 五月玫瑰六月丁香| 97热精品久久久久久| 熟女人妻精品中文字幕| 观看免费一级毛片| 色综合婷婷激情| 国产乱人伦免费视频| 日本色播在线视频| 亚洲av日韩精品久久久久久密| 日韩高清综合在线| 高清在线国产一区| 亚洲欧美日韩东京热| 国产乱人视频| 日本色播在线视频| 美女黄网站色视频| 性插视频无遮挡在线免费观看| 久99久视频精品免费| 久久精品国产99精品国产亚洲性色| 国产av一区在线观看免费| a级一级毛片免费在线观看| 国产白丝娇喘喷水9色精品| 日本五十路高清| 国产精品一区二区三区四区久久| 欧美日韩乱码在线| 性色avwww在线观看| 久久久国产成人精品二区| 日韩欧美在线乱码| 国产午夜福利久久久久久| 国产成年人精品一区二区| 最近最新免费中文字幕在线| 欧美极品一区二区三区四区| 伦精品一区二区三区| 国产高清视频在线播放一区| 97超视频在线观看视频| 亚洲av一区综合| 国产淫片久久久久久久久| 午夜久久久久精精品| 国产成人一区二区在线| 三级国产精品欧美在线观看| 国产精品一及| a在线观看视频网站| 校园春色视频在线观看| 特级一级黄色大片| 无遮挡黄片免费观看| 国国产精品蜜臀av免费| 午夜精品在线福利| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久久久久免费视频| 日日摸夜夜添夜夜添av毛片 | 天堂√8在线中文| 嫩草影院入口| 尤物成人国产欧美一区二区三区| 麻豆成人av在线观看| 人人妻人人澡欧美一区二区| 露出奶头的视频| 美女高潮的动态| 天天一区二区日本电影三级| 午夜福利欧美成人| 91在线精品国自产拍蜜月| 国产人妻一区二区三区在| 欧美3d第一页| 91狼人影院| 在线观看66精品国产| 欧美bdsm另类| 免费人成视频x8x8入口观看| 精品99又大又爽又粗少妇毛片 | 岛国在线免费视频观看| 日韩欧美在线二视频| 变态另类成人亚洲欧美熟女| 国产精品国产三级国产av玫瑰| 18禁裸乳无遮挡免费网站照片| 欧洲精品卡2卡3卡4卡5卡区| 麻豆一二三区av精品| 麻豆成人av在线观看| 在线观看av片永久免费下载| 国产私拍福利视频在线观看| 久久久久九九精品影院| 亚洲三级黄色毛片| 99久久久亚洲精品蜜臀av| 国产精品日韩av在线免费观看| 桃色一区二区三区在线观看| 欧美日韩国产亚洲二区| 少妇人妻一区二区三区视频| 女人十人毛片免费观看3o分钟| 一卡2卡三卡四卡精品乱码亚洲| 一个人看的www免费观看视频| 波多野结衣高清作品| 成人av在线播放网站| 18+在线观看网站| 18禁黄网站禁片免费观看直播| 亚洲黑人精品在线| 亚洲性久久影院| 国产精品一区二区三区四区久久| 又爽又黄a免费视频| 村上凉子中文字幕在线| 国产成人aa在线观看| 国产精品亚洲一级av第二区| 欧美日韩乱码在线| 搡老熟女国产l中国老女人| 五月玫瑰六月丁香| 久久精品国产亚洲av天美| 午夜影院日韩av| 午夜免费激情av| 伊人久久精品亚洲午夜| 免费看美女性在线毛片视频| 免费av观看视频| 亚洲18禁久久av| 免费黄网站久久成人精品| 国产av在哪里看| 最近在线观看免费完整版| 波野结衣二区三区在线| 日韩av在线大香蕉| 人妻制服诱惑在线中文字幕| 少妇裸体淫交视频免费看高清| 国产欧美日韩精品一区二区| 久久精品国产自在天天线| 国内精品宾馆在线| 天堂影院成人在线观看| 免费观看在线日韩| 一级a爱片免费观看的视频| 亚洲精品国产成人久久av| 久久99热6这里只有精品| 日韩精品中文字幕看吧| 日韩欧美三级三区| 欧美xxxx性猛交bbbb| 日韩国内少妇激情av| 国产高清激情床上av| 欧美色视频一区免费| 91在线观看av| 男人舔女人下体高潮全视频| 高清毛片免费观看视频网站| 欧美潮喷喷水| 村上凉子中文字幕在线| 性插视频无遮挡在线免费观看| 成人综合一区亚洲| av天堂在线播放| 国产精品久久久久久av不卡| 联通29元200g的流量卡| 色在线成人网| 熟女人妻精品中文字幕| 22中文网久久字幕| 可以在线观看的亚洲视频| 在线天堂最新版资源| 一个人观看的视频www高清免费观看| 国产伦一二天堂av在线观看| 日韩强制内射视频| 成人欧美大片| 久久久久久国产a免费观看| 真实男女啪啪啪动态图| 99热只有精品国产| 老师上课跳d突然被开到最大视频| 一a级毛片在线观看| 色视频www国产| 久久久久久九九精品二区国产| 国产精华一区二区三区| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 亚洲国产精品合色在线| 久久久久久久精品吃奶| 欧美激情国产日韩精品一区| 国产伦人伦偷精品视频| 亚洲最大成人中文| 欧美激情在线99| 熟妇人妻久久中文字幕3abv| 制服丝袜大香蕉在线| 中文字幕高清在线视频| 舔av片在线| 精品久久久久久久人妻蜜臀av| 欧美精品国产亚洲| 亚洲欧美日韩高清在线视频| 天堂√8在线中文| 亚洲人成网站高清观看| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 国内久久婷婷六月综合欲色啪| 一a级毛片在线观看| 91精品国产九色| 国产一区二区三区av在线 | 中亚洲国语对白在线视频| 婷婷亚洲欧美| 亚洲五月天丁香| 老熟妇乱子伦视频在线观看| 中亚洲国语对白在线视频| 欧美不卡视频在线免费观看| 日本 欧美在线| 看十八女毛片水多多多| 最近最新免费中文字幕在线| 国产精品亚洲一级av第二区| 99riav亚洲国产免费| 国产精品亚洲一级av第二区| 久久九九热精品免费| 亚洲国产高清在线一区二区三| a在线观看视频网站| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 日本黄色视频三级网站网址| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线观看免费| 成人无遮挡网站| 一个人看的www免费观看视频| 色播亚洲综合网| 伦理电影大哥的女人| а√天堂www在线а√下载| 国产免费男女视频| 91久久精品国产一区二区成人| 国产毛片a区久久久久| 一区二区三区四区激情视频 | 午夜日韩欧美国产| 亚洲无线在线观看| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| 免费观看人在逋| 日韩国内少妇激情av| а√天堂www在线а√下载| 少妇熟女aⅴ在线视频| 一区二区三区四区激情视频 | 亚洲无线观看免费| 欧美一区二区亚洲| 免费看光身美女| 欧美成人性av电影在线观看| 波野结衣二区三区在线| 美女cb高潮喷水在线观看| 国产午夜福利久久久久久| 婷婷丁香在线五月| 精品久久久久久久久亚洲 | 国产精品一区二区三区四区久久| 成人欧美大片| 久久久精品欧美日韩精品| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 熟妇人妻久久中文字幕3abv| 一个人观看的视频www高清免费观看| 成年版毛片免费区| 在线观看午夜福利视频| 国产精品亚洲一级av第二区| 欧美色视频一区免费| 国产爱豆传媒在线观看| aaaaa片日本免费| 亚洲成人精品中文字幕电影| 欧美日韩亚洲国产一区二区在线观看| 国产精品av视频在线免费观看| 亚洲国产精品久久男人天堂| 美女高潮喷水抽搐中文字幕| 麻豆一二三区av精品| 亚洲狠狠婷婷综合久久图片| 不卡视频在线观看欧美| 日本 欧美在线| 九九在线视频观看精品| 日韩亚洲欧美综合| 1000部很黄的大片| 在线看三级毛片| 美女大奶头视频| 欧美极品一区二区三区四区| 国产蜜桃级精品一区二区三区| 久久国产精品人妻蜜桃| 亚洲熟妇熟女久久| 色5月婷婷丁香| 校园人妻丝袜中文字幕| 亚洲专区国产一区二区| 男人和女人高潮做爰伦理| 精品人妻偷拍中文字幕| 国产亚洲欧美98| 欧美日韩精品成人综合77777| 亚洲精品乱码久久久v下载方式| 精品一区二区免费观看| 国产探花在线观看一区二区| 国产精品自产拍在线观看55亚洲| netflix在线观看网站| 日本精品一区二区三区蜜桃| 内地一区二区视频在线| av中文乱码字幕在线| 嫁个100分男人电影在线观看| 国产伦一二天堂av在线观看| 国产精品爽爽va在线观看网站| 真实男女啪啪啪动态图| 欧美日本亚洲视频在线播放| 欧美xxxx性猛交bbbb| 国内精品美女久久久久久| 嫩草影院入口| 亚洲成a人片在线一区二区| 国产精品一区二区性色av| 午夜老司机福利剧场| av在线亚洲专区| 琪琪午夜伦伦电影理论片6080| 中文在线观看免费www的网站| 国产一区二区在线观看日韩| 可以在线观看的亚洲视频| 成年女人看的毛片在线观看| 99久久无色码亚洲精品果冻| 国产精品99久久久久久久久| 国产精品爽爽va在线观看网站| 成人综合一区亚洲| 国产精品国产三级国产av玫瑰| 直男gayav资源| 成人国产麻豆网| 特大巨黑吊av在线直播| 久久欧美精品欧美久久欧美| 午夜精品久久久久久毛片777| 日本黄色视频三级网站网址| 一区二区三区免费毛片| 白带黄色成豆腐渣| 午夜亚洲福利在线播放| 91在线观看av| a级毛片免费高清观看在线播放| 人妻丰满熟妇av一区二区三区| 丰满人妻一区二区三区视频av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人中文字幕在线播放| 久久草成人影院| 无遮挡黄片免费观看| 久久精品国产亚洲av天美| 十八禁网站免费在线| 日日摸夜夜添夜夜添av毛片 | 国产精品98久久久久久宅男小说| 简卡轻食公司| 日韩中字成人| 国内精品久久久久精免费| 亚洲性夜色夜夜综合| 国语自产精品视频在线第100页| 成年人黄色毛片网站| 婷婷精品国产亚洲av在线| 亚洲人成网站在线播| 88av欧美| 欧美+日韩+精品| www.www免费av| 成人国产一区最新在线观看| 校园春色视频在线观看| 国产av不卡久久| 熟妇人妻久久中文字幕3abv| 春色校园在线视频观看| 99国产极品粉嫩在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品久久国产蜜桃| 黄色一级大片看看| 1024手机看黄色片| 少妇人妻一区二区三区视频| 亚洲人成网站高清观看| 18禁在线播放成人免费| 国产高清不卡午夜福利| 国产中年淑女户外野战色| 精品国产三级普通话版| 99久久精品一区二区三区| 亚洲七黄色美女视频| 国产极品精品免费视频能看的| 免费在线观看成人毛片| 一进一出好大好爽视频| 午夜日韩欧美国产| videossex国产| 91午夜精品亚洲一区二区三区 | 国产国拍精品亚洲av在线观看| 欧美绝顶高潮抽搐喷水| 久久午夜亚洲精品久久| 亚洲精品日韩av片在线观看| 精品久久久噜噜| 搡老妇女老女人老熟妇| 国产精品一区二区性色av| 免费看日本二区| 精品午夜福利在线看| 亚洲,欧美,日韩| 午夜免费激情av| 2021天堂中文幕一二区在线观| 午夜久久久久精精品| 亚洲精品国产成人久久av| 亚洲av二区三区四区| 欧美xxxx性猛交bbbb| 国产精品av视频在线免费观看| 看片在线看免费视频| 久久久久久久精品吃奶| 亚洲专区中文字幕在线| 别揉我奶头 嗯啊视频| 国产一区二区激情短视频| 欧美日韩瑟瑟在线播放| 国产人妻一区二区三区在| 久久精品久久久久久噜噜老黄 | 日韩av在线大香蕉| 在线观看av片永久免费下载| 欧美国产日韩亚洲一区| 久久精品综合一区二区三区| 亚洲专区国产一区二区| 韩国av在线不卡| 国产91精品成人一区二区三区| 亚洲精华国产精华液的使用体验 | 久久6这里有精品| 日韩欧美 国产精品| 亚洲男人的天堂狠狠| 欧美日韩中文字幕国产精品一区二区三区| 婷婷精品国产亚洲av| 热99re8久久精品国产| 国产精品久久久久久精品电影| 亚洲成av人片在线播放无| 中国美女看黄片| 亚洲av中文av极速乱 | a级毛片a级免费在线| 亚洲专区国产一区二区| 国产亚洲欧美98| 五月玫瑰六月丁香| 又黄又爽又免费观看的视频| 国产美女午夜福利| 一进一出好大好爽视频| 欧美+日韩+精品| 麻豆国产av国片精品| 一区福利在线观看| 国产亚洲精品综合一区在线观看| 男人的好看免费观看在线视频| 少妇人妻一区二区三区视频| 日日夜夜操网爽| 男女视频在线观看网站免费| 午夜免费成人在线视频| 真实男女啪啪啪动态图| 国产高清有码在线观看视频| 国产乱人视频| 国产aⅴ精品一区二区三区波| 国产伦人伦偷精品视频| 欧美一区二区精品小视频在线| 日本在线视频免费播放| 99精品久久久久人妻精品| 色噜噜av男人的天堂激情| 看黄色毛片网站| 午夜福利视频1000在线观看| 精品久久久久久久久亚洲 | 免费搜索国产男女视频| 亚洲精品在线观看二区| 久久久精品大字幕| 一个人观看的视频www高清免费观看| 波多野结衣高清无吗| 又黄又爽又刺激的免费视频.| 男女之事视频高清在线观看| 精品人妻偷拍中文字幕| 国产大屁股一区二区在线视频| 丝袜美腿在线中文| 国产精品久久久久久亚洲av鲁大| 中文在线观看免费www的网站| 日本a在线网址| 91午夜精品亚洲一区二区三区 | 99九九线精品视频在线观看视频| 精品一区二区三区人妻视频| 亚洲精品乱码久久久v下载方式| 美女xxoo啪啪120秒动态图| 两个人的视频大全免费| 国产av一区在线观看免费| 国产伦人伦偷精品视频| 男女啪啪激烈高潮av片| 久久人妻av系列| 99国产精品一区二区蜜桃av| 亚洲av中文av极速乱 | 日韩欧美免费精品| 91久久精品电影网| 亚洲天堂国产精品一区在线| 一区二区三区四区激情视频 | 国产精品久久久久久久电影|