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

響應面法優化蠶豆多酚 超聲輔助提取工藝

2018-05-30 19:02:16李春陽劉玉皎
食品工業科技 2018年9期
關鍵詞:優化影響實驗

錢 敏,李春陽,*,劉玉皎,王 帆,陳 新,吳 寒

(1.南京農業大學食品科技學院,江蘇南京 210095;2.江蘇省農業科學院農產品加工研究所,江蘇南京 210014;3.青海省農林科學院,青海西寧 810000;4.江蘇省農業科學院經濟作物研究所,江蘇南京 210014)

蠶豆在我國已有二千年栽培史[1],是我國目前種植面積和產量位居第三的食用豆類作物,僅次于大豆和花生[2],分別占世界蠶豆種植面積和總產量的53.8%和64.1%[3]。中醫認為蠶豆的花、果實、種殼,及葉均可入藥,有止血,利尿解毒,消腫之功效,并且國外已有用蠶豆提取抗癌物質的報道[4]。有研究證明蠶豆中含有豐富的酚類化合物,具有降低小鼠心血管疾病的風險,可作為很好的天然抗氧化食物和功能食物的來源[5]。

目前,國內對豆類酚類物質研究較多,主要是比較不同品種豆類多酚含量及其抗氧化活性[6-8],李思荻等[7]研究發現,合豐大豆和白蕓豆對羥基自由基的清除能力強于其他豆類,提取液中的多酚含量與其抗氧化性成正相關關系。國外對于蠶豆多酚的研究主要集中在其抗氧化活性[9-11]及酚類物質組成成分分析[12-14],Marathe等[9]研究表明,表皮顏色較深的雜豆抗氧化性較強,且抗氧化性強弱與豆中所含酚類物質有關,Journi等[12]從蠶豆多酚中鑒定出槲皮素糖苷、表兒茶素糖苷等五種成分,Merghem等[13]利用電子噴霧質譜及高效液相色譜技術從蠶豆原花青素中分離鑒定出原花青素二聚物B1、B2和B3及原花青素三聚物等六種主要化合物。

由于超聲提取具有操作簡單方便、提取率高、提取溫度低、不破環提取物結構等特點,已廣泛用于提取天然植物有效成分,如橙皮多酚的提取[15]、海棠葉中總黃酮的提取[16]、茶多酚的提取[17]、獼猴桃根熊果酸的提取[18]、大豆異黃酮的提取[19]等。目前有關蠶豆多酚的提取主要還是溶劑浸提法,關于超聲波輔助提取蠶豆多酚工藝優化的報道尚不多見。因此本文將系統優化蠶豆多酚超聲提取工藝,并在已有報道的基礎上考慮了提取液的pH對多酚提取量的影響,以期獲得適宜的提取方法,為將來工業化高效提取蠶豆中酚類物質提供參考,以便進行蠶豆的保健品開發和蠶豆抗氧化劑制品的研究。

1 材料與方法

1.1 材料與儀器

蘇蠶三號 由江蘇省農業科學院蔬菜所提供;甲醇、乙醇、丙酮、無水碳酸鈉 均為國產分析純;沒食子酸標準品與Folin-Ciocalteu試劑 購自Sigma公司。

AL104型電子天平 梅特勒托利多儀器有限公司;KQ-500DE型數控超聲波清洗器 昆山市超聲儀器有限公司;Sigma 3K15離心機 德國Sigma公司;RE-5203旋轉蒸發儀 上海亞榮生化儀器廠;萬能高速粉碎機 歐凱萊芙(香港)寶業公司;PHS-2C型數顯pH計 上海智光儀器儀表公司;752S型紫外可見分光光度計 上海棱光技術有限公司。

1.2 實驗方法

