999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

不同產地絨柄牛肝菌紫外指紋圖譜鑒別分析

2015-02-16 01:54:10楊天偉李杰慶劉鴻高王元忠
食品工業科技 2015年9期

楊天偉,李 濤,張 霽,李杰慶,劉鴻高,王元忠,*

(1.云南農業大學農學與生物技術學院,云南昆明 650201;2.云南省農業科學院藥用植物研究所,云南昆明 650200;3.玉溪師范學院資源環境學院,云南玉溪 653100)

?

不同產地絨柄牛肝菌紫外指紋圖譜鑒別分析

楊天偉1,2,李 濤3,張 霽2,李杰慶1,劉鴻高1,王元忠2,*

(1.云南農業大學農學與生物技術學院,云南昆明 650201;2.云南省農業科學院藥用植物研究所,云南昆明 650200;3.玉溪師范學院資源環境學院,云南玉溪 653100)

應用紫外光譜技術建立了云南9個不同地區絨柄牛肝菌的紫外指紋圖譜,采用歐氏距離和主成分分析法對230~450nm的紫外光譜數據進行分析。結果顯示該方法的精密度、重現性及10h內穩定性的RSD分別在0~2.87%、0.03%~0.63%、0.04%~1.73%之間;不同產地樣品間的歐氏距離值在0.26~6.52之間,樣品間的歐氏距離明顯;主成分分析的前三個主成分累積貢獻率達到86.848%,能夠表達樣品主要信息,前兩個主成分的二維投影圖能夠較好地區分不同產地絨柄牛肝菌樣品。紫外光譜結合歐氏距離、主成分分析法能夠快速鑒別不同產地絨柄牛肝菌。

紫外光譜,歐氏距離,主成分分析,絨柄牛肝菌,鑒別

絨柄牛肝菌(BoletustomentipesEarle)又名黑牛肝、大巴菌、毛腳牛肝菌,是牛肝菌科(Boletaceae)的一種野生食用菌[1],其子實體富含蛋白質、氨基酸、多糖、甘露醇、脂肪及多種礦質元素[2-3],其中多糖等有增強免疫、抗癌、抗病毒等功能,具有極大的食藥用價值[4-6],在我國主要分布于云南[6]。云南特殊的生境資源,孕育了豐富的大型真菌資源;野生牛肝菌種類繁多,分布廣泛,受地理因素影響,不同產地野生食用菌,化學成分、營養物質積累、藥效等不盡相同[7-8];對其進行產地鑒別是食用菌開發利用和深入研究的基礎,同時有助于野生食用菌質量控制。傳統的野生食用菌的分類主要根據食用菌子實體的外觀形貌、顯微結構、生長特性等進行鑒別[9],但這種方法難以鑒別出產地來源。目前,對食用菌鑒別分類的研究主要用紅外光譜法;周在進等[10]用傅里葉變換紅外光譜(FTIR)技術鑒別了不同產地雙色牛肝菌;孫素琴等[11]用傅里葉變換紅外光譜鑒別了不同廠家生產的36種靈芝產品;劉剛等[9]根據青頭菌、大紅菇等不同食用菌紅外吸收光譜的峰形、峰高對其進行鑒別,但這種方法需要豐富的經驗,難以建立完整的鑒別體系。

紫外光譜法具有簡便快速、靈敏可靠的特點,根據不同物質體系的紫外吸收圖譜峰形、峰高、峰面積等差異,可以鑒別不同物質體系[12-13]。紫外指紋圖譜技術已被廣泛用于中藥材鑒別及含量測定[14-15],食品摻假鑒別及檢測[16-17]等多個領域。本研究采用紫外光譜技術結合歐氏距離和主成分分析法對9個采自云南不同產地的絨柄牛肝菌進行鑒別分析,該方法能夠快速鑒別出不同產地絨柄牛肝菌,為野生食用菌產地鑒別和質量控制提供參考。

1 材料與方法

1.1 材料與儀器

絨柄牛肝菌樣品均采自2011年,并由云南農業大學劉鴻高教授鑒定為絨柄牛肝菌BoletustomentipesEarle,來源見表1。

表1 絨柄牛肝菌樣品來源Table1 Sources of Boletus tomentipes samples

氯仿(分析純) 西隴化工股份有限公司;石油醚(分析純) 天津市風船化學試劑科技有限公司;氫氧化鈉(分析純) 西隴化工股份有限公司;無水乙醇(分析純) 云南汕滇藥業有限公司;蒸餾水。

