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

    基于響應面法對皇冠牡丹(Paeonia×lemoinei ‘Yellow Crown’)增殖培養(yǎng)基大量元素配方的優(yōu)化

    2024-09-24 00:00:00李勝皓張曉芝潘月鄭冰雁文書生
    江蘇農(nóng)業(yè)學報 2024年8期

    收稿日期:2023-08-09

    基金項目:國家自然科學基金項目(32001359);南京林業(yè)大學大學生創(chuàng)新訓練計劃項目(2022NFUSPITP0436)

    作者簡介:李勝皓(1999-),女,重慶人,碩士研究生,研究方向為園林植物與觀賞園藝。(E-mail)lynel1104@163.com

    通訊作者:文書生,(E-mail)shusheng0507@126.com

    摘要: 增殖系數(shù)低一直是制約牡丹微繁殖技術(shù)推廣應用的障礙之一,但目前的研究多集中于植物生長調(diào)節(jié)劑的使用與組合,較少關(guān)注培養(yǎng)基配方對牡丹試管苗增殖效果的影響。本研究以皇冠牡丹(Paeonia×lemoinei ‘Yellow Crown’)為試驗材料,使用Design-Expert響應面法設(shè)計優(yōu)化試驗,探究WPM中5種大量元素[K2SO4、MgSO4·7H2O、KH2PO4、NH4NO3、Ca(NO3)2·4H2O]對皇冠牡丹試管苗增殖效果的影響,并建立數(shù)學模型分析預測最佳優(yōu)化培養(yǎng)基配方。單因素效應分析結(jié)果表明,KH2PO4、NH4NO3與Ca(NO3)2·4H2O對牡丹試管苗的增殖效果有顯著或極顯著影響,而MgSO4·7H2O與K2SO4對牡丹試管苗的增殖效果無明顯影響。交互作用分析結(jié)果表明,K2SO4與KH2PO4的交互作用對株高、莖長有顯著或極顯著影響,K2SO4與MgSO4·7H2O、KH2PO4與NH4NO3的交互作用對莖長有顯著或極顯著影響。根據(jù)試驗數(shù)據(jù)構(gòu)建模型進行分析預測,結(jié)果顯示,5種大量元素優(yōu)化配方為:K2SO4 2 257 mg/L、MgSO4·7H2O 907 mg/L、KH2PO4 481 mg/L、NH4NO3 200 mg/L與Ca(NO3)2·4H2O 1 963 mg/L。在該條件下開展驗證試驗,獲得試管苗的增殖系數(shù)、株高、莖長與葉片數(shù)分別為2.43、4.62 cm、1.30 cm、6.75張,略高于Design-Expert試驗得出的預測值,顯著優(yōu)于對照組,說明本研究成功優(yōu)化了皇冠牡丹的專用增殖培養(yǎng)基配方。研究結(jié)果不僅能大幅提高皇冠牡丹的增殖效率,也為其他牡丹品種增殖培養(yǎng)體系的構(gòu)建與改良提供重要的技術(shù)參考。

    關(guān)鍵詞: 牡丹;微繁殖;響應面法;增殖;WPM培養(yǎng)基

    中圖分類號: S685.11 文獻標識碼: A 文章編號: 1000-4440(2024)08-1493-14

    Optimization of macro-elements composition in the multiplication medium of Paeonia×lemoinei ‘Yellow Crown’ based on response surface methodology

    LI Shenghao, ZHANG Xiaozhi, PAN Yue, ZHENG Bingyan, WEN Shusheng

    (College of Landscape Architecture, Nanjing Forestry University, Nanjing 210000, China)

    Abstract: Low multiplication coefficient has been one of the obstacles limiting the industrial application of tree peony (Paeonia sect. Moutan) micropropagation technology, but most of the current studies have been focused on the use and combination of plant growth regulators, and less attention has been paid to the effect of medium formulation on the multiplication of test-tube seedlings of tree peony. In this study, we used Paeonia×lemoinei ‘Yellow CrowECVEQI0XsugZ00QQ/y8GO951BQ0uRBlPk+pQUPBaqNE=n’ as the experimental material and designed an optimization test by using Design-Expert response surface method to investigate the effects of five macro-elements, such as K2SO4, MgSO4·7H2O, KH2PO4, NH4NO3, Ca(NO3)2·4H2O in WPM medium on the multiplication effect of the test-tube seedlings of Paeonia sect. Moutan, and developed a mathematical model to predict the best optimized medium formulation. Results of single factor effect analysis showed that, KH2PO4, NH4NO3 and Ca(NO3)2·4H2O had significant or highly significant effects on the multiplication effect of the test-tube seedlings, while MgSO4·7H2O and K2SO4 had no significant effects. Results of interaction analysis showed that, the interaction between K2SO4 and KH2PO4 had significant or highly significant effects on plant height, and stem length, the interactions between K2SO4 and MgSO4·7H2O, KH2PO4 and NH4NO3 had significant or hightly significant effects on stem length. Based on the experimental data, the model analysis results predicted that the optimized formulations for five macro-elements were K2SO4 2 257 mg/L, MgSO4·7H2O 907 mg/L, KH2PO4 481 mg/L, NH4NO3 200 mg/L and Ca(NO3)2·4H2O 1 963 mg/L. The validation test was conducted, and the multiplication coefficient, plant height, stem length and leaf number of the test-tube seedlings were 2.43, 4.62 cm, 1.30 cm and 6.75 respectively, which were slightly higher than the predicted values obtained by Design-Expert experiment and were significantly better than the control group, indicating the special multiplication medium for P.×lemoinei ‘Yellow Crown’ was constructed successfully. The results of this study can not only greatly improve the multiplication efficiency of P.×lemoinei ‘Yellow Crown’, but can also provide an important technical reference for the construction and improvement of multiplication culture system of other tree peony species.

    Key words: tree peony;micropropagation;response surface method;multiplication;WPM medium

    牡丹(Paeonia sect. Moutan)為芍藥科芍藥屬落葉灌木,是原產(chǎn)于中國的傳統(tǒng)名花和藥用植物,近年來被發(fā)現(xiàn)作為新型油料植物也有巨大的市場潛力[1]?;使谀档ぃ≒.×lemoinei ‘Yellow Crown’)為20世紀70年代由日本育種學家選育的牡丹遠緣雜交品種,花乳黃色,花頭直立,芳香馥郁,生長勢強,品質(zhì)優(yōu)良,具有很高的觀賞與經(jīng)濟價值。但長期以來該品種的傳統(tǒng)繁殖方式(如分株、嫁接等)存在繁殖系數(shù)低、周期長、易受環(huán)境條件限制等問題,嚴重制約了其規(guī)?;a(chǎn)。微繁殖技術(shù)即植物的離體無性繁殖技術(shù),在保留母株優(yōu)良性狀的同時,具有繁殖快速、繁殖系數(shù)高等特點,能彌補傳統(tǒng)繁殖方法的不足,因此,構(gòu)建皇冠牡丹的微繁殖技術(shù)體系是推動其規(guī)模化生產(chǎn)的必然趨勢。

    增殖培養(yǎng)是牡丹微繁殖技術(shù)的第2個階段,其中試管苗能否高效地增殖和正常生長是決定微繁殖技術(shù)成敗的關(guān)鍵。國內(nèi)自1984年首次開展牡丹微繁殖技術(shù)研究以來[2],迄今已初步建立了菱花湛露(P. suffruticosa‘Ling Hua Zhan Lu’)、洛陽紅(P.×suffruticosa‘Luo Yang Hong’)、正午S9Weqz4lRk5OmoyZwW1O2Q==(P.×lemoinei ‘High Noon’)等40多個牡丹品種的增殖培養(yǎng)體系[3]。但目前關(guān)于皇冠牡丹增殖培養(yǎng)體系的研究仍鮮見報道,僅黃素姣[4]就植物生長調(diào)節(jié)劑對皇冠牡丹增殖生長的影響進行了探究,并獲得了適合皇冠牡丹增殖的植物生長調(diào)節(jié)劑組合,即6-BA(6-芐氨基嘌呤) 0.20mg/L+GA3(赤霉素) 0.10mg/L+ZT(玉米素) 0.05mg/L,但并未探究培養(yǎng)基中基礎(chǔ)成分對其增殖效果的影響。此外,已有研究發(fā)現(xiàn)牡丹增殖培養(yǎng)過程中的培養(yǎng)基成分與試管苗增殖效果關(guān)系密切,如Beruto等[5]發(fā)現(xiàn),培養(yǎng)基中的Ca2+濃度會對試管苗產(chǎn)生影響,并且提高WPM(木本植物專用培養(yǎng)基)中的Ca2+濃度,不僅能有效減少頂芽壞死的現(xiàn)象,還能提高增殖系數(shù)和莖長[6];Li等[7]則發(fā)現(xiàn)降低培養(yǎng)基中NH+4與NO-3濃度的比值可以促進試管苗的生長,并降低試管苗玻璃化率;在此基礎(chǔ)上,有研究將WPM中Ca(NO3)2濃度分別提高至原始濃度的3.0倍與4.0倍,從而顯著提高了正午牡丹(P.×lemoinei ‘High noon’)和鳳丹牡丹(P. ostii)試管苗的增殖系數(shù),且子代腋芽較為粗壯,而Ca(NO3)2濃度過低會導致試管苗生長不良、易玻璃化,Ca(NO3)2濃度過高則會降低增殖系數(shù),使葉片變小卷曲[8-9]。綜上,前人關(guān)于牡丹增殖培養(yǎng)階段的研究多集中在植物生長調(diào)節(jié)劑的篩選,較少對培養(yǎng)基配方進行研究,且培養(yǎng)基構(gòu)成組分復雜,各因素間會產(chǎn)生相互作用,使用常規(guī)試驗方法將耗費大量時間和精力。

    響應面分析法(RSM)通過結(jié)合數(shù)學方法和統(tǒng)計分析,對響應過程變量進行數(shù)學建模和分析,尋求最優(yōu)工藝參數(shù),是解決多變量問題的一種統(tǒng)計學方法[10],主要利用二次回歸模型對各因素和響應值之間的關(guān)系進行分析和擬合,試驗精度高、次數(shù)少、周期短,克服了正交試驗與單因素試驗法的不足,已被廣泛應用于食品學、生態(tài)學等領(lǐng)域,并常被用來改良植物組織培養(yǎng)過程中培養(yǎng)基成分與培養(yǎng)條件[11-14]。在牡丹微繁殖技術(shù)研究中,響應面法也有應用,黃弄璋[15]使用響應面法中的Plackett-Burman試驗,以WPM中的大量元素、微量元素、肌醇、鈣、有機物質(zhì)、鐵鹽6大組分進行設(shè)計,篩選發(fā)現(xiàn)鐵鹽、鈣和肌醇是影響正午牡丹增殖效果的主要影響因子,并選取這3個因素設(shè)計Box-Behnken(效應面法)試驗對培養(yǎng)基進行優(yōu)化,將增殖系數(shù)提高了10%左右。這表明響應面法在牡丹微繁殖技術(shù)培養(yǎng)基配方改良中具有較大的應用潛力,但目前尚未見使用響應面法對牡丹微繁殖技術(shù)培養(yǎng)基中的大量元素濃度進行改良調(diào)整的研究。鑒于此,本研究擬以皇冠牡丹試管苗為試驗材料,利用響應面法對WPM中5種大量元素配方進行優(yōu)化,以期篩選獲得適宜皇冠牡丹的特定增殖培養(yǎng)基,為該品種規(guī)?;a(chǎn)提供技術(shù)支持,并為其他牡丹品種微繁殖技術(shù)的研發(fā)提供參考與借鑒。

    1 材料與方法

    1.1 試驗材料

    試驗于2022年2月底于山東省菏澤市古韻牡丹基地進行,選取長勢良好、無病蟲害的皇冠牡丹母株,剪取帶有飽滿鱗芽的枝,剝下腋芽作為外植體,并按照文獻[16]的方法對外植體進行消毒后進行接種,啟動培養(yǎng)基為改良WPM[Ca(NO3)2·4H2O 1 668 mg/L+BA 0.5 mg/L+GA3 0.1 mg/L,附加蔗糖30.0 g/L、瓊脂6.5 g/L,pH=5.9],后續(xù)如無特殊說明,培養(yǎng)基同上。在第5次繼代培養(yǎng)末期,切取芽叢上莖長約為1~2 cm的健壯單芽,去除全部葉片作為試驗材料。

    1.2 試驗方法

    以WPM中5種大量元素的無機鹽[K2SO4、MgSO4·7H2O、KH2PO4、NH4NO3、Ca(NO3)2·4H2O]作為設(shè)計因素,前4種大量元素的無機鹽(K2SO4、MgSO4·7H2O、KH2PO4、NH4NO3)的濃度范圍設(shè)置為WPM原濃度[K2SO4為990 mg/L,MgSO4·7H2O為370 mg/L,KH2PO4為170 mg/L,NH4NO3為400 mg/L,Ca(NO3)2·4H2O為556 mg/L]的0.5~3.0倍,Ca(NO3)2·4H2O的濃度范圍設(shè)置為WPM原濃度的0.5~4.0倍,使用Design-Expert Optimal軟件中的響應面法進行最佳優(yōu)化試驗設(shè)計,得到30個處理,并以WPM原濃度作為對照(表1)。上述各處理培養(yǎng)基均附加6-BA 0.5 mg/L、GA3 0.1 mg/L、蔗糖30.0 g/L、瓊脂6.5 g/L,pH 5.9。

    1.3 培養(yǎng)條件與數(shù)據(jù)統(tǒng)計分析

    本研究中啟動和增殖培養(yǎng)條件均為培養(yǎng)溫度(24±1) ℃,光照時間14 h/d,熒光燈光照度為32.4 μmol/(m2·s)。

    增殖培養(yǎng)50 d后,統(tǒng)計不同處理試管苗的增殖系數(shù)(莖長為1~2 cm健壯單芽數(shù)/接種外植體個數(shù))、株高、莖長和葉片數(shù),并以上述指標為目標函數(shù)建立模型,根據(jù)模型計算預測目標值對應各因素的置信區(qū)間,獲得優(yōu)化配方,從而確定優(yōu)化培養(yǎng)基的組成比例,并進行驗證試驗。

    2 結(jié)果與分析

    將原始數(shù)據(jù)導入SPSS22.0,可得增殖系數(shù)、平均株高、平均莖長與單株葉片數(shù)(表2)。由表2可知,增殖系數(shù)較高的前3個處理依次為第13組、第28組、第15組;平均株高較高的前3個處理依次為第19組、第28組、第6組;平均莖長較高的前3個處理依次為第28組、第15組、第14組;單株葉片數(shù)較高的前3個處理依次為第15組、第19組、第13組。綜合來看,第28組為最優(yōu)處理,其增殖系數(shù)、平均株高、平均莖長與單株葉片數(shù)均優(yōu)于對照,說明培養(yǎng)基大量元素濃度對試管苗增殖效果會產(chǎn)生顯著影響。

    將SPSS22.0處理后的數(shù)據(jù)通過軟件Design-Expert 12.0進行數(shù)據(jù)擬合和顯著性檢驗,對試驗結(jié)果與試驗因子進行三維建模,預測增殖系數(shù)和平均莖長處于較高水平時培養(yǎng)基成分的最優(yōu)配比。其中,5項變量因子為K2SO4、MgSO4·7H2O、KH2PO4、NH4NO3、Ca(NO3)2·4H2O的濃度倍數(shù),分別用x1、x2、x3、x4、x5表示。線性項顯著表示該大量元素濃度對響應值影響顯著,二次項顯著表示2種大量元素濃度交互作用顯著,系數(shù)為正、負分別表示濃度增大時影響方向的正、負。

    2.1 平均增殖系數(shù)模型建立與分析

    以增殖系數(shù)為Y1,根據(jù)Design-Expert 12.0軟件擬合試驗處理,得到增殖系數(shù)的二次回歸方程:

    Y1=1.414 280+0.001 072x1-0.030 090x2+0.398 823x3-0.016 274x4+0.224 442x5+0.037 415x1x2+0.035 269x1x3+0.036 119x1x4-0.019 207x1x5+0.016 060x2x3-0.049 592x2x4-0.002 236x2x5-0.024 230x3x4+0.022 937x3x5-0.002 875x4x5-0.038 994x21+0.001 120x22-0.113 775x23-0.003 998x24-0.035 254x25。

    分析結(jié)果(表3)表明,線性項x3、x4對增殖系數(shù)具有極顯著影響,線性項x5對增殖系數(shù)具有顯著影響,由F(x3)>F(x4)>F(x5)>F(x2)>F(x1)可知,5種大量元素濃度對增殖系數(shù)的影響程度為KH2PO4>NH4NO3>Ca(NO3)2·4H2O>MgSO4·7H2O>K2SO4。具體影響(圖1)如下: KH2PO4、Ca(NO3)2·4H2O濃度倍數(shù)與增殖系數(shù)呈線性正相關(guān)(圖1C;圖1E),相關(guān)系數(shù)分別為0.398 823與0.224 442;NH4NO3濃度倍數(shù)與增殖系數(shù)呈線性負相關(guān)(圖1D),相關(guān)系數(shù)為-0.016 274;K2SO4、MgSO4·7H2O濃度倍數(shù)與增殖系數(shù)無顯著相關(guān)性(圖1A、圖1B)。

    取F值>1.5的二元二次項分析各因素間的交互作用與對應響應曲面(圖2)。比較F值大小可知,x1x2、x1x3、x1xiQpEOpm0Pdu9NRAwmbSco02ONOfP4j2cipENJgXQAQY=4和x2x4這4項對增殖系數(shù)的影響程度為:x2x4>x1x3>x1x4>x1x2(表3),即MgSO4·7H2O與NH4NO3的交互作用>K2SO4與KH2PO4的交互作用>K2SO4與NH4NO3的交互作用>K2SO4與MgSO4·7H2O的交互作用。影響如下:(1)K2SO4與MgSO4·7H2O濃度比保持為1∶1時,試管苗增殖系數(shù)未見明顯變化;但低濃度K2SO4與高濃度MgSO4·7H2O相互作用會導致增殖系數(shù)下降(圖2A)。(2)K2SO4與KH2PO4濃度同時升高,增殖系數(shù)先上升后下降,且增殖系數(shù)受KH2PO4濃度影響較大(圖2B)。(3)K2SO4與NH4NO3兩者濃度同時升高,增殖系數(shù)逐漸下降且受NH4NO3濃度影響較大;低濃度NH4NO3可以使增殖系數(shù)保持較高水平(圖2C)。(4)當MgSO4·7H2O與NH4NO3濃度同時升高時,試管苗增殖系數(shù)逐漸降低,而高濃度MgSO4·7H2O與低濃度NH4NO3相互作用時增殖系數(shù)可達最高值(圖2D)。

    2.2 平均株高模型建立與分析

    以平均株高為Y2,根據(jù)Design-Expert 12.0軟件擬合試驗處理,得到株高的二次回歸方程:

    Y2=2.998 730-0.263 603x1+0.588 062x2+0.264 487x3-0.161 580x4+0.462 537x5-0.007 257x1x2+0.135 470x1x3+0.037 222x1x4-0.031 871x1x5+0.097 298x2x3-0.021 955x2x4-0.046 501x2x5-0.119 603x3x4-0.031 871x3x5-0.042 601x4x5-0.027 184x21-0.178 435x22-0.053 012x23+0.028 309x24+0.000 039x25。

    分析結(jié)果(表4)表明,線性項x3、x4、x5對株高均有極顯著影響,二次項x1x3對株高有顯著影響。由F(x4)>F(x5)>F(x3)>F(x1)>F(x2)可知,5種大量元素濃度對平均株高的影響程度為NH4NO3>Ca(NO3)2·4H2O>KH2PO4>K2SO4>MgSO4·7H2O。具體影響如下:(1)KH2PO4、Ca(NO3)2·4H2O濃度倍數(shù)與株高呈線性正相關(guān)(圖3C;圖3E),相關(guān)系數(shù)分別為0.264 487與0.462 537;NH4NO3濃度倍數(shù)與株高呈線性負相關(guān),相關(guān)系數(shù)為-0.161 580(圖3D);K2SO4與MgSO4·7H2O濃度倍數(shù)與株高無顯著相關(guān)性(圖3A、圖3B)。(2)K2SO4與KH2PO4濃度倍數(shù)的交互作用與株高呈線性正相關(guān)(圖3F),相關(guān)系數(shù)為0.135 470。

    取F值>1.0的二元二次項分析交互作用與響應曲面。比較F值大小可知,x1x3、x2x3、x2x5和x3x4這4項對株高的影響程度為:x1x3>x3x4>x2x3>x2x5(表4),即K2SO4與KH2PO4的交互作用>KH2PO4與NH4NO3的交互作用>MgSO4·7H2O與KH2PO4的交互作用>MgSO4·7H2O與Ca(NO3)2·4H2O的交互作用。影響如下:(1)當K2SO4與KH2PO4濃度同時升高時,株高無明顯變化,當K2SO4濃度較高且KH2PO4濃度較低時,株高呈下降趨勢(圖4A)。(2)當MgSO4·7H2O濃度倍數(shù)為0.5~1.8倍時,株高隨KH2PO4濃度升高而上升;當濃度倍數(shù)在1.8~3.0倍時,株高隨KH2PO4濃度升高而下降(圖4B)。(3)MgSO4·7H2O濃度倍數(shù)在0.5~2.0倍時,株高隨Ca(NO3)2·4H2O濃度升高而上升;當濃度倍數(shù)在2.0~3.0倍時,株高隨Ca(NO3)2·4H2O濃度升高而下降(圖4C)。(4)當KH2PO4濃度升高、NH4NO3濃度下降時,株高呈上升趨勢并在KH2PO4濃度倍數(shù)為3.0倍、NH4NO3濃度倍數(shù)為0.5倍時達到最高值(圖4D)。

    2.3 平均莖長模型建立與分析

    以平均莖長為Y3,根據(jù)Design-Expert 12.0軟件擬合試驗處理,得到莖長的二次回歸方程:

    Y3=1.301 670+0.148 644x1-0.196 973x2+0.071 288x3-0.193 460x4+0.116 001x5+0.060 904x1x2+0.062 311x1x3+0.005 765x1x4+0.009 485x1x5+0.020 522x2x3-0.023 159x2x4+0.021 980x2x5-0.050 456x3x4+0.023 133x3x5+0.002 436x4x5-0.027 867x12+0.008 299x22-0.028 877x32-0.035 352x42-0.042 483x52。

    分析結(jié)果(表5)表明,線性項x3、二次項x1x2和x1x3對平均莖長具有極顯著影響,線性項x4、二次項x3x4、x25對莖長具有顯著影響。由F(x3)>F(x4)>F(x5)>F(x2)>F(x1)可知,5種大量元素濃度對平均莖長的影響程度為KH2PO4>NH4NO3>Ca(NO3)2·4H2O>MgSO4·7H2O>K2SO4。具體影響如下:(1)KH2PO4濃度倍數(shù)與莖長呈正線性相關(guān)(圖5C),相關(guān)系數(shù)為0.071 288;NH4NO3濃度倍數(shù)與莖長呈負線性相關(guān)(圖5D),相關(guān)系數(shù)為-0.193 460;K2SO4、MgSO4·7H2O濃度倍數(shù)與莖長無顯著相關(guān)性(圖5A、圖5B)。(2)K2SO4與MgSO4·7H2O濃度倍數(shù)的交互作用與莖長呈正線性相關(guān)(圖5F),相關(guān)系數(shù)為0.060 904;K2SO4與KH2PO4濃度倍數(shù)的交互作用與莖長呈正線性相關(guān)(圖5G),相關(guān)系數(shù)為0.062 311;NH4NO3與KH2PO4濃度倍數(shù)的交互作用與莖長呈負線性相關(guān)(圖5H),相關(guān)系數(shù)為-0.050 456。(3)Ca(NO3)2·4H2O濃度倍數(shù)對莖長無顯著影響(圖5E),但其平方值對莖長具有顯著影響,相關(guān)系數(shù)為-0.042 483,莖長隨Ca(NO3)2·4H2O濃度倍數(shù)增大呈現(xiàn)先上升后下降的趨勢。

    取F值>3.0的二元二次項分析交互作用與響應曲面。比較F值大小可知,x1x2、x1x3、x3x4和x3x5這4項交互作用對莖長的影響程度順序為:x1x3>x1x2>x3x4>x3x5,即K2SO4與KH2PO4的交互作用>K2SO4與MgSO4·7H2O的交互作用>KH2PO4與NH4NO3的交互作用>KH2PO4與Ca(NO3)2·4H2O的交互作用。具體影響如下:(1)K2SO4與MgSO4·7H2O濃度為1∶1且同時升高時,莖長先下降后升高(圖6A)。(2)K2SO4與KH2PO4濃度同時增加時,莖長上升,在兩者同時達到3.0倍時莖長達到最高值(圖6B)。(3)KH2PO4與NH4NO3濃度同時升高時,莖長先升高后下降;當KH2PO4處于高濃度且NH4NO3處于低濃度時,試管苗莖長最高(圖6C)。(4)KH2PO4與Ca(NO3)2·4H2O濃度同時升高時,莖長隨之升高;當KH2PO4濃度較高且Ca(NO3)2·4H2O濃度倍數(shù)約為3.0倍時,莖長達到最高值(圖6D)。

    2.4 WPM中的5種大量元素配方優(yōu)化及驗證

    通過Design-Expert 12.0建立模型,將皇冠牡丹試管苗增殖系數(shù)目標設(shè)為max,重要程度為5(+++++);莖長目標設(shè)為max,重要程度為4(++++);株高目標設(shè)為max,重要程度為3(+++);葉片數(shù)目標設(shè)為max,重要程度為2(++)。根據(jù)模型分析結(jié)果,可得出適宜皇冠牡丹增殖培養(yǎng)的WPM中的5種大量元素的最佳濃度倍數(shù)分別為x1=2.28、x2=2.45、x3=2.83、x4=0.50、x5=3.53,即K2SO4、MgSO4·7H2O、KH2PO4、NH4NO3與Ca(NO3)2·4H2O的質(zhì)量濃度分別為2 257 mg/L、907 mg/L、481 mg/L、200 mg/L與1 963 mg/L,此時培養(yǎng)基中離子濃度分別為NO-3 19.13 mmol/L,NH+4 2.50 mmol/L,Ca2+ 8.31 mmol/L,K+ 29.44 mmol/L,Mg2+ 3.68 mmol/L,PO3-4 3.53 mmol/L,SO2-4 16.63 mmol/L,并通過模型預測出在該條件下皇冠牡丹的增殖系數(shù)、株高、莖長與單株葉片數(shù)分別為2.24、4.32 cm、1.64 cm、5.31張。使用該配方進行驗證試驗(圖7),得到皇冠牡丹的增殖系數(shù)、株高、莖長與單株葉片數(shù)分別為2.43、4.62 cm、1.30 cm、6.75張,除莖長稍低外,其余指標均高于預測值,并顯著優(yōu)于對照組(1.80、3.84 cm、1.32 cm、4.32張),這說明該模型可信,本研究得出的5種大量元素的最佳濃度可以用于皇冠牡丹的增殖培養(yǎng)。

    3 討論

    3.1 單因素條件對皇冠牡丹增殖和生長的影響

    在對植物進行離體培養(yǎng)時,植物生長發(fā)育所需養(yǎng)分主要來源于培養(yǎng)基,基本培養(yǎng)基的成分與濃度對植物有重要影響[17]。本研究結(jié)果表明,WPM中的5種大量組分中,KH2PO4、NH4NO3與Ca(NO3)2·4H2O對牡丹試管苗的增殖效果有顯著或極顯著影響,而MgSO4·7H2O與K2SO4對牡丹試管苗的增殖效果無明顯影響。

    KH2PO4可同時為植物提供磷、鉀2種大量元素,在WPM、MS與DKW等植物培養(yǎng)基中均有添加,增加KH2PO4用量會促進植物生物量的累積[18]。本研究發(fā)現(xiàn),隨著KH2PO4濃度上升,試管苗的增殖系數(shù)、株高與莖長明顯提高,481 mg/L(2.83倍)KH2PO4為皇冠牡丹試管苗增殖最適濃度,而高濃度KH2PO4雖然會促進牡丹試管苗株高與莖長增加,但同時會降低增殖系數(shù)與葉片數(shù)。多數(shù)研究結(jié)果表明,不同植物適宜的KH2PO4濃度差異顯著,蘋果矮化砧木SH6繼代培養(yǎng)的最適濃度為175 mg/L[19],杏(Prunus armeniaca)試管苗培養(yǎng)的最適KH2PO4濃度為125 mg/L[12],而油松(Pinus tabuliformis)成熟離體胚不定芽誘導所需最適KH2PO4濃度為340 mg/L[20]。此外,同種植物的不同品種適宜的KH2PO4濃度也存在較大差異,如在使用MS的基礎(chǔ)上,KH2PO4濃度對白鶴芋屬(Spathiphyllum)不同品種增殖系數(shù)的影響存在差異[21]。因此,皇冠牡丹試管苗增殖階段的最佳KH2PO4質(zhì)量濃度為481 mg/L(2.83倍),但對于其他牡丹品種還需要根據(jù)實際情況具體分析。

    鈣元素在細胞分裂、細胞壁形成和分生組織生長過程中是非常重要的[22],前人研究發(fā)現(xiàn)高濃度Ca2+對牡丹試管苗的增殖效果有促進作用,將MS培養(yǎng)基中CaCl2濃度提升至原濃度的2倍可顯著提高Mme de Vatry牡丹試管苗的增殖系數(shù),并減輕莖尖壞死和玻璃化程度[6]。然而,不同牡丹品種的最適鈣鹽濃度具有較大差異,同樣以WPM作為增殖培養(yǎng)基,正午牡丹的最適Ca(NO3)2·4H2O濃度為1 668 mg/L(3.00倍),此時試管苗的增殖系數(shù)可達到3.0,子代腋芽較為粗壯[8];而鳳丹牡丹的最適Ca(NO3)2濃度為1 544 mg/L(4.00倍),此時增殖系數(shù)可達3.07[9],這可能是由牡丹品種間基因型不同所致。本研究發(fā)現(xiàn),隨著Ca(NO3)2·4H2O濃度升高,皇冠牡丹試管苗的增殖系數(shù)與莖長呈先升高后下降的趨勢,而株高則持續(xù)上升,最佳Ca(NO3)2·4H2O質(zhì)量濃度為1 963 mg/L(3.53倍),在該條件下試管苗較為健壯,分枝與莖葉數(shù)量較對照顯著提高,且莖尖壞死現(xiàn)象得到有效緩解,與前人研究結(jié)果一致。因此,在一定范圍內(nèi)提高WPM中Ca(NO3)2·4H2O濃度是提高牡丹試管苗增殖效果的有效途徑,而其具體的作用機理也值得后續(xù)研究深入探索。

    3.2 交互作用對皇冠牡丹試管苗有增殖和生長的影響

    本研究結(jié)果顯示,K2SO4與KH2PO4、K2SO4與MgSO4·7H2O、KH2PO4與NH4NO3的交互作用對試管苗增殖效果影響顯著或極顯著。因電荷存在動態(tài)平衡,培養(yǎng)基中大量元素間的交互作用實則為離子間的交互作用[23],因此本研究以分析離子間的交互作用與元素對植物的影響為主。

    K2SO4與KH2PO4的交互作用對株高、莖長均有顯著或極顯著影響,當二者濃度同時增大時,增殖系數(shù)呈先升高后下降的趨勢,株高與莖長呈升高趨勢。因二者的陽離子相同,其交互作用實則為SO2-4與PO3-4間的交互作用,SO2-4提供硫元素,含硫氨基酸不僅是構(gòu)成蛋白質(zhì)的主要成分,還是細胞內(nèi)某些化合物合成的介質(zhì)之一[24];PO3-4提供磷元素,磷是植物必需的大量元素,參與磷脂、核酸、腺嘌呤核苷三磷酸(ATP)和輔酶等大分子合成[25]。已有研究結(jié)果表明,施磷能顯著提高牡丹植株葉片的葉綠素含量并促進光合作用[26],而硫是植物光合作用過程中的重要組分,二者均與光合作用有關(guān);Misson等[27]發(fā)現(xiàn),在低磷環(huán)境下植物會加大對硫酸鹽的吸收轉(zhuǎn)運以合成硫脂,而硫脂可以在一定程度上代替磷脂,說明硫、磷元素在功能上存在重疊性,存在相互影響。目前關(guān)于SO2-4與PO3-4在植物體內(nèi)的交互作用尚未見報道,但從本研究結(jié)果可以看出二者存在協(xié)同的交互作用,適宜的SO2-4與PO3-4濃度比會促進試管苗的增殖。但本研究也發(fā)現(xiàn),當二者濃度同時超過2.50倍后,增殖系數(shù)與葉片數(shù)會出現(xiàn)下降趨勢,這可能是過高的離子濃度對植物產(chǎn)生了毒害作用。因此,在提高培養(yǎng)基內(nèi)營養(yǎng)物質(zhì)濃度以促進試管苗增殖的過程中,應注意保持培養(yǎng)基中離子濃度平衡,避免過高的濃度對植物生長產(chǎn)生負面影響。

    本研究結(jié)果顯示,K2SO4與MgSO4·7H2O的交互作用對莖長有極顯著影響,當兩者濃度同時上升時,莖長呈先下降后上升的趨勢。已有研究結(jié)果表明鉀元素與鎂元素存在交互作用,高濃度鉀會抑制植物對鎂的吸收,并導致植物因缺鎂而生長不良[28]。這可能與二者相似的化學性質(zhì)產(chǎn)生的拮抗作用有關(guān)[29],而對細胞膜上相同轉(zhuǎn)運蛋白的競爭是二者產(chǎn)生拮抗作用的主要原因[30],且鉀、鎂間的交互作用會影響植物對其他養(yǎng)分如氮、磷、鈣等的積累[31-32]。本研究中當K2SO4與MgSO4·7H2O濃度倍數(shù)介于0.5~2.0時,莖長呈下降趨勢,這與前人研究發(fā)現(xiàn)的鉀、鎂元素間的拮抗作用相符合;而當二者濃度倍數(shù)大于2.0后,莖長呈上升趨勢,這可能是由于更高濃度的鉀、鎂元素能抵消一部分拮抗作用對試管苗產(chǎn)生的負面影響。因此,在改良培養(yǎng)基中大量元素配方時應重視礦質(zhì)元素間的平衡關(guān)系,盡量避免鉀、鎂元素間產(chǎn)生拮抗作用,而關(guān)于鉀、鎂元素在牡丹體內(nèi)的作用機理則需要進一步研究。

    本研究結(jié)果顯示,KH2PO4與NH4NO3間的交互作用對莖長有顯著影響,當KH2PO4濃度上升、NH4NO3濃度下降時,莖長呈升高趨勢。前人研究發(fā)現(xiàn)植物對二者需求的差異可能與基因型、統(tǒng)計方式和使用濃度范圍有關(guān)[33-35]。Akin等[11]在對榛子(Corylus heterophylla)所使用的DKW培養(yǎng)基中的離子濃度進行響應面法優(yōu)化時也發(fā)現(xiàn)相似情況,即當KH2PO4為高濃度、NH4NO3為低濃度時,榛子Wepster能得到較高的芽質(zhì)量,榛子Dorris和Wepster能達到理想芽長。這可能是由于NH4NO3與KH2PO4在產(chǎn)生交互作用時,隨著兩者濃度升高,NH4NO3與Ca(NO3)2·4H2O同時產(chǎn)生了反應,而Ca(NO3)2·4H2O和NH4NO3同時為試管苗提供氮元素,降低NH4+與NO-3比值能夠顯著促進牡丹試管苗增殖[7-9],本研究據(jù)此優(yōu)化得出結(jié)果,將NH4NO3濃度降低為基本W(wǎng)PM中濃度的0.50倍后,增殖效果得到顯著提升,與前人研究結(jié)果相似。

    3.3 銨態(tài)氮(NH+4-N)和硝態(tài)氮(NO-3-N)濃度比值對皇冠牡丹增殖的影響

    氮是植物生長過程中需求量最大的礦質(zhì)元素,是組成蛋白質(zhì)和核酸的必要成分。WPM為木本植物專用培養(yǎng)基[36],供給植物的主要氮源為銨態(tài)氮(NH+4-N)和硝態(tài)氮(NO-3-N),其中NH4NO3提供NH+4-N,NH4NO3和Ca(NO3)2·4H2O共同提供NO-3-N。二者雖共同提供氮元素,但對植物生長的作用存在差異,關(guān)于培養(yǎng)基中NH+4-N與NO-3-N的濃度比值對試管苗增殖效果的影響已有諸多研究報道,結(jié)果顯示,較低的NH+4-N與NO-3-N的濃度比會促進試管苗的增殖與生長[11, 37-40]。在牡丹研究中也發(fā)現(xiàn)類似現(xiàn)象,例如Li等[7]通過降低MS培養(yǎng)基中NH+4-N與NO-3-N濃度比值以促進牡丹試管苗增殖,而后有試驗通過增加培養(yǎng)基中Ca(NO3)2·4H2O濃度,在NH+4-N與NO-3-N的濃度比值降低的同時提高培養(yǎng)基中Ca2+的濃度,顯著提高了牡丹試管苗的增殖效果[8-9]。這表明降低培養(yǎng)基中NH+4-N與NO-3-N的濃度比值有利于牡丹試管苗的增殖,這可能是由于高濃度NH4+會促進植物吸收過多水分,從而導致試管苗枝條壞死、存活率降低[40-41],同時較低的NH+4-N與NO-3-N的濃度比值可以更好地緩沖培養(yǎng)基的pH,減少培養(yǎng)基pH波動對試管苗的影響[42]。本研究也得到類似結(jié)果,將Ca(NO3)2·4H2O濃度提升至3.50倍,NH4NO3降低至0.50倍后,試管苗增殖系數(shù)得到顯著提升,并且莖尖壞死現(xiàn)象得到緩解。綜上可見,WPM中原有的NH+4-N與NO-3-N的濃度比值無法滿足牡丹試管苗的增殖需要,在一定范圍內(nèi)提高培養(yǎng)基中Ca(NO3)2·4H2O濃度并降低NH4NO3濃度,是促進皇冠牡丹試管苗增殖生長的有效途徑。當然,不同牡丹品種的適宜濃度存在較大差異,關(guān)于NH+4-N與NO-3-N的濃度比值對牡丹試管苗產(chǎn)生影響的機理研究還尚在探索階段,還需進一步研究。

    4 結(jié)論

    本研究基于響應面法,優(yōu)化得到皇冠牡丹增殖培養(yǎng)專用WPM大量元素配方[K2SO4 2 257 mg/L、MgSO4·7H2O 907 mg/L、KH2PO4 481 mg/L、NH4NO3 200 mg/L、Ca(NO3)2·4H2O 1 963 mg/L],建立了高效的皇冠牡丹增殖培養(yǎng)體系,在該條件下試管苗的增殖效率(增殖系數(shù)2.43、株高4.62 cm、莖長1.30 cm、葉片數(shù)6.75張)得到顯著提高。

    參考文獻:

    [1] LI S S, YUAN R Y, CHEN L G, et al. Systematic qualitative and quantitative assessment of fatty acids in the seeds of 60 tree peony (Paeonia section Moutan DC.) cultivars by GC-MS[J]. Food Chemistry,2015,173(15):133-140.

    [2] 李玉龍,吳德玉,潘淑龍,等. 牡丹試管苗繁殖技術(shù)的研究[J]. 科學通報,1984(8):500-502.

    [3] WEN S S, CHEN L, TIAN R N. Micropropagation of tree peony (Paeonia sect. Moutan):a review[J]. Plant Cell Tissue and Organ Culture,2020,141(4):1-14.

    [4] 黃素姣. 三個牡丹品種的離體快繁技術(shù)研究[D]. 北京:北京林業(yè)大學,2020.

    [5] BERUTO M, CURIR P. In vitro culture of tree peony through axillary budding[M]. Berlin:Springer Netherlands,2007:477-497.

    [6] BOUZA L, JACQUES M, MIGINIAC E. In vitro propagation of Paeonia suffruticosa Andr. cv. ‘Mme de Vatry’: developmental effects of exogenous hormones during the multiplication phase[J]. Scientia Horticulturae,1994,3(57):241-251.

    [7] LI P, CHENG P Y. Basal medium with modified calcium source and other factors influence on shoots culture of tree peony[J]. Acta Hortic,2008(766):383-390.

    [8] 文書生,成仿云,鐘 原,等. ‘正午’牡丹微繁殖體系的建立[J]. 植物科學學報,2016,34(1):143-150.

    [9] 王 新,成仿云,鐘 原,等. 鳳丹牡丹鱗芽離體培養(yǎng)與快繁技術(shù)[J]. 林業(yè)科學,2016,52(5):101-110.

    [10]武愛龍,何 冰,吳建陽,等. 響應面法在植物組織培養(yǎng)中的應用進展[J]. 安徽農(nóng)學通報,2020,26(17):21-22.

    [11]AKIN M, EYDURAN E, REED B M. Use of RSM and CHAID data mining algorithm for predicting mineral nutrition of hazelnut[J]. Plant Cell, Tissue and Organ Culture,2017,128(2):303-316.

    [12]KOVALCHUK I Y, MUKHITDINOVA Z, TURDIYEV T, et al. Modeling some mineral nutrient requirements for micropropagated wild apricot shoot cultures[J]. Plant Cell, Tissue and Organ Culture,2017,129(2):325-335.

    [13]曹 劍,黃志偉,李華,等. 利用響應面法優(yōu)化毛葉木姜子莖段腋芽誘導條件[J]. 分子植物育種,2022,20(21):7209-7215.

    [14]蔡正禹,文書生,田如男. 響應面法優(yōu)化歐洲冬青‘Ferox Argentea’增殖培養(yǎng)基大量元素配方[J]. 東北林業(yè)大學學報,2023,51(6):53-63.

    [15]黃弄璋. ‘正午’與‘鳳丹’牡丹的離體增殖與生根移栽技術(shù)研究[D]. 北京:北京林業(yè)大學,2017.

    [16]邱金梅. 牡丹離體快繁技術(shù)的研究[D]. 北京:北京林業(yè)大學,2010.

    [17]文書生,何絨絨,鄭佳康,等. 牡丹組織培養(yǎng)技術(shù)研究進展[J]. 林業(yè)科學,2018,54(10):143-155.

    [18]EL-HAWAZ R F, BRIDGES W C, ADELBERG J W, et al. In Vitro growth of Curcuma longa L. in response to five mineral elements and plant density in fed-batch culture systems[J]. PLoS One,2015,10(4):e0118912.

    [19]楊雨璋,周貝貝,李民吉,等. 蘋果矮化砧木‘SH6’組培快繁培養(yǎng)基大量元素配方的優(yōu)化[J]. 果樹學報,2020,37(1):40-49.

    [20]董麗芬,肖 穎,邵崇斌. 氮、磷、鉀元素形態(tài)配比及濃度對油松胚培養(yǎng)的影響[J]. 西北林學院學報,2006,21(3):64-66.

    [21]朱根發(fā). 白鶴芋屬觀賞植物的組織培養(yǎng)和快速繁殖技術(shù)研究[J]. 中國農(nóng)學通報,2003,19(3):75-76.

    [22]HIRSCHI K D. The calcium conundrum. Both versatile nutrient and specific signal[J]. Plant Physiology,2004,136(1):2438-2442.

    [23]NIEDZ R P, EVENS T J .The effects of nitrogen and potassium nutrition on the growth of nonembryogenic and embryogenic tissue of sweet orange (Citrus sinensis (L.) Osbeck)[J]. Bmc Plant Biology,2008,8(1):126.

    [24]王慶仁,林 葆. 植物硫營養(yǎng)研究的現(xiàn)狀與展望[J]. 中國土壤與肥料,1996(3):16-19.

    [25]PANG J Y, RYAN H M, HANS L, et al. Phosphorus acquisition and utilisation in crop legumes under global change[J]. Current Opinion in Plant Biology,2018,45:248-254.

    [26]單佩佩. 牡丹體內(nèi)外磷循環(huán)和根際微生物種群構(gòu)成對外源磷素的響應[D]. 泰安:山東農(nóng)業(yè)大學,2016.

    [27]MISSON J, RAGHOTHAMA K G, JAIN A, et al. A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation[J]. Proceedings of the National Academy of Sciences,2005,102(33):11934-11939.

    [28]LI H X, CHEN Z J, ZHOU T, et al. High potassium to magnesium ratio affected the growth and magnesium uptake of three tomato (Solanum lycopersicum L.) cultivars[J]. Journal of Integrative Agriculture,2018,17(12):2813-2821.

    [29]VIAD A, FERNNDEZ-MARCOS M L, HERNNDEZ-NISTAL J, et al. Effect of particle size of limestone on Ca, Mg and K contents in soil and in sward plants[J]. Scientia Agricola,2011,68:200-208.

    [30]KAMIYA T, YAMAGAMI M, HIRAI M Y, et al. Establishment of an in planta magnesium monitoring system using CAX3 promoter-luciferase in Arabidopsis[J]. Journal of Experimental Botany,2012,63(1):355-363.

    [31]薛欣欣,魏云霞,王文斌,等. 鉀、鎂交互作用對橡膠幼苗生長及養(yǎng)分吸收的影響[J]. 植物營養(yǎng)與肥料學報,2020,26(10):1870-1878.

    [32]FARHAT N, ELKHOUNI A, ZORRIG W, et al. Effects of magnesium deficiency on photosynthesis and carbohydrate partitioning[J]. Acta Physiologiae Plantarum,2016,38(6):145.

    [33]HAND C, MAKI S, REED B M. Modeling optimal mineral nutrition for hazelnut micropropagation[J]. Plant Cell, Tissue and Organ Culture,2014,119(2):411-425.

    [34]AKIN M, EYDURAN S P, EYDURAN E, et al. Analysis of macro nutrient related growth responses using multivariate adaptive regression splines[J]. Plant Cell, Tissue and Organ Culture,2020,140(3):661-670.

    [35]REED B M, WADA S, DENOMA, J. et al. Improving in vitro mineral nutrition for diverse pear germplasm.[J]. In Vitro Cellular & Developmental Biology-Plant,2013,49(3):343-355.

    [36]LLOYD G, MCCOWN B. Commercially-feasible micropropagation of mountain laurel, Kalmia latifolia, by use of shoot-tip culture[J]. Combined Proceedings - International Plant Propagators’ Society (USA),1980,30:421-427.

    [37]ALANAGH E N, GAROOSI G, HADDAD R, et al. Design of tissue culture media for efficient Prunus rootstock micropropagation using artificial intelligence models[J]. Plant Cell, Tissue and Organ Culture,2014,117(3):349-359.

    [38]PARK K, JANG B K, LEE H M, et al. Effective Propagation of Selaginella tamariscina through optimized medium composition[J]. Agronomy,2021,11(3):578.

    [39]JAMSHIDI S, YADOLLAHI A, AHMADI H, et al. Predicting In vitro culture medium macro-nutrients composition for pear rootstocks using regression analysis and neural network models[J]. Frontiers in Plant Science,2016,7. DOI:10.3389/fpls.2016.00274.

    [40]POOTHONG S, REED B M. Optimizing shoot culture media for Rubus germplasm: the effects of NH+4, NO-3, and total nitrogen[J]. In Vitro Cellular & Developmental Biology - Plant, 2016, 52(3):265-275.

    [41]GAGO J, PREZ-TORNERO O, LANDN M, et al. Improving knowledge of plant tissue culture and media formulation by neurofuzzy logic: a practical case of data mining using apricot databases[J]. Journal of Plant Physiology,2011,168(15):1858-1865.

    [42]MARTIN S M, ROSE D, HUI V. Growth of plant cell suspension cultures with ammonium as the sole sou. [J]. Canadian Journal of Botany,1977,55(22):2838-2843.

    (責任編輯:陳海霞)

    久久中文看片网| av视频在线观看入口| 黑人巨大精品欧美一区二区mp4| 久久久久国产精品人妻aⅴ院| 亚洲色图综合在线观看| 伦理电影免费视频| 免费看十八禁软件| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 人成视频在线观看免费观看| 久久精品影院6| 久久天堂一区二区三区四区| 日韩中文字幕欧美一区二区| 少妇熟女aⅴ在线视频| 亚洲av片天天在线观看| 日韩欧美国产在线观看| 可以免费在线观看a视频的电影网站| 国产高清有码在线观看视频 | 欧美日本中文国产一区发布| 岛国视频午夜一区免费看| 精品卡一卡二卡四卡免费| 久久人人爽av亚洲精品天堂| 18美女黄网站色大片免费观看| 黄片小视频在线播放| 老司机靠b影院| 国产伦一二天堂av在线观看| 色尼玛亚洲综合影院| 免费在线观看黄色视频的| 国产精品精品国产色婷婷| 国产又色又爽无遮挡免费看| 亚洲熟女毛片儿| 色播在线永久视频| 欧美中文综合在线视频| 波多野结衣一区麻豆| 国内毛片毛片毛片毛片毛片| 久久人妻熟女aⅴ| 巨乳人妻的诱惑在线观看| 免费无遮挡裸体视频| 亚洲av成人一区二区三| 欧美性长视频在线观看| 黄色成人免费大全| 国产精品免费视频内射| 9色porny在线观看| 国产成人av教育| 18禁美女被吸乳视频| а√天堂www在线а√下载| 亚洲五月天丁香| 久久精品成人免费网站| 十八禁人妻一区二区| 叶爱在线成人免费视频播放| 国产一卡二卡三卡精品| 日韩欧美三级三区| 国产精品九九99| 搡老妇女老女人老熟妇| 男女之事视频高清在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区精品视频观看| 久久国产亚洲av麻豆专区| 成人精品一区二区免费| 无限看片的www在线观看| 国产乱人伦免费视频| 正在播放国产对白刺激| 精品无人区乱码1区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲无线在线观看| 91麻豆精品激情在线观看国产| 丁香欧美五月| 嫩草影视91久久| 国产视频一区二区在线看| 激情在线观看视频在线高清| 777久久人妻少妇嫩草av网站| 亚洲在线自拍视频| 国产av一区在线观看免费| 九色国产91popny在线| 久久人妻av系列| 亚洲成av人片免费观看| 两性夫妻黄色片| 日本撒尿小便嘘嘘汇集6| 青草久久国产| 两性夫妻黄色片| 无人区码免费观看不卡| 精品无人区乱码1区二区| 欧美另类亚洲清纯唯美| 人人澡人人妻人| 色av中文字幕| 欧美日本中文国产一区发布| 一级黄色大片毛片| 黄网站色视频无遮挡免费观看| 精品国产超薄肉色丝袜足j| 身体一侧抽搐| 国产精品综合久久久久久久免费 | 黄网站色视频无遮挡免费观看| 精品国产亚洲在线| 国产视频一区二区在线看| 变态另类成人亚洲欧美熟女 | 啦啦啦 在线观看视频| 国产精品久久久人人做人人爽| 1024香蕉在线观看| 久久精品影院6| 国产伦人伦偷精品视频| 99在线人妻在线中文字幕| 97人妻精品一区二区三区麻豆 | 中文字幕色久视频| 麻豆国产av国片精品| 欧美在线一区亚洲| 叶爱在线成人免费视频播放| 99国产精品免费福利视频| 久久中文看片网| av有码第一页| 日本a在线网址| 黄色女人牲交| 国产xxxxx性猛交| 91字幕亚洲| 在线视频色国产色| 亚洲激情在线av| 亚洲九九香蕉| 精品欧美国产一区二区三| 亚洲国产精品成人综合色| 国产精品久久视频播放| 欧美日韩一级在线毛片| 中出人妻视频一区二区| 亚洲情色 制服丝袜| 搡老熟女国产l中国老女人| 日日夜夜操网爽| 久久伊人香网站| 国产aⅴ精品一区二区三区波| 亚洲 欧美一区二区三区| 欧美绝顶高潮抽搐喷水| 亚洲人成网站在线播放欧美日韩| 午夜日韩欧美国产| 中文字幕高清在线视频| 啦啦啦韩国在线观看视频| 免费在线观看日本一区| 欧美久久黑人一区二区| 亚洲情色 制服丝袜| 看片在线看免费视频| 久久精品亚洲熟妇少妇任你| 搞女人的毛片| 侵犯人妻中文字幕一二三四区| www.自偷自拍.com| 国产主播在线观看一区二区| 香蕉丝袜av| 色精品久久人妻99蜜桃| 国产成人影院久久av| 国产欧美日韩综合在线一区二区| 精品日产1卡2卡| www.精华液| 成人18禁在线播放| 亚洲色图综合在线观看| 国产乱人伦免费视频| 日韩精品青青久久久久久| 国产精品,欧美在线| 国产91精品成人一区二区三区| 午夜福利一区二区在线看| 国产高清视频在线播放一区| av天堂在线播放| 免费av毛片视频| 国产一区二区三区在线臀色熟女| 这个男人来自地球电影免费观看| 日本vs欧美在线观看视频| 成人国语在线视频| 国产av一区在线观看免费| 国产精品野战在线观看| 中文字幕av电影在线播放| 国产精品亚洲美女久久久| 动漫黄色视频在线观看| 午夜福利高清视频| 自线自在国产av| av在线播放免费不卡| 男女下面进入的视频免费午夜 | 久久国产亚洲av麻豆专区| 亚洲av五月六月丁香网| 成人永久免费在线观看视频| 搡老熟女国产l中国老女人| 成人永久免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成国产人片在线观看| 成人国语在线视频| 午夜福利成人在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品国产一区二区精华液| 多毛熟女@视频| 黑人欧美特级aaaaaa片| 精品国内亚洲2022精品成人| 成人精品一区二区免费| 精品国产国语对白av| 999久久久精品免费观看国产| 99re在线观看精品视频| 免费高清在线观看日韩| 亚洲人成电影观看| 亚洲欧美一区二区三区黑人| 色精品久久人妻99蜜桃| 嫁个100分男人电影在线观看| 免费不卡黄色视频| 99精品在免费线老司机午夜| 欧美成人午夜精品| 在线观看免费视频网站a站| 亚洲精品粉嫩美女一区| 性色av乱码一区二区三区2| 麻豆国产av国片精品| 午夜久久久久精精品| 国产精品国产高清国产av| 久久人人精品亚洲av| 久久久水蜜桃国产精品网| 一个人免费在线观看的高清视频| 88av欧美| 国内毛片毛片毛片毛片毛片| 18禁黄网站禁片午夜丰满| 怎么达到女性高潮| 一区二区三区高清视频在线| 亚洲视频免费观看视频| 亚洲精品中文字幕一二三四区| 欧美成人性av电影在线观看| 亚洲国产欧美一区二区综合| 99在线视频只有这里精品首页| 精品无人区乱码1区二区| 国产精品久久电影中文字幕| 窝窝影院91人妻| 两性午夜刺激爽爽歪歪视频在线观看 | 精品福利观看| 亚洲精品在线观看二区| 午夜福利,免费看| 一级,二级,三级黄色视频| 香蕉国产在线看| 99久久综合精品五月天人人| 久久青草综合色| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 岛国视频午夜一区免费看| 免费看美女性在线毛片视频| 精品一区二区三区四区五区乱码| 人妻丰满熟妇av一区二区三区| 一级a爱视频在线免费观看| 国产精品久久久av美女十八| 纯流量卡能插随身wifi吗| 日韩免费av在线播放| 成人av一区二区三区在线看| 婷婷精品国产亚洲av在线| 国产真人三级小视频在线观看| 久久久精品欧美日韩精品| 精品熟女少妇八av免费久了| 国产精品一区二区免费欧美| 老司机在亚洲福利影院| 日韩 欧美 亚洲 中文字幕| 欧美黄色片欧美黄色片| 亚洲性夜色夜夜综合| 十八禁人妻一区二区| 最近最新中文字幕大全电影3 | 欧美成人一区二区免费高清观看 | 午夜老司机福利片| 女人精品久久久久毛片| 国产精品野战在线观看| 国产一级毛片七仙女欲春2 | 日韩中文字幕欧美一区二区| 亚洲专区中文字幕在线| 一区福利在线观看| 国产欧美日韩精品亚洲av| 免费女性裸体啪啪无遮挡网站| 午夜福利在线观看吧| 一进一出抽搐动态| 亚洲精品国产一区二区精华液| 啦啦啦观看免费观看视频高清 | 亚洲精品国产精品久久久不卡| 国产高清有码在线观看视频 | 国产av一区在线观看免费| 看免费av毛片| 国产精品99久久99久久久不卡| 丝袜在线中文字幕| 自线自在国产av| 亚洲美女黄片视频| 国产视频一区二区在线看| АⅤ资源中文在线天堂| 满18在线观看网站| 欧美成人免费av一区二区三区| 又黄又爽又免费观看的视频| 男女午夜视频在线观看| 电影成人av| 精品国产乱子伦一区二区三区| 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产精品综合久久久久久久免费 | 国产91精品成人一区二区三区| 母亲3免费完整高清在线观看| 国产亚洲欧美精品永久| 色播在线永久视频| 身体一侧抽搐| 天堂动漫精品| 亚洲精品美女久久久久99蜜臀| 午夜福利视频1000在线观看 | av在线播放免费不卡| av网站免费在线观看视频| 欧美日韩亚洲国产一区二区在线观看| 老司机在亚洲福利影院| 免费少妇av软件| 中国美女看黄片| xxx96com| 好男人电影高清在线观看| 俄罗斯特黄特色一大片| 成年女人毛片免费观看观看9| 激情在线观看视频在线高清| 99re在线观看精品视频| 亚洲精华国产精华精| 午夜免费鲁丝| 亚洲全国av大片| 精品一品国产午夜福利视频| 69av精品久久久久久| 精品久久久久久成人av| 国语自产精品视频在线第100页| 日韩中文字幕欧美一区二区| 久久婷婷人人爽人人干人人爱 | 少妇熟女aⅴ在线视频| 一a级毛片在线观看| 夜夜爽天天搞| 九色国产91popny在线| 在线观看66精品国产| 国产精品影院久久| or卡值多少钱| 十分钟在线观看高清视频www| 91在线观看av| www.自偷自拍.com| 亚洲人成77777在线视频| av片东京热男人的天堂| 又黄又粗又硬又大视频| 国产精品香港三级国产av潘金莲| 久久精品成人免费网站| 亚洲色图综合在线观看| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 久99久视频精品免费| 黄色毛片三级朝国网站| av中文乱码字幕在线| 51午夜福利影视在线观看| 亚洲国产中文字幕在线视频| 亚洲五月色婷婷综合| 长腿黑丝高跟| 一级作爱视频免费观看| 热re99久久国产66热| 最好的美女福利视频网| 日韩欧美免费精品| 亚洲avbb在线观看| 免费在线观看视频国产中文字幕亚洲| 精品人妻1区二区| 国产成人精品在线电影| 日本一区二区免费在线视频| 桃色一区二区三区在线观看| 欧美国产日韩亚洲一区| av免费在线观看网站| 国产精品久久久人人做人人爽| 黄片大片在线免费观看| 国产精品一区二区免费欧美| 此物有八面人人有两片| 国产精品精品国产色婷婷| 国产精品秋霞免费鲁丝片| 脱女人内裤的视频| 999久久久精品免费观看国产| 国产精品免费一区二区三区在线| 丝袜人妻中文字幕| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 久热这里只有精品99| 亚洲欧美一区二区三区黑人| 在线av久久热| 欧美乱色亚洲激情| 搡老岳熟女国产| 亚洲一码二码三码区别大吗| 免费久久久久久久精品成人欧美视频| 九色国产91popny在线| 伊人久久大香线蕉亚洲五| 精品免费久久久久久久清纯| 夜夜躁狠狠躁天天躁| 国产av在哪里看| 狠狠狠狠99中文字幕| 可以免费在线观看a视频的电影网站| 在线观看舔阴道视频| 亚洲国产精品999在线| cao死你这个sao货| 久久人人爽av亚洲精品天堂| 亚洲少妇的诱惑av| 色播亚洲综合网| 脱女人内裤的视频| 久久国产精品男人的天堂亚洲| 久久性视频一级片| 91精品国产国语对白视频| 大型黄色视频在线免费观看| 极品教师在线免费播放| 午夜成年电影在线免费观看| 999久久久国产精品视频| 18禁美女被吸乳视频| 久久热在线av| 久久国产乱子伦精品免费另类| a在线观看视频网站| 成人精品一区二区免费| 女人高潮潮喷娇喘18禁视频| 成人亚洲精品一区在线观看| 精品欧美一区二区三区在线| 中文亚洲av片在线观看爽| 色综合亚洲欧美另类图片| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 午夜影院日韩av| 在线观看午夜福利视频| 看免费av毛片| 国产主播在线观看一区二区| 女警被强在线播放| 9热在线视频观看99| 欧美激情极品国产一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 欧美激情 高清一区二区三区| 久久精品国产99精品国产亚洲性色 | 精品日产1卡2卡| 在线观看66精品国产| 十八禁网站免费在线| aaaaa片日本免费| 又紧又爽又黄一区二区| 在线天堂中文资源库| 中文字幕色久视频| 午夜两性在线视频| 亚洲片人在线观看| 电影成人av| 我的亚洲天堂| 99久久久亚洲精品蜜臀av| 亚洲性夜色夜夜综合| 精品第一国产精品| 午夜影院日韩av| 十八禁人妻一区二区| netflix在线观看网站| 国产不卡一卡二| 日韩高清综合在线| 成人精品一区二区免费| 久久国产乱子伦精品免费另类| 夜夜爽天天搞| 国产精品 国内视频| 脱女人内裤的视频| 亚洲片人在线观看| 老司机午夜十八禁免费视频| 一级a爱视频在线免费观看| 免费在线观看完整版高清| 国产熟女午夜一区二区三区| 久久久久久久久中文| 免费少妇av软件| 久久久久九九精品影院| 日日爽夜夜爽网站| 日韩成人在线观看一区二区三区| 日本 欧美在线| 国产99久久九九免费精品| av免费在线观看网站| 国产精品久久久久久精品电影 | 后天国语完整版免费观看| 亚洲午夜精品一区,二区,三区| videosex国产| 黑人欧美特级aaaaaa片| 国产精品久久久av美女十八| 99riav亚洲国产免费| 日日爽夜夜爽网站| 最好的美女福利视频网| 青草久久国产| 欧美激情久久久久久爽电影 | 男男h啪啪无遮挡| 亚洲,欧美精品.| 久久草成人影院| 人人妻,人人澡人人爽秒播| 两性夫妻黄色片| 国产精品1区2区在线观看.| 18禁美女被吸乳视频| 亚洲欧美日韩高清在线视频| 中文字幕人成人乱码亚洲影| 亚洲成人国产一区在线观看| 欧美日韩一级在线毛片| 九色亚洲精品在线播放| 久久人妻福利社区极品人妻图片| 啦啦啦韩国在线观看视频| 1024香蕉在线观看| 亚洲国产精品久久男人天堂| 午夜免费成人在线视频| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 国产一区在线观看成人免费| 欧美成狂野欧美在线观看| 亚洲国产欧美网| 日韩欧美三级三区| 免费久久久久久久精品成人欧美视频| 91成人精品电影| 国产麻豆成人av免费视频| 亚洲成人免费电影在线观看| 免费搜索国产男女视频| 国产成+人综合+亚洲专区| 中文字幕人妻丝袜一区二区| 成熟少妇高潮喷水视频| 久久久久久人人人人人| 制服诱惑二区| 亚洲五月色婷婷综合| 色播亚洲综合网| 久久久久久久久中文| 黄片小视频在线播放| 大码成人一级视频| 人人妻人人澡欧美一区二区 | 国产精品免费视频内射| а√天堂www在线а√下载| 91在线观看av| 国产精品日韩av在线免费观看 | 国产真人三级小视频在线观看| 一区在线观看完整版| 欧美精品啪啪一区二区三区| 美女高潮到喷水免费观看| 久久久久久久久免费视频了| 久久精品aⅴ一区二区三区四区| 久久精品国产99精品国产亚洲性色 | 色婷婷久久久亚洲欧美| 欧美午夜高清在线| 国产国语露脸激情在线看| 中文亚洲av片在线观看爽| 少妇裸体淫交视频免费看高清 | 久久精品91蜜桃| 99国产精品一区二区蜜桃av| 午夜日韩欧美国产| 99精品欧美一区二区三区四区| 色尼玛亚洲综合影院| 国产不卡一卡二| 大码成人一级视频| www.自偷自拍.com| 国产欧美日韩综合在线一区二区| 欧美精品亚洲一区二区| 亚洲黑人精品在线| 欧美性长视频在线观看| 99热只有精品国产| 久久伊人香网站| 色在线成人网| 日韩有码中文字幕| 一区二区三区激情视频| 99国产精品一区二区蜜桃av| 亚洲午夜理论影院| 免费高清视频大片| 亚洲欧洲精品一区二区精品久久久| av天堂在线播放| 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 亚洲狠狠婷婷综合久久图片| 热re99久久国产66热| 国产亚洲欧美在线一区二区| 欧美一级a爱片免费观看看 | 美女 人体艺术 gogo| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 三级毛片av免费| 91字幕亚洲| 日韩有码中文字幕| 自线自在国产av| 亚洲成a人片在线一区二区| 精品少妇一区二区三区视频日本电影| 国产精品日韩av在线免费观看 | 婷婷六月久久综合丁香| 国产成人av教育| 久久香蕉精品热| 色精品久久人妻99蜜桃| 一边摸一边抽搐一进一出视频| 18禁美女被吸乳视频| www日本在线高清视频| 免费搜索国产男女视频| 黑人巨大精品欧美一区二区mp4| 在线观看午夜福利视频| 国产亚洲av高清不卡| 丝袜人妻中文字幕| 精品不卡国产一区二区三区| 69av精品久久久久久| 99国产精品一区二区三区| 日本免费a在线| 91九色精品人成在线观看| 岛国视频午夜一区免费看| www.熟女人妻精品国产| avwww免费| 婷婷丁香在线五月| 久久狼人影院| www.熟女人妻精品国产| 日韩欧美三级三区| 啦啦啦韩国在线观看视频| 国产亚洲欧美在线一区二区| av电影中文网址| av视频在线观看入口| 亚洲av五月六月丁香网| 麻豆av在线久日| 国产99白浆流出| 欧美色视频一区免费| 久久婷婷成人综合色麻豆| 给我免费播放毛片高清在线观看| 亚洲熟女毛片儿| 久久午夜亚洲精品久久| 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 一进一出抽搐gif免费好疼| 香蕉国产在线看| 亚洲中文字幕一区二区三区有码在线看 | 久久亚洲精品不卡| 国产av一区在线观看免费| 国产精品免费一区二区三区在线| 9色porny在线观看| 欧美黄色片欧美黄色片| 黄色a级毛片大全视频| 人人妻人人爽人人添夜夜欢视频| 国产av一区在线观看免费| 久久中文看片网| 免费在线观看日本一区| 色综合站精品国产| 91在线观看av| 日本免费a在线| 亚洲欧美精品综合一区二区三区| 青草久久国产| 一进一出抽搐gif免费好疼| 精品熟女少妇八av免费久了| 午夜福利高清视频| 无人区码免费观看不卡| 国产男靠女视频免费网站|