• <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免费在线观看| 欧美成狂野欧美在线观看| avwww免费| 成人av一区二区三区在线看| 一级毛片高清免费大全| 国产高清三级在线| 麻豆av在线久日| 91av网站免费观看| 亚洲精品一区av在线观看| 高潮久久久久久久久久久不卡| 热99re8久久精品国产| 国内毛片毛片毛片毛片毛片| 黄片小视频在线播放| 国产av在哪里看| 无限看片的www在线观看| 窝窝影院91人妻| 国产v大片淫在线免费观看| 久久久国产成人免费| 亚洲精品中文字幕一二三四区| 欧美性猛交黑人性爽| 亚洲成人久久性| 国产熟女xx| 日韩欧美在线二视频| 成人av在线播放网站| 男女之事视频高清在线观看| 国产高清视频在线观看网站| av国产免费在线观看| 人妻丰满熟妇av一区二区三区| 99riav亚洲国产免费| 日本免费a在线| 深夜精品福利| 亚洲 国产 在线| 国产毛片a区久久久久| 黄色丝袜av网址大全| 亚洲欧美激情综合另类| 嫁个100分男人电影在线观看| 听说在线观看完整版免费高清| 亚洲成人精品中文字幕电影| 国内精品美女久久久久久| 激情在线观看视频在线高清| 小蜜桃在线观看免费完整版高清| 日本三级黄在线观看| 非洲黑人性xxxx精品又粗又长| a级毛片在线看网站| 免费大片18禁| 久久久久久国产a免费观看| 国产免费av片在线观看野外av| 亚洲精品在线观看二区| 日韩精品中文字幕看吧| 99热精品在线国产| 可以在线观看毛片的网站| 淫妇啪啪啪对白视频| 精品99又大又爽又粗少妇毛片 | 一区福利在线观看| 神马国产精品三级电影在线观看| 亚洲中文字幕日韩| 好看av亚洲va欧美ⅴa在| 国产综合懂色| 欧美国产日韩亚洲一区| 国产成人av教育| 一边摸一边抽搐一进一小说| 偷拍熟女少妇极品色| 老汉色av国产亚洲站长工具| 一二三四社区在线视频社区8| av欧美777| 嫩草影视91久久| 国产三级中文精品| 欧美又色又爽又黄视频| 九九在线视频观看精品| 国产精品99久久久久久久久| 看免费av毛片| 在线播放国产精品三级| av片东京热男人的天堂| 欧美成狂野欧美在线观看| 俺也久久电影网| 国内毛片毛片毛片毛片毛片| 日韩欧美三级三区| 岛国视频午夜一区免费看| 亚洲美女黄片视频| 亚洲真实伦在线观看| 国产精品久久久久久亚洲av鲁大| 国产精品98久久久久久宅男小说| 亚洲午夜理论影院| 久久久久精品国产欧美久久久| 两个人看的免费小视频| 观看免费一级毛片| 国产成人欧美在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲人成电影免费在线| 国产精华一区二区三区| 国产精华一区二区三区| 后天国语完整版免费观看| 精品无人区乱码1区二区| 99久久综合精品五月天人人| 国产在线精品亚洲第一网站| 超碰成人久久| 亚洲va日本ⅴa欧美va伊人久久| 午夜日韩欧美国产| 日日夜夜操网爽| 国产1区2区3区精品| 看片在线看免费视频| 欧美日本视频| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩东京热| 国产精品亚洲美女久久久| 人人妻人人看人人澡| 精品久久久久久久久久免费视频| 国产精品久久电影中文字幕| 最新在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 两性夫妻黄色片| 免费看十八禁软件| 国产精品 国内视频| 91在线精品国自产拍蜜月 | 日本黄色视频三级网站网址| 久久精品亚洲精品国产色婷小说| 亚洲人成网站高清观看| 这个男人来自地球电影免费观看| 久久久久久大精品| 日本精品一区二区三区蜜桃| 19禁男女啪啪无遮挡网站| 成人特级av手机在线观看| 国内毛片毛片毛片毛片毛片| 俺也久久电影网| 国产精品亚洲av一区麻豆| 国产精品 欧美亚洲| 最新美女视频免费是黄的| 欧美最黄视频在线播放免费| 欧美乱色亚洲激情| 亚洲人成网站在线播放欧美日韩| 99视频精品全部免费 在线 | 免费观看人在逋| 午夜免费激情av| 男人和女人高潮做爰伦理| 精品久久久久久久末码| 亚洲成av人片免费观看| 亚洲人成伊人成综合网2020| 91麻豆av在线| 伦理电影免费视频| 亚洲欧美日韩高清在线视频| 亚洲无线观看免费| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美网| 国产精品99久久久久久久久| 亚洲成av人片在线播放无| 亚洲成av人片在线播放无| 麻豆一二三区av精品| 国产伦精品一区二区三区四那| 国产午夜精品论理片| 少妇的逼水好多| h日本视频在线播放| 久久精品91蜜桃| 村上凉子中文字幕在线| 久久久成人免费电影| 日韩成人在线观看一区二区三区| 麻豆成人av在线观看| 日韩欧美 国产精品| 嫩草影视91久久| avwww免费| 深夜精品福利| 精品久久久久久久毛片微露脸| 欧美成人性av电影在线观看| 中国美女看黄片| 可以在线观看的亚洲视频| 久久久久国产一级毛片高清牌| 女人被狂操c到高潮| 国产久久久一区二区三区| 人妻丰满熟妇av一区二区三区| 欧美乱码精品一区二区三区| 亚洲熟女毛片儿| 亚洲电影在线观看av| 999久久久精品免费观看国产| 亚洲专区中文字幕在线| 嫩草影视91久久| 久久精品aⅴ一区二区三区四区| 国产一区二区在线av高清观看| 欧美激情在线99| 国产精华一区二区三区| 久久精品国产99精品国产亚洲性色| 欧美中文日本在线观看视频| 亚洲一区高清亚洲精品| 99久久综合精品五月天人人| 成人三级做爰电影| 久99久视频精品免费| 夜夜爽天天搞| 亚洲成av人片免费观看| 亚洲欧美日韩高清专用| 国内毛片毛片毛片毛片毛片| 99国产综合亚洲精品| 午夜福利免费观看在线| 欧美zozozo另类| 欧美日韩一级在线毛片| 最近在线观看免费完整版| 最新在线观看一区二区三区| 国产久久久一区二区三区| 黑人操中国人逼视频| 一区二区三区国产精品乱码| 国产亚洲精品一区二区www| 成年版毛片免费区| 999久久久精品免费观看国产| 男人舔女人的私密视频| 国产av在哪里看| 黄色女人牲交| 欧美黑人欧美精品刺激| 九色成人免费人妻av| 听说在线观看完整版免费高清| 欧美日韩乱码在线| 夜夜爽天天搞| 美女cb高潮喷水在线观看 | 国产成人av激情在线播放| 真人做人爱边吃奶动态| 网址你懂的国产日韩在线| 一进一出抽搐gif免费好疼| 成人三级黄色视频| 免费高清视频大片| 国内毛片毛片毛片毛片毛片| 99久久无色码亚洲精品果冻| 熟女人妻精品中文字幕| 最新美女视频免费是黄的| 少妇的逼水好多| 曰老女人黄片| 亚洲中文av在线| 国产精品av久久久久免费| 成人无遮挡网站| 日本撒尿小便嘘嘘汇集6| 日本在线视频免费播放| 国产在线精品亚洲第一网站| 一本综合久久免费| 国产毛片a区久久久久| 三级男女做爰猛烈吃奶摸视频| 嫩草影视91久久| 久久精品夜夜夜夜夜久久蜜豆| 久久天堂一区二区三区四区| 免费观看人在逋| a级毛片在线看网站| xxxwww97欧美| 精品久久久久久久久久久久久| 搞女人的毛片| 我的老师免费观看完整版| 国产主播在线观看一区二区| 一个人免费在线观看的高清视频| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 国产伦一二天堂av在线观看| 夜夜看夜夜爽夜夜摸| 国产精品免费一区二区三区在线| 午夜福利视频1000在线观看| 免费观看精品视频网站| 午夜两性在线视频| 天堂√8在线中文| 亚洲人与动物交配视频| 淫妇啪啪啪对白视频| ponron亚洲| 国产精品一及| 午夜福利视频1000在线观看| 精品久久久久久久人妻蜜臀av| 午夜免费成人在线视频| 天天添夜夜摸| xxxwww97欧美| 欧美绝顶高潮抽搐喷水| 草草在线视频免费看| 亚洲七黄色美女视频| 国产人伦9x9x在线观看| 国内精品美女久久久久久| 90打野战视频偷拍视频| 女人被狂操c到高潮| 91九色精品人成在线观看| 老熟妇仑乱视频hdxx| 国产av不卡久久| 精品熟女少妇八av免费久了| 看黄色毛片网站| 男女午夜视频在线观看| 观看美女的网站| 少妇的丰满在线观看| 中文资源天堂在线| 久久国产乱子伦精品免费另类| 久久久精品欧美日韩精品| x7x7x7水蜜桃| 国产亚洲精品av在线| 好男人电影高清在线观看| 亚洲熟妇中文字幕五十中出| 亚洲五月婷婷丁香| 亚洲欧洲精品一区二区精品久久久| 一级黄色大片毛片| 在线免费观看的www视频| 欧美日韩福利视频一区二区| 99久久99久久久精品蜜桃| 欧美在线黄色| 色精品久久人妻99蜜桃| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| 免费av毛片视频| 亚洲精品中文字幕一二三四区| 校园春色视频在线观看| 国产精品国产高清国产av| 色噜噜av男人的天堂激情| 搡老妇女老女人老熟妇| 夜夜夜夜夜久久久久| 久久精品亚洲精品国产色婷小说| 亚洲自拍偷在线| 成人午夜高清在线视频| 99久久综合精品五月天人人| 久久久久国产精品人妻aⅴ院| 国内揄拍国产精品人妻在线| 十八禁人妻一区二区| 中文字幕久久专区| 日韩精品青青久久久久久| 成人特级av手机在线观看| 亚洲人成伊人成综合网2020| 国产毛片a区久久久久| 长腿黑丝高跟| 桃红色精品国产亚洲av| 女警被强在线播放| 国产探花在线观看一区二区| 国产精品1区2区在线观看.| 成人三级黄色视频| 天堂动漫精品| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| ponron亚洲| 波多野结衣高清作品| 成人性生交大片免费视频hd| 小蜜桃在线观看免费完整版高清| 最好的美女福利视频网| 听说在线观看完整版免费高清| 女警被强在线播放| 亚洲 欧美一区二区三区| 全区人妻精品视频| 黄片小视频在线播放| 精品久久久久久久末码| 欧美日韩亚洲国产一区二区在线观看| 成人高潮视频无遮挡免费网站| 精品国产美女av久久久久小说| 欧美黑人巨大hd| 亚洲第一电影网av| 久久国产精品影院| 久久中文字幕人妻熟女| 九色国产91popny在线| 男女午夜视频在线观看| 后天国语完整版免费观看| 小蜜桃在线观看免费完整版高清| 18禁国产床啪视频网站| 国内精品久久久久久久电影| 国内精品一区二区在线观看| 国产一区二区激情短视频| 欧美极品一区二区三区四区| 悠悠久久av| 嫩草影院入口| 国产亚洲精品av在线| 免费电影在线观看免费观看| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 精品国产乱码久久久久久男人| 久久人妻av系列| 成年免费大片在线观看| 日韩欧美在线乱码| 国产一区在线观看成人免费| 精品一区二区三区视频在线 | 级片在线观看| 人人妻人人澡欧美一区二区| 国产成人影院久久av| 好看av亚洲va欧美ⅴa在| 久久亚洲精品不卡| 三级毛片av免费| 国产精品1区2区在线观看.| 九九热线精品视视频播放| 99久久成人亚洲精品观看| www.自偷自拍.com| 久久精品91无色码中文字幕| 人人妻人人看人人澡| 他把我摸到了高潮在线观看| av片东京热男人的天堂| 国产99白浆流出| 岛国在线免费视频观看| 国产精品野战在线观看| 久久99热这里只有精品18| x7x7x7水蜜桃| 中文亚洲av片在线观看爽| 制服丝袜大香蕉在线| 中文字幕久久专区| 精品国产三级普通话版| 很黄的视频免费| 亚洲电影在线观看av| 亚洲七黄色美女视频| 国产av麻豆久久久久久久| 久久精品国产亚洲av香蕉五月| 亚洲成人久久爱视频| 变态另类成人亚洲欧美熟女| 欧美性猛交╳xxx乱大交人| 日本一本二区三区精品| 亚洲av成人不卡在线观看播放网| 亚洲美女视频黄频| 黄片大片在线免费观看| 国产熟女xx| 在线免费观看的www视频| 亚洲国产中文字幕在线视频| 欧美又色又爽又黄视频| 嫁个100分男人电影在线观看| 成年版毛片免费区| 淫秽高清视频在线观看| 男女那种视频在线观看| 欧美av亚洲av综合av国产av| 日本 欧美在线| 无限看片的www在线观看| 国产黄a三级三级三级人| 法律面前人人平等表现在哪些方面| 麻豆久久精品国产亚洲av| 久久久久久久精品吃奶| 小蜜桃在线观看免费完整版高清| 91麻豆精品激情在线观看国产| 成年女人毛片免费观看观看9| 日韩人妻高清精品专区| 亚洲欧美日韩卡通动漫| 国产亚洲精品av在线| 老汉色av国产亚洲站长工具| 亚洲精品乱码久久久v下载方式 | 啦啦啦免费观看视频1| bbb黄色大片| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕大全免费视频| 日日夜夜操网爽| 美女黄网站色视频| 成年女人毛片免费观看观看9| 久久久久久九九精品二区国产| 18美女黄网站色大片免费观看| 欧美乱色亚洲激情| 成年女人永久免费观看视频| 日韩欧美国产在线观看| 啦啦啦免费观看视频1| 日本撒尿小便嘘嘘汇集6| 精品国产超薄肉色丝袜足j| 69av精品久久久久久| 久久久久免费精品人妻一区二区| 日本精品一区二区三区蜜桃| 国产成人欧美在线观看| 男人舔奶头视频| 1024手机看黄色片| 亚洲欧美精品综合一区二区三区| 亚洲欧美激情综合另类| 少妇的丰满在线观看| 黄色日韩在线| 中文字幕熟女人妻在线| 天堂网av新在线| 热99在线观看视频| 久久久久久久久免费视频了| 午夜福利在线在线| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 两个人看的免费小视频| 久久人人精品亚洲av| 一级作爱视频免费观看| 免费看光身美女| 免费av毛片视频| 熟女少妇亚洲综合色aaa.| 伊人久久大香线蕉亚洲五| 国产精华一区二区三区| 色老头精品视频在线观看| 综合色av麻豆| 脱女人内裤的视频| 中文字幕人成人乱码亚洲影| 国产成人一区二区三区免费视频网站| 精品不卡国产一区二区三区| 成在线人永久免费视频| 国产av麻豆久久久久久久| 久久久水蜜桃国产精品网| 黄色丝袜av网址大全| 日韩欧美精品v在线| 国产成人一区二区三区免费视频网站| 成年女人看的毛片在线观看| 啪啪无遮挡十八禁网站| 一个人看视频在线观看www免费 | 国产不卡一卡二| 亚洲一区高清亚洲精品| 国产精品影院久久| 欧美日韩一级在线毛片| 国产99白浆流出| 久久久久九九精品影院| 日本成人三级电影网站| 最近最新免费中文字幕在线| 国产精品国产高清国产av| 亚洲av第一区精品v没综合| 国产伦在线观看视频一区| 舔av片在线| 亚洲国产欧美网| 90打野战视频偷拍视频| 国产激情偷乱视频一区二区| 男女床上黄色一级片免费看| 制服人妻中文乱码| 男女下面进入的视频免费午夜| 91麻豆av在线| 成人三级黄色视频| 国产一区二区激情短视频| 国产激情偷乱视频一区二区| 久久久久久久精品吃奶| 免费在线观看影片大全网站| 久久草成人影院| 97超视频在线观看视频| 亚洲午夜理论影院| 欧美xxxx黑人xx丫x性爽| 国产爱豆传媒在线观看| 欧美一级毛片孕妇| 天堂动漫精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产亚洲av嫩草精品影院| 婷婷精品国产亚洲av在线| 欧美中文综合在线视频| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| tocl精华| 18禁黄网站禁片免费观看直播| 久久国产精品人妻蜜桃| 俄罗斯特黄特色一大片| 免费av不卡在线播放| 18禁黄网站禁片午夜丰满| 色噜噜av男人的天堂激情| 怎么达到女性高潮| 91在线精品国自产拍蜜月 | 精品国产美女av久久久久小说| 亚洲av中文字字幕乱码综合| 久久亚洲真实| 国产乱人视频| 国产成人啪精品午夜网站| 国产一区二区在线av高清观看| 国内精品一区二区在线观看| 国产淫片久久久久久久久 | 亚洲va日本ⅴa欧美va伊人久久| 禁无遮挡网站| 久久久国产欧美日韩av| 亚洲人成网站高清观看| 国产黄色小视频在线观看| 国产极品精品免费视频能看的| 桃色一区二区三区在线观看| 亚洲,欧美精品.| 韩国av一区二区三区四区| 久久久久国内视频| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 国产午夜精品久久久久久| 婷婷亚洲欧美| 国产精品美女特级片免费视频播放器 | 免费电影在线观看免费观看| 亚洲无线观看免费| 免费看十八禁软件| 日本黄大片高清| bbb黄色大片| 国产三级在线视频| 亚洲第一欧美日韩一区二区三区| 午夜福利欧美成人| 亚洲成a人片在线一区二区| 2021天堂中文幕一二区在线观| 中出人妻视频一区二区| 国产伦精品一区二区三区视频9 | av女优亚洲男人天堂 | 黄色片一级片一级黄色片| 女人高潮潮喷娇喘18禁视频| 亚洲第一欧美日韩一区二区三区| 男女下面进入的视频免费午夜| 精品久久久久久久久久久久久| 久久久久久久久久黄片| 成人一区二区视频在线观看| 露出奶头的视频| 两性午夜刺激爽爽歪歪视频在线观看| 免费人成视频x8x8入口观看| 国产精品99久久99久久久不卡| 男人的好看免费观看在线视频| 两个人的视频大全免费| 亚洲人成伊人成综合网2020| 精品久久久久久,| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 久久人妻av系列| 国产亚洲欧美在线一区二区| 在线观看舔阴道视频| 精品免费久久久久久久清纯| 超碰成人久久| 在线观看免费午夜福利视频| 国产一区二区在线观看日韩 | 国产毛片a区久久久久| 一进一出抽搐动态| 91麻豆精品激情在线观看国产| 国产伦在线观看视频一区| 后天国语完整版免费观看| 国产精品香港三级国产av潘金莲| 色综合婷婷激情| 男插女下体视频免费在线播放| 亚洲成人中文字幕在线播放| 麻豆一二三区av精品| 99riav亚洲国产免费| 少妇人妻一区二区三区视频| 精品国内亚洲2022精品成人| 免费搜索国产男女视频| av中文乱码字幕在线| 天堂网av新在线| 啪啪无遮挡十八禁网站| 国产av在哪里看| 国内精品久久久久精免费| 天天躁日日操中文字幕| av天堂在线播放| 天堂影院成人在线观看| 亚洲精品中文字幕一二三四区| 男人舔女人的私密视频| 久久人妻av系列| 久久久久久久精品吃奶| 婷婷精品国产亚洲av在线| 欧美一级a爱片免费观看看| 母亲3免费完整高清在线观看| 亚洲精品在线观看二区| 高潮久久久久久久久久久不卡| 老鸭窝网址在线观看| 精品人妻1区二区| 免费观看精品视频网站| 成人无遮挡网站| 亚洲成人中文字幕在线播放|