1.2.1 沒食子酸標準曲線的制作 配制濃度為1 mg/mL的沒食子酸標準溶液,分別吸取1.00、2.00、3.00、4.00、5.00、10.00 mL于100 mL容量瓶中定容,得10、20、30、40、50、100 μg/mL的沒食子酸標準溶液,移取1 mL沒食子酸標準溶液于試管中,加入5 mL蒸餾水,1 mL 10%福林酚,反應3 min后加3 mL 7.5% Na2CO3,室溫下放置2 h,于765 nm處測吸光值[20]。以濃度(X)對吸光度(Y)進行線性回歸分析,繪制出沒食子酸標準曲線。

求得回歸方程為Y=0.0088X+0.0558(0~100 μg/mL,R2=0.9996),其中X為沒食子酸濃度,單位為μg/mL,Y為吸光度。

1.2.2 蠶豆多酚的提取及含量測定 將蠶豆冷凍干燥(-20 ℃預凍6 h,-80 ℃冷凍20 h,放入凍干機中,在冷阱溫度為-45 ℃,真空度為13.3 Pa、加熱隔板溫度為30 ℃左右冷凍干燥24 h),用萬能粉碎機將蠶豆粉碎,過60目篩,制得蠶豆粉(水分含量2.98%±0.06%,w/w),-20 ℃保存備用。準確稱取1 g蠶豆粉,按料液比1∶15 (g/mL)加入60%乙醇溶液進行超聲提取,超聲功率300 W,超聲溫度35 ℃,提取20 min后4500 r/min離心10 min,收集上清液,殘渣用同樣的方法再次提取,合并兩次上清液,45 ℃減壓旋轉蒸干,甲醇定容至10 mL,-20 ℃避光儲存,用于多酚含量的測定。多酚提取量以每克提取物(干基)中所含相當于沒食子酸的量mg表示。

式中:C為測量液總酚濃度(μg/mL);V為提取液的總體積(mL);N為稀釋倍數;M為樣品除去水分質量(g);1000為質量轉化單位(mg/μg)。

1.2.3 提取劑種類的選擇 準確稱取1 g蠶豆粉,按料液比1∶15 (g/mL)的比例分別加入體積分數均為60%的甲醇、乙醇、丙酮,用1 mol/L HCl調節至pH4(下文同),在300 W、35 ℃下超聲提取20 min,提取兩次,研究不同提取劑種類對蠶豆多酚提取效果的影響。

1.2.4 單因素實驗

1.2.4.1 乙醇體積分數對提取效果的影響 準確稱取1 g蠶豆粉,按料液比1∶15 (g/mL)分別加入pH4,不同濃度的乙醇水溶液(20%、40%、60%、80%、100%),300 W、35 ℃下超聲提取20 min,提取兩次,研究乙醇體積分數對蠶豆多酚提取效果的影響。

1.2.4.2 料液比對提取效果的影響 準確稱取1 g蠶豆粉,分別加入10、15、20、25、30 mL體積分數為60%、pH4的乙醇水溶液,300 W、35 ℃下超聲提取20 min,提取兩次,研究料液比對蠶豆多酚提取效果的影響。

1.2.4.3 提取時間對提取效果的影響 準確稱取1 g蠶豆粉5份,分別加入15 mL體積分數為60%、pH4的乙醇溶液。300 W、35 ℃下超聲提取10、15、20、25、30 min,提取兩次,研究提取時間對蠶豆多酚提取效果的影響。

1.2.4.4 提取液pH對提取效果的影響 準確稱取1 g蠶豆粉,分別加入15 mL體積分數為60%不同pH3、4、5、6、7、8、9的乙醇溶液,300 W、35 ℃下超聲提取20 min,提取兩次,研究pH對蠶豆多酚提取效果的影響。

1.2.5 響應面實驗設計 根據Box-Benhnken中心組合設計原則,以乙醇體積分數、料液比、提取時間、pH四個因素為自變量,多酚提取量為響應值,設計了四因素三水平的響應面分析實驗,各因素水平見表1。

表1 Box-Benhnken響應面實驗因子與水平Table 1 Variables and levels in response surface design

1.3 數據處理

