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

響應面法優化植物乳桿菌代謝產細菌素的發酵條件

2011-03-30 10:00:30孟祥晨
食品科學 2011年3期
關鍵詞:植物優化模型

陳 琳,孟祥晨

(東北農業大學 乳品科學教育部重點實驗室,黑龍江 哈爾濱 150030)

響應面法優化植物乳桿菌代謝產細菌素的發酵條件

陳 琳,孟祥晨*

(東北農業大學 乳品科學教育部重點實驗室,黑龍江 哈爾濱 150030)

為提高一株分離自內蒙古傳統發酵稀奶油“焦克”中的植物乳桿菌KLDS1.0391代謝產細菌素量,以中性蛋白酶水解脫脂乳為培養基,以枯草芽孢桿菌為指示菌,以抑菌圈直徑為考察指標,在單因素試驗的基礎上,采用響應面法優化發酵pH值、溫度以及接種量。結果表明:對該菌代謝產細菌素的活性影響大小依次為:發酵pH值>接種量>發酵溫度;最優發酵條件為:pH5.1、發酵溫度33℃、接種量1%。在此條件下,發酵液的抑菌圈直徑為15.00mm,細菌素的效價為601.32IU/mL,較優化前提高了43.08%。在最優發酵條件下獲得的實驗結果與模型預測值吻合,說明所建立的模型是切實可行的。

響應面;細菌素;植物乳桿菌;焦克

“焦克”是內蒙古的一種傳統發酵稀奶油產品,屬于自然發酵乳制品,保質期較長,其中蘊含豐富的乳酸菌資源。乳酸菌發酵乳糖產酸形成低pH值是其能夠較長時間保藏的主要原因,其次,乳酸菌代謝產生的具有抑菌活性的物質,如細菌素,也有利于這種產品的保藏[1]。在前期研究中筆者分析了從“焦克”中分離的67株乳酸菌產抑菌物質的能力,發現4株乳桿菌具有較強的產抑菌物質的能力,經排除干擾因素后確定主要的抑菌物質為細菌素。其中,植物乳桿菌KLDS1.0391代謝產生的植物乳桿菌素MG不僅可以抑制單細胞增生利斯特氏菌等革蘭氏陽性細菌,還可抑制沙門氏菌等革蘭氏陰性細菌,該細菌素是具有良好熱、酸穩定性的蛋白質性物質[2],在食品中有良好的應用前景。

乳酸菌代謝產生細菌素的能力和產量,除了受菌株自身特性的影響,還受發酵條件的影響,如培養基、發酵pH值[3]、發酵溫度[4]以及接種量[5]等。本研究以分離自“焦克”的植物乳桿菌KLDS1.0391為研究對象,以前期優化的酶水解脫脂乳為培養基,通過單因素試驗和響應面法優化發酵條件,獲得該菌代謝產生細菌素的最優發酵參數,為該細菌素的開發利用提供參考。

1 材料與方法

1.1 菌株、培養基與試劑

植物乳桿菌(Lactobacillus plantarum)KLDS1.0391分離自內蒙古自治區傳統發酵稀奶油,由乳品科學教育部重點實驗室工業微生物菌種保藏中心凍干保藏。指示菌株:枯草芽孢桿菌(Bacillus subtilis)ATCC6633購自中國藥品生物制品檢定所。

MRS培養基 英國Oxoid公司;營養肉湯培養基北京奧博星生物技術有限責任公司。

脫脂乳(蛋白質含量為35%) 新西蘭Karicare公司;中性蛋白酶(酶活力為0.8AU/g) 北京諾維信公司;乳鏈球菌素Nisin(含Nisin 2.5%) 美國Sigma公司;其他試劑均為分析純。

1.2 儀器與設備

KLF2000 3.7L發酵罐 瑞士比歐生物工程公司;SPX-150B生化培養箱 上海智城分析儀器制造有限公司;ES-2030冷凍干燥機 日本日立公司;HVE-50型高壓滅菌器 日本Hirayama公司;BCN1360型生物潔凈工作臺 上海佳勝實驗設備有限公司。

1.3 方法

1.3.1 菌種的活化與培養

植物乳桿菌KLDS1.0391經活化傳代恢復活力后,以1%接種量接入MRS培養基中,37℃靜置培養24h。