UV-2550雙通道紫外可見分光光度計配有UV probe工作站 日本島津公司;KQ5200型超聲波清洗機 昆明市超聲儀器有限公司;FW-100型高速粉碎機 天津市華鑫儀器廠;100目標準篩盤 浙江上虞市道墟五四儀器廠。

1.2 樣品處理及紫外光譜測定

樣品采集后,清洗干凈,50℃烘干,取同一產地的10個絨柄牛肝菌子實體混合粉碎,過100目篩,備用;準確稱取0.10g絨柄牛肝菌樣品于25mL比色管中,加入10mL氯仿,超聲提取30min,過濾得氯仿提取液。以氯仿作為參比液和測試液進行基線校正及空白測定,以消除基線漂移;以氯仿為參比液測定牛肝菌樣品提取液的紫外指紋圖譜,設定掃描波長為190~450nm,重復掃描3次,狹縫1.0nm,采樣間隔0.2nm。對氯仿提取液所測得的原始光譜進行3組總平均,減去空白,以消除溶劑干擾,提高光譜測定的準確度[18]。

1.3 牛肝菌樣品的提取條件

1.3.1 最佳提取溶劑 以樣品B4為考察對象,準確稱取0.10g樣品于25mL比色管中,分別加入10mL蒸餾水、無水乙醇、氯仿、石油醚和0.5mol/L氫氧化鈉,分別平行3次;超聲提取30min,過濾,以對應溶劑為參比液測定紫外光譜,根據圖譜吸收峰數考察不同溶劑對絨柄牛肝菌提取率的影響。

1.3.2 最佳提取時間 以樣品B4為考察對象,準確稱取0.10g樣品,加入10mL氯仿,分別超聲提取20、30、40、50、60min,過濾,提取液進行紫外光譜測定,根據吸收峰數確定最佳提取時間。

1.4 方法學實驗

以樣品B4為考察對象,稱取0.10g樣品5份,分別加入10mL氯仿,超聲提取30min過濾,經190~450nm紫外光譜測定,所得光譜數據轉置后計算共有峰吸收波長的平均相對標準偏差(RSD%),考察重現性;其中共有峰的確定參照鄒華彬的共有峰識別方法[19]。

樣品B4的氯仿提取液在190~450nm重復測定6次,計算吸收波長的RSD%,考察精密度。

樣品B4的氯仿提取液分別在2、4、6、8、10h時進行紫外光譜測定,計算吸收波長的RSD%,考察穩定性。

1.5 數據處理

由于190~230nm紫外光譜的干擾嚴重,而絨柄牛肝菌紫外光譜的特征吸收峰在230~450nm波長范圍內,因此將230~450nm波長范圍內的絨柄牛肝菌紫外光譜數據(波長、吸光度)轉置后,用SPSS17.0軟件計算絨柄牛肝菌樣品間的歐氏距離及進行主成分分析,以直觀表征不同產地絨柄牛肝菌樣品間的相似性。

2 結果與分析

2.1 絨柄牛肝菌提取條件

2.1.1 最佳提取溶劑 不同溶劑提取絨柄牛肝菌的紫外光譜見圖1,由圖1可以看出,蒸餾水、乙醇和0.5mol/L氫氧化鈉提取液幾乎沒有紫外吸收峰;石油醚提取液吸收峰數少且吸收峰不明顯;用氯仿提取的絨柄牛肝菌在260~300nm波長范圍有明顯的吸收峰,吸收峰數多于其它溶劑提取的吸收峰數,因此用氯仿作為提取溶劑。

圖1 不同溶劑提取絨柄牛肝菌的紫外指紋圖譜Fig. 1 UV fingerprint spectra of Boletus tomentipes extracted by different solvents

表2 實驗方法的精密度Table2 The precision of the method

注:表中的同一列為一組共有峰的吸收波長,下同。

表3 實驗方法的重現性Table3 The repeatability of experimental methods

表4 實驗方法的穩定性Table4 The stability of the experimental method

2.1.2 最佳提取時間 樣品B4不同提取時間的紫外光譜見圖2,由圖2可看出提取時間為20min時絨柄牛肝菌提取液的紫外吸收峰不明顯;當超聲提取時間為30min時出現了明顯的紫外吸收峰,此后延長提取時間對絨柄牛肝菌的紫外指紋圖譜影響不明顯,故選30min為提取時間。