所有實驗均重復三次,實驗數據以平均值±標準差(Means±SD)表示,采用Excel 2003軟件作圖,使用SAS 9.3 軟件對數據結果進行處理分析,設置顯著水平為p<0.05,極顯著水平為p<0.01。中心組合實驗數據處理使用Design Expert 8.0.6軟件進行。

2 結果與分析

2.1 提取劑種類對提取效果的影響

從圖1可以看出,60%不同提取劑的提取效果依次為:乙醇>甲醇>丙酮,且乙醇的成本較低、環境友好,所以本文選擇乙醇作為提取溶劑。在楊希鵑等[21]的蠶豆多酚提取實驗中,60%甲醇的提取率高于60%乙醇,這可能是因為不同品種蠶豆中的酚類物質組成不同,導致極性不同,對甲醇和乙醇的親和性不同,從而導致最佳提取溶劑種類有差異。

圖1 提取劑種類對多酚提取量的影響Fig.1 The influence of the type of extracting solution on extraction volume

2.2 單因素實驗結果與分析

2.2.1 乙醇體積分數對提取效果的影響 由圖2可以看出,隨著乙醇濃度增大,多酚提取量呈現先增加后降低的趨勢,當乙醇濃度達到40%時,提取量最大。分析原因可能是蠶豆中多酚與多糖、蛋白質等以氫鍵和疏水鍵形成穩定的化合物,隨著乙醇濃度升高,對氫鍵斷裂作用越強,多酚提取量升高[22]。然而,當乙醇濃度超過40%時,多酚提取量下降,根據相似相溶原理推測可能是由于溶劑與多酚極性差異增大,多酚類物質得不到充分溶解,提取量反而下降。本文中多酚提取量隨乙醇體積分數的變化趨勢與王文昕等[23]的響應面法優化花生紅衣多酚提取的報道是相符的。因此本文選擇乙醇體積分數為40%。

圖2 乙醇體積分數對多酚提取量的影響Fig.2 The influence of the concentration of the extracting solution on extraction volume

2.2.2 料液比對提取效果的影響 由圖3可以看出,當料液比小于1∶25 (g/mL)時,隨料液比增大,蠶豆多酚提取量增大,這可能是因為,蠶豆組織中多酚與蛋白質、生物堿、多糖等以氫鍵等形式結合在一起,阻礙多酚類物質的提取,而有機溶劑和水的混合液能有效的抑制這種結合,因此提取液增大有利于蠶豆多酚的提取。但當料液比達到1∶25 (g/mL)以后,再增大溶劑的用量,多酚提取量下降,而且溶劑用量太大,會使熱負荷增大,提取完全所需要的時間增大[24]。在郭彩霞等[25]的響應面優化超聲波輔助提取獼猴桃果皮多酚實驗中也觀察到類似現象。因此本文選擇的料液比為1∶25 (g/mL)。

圖3 料液比對提取量的影響Fig.3 The influence of the solid-liquid ratio on extraction volume

2.2.3 提取時間對提取效果的影響 由圖4可以看出,當提取時間從5 min增加到20 min時,多酚提取量由1.31 mg/g d.w.增加至1.74 mg/g d.w.。這可能是因為超聲波提取能加速介質質點運動,同時產生空化作用,使提取物質加速浸出[26],當超聲時間大于20 min時,多酚提取量呈下降趨勢,這可能與超聲的時間過長,產生熱效應使溫度升高,從而使得多酚類物質結構被破環有關[27]。在Yang等[28]的響應面法優化余甘子樹皮多酚及Sun等[29]的響應面法優化苦丁茶中酚類物質提取的實驗中也觀察到類似現象。因此,選擇超聲提取時間為20 min。

圖4 提取時間對多酚提取量的影響Fig.4 The influence of the extracting time on extraction volume