枯草芽孢桿菌ATCC6633經活化傳代恢復活力后,以0.3%接種量接種到營養肉湯培養基中,37℃搖床培養18h。

1.3.2 發酵液抑菌活性的測定

無細胞發酵上清液的制備:取發酵液于4℃、12000×g離心15min,收集上清液,用6mol/L NaOH調pH6.5,以中和有機酸的干擾,經0.22μm濾膜過濾,除去菌體及其他雜質,然后冷凍干燥,凍干后取1g凍干樣品加入1mL滅菌雙蒸水,溶解后放入-20℃冰箱中備用。

采用雙層平板打孔法測定抑菌活性[6]。配制質量濃度1.2g/100mL瓊脂,滅菌后按每平皿10mL傾倒平板,備用。制備含質量濃度0.7g/100mL瓊脂的100mL肉湯培養基,滅菌后冷卻至 50℃左右,接入0.6mL枯草芽孢桿菌培養液,混勻后立即取6mL倒入含有瓊脂的平皿上,冷卻后,用打孔器打孔,孔直徑為6 m m。在孔中加入50μL無細胞發酵上清液,在超凈工作臺上放置3h后,于30℃靜置培養24h,取出平皿后用游標卡尺測量抑菌圈直徑(不含孔徑)。

細菌素效價的測定:根據Cabo等[7]的方法繪制標準效價曲線,以枯草芽孢桿菌ATCC6633為指示菌,以Nisin為陽性對照。將上述獲得的抑菌圈直徑代入效價回歸方程,計算發酵液中細菌素的效價。本實驗測得回歸方程為y=0.2425x+0.6111,y表示效價的對數值;x表示抑菌圈直徑(不含孔徑,單位為mm),R2=0.9868。

1.3.3 酶水解脫脂乳培養基的制備

配制質量濃度4g/100mL的脫脂乳溶液2L,調pH7.0,水浴中升溫至50℃,添加0.05g/100mL的中性蛋白酶,恒溫水解,水解過程中加入2mol/L的氫氧化鈉維持pH 7.0。根據pH-stat方法[8]制備成水解度為15%的脫脂乳溶液,然后添加質量濃度2g/100mL葡萄糖和0.1g/100mL吐溫-80制成用于發酵的酶水解脫脂乳培養基。

1.3.4 單因素試驗

按照1.3.3節方法制備酶水解脫脂乳培養基,裝入3.7L發酵罐中,裝罐量為2L,115℃滅菌15min。選擇發酵pH值、發酵溫度、接種量作為影響抑菌活性的主要因素,通過單因素試驗選取響應面試驗的因素和水平。

1.3.4.1 發酵pH值對植物乳桿菌代謝產細菌素的影響

以1%接種量接種植物乳桿菌KLDS1.0391培養液,發酵pH值分別設為4.5、5.0、5.5、6.0、6.5、7.0,發酵溫度30℃,發酵時間24h。發酵結束后,放罐、取樣,按照1.3.2節所述方法測定發酵液的抑菌圈直徑。

1.3.4.2 接種量對植物乳桿菌代謝產細菌素的影響

分別以1%、2%、3%、4%的接種量接種植物乳桿菌KLDS1.0391培養液,發酵pH值為5.0,發酵溫度30℃,發酵時間24h。發酵結束后,放罐、取樣,按照1.3.2節所述方法測定發酵液的抑菌圈直徑。

1.3.4.3 發酵溫度對植物乳桿菌代謝產細菌素的影響

以1%接種量接種植物乳桿菌KLDS1.0391培養液,發酵pH值為5.0,發酵溫度分別設為30、35、40℃,發酵時間24h。發酵結束后,放罐、取樣,按照1.3.2節所述方法測定發酵液的抑菌圈直徑。

1.3.5 響應面法優化植物乳桿菌代謝產細菌素的發酵條件

根據單因素試驗結果,選取發酵pH值(X1)、發酵溫度(X2)和接種量(X3)3個因素為自變量,以抑菌圈直徑為響應值,進行中心組合試驗設計[9]。每個試驗重復3次,取其平均值,試驗因素水平及編碼見表1。

表1 試驗因素水平及編碼Table 1 Factors, levels and codes of response surface tests