圖2 絨柄牛肝菌不同提取時間的紫外指紋圖譜Fig. 2 UV fingerprint spectra of Boletus tomentipes with different extraction times

2.1.3 方法學考察 以樣品B4為方法學考察對象,其結果見表2~表4。由表2可知樣品B4提取液重復測定6次,用6次測定的共有峰的波長數據計算平均相對標準偏差(RSD)在0~2.87%之間,表明該方法精密度好。稱取5份樣品,分別提取、測定紫外光譜,以共有峰波長數據計算的RSD介于0.03%~0.63%之間,表明該方法重現性好。將樣品B4提取液在不同時間段測定后,以共有峰吸收波長計算的RSD在0.04%~1.73%之間,表明樣品制備液在10h內穩定。

2.2 不同產地絨柄牛肝菌紫外指紋圖譜

圖3為9個不同產地絨柄牛肝菌紫外指紋圖譜,由圖可以看出不同產地絨柄牛肝菌紫外指紋圖譜較為相似,276、281、295nm等是樣品的共有峰,表明不同產地絨柄牛肝菌主要的化學組分相似;但不同產地樣品的紫外指紋圖譜的吸收峰數目、峰高、峰位等存在一定的差異,具有明顯的指紋特性,如采自南華天申堂的絨柄牛肝菌(圖中B5)紫外吸收強度明顯高于其他產地樣品的紫外吸光度,而采自曲靖澤州桂花樹的樣品(圖中B9)是所有樣品中吸光度最小的。

圖3 不同產地絨柄牛肝菌紫外指紋圖譜Fig. 3 UV fingerprint spectra of Boletus tomentipes from different regions

2.3 不同產地絨柄牛肝菌歐氏距離

歐氏距離能用來描述樣品間的相似或不相似程度,距離越小,說明兩樣品越相似,距離越大兩樣品的差別越大[20]。根據不同產地絨柄牛肝菌在230~450nm的紫外吸光度計算樣品間的歐氏距離,結果見表5。由表5可知,不同產地絨柄牛肝菌樣品間距離最大的是B5∶B9(B5∶B9為樣品B5與樣品B9之間的歐氏距離),為6.52,表明兩樣品間相似度最小;采自楚雄南華天申堂的樣品B5及采自曲靖澤州桂花樹的樣品B9與其他產地樣品間的歐氏距離比較明顯,表明這兩個產地的絨柄牛肝菌與其他產地絨柄牛肝菌最易于區分;9個不同產地絨柄牛肝菌間的歐氏距離介于0.26~6.52之間,表明同一種牛肝菌在不同的生長環境下營養成分及含量的積累具有差別,這與絨柄牛肝菌紫外指紋圖譜反映的信息一致。歐氏距離能夠直觀地反映出不同產地樣品間的差異性,可用于不同產地絨柄牛肝菌鑒別。

表5 不同產地絨柄牛肝菌樣品間的歐氏距離Table5 The euclidean distance among differentBoletus tomentipes samples

2.4 不同產地絨柄牛肝菌主成分分析

主成分分析(PCA)是以降維的思想將多個指標綜合為少數幾個指標,用幾個綜合因子替代多個原始變量,使綜合因子盡可能的反映原始變量的信息,以簡化分析過程[20]。對9個不同產地絨柄牛肝菌進行主成分分析,前三個主成分的特征值和累積貢獻率見表6,由表6可知主成分1占總方差貢獻率的65.120%,是絨柄牛肝菌最重要的成分,主成分2、主成分3分別占總方差貢獻率的11.741%和9.987%。前三個主成分的累積貢獻率達到86.848%,能夠反映絨柄牛肝菌紫外指紋圖譜全部信息的86.848%(大于85%),僅有13.152%的信息丟失,能用于樣品間的相似性分析。圖4為前三個主成分的三維投影圖,由圖4可以看出不同產地絨柄牛肝菌樣品間的相似程度,其中樣品B9與其他樣品的相似度最小;樣品B1、B6、B7、B8之間的相似度較高,從三維圖中不易區分所有樣品。比較發現用PC1、PC2構成的二維圖(圖5)能更好的區分不同產地絨柄牛肝菌,從圖5可以看出不同產地絨柄牛肝菌樣品的主成分圖產生明顯的離散現象,反映出不同產地樣品的化學成分積累不相同。