2.2.4 提取液pH對提取效果的影響 由圖5可以看出,不同的pH對提取效果有顯著的影響,當pH為5時,多酚提取量最大。這可能是由于在酸性條件下多酚與酸形成佯鹽,以離子形式存在,容易析出;而且在酸性條件下多酚氧化酶的活性受到抑制,從而抑制了多酚在提取過程中的酶促氧化,進而提高了多酚的提取率。與趙謀明等[30]報道的多酚提取量隨pH變化趨勢相符。因此本文選擇pH5為適宜pH。

圖5 提取液pH對多酚提取量的影響Fig.5 The influence of the pH on extraction volume

2.3 響應面優化實驗

2.3.1 二次響應面回歸模型的建立與分析 響應面設計與結果見表2。應用Design Expert 進行回歸擬合分析,得到提取條件與多酚提取量之間的二次多項式模型為:

表2 響應面實驗設計與結果Table 2 Response surface experimental design and results

Y=2.17-(9.167E-003)A-0.028B+0.064C-0.026D-0.057AB-0.053AC-0.057AD-0.013BC+(2.500E-003)BD+0.047CD-0.13A2-0.19B2-0.22C2-0.20D2。

表3 回歸方程系數顯著性檢驗表Table 3 Test of significance for regression equation coefficients

2.3.2 兩因子間交互作用分析 響應面分析圖見圖6~圖11。

由圖6可知,當料液比一定時,隨著乙醇體積分數的增大,提取量先增大后減小;當乙醇體積分數一定時,隨料液比的變化,提取量也呈先增大后減小趨勢,且響應面坡度較陡峭,說明乙醇體積分數和料液比的交互作用極顯著。

圖6 乙醇體積分數和料液比及其相互作用對蠶豆多酚提取量影響的響應面和等高線圖Fig.6 Response surface and contour plots for extraction volume under the concentration of the extracting solution and solid-liquid ratio

圖7表示pH為5、料液比為1∶25 (g/mL)時乙醇體積分數和時間對蠶豆多酚提取量的影響。當時間一定時,隨乙醇體積分數變化,多酚提取量變化不明顯,趨勢較平緩;當乙醇體積分數一定時,隨時間變化,提取量先增后減,且弧度較陡,說明時間的主效應大于乙醇體積分數。

圖7 乙醇體積分數和時間及其相互作用對蠶豆多酚提取量影響的響應面和等高線圖Fig.7 Response surface and contour plots for extraction volume under the concentration of the extracting solution and extracting time

圖8表示當料液比為1∶25 (g/mL)、提取時間為20 min時,乙醇體積分數和pH對蠶豆多酚提取量的影響。當pH一定時,隨乙醇體積分數的變化,提取量先增后減,趨勢較平緩;當乙醇體積分數一定時,隨pH變化,提取量先增大后減小,幅度較陡峭,說明pH的主效應大于乙醇體積分數,又等高線圖呈橢圓形,說明兩者間交互作用顯著。

圖8 乙醇體積分數和pH及其相互作用對蠶豆多酚提取量影響的響應面和等高線圖Fig.8 Response surface and contour plots for extraction volume under the concentration of the extracting solution and pH

圖9表示當乙醇體積分數為45%、pH為5時,料液比和時間對蠶豆多酚提取量的影響。當時間一定時,隨料液比的增大,蠶豆提取量先增大后減小;當料液比一定時,隨時間變化,料液比先增后減。由等高線為圓形,可知兩者交互作用不顯著。

圖9 料液比和時間及其相互作用對蠶豆多酚提取量影響的響應面和等高線圖Fig.9 Response surface and contour plots for extraction volume under solid-liquid ratio and extracting time

圖10表示乙醇體積分數為45%、提取時間為20 min時,料液比和pH對蠶豆多酚提取量的影響。當pH一定時,隨料液比升高,提取量先增后減;當料液比一定時,隨pH變化,提取量先升高后降低。由其等高線呈圓形,可知兩者交互作用不顯著。

