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

    基于氮穩(wěn)定同位素比例質(zhì)譜和液相色譜-高分辨質(zhì)譜的有機番茄鑒別

    2021-07-29 02:10:14王世成李國琛高凝軒王顏紅
    食品科學(xué) 2021年14期
    關(guān)鍵詞:同位素比值質(zhì)譜

    王世成,李國琛,王 瑩,李 波,高凝軒,孫 辭,宋 陽,王顏紅

    (中國科學(xué)院沈陽應(yīng)用生態(tài)研究所,遼寧省農(nóng)產(chǎn)品環(huán)境與質(zhì)量安全控制工程技術(shù)研究中心,遼寧 沈陽 110016)

    近年來,農(nóng)業(yè)生產(chǎn)從數(shù)量滿足型向質(zhì)量需求型轉(zhuǎn)變,環(huán)境友好型的循環(huán)農(nóng)業(yè)也日益凸顯優(yōu)勢。隨著生活水平的提高,消費者需求也逐漸從傳統(tǒng)食品升級為綠色食品、有機食品。有機食品在生產(chǎn)過程中拒絕使用化學(xué)性肥料和農(nóng)藥,因而具有較高的安全性和經(jīng)濟價值。受利益驅(qū)使,一些不法商家隨意標(biāo)注“有機”字樣,假冒“有機食品”坑害消費者。歐洲議會“公共衛(wèi)生和食品安全委員會”將有機食品列為10 個風(fēng)險較高的摻假食品的第3位[1]。因此,鑒別有機產(chǎn)品的真實性對于保障消費者安全和市場誠信,保護消費者和有機品牌的經(jīng)濟利益至關(guān)重要。

    在對有機與非有機種植模式的農(nóng)產(chǎn)品鑒定研究中,有近紅外光譜[2]、高效液相色譜[3-4]和氮同位素比例質(zhì)譜[1,5]等技術(shù)得到應(yīng)用,但僅通過一種檢測手段或只對單一目標(biāo)物檢測的技術(shù)很難滿足鑒別需求,將光譜、色譜以及質(zhì)譜結(jié)合的組學(xué)及化學(xué)計量學(xué)分析技術(shù)已逐漸成為研究熱點[6]。

    有機食品是指根據(jù)有機農(nóng)業(yè)生產(chǎn)的規(guī)范生產(chǎn)加工并經(jīng)過認證的農(nóng)產(chǎn)品和加工品。有機食品禁用化肥,允許施用有機肥,土壤有機肥的施用是種植業(yè)有機產(chǎn)品生產(chǎn)的主要控制要素[7]。研究表明[8-10],化學(xué)肥料中δ15N的豐度變化范圍小,平均值為0.2‰,而有機肥料中δ15N的豐度變化范圍大,平均值為8.5‰,有機肥中15N同位素豐度高于化肥中15N豐度。有機模式種植的農(nóng)產(chǎn)品主要施用有機肥,不施用化肥,因此致使農(nóng)產(chǎn)品中氮同位素比值(δ15N)存在差異,通常表現(xiàn)為有機產(chǎn)品中δ15N高于施用化肥的普通產(chǎn)品。農(nóng)產(chǎn)品中δ15N豐度依次為施用有機肥>不施肥(土壤內(nèi)源氮)>施用化學(xué)合成氮肥[11]。對有機番茄的化學(xué)成分研究表明[12-13],有機番茄中山柰酚、槲皮素等類黃酮化合物高于普通番茄,且有機種植的番茄中的類黃酮含量隨種植時間延長而增加,而傳統(tǒng)種植的類黃酮含量并無顯著差異。因此,利用氮穩(wěn)定同位素結(jié)合化學(xué)成分的變化,可以甄別有機食品與普通食品。本實驗對有機和普通種植番茄的穩(wěn)定同位素比值特征與化學(xué)組成成分進行對比研究。

    1 材料與方法

    1.1 材料與試劑

    番茄:番茄樣品共計94 個,包括寶石、大紅、3944、安哥拉等大果番茄品種,其中有機番茄樣品50 個,普通番茄樣品44 個。采集于2017年6月至10月期間。地點分別位于遼寧沈陽新民、沈陽法庫、江蘇南京、山東濰坊、四川西昌的有機番茄及附近的普通番茄生產(chǎn)基地,其中有機番茄采自當(dāng)?shù)匕从袡C種植方式生產(chǎn)的有機農(nóng)場或企業(yè)。每個產(chǎn)地樣品的采摘方式為隨機選擇10 個有機種植模式的番茄植株和普通種植模式的番茄植株,在植株的不同部位采摘成熟度一致的果實。為增大采集樣品數(shù)量,每株植株隨機摘取4~6 枚果實,所摘取的番茄大小均勻,具有代表性。每個植株的番茄為1 組。樣品的來源及相關(guān)信息見表1。

    表1 有機與普通番茄樣品的產(chǎn)地來源Table 1 Geographical origins of organic and common tomato samples tested in this study

    甲醇、乙腈(均為色譜級) 德國默克公司;乙醇、甲酸(均為色譜級) 天津市科密歐化學(xué)試劑有限公司;超純水經(jīng)Millipore超純水系統(tǒng)純化。

    1.2 儀器與設(shè)備

    HPLC-Q Exactive高效液相色譜-四極桿-靜電場軌道阱高分辨質(zhì)譜(high performance liquid chromatographyhigh resolution mass spectrometry,HPLC-HRMS)聯(lián)用儀美國Thermo Scientific公司;Deltaplus XP穩(wěn)定同位素比例質(zhì)譜(isotope ratio mass spectrometry,IRMS)儀、Flash EA1112元素分析儀、ConFlo III元素分析儀-IRMS連續(xù)流接口 美國Thermo-Finnigan公司;SQP電子天平賽多利科學(xué)儀器(北京)有限公司;KQ-250B型超聲波清洗器 昆山市超聲儀器有限公司;TTL-DC氮吹儀II型 北京同泰聯(lián)科技發(fā)展有限公司;CT15RE冷凍離心機 日本日立公司。

    1.3 方法

    1.3.1 樣品的制備

    采集獲得的番茄樣品洗凈切片,混勻后,用勻漿機打碎,冷凍干燥成粉末,過20 目篩放于-80 ℃?zhèn)溆谩?/p>

    1.3.2 HPLC-HRMS分析

    凍干后的樣品,稱量約150 mg,加入3 mL水-乙醇(1∶4,V/V),超聲30 min后,4 000 r/min離心15 min,取2 mL上清液于氮氣流下?lián)]發(fā)(45 ℃)至近干。殘渣加水-甲醇(1∶4,V/V)2 mL進行復(fù)溶,超聲30 min,4 ℃、14 000 r/min離心15 min,過0.22 μm濾膜,分取1.5 mL上清液于液相進樣瓶中,進行分析。

    1.3.3 氮同位素比值測定

    凍干后的樣品,分別稱量約2 mg,并用錫杯包裹壓實,利用元素分析儀-同位素比值質(zhì)譜儀測定穩(wěn)定氮同位素比值。每份樣品重復(fù)3 次,并根據(jù)下式計算δ15N值[14-17]:

    式中:(15N/14N)樣品為樣品穩(wěn)定氮同位素比值;(15N/14N)標(biāo)準為標(biāo)準物質(zhì)穩(wěn)定氮同位素比值。

    1.3.4 液相色譜條件

    采用Waters Atlantis T3色譜柱(2.1 mm×150 mm,3 μm);流動相為0.1%甲酸-乙腈溶液(A)和0.1%甲酸溶液(B);流速0.3 mL/min;柱溫45 ℃;梯度洗脫條件:0~0.5 min,5% A,95% B;0.5~5 min,5%~30% A,95%~70% B;5~9.5 min,30%~90% A,70%~10% B;9.5~9.75 min,90%~5% A,10%~95% B;9.75~12 min,5% A,95% B。進樣體積5 μL。

    1.3.5 高分辨質(zhì)譜條件

    電噴霧離子源;噴霧電壓3.2 kV;蒸汽溫度350 ℃;離子傳輸管溫度320 ℃;鞘氣流速35 arb;輔助氣流速15 arb;Full MS-ddMS2掃描模式,F(xiàn)ull MS質(zhì)量掃描范圍m/z100~900,一級分辨率70 000,二級分辨率17 500。

    1.3.6 同位素比值測試條件

    元素分析系統(tǒng)氧化柱填料Cr2O3/Co3O4,溫度為960 ℃,還原柱填料為高純Cu/CuO,溫度為640 ℃,載氣He流量為90~100 mL/min,氧噴條件為175 mL/min。

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

    HPLC-HRMS數(shù)據(jù)采用Thermo Scientific公司的Xcalibur4.0數(shù)據(jù)軟件進行原始數(shù)據(jù)采集。Compound Discover軟件(version 3.0, Thermo Scientific)對總離子流圖中的色譜峰進行時間和峰對齊校正、濾噪,數(shù)據(jù)處理分析,與數(shù)據(jù)庫進行計算匹配。處理結(jié)果導(dǎo)出樣品名稱、保留時間、分子式、分子質(zhì)量、質(zhì)荷比以及對應(yīng)的離子峰面積強度組成的數(shù)據(jù)集。主要參數(shù)設(shè)置如下:保留時間范圍為1~12 min;質(zhì)荷比范圍為m/z100~900;偏差為5×10-6;信噪比為3;保留時間漂移值為0.2 min。得到的峰面積利用SIMCA-P 11.5軟件(Umetrics, Umea, Sweden)進行化學(xué)計量學(xué)分析。采用主成分分析(principal component analysis,PCA)非監(jiān)督分析方法進行探索性數(shù)據(jù)分析,實現(xiàn)樣本間分析聯(lián)系的可視化。使用SIMCA-P剔除異常值樣品,使用有監(jiān)督的偏最小二乘-判別分析(partial least squares-discriminant analysis,PLS-DA)方法對樣本及已確定的類別進行建模[18-22]。采用多變量和單變量分析選擇法選定差異顯著的變量,篩選出統(tǒng)計學(xué)意義顯著的化合物[23-26]。化合物的初步鑒定是在準確的物質(zhì)組成基礎(chǔ)上進行,從精確質(zhì)量數(shù)、同位素組成及分布與預(yù)測分子式的吻合度、一級質(zhì)譜/二級質(zhì)譜碎片和數(shù)據(jù)庫匹配等方面進行推測。

    IRMS數(shù)據(jù)利用SPSS 19.0軟件進行統(tǒng)計分析,并利用Excel和Origin 8.5軟件進行繪圖統(tǒng)計處理。

    2 結(jié)果與分析

    2.1 穩(wěn)定同位素比值分析

    2.1.1 有機與普通番茄穩(wěn)定氮同位素比值的差異分析

    利用元素分析儀-IRMS測得不同地區(qū)50 種有機種植番茄樣品以及44 種普通種植番茄樣品,由圖1A可以看出,有機番茄δ15N值的中位值在8.68‰左右,而普通方式種植的番茄δ15N值的中位值則為3.61‰,有機番茄氮同位素比值明顯高于普通番茄。有機番茄的肥料主要為有機肥,普通番茄主要施用化肥。有機肥以畜禽糞便、動植物殘體等富含有機質(zhì)的副產(chǎn)品資源為主要原料,Peterson等[10]研究發(fā)現(xiàn),植物氮源主要來自于土壤中動植物腐爛物質(zhì)及動物排泄物。由于15N與14N在生物體內(nèi)反應(yīng)所需能量不同,反應(yīng)速率存在差異,造成15N在生物體內(nèi)產(chǎn)生富集。因而有機肥15N豐度高于化肥,有機肥中高含量的15N使有機作物果實中15N比例增高[8,11]。

    圖1 所有樣品(A)及不同產(chǎn)地來源(B)的有機與普通番茄穩(wěn)定氮同位素比值Fig.1 Stable nitrogen isotopic ratios of organic and common tomato samples (A) and samples from different geographical origins (B)

    所有有機和普通番茄穩(wěn)定同位素比值在一定范圍內(nèi)互有重疊(圖1A),但按4 個地區(qū)分組比較發(fā)現(xiàn),同一地區(qū)所產(chǎn)有機番茄的δ15N均高于普通番茄,并且存在顯著差異(P<0.01),且無重疊(圖1B)。因此,可以通過比較同一地區(qū)番茄的δ15N初步判斷是否為有機番茄。進一步比較有機番茄和普通番茄種植土壤的同位素比值發(fā)現(xiàn),除四川西昌外,其他4 個地區(qū)的有機番茄土壤的δ15N值也均高于普通土壤,多年有機種植模式下土壤中氮同位素比值有增加的趨勢。但統(tǒng)計分析發(fā)現(xiàn)除四川土壤外,其他地區(qū)的有機和普通土壤二者無顯著差異(圖2)。

    圖2 土壤及對應(yīng)番茄穩(wěn)定氮同位素比值Fig.2 Sable nitrogen isotope ratios of soil and corresponding tomatoes

    2.1.2 番茄與對應(yīng)的土壤氮同位素變化特征

    番茄與對應(yīng)土壤的δ15N值對比發(fā)現(xiàn),除遼寧一個產(chǎn)地(法庫)外,其他4 個產(chǎn)地有機番茄的氮同位素比值均高于對應(yīng)的土壤同位素比值(圖2);除山東外,普通番茄δ15N值均低于對應(yīng)的土壤(圖2),在土壤氮源中δ15N豐度差異不顯著的情況下,施用外源肥料后,有機番茄主要吸收有機肥中的高δ15N的有效氮,導(dǎo)致有機番茄中δ15N高于土壤的δ15N,普通番茄吸收無機肥的低δ15N有效氮,導(dǎo)致普通番茄的同位素比值低于土壤同位素比值[8]。

    2.2 液相色譜高分辨質(zhì)譜分析番茄化學(xué)組成成分

    2.2.1 液相色譜分析結(jié)果

    分別采集番茄樣品的HPLC-HRMS正離子和負離子模式的總離子流色譜圖。圖3為具有代表性的有機番茄樣品正離子模式指紋圖譜。在正離子模式下共有約25 個明顯可見的色譜峰,而在負離子模式下,樣品出峰較少,在扣除試劑空白后,僅有5 個明顯的色譜峰存在。

    圖3 番茄樣品的HPLC-HRMS總離子流色譜圖(正離子模式)Fig.3 HPLC-HRMS total ion current chromatogram of tomato samples (in the positive ion mode)

    2.2.2 液相色譜-高分辨質(zhì)譜的多元統(tǒng)計判別分析

    94 個樣品的色譜峰面積經(jīng)SIMCA-P軟件剔除異常值樣品,獲得符合要求的樣品91 個,隨機取56 個作為訓(xùn)練集,其中有機樣品28 個,普通樣品28 個,35 個作為驗證集樣品,其中有機樣品20 個,普通樣品15 個。

    應(yīng)用Compound Discover軟件,對色譜峰進行峰對齊、保留時間校正和濾噪處理后,提取其中的所有質(zhì)譜信息。使用軟件的化合物篩查過濾模式,提取化合物峰面積大于105,且60%以上樣品都有響應(yīng)的化合物組成成分,提取其中的875 個化合物的峰面積數(shù)據(jù)。采用SIMCA-P軟件進行PCA和PLS-DA,PCA不能將樣品很好分開,采用PLS-DA提取3 個PC,模型解釋的訓(xùn)練集累計貢獻率(R2Y)達到0.915(通常R2越接近1,模型越良好),模型交叉驗證對訓(xùn)練集預(yù)測的累計貢獻率(Q2)為0.539(通常Q2>0.5表明模型對新數(shù)據(jù)的預(yù)測能力良好)。采用Permutation校驗方法[21]驗證模型是否存在過擬合現(xiàn)象(圖4),通過200 次迭代的排列測試,模型結(jié)果顯示Q2點散射回歸線與垂直軸的截距小于零,說明PLS-DA模型驗證結(jié)果良好,數(shù)據(jù)模型擬合成功,模型可靠。建立的模型能夠?qū)崿F(xiàn)2 組樣本完全分離,模型的總預(yù)測正確率為92.3%(表2)。

    圖4 PLS-DA模型的200 次permutation校驗結(jié)果Fig.4 Results of verification of PLS-DA model with 200 permutation tests

    表2 模型對有機和普通番茄的判定結(jié)果Table 2 Results of the model for organic and common tomatoes

    2.2.3 有機番茄和普通番茄組成成分的單變量差異性分析

    采用Compound Discover軟件,同2.2.1節(jié)處理,對番茄樣品的HPLC-HRMS組間數(shù)據(jù)峰面積均值進行比對,形成單變量統(tǒng)計分析結(jié)果的火山圖(圖5)。圖中框圖標(biāo)注區(qū)域為有機番茄和普通番茄單變量統(tǒng)計差異顯著(P<0.05),并且有機番茄相比于普通番茄(普通/有機),面積比值增加(log2fold change<0.5)或降低(log2fold change>2)的特征化合物,共有18 個,這些化合物成分在有機和普通番茄中存在較大差異,需進一步確證。

    圖5 不同化合物的有機與普通番茄組間的火山圖Fig.5 Volcano plots showing differential compounds between organic and common tomatoes

    2.2.4δ15N與HPLC-HRMS組合的化學(xué)計量學(xué)判別分析

    為了提高模型的預(yù)測能力,將δ15N數(shù)據(jù)與875 個色譜峰數(shù)據(jù)合并形成新的數(shù)據(jù)集,采用SIMCA-P軟件進行PLS-DA,提取3 個PC,模型解釋的訓(xùn)練集累計貢獻率(R2Y)達到0.918,交叉驗證模型預(yù)測的訓(xùn)練集累計貢獻率(Q2)為0.559,模型預(yù)測的正確率為93.4%(表2)。與δ15N變量加入前相比,信息提取率、預(yù)測能力和準確率均有一定提高。前3 個PC的得分圖上,2 組樣本實現(xiàn)了完全分離(圖6a)。在前2 個PC的載荷圖(圖6b)中,距離原點最遠的變量為318.07和δ15N,說明這2 個指標(biāo)對判定有機番茄的貢獻最大,其次為分子質(zhì)量192.03、208.15、274.19、159.07、609.39、305.27、453.29 Da等變量,也有可能成為有機和普通番茄產(chǎn)生差異的主要標(biāo)志化合物。

    圖6 有機和普通番茄的得分圖(a)及載荷圖(b)Fig.6 Score plot (a) and loading plot (b) for organic and common tomatoes

    2.2.5 特征變量因子篩查與確證分析

    將SIMCA-P軟件多變量處理結(jié)果與Compound Discover軟件單變量分析的結(jié)果相結(jié)合,進一步篩選特征化合物,以更全面闡明有機和普通樣品之間的差異化合物成分。單變量處理的組間存在顯著差異(P<0.05)的219 個變量中,滿足多變量處理變量重要性投影(variable importance in projection,VIP)大于1.0的變量153 個,其中VIP值大于1.5的變量53 個,這53 個化合物可能成為有機番茄和普通番茄具有顯著差異的特征化合物。這些化合物與采用Compound Discover處理的火山圖結(jié)果(圖5)中的18 個化合物中有14 個一致,這14 個成分中有9 個化合物在有機番茄中含量低于普通番茄,有5 個高于普通番茄,說明通過簡單的化合物峰面積比較可初步發(fā)現(xiàn)存在含量差異的化合物。

    針對篩選出的除δ15N以外的152 個變量,通過精確分子質(zhì)量、二級質(zhì)譜分析,并檢索mzclound、chemspider、massbank、KEGG等數(shù)據(jù)庫,結(jié)合文獻[27-31],初步確定mzclound數(shù)據(jù)庫匹配度大于70%,在有機與普通番茄中存在差異的13 個化合物,見表3。在這些化合物中,主要為有機酸、黃酮類化合物,這兩類化合物是番茄中的主要抗氧化成分。與普通番茄組相比較,有機番茄中的小分子有機酸含量降低,黃酮類化合物槲皮素-3-O-三糖苷的含量增加。這與已有研究結(jié)果基本一致,環(huán)境脅迫可激活番茄植株體內(nèi)苯丙醇類化合物的生物合成,有機處理的番茄中類黃酮含量隨種植年限增加而增加,而傳統(tǒng)處理的類黃酮含量并無顯著差異[5,12]。

    表3 初步確定兩組間差異顯著的化合物信息Table 3 Information about significantly differential compounds between organic and common tomatoes

    3 結(jié) 論

    本實驗采用IRMS和HPLC-HRMS分析技術(shù),借助PLS-DA技術(shù)建立一種有機番茄鑒別方法。有機和普通種植方式所生產(chǎn)的番茄,其果實中穩(wěn)定氮同位素比值差異顯著??山Y(jié)合地區(qū)特征,通過穩(wěn)定氮同位素比值大小在一定程度上判定樣品是否為有機番茄。為進一步提高判別效果,將液相色譜高分辨質(zhì)譜分析數(shù)據(jù)與穩(wěn)定碳同位素比值相結(jié)合,構(gòu)建了有機番茄判別模型,篩查出特征變量因子,包括穩(wěn)定氮同位素比值,50多個可能的化學(xué)成分,鑒定出其中的13 個化合物,可作為有機和普通番茄區(qū)別的特征成分。所建立的模型對有機番茄的判別正確率為93.4%,基本可以滿足有機番茄判別的需要。對有機產(chǎn)品的判定,穩(wěn)定氮同位素差異分析更為簡潔方便,但由于有機與普通產(chǎn)品的穩(wěn)定氮同位素比值界限時有重疊,給判定增加了不確定性,因此應(yīng)該進一步結(jié)合液相色譜-串聯(lián)質(zhì)譜的組學(xué)技術(shù),提高判別的可靠性。

    猜你喜歡
    同位素比值質(zhì)譜
    氣相色譜質(zhì)譜聯(lián)用儀在農(nóng)殘檢測中的應(yīng)用及維護
    比值遙感蝕變信息提取及閾值確定(插圖)
    河北遙感(2017年2期)2017-08-07 14:49:00
    吹掃捕集-氣相色譜質(zhì)譜聯(lián)用測定水中18種揮發(fā)性有機物
    不同應(yīng)變率比值計算方法在甲狀腺惡性腫瘤診斷中的應(yīng)用
    深空探測用同位素電源的研究進展
    《同位素》(季刊)2015年征訂通知
    同位素(2014年3期)2014-06-13 08:22:28
    硼同位素分離工藝與生產(chǎn)技術(shù)
    同位素(2014年2期)2014-04-16 04:57:15
    穩(wěn)定同位素氘標(biāo)記蘇丹紅I的同位素豐度和化學(xué)純度分析
    同位素(2014年2期)2014-04-16 04:57:12
    棗霜化學(xué)成分的色譜質(zhì)譜分析
    Modeled response of talik development under thermokarst lakes to permafrost thickness on the Qinghai-Tibet Plateau
    国产精品电影一区二区三区| 老司机影院成人| 国产精品一二三区在线看| 国产成人精品一,二区 | 国产黄色小视频在线观看| 日本一二三区视频观看| 欧美另类亚洲清纯唯美| 国产探花极品一区二区| 国产精品一区二区性色av| 中文字幕av成人在线电影| 欧美一区二区精品小视频在线| 亚洲国产高清在线一区二区三| 变态另类成人亚洲欧美熟女| 成人特级黄色片久久久久久久| 亚洲经典国产精华液单| 成年版毛片免费区| 少妇高潮的动态图| av.在线天堂| 日韩欧美 国产精品| 18禁黄网站禁片免费观看直播| 久久精品久久久久久噜噜老黄 | 欧美日韩一区二区视频在线观看视频在线 | 26uuu在线亚洲综合色| 久久99热这里只有精品18| 亚洲av免费在线观看| 综合色av麻豆| 亚洲欧美日韩卡通动漫| 性色avwww在线观看| 51国产日韩欧美| 男人狂女人下面高潮的视频| 又黄又爽又刺激的免费视频.| 小蜜桃在线观看免费完整版高清| 亚洲av中文av极速乱| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 日韩大尺度精品在线看网址| 日韩精品青青久久久久久| 国产一区二区亚洲精品在线观看| a级毛色黄片| 免费av观看视频| 国产乱人视频| 黄色视频,在线免费观看| 少妇熟女欧美另类| 搞女人的毛片| 日韩精品有码人妻一区| 午夜免费激情av| 在线播放国产精品三级| 欧美最新免费一区二区三区| 亚洲真实伦在线观看| 97超视频在线观看视频| 美女xxoo啪啪120秒动态图| 51国产日韩欧美| 亚洲人成网站在线观看播放| 12—13女人毛片做爰片一| 亚洲欧洲国产日韩| 午夜亚洲福利在线播放| 国产美女午夜福利| 欧美极品一区二区三区四区| 国产欧美日韩精品一区二区| 中文欧美无线码| 超碰av人人做人人爽久久| 欧美色欧美亚洲另类二区| 欧美高清成人免费视频www| 久久欧美精品欧美久久欧美| 一夜夜www| 国产精品蜜桃在线观看 | 亚洲国产欧洲综合997久久,| 秋霞在线观看毛片| 日韩成人伦理影院| 国内揄拍国产精品人妻在线| 夜夜爽天天搞| 日本在线视频免费播放| av在线蜜桃| 久久久久久久久大av| 欧美日韩乱码在线| 日本欧美国产在线视频| 欧美xxxx性猛交bbbb| 乱系列少妇在线播放| 国产69精品久久久久777片| 老司机福利观看| 级片在线观看| 级片在线观看| 日韩大尺度精品在线看网址| 不卡视频在线观看欧美| 中文亚洲av片在线观看爽| 三级国产精品欧美在线观看| 91午夜精品亚洲一区二区三区| 久久99精品国语久久久| 亚洲美女视频黄频| 最近视频中文字幕2019在线8| 久久久久九九精品影院| 非洲黑人性xxxx精品又粗又长| 又爽又黄a免费视频| 亚洲,欧美,日韩| 国产一级毛片七仙女欲春2| 91av网一区二区| 欧美一区二区国产精品久久精品| 国产精品久久久久久av不卡| 久久人人爽人人片av| 久久精品国产亚洲av天美| 美女被艹到高潮喷水动态| 欧美+日韩+精品| 波多野结衣高清无吗| 久久精品国产清高在天天线| 美女xxoo啪啪120秒动态图| 亚洲va在线va天堂va国产| 一进一出抽搐动态| 国产单亲对白刺激| 国产 一区 欧美 日韩| 免费av观看视频| 国产一级毛片七仙女欲春2| 少妇丰满av| 国产精品一区二区性色av| 亚洲精品亚洲一区二区| 国产精品国产三级国产av玫瑰| 精品免费久久久久久久清纯| 99在线视频只有这里精品首页| 国产精品永久免费网站| 亚洲丝袜综合中文字幕| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| 性插视频无遮挡在线免费观看| 久久久久久久久久久丰满| 中文亚洲av片在线观看爽| АⅤ资源中文在线天堂| 国产一级毛片在线| 九九在线视频观看精品| 99九九线精品视频在线观看视频| 蜜桃亚洲精品一区二区三区| 日本黄大片高清| 国国产精品蜜臀av免费| 亚洲av电影不卡..在线观看| 国产免费一级a男人的天堂| 嫩草影院精品99| 一级av片app| 婷婷精品国产亚洲av| 亚洲精品色激情综合| 人妻少妇偷人精品九色| 网址你懂的国产日韩在线| 国产精品国产三级国产av玫瑰| 国产精品久久久久久精品电影小说 | 国产一级毛片在线| 亚洲欧美精品自产自拍| 久久精品综合一区二区三区| 性欧美人与动物交配| 看十八女毛片水多多多| 久久久久久久午夜电影| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 97在线视频观看| 免费看日本二区| 小说图片视频综合网站| 一级毛片久久久久久久久女| 国产 一区 欧美 日韩| 日韩一本色道免费dvd| 欧美日韩乱码在线| 免费观看精品视频网站| 国产一区二区三区av在线 | 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 老司机福利观看| av.在线天堂| 日韩强制内射视频| 国产高清三级在线| av天堂在线播放| 亚洲欧美日韩卡通动漫| 久久久精品大字幕| 少妇人妻一区二区三区视频| 久久午夜亚洲精品久久| 久久久久国产网址| 成年女人永久免费观看视频| 又爽又黄a免费视频| 久久热精品热| 观看美女的网站| 国产av一区在线观看免费| 22中文网久久字幕| 狠狠狠狠99中文字幕| 国产 一区 欧美 日韩| 男人舔奶头视频| av视频在线观看入口| 国产成人精品一,二区 | 免费观看精品视频网站| 日韩精品有码人妻一区| 国产人妻一区二区三区在| 特级一级黄色大片| 99热网站在线观看| 男人狂女人下面高潮的视频| 22中文网久久字幕| 99riav亚洲国产免费| 波多野结衣巨乳人妻| 少妇丰满av| 1024手机看黄色片| 国产黄片美女视频| 99九九线精品视频在线观看视频| 一本精品99久久精品77| 此物有八面人人有两片| 性插视频无遮挡在线免费观看| 美女被艹到高潮喷水动态| 精品久久久久久久久亚洲| 99久久精品国产国产毛片| 欧美高清成人免费视频www| 国产高清激情床上av| 久久久国产成人免费| 国产免费男女视频| 国产一区二区激情短视频| 亚洲av电影不卡..在线观看| 日产精品乱码卡一卡2卡三| 欧美xxxx黑人xx丫x性爽| 少妇的逼水好多| 国产伦精品一区二区三区视频9| 可以在线观看的亚洲视频| 看免费成人av毛片| 国产精品爽爽va在线观看网站| 国内久久婷婷六月综合欲色啪| 22中文网久久字幕| 午夜免费激情av| 欧美日韩国产亚洲二区| 国产精品三级大全| 级片在线观看| 少妇被粗大猛烈的视频| 国产成人aa在线观看| 免费无遮挡裸体视频| 成年女人永久免费观看视频| 国产精品麻豆人妻色哟哟久久 | 国产女主播在线喷水免费视频网站 | 1000部很黄的大片| 哪里可以看免费的av片| 色播亚洲综合网| 秋霞在线观看毛片| 久久久精品大字幕| 久久精品91蜜桃| 亚洲精品成人久久久久久| 欧美一区二区亚洲| 麻豆乱淫一区二区| 国产色婷婷99| 老熟妇乱子伦视频在线观看| 日韩欧美精品v在线| 日韩成人av中文字幕在线观看| 久久久久性生活片| 日韩一区二区视频免费看| 国内精品美女久久久久久| www.色视频.com| 国产精品久久视频播放| 国产av一区在线观看免费| 永久网站在线| 久久精品影院6| 国产精品美女特级片免费视频播放器| 成人欧美大片| 亚洲真实伦在线观看| 嫩草影院精品99| 女人被狂操c到高潮| 免费看光身美女| 亚洲欧美成人综合另类久久久 | 人人妻人人看人人澡| 女同久久另类99精品国产91| av国产免费在线观看| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 内射极品少妇av片p| 哪个播放器可以免费观看大片| 亚洲av熟女| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| 久久99热这里只有精品18| 春色校园在线视频观看| 久久久久久久亚洲中文字幕| 国内精品美女久久久久久| 久久人人爽人人爽人人片va| 人妻少妇偷人精品九色| 久久久久免费精品人妻一区二区| 色综合亚洲欧美另类图片| 国产日本99.免费观看| 高清午夜精品一区二区三区 | 三级经典国产精品| 嫩草影院新地址| 成人午夜精彩视频在线观看| 国产高清视频在线观看网站| 亚洲精品乱码久久久v下载方式| 久久欧美精品欧美久久欧美| 久久午夜亚洲精品久久| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美中文字幕日韩二区| 寂寞人妻少妇视频99o| 少妇猛男粗大的猛烈进出视频 | 欧美zozozo另类| 小蜜桃在线观看免费完整版高清| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 五月伊人婷婷丁香| 国产视频内射| 久久热精品热| 可以在线观看的亚洲视频| 你懂的网址亚洲精品在线观看 | 成熟少妇高潮喷水视频| 91aial.com中文字幕在线观看| 久99久视频精品免费| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产亚洲网站| 日韩精品青青久久久久久| 一本一本综合久久| 不卡一级毛片| 禁无遮挡网站| 日本-黄色视频高清免费观看| 免费搜索国产男女视频| 国产精品乱码一区二三区的特点| 国产精品女同一区二区软件| 欧美日韩国产亚洲二区| 一级二级三级毛片免费看| 免费一级毛片在线播放高清视频| 亚洲av成人av| 精品国内亚洲2022精品成人| 天堂√8在线中文| 在线播放国产精品三级| 精品一区二区免费观看| 国产乱人偷精品视频| 少妇丰满av| 男人狂女人下面高潮的视频| 久久国产乱子免费精品| 如何舔出高潮| 久久精品国产自在天天线| 国产69精品久久久久777片| 91午夜精品亚洲一区二区三区| 久久久色成人| 全区人妻精品视频| 99久久精品国产国产毛片| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 国产午夜精品论理片| 国产视频首页在线观看| 少妇的逼水好多| 麻豆国产97在线/欧美| 久久久a久久爽久久v久久| 久久久成人免费电影| 99久国产av精品| 国产极品天堂在线| 国产视频内射| 国产一区二区三区av在线 | 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 全区人妻精品视频| 国产精品,欧美在线| 日产精品乱码卡一卡2卡三| 日本av手机在线免费观看| 欧美潮喷喷水| 搡女人真爽免费视频火全软件| 国产69精品久久久久777片| 边亲边吃奶的免费视频| 精品久久久久久久末码| 亚洲性久久影院| 精品久久久久久久久av| 日本在线视频免费播放| 自拍偷自拍亚洲精品老妇| 麻豆国产av国片精品| 三级男女做爰猛烈吃奶摸视频| 日本爱情动作片www.在线观看| 成人欧美大片| 丰满人妻一区二区三区视频av| 午夜精品国产一区二区电影 | 国产高清不卡午夜福利| 国产av不卡久久| 精品久久久久久久久久久久久| 欧美日本亚洲视频在线播放| videossex国产| 免费人成在线观看视频色| 91麻豆精品激情在线观看国产| 国产真实乱freesex| 国产黄片美女视频| 一夜夜www| 国内精品久久久久精免费| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 国产精品一区二区性色av| 亚洲精品粉嫩美女一区| 中文字幕免费在线视频6| 国产 一区精品| 成年免费大片在线观看| 黄片无遮挡物在线观看| 午夜福利成人在线免费观看| 久久久成人免费电影| 又爽又黄a免费视频| 91麻豆精品激情在线观看国产| 亚洲天堂国产精品一区在线| 精品久久久久久久久亚洲| 午夜福利成人在线免费观看| 精品一区二区免费观看| 国产成人a区在线观看| 黑人高潮一二区| 日韩视频在线欧美| 国国产精品蜜臀av免费| 国产伦理片在线播放av一区 | 精品久久久噜噜| 又粗又硬又长又爽又黄的视频 | 久久精品国产清高在天天线| 久久久精品大字幕| a级毛片免费高清观看在线播放| 人人妻人人看人人澡| 久久久精品欧美日韩精品| av在线亚洲专区| 嫩草影院新地址| 九草在线视频观看| 在线观看午夜福利视频| 国产成人精品久久久久久| 在现免费观看毛片| 日韩精品青青久久久久久| 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 内地一区二区视频在线| 97超视频在线观看视频| a级一级毛片免费在线观看| 久久久精品94久久精品| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 69人妻影院| 久久久久久久久久久免费av| 亚洲va在线va天堂va国产| 亚洲一区二区三区色噜噜| 国产精品一区www在线观看| 国产三级中文精品| 性欧美人与动物交配| 欧美bdsm另类| 真实男女啪啪啪动态图| 床上黄色一级片| 97超碰精品成人国产| 欧美一区二区亚洲| 网址你懂的国产日韩在线| 午夜老司机福利剧场| 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放| 成人午夜精彩视频在线观看| 亚洲欧美日韩无卡精品| 欧美又色又爽又黄视频| 夜夜爽天天搞| 特级一级黄色大片| 五月玫瑰六月丁香| 可以在线观看的亚洲视频| av.在线天堂| 亚洲国产精品成人综合色| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 亚洲第一电影网av| 欧美+亚洲+日韩+国产| 最近中文字幕高清免费大全6| 国产成人a区在线观看| 男人舔奶头视频| 99国产精品一区二区蜜桃av| 又黄又爽又刺激的免费视频.| 日韩制服骚丝袜av| 亚洲图色成人| 精品久久久久久久久亚洲| 嘟嘟电影网在线观看| 精品一区二区三区人妻视频| 成人无遮挡网站| 亚洲最大成人中文| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 日韩精品有码人妻一区| 中文字幕制服av| 97人妻精品一区二区三区麻豆| 插阴视频在线观看视频| 97超碰精品成人国产| 国产伦精品一区二区三区视频9| 国产老妇伦熟女老妇高清| 欧美+日韩+精品| 日韩欧美三级三区| 久久午夜福利片| 欧美极品一区二区三区四区| 亚洲在久久综合| 美女大奶头视频| 观看美女的网站| 亚洲精品456在线播放app| 国产免费男女视频| 91aial.com中文字幕在线观看| 嫩草影院精品99| 免费在线观看成人毛片| 在线国产一区二区在线| 亚洲精品亚洲一区二区| 欧美性猛交黑人性爽| 欧美xxxx性猛交bbbb| 成人午夜精彩视频在线观看| 在线观看av片永久免费下载| 亚洲精品粉嫩美女一区| 久99久视频精品免费| 国产午夜精品一二区理论片| 亚洲av中文av极速乱| 淫秽高清视频在线观看| 国产视频内射| 国产探花极品一区二区| 日韩中字成人| 精品久久久久久久末码| 精品人妻一区二区三区麻豆| 亚洲国产高清在线一区二区三| 亚洲av免费在线观看| 亚洲国产色片| 麻豆一二三区av精品| 国产av在哪里看| 欧美日韩乱码在线| 中文字幕免费在线视频6| av天堂中文字幕网| 国产一区亚洲一区在线观看| 久久九九热精品免费| av专区在线播放| 只有这里有精品99| 国产黄色小视频在线观看| 亚洲人成网站在线播| 变态另类成人亚洲欧美熟女| 免费不卡的大黄色大毛片视频在线观看 | 免费不卡的大黄色大毛片视频在线观看 | 欧美xxxx黑人xx丫x性爽| 又黄又爽又刺激的免费视频.| 两个人的视频大全免费| 综合色av麻豆| 99久国产av精品| 插逼视频在线观看| 国产黄a三级三级三级人| 一级黄片播放器| 搞女人的毛片| 99久久精品热视频| 欧美色视频一区免费| 亚洲欧美成人精品一区二区| 美女大奶头视频| 在线观看午夜福利视频| 一级毛片aaaaaa免费看小| 只有这里有精品99| 夜夜爽天天搞| 国内精品久久久久精免费| 天天躁夜夜躁狠狠久久av| 午夜亚洲福利在线播放| 国内久久婷婷六月综合欲色啪| 亚洲一区二区三区色噜噜| 黄色视频,在线免费观看| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 欧美性感艳星| 级片在线观看| 波多野结衣高清无吗| 久久韩国三级中文字幕| 国产色婷婷99| 日韩欧美精品v在线| 日韩高清综合在线| 少妇的逼水好多| 午夜a级毛片| 成人av在线播放网站| 一个人看视频在线观看www免费| 亚洲欧洲国产日韩| 麻豆成人午夜福利视频| 久久鲁丝午夜福利片| 又粗又硬又长又爽又黄的视频 | 午夜久久久久精精品| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 男人狂女人下面高潮的视频| 高清午夜精品一区二区三区 | 日本在线视频免费播放| 日韩人妻高清精品专区| 夜夜爽天天搞| 久久久欧美国产精品| 99热网站在线观看| 亚洲成a人片在线一区二区| 99视频精品全部免费 在线| 嫩草影院精品99| 观看免费一级毛片| 精品人妻视频免费看| 精品欧美国产一区二区三| 国内精品一区二区在线观看| 国产av不卡久久| 亚洲人成网站在线播| 欧美+亚洲+日韩+国产| videossex国产| av在线亚洲专区| 深夜精品福利| 中文字幕制服av| 亚洲精品亚洲一区二区| 插阴视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品一区二区三区四区久久| 国产精品蜜桃在线观看 | eeuss影院久久| 中文字幕精品亚洲无线码一区| 免费av不卡在线播放| 一级毛片久久久久久久久女| 亚洲在线观看片| 一区福利在线观看| 青春草国产在线视频 | 黄片无遮挡物在线观看| 国产真实乱freesex| 欧美日韩综合久久久久久| av在线天堂中文字幕| 精品国内亚洲2022精品成人| 99久久精品热视频| 国产一级毛片在线| 此物有八面人人有两片| 亚洲四区av| 小说图片视频综合网站| av视频在线观看入口| kizo精华| 国产一区二区在线av高清观看| 精品久久久久久久久久久久久| 免费观看a级毛片全部| 精品无人区乱码1区二区| 色哟哟哟哟哟哟| 高清毛片免费观看视频网站| 白带黄色成豆腐渣| 日本av手机在线免费观看| 国产成人午夜福利电影在线观看| 欧美精品国产亚洲| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 欧美3d第一页| av天堂在线播放| 熟女人妻精品中文字幕| 精品一区二区三区人妻视频| 一本久久精品| 精品不卡国产一区二区三区| 欧美在线一区亚洲| 亚洲欧美日韩东京热|