表6 主成分特征值及貢獻率Table6 Characteristic value and contributing ratios of principal components

圖4 不同產地絨柄牛肝菌主成分分析的三維圖Fig.4 Three-dimensional diagram of Boletus tomentipes from different regions by PCA

圖5 不同產地絨柄牛肝菌樣品主成分分析的二維圖Fig.5 Two-dimensional diagram of Boletus tomentipes from different regions by PCA

3 討論

絨柄牛肝菌在不同生長環境或不同生長階段會出現大量的表形變異,如其菌蓋有暗褐色和棕褐色;幼嫩時邊緣內卷,老后邊緣一般變平[1];此外,受生長環境影響,不同產地野生食用菌的化學成分、營養價值等不盡相同[21],傳統的形態學鑒定方法難以鑒別產地來源[9]。

本文采用紫外光譜結合歐氏距離和主成分分析法對不同產地絨柄牛肝菌進行鑒別研究,結果顯示不同產地絨柄牛肝菌紫外指紋圖譜的峰形比較相似,表明絨柄牛肝菌的氯仿提取液中主要化學成分比較相似,即絨柄牛肝菌積累的化學組分基本相同;而絨柄牛肝菌氯仿提取液在同一波長處紫外吸光度差異較明顯,表明不同生長環境下絨柄牛肝菌對營養物質的積累量不同,這與Svoboda和Chrastny[22],Falandysz[21]等研究的不同生態環境下食用菌中富集的微量元素有差異的結論相符。

Chojnacka等測定了12個地區絨蓋牛肝菌屬(Xerocomus subtomentosus)及相應土壤中Hg元素含量,發現不同地區絨蓋牛肝菌屬對Hg元素的積累不同,采自礦產周圍的樣品Hg元素含量最高,同時發現絨蓋牛肝菌屬中Hg的含量與土壤中Hg的含量呈負相關關系[23]。本實驗將不同產地絨柄牛肝菌紫外光譜數據轉置后進行歐氏距離和主成分分析,結果顯示,樣品間的歐氏距離較明顯,其中楚雄南華天申堂的樣品B5及采自曲靖澤州桂花樹的樣品B9間的歐氏距離最大,這與圖3反映的結果一致,表明這兩個產地的樣品對化學成分的積累差異較大。主成分分析結果顯示,前三個主成分能夠反映絨柄牛肝菌樣品的主要信息,9個產地樣品在主成分分析的二維圖中產生離散分布,表明每個產地絨柄牛肝菌對化學成分的積累具有差異。這可能是由于不同產地氣候條件、土壤環境、及地質地貌等的差異,導致絨柄牛肝菌對營養物質的積累具有差異性和特殊性。

4 結論

應用紫外光譜結合歐氏距離和主成分分析法分析了9個不同產地絨柄牛肝菌樣品。首先通過單因素實驗考察不同溶劑,不同提取時間對絨柄牛肝菌提取效果的影響,確定提取絨柄牛肝菌特征成分的溶劑為氯仿,超聲提取時間為30min;然后建立不同產地絨柄牛肝菌紫外指紋圖譜。不同產地絨柄牛肝菌紫外指紋圖譜較為相似,而紫外圖譜的峰強、峰位等存在一定的差異,表明絨柄牛肝菌的主要化學組分相似,但含量不同。將230~450nm的紫外光譜數據轉置后,分別用歐氏距離和主成分分析法分析不同產地絨柄牛肝菌樣品的相似度,結果顯示兩種分析方法均能反映絨柄牛肝菌樣品間的差異性;其中,樣品間的歐氏距離在0.26~6.52之間,根據歐氏距離大小,可以反映樣品間相似度的從而鑒別不同產地絨柄牛肝菌;主成分分析的前三個主成分能夠反映樣品的主要信息,前兩個主成分的二維投影圖產生離散現象。紫外光譜技術結合歐氏距離和主成分分析法能夠鑒別不同產地絨柄牛肝菌,為野生食用菌的產地鑒別和質量控制提供參考。

[1]卯曉嵐. 中國大型真菌[M]. 河南:科學技術出版社,2000,322.

[2]李濤,王元忠,李榮春. 絨柄牛肝菌子實體營養成分分析[J]. 食用菌學報,2009,15(4):70-72.