圖10 料液比和pH及其相互作用 對蠶豆多酚提取量影響的響應面和等高線圖Fig.10 Response surface and contour plots for extraction volume under solid-liquid ratio and pH

圖11表示乙醇體積分數為45%、料液比為1∶25 (g/mL)時,時間和pH對蠶豆多酚提取量的影響。當pH一定時,隨時間升高,提取量先增后減;當時間一定時,隨pH變化,提取量先升高后降低。

圖11 時間和pH及其相互作用 對蠶豆多酚提取量影響的響應面和等高線圖Fig.11 Response surface and contour plots for extraction volume under extracting time and pH

2.3.3 最佳條件的預測及驗證實驗 通過響應面優化得到超聲波提取蠶豆多酚的最佳提取工藝條件為:乙醇體積分數44.71%、料液比1∶23.89 (g/mL)、時間20.09 min、pH4.8,此時蠶豆多酚提取量達到最大值(2.06±0.23) mg/g d.w.。為方便實際操作,各條件調整為:乙醇體積分數45%、料液比1∶25 (g/mL)、時間20 min、pH5,在此條件下進行驗證,測得蠶豆多酚提取量為(2.03±0.16) mg/g d.w.,與預測值(2.06±0.23) mg/g d.w.誤差值為1%,證實了該模型的有效性。本實驗多酚提取量與楊希鵑等[21]在正交優化條件甲醇體積分數60%,料液比1∶15 (g∶mL),超聲波功率300 W 的條件下提取 20 min,得到的蠶豆多酚提取量47.78 mg/100 g d.w.相比較高,推測可能是因為本實驗考慮了提取溶劑的pH,有機酸的添加促進了結合酚的釋放,從而提高了總酚的含量;該值與鞏藹[8]報道的用60%乙醇溶液水浴浸提3 h得到的蠶豆多酚含量值0.25%相近,但與Siah[31]用70%丙酮提取得到的多酚提取量6.2~10.7 mg/g d.w.、趙艷等[6]用70%丙酮(丙酮∶水∶乙酸=70∶29.5∶0.5,V/V/V)在25 ℃下振蕩3 h得到的總酚提取量5.80~8.37 mg/g d.w.實驗值相比還是較低,結果的差異可能與蠶豆品種、生長環境、栽培條件以及實驗方法等的不同都有很大關系。

3 結論

本實驗采用單因素實驗及響應面Box-Benhnken實驗設計,對超聲輔助提取蠶豆多酚工藝進行優化,得到的最佳工藝條件為:乙醇體積分數45%、料液比1∶25 (g/mL)、提取時間20 min、pH5,在此條件下蠶豆多酚提取量為(2.03±0.16) mg/g d.w.,與預測值(2.06±0.23) mg/g d.w.高度相符,模型可靠。

[1]游修齡. 蠶豆的起源和傳播問題[J]. 自然科學史研究,1993,12(2):166-167.

[2]Xiaohong T H T B T,Ming L. Chemicophysical properties of twenty broad bean varieties and their correlations[J]. Journal of the Chinese Cereals and Oils Association,2009,24(12):153-157.

[3]Tian X H,Tan H Z,Tan B,et al. Analyse physicochemical properties of broad beans in main producing areas in China[J]. Science and Technology of Cereals Oils and Foods,2009,17(2):7-12.

[4]李雪琴,裘愛泳. 蠶豆生理活性物質研究進展[J]. 糧食與油脂,2002(7):34-35.

[5]Macarulla M T,Medina C,De Diego M A,et al. Effects of the whole seed and a protein isolate of faba bean(Viciafaba)on the cholesterol metabolism of hypercholesterolaemic rats[J]. British Journal of Nutrition,2001,85(5):607-614.

[6]趙艷,杜雙奎,王晗欣,等. 5 種雜豆體外抗氧化性研究[J].中國糧油學報,2015,30(10):6-10.

[7]李思荻,趙永煥. 東北幾種豆類提取物的體外抗氧化活性研究[J]. 糧油加工,2010(4):54-57.