1.3.6 模型的驗證

通過響應面法優化植物乳桿菌代謝產細菌素的發酵條件。在優化條件下進行發酵實驗,通過比較預測值和實驗值驗證模型的有效性。

1.3.7 數據處理與分析

每個試驗重復3次,取其平均值。采用SAS 8.1及Microsoft Excel 2003統計軟件進行單因素方差分析和差異顯著性分析,采用Design Expert 7.1 軟件進行響應面分析。

2 結果與分析

2.1 單因素試驗

2.1.1 發酵pH值對植物乳桿菌KLDS1.0391代謝產細菌素的影響

圖1 發酵pH值對植物乳桿菌KLDS1.0391代謝產細菌素的影響Fig.1 Effect of fermentation pH on bacteriocin produced by Lactobacillus plantarum KLDS1.0391

如圖1所示,發酵pH值對植物乳桿菌發酵產細菌素的影響顯著(P<0.05)。當發酵pH值在4.5~5.0范圍內,菌株代謝產生細菌素的抑菌活性迅速增加,當發酵pH值為5.0時抑菌圈直徑達到最大,隨發酵pH值繼續增加,抑菌活性呈逐漸下降趨勢。根據抑菌圈直徑,選擇發酵pH5.0為響應面試驗中心試驗點。

2.1.2 接種量對植物乳桿菌KLDS1.0391發酵產細菌素的影響

圖2 接種量對植物乳桿菌KLDS1.0391代謝產細菌素的影響Fig.2 Effect of inoculation amount on bacteriocin produced by Lactobacillus plantarum KLDS1.0391

接種量的大小影響乳桿菌的生長周期,接種量對該菌代謝產細菌素的影響如圖2所示。抑菌活性隨接種量增加而逐漸下降,接種量對該菌產細菌素的影響顯著(P<0.05)。根據抑菌圈直徑,選擇接種量為1%進行響應面試驗。

2.1.3 發酵溫度對植物乳桿菌KLDS1.0391代謝產細菌素的影響

圖3 發酵溫度對植物乳桿菌KLDS1.0391代謝產細菌素的影響Fig.3 Effect of fermentation temperature on bacteriocin produced by Lactobacillus plantarum KLDS1.0391

發酵溫度影響菌體的生長和細菌素的合成。如圖3所示,發酵溫度為30℃時,抑菌活性最高,抑菌活性隨發酵溫度的升高呈逐漸下降趨勢,但差異不顯著(P>0.05)。根據抑菌圈直徑,選擇發酵溫度為30℃進行響應面試驗。

2.2 響應面模型的建立及其顯著性

響應面試驗設計及結果見表2。利用Design Expert7.1軟件對表2數據進行多元回歸擬合,獲得的響應值抑菌圈直徑(Y)對編碼自變量x1、x2和x3的二次多項式回歸模

對上述回歸模型進行方差分析,結果見表3。結果表明:模型是顯著的(P<0.0001),失擬項不顯著,回歸模型的決定系數為0.9905,說明該模型能夠解釋99.05%的變化,僅有總變異的0.95%不能用此模型來解釋,該模型擬合程度良好,試驗誤差小,因此,可用此模型對植物乳桿菌KLDS1.0391代謝產細菌素的發酵條件進行分析和預測。

表2 響應面試驗設計及結果Table 2 Design and results of response surface tests

表3 回歸模型方差分析表Table 3 Variance analysis for regression equation

表4 回歸方程系數顯著性檢驗表Table 4 Significance test for regression coefficients

回歸模型的顯著性分析見表4。由回歸方程系數顯著性檢驗可知:模型一次項x1、x3極顯著,x2不顯著;二次項x12極顯著,x22、x32均處于顯著水平;交互項 x2x3顯著,x1x2、x1x3均不顯著。對回歸方程進行中心標準化處理,通過回歸系數絕對值大小分析各個因素的改變對該植物乳桿菌KLDS1.0391代謝產細菌素的影響程度。回歸方程一次項的回歸系數絕對值依次為x1、x3、x2,結果表明:3個因素對該菌代謝產細菌素的影響大小依次為:發酵pH值>接種量>發酵溫度。