[3]李濤,王元忠,張霽,等. 絨柄牛肝菌微量元素含量測定及相關性分析[J]. 光譜學與光譜分析,2011,31(7):1939-1941.

[4]倪宗耀. 食用菌的食療作用與保健功效[J]. 食藥用菌,2013,21(1):22-25.

[5]Yaltirak T,Aslim B,Ozturk S,et al. Antimicrobial and antioxidant activities ofRussuladelicaFr[J]. Food and Chemical Toxicology,2009,47(8):2052-2056.

[6]趙云霞,陶明煊,程光宇,等. 黑牛肝菌多糖對小鼠酒精性腎損傷保護作用的研究[J]. 食品工業科技,2013,34(20):368-371.

[7]李泰輝,宋斌. 中國食用牛肝菌的種類及其分布[J]. 食用菌學報,2002,9(2):22-30.

[8]李樹紅,趙永昌,于富強,等. 云南商品牛肝菌中易混淆毒牛肝系統學研究[J]. 中國食用菌,2011,30(5):34-36.

[9]劉剛,劉劍虹,楊愛明,等. 食用菌的傅里葉變換紅外光譜鑒別[J]. 光譜學與光譜分析,2004,24(8):941-945.

[10]周在進,劉剛,任先培. 不同產地小美牛肝菌的紅外光譜聚類分析研究[J]. 激光與紅外,2009,39(11):1158-1162.

[11]孫素琴,杜德國,梁曦云,等. 36種靈芝產品傅里葉變換紅外光譜快速鑒別研究[J]. 分析化學,2001,29(3):309-312.

[12]Casale M,Oliveri P,Armanino C,et al. NIR and UV-vis spectroscopy,artificial nose and tongue:Comparison of four fingerprinting techniques for the characterisation of Italian red wines[J]. Analytica chimica acta,2010,668(2):143-148.

[13]賈萍. 指紋圖譜與中藥質量控制[J]. 中國藥業,2004,13(2):79-80.

[14]Baciu A,Ranga F,Fetea F,et al. Fingerprinting food supplements and their botanical ingredients by coupled UV/Vis/FTIR spectrometry[J]. Bulletin of university of agricultural sciences and veterinary medicine Cluj-Napoca. Food Science and Technology,2013,70(1):8-15.

[15]查圣華,姜水紅,張宏. 燕窩及其偽品鑒定方法研究[J]. 食品科技,2010(4):281-283.

[16]Dost K,Ideli C. Determination of polycyclic aromatic hydrocarbons in edible oils and barbecued food by HPLC/UV-Vis detection[J]. Food Chemistry,2012,133(1):193-199.

[17]楊晨,于修燭,王昕,等. 基于紫外光譜的花生油摻偽檢測[J]. 食品科學,2012,33(6):186-189.

[18]袁天軍,王元忠,趙艷麗,等. 滇龍膽紫外指紋圖譜共有峰率和變異峰率雙指標序列分析法[J]. 光譜學與光譜分析,2011,31(8):2161-2165.

[19]鄒華彬. 中藥質量控制多維指紋圖譜共有峰率和變異峰率雙指標序列分析法研究[D]. 濟南:山東中醫藥大學,2004.

[20]王力賓,顧光同. 多元統計分析:模型、案例及SPPS應用[M].北京:經濟科學出版社,2010:135-137.

[21]Falandysz J,Kojta A K,Jarzyńska G,et al. Mercury in bay bolete(Xerocomus badius):bioconcentration by fungus and assessment of element intake by humans eating fruiting bodies[J].Food Additives & Contaminants:Part A,2012,29(6):951-961.

[22]Svoboda L,Chrastny V. Levels of eight trace elements in edible mushrooms from a rural area[J]. Food Additives and Contaminants,2008,25(1):51-58.

[23]Chojnacka A,Drewnowska M,Jarzyńska G,et al. Mercury in Yellow-cracking Boletes Xerocomus subtomentosus mushrooms and soils from spatially diverse sites:Assessment of bioconcentration potential by species and human intake[J]. Journal of Environmental Science and Health,Part A,2012,47(13):2094-2100.

Ultraviolet spectrum identification ofBoletustomentipesfrom different regions

YANG Tian-wei1,2,LI Tao3,ZHANG Ji2,LI Jie-qing1,LIU Hong-gao1,WANG Yuan-zhong2,*