[8]鞏藹,任順成,常云彩,等. 云南特色豆類的酚類含量及其抗氧化活性[J]. 中國糧油學報,2015,30(4):1-5.

[9]Marathe S A,Rajalakshmi V,Jamdar S N,et al. Comparative study on antioxidant activity of different varieties of commonly consumed legumes in India[J]. Food and chemical toxicology,2011,49(9):2005-2012.

[10]Turco I,Bacchetti T,Bender C,et al. Polyphenol content and glycemic load of pasta enriched with Faba bean flour[J]. Functional Foods in Health and Disease,2016,6(5):291-305.

[11]Okumura K,Hosoya T,Kawarazaki K,et al. Antioxidant activity of phenolic compounds from fava bean sprouts[J]. Journal of Food Science,2016,81(6):C1394-C1398.

[12]Journi M,Hammouda H,Trabelsi-Ayadi M,et al. Quantitative determination and identification of phenolic compounds of three Tunisian legumes:Viciafaba,LensculinarisandPhaseolusvulgaris[J]. Chemistry and Biochemistry Sciences,2015,2(3):1-12.

[13]Merghem R,Jay M,Brun N,et al. Qualitative analysis and HPLC isolation and identification of procyanidins fromViciafaba[J]. Phytochemical Analysis,2004,15(2):95-99.

[14]Baginsky C,Pea-Neira á,Cáceres A,et al. Phenolic compound composition in immature seeds of fava bean(ViciafabaL.)varieties cultivated in Chile[J]. Journal of food composition and analysis,2013,31(1):1-6.

[15]Khan M K,Abert-Vian M,Fabiano-Tixier A S,et al.

Ultrasound-assisted extraction of polyphenols(flavanone glycosides)from orange(CitrussinensisL.)peel[J]. Food Chemistry,2010,119(2):851-858.

[16]喬孟,屈曉清,丁之恩. 響應面法優化超聲波輔助提取湖北海棠葉中總黃酮工藝[J]. 食品科學,2013,34(2):143-147.

[17]Guo S Q,Sheng-Ju W U,Dai L I. Optimization of ultrasonic extraction of tea polyphenols from green tea by response surface methodology[J]. Chinese Journal of Bioprocess Engineering,2009,7(1):39-43.

[18]楊艷,任亞梅,馬婷,等. 響應面優化超聲波提取獼猴桃根熊果酸工藝[J]. 食品科學,2014,35(4):44-49.

[19]Rostagno M A,Palma M,Barroso C G. Ultrasound-assisted extraction of soy isoflavones[J]. Journal of Chromatography A,2003,1012(2):119-128.

[20]張清安,范學輝,張志琪,等. 沙苑子提取物中多酚含量的測定[J]. 食品科學,2010,31(14):178-181.

[21]楊希娟. 青海蠶豆多酚的提取工藝優化及不同品種(系)酚類含量及抗氧化活性[J]. 食品與發酵工業,2016,42(11):237-243.

[22]孫曉雪. 甘蔗梢中酚類物質的提取與應用研究[D]. 南寧:廣西大學,2007.

[23]王文昕,董全. 響應面法優化超聲輔助提取花生紅衣多酚工藝[J]. 食品科學,2012,33(22):1-5.

[24]郭新竹,寧正祥. 天然酚類化合物及其保健作用[J]. 食品工業,2002,23(3):28-29.

[25]郭彩霞,任曉婷,張生萬,等. 響應面法優化超聲波輔助提取獼猴桃果皮多酚工藝研究[J]. 食品工業科技,2017(11):38.

[26]吳棉國,林彥銓. 大力推進甘蔗標準化體系建設的對策及建議[J]. 中國農學通報,2011,27(1):456-460.

[27]鄭利琴. 楊梅多酚提取及應用研究[D]. 無錫:江南大學,2011.