利用Design Expert軟件對表2數據進行二次多元回歸擬合,所得到的二次回歸方程的響應面分別見圖4~6。

圖4 發酵溫度和pH值對植物乳桿菌KLDS1.0391所產細菌素抑菌活性影響的響應面圖Fig.4 Response surface and contour plot for the effect of cross-interaction between fermentation temperature and pH on antibacterial activity

由圖4可見,兩個因素的交互項不顯著,由響應面可以看出,發酵溫度不變,隨發酵pH值的增加,抑菌圈直徑先增大,當發酵pH值超過5.1后,抑菌活性呈逐漸下降的趨勢。發酵溫度在30~35℃,發酵pH值在4.8~5.5范圍內,抑菌活性較高。乳酸菌在一定pH值條件下合成細菌素,菌體生長的最適pH值與細菌素合成的最適pH值不一定相同。Izquierdo等[10]在對屎腸球菌WHE 81產細菌素的研究中發現,菌體最適生長pH值為7.0,而細菌素生成的最適pH值為5.0。這與Nel等[11]的研究結果相符,細菌素一般在弱酸性條件下生成,發酵pH值低于菌體生長的最適pH值,細菌素產生的最適pH值具有菌株特異性。

圖5 接種量和發酵pH值對植物乳桿菌KLDS1.0391所產細菌素抑菌活性影響的響應面Fig.5 Response surface and contour plot for the effect of crossinteraction between inoculation amount and pH on antibacterial activity

由圖5可見,兩個因素的交互項不顯著,由響應面可以看出,發酵pH值不變,抑菌活性隨接種量的增加而迅速增加,當接種量超過2%時,抑菌活性迅速下降。接種量在1%~2%,發酵pH值在4.7~5.5范圍內,抑菌活性較高。細菌素的產量多在細菌對數生長后期達到最大值,接種量過大,縮短合成時間,產量有很大程度的下降,細菌素的活性也會降低,這可能也與細菌素吸附到產生菌菌體表面有關[12]。

圖6 接種量和發酵溫度對植物乳桿菌KLDS1.0391所產細菌素抑菌活性影響的響應面Fig.6 Response surface and contour plot for the effect of crossinteraction between inoculation amount and fermentation temperature on antibacterial activity

由圖6可見,兩個因素的交互項顯著,由響應面可以看出,接種量不變,抑菌活性隨發酵溫度的升高而逐漸增加,當發酵溫度超過35℃時,抑菌活性迅速下降。發酵溫度在30~37℃,接種量在1%~2%范圍內,抑菌活性較高,抑菌圈直徑可達到14.95mm。菌體生長的最適溫度不一定適合細菌素的生成。相關研究表明,細菌素在低于菌體最適生長溫度的條件下產量較高,這與細菌素的合成機制有關,而且不同菌株產細菌素的最適發酵溫度也不同[13],擴展短桿菌ATCC9175在28℃時細菌素的產量最高,溫度過高產量下降[14],腸球菌RZSC5產細菌素的溫度必須控制在35℃,溫度過低或過高細菌素產量明顯減少[15]。

2.3 回歸模型的驗證

通過上述回歸模型,采用 Design-Expert 7.1軟件優化發酵條件,獲得植物乳桿菌KLDS1.0391代謝產細菌素的最適發酵條件為:發酵pH5.1、發酵溫度33℃、接種量1%,在此條件下抑菌圈直徑可達到15.08mm。采用上述條件進行驗證實驗,結果獲得的抑菌圈直徑為15.00mm,實驗值與預測值的相對誤差為0.5%,證明采用響應面法優化得到的發酵條件準確可靠,具有實用價值。未進行發酵條件優化時,細菌素的效價為420.26IU/mL,優化后細菌素的效價達到601.32 IU/mL,較優化前提高了43.08%。

3 結 論

3.1 在單因素試驗的基礎上,利用響應面法對植物乳桿菌KLDS1.0391代謝產細菌素的發酵條件進行了優化,建立了抑菌圈直徑與發酵pH值、發酵溫度、接種量3個因素的二次多項式回歸模型,經驗證實驗證明該模型合理可靠。