(1. College of Agronomy and Biotechnology,Yunnan Agricultural University,Kunming 650201,China;2. Institute of Medicinal Plants,Yunnan Academy of Agricultural Sciences,Kunming 650200,China 3. College of Resources and Environment,Yuxi Normal University,Yuxi 653100,China)

The ultraviolet fingerprints ofBoletustomentipesfrom nine different areas in Yunnan were established. UV spectral data in the range of 230~450nm were analyzed by Euclidean distance and principal component analysis(PCA). The results showed that,the RSDs of stability within 10h,accuracy and repeatability of the method,were 0.04%~1.73%,0~2.87% and 0.03%~0.63%,respectively. The Euclidean distance value of samples among different origins differed significantly from 0.26 to 6.52. PCA showed that the cumulative contribution rate of the first three factors was 86.848%,which could reflect the most information of the samples,and the two-dimensional projection of the first two principal components could distinguishBoletustomentipessamples from different regions. Ultraviolet spectrum combined with Euclidean distance and principal component analysis(PCA)can rapidly identify different origin ofBoletustomentipes.

ultraviolet spectrum;euclidean distance;principal component analysis(PCA);Boletustomentipes;identification

2014-07-23

楊天偉(1989-),男,碩士研究生,研究方向:食用真菌資源研究。

*通訊作者:王元忠(1981-),男,碩士,助理研究員,研究方向:藥用植物與藥用真菌資源研究。

國家自然科學基金項目(31260496,31160409);云南省自然科學基金項目(2011FB053,2011FZ195);云南省教育廳科學研究基金項目(2013Z074)。

TS201.2

A

:1002-0306(2015)09-0301-06

10.13386/j.issn1002-0306.2015.09.057

主站蜘蛛池模板: 精品国产香蕉伊思人在线| 午夜福利网址| 亚洲精品福利视频| 亚洲欧洲天堂色AV| 色AV色 综合网站| 久久国产精品麻豆系列| 人妻丰满熟妇av五码区| 中文国产成人精品久久| 欧美成人手机在线观看网址| 色网在线视频| 国产真实自在自线免费精品| 中文无码精品a∨在线观看| 欧美成人综合在线| 国产免费自拍视频| 日本成人不卡视频| 成人在线亚洲| 一本大道无码日韩精品影视| 天堂av高清一区二区三区| 国产精品视频导航| 色偷偷综合网| 久爱午夜精品免费视频| 国产精品99久久久| 欧美伊人色综合久久天天| 国产白丝av| 精品国产www| 国产成人亚洲综合a∨婷婷| 色亚洲成人| 国产一区二区三区在线观看视频| 伊人久久福利中文字幕| 激情五月婷婷综合网| 婷婷亚洲综合五月天在线| 亚洲一区色| 一级毛片免费观看久| 国产精品999在线| 国产一区二区三区在线精品专区| 精品無碼一區在線觀看 | 国产又色又刺激高潮免费看| 久久性视频| 国产精品久线在线观看| 国产一区二区三区精品久久呦| 亚洲成人一区二区| 免费无遮挡AV| 中美日韩在线网免费毛片视频| 亚洲最大福利视频网| 欧美一级高清免费a| 亚洲精品少妇熟女| 国产成+人+综合+亚洲欧美| 青青久在线视频免费观看| 国产成人综合亚洲欧美在| a级毛片免费网站| 一区二区偷拍美女撒尿视频| 波多野结衣视频网站| 精品免费在线视频| 中文字幕2区| 97在线免费| 亚洲av日韩综合一区尤物| 久久中文字幕2021精品| a免费毛片在线播放| 久久6免费视频| 免费jjzz在在线播放国产| 91人人妻人人做人人爽男同| 中文字幕亚洲精品2页| www.狠狠| 伊人激情综合网| 91免费片| 国产视频久久久久| 日韩高清一区 | 人人91人人澡人人妻人人爽| 国产欧美日本在线观看| 欧美午夜在线观看| 永久免费AⅤ无码网站在线观看| 亚洲综合二区| 制服丝袜一区| 日韩在线影院| 亚洲色图欧美激情| 久久福利网| 日韩欧美国产精品| 免费毛片全部不收费的| 午夜少妇精品视频小电影| 国产免费网址| 精品国产中文一级毛片在线看| 日本五区在线不卡精品|