[28]Yang L,Jiang J G,Li W F,et al. Optimum extraction process of polyphenols from the bark ofPhyllanthusemblicaL. based on the response surface methodology[J]. Journal of Separation Science,2009,32(9):1437.

[29]Sun Y,Xu W,Zhang W,et al. Optimizing the extraction of phenolic antioxidants from kudingcha made from Ilex kudingcha C.J. Tseng by using response surface methodology[J]. Separation & Purification Technology,2011,78(3):311-320.

[30]趙謀明,劉曉麗,崔春,等. 余甘子多酚響應面法優化提取及其抗氧化活性研究[J]. 食品工業科技,2007,28(6):117-120.

[31]Siah S,Konczak I,Wood J A,et al. Effects of roasting on phenolic composition andinvitroantioxidant capacity of Australian grown faba beans(ViciafabaL.)[J]. Plant Foods for Human Nutrition,2014,69(1):85-91.

歡迎訂閱《食品工業科技》

猜你喜歡
優化影響實驗
記一次有趣的實驗
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
是什么影響了滑動摩擦力的大小
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
哪些顧慮影響擔當?
當代陜西(2021年2期)2021-03-29 07:41:24
做個怪怪長實驗
NO與NO2相互轉化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
主站蜘蛛池模板: 91美女视频在线| 中国一级特黄视频| 亚洲乱码精品久久久久..| 国产精品福利导航| 青青青国产视频| 亚洲AV无码久久精品色欲| 54pao国产成人免费视频| 久久久久国产一区二区| 一级毛片在线播放| 91视频99| 在线播放精品一区二区啪视频| 狠狠做深爱婷婷综合一区| 成人在线欧美| 国产黄色视频综合| 99在线国产| 日本午夜精品一本在线观看| 色135综合网| 99久久精品免费看国产免费软件| 亚洲无码高清一区| 国产91丝袜在线播放动漫| 国产成人精品高清在线| 日韩免费成人| 亚洲av无码专区久久蜜芽| 欧美成人精品高清在线下载| 91精品国产自产在线观看| 国产香蕉在线视频| 野花国产精品入口| 国产免费网址| 欧美在线黄| 欧美三级视频在线播放| 欧美a级在线| 日韩a在线观看免费观看| 色一情一乱一伦一区二区三区小说| 欧美特级AAAAAA视频免费观看| 国产欧美日韩在线一区| www.亚洲天堂| 久久国产毛片| 欧美日韩精品一区二区视频| 超碰色了色| 亚洲精品欧美日韩在线| 欧美不卡视频一区发布| 欧美啪啪网| 久久久久久久97| 国产乱子伦一区二区=| 国产JIZzJIzz视频全部免费| 91啦中文字幕| 影音先锋亚洲无码| 午夜一级做a爰片久久毛片| 成年女人a毛片免费视频| 国产白丝av| 欧美亚洲综合免费精品高清在线观看| 手机永久AV在线播放| 久久综合亚洲色一区二区三区| 五月婷婷激情四射| 亚洲a免费| 国产欧美亚洲精品第3页在线| 9cao视频精品| 婷婷在线网站| 亚洲成a人片77777在线播放| 亚洲一级色| 亚洲无线一二三四区男男| 欧美视频在线播放观看免费福利资源| 久久午夜夜伦鲁鲁片不卡| 国产精品久久久久久久久久久久| 天堂岛国av无码免费无禁网站 | 国产99视频在线| 亚洲日韩久久综合中文字幕| 91在线一9|永久视频在线| 日韩精品久久无码中文字幕色欲| 国产网友愉拍精品视频| 国产乱人乱偷精品视频a人人澡| 伊人精品成人久久综合| 欧美 亚洲 日韩 国产| 国内精品九九久久久精品| 国产第一福利影院| 狠狠ⅴ日韩v欧美v天堂| 亚洲日本韩在线观看| 国产精品网拍在线| 九九视频在线免费观看| 国产91熟女高潮一区二区| 中文国产成人久久精品小说| 九九视频在线免费观看|