3.23 個因素對該菌代謝產細菌素的影響大小依次為:發酵pH值>接種量>發酵溫度;最優發酵條件為:pH 5.1、發酵溫度33℃、接種量1%。在此條件下,發酵液的抑菌圈為15.00mm,細菌素的效價為601.32IU/mL,較優化前提高了43.08%。

[1]REVIRIEGO C, FERNANDEZ L, RODRIGUEZ J M. A food-grade system for production of pediocin PA-1 in nisin-producing and nonnisin-producing Lactococcus lactis strains: Application to inhibit Listeria growth in a cheese model system[J]. Food Protection, 2007, 70 (11): 2512-2517.

[2]GONG Hansheng, MENG Xiangchen, WANG Hui. Plantaricin MG active against Gram-negative,,,, by Lactobacillus plantarum KLDS1.0391 isolated from Jiaoke , a traditional fermented cream from China[J]. Food Control, 2010, 21(1): 89-96.

[3]呂燕妮, 李平蘭, 周偉. 戊糖乳桿菌31-1菌株產細菌素發酵條件的優化[J]. 微生物學通報, 2005, 32(3): 13-19.

[4]GRATTEPANCHE F, AUDET P, LACROIX C. Milk fermentation by functional mixed culture producing nisin Z and exopolysaccarides in a fresh cheese model[J]. International Dairy Journal, 2007, 17(2): 123-132.

[5]張大為, 張潔, 鄺仁亮. 產細菌素的嗜酸乳桿菌WS發酵動力學模型的建立[J]. 農業工程學報, 2008, 24(3): 295-297.

[6]NIETO-LOZANO J C, REGUERA-USEROS J I, PELAEZ-MARTINEZ M C, et al. Effect of a bacteriocin produced by Pediococcus acidilactici against Listeria monocytogenes and Clostridium perfringens on Spanish raw meat[J]. Meat Science, 2006, 72(1): 57-61.

[7]CABO M L, MURADO M A, GONZALEZ M P, et al. A method for bacteriocin quantification[J]. Applied Microbiology, 1999, 87(6): 907-914.

[8]CHEISON S C, ZHANG Shaobing, WANG Zhang, et al. Comparison of a modified spectrophotmetric and the pH-stat methods for determination of the degree of hydrolysis of whey proteins hydrolysed in a tangential-fiow filter membrane reactor[J]. Food Research International, 2009, 42(1): 91-97.

[9]董周永, 胡青霞, 郭松年, 等. 石榴果皮中抑菌活性物質提取工藝優化[J]. 農業工程學報, 2008, 24(3): 274-277.

[10]IZQUIERDO E, MARCHIONI E, AOUDE-WERNER D, et al. Smearing of soft cheese with Enterococcus faecium WHE 81, a multi-bacteriocin producer, against Listeria monocytogenes[J]. Food Microbiology, 2009, 26(1): 16-20.

[11]NEL H A, BAUER R, VANDAMME E J, et al. Growth optimization of Pediococcus damnosus NCFB1832 and the influence of pH and nutrients on the production of pediocin PD-1[J]. Applied Microbiology, 2001, 91 (6): 1131-1138.

[12]ASEN I M, MORETRO T, KATLA T, et al. Influence of complex nutrients, temperature and pH on bacteriocin production by Lactobacillus sakei CCUG42687[J]. Applied Microbiology and Biotechnology, 2000, 53(2): 159-166.

[13]ZAMFIR M, CALLEWAERT R, CORNEA P C, et al. Purification and characterization of a bacteriocin produced by Lactobacillus acidophilus IBB801[J]. Applied Microbiology, 2010, 87(6): 923-931.

[14]MOTTA A S, BRANDELLI A. Influence of growth conditions on bacteriocin production on by Brevibacterium linens[J]. Applied Microbiology and Biotechnology, 2003, 62(2/3): 163-167.

[15]LEROY F, VANKRUNKELSVEN S, GREEF J D. The stimulating effect of a harsh environment on the bacteriocin activity by Enterococcus faecium RZSC5 and dependency on the environmental stress factor used [J]. Food Microbiology, 2003, 83(1): 27-38.

Optimization of Fermentation Conditions of Lactobacillus plantarum for Bacteriocin Production by Response Surface Methodology

CHEN Lin,MENG Xiang-chen*
(Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, China)

The fermentation condition of Lactobacillus plantarum KLDS1.0391 isolated from ,, a traditional fermented cream from Inner Mongolia in China, were optimized for bacteriocin production. Skim milk hydrolyzed by neutral protease was used as the medium; Baicillus subtilis was used as indicator bacteria and diameter of inhibition zone was used as evaluation index of antibacterial activity. Based on the single factor tests, the optimal fermentation conditions for producing antibacterial components with high yield in this strain was explored by response surface methodology with three variables of fermentation pH, fermentation temperature and inoculation amount. Results indicated that the order for affecting the yield of antibacterial components from strong to weak was fermentation pH, inoculation amount and fermentation temperature. The optimal fermentation conditions were pH 5.1, inoculation amount of 1% and fermentation temperature of 33 ℃. Under these optimal conditions, the diameter of inhibition zone reached up to 15.00 mm. The antibacterial activity was increased by 43.08% and reached up to 601.32 IU/mL. Therefore, it is feasible for the established model due to the consistent results between the prediction and experiments.

response surface;bacteriocin;Lactobacillus plantarum;Jiaoke

TQ920.6

A

1002-6630(2011)03-0176-05

2010-04-16

黑龍江省留學歸國科學基金項目(LC2009C30);國家“863”計劃項目(2008AA10Z335)

陳琳(1985—),女,碩士研究生,主要從事食品微生物研究。E-mail:chenlinlucky@163.com

*通信作者:孟祥晨(1970—),女,教授,博士,主要從事乳品科學及食品發酵研究。E-mail:xchmeng@163.com

猜你喜歡
植物優化模型
一半模型
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
哦,不怕,不怕
將植物穿身上
3D打印中的模型分割與打包
主站蜘蛛池模板: 特级精品毛片免费观看| 中文毛片无遮挡播放免费| 国产精品人成在线播放| 久久不卡国产精品无码| 另类综合视频| 97国产一区二区精品久久呦| 欧美高清国产| AV网站中文| 欧美亚洲第一页| www.亚洲天堂| 婷婷色狠狠干| 99re在线观看视频| 欧美色伊人| 国产成人91精品免费网址在线| 成年A级毛片| 国产丝袜丝视频在线观看| 日韩黄色精品| 青青国产视频| 国产毛片基地| 午夜视频免费一区二区在线看| 国产精品视频导航| 国产精品三级av及在线观看| 9久久伊人精品综合| 亚洲性影院| aaa国产一级毛片| 毛片网站免费在线观看| 亚洲香蕉久久| 久无码久无码av无码| 日本在线国产| 午夜精品福利影院| 日韩精品一区二区三区中文无码| 蜜臀av性久久久久蜜臀aⅴ麻豆| 成人午夜久久| 国产成人一区二区| 色丁丁毛片在线观看| 九九视频免费看| 日韩欧美一区在线观看| 精品在线免费播放| 免费午夜无码18禁无码影院| 日本人妻丰满熟妇区| 国产欧美精品专区一区二区| 国产精品v欧美| V一区无码内射国产| 国产极品美女在线| 中文字幕无码av专区久久| 波多野结衣在线一区二区| 高清欧美性猛交XXXX黑人猛交| 91色综合综合热五月激情| 国产极品美女在线观看| a欧美在线| 亚洲人成在线免费观看| 久久成人国产精品免费软件| 久久久波多野结衣av一区二区| 超碰精品无码一区二区| 波多野结衣一区二区三区AV| 青青草一区| 激情無極限的亚洲一区免费| 丰满人妻久久中文字幕| 亚卅精品无码久久毛片乌克兰 | 色天天综合久久久久综合片| 精品伊人久久久久7777人| 久久伊人色| 欧美区一区| 日本91在线| 亚洲中文字幕国产av| 国产精品成人一区二区| 日本91在线| 九一九色国产| 亚洲天堂免费在线视频| 国产激情无码一区二区三区免费| 欧美亚洲欧美区| 国产成人精品免费av| 5555国产在线观看| 国产香蕉国产精品偷在线观看| 国产69精品久久久久妇女| 久久国产高潮流白浆免费观看| 深爱婷婷激情网| 久久免费成人| 亚洲第一区在线| 999精品在线视频| 欧美视频在线播放观看免费福利资源 | 国产精品久线在线观看|