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

    基于GC-MS技術(shù)對(duì)不同產(chǎn)地稻米的代謝組學(xué)分析

    2019-07-05 02:13:14富天昕馮玉超張麗媛王長(zhǎng)遠(yuǎn)
    食品科學(xué) 2019年12期
    關(guān)鍵詞:差異

    富天昕,馮玉超,張麗媛,2,李 雪,王長(zhǎng)遠(yuǎn),2,*

    (黑龍江八一農(nóng)墾大學(xué)食品學(xué)院,黑龍江省農(nóng)產(chǎn)品加工與質(zhì)量安全重點(diǎn)實(shí)驗(yàn)室,黑龍江 大慶 163319)

    代謝組學(xué)這一概念由Nicholson在2002年提出,是繼基因組學(xué)和蛋白質(zhì)組學(xué)之后新發(fā)展起來(lái)的一門學(xué)科,是系統(tǒng)生物學(xué)的重要組成部分[1]。概念提出后,該領(lǐng)域的研究迅速發(fā)展并滲透到多個(gè)領(lǐng)域[2-3]。植物代謝組學(xué)是代謝組學(xué)的重要分支,它主要研究的是作為各種代謝路徑的底物和產(chǎn)物的小分子代謝物(mw<1 000 Da),結(jié)合基因、蛋白質(zhì)、表型等綜合分析植物的生理生化過(guò)程[4-5]。而植物體的最終表現(xiàn)方式則由代謝產(chǎn)物體現(xiàn),代謝產(chǎn)物在生物體信息傳遞鏈的末端,是植物代謝的最終產(chǎn)物[6]。

    近年來(lái)植物代謝得到迅速發(fā)展,而稻米作為主要糧食作物之一,也成為研究的熱點(diǎn)[7-11]。稻米中富含豐富的蛋白質(zhì)、糖類、脂肪和膳食纖維等營(yíng)養(yǎng)成分且氨基酸組成比較完全,有助于人體的吸收[12]。稻米的品質(zhì)不僅與品種有關(guān),還與其生長(zhǎng)環(huán)境密切相關(guān)。不同產(chǎn)地的稻米因其種植的土壤、氣候、水源不同,其營(yíng)養(yǎng)品質(zhì)和代謝產(chǎn)物的種類、含量也不盡相同[13]。代謝組學(xué)的研究方法,可以從系統(tǒng)生物學(xué)的角度解釋代謝物變化與生命活動(dòng)現(xiàn)象的關(guān)系,近幾年代謝組學(xué)在農(nóng)產(chǎn)品鑒別與產(chǎn)地區(qū)分方向應(yīng)用廣泛,為農(nóng)產(chǎn)品的溯源提供了新的途徑[14-18]。氣相色譜-質(zhì)譜(gas chromatography-mass spectrometry,GC-MS)聯(lián)用技術(shù)具有重復(fù)性好、靈敏度高等特點(diǎn),是目前代謝組學(xué)研究中廣泛運(yùn)用的分析技術(shù)并且有很長(zhǎng)的研究歷史[18-22]。田福林等[23]利用GC-MS技術(shù)對(duì)4 種不同產(chǎn)地大米中的脂溶性成分進(jìn)行分析,結(jié)合主成分分析(principal component analysis,PCA)和聚類分析,建立一種對(duì)不同產(chǎn)地大米進(jìn)行快速鑒定的方法。本實(shí)驗(yàn)利用GC-MS聯(lián)用技術(shù)對(duì)黑龍江省建三江地區(qū)和響水地區(qū)的稻米進(jìn)行差異代謝物的研究,探討產(chǎn)地對(duì)水稻代謝物的影響,并從代謝物的角度著手,篩選對(duì)區(qū)分2 個(gè)產(chǎn)地具有重要作用的差異代謝物,為通過(guò)代謝物對(duì)稻米進(jìn)行產(chǎn)地溯源提供理論依據(jù)。

    1 材料與方法

    1.1 材料與試劑

    水稻采自黑龍江省建三江和響水地理標(biāo)志稻米保護(hù)區(qū),建三江采樣地點(diǎn)選擇為12 個(gè)農(nóng)場(chǎng)布點(diǎn);響水采樣地點(diǎn)選擇5 個(gè)鄉(xiāng)鎮(zhèn)布點(diǎn),均選擇主產(chǎn)縣、主產(chǎn)鄉(xiāng)(鎮(zhèn)或管理區(qū))、主產(chǎn)村(屯)的大面積種植地塊進(jìn)行采樣,采樣點(diǎn)依照種植范圍大小設(shè)置,每塊地塊隨機(jī)設(shè)置分點(diǎn)5 個(gè)重復(fù),采樣地點(diǎn)按照東、西、南、北、中5 個(gè)區(qū)域設(shè)計(jì),依據(jù)代表性采樣原則,每個(gè)采樣點(diǎn)按照不同方位采集1~2 kg稻穗,每個(gè)產(chǎn)區(qū)選取30 個(gè)樣本。響水地區(qū)水稻品種為稻花香,建三江地區(qū)品種主要為龍粳系列,龍粳26、龍粳29、龍粳31、龍粳39、龍粳40、龍粳43、龍粳46,以及綏粳4、松粳9等。

    吡啶(純度≥99.9%)、甲醇(均為色譜純) 美國(guó)阿拉丁試劑有限公司;2-氯苯丙氨酸(純度98.5%)、甲氧基胺鹽酸鹽(純度98%)、99% N,O-雙(三甲基硅)三氟乙酰胺(含1%三甲基氯硅烷) 上海麥克林試劑有限公司。

    1.2 儀器與設(shè)備

    7890A/5975C GC-MS聯(lián)用儀、HP-5ms色譜柱(30 m×0.25 mm,0.25 μm) 美國(guó)Agilent公司;QP2010 GC-MS聯(lián)用儀 日本島津公司;DHG-9123A型電熱恒溫鼓風(fēng)干燥箱 上海精宏實(shí)驗(yàn)設(shè)備有限公司;MTN-2800D氮吹濃縮裝置 天津奧特賽恩斯儀器有限公司;恒溫箱 上海森信實(shí)驗(yàn)儀器有限公司;TGL-16B臺(tái)式離心機(jī) 上海安亭科學(xué)儀器廠。

    1.3 方法

    樣品制備根據(jù)王玲[24]的方法加以改進(jìn)。稻種子液氮研磨后稱取50 mg粉末于2 mL的EP管中,加入800 μL甲醇和10 μL內(nèi)標(biāo)(2-氯苯丙氨酸),渦旋混勻30 s,置于4 ℃離心機(jī)中,12 000 r/min離心15 min,吸取200 μL上清液,轉(zhuǎn)入進(jìn)樣小瓶中氮?dú)獯蹈伞Q苌喝?0 μL甲氧氨鹽酸吡啶溶液至濃縮后的樣品中,渦旋混勻至完全溶解,置于37 ℃恒溫箱90 min,取出后加入30 μL N,O-雙(三甲基硅)三氟乙酰胺,置于70 ℃烘箱60 min,取出后待檢測(cè)。

    1.4 數(shù)據(jù)處理

    數(shù)據(jù)的特征峰提取和預(yù)處理在R軟件平臺(tái)采用XCMS軟件包條件下進(jìn)行,將處理的數(shù)據(jù)矩陣導(dǎo)入SIMCA-P軟件,進(jìn)行PCA[24]、正交偏最小二乘-判別分析(orthogonal partial least squares-discriminant analysis,OPLS-DA)等多元統(tǒng)計(jì)分析,再根據(jù)變量對(duì)分組貢獻(xiàn)值(variable importance in the projection,VIP)的大小和組間變化的顯著性(P<0.05)進(jìn)行差異性代謝產(chǎn)物的篩選。差異性代謝物的定性方法為:根據(jù)保留時(shí)間及質(zhì)荷比(m/z)在NIST商業(yè)數(shù)據(jù)庫(kù)和Wiley Registry代謝組數(shù)據(jù)庫(kù)進(jìn)行注釋。

    2 結(jié)果與分析

    2.1 GC-MS分析

    2.1.1 典型離子流色譜分析

    由圖1可以看出,色譜峰基線平穩(wěn),峰分離效果良好,儀器穩(wěn)定性良好。

    圖1 典型樣本的總離子流色譜圖Fig. 1 Typical total ion chromatogram of rice metabolites

    2.1.2 代謝組分的質(zhì)譜分析

    通過(guò)對(duì)GC-MS原始數(shù)據(jù)的處理,兩產(chǎn)地60 個(gè)稻米種子樣品中共掃描到143 個(gè)峰,根據(jù)保留時(shí)間及m/z在NIST商業(yè)數(shù)據(jù)庫(kù)和Wiley Registry代謝組數(shù)據(jù)庫(kù)進(jìn)行注釋,共注釋到39 個(gè)代謝產(chǎn)物,其中建三江地區(qū)有29 個(gè)代謝產(chǎn)物,響水地區(qū)有26 個(gè)代謝產(chǎn)物,如表1所示。

    表1中1~16號(hào)為建三江地區(qū)和響水地區(qū)共同存在的代謝物,17~29號(hào)為建三江地區(qū)稻米樣品相對(duì)于響水地區(qū)樣品具有的代謝物,30~39號(hào)為響水地區(qū)稻米樣品相對(duì)于建三江地區(qū)樣品具有的代謝物,除共同存在的代謝物外,2 個(gè)產(chǎn)地單獨(dú)具有的代謝物并不一定是相應(yīng)產(chǎn)地特有的代謝物,可能還與稻米樣品的品種有關(guān),由此可知,產(chǎn)地不同,稻米的代謝物的數(shù)量不也相同。如表2所示,建三江地區(qū)與響水地區(qū)相比,多元醇、有機(jī)酸類物質(zhì)含量要高于響水地區(qū),而響水地區(qū)的糖類、脂肪酸、氨基酸和其他類物質(zhì)要高于建三江地區(qū),核苷酸類物質(zhì)含量相差不多,可見(jiàn)產(chǎn)地不同對(duì)稻米代謝物的種類也會(huì)有影響。

    表2 建三江地區(qū)和響水地區(qū)稻米代謝物分類Table 2 Classification of metabolites in rice samples from Jiansanjiang and Xiangshui%

    2.2 兩產(chǎn)地稻米的代謝組學(xué)差異分析

    2.2.1 PCA結(jié)果

    圖2 2 種樣品的PCA得分圖Fig. 2 PCA scores plot of two rice samples

    為判別建三江地區(qū)和響水地區(qū)稻米樣品間的差異,首先采用PCA方法進(jìn)行分析,共得到3 個(gè)主成分,累計(jì)R2X=0.675,Q2=0.52,2 個(gè)參數(shù)數(shù)值均大于0.5,且二者數(shù)值相差不大,說(shuō)明PCA效果較好。如圖2所示,兩產(chǎn)地稻米樣品除個(gè)別異常樣品外絕大部分的樣品均處于在置信區(qū)間內(nèi),但兩產(chǎn)地樣品并未完全分開(kāi),還存在個(gè)別的重疊樣品,說(shuō)明這些樣本中所含有的分子的組成和濃度越接近,而分子組成和濃度接近可能是相同品種間的聚集,也有可能是不同品種但是組分含量上相近。PCA這種無(wú)監(jiān)督分析方法不能忽略組內(nèi)誤差、消除與研究目的無(wú)關(guān)的隨機(jī)誤差,所以不利于發(fā)現(xiàn)組間差異,之后將采用有監(jiān)督的方法進(jìn)行進(jìn)一步研究。

    2.2.2 OPLS-DA結(jié)果

    通過(guò)OPLS-DA可以過(guò)濾掉代謝物中,與分類變量不相關(guān)的正交變量,并對(duì)非正交變量和正交變量分別分析,從而獲取更加可靠的代謝物的組間差異與實(shí)驗(yàn)組的相關(guān)程度信息[26]。如圖3所示,建三江地區(qū)的稻米樣品主要分布在置信區(qū)間的左側(cè),響水地區(qū)的稻米樣品主要分布在置信區(qū)間的右側(cè),兩產(chǎn)地稻米樣品區(qū)分效果很好。本次OPLS-DA得到2 個(gè)主成分,R2X=0.575,R2Y=0.63,Q2=0.479,3 個(gè)參數(shù)數(shù)值接近或大于0.5,且相差不多,說(shuō)明OPLS-DA模型可解釋度較好,相比于PCA模型效果更好。

    圖3 兩產(chǎn)地樣本OPLS-DA得分圖Fig. 3 OPLS-DA score plot of two rice samples

    2.2.3 OPLS-DA置換檢驗(yàn)結(jié)果

    置換檢驗(yàn)通過(guò)隨機(jī)改變分類變量Y的排列順序,多次建立對(duì)應(yīng)的OPLS-DA模型以獲取隨機(jī)模型的R2和Q2值,在避免檢驗(yàn)?zāi)P偷倪^(guò)擬合以及評(píng)估模型的統(tǒng)計(jì)顯著性上有重要作用。對(duì)OPLS-DA模型的置換檢驗(yàn)結(jié)果如4圖所示。

    圖4 建三江地區(qū)與響水地區(qū)稻米樣品OPLS-DA置換檢驗(yàn)結(jié)果圖Fig. 4 OPLS-DA replacement test results of rice samples collected from Jiansanjiang and Xiangshui

    由圖4可知,所有R2點(diǎn)從左到右均低于最右的原始的R2點(diǎn),所有Q2點(diǎn)均低于右側(cè)的原始Q2點(diǎn),且點(diǎn)的回歸線與橫坐標(biāo)交叉或者小于0,說(shuō)明OPLS-DA模型是可靠的,可以較好地解釋兩組樣本之間的差異。同時(shí)隨著置換保留度逐漸降低,置換的Y變量比例增大,隨機(jī)模型的R2和Q2均逐漸下降,說(shuō)明原模型不存在過(guò)擬合現(xiàn)象,模型穩(wěn)健性較好。

    2.2.4 兩地區(qū)稻米樣品的差異代謝物鑒定

    采用OPLS-DA模型的VIP值(閾值>1),并結(jié)合學(xué)生氏t檢驗(yàn)(t-test)的P值(閾值<0.05)來(lái)尋找差異性表達(dá)代謝物。差異性代謝物的定性方法為:根據(jù)保留時(shí)間及m/z在NIST庫(kù)和Wiley Registry代謝組數(shù)據(jù)庫(kù)中查找與其匹配的物質(zhì)。兩地區(qū)樣品差異代謝物的定性結(jié)果見(jiàn)表3。由表3可知,兩產(chǎn)地水稻樣品共篩選出16 種差異代謝產(chǎn)物,由log2(fc-XS/JSJ)值可知,篩選出的16 種差異代謝物中除α-D-甲基呋喃果糖苷外,其余15 種差異代謝物含量均是響水地區(qū)稻米樣品高于建三江地區(qū)的稻米樣品,且倍數(shù)在1.54~3.05之間,含量差異顯著,雖然α-D-甲基呋喃果糖苷的含量是建三江地區(qū)高于響水地區(qū),但是可以看出高出的倍數(shù)為0.02,數(shù)值很小。整體看來(lái)響水地區(qū)稻米中差異代謝物的含量要高于建三江地區(qū),說(shuō)明稻米的產(chǎn)地不同,相同代謝產(chǎn)物的含量也有明顯差異,稻米代謝產(chǎn)物受產(chǎn)地的影響顯著。

    表3 兩地區(qū)樣品差異代謝物的定性結(jié)果Table 3 Qualitative analysis of differential metabolites in rice samples collected from two regions

    2.3 差異性代謝產(chǎn)物分析

    2.3.1 差異代謝物熱圖分析

    圖5 差異代謝物熱圖Fig. 5 Heatmap of the differential metabolites

    通過(guò)R(v3.3.2)中pheatmap程序包對(duì)數(shù)據(jù)集進(jìn)行數(shù)據(jù)縮放,同時(shí)對(duì)樣本和代謝物進(jìn)行雙向聚類。由圖5可知,圖中紅色區(qū)域表示差異代謝物的高表達(dá)區(qū),綠色部分表示差異代謝產(chǎn)物的低表達(dá)區(qū)。通過(guò)篩選出的差異代謝物將兩產(chǎn)地的樣品進(jìn)行聚類分析,響水地區(qū)的稻米樣品主要集中在紅色的高表達(dá)區(qū)域,而建三江地區(qū)的稻米樣品則主要集中在綠色低表達(dá)區(qū),說(shuō)明響水地區(qū)的稻米樣品中差異代謝物的含量普遍高于建三江地區(qū),響水稻米樣品中的差異代謝物的代謝機(jī)能更為活躍,同時(shí)也可以看出同產(chǎn)地樣品的代謝物具有相近的表達(dá)程度,可見(jiàn)產(chǎn)地對(duì)稻米的代謝具有很大的影響。影響水稻生長(zhǎng)的因素非常復(fù)雜,其中溫度與光照是最為重要的因素[27]。響水地區(qū)所處的地理位置位于火山爆發(fā)出的巖漿凝結(jié)的熔巖臺(tái)地,地溫、水溫比一般的稻田地高2~3 ℃左右,在晝夜溫差大的北方,形成了有利的稻米生長(zhǎng)環(huán)境,使其具備更好的溫度和濕度條件,較高的溫度條件能增強(qiáng)水稻的光合作用和根葉的呼吸作用,從而提高了稻米的代謝速率[28],這可能是導(dǎo)致響水地區(qū)的稻米代謝產(chǎn)物高表達(dá)的原因。

    2.3.2 差異代謝物通路分析

    京都基因與基因組百科全書(shū)(Kyoto Encyclopedia of Genes and Gnomes,KEGG)注釋分析僅找到所有差異代謝物參與的通路,MetPA數(shù)據(jù)庫(kù)通過(guò)代謝通路濃縮和拓?fù)浞治?,能識(shí)別出可能的受生物擾動(dòng)的代謝通路,可以對(duì)通路進(jìn)行進(jìn)一步的篩選,找到與代謝物差異相關(guān)性最高的關(guān)鍵通路,結(jié)果如表4所示。

    表4 差異代謝物參與的代謝途徑明細(xì)Table 4 List of metabolic pathways with the involvement of differential metabolites

    富集分析和拓?fù)浞治龉矘?biāo)注到4 個(gè)較為重要的代謝途徑,包括脂肪酸生物合成、精氨酸和脯氨酸代謝、嘧啶代謝、氨酰-tRNA生物合成。脂肪酸可為稻米提供能量,且對(duì)植物本身的抗逆性起重要的作用,使其可以適應(yīng)寒冷的生長(zhǎng)環(huán)境[29]。精氨酸和脯氨酸是稻米氨基酸重要組成,脯氨酸是稻米代謝產(chǎn)物中重要的滲透調(diào)節(jié)物質(zhì),它起到保持原生質(zhì)和環(huán)境滲透平衡、阻止水分喪失的作用。此外,脯氨酸還直接影響蛋白質(zhì)的穩(wěn)定性,還可作為含氮的儲(chǔ)藏物質(zhì)和恢復(fù)生長(zhǎng)的能源,還有研究表明,低溫下水稻幼苗脯氨酸含量會(huì)隨著脅迫時(shí)間的延長(zhǎng)而增加[30-31]。不同產(chǎn)地環(huán)境特點(diǎn)不同,盡管建三江地區(qū)與響水地區(qū)位于同省份內(nèi),且均適合水稻的生長(zhǎng),但建三江地處三江平原腹地,屬寒溫帶濕潤(rùn)季風(fēng)氣候,年平均氣溫1~2 ℃,平均降水量為550~600 mm,而響水大米產(chǎn)區(qū)寧安市屬溫帶大陸性季風(fēng)氣候,年平均氣溫3.5 ℃,平均降水量為400~600 mm,氣溫、降水量均有較大差異,所以其產(chǎn)地生長(zhǎng)的稻米為適應(yīng)環(huán)境會(huì)通過(guò)調(diào)節(jié)代謝水平來(lái)適應(yīng)環(huán)境,以便更好地生長(zhǎng),從而其各種代謝物的含量會(huì)出現(xiàn)差異,同時(shí)也會(huì)賦予稻米不同的品質(zhì)。氨酰-tRNA生物合成是形成蛋白質(zhì)的重要通路。在蛋白質(zhì)的生物合成中,氨酰-tRNA作為氨基酸轉(zhuǎn)移到肽鏈的羧末端過(guò)程的中間物而發(fā)揮重要的作用[32-33]。氨酰-tRNA合成酶在識(shí)別它所催化的氨基酸時(shí)需要ATP,所以光照是氨酰-tRNA生物合成途徑的關(guān)鍵因素[34],響水大米的光照周期明期與暗期比為7∶17,建三江大米的光照周期明期與暗期比為6∶18,兩地區(qū)均有充足的光照可滿足氨酰-tRNA生物合成途徑,但光照情況還是有所差異。產(chǎn)地的差異對(duì)水稻機(jī)體的多種終端代謝途徑造成了一定影響,但產(chǎn)地影響稻米代謝途徑及其機(jī)理尚未清楚,有待進(jìn)一步研究。

    3 結(jié) 論

    本研究基于GC-MS的代謝組學(xué)技術(shù)對(duì)黑龍江省響水地區(qū)和建三江地區(qū)的60 個(gè)的稻米樣品進(jìn)行代謝組學(xué)分析,考察不同產(chǎn)地對(duì)稻米代謝產(chǎn)物的影響。實(shí)驗(yàn)共鑒定出39 個(gè)代謝產(chǎn)物,其中響水地區(qū)包括26 個(gè)代謝產(chǎn)物,建三江地區(qū)包括29 個(gè)代謝產(chǎn)物。對(duì)樣品進(jìn)行PCA和OPLSDA等多元分析,兩產(chǎn)地共篩選出16 種差異代謝物,響水地區(qū)稻米差異代謝物含量要高于建三江地區(qū),且倍數(shù)在1.54~3.05 倍之間,通過(guò)差異代謝物可以對(duì)兩產(chǎn)地的稻米樣品進(jìn)行較好的區(qū)分,說(shuō)明不同產(chǎn)地生長(zhǎng)的稻米其代謝產(chǎn)物的種類和含量均不相同,產(chǎn)地對(duì)稻米的代謝具有重要的影響,而稻米的代謝產(chǎn)物可能攜帶著其產(chǎn)地的信息。綜上所述,基于GC-MS技術(shù)的代謝組學(xué)分析方法用于分析產(chǎn)地對(duì)水稻代謝產(chǎn)物的影響具有一定的可行性,可為因環(huán)境不同而導(dǎo)致的稻米品質(zhì)差別提供可視化的理論基礎(chǔ)。實(shí)驗(yàn)證明利用基于GC-MS聯(lián)用技術(shù),結(jié)合多元統(tǒng)計(jì)分析的方法,可用于分析環(huán)境對(duì)水稻代謝產(chǎn)物的影響,可為稻米的溯源提供理論基礎(chǔ)。

    猜你喜歡
    差異
    “再見(jiàn)”和bye-bye等表達(dá)的意義差異
    JT/T 782的2020版與2010版的差異分析
    相似與差異
    關(guān)于中西方繪畫(huà)差異及對(duì)未來(lái)發(fā)展的思考
    收藏界(2019年3期)2019-10-10 03:16:40
    找句子差異
    DL/T 868—2014與NB/T 47014—2011主要差異比較與分析
    生物為什么會(huì)有差異?
    法觀念差異下的境外NGO立法效應(yīng)
    構(gòu)式“A+NP1+NP2”與“A+NP1+(都)是+NP2”的關(guān)聯(lián)和差異
    論言語(yǔ)行為的得體性與禮貌的差異
    亚洲精品aⅴ在线观看| 亚洲成色77777| 精品一区二区三区视频在线| 午夜福利视频精品| 大码成人一级视频| 中文欧美无线码| 亚洲成人av在线免费| 久久婷婷青草| 韩国高清视频一区二区三区| 久久精品夜色国产| 国产精品一区www在线观看| 最后的刺客免费高清国语| 王馨瑶露胸无遮挡在线观看| 一区二区三区四区激情视频| 91在线精品国自产拍蜜月| 亚洲精品色激情综合| 桃花免费在线播放| 国产麻豆69| 欧美日韩视频精品一区| 亚洲精品久久久久久婷婷小说| 久久国产精品大桥未久av| 成人午夜精彩视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 又黄又爽又刺激的免费视频.| 香蕉精品网在线| 亚洲人与动物交配视频| 久久精品夜色国产| 免费看av在线观看网站| 午夜福利视频在线观看免费| 日本91视频免费播放| av在线app专区| 精品亚洲成a人片在线观看| 男女午夜视频在线观看 | 国产熟女欧美一区二区| 久久久精品免费免费高清| 国产精品蜜桃在线观看| 亚洲综合精品二区| 久久青草综合色| 精品一区二区免费观看| 午夜日本视频在线| 全区人妻精品视频| 十八禁网站网址无遮挡| 日韩视频在线欧美| 国产在视频线精品| 男女高潮啪啪啪动态图| 国产国语露脸激情在线看| 久久久久久久久久人人人人人人| 熟女av电影| 内地一区二区视频在线| 久久精品国产综合久久久 | 乱人伦中国视频| 亚洲五月色婷婷综合| 国产成人a∨麻豆精品| 国国产精品蜜臀av免费| 精品久久久久久电影网| 婷婷色麻豆天堂久久| 大陆偷拍与自拍| 午夜福利,免费看| av在线老鸭窝| 欧美成人精品欧美一级黄| 人人澡人人妻人| 大码成人一级视频| 日韩中字成人| av.在线天堂| 欧美激情国产日韩精品一区| 在线观看一区二区三区激情| 黄色 视频免费看| 午夜激情av网站| 亚洲婷婷狠狠爱综合网| 午夜福利视频精品| 免费av中文字幕在线| 亚洲国产av影院在线观看| 国产精品国产三级专区第一集| 日本-黄色视频高清免费观看| 天天操日日干夜夜撸| 1024视频免费在线观看| 亚洲欧美成人综合另类久久久| 精品人妻偷拍中文字幕| 日韩一本色道免费dvd| 久久久亚洲精品成人影院| 亚洲精品自拍成人| 亚洲国产精品一区二区三区在线| 亚洲成人av在线免费| 久久久欧美国产精品| 国产深夜福利视频在线观看| 久久人人97超碰香蕉20202| 青春草国产在线视频| 大香蕉97超碰在线| 美女福利国产在线| 天美传媒精品一区二区| av一本久久久久| 在线观看国产h片| 久久久久视频综合| 各种免费的搞黄视频| 午夜激情久久久久久久| 久久久久精品久久久久真实原创| 免费人妻精品一区二区三区视频| 亚洲一区二区三区欧美精品| 人人妻人人澡人人爽人人夜夜| 成人漫画全彩无遮挡| 久久ye,这里只有精品| 伦理电影免费视频| 国国产精品蜜臀av免费| 国国产精品蜜臀av免费| 亚洲国产欧美在线一区| 综合色丁香网| 亚洲欧洲精品一区二区精品久久久 | 久久久精品免费免费高清| 国产成人精品一,二区| 亚洲伊人色综图| 免费黄网站久久成人精品| 18禁观看日本| 欧美日韩亚洲高清精品| 麻豆精品久久久久久蜜桃| 97精品久久久久久久久久精品| 麻豆乱淫一区二区| 大香蕉97超碰在线| 久久久久精品性色| 一本久久精品| 亚洲高清免费不卡视频| 久久久久久久精品精品| 成人影院久久| 国国产精品蜜臀av免费| 日本av手机在线免费观看| 美女内射精品一级片tv| 欧美亚洲日本最大视频资源| 日韩制服丝袜自拍偷拍| 国产av一区二区精品久久| 最近最新中文字幕大全免费视频 | av播播在线观看一区| 精品国产一区二区久久| √禁漫天堂资源中文www| 一区二区三区四区激情视频| 成年人午夜在线观看视频| 蜜桃国产av成人99| 国产又爽黄色视频| 亚洲,欧美,日韩| 亚洲第一区二区三区不卡| 香蕉国产在线看| 一本色道久久久久久精品综合| 制服人妻中文乱码| av电影中文网址| 啦啦啦啦在线视频资源| 丝袜喷水一区| 久久av网站| 免费女性裸体啪啪无遮挡网站| 日本与韩国留学比较| 亚洲成国产人片在线观看| 国产精品女同一区二区软件| av国产久精品久网站免费入址| 久久精品夜色国产| 亚洲少妇的诱惑av| 色5月婷婷丁香| 最近中文字幕2019免费版| 亚洲美女黄色视频免费看| 免费人妻精品一区二区三区视频| 日本vs欧美在线观看视频| 国产亚洲最大av| 人妻人人澡人人爽人人| 久久久精品区二区三区| 制服丝袜香蕉在线| 国产日韩欧美在线精品| 精品一区二区三区四区五区乱码 | 一区二区日韩欧美中文字幕 | 国产一级毛片在线| 国产 精品1| 亚洲图色成人| 久久精品熟女亚洲av麻豆精品| 久久久精品区二区三区| 日韩av免费高清视频| 精品午夜福利在线看| 2018国产大陆天天弄谢| 午夜久久久在线观看| 免费av中文字幕在线| 九色成人免费人妻av| xxx大片免费视频| 国产亚洲午夜精品一区二区久久| 在线观看免费高清a一片| 一区二区三区乱码不卡18| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜爱| 久久99热这里只频精品6学生| 最近中文字幕高清免费大全6| 日韩成人av中文字幕在线观看| 在线观看美女被高潮喷水网站| 观看美女的网站| 久久精品久久精品一区二区三区| 亚洲,一卡二卡三卡| 大片免费播放器 马上看| 丰满少妇做爰视频| 精品久久久精品久久久| 亚洲av免费高清在线观看| 美女国产高潮福利片在线看| 最近的中文字幕免费完整| 天堂俺去俺来也www色官网| 丰满饥渴人妻一区二区三| 国产一级毛片在线| 亚洲精品一二三| 久久精品aⅴ一区二区三区四区 | av不卡在线播放| 卡戴珊不雅视频在线播放| 婷婷色综合www| 黄网站色视频无遮挡免费观看| 国产亚洲精品第一综合不卡 | 免费av中文字幕在线| 2021少妇久久久久久久久久久| 最近2019中文字幕mv第一页| 啦啦啦啦在线视频资源| 久久国产精品男人的天堂亚洲 | 两个人看的免费小视频| av女优亚洲男人天堂| 欧美精品高潮呻吟av久久| 人妻 亚洲 视频| 久久久精品免费免费高清| 丝袜喷水一区| 免费在线观看完整版高清| 中文字幕人妻丝袜制服| 少妇人妻 视频| 91久久精品国产一区二区三区| 在线 av 中文字幕| 婷婷色av中文字幕| 韩国高清视频一区二区三区| 亚洲精品一区蜜桃| 丝袜喷水一区| 18禁国产床啪视频网站| 亚洲人成77777在线视频| 欧美精品一区二区免费开放| 男女国产视频网站| 亚洲欧美中文字幕日韩二区| 丝瓜视频免费看黄片| 巨乳人妻的诱惑在线观看| 国产探花极品一区二区| 亚洲国产精品一区二区三区在线| freevideosex欧美| 色婷婷av一区二区三区视频| 男女高潮啪啪啪动态图| 高清欧美精品videossex| 亚洲图色成人| 90打野战视频偷拍视频| 国产一区二区在线观看日韩| 日韩大片免费观看网站| 久久这里只有精品19| 亚洲成国产人片在线观看| 亚洲 欧美一区二区三区| 亚洲国产精品成人久久小说| 久久久a久久爽久久v久久| 成年女人在线观看亚洲视频| 18禁国产床啪视频网站| 99国产综合亚洲精品| 久久久久精品久久久久真实原创| 大片电影免费在线观看免费| av不卡在线播放| 成人二区视频| 精品视频人人做人人爽| 国产精品.久久久| a级毛片在线看网站| 精品一区二区三区四区五区乱码 | 亚洲第一区二区三区不卡| 男女免费视频国产| 1024视频免费在线观看| 人人妻人人澡人人看| 777米奇影视久久| 日韩成人av中文字幕在线观看| 2021少妇久久久久久久久久久| 欧美亚洲日本最大视频资源| 自线自在国产av| 久久久精品免费免费高清| 亚洲经典国产精华液单| 边亲边吃奶的免费视频| 五月伊人婷婷丁香| 国产精品欧美亚洲77777| av在线播放精品| 大片免费播放器 马上看| 91国产中文字幕| 欧美精品一区二区免费开放| 男女边吃奶边做爰视频| 18禁在线无遮挡免费观看视频| 国产精品蜜桃在线观看| 亚洲性久久影院| 成人手机av| 啦啦啦中文免费视频观看日本| 春色校园在线视频观看| 日韩一区二区三区影片| 国产色婷婷99| 边亲边吃奶的免费视频| 亚洲一级一片aⅴ在线观看| 性色avwww在线观看| 丰满饥渴人妻一区二区三| 久久久久久人人人人人| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品999| 亚洲第一区二区三区不卡| 少妇人妻 视频| 午夜免费鲁丝| 免费看av在线观看网站| 午夜视频国产福利| 老司机影院毛片| 91午夜精品亚洲一区二区三区| 亚洲成国产人片在线观看| 嫩草影院入口| av不卡在线播放| 国产爽快片一区二区三区| 成年动漫av网址| 精品酒店卫生间| 中文欧美无线码| 国产成人91sexporn| 在线观看国产h片| 亚洲精品日本国产第一区| 男人添女人高潮全过程视频| 国产在线视频一区二区| 中文天堂在线官网| 国产成人精品婷婷| 欧美人与性动交α欧美软件 | 黄片无遮挡物在线观看| 美女xxoo啪啪120秒动态图| 性色avwww在线观看| 亚洲国产精品成人久久小说| 亚洲精品中文字幕在线视频| 男人舔女人的私密视频| 黄色配什么色好看| 亚洲国产av影院在线观看| 久久国内精品自在自线图片| 久久99热6这里只有精品| 亚洲少妇的诱惑av| 看免费成人av毛片| av播播在线观看一区| 91久久精品国产一区二区三区| 在线 av 中文字幕| 精品亚洲乱码少妇综合久久| 黄色一级大片看看| 国产 一区精品| 亚洲欧洲精品一区二区精品久久久 | 久久精品aⅴ一区二区三区四区 | 亚洲国产色片| tube8黄色片| 色网站视频免费| 哪个播放器可以免费观看大片| 亚洲综合色网址| 国产精品久久久久成人av| 亚洲精品国产av蜜桃| 色5月婷婷丁香| 一级,二级,三级黄色视频| 色吧在线观看| 亚洲,一卡二卡三卡| 999精品在线视频| 成人午夜精彩视频在线观看| 成年美女黄网站色视频大全免费| 国产一区二区在线观看日韩| www.熟女人妻精品国产 | 成人免费观看视频高清| 国产一区二区三区av在线| 18禁裸乳无遮挡动漫免费视频| 国产免费视频播放在线视频| 女的被弄到高潮叫床怎么办| 成人免费观看视频高清| 桃花免费在线播放| 国产老妇伦熟女老妇高清| 国产又爽黄色视频| 久久99一区二区三区| 国产亚洲午夜精品一区二区久久| av女优亚洲男人天堂| 欧美精品高潮呻吟av久久| 熟女电影av网| 国产麻豆69| 大香蕉97超碰在线| 日韩欧美一区视频在线观看| 久久99一区二区三区| 欧美精品av麻豆av| 国产69精品久久久久777片| 成年女人在线观看亚洲视频| 久久人妻熟女aⅴ| 有码 亚洲区| 狂野欧美激情性bbbbbb| 人妻一区二区av| 亚洲av成人精品一二三区| 中国国产av一级| 亚洲av免费高清在线观看| 内地一区二区视频在线| 亚洲性久久影院| 国产福利在线免费观看视频| 久久久久视频综合| 日韩 亚洲 欧美在线| 伦精品一区二区三区| 美女国产视频在线观看| 亚洲精品av麻豆狂野| 在线免费观看不下载黄p国产| 久久婷婷青草| 在线观看免费日韩欧美大片| 青青草视频在线视频观看| 国产乱人偷精品视频| 看十八女毛片水多多多| 一本—道久久a久久精品蜜桃钙片| 欧美国产精品va在线观看不卡| 午夜视频国产福利| 69精品国产乱码久久久| 一本久久精品| 午夜精品国产一区二区电影| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 在线观看美女被高潮喷水网站| 下体分泌物呈黄色| 18禁国产床啪视频网站| 日韩制服骚丝袜av| 三级国产精品片| videosex国产| 精品久久国产蜜桃| 亚洲美女搞黄在线观看| 国产精品久久久久久久电影| 亚洲色图 男人天堂 中文字幕 | 91精品伊人久久大香线蕉| 高清毛片免费看| 亚洲欧美一区二区三区黑人 | 国产日韩一区二区三区精品不卡| 9热在线视频观看99| 久久精品国产亚洲av涩爱| 亚洲国产av新网站| av不卡在线播放| 精品视频人人做人人爽| 香蕉国产在线看| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 欧美日韩视频精品一区| 最近手机中文字幕大全| 一级爰片在线观看| 黑丝袜美女国产一区| 久久久久久久久久人人人人人人| 最近最新中文字幕免费大全7| 国产成人av激情在线播放| 99热全是精品| www.av在线官网国产| 观看av在线不卡| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 大香蕉久久成人网| 91精品伊人久久大香线蕉| 在线观看一区二区三区激情| 日韩av在线免费看完整版不卡| xxxhd国产人妻xxx| 国产亚洲精品第一综合不卡 | 精品少妇久久久久久888优播| 亚洲精品国产av蜜桃| 少妇高潮的动态图| 女人精品久久久久毛片| 99久国产av精品国产电影| 看免费成人av毛片| 大片电影免费在线观看免费| 青春草国产在线视频| 久久久欧美国产精品| 精品福利永久在线观看| 天美传媒精品一区二区| 极品少妇高潮喷水抽搐| 精品人妻偷拍中文字幕| 观看美女的网站| 美女中出高潮动态图| 精品久久久久久电影网| 久久99热6这里只有精品| 国产精品久久久久久久电影| av又黄又爽大尺度在线免费看| 国产免费一区二区三区四区乱码| 男女无遮挡免费网站观看| 黄片无遮挡物在线观看| 欧美国产精品va在线观看不卡| 在线精品无人区一区二区三| 午夜久久久在线观看| 在线观看人妻少妇| 全区人妻精品视频| 欧美日韩视频精品一区| 最近中文字幕高清免费大全6| 免费黄网站久久成人精品| 全区人妻精品视频| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| 国产黄频视频在线观看| 国产又爽黄色视频| 成人手机av| av卡一久久| 久久久久精品人妻al黑| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 男人操女人黄网站| 国产亚洲一区二区精品| 人体艺术视频欧美日本| 啦啦啦中文免费视频观看日本| 考比视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 三级国产精品片| 日本免费在线观看一区| 91在线精品国自产拍蜜月| 男女无遮挡免费网站观看| 久久久久久伊人网av| 久久久久久人人人人人| 最新的欧美精品一区二区| 在现免费观看毛片| 久久国产精品大桥未久av| av国产精品久久久久影院| 美女中出高潮动态图| 欧美日韩精品成人综合77777| 下体分泌物呈黄色| 国产高清不卡午夜福利| 久久午夜综合久久蜜桃| 免费少妇av软件| 最新中文字幕久久久久| 亚洲精品第二区| 久久人妻熟女aⅴ| 一二三四中文在线观看免费高清| 国产色婷婷99| 精品国产一区二区久久| 欧美精品亚洲一区二区| 国产精品久久久久久久电影| 久久韩国三级中文字幕| 免费在线观看完整版高清| 自拍欧美九色日韩亚洲蝌蚪91| 狠狠精品人妻久久久久久综合| 日韩av在线免费看完整版不卡| 精品久久久精品久久久| 亚洲精品美女久久久久99蜜臀 | 成人手机av| 狂野欧美激情性bbbbbb| 国产熟女欧美一区二区| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 久久久精品94久久精品| 成年女人在线观看亚洲视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品国产精品| 搡女人真爽免费视频火全软件| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 久久精品国产综合久久久 | 欧美少妇被猛烈插入视频| 精品视频人人做人人爽| 在线观看三级黄色| 亚洲欧美色中文字幕在线| www.色视频.com| 老熟女久久久| 一区二区av电影网| 免费黄网站久久成人精品| 亚洲色图 男人天堂 中文字幕 | 免费少妇av软件| 18禁裸乳无遮挡动漫免费视频| 乱码一卡2卡4卡精品| 国产成人91sexporn| 99久久人妻综合| 国产国语露脸激情在线看| 欧美日韩亚洲高清精品| 香蕉国产在线看| 久久人人爽人人片av| 男女边摸边吃奶| 亚洲av免费高清在线观看| 最近中文字幕2019免费版| 高清不卡的av网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲性久久影院| 观看美女的网站| 亚洲国产日韩一区二区| 91精品国产国语对白视频| 成人亚洲精品一区在线观看| 观看美女的网站| 亚洲精华国产精华液的使用体验| 免费av中文字幕在线| 久久精品熟女亚洲av麻豆精品| 精品国产乱码久久久久久小说| 午夜激情av网站| 日本黄大片高清| 久久久久网色| 人成视频在线观看免费观看| 女性被躁到高潮视频| 天天操日日干夜夜撸| 最近2019中文字幕mv第一页| 日韩一本色道免费dvd| 欧美xxxx性猛交bbbb| 亚洲精品乱码久久久久久按摩| 国产精品一区www在线观看| 男的添女的下面高潮视频| 人妻少妇偷人精品九色| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 亚洲精品aⅴ在线观看| 国产黄色视频一区二区在线观看| 亚洲av国产av综合av卡| 99热这里只有是精品在线观看| 丝袜人妻中文字幕| 丰满饥渴人妻一区二区三| 国产极品天堂在线| 国产亚洲一区二区精品| 波多野结衣一区麻豆| 熟妇人妻不卡中文字幕| 少妇人妻久久综合中文| 丁香六月天网| 国产国拍精品亚洲av在线观看| 考比视频在线观看| 国产日韩一区二区三区精品不卡| 十八禁网站网址无遮挡| 久久青草综合色| 五月开心婷婷网| 美女福利国产在线| 国产精品一区二区在线观看99| 91精品伊人久久大香线蕉| 九九在线视频观看精品| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 国产日韩欧美视频二区| 国产 一区精品| 自线自在国产av| 免费不卡的大黄色大毛片视频在线观看| 精品一区二区三卡| 久久久国产一区二区| 成人综合一区亚洲| 多毛熟女@视频| 欧美激情 高清一区二区三区| 丰满迷人的少妇在线观看| 亚洲精品美女久久久久99蜜臀 | 搡老乐熟女国产| 人妻一区二区av| 一区二区三区四区激情视频|