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

    基于葉片電信號邊際譜熵的玉米耐鹽堿性無損評價方法

    2018-02-28 06:31:53趙燕燕楊運經(jīng)杜光源
    農(nóng)業(yè)工程學(xué)報 2018年2期
    關(guān)鍵詞:鹽堿邊際電信號

    劉 鍇,趙燕燕,習(xí) 崗,楊運經(jīng),杜光源

    ?

    基于葉片電信號邊際譜熵的玉米耐鹽堿性無損評價方法

    劉 鍇1,趙燕燕2,習(xí) 崗1※,楊運經(jīng)3,杜光源3

    (1. 西安理工大學(xué)理學(xué)院,西安 710054; 2. 鄭州工業(yè)應(yīng)用技術(shù)學(xué)院基礎(chǔ)教學(xué)部,鄭州 451100; 3. 西北農(nóng)林科技大學(xué)理學(xué)院,楊凌 712100)

    為了探索能夠早期、靈敏、在位和無損檢測與評價植物耐鹽堿性的方法,將NaCl、Na2SO4、NaHCO3和Na2CO3配置成復(fù)合鹽堿溶液對耐鹽堿性較弱的玉米品種鄭單958和耐鹽堿性較強的玉米品種名玉20的玉米幼苗進行鹽堿脅迫,采集了鹽堿脅迫過程中鄭單958和名玉20幼苗葉片電信號,應(yīng)用HHT(Hilbert-Huang transformation)方法得到了2種玉米葉片電信號的邊際譜,分析了鹽堿脅迫過程中2個玉米品種葉片電信號邊際譜熵MSE(marginal spectrum entropy)變化的差異和生物學(xué)意義。結(jié)果顯示:鹽堿脅迫過程中,鄭單958葉片電信號的MSE表現(xiàn)出不斷下降的趨勢,葉片中丙二醛MDA(malondialdehyde)含量迅速升高;名玉20的MSE表現(xiàn)出波動性的變化,MDA含量變化不大,表明鄭單958葉片細(xì)胞的離子跨膜運輸被抑制,名玉20的葉片細(xì)胞存在著復(fù)雜的代謝調(diào)節(jié),鹽堿脅迫造成的葉片細(xì)胞膜脂過氧化可能是葉片電信號MSE變化的原因。研究發(fā)現(xiàn),鹽堿脅迫下耐鹽堿性不同的2個玉米品種的葉片電信號響應(yīng)指數(shù)RI(response index)差異明顯,在脅迫2、3和4 d時鄭單958的RI值分別比名玉20增長了42%、193%和332%。根據(jù)RI值的大小有可能對鹽堿脅迫下玉米葉片細(xì)胞離子運輸和細(xì)胞膜傷害的影響程度進行靈敏和早期的定量診斷,進而實現(xiàn)對玉米幼苗期耐鹽堿性強弱的在位和無損傷的評價。

    脅迫;作物;無損檢測;玉米;鹽堿脅迫;電信號;邊際譜熵;評價

    0 引 言

    目前,全球范圍的土壤鹽堿化不斷擴大[1]。土壤鹽堿化對植物會形成鹽脅迫導(dǎo)致滲透脅迫和離子毒害,還會形成堿脅迫下pH值升高造成營養(yǎng)脅迫[2-4],嚴(yán)重影響植物生長[5]。挖掘植物的耐鹽堿能力, 開發(fā)和選育耐鹽堿植物品種正成為研究的熱點[6]。

    要開發(fā)和選育耐鹽堿性植物,必須首先建立能夠?qū)χ参锬望}堿性強弱進行準(zhǔn)確判斷的評價方法。在現(xiàn)有研究中,植物耐鹽堿性評價采用的是基于多個形態(tài)指標(biāo)和生理生化指標(biāo)的隸屬函數(shù)法、聚類分析法和灰色關(guān)聯(lián)度分析法等綜合數(shù)量分析方法[7-11]。這些評價方法需要測量的指標(biāo)多,樣品量大,周期長,不能進行早期診斷,而且許多指標(biāo)必須通過試管試驗的破壞性測量來獲取,無法做到無損檢測。

    由于植物抗逆性是由多因素控制的,在細(xì)胞層面上體現(xiàn)的綜合性狀,是植物在逆境下通過代謝調(diào)節(jié)維持細(xì)胞功能狀態(tài)的能力[12-13],植物耐鹽堿性應(yīng)該基于鹽堿脅迫下植物細(xì)胞發(fā)出的生命信息,根據(jù)細(xì)胞代謝調(diào)節(jié)和功能狀態(tài)的變化程度來評價,其評價指標(biāo)必須能夠?qū)崿F(xiàn)活細(xì)胞或組織的無損和在位檢測,此外,還要求做到反應(yīng)靈敏,可以進行早期診斷。顯然,現(xiàn)有評價方法難以滿足這些要求。

    植物電信號是來自活細(xì)胞的生命信息,大量研究表明,植物電信號與化學(xué)信號、力學(xué)信號和水信號等構(gòu)成了植物的信息系統(tǒng)[14]。植物電信號一方面在組織和器官中傳遞信息流,另一方面直接參與代謝控制[15-19],調(diào)節(jié)呼吸代謝、光合作用、水分吸收和氣孔導(dǎo)度變化等核心生理過程,維持細(xì)胞的功能狀態(tài)[20-22]。分析和解讀植物電信號能夠?qū)χ参锶~片細(xì)胞代謝調(diào)節(jié)和功能狀態(tài)進行定量診斷和評價,并做到實時、在位和無損檢測[22-25]。因此,研究鹽堿脅迫下植物葉片電信號的變化規(guī)律及其生物學(xué)意義有可能為植物耐鹽堿性的無損檢測與評價提供一種新方法,然而有關(guān)研究未見報道。

    植物葉片電信號邊際譜是基于HHT(Hilbert-Huang transformation)方法獲得的信號在頻域上的分布[26],邊際譜熵MSE(marginal spectrum entropy)是邊際譜分布混亂程度的量度,其綜合反映了葉片細(xì)胞群跨膜電位和細(xì)胞之間的電偶聯(lián)狀況[27]。玉米是全球主要的糧食作物,其對鹽堿脅迫十分敏感[7,28],研究玉米耐鹽堿性評價方法具有重要的現(xiàn)實意義。鑒于此,本文對比研究了鹽堿脅迫下2個玉米品種葉片電信號邊際譜MSE變化的差異及其生物學(xué)意義,探討了基于葉片電信號邊際譜MSE評價植物耐鹽堿性的可行性,為開發(fā)能夠早期、靈敏、快速和無損傷在位檢測和評價植物耐鹽堿性的技術(shù)提供參考。

    1 材料與方法

    1.1 材料選擇與培養(yǎng)

    玉米品種鄭單958和名玉20為種子市場購得。試驗前選取2個品種外觀一致、顆粒飽滿的玉米種子各200粒,用蒸餾水洗掉外層紅色包衣,用質(zhì)量分?jǐn)?shù)為0.2%的HgCl2溶液消毒并用大量蒸餾水洗滌,隨后放置在溫度35 ℃,濕度50%和光照3 000 lx的PRX-1000A型人工培養(yǎng)箱中進行催芽,待種子萌發(fā)后選取發(fā)芽一致的種子進行培養(yǎng)。

    1.2 鹽堿處理

    由于現(xiàn)有玉米耐鹽堿性評價研究一般采用的脅迫方式是NaCl脅迫或NaCl和Na2CO3混合脅迫[7],脅迫方式較為簡單。而實際鹽堿地多是中性鹽(NaCl、Na2SO4)和堿性鹽(Na2CO3、NaHCO3)的復(fù)合鹽堿地,其pH值較高[8],這種復(fù)合鹽堿脅迫下玉米的研究報道較少?;诖耍疚膮⒖嘉墨I[29]的方法,將2種中性鹽NaCl、Na2SO4和2種堿性鹽NaHCO3、Na2CO3按照1:9:9:1摩爾質(zhì)量比混合后加入蒸餾水,定容成pH值為9.09±0.05的重度復(fù)合鹽堿溶液。在試驗中證實,在該復(fù)合鹽堿溶液處理下耐鹽堿性不同的玉米可以被有效地區(qū)分開來。因此,本研究采用此種復(fù)合鹽堿溶液對玉米幼苗進行鹽堿處理。試驗時待2個玉米品種的幼苗長至3片真葉(株高約15 cm)時,將培養(yǎng)皿中的蒸餾水(pH值為7.08)換為復(fù)合鹽堿溶液對玉米進行鹽堿脅迫處理。

    1.3 玉米葉片干質(zhì)量的測量

    在鹽堿脅迫的不同時間,每個品種各取15株長勢均勻的玉米幼苗,分為3組,每組5株,剪下全部葉片部分,用濾紙清理和擦干玉米葉片表面雜質(zhì)和水分,將其放入樣品杯子中置于烘箱內(nèi)以105 ℃殺青15 min,80 ℃烘干至質(zhì)量恒定。然后取出樣品,冷卻至室溫,用電子天平分別測量葉片干質(zhì)量,取每個品種3組葉片干質(zhì)量的平均值。

    1.4 玉米葉片電信號的采集

    參照文獻[23,27]的方法在鹽堿脅迫開始后每天上午10:00測量葉片干質(zhì)量的同時,分別采集2個玉米品種另外各3株玉米相同葉位的葉片電信號。采集儀器為BL420S生物機能試驗系統(tǒng)(成都泰盟科技有限公司生產(chǎn)),該系統(tǒng)具備高輸入阻抗(>1010Ω)、高共模抑制比(>120 dB)和低噪聲(<1V)、低漂移等特點。研究已經(jīng)證明,該系統(tǒng)能夠滿足采集植物葉片電信號的要求[23,27-30]。測試電極為丹麥Ambu公司生產(chǎn)的P-00-S型醫(yī)用高靈敏度Ag/AgCl心電電極。試驗之前對電極進行的測試表明,溫度在20~45 ℃范圍內(nèi)變化時,電極采集到的0~50 Hz頻段內(nèi)的信號振幅波動幅度小于0.5%,表明電極的穩(wěn)定性好[27]。

    采集葉片電信號時,1片采集電極貼于玉米葉片正面靠近葉尖位置,另1片電極貼于玉米葉片反面根部位置,電極間距15 cm,參考電極通過導(dǎo)線直接連接至培養(yǎng)槽中放置的銅片上。采集電極通過導(dǎo)電膠與葉片保持良好接觸,系統(tǒng)設(shè)置的采樣頻率為2 kHz,開啟50 Hz工頻抑制[23,27]。同時將被采集的玉米植株和采集系統(tǒng)放置在Faraday籠內(nèi)屏蔽環(huán)境電磁干擾。此外,植物生長和每次采集電信號時的環(huán)境條件相同,環(huán)境溫度為20 ℃,濕度為50%,光照強度為3 000 lx。為了減少對葉片細(xì)胞活動的影響,每次信號采集完畢后,及時拆去采集電極。

    1.5 玉米葉片電信號邊際譜MSE的計算

    研究發(fā)現(xiàn),當(dāng)采集到的玉米葉片電信號時域波形的數(shù)據(jù)時長在200 s以上時,信號的復(fù)雜度是穩(wěn)定的[27],因此,本文在采集到的葉片電信號時域波形的數(shù)據(jù)中,取數(shù)據(jù)時長為200 s的信號進行HHT變換,得到鹽堿脅迫過程中2個玉米品種各3個處理組葉片電信號的希爾伯特譜,再由希爾伯特譜得到邊際譜,然后分別計算各邊際譜的相對譜熵MSE,分別取2個品種各3個處理組的平均值進行動態(tài)分析。

    式中為()的序列長度。

    1.6 玉米葉片電信號相應(yīng)指數(shù)RI的計算

    在鹽堿脅迫下,玉米葉片電信號MSE將發(fā)生改變,其變化程度可以通過基于邊際譜熵的電信號響應(yīng)指數(shù)RI(response index)來表征

    式中MSE0為脅迫發(fā)生前的邊際譜熵,MSE為脅迫過程中的邊際譜熵。

    1.7 玉米葉片MDA的測量

    在上述采集玉米葉片電信號的同時,測量鹽堿脅迫過程中其他玉米植株相同葉位的葉片中MDA(Malondialdehyde)含量,MDA含量的測量方法見參考文獻[31]。將葉片先用濾紙吸干表面培養(yǎng)液,然后放入電子天平稱量鮮質(zhì)量并記錄。將樣品剪碎后,加入2 mL預(yù)冷的5%三氯乙酸,加入少量石英砂,在經(jīng)過冰浴的研缽內(nèi)研磨至勻漿,轉(zhuǎn)移到5 mL刻度離心試管,將研缽用緩沖液洗凈,清洗液移入離心管中,用5%三氯乙酸定容至5 mL,4 000 r/min下離心10 min,吸取上清液2 mL,加入0.67%硫代巴比妥酸的10%三氯乙酸溶液3 mL,于沸水浴中加熱30 min,迅速冷卻,再于4 000 r/min離心機中離心10 min,取上清液,以蒸餾水為空白調(diào)透光率,測定450、532、600 nm波長下的吸光度,按照下式計算MDA的濃度

    (mol/L)=6.45(532-600)-0.56450(4)

    再根據(jù)葉片鮮質(zhì)量計算測定樣品中MDA的含量

    每次試驗重復(fù)3次,取平均值。

    1.8 統(tǒng)計分析

    各測量結(jié)果均用平均值±標(biāo)準(zhǔn)差表示,差異顯著性采用SPSS軟件進行分析,<0.05表示差異到達顯著水平,<0.01表示差異達到極顯著水平。數(shù)值擬合采用origin軟件進行。

    2 結(jié)果與分析

    2.1 鹽堿脅迫對2個玉米品種幼苗葉片干質(zhì)量的影響

    葉片是植物進行物質(zhì)累積的主要場所,葉片干質(zhì)量是葉片物質(zhì)生產(chǎn)量的最終體現(xiàn)。葉片干質(zhì)量的變化是植物葉片對鹽堿脅迫的綜合反映,也是判斷植物耐鹽堿性強弱最可靠的指標(biāo)[7,32]。為了確定試驗的2個玉米品種耐鹽堿性的強弱,首先測量了鹽堿脅迫下鄭單958和名玉20玉米品種葉片干質(zhì)量變化的差異,結(jié)果見圖1。由圖1可見,在鹽堿脅迫過程中鄭單958品種葉片干質(zhì)量的增長緩慢,在脅迫4 d時干質(zhì)量反而有所下降;而名玉20品種的葉片干質(zhì)量則一直明顯增長,在脅迫4 d時干質(zhì)量才停止增長,此后保證基本不變的趨勢。在表觀上,在脅迫3 d時鄭單958品種的葉片開始卷曲,葉色變黃,名玉20品種的葉片形態(tài)和葉色變化不明顯。結(jié)果表明:鹽堿脅迫抑制了鄭單958品種葉片物質(zhì)的積累(<0.01),對名玉20品種物質(zhì)積累的影響相對較小。由此可見鄭單958品種的耐鹽堿性較差,這個結(jié)論與文獻報道是一致的[9]。

    圖1 鹽堿脅迫對玉米幼苗葉片干質(zhì)量的影響

    2.2 鹽堿脅迫下玉米葉片電信號的時域波形

    圖2和圖3分別為電極采集到的鹽堿脅迫下鄭單958和名玉20玉米各3株相同葉位的葉片電信號平均值的時域波形,其具有非平穩(wěn)和隨機的特征。在正常情況下(未脅迫時),振幅波動的幅度在1 000~2 000V以內(nèi),與以往報道的結(jié)果是一致的[23,27,30]。鹽堿脅迫對玉米葉片電信號有明顯影響,隨著脅迫時間的延長,葉片電信號的振幅呈現(xiàn)出減小的趨勢;在脅迫第5天時,鄭單958玉米葉片電信號的振幅降至20V以內(nèi),已經(jīng)被噪聲所湮沒;與此對應(yīng),名玉20品種還存在振幅在100V以內(nèi)的微弱電信號。圖2和圖3中測量開始后電信號出現(xiàn)的漸變過程,可能是采集電極貼附在葉片上時細(xì)胞的適應(yīng)過程。

    圖2 鹽堿脅迫下鄭單958玉米葉片電信號時域波形

    圖3 鹽堿脅迫下名玉20玉米葉片電信號的時域波形

    2.3 鹽堿脅迫下玉米葉片電信號MSE的變化

    電極采集到的玉米葉片電信號是一種微弱和非穩(wěn)定的隨機信號,HHT方法特別適合于這種信號的信息提取[27,33]?;贖HT方法分別計算出2個玉米品種葉片電信號的邊際譜如圖4和圖5所示。貌似無序的玉米葉片電信號的邊際譜具有連續(xù)的特征,經(jīng)歷不同時間的鹽堿脅迫,邊際譜的復(fù)雜度發(fā)生了變化。

    圖4 鹽堿脅迫下鄭單958玉米葉片電信號邊際譜

    圖5 鹽堿脅迫下名玉20玉米葉片電信號的邊際譜

    邊際譜的復(fù)雜度通過邊際譜熵MSE來描述,邊際譜越狹窄,MSE越小,表示信號中存在明顯的振蕩節(jié)律,復(fù)雜度??;反之,邊際譜越平坦,MSE越大,表明信號的復(fù)雜程度越高[27]。根據(jù)圖4和圖5計算得到的鹽堿脅迫下鄭單958和名玉20玉米葉片電信號邊際譜熵MSE的變化如圖6a所示。由圖6a可見,鄭單958玉米葉片正常生長時的MSE在0.53附近,在鹽堿脅迫過程中MSE出現(xiàn)了單調(diào)下降的趨勢,在脅迫3、4和5 d時,MSE分別下降了39.3%、48.9%和51.2%。與此相對應(yīng),在鹽堿脅迫過程中,名玉20品種的MSE呈現(xiàn)出波動的變化,在脅迫4 d時,MSE才開始迅速下降,說明此時名玉20品種葉片細(xì)胞膜的結(jié)構(gòu)才發(fā)生了不可逆的變化。

    圖6 鹽堿脅迫下2種玉米品種葉片電信號邊際譜熵MSE和丙二醛MDA含量的變化

    2.4 鹽堿脅迫下2個玉米品種葉片MDA含量的變化

    為了說明玉米葉片電信號邊際譜MSE變化的生物學(xué)意義,本文在鹽堿脅迫的不同時間采集玉米葉片電信號的同時,對比研究了同批次復(fù)合鹽堿脅迫處理的2個玉米品種相同葉位的葉片中MDA含量的變化,如圖6 b所示。研究MDA的原因在于MDA含量常作為植物細(xì)胞膜受損的主要指示物,通過MDA含量的變化可以了解細(xì)胞膜脂過氧化的程度,進而對膜系統(tǒng)受損程度進行評價[34-35]。在植物耐鹽堿性的研究中,MDA作為最重要的生化評價指標(biāo),被廣泛的采用[8,36]。由圖6b可見,在復(fù)合鹽堿脅迫2 d時,鄭單958品種的葉片MDA含量開始迅速上升,在脅迫3、4和5 d時MDA的相對增長率分別為72.7%、125.8%和143.8%,差異均達到了極顯著的水平(<0.01),MDA的這種變化趨勢與一些有關(guān)鹽脅迫的研究報道是類似的[37-39]。與此相對應(yīng),在脅迫1~4 d期間,名玉20品種的葉片中MDA含量變化不大,從4d開始才顯著升高,在脅迫5 d時相對增長率達到了114.0%,差異達到了極顯著的水平(<0.01)。

    比較圖6a和圖6b可知,鄭單958和名玉20品種的MSE和MDA之間可能具有一定的相關(guān)性,擬合結(jié)果如圖7所示。

    圖7 鹽堿脅迫下2種玉米品種葉片邊際譜熵MSE和丙二醛MDA的相關(guān)性

    2.5 鹽堿脅迫下2個玉米品種葉片電信號響應(yīng)指數(shù)RI的變化

    由于在相同的鹽堿脅迫過程中,耐鹽堿性不同的鄭單958和名玉20葉片電信號MSE的變化程度有所差異。為了定量表示MSE的變化程度,本文定義了基于MSE的葉片電信號響應(yīng)指數(shù)RI。RI越大表明葉片電信號MSE的變化越大,即葉片電信號邊際譜的復(fù)雜程度變化越大。鹽堿脅迫過程中鄭單958和名玉20品種的RI值如圖8所示。在相同的鹽堿脅迫條件下,耐鹽堿性較弱的鄭單958品種的RI值在脅迫開始后就呈現(xiàn)出單調(diào)上升的趨勢,在脅迫2 d后迅速增大,而耐鹽堿性較強的名玉20品種的RI值在脅迫第4 d以后才開始明顯增大,從脅迫2 d開始,鄭單958的RI值就一直明顯高于名玉20。在脅迫2、3和4 d時鄭單958的RI值分別比名玉20的RI值相對增長了42%、193%和332%,差異十分顯著。

    圖8 鹽堿脅迫下2個玉米品種電信號響應(yīng)指數(shù)RI的變化

    3 討 論

    植物電信號的來源非常復(fù)雜,包括系統(tǒng)電位(system potential,SP)、動作電位(action potential, AP)和變異電位(variation potential, VP),AP與細(xì)胞膜上Ca2+、Cl-和K+等離子通道的活性有關(guān),VP與細(xì)胞膜上H+泵的失活有關(guān),SP與H+泵的激活有關(guān)[40]。本文的研究發(fā)現(xiàn),在正常的生理狀態(tài)下(脅迫開始時),葉片電信號的MSE數(shù)值相對最大,表明葉片細(xì)胞生理電活動的復(fù)雜度較高;在鹽堿脅迫下,MSE表現(xiàn)出下降的趨勢。由于在鹽堿脅迫過程中,植物首先遭遇滲透脅迫和離子毒害[4,41],堿脅迫下pH 的升高會進一步破壞細(xì)胞的離子吸收[42],使根系細(xì)胞膜的H+泵失活和K+等離子通道活性變化,離子外流增強[43-44],其結(jié)果導(dǎo)致根系細(xì)胞膜上的SP發(fā)生變化、并誘發(fā)AP產(chǎn)生[45-47]。當(dāng)AP傳遞到葉片時葉片細(xì)胞膜上的H+泵將受到抑制,導(dǎo)致葉片電信號中SP信號的振幅和復(fù)雜度減少,影響葉片細(xì)胞的光合作用[22]、呼吸作用[48-49]和ATP含量產(chǎn)生[50]等重要的生理過程。因此,鹽堿脅迫下葉片細(xì)胞膜上離子跨膜運輸?shù)囊种瓶赡苁侨~片電信號MSE減小的原因。在鹽堿脅迫的持續(xù)進行下,葉片細(xì)胞膜會發(fā)生嚴(yán)重?fù)p傷,造成細(xì)胞膜上各種離子通道的破壞和離子泵的瓦解,大量電解質(zhì)外滲[51],結(jié)果使葉片電信號的節(jié)律性大為增強,MSE進一步減小。由此看來,鹽堿脅迫下植物葉片電信號MSE發(fā)生變化意味著葉片細(xì)胞膜上的離子運輸受到了影響,而MSE的快速下降則表明細(xì)胞膜受到了損傷。由于MSE的變化發(fā)生在細(xì)胞功能狀態(tài)變化之前,采集電極所涵蓋的任何一個細(xì)胞生理電活動的變化都會導(dǎo)致葉片電信號MSE的變化,因此,依據(jù)MSE的變化判斷葉片細(xì)胞功能狀態(tài)的變化具有早期和靈敏的特點;又由于貼片電極是將表面電極直接貼附在植物葉片上活體測量,不會對植物造成傷害,所以還具有在位和無損傷的優(yōu)點。

    本研究發(fā)現(xiàn),鹽堿脅迫下耐鹽堿性不同的鄭單958和名玉20品種葉片電信號MSE和丙二醛MDA的變化規(guī)律不同。由于2個玉米品種葉片中MDA含量的變化與MSE的變化呈現(xiàn)負(fù)相關(guān),說明鹽堿脅迫造成的葉片細(xì)胞膜脂過氧化可能是MSE下降的原因。

    基于鹽堿脅迫下耐鹽堿性不同的2個玉米品種MSE的變化程度不同,本文定義了電信號響應(yīng)指數(shù)RI,用于表征MSE的變化程度。RI值越大,表明葉片細(xì)胞受到的傷害越大,耐鹽堿性越差。因此,依據(jù)RI的大小可以對鹽堿脅迫下不同品種的植物葉片細(xì)胞傷害的影響程度進行定量診斷,進而對耐鹽堿性強弱做出評價。在本文中,相同的鹽堿脅迫下耐鹽堿性不同的2個玉米品種的RI值差異明顯,根據(jù)這種差異可以方便的確定出玉米耐鹽堿性的強弱。將RI作為評價指標(biāo)的優(yōu)點在于RI值始終處于0和1之間,與作物種類和品種無關(guān),是統(tǒng)一的量化標(biāo)準(zhǔn),易于計算機定量分析與自動評判。

    應(yīng)該指出,本文中電信號是在植物葉片上無損和在位采集的,MSE和RI的數(shù)值大小反映的是鹽堿脅迫下葉片細(xì)胞群離子運輸和細(xì)胞膜發(fā)生的變化,評價指標(biāo)RI值與物種無關(guān),因此,本文提出的評價植物耐鹽堿性的方法不但具有早期、靈敏和無損傷的特點,同時還具有普適性。

    4 結(jié) 論

    1)在pH值為9.09的復(fù)合鹽堿溶液脅迫過程中,耐鹽堿性較差的鄭單958品種葉片電信號MSE(marginal spectrum entropy)表現(xiàn)出不斷下降的趨勢,在脅迫3、4和5 d時,MSE分別下降了39.3%、48.9%和51.2%。耐鹽堿性較強的名玉20品種的MSE呈現(xiàn)出波動的變化,表明鹽堿脅迫下鄭單958品種葉片細(xì)胞的離子跨膜運輸被抑制,名玉20品種的葉片細(xì)胞內(nèi)存在復(fù)雜的代謝調(diào)節(jié),以維持葉片細(xì)胞離子運輸?shù)膭討B(tài)平衡和正常的功能狀態(tài)。

    2)在復(fù)合鹽堿脅迫2 d時,鄭單958品種的葉片MDA(malondialdehyde)含量迅速上升,在脅迫3、4和5 d時MDA的相對增長率分別為72.7%、125.8%和143.8%,名玉20品種的MDA含量在脅迫4 d后才開始明顯上升,表明鄭單958品種在脅迫初期就出現(xiàn)了膜脂過氧化現(xiàn)象,隨著脅迫過程的進行,膜脂過氧化現(xiàn)象越來越嚴(yán)重,名玉20品種在脅迫4 d以后才發(fā)生顯著的膜脂過氧化現(xiàn)象,鹽堿脅迫造成的葉片細(xì)胞膜脂過氧化可能是葉片電信號MSE下降的原因。

    3)從鹽堿脅迫2 d開始,鄭單958玉米品種葉片電信號響應(yīng)指數(shù)RI(response index)的數(shù)值一直顯著高于名玉20。在脅迫2、3和4 d時鄭單958的RI值分別比名玉20的RI值相對增長了42%、193%和332%,根據(jù)RI值的大小可以對鹽堿脅迫下玉米葉片細(xì)胞離子運輸和細(xì)胞膜傷害的影響程度進行靈敏和早期的定量診斷,進而實現(xiàn)對玉米幼苗耐鹽堿性強弱的在位和無損傷的評價。

    [1] Wang W X, Vinocur B, Altman A. Plant responses to drought, salinity and extreme temperatures: Towards genetic engineering for stress tolerance[J]. Planta, 2003, 218(1): 1-14.

    [2] Yang J Y, Zhang W, Tian Y, et al. Effects of various mixed salt-alkaline stresses on growth, photosynthesis, and photosynthetic pigment concentrations ofseedlings[J]. Photosynthetica, 2011, 49(2): 275-284.

    [3] Gong B, Wang X F, Wei M, et al. Overexpression of S-adenosylmethionine synthetase 1 enhances tomato callus tolerance to alkali stress through polyamine and hydrogen peroxide cross-linked networks[J]. Plant Cell, Tissue and Organ Culture, 2016, 124(2): 377-391.

    [4] Wang X P, Jiang P, Ma Y, et al. Physiological strategies of sunflower exposed to salt or alkali stresses: Restriction of ion transport in the cotyledon node zone and solute accumulation[J]. Agronomy Journal, 2015, 107(6): 2181-2192.

    [5] Wang J C, Yao L R, Li B C, et al. Comparative proteomic analysis of cultured suspension cells of the halophyteby iTRAQ provides insights into response mechanisms to salt stress[J]. Frontiers in Plant Science, 2016, 7(30): 1-12.

    [6] 王佺珍,劉倩,高婭妮,等. 植物對鹽堿脅迫的響應(yīng)機制研究進展[J]. 生態(tài)學(xué)報,2017,37(16):5565-5577.

    Wang Quanzhen, Liu Qian, Gao Yani, et al. Review on the mechanism of the response to salinity-alkalinity stress in plants[J]. Acta Ecologica Sinica, 2017, 37(16): 5565-5577. (in Chinese with English abstract)

    [7] 張春宵,袁英,劉文國,等. 玉米雜交種苗期耐鹽堿篩選與大田鑒定的比較分析[J]. 玉米科學(xué),2010,18(5):14-18.

    Zhang Chunxiao, Yuan Ying, Liu Wenguo, et al. Comparative analysis between salt-alkali tolerance in seedling stage and production in the field of maize hybrid[J]. Corn Science, 2010, 18(5): 14-18. (in Chinese with English abstract)

    [8] 趙俊香,任翠梅,吳鳳芝,等. 16份菊芋種質(zhì)苗期耐鹽堿性篩選與綜合鑒定[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報,2015,23(5):620-627.

    Zhao Junxiang, Ren Cuimei, Wu Fengzhi, et al. Comprehensive identification of saline-alkaline tolerance of 16 Jerusalem artichoke accessions at seedling stage[J]. Chinese Journal of Eco-Agriculture, 2015, 23(5): 620-627. (in Chinese with English abstract)

    [9] 孫浩,張保望,李宗新,等. 夏玉米品種鹽堿脅迫耐受能力評價[J]. 玉米科學(xué),2016,24(1):81-87.

    Sun Hao, Zhang Baowang, Li Zongxin, et al. Different tolerance ability of maize varieties to saline-alkali soil conditions [J]. Journal of Maize Sciences, 2016, 24(1): 81-87. (in Chinese with English abstract)

    [10] 張婷婷,于崧,于立河,等. 松嫩平原春小麥耐鹽堿性鑒定及品種(系)篩選[J]. 麥類作物學(xué)報,2016,36(8):1008-1019.

    Zhang Tingting,Yu Song,Yu Lihe,et al. Saltine-alkaline tolerance identification and varieties (Line) screening of spring wheat in songnen plain[J]. Journal of Triticeae Crops, 2016,36(8):1008?1019. (in Chinese with English abstract)

    [11] 王晨,陳吉寶,龐振凌,等. 甜高粱對混合鹽堿脅迫的響應(yīng)及耐鹽堿種質(zhì)鑒定[J]. 作物雜志,2016(1):56-61. Wang Chen, Chen Jibao, Pang Zhenling, et al. The Response and screening of germplasm tolerant to mixed saline-alkali stress in sweet sorghum[J]. Crops, 2016(1): 56-61. (in Chinese with English abstract)

    [12] Zhu J K. Salt and drought stress signal transduction in plant[J]. Annual Review of Plant Biology, 2002, 53: 247-273.

    [13] 張燦軍,冀天會,楊子光,等. 小麥抗旱性鑒定方法及評價指標(biāo)研究Ⅰ,鑒定方法及評價指標(biāo)[J]. 中國農(nóng)學(xué)通報,2007,23(9):226-230.

    Zhang Canjun, Ji Tianhui, Yang Ziguang, et al. Study on resistance drought identify method and evaluation index of wheat Ⅰ, identify method and evaluation index[J]. Chinese Agricultural Science Bulletin, 2007, 23(9): 226-230. (in Chinese with English abstract)

    [14] Fromm J, Lautner S. Electrical signals and their physiological significance in plants[J]. Plant, Cell & Environment, 2007, 30(3): 249-257.

    [15] Fromm J, Fei H. Electrical signaling and gas exchange in maize plants of drying soil[J]. Plant Science, 1998, 132(2): 203-213.

    [16] Lautner S, Erhard T, Matyssek R, et al. Characteristics of electrical signals in poplar and responses in photosynthesis[J]. Plant Physiology, 2005, 138(4): 2200-2209.

    [17] Volkov A, Adesina T, Markin V, et al. Kinetics and mechanism of Dionaea muscipula trap closing[J]. Plant Physiology, 2008, 146(2): 694-702.

    [18] Gil P M, Gurovich L, Schaffer B, et al. Root to leaf electrical signaling in avocado in response to light and soil water content[J]. Journal of Plant Physiology, 2008, 165(10): 1070-1078.

    [19] Gurovich L,Hermosilla P.Electric signaling in fruit trees in response to water applications and light-darkness conditions[J]. Journal of Plant Physiology, 2009, 166(3): 290-300.

    [20] Koziolek C, Grams T E E, Schreiber U, et al. Transient knockout of photosynthesis mediated by electrical signals[J]. New Phytologist, 2004, 161(3): 715-722.

    [21] Masi E, Ciszak M, Stefano G, et al. Spatiotemporal dynamics of the electrical network activity in the root apex[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(10): 4048-4053.

    [22] Oyarce P, Gurovich L. Evidence for the transmission of information through electric potentials in injured avocado trees[J]. Journal of Plant Physiology, 2011, 168(2): 103-108.

    [23] Zhang X H, Yu N M, Xi G, et al. Changes in the power spectrum of electrical signals in maize leaf induced by osmotic stress[J]. Chinese Science Bulletin, 2012, 57(4): 413-420.

    [24] Gil P M, Saavedra J, Schaffer B, et al. Quantifying effects of irrigation and soil water content on electrical potentials in grapevines (Vitis vinifera) using multivariate statistical methods[J]. Scientia Horticulturae, 2014, 173(3): 71-78.

    [25] Gallé A, Lautner S, Flexas J, et al. Environmental stimuli and physiological responses: The current view on electrical signalling[J]. Environmental & Experimental Botany, 2015, 114: 15-21.

    [26] 董紅生,邱天爽,張愛華,等. 基于HHT邊際譜熵和能量譜熵的心率變異信號的分析方法[J]. 中國生物醫(yī)學(xué)工程學(xué)報,2010,29(3):336-344.

    Dong Hongsheng, Qiu Tianshuang, Zhang Aihua, et al. The analysis method of heart rate variability signal based on the HHT marginal spectrum entropy and energy spectrum entropy[J]. Chinese Journal of Biomedical Engineering, 2010, 29(3): 336-344. (in Chinese with English abstract)

    [27] 劉鍇,習(xí)崗,賀瑞瑞,等. 滲透脅迫下玉米葉片電位波動邊際譜的變化與意義[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(1):199-205.

    Liu Kai, Xi Gang, He Ruirui, et al. Changes and significance of marginal spectrum on maize leaves potential fluctuations under osmotic stress[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(1): 199-205. (in Chinese with English abstract)

    [28] 馮艷,江鵬,王亞峰,等. 玉米種質(zhì)資源親緣關(guān)系的分子標(biāo)記鑒定及其耐鹽性評價[J]. 中國農(nóng)學(xué)通報,2013,29(33):79-84.

    Feng Yan, Jiang Peng, Wang Yafeng, et al. Identification of genetic relationship of corn germplasms with molecular markers and evaluation of their salt tolerance[J]. Chinese Agricultural Science Bulletin, 2013, 29(33): 79-84. (in Chinese with English abstract)

    [29] 陳金元,陳學(xué)林,滿吉琳,等. 混合鹽堿脅迫對紅砂種子萌發(fā)的影響[J]. 西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2016,44(5):113-119. Chen Jinyuan, Chen Xuelin, Man Jilin, et al. Effects of mixed saline-alkali stress on germination of Reaumuria songarica seeds[J]. Journal of Northwest A&F University: Nat. Sci.Ed., 2016, 44(5): 113-119. (in Chinese with English abstract)

    [30] Zhang X H, Sun H Q, Ma H B, et al. The changes of power spectrum of electrical signal in aloe under several damage conditions[J]. International Journal of Control and Automation, 2016, 9(6): 161-172.

    [31] 張玉榮,周顯青,張勇. 儲存玉米膜脂過氧化與生理指標(biāo)的研究 [J]. 中國農(nóng)業(yè)科學(xué),2008,41(10):3410-3414.

    Zhang Yurong, Zhou Xianqing, Zhang Yong. Research on membrane lipid peroxidation and physiological parameters of storage maize[J]. Scientia Agricultura Sinica, 2008, 41(10): 3410-3414. (in Chinese with English abstract)

    [32] Gong B, Wen D, Vanden L K, et al. Comparative effects of NaCl and NaHCO3stress on photosynthetic parameters, nutrient metabolism, and the antioxidant system in tomato leaves[J]. Scientia Horticulturae, 2013, 157(3): 1-12.

    [33] 關(guān)海鷗,李偉凱,杜松懷,等. 基于Hilbert-Huang 變換的生物觸電電流檢測模型[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(14):202-209.

    Guan Haiou, Li Weikai, Du Songhuai, et al. Detection model of biological electric shock current based on Hilbert-Huang transform[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(14): 202-209. (in Chinese with English abstract)

    [34] Sorkheh K, Shiran B, Rouhi V, et al. Salt stress induction of some key antioxidant enzymes and metabolites in eight Iranian wild almond species[J]. Acta Physiologiae Plantarum, 2012, 34(1): 203-213.

    [35] 習(xí)崗,賀瑞瑞,劉鍇,等. 應(yīng)用超弱光子輻射評價菠菜葉片衰老方法可行性[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(17):268-275.

    Xi Gang, He Ruirui, Liu Kai, et al. Feasibility of evaluation method for spinach leaf senescence based on biological ultraweak photon emission[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(17): 268-275. (in Chinese with English abstract)

    [36] 劉濱碩,康春莉,王鑫,等. 羊草對鹽堿脅迫的生理生化響應(yīng)特征[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(23):166-173.

    Liu Binshuo, Kang Chunli, Wang Xin, et al. Physiological and biochemical response characteristics of Leymus chinensis to saline-alkali stress[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(23): 166-173. (in Chinese with English abstract)

    [37] Tavakkoli E, Fatehi F, Coventry S, et al. Additive effects of Na+and Cl-ions on barley growth under salinity stress[J]. Journal of Experimental Botany, 2011, 62(6): 2189-2203.

    [38] 閆永慶,劉興亮,王崑,等. 白刺對不同濃度混合鹽堿脅迫的生理響應(yīng)[J]. 植物生態(tài)學(xué)報,2010,34(10):1213-1219.

    Yan Yongqing, Liu Xingliang, Wang Kun, et al. Effect of complex saline-alkali stress on physiological parameters of Nitratia tangutorum[J]. Chinese Journal of Plant Ecology, 2010, 34(10): 1213-1219. (in Chinese with English abstract)

    [39] Hazman M, Hause B, Eiche E, et al. Increased tolerance to salt stress in Opda-deficient rice Allene Oxide Cyclase mutants is linked to an increased Ros-scavenging activity[J]. Journal of Experimental Botany, 2015, 66(11): 3339-3352.

    [40] Zimmermann M R, Maischak H, Mithofer A, et al. System potentials,a novel electrical long-distance apoplastic signal in plants,induced by wounding[J]. Plant Physiology, 2009, 149(3): 1593-1600.

    [41] Munns R. Comparative physiology of salt and water stress[J]. Plant, Cell & Environment, 2002, 25(2): 239-250.

    [42] Ruiz K B, Biondi S, Martinez E A, et al. Quinoa-a model crop for understanding salt-tolerance mechanisms in halophytes[J]. Plant Biosystems, 2016, 150(2): 357-371.

    [43] Chen Z, Newman I, Zhou M, et al. Screening plants for salt tolerance by measuring K+flux: A case study for barley[J]. Plant, Cell & Environment, 2005, 28(10): 1230-1246.

    [44] Chen Z, Pottosin I I, Cuin T A, et al. Root plasma membrane transporters controlling K+/Na+homeostasis in salt-stressed barky[J]. Plant Physiology, 2007, 145(4): 1714-1725.

    [45] Felle H H, Zimmermann M R. Systemic signaling in barley through action potentials[J]. Planta, 2007, 226(1): 203-214.

    [46] Favre P, Greppin H, Agosli H D. Accession-dependent action potentials in Arabidopsis [J]. Journal of Plant Physiology, 2011, 168(7): 653-660.

    [47] Maffei M E, Mithofer A, Boland W. Before gene express ion:Early events in Plant-insect interaction[J]. Trends in Plant Science, 2007, 12(7): 310-316.

    [48] Gro?kinsky K D, Svensgaard J, Christensen S, et al. Plant phenomics and the need for physiological phenotyping across scales to narrow the genotype-to-phenotype knowledge gap[J]. Journal of Experimental Botany, 2015, 66(18): 5429-5440.

    [49] Hedrich R, Salvador-Recatal V, Dreyer I. Electrical wiring and long-distance plant communication[J]. Trends in Plant Science, 2016, 21(5): 376-387.

    [50] Fiorani F, Schurr U. Future scenarios for plant phenotyping[J]. Annual Review of Plant Biology, 2013, 64(1): 267-291.

    [51] Dinneny J R. Traversing organizational scales in plant salt-stress responses[J]. Current Opinion in Plant Biology, 2015, 23: 70-75.

    Nondestructive evaluation method for saline-alkaline tolerance of maize based on marginal spectral entropy of electric signal in leaf

    Liu Kai1, Zhao Yanyan2, Xi Gang1※, Yang Yunjing3, Du Guangyuan3

    (1.,’710054,; 2.,,451100,; 3.,712100,)

    Saline-alkaline stress (SAS) is one of the major abiotic stresses affecting the growth of plants. It has been a severe problem that restricts plant production and even the development of the ecological environment. The improvement of plant saline-alkaline tolerance and selection of saline-alkaline tolerance plant varieties are becoming hot spots for research. To develop and select saline-alkaline tolerance plants, an evaluation method that can accurately judge the plant saline-alkaline tolerance must be first established. In the present study, the evaluation of saline-alkaline tolerance of plants is generally based on morphological indicators and physiological and biochemical indicators. These evaluation methods require a large number of samples and long cycle, and cannot be early diagnosed. Moreover, many of the indicators must be obtained through the destructive measurement of test-tube experiments, which are not nondestructive testing. Therefore, the traditional evaluation method has many disadvantages. In order to explore the methods of early, sensitive, in situ and nondestructive testing saline-alkaline tolerance of plants, a complex solution consisting of NaCl, Na2SO4, NaHCO3and Na2CO3with pH value of 9.09 was used to stress 2 kinds of maize variety seedlings of Zhengdan 958 with poor saline-alkaline tolerance and Mingyu 20 with strong saline-alkaline tolerance. Time-domain waveforms of leaf electrical signals of Zhengdan 958 and Mingyu 20 seedlings during saline-alkaline stress were collected. The marginal spectra of 2 kinds of maize leaf electrical signals were obtained by Hilbert-Huang transformation (HHT). The changes of marginal spectrum entropy (MSE) of maize leaf electrical signals and the biological significance were analyzed. The results showed that the MSE of leaf electrical signals from maize variety Zhengdan 958 continued to decline in the process of saline-alkaline stress, while the MSE of leaf electrical signals from maize variety Mingyu 20 changed in volatility. It was indicated that the ion transport of leaf cells in Zhengdan 958 was inhibited under saline-alkaline stress, and there was complex metabolic regulation in leaf cells of Mingyu 20 to maintain the dynamic balance of ion transport and normal functional status of leaf cell. The study also found that the malondialdehyde (MDA) content in leaf of Zhengdan 958 was increasing in the process of saline-alkaline stress, and the MDA content in leaves of Mingyu 20 began to increase significantly after 4 days of stress. This phenomenon suggested that there was membrane lipid peroxidation in leaves of Zhengdan 958 in the early stages of saline-alkaline stress, and it was more and more serious with the process of the stress, however, there was significant membrane lipid peroxidation in leaves of Mingyu 20 after 4 days of stress. The membrane lipid peroxidation of leaf cells caused by saline-alkaline stress could be the reason to the decrease of MSE about maize leaf electrical signals. Due to that the variation of the MSE from 2 maize varieties with different saline-alkaline tolerance under saline-alkaline stress was different, the response index (RI) of electrical signal based on the MSE was defined in this paper. The results showed that the RI values of maize varieties Zhengdan 958 and Mingyu 20 were obviously different in the processes of saline-alkaline stress. The influence of saline-alkaline stress on the ion transport and cell membrane injury of maize leaf cell could be sensitive and early quantitatively diagnosed according to the size of RI, and then to achieve in situ measurement and nondestructive evaluating saline-alkaline tolerance of maize seedlings. Since the RI based on plant electrical signals has not relationship with the species, the method proposed in this paper to evaluate the saline-alkaline tolerance of maize seedlings may also have a wide range of applicability. It is expected that the evaluation method about saline-alkaline tolerance of plant proposed in this paper can be verified through a large number of experiments.

    stresses; crops; nondestructive testing; maize; saline-alkaline stress; electrical signal; marginal spectrum entropy; evaluation

    10.11975/j.issn.1002-6819.2018.02.027

    Q64

    A

    1002-6819(2018)-02-0197-08

    2017-09-08

    2018-01-08

    國家自然科學(xué)基金資助項目(31471412,51277151);陜西省教育廳科學(xué)研究計劃項目(15JK1515)

    劉鍇,講師,博士,主要從事生物光電信息的分析與應(yīng)用研究。Email:leaukai@gmail.com

    習(xí) 崗,陜西楊凌人,教授,研究方向為生物光電信息的分析與應(yīng)用。Email:xig@xaut.edu.cn

    劉 鍇,趙燕燕,習(xí) 崗,楊運經(jīng),杜光源. 基于葉片電信號邊際譜熵的玉米耐鹽堿性無損評價方法[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(2):197-204. doi:10.11975/j.issn.1002-6819.2018.02.027 http://www.tcsae.org

    Liu Kai, Zhao Yanyan, Xi Gang, Yang Yunjing, Du Guangyuan. Nondestructive evaluation method for saline-alkaline tolerance of maize based on marginal spectral entropy of electric signal in leaf[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 197-204. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.02.027 http://www.tcsae.org

    猜你喜歡
    鹽堿邊際電信號
    隨身新配飾
    基于聯(lián)合聚類分析的單通道腹部心電信號的胎心率提取
    鹽堿荒灘的“底色之變”
    金橋(2021年11期)2021-11-20 06:37:14
    鹽堿地區(qū)南美白對蝦、青蝦混養(yǎng)技術(shù)
    基于Code Composer Studio3.3完成對心電信號的去噪
    科技傳播(2019年24期)2019-06-15 09:29:28
    追求騎行訓(xùn)練的邊際收益
    社會治理的邊際成本分析
    基于隨機森林的航天器電信號多分類識別方法
    鹽堿脅迫對三種西藏牧草種子萌發(fā)及幼苗影響
    西藏科技(2016年8期)2016-09-26 09:00:55
    鹽堿脅迫對亞麻苗期生長及陽離子吸收和分配的影響
    午夜a级毛片| 久久久国产成人免费| 免费在线观看成人毛片| 亚洲三级黄色毛片| 动漫黄色视频在线观看| 久久99热这里只有精品18| 国产一区二区三区在线臀色熟女| 日本色播在线视频| 欧美日韩乱码在线| 亚洲成a人片在线一区二区| 人妻丰满熟妇av一区二区三区| 亚洲三级黄色毛片| 国产免费男女视频| 午夜福利在线在线| av在线天堂中文字幕| 日韩欧美免费精品| or卡值多少钱| 亚洲综合色惰| 天堂动漫精品| 最近视频中文字幕2019在线8| 久久精品影院6| 乱人视频在线观看| 欧美激情在线99| 亚洲欧美清纯卡通| 国产精品一区www在线观看 | 乱码一卡2卡4卡精品| 啦啦啦观看免费观看视频高清| 麻豆一二三区av精品| 亚洲电影在线观看av| 久久婷婷人人爽人人干人人爱| 听说在线观看完整版免费高清| 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 婷婷精品国产亚洲av| www日本黄色视频网| 俄罗斯特黄特色一大片| 中文字幕熟女人妻在线| 精品久久久久久久久av| 国产精品免费一区二区三区在线| 亚洲精品成人久久久久久| 久久久久久久久久黄片| 综合色av麻豆| 超碰av人人做人人爽久久| 99热只有精品国产| 淫秽高清视频在线观看| 亚洲性夜色夜夜综合| 亚洲欧美日韩卡通动漫| 春色校园在线视频观看| 99热6这里只有精品| 成年女人毛片免费观看观看9| 桃色一区二区三区在线观看| 免费黄网站久久成人精品| 久久国产精品人妻蜜桃| 窝窝影院91人妻| 波多野结衣巨乳人妻| 黄色一级大片看看| 久久国产精品人妻蜜桃| 精品久久久噜噜| 在线国产一区二区在线| 免费无遮挡裸体视频| 老熟妇仑乱视频hdxx| 在线观看一区二区三区| 久久人人爽人人爽人人片va| 在线观看美女被高潮喷水网站| av视频在线观看入口| 午夜精品在线福利| 啦啦啦啦在线视频资源| 国产综合懂色| 国产亚洲精品久久久久久毛片| 五月伊人婷婷丁香| 一区二区三区免费毛片| 老女人水多毛片| 在线观看免费视频日本深夜| 精品久久久久久久久av| 韩国av在线不卡| 久久久久久久久中文| 国产精品久久久久久亚洲av鲁大| 成人特级av手机在线观看| 国产美女午夜福利| 人妻制服诱惑在线中文字幕| 亚洲av.av天堂| www.www免费av| 老女人水多毛片| 美女被艹到高潮喷水动态| 久久久久性生活片| 久久热精品热| 日本a在线网址| 国产精品福利在线免费观看| 亚洲乱码一区二区免费版| 免费观看在线日韩| 欧美日本亚洲视频在线播放| 亚州av有码| 午夜亚洲福利在线播放| 久久九九热精品免费| 男人狂女人下面高潮的视频| 久久久久久久亚洲中文字幕| 色5月婷婷丁香| 成年女人看的毛片在线观看| 国产一区二区三区av在线 | 久久草成人影院| 18+在线观看网站| 国产av在哪里看| 亚洲av第一区精品v没综合| 欧美成人免费av一区二区三区| netflix在线观看网站| 人人妻人人澡欧美一区二区| 亚洲七黄色美女视频| 蜜桃久久精品国产亚洲av| .国产精品久久| 国产精品一及| 久久久国产成人免费| 看黄色毛片网站| 自拍偷自拍亚洲精品老妇| 日本免费a在线| 欧美另类亚洲清纯唯美| 成年女人看的毛片在线观看| 国产一区二区亚洲精品在线观看| 国产欧美日韩一区二区精品| av专区在线播放| 此物有八面人人有两片| 69av精品久久久久久| 直男gayav资源| 动漫黄色视频在线观看| 琪琪午夜伦伦电影理论片6080| 成人av一区二区三区在线看| 天堂动漫精品| 夜夜夜夜夜久久久久| 亚洲va日本ⅴa欧美va伊人久久| 美女大奶头视频| 丝袜美腿在线中文| 精品久久久久久久久久免费视频| 国产高清视频在线播放一区| 日本免费a在线| 国内久久婷婷六月综合欲色啪| av福利片在线观看| 麻豆国产97在线/欧美| 成年免费大片在线观看| 成人特级黄色片久久久久久久| 国产91精品成人一区二区三区| 精品久久久久久久末码| 国产av麻豆久久久久久久| 久久欧美精品欧美久久欧美| 亚洲精品影视一区二区三区av| 蜜桃久久精品国产亚洲av| 成年女人看的毛片在线观看| 午夜亚洲福利在线播放| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 三级男女做爰猛烈吃奶摸视频| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| 亚洲欧美精品综合久久99| 波野结衣二区三区在线| 我要看日韩黄色一级片| 国产综合懂色| 亚洲欧美日韩无卡精品| 日本 欧美在线| 亚洲精品乱码久久久v下载方式| 日韩欧美在线二视频| 免费在线观看日本一区| 波多野结衣巨乳人妻| 午夜福利高清视频| 精品午夜福利在线看| 日本 欧美在线| 午夜免费激情av| 99热这里只有是精品在线观看| 窝窝影院91人妻| 久久久午夜欧美精品| 女人十人毛片免费观看3o分钟| 欧美日韩亚洲国产一区二区在线观看| 黄色配什么色好看| 美女高潮喷水抽搐中文字幕| 久久久国产成人精品二区| 亚洲色图av天堂| 在线免费观看的www视频| 免费在线观看成人毛片| 精品一区二区三区视频在线观看免费| 我的女老师完整版在线观看| 午夜激情福利司机影院| 欧美区成人在线视频| 精品久久国产蜜桃| 国产在视频线在精品| 久久久久久九九精品二区国产| 久久精品综合一区二区三区| 久久精品国产鲁丝片午夜精品 | 亚洲最大成人手机在线| 久久久久久久久久成人| 国产欧美日韩精品亚洲av| 麻豆精品久久久久久蜜桃| 少妇熟女aⅴ在线视频| 久久中文看片网| 麻豆精品久久久久久蜜桃| 在线观看66精品国产| 成年女人永久免费观看视频| 黄色丝袜av网址大全| 日本免费一区二区三区高清不卡| 国产伦人伦偷精品视频| 免费观看精品视频网站| 99久久精品热视频| 色视频www国产| 变态另类成人亚洲欧美熟女| 偷拍熟女少妇极品色| 精品久久久久久久久久久久久| 久久精品91蜜桃| 超碰av人人做人人爽久久| 最近最新中文字幕大全电影3| 日本-黄色视频高清免费观看| 国产一区二区亚洲精品在线观看| 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| 国产精品不卡视频一区二区| 国产精品一区二区三区四区久久| 99精品久久久久人妻精品| 人妻制服诱惑在线中文字幕| 变态另类丝袜制服| 亚洲第一电影网av| 日韩欧美精品v在线| 一本一本综合久久| 精品不卡国产一区二区三区| 在线观看美女被高潮喷水网站| 别揉我奶头 嗯啊视频| 日韩欧美国产一区二区入口| 久久午夜福利片| 国产在线男女| 亚洲av美国av| 成人二区视频| 伦精品一区二区三区| 三级毛片av免费| 亚洲人成伊人成综合网2020| 国产亚洲精品综合一区在线观看| 久久久久久久久久成人| 日本黄色视频三级网站网址| 搡女人真爽免费视频火全软件 | 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看| 99久久九九国产精品国产免费| 国产爱豆传媒在线观看| 亚洲精品色激情综合| 99久久精品国产国产毛片| 成人永久免费在线观看视频| 99在线人妻在线中文字幕| 亚洲人成网站在线播| 中文字幕av成人在线电影| 成人无遮挡网站| 亚洲成av人片在线播放无| 久久精品国产清高在天天线| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 人妻少妇偷人精品九色| 久久香蕉精品热| 日本免费a在线| 亚洲人成网站高清观看| 黄片wwwwww| 乱系列少妇在线播放| 免费观看精品视频网站| 亚洲电影在线观看av| ponron亚洲| 亚洲av五月六月丁香网| 搡女人真爽免费视频火全软件 | 国产三级中文精品| av视频在线观看入口| 日韩大尺度精品在线看网址| 日韩欧美国产一区二区入口| 91久久精品国产一区二区成人| 久久久久九九精品影院| 久久久精品大字幕| 一区二区三区高清视频在线| 精品久久国产蜜桃| 俺也久久电影网| 最近最新中文字幕大全电影3| 国产乱人伦免费视频| 国产高清激情床上av| 欧美xxxx黑人xx丫x性爽| 男人狂女人下面高潮的视频| 人妻制服诱惑在线中文字幕| 久久这里只有精品中国| 18禁裸乳无遮挡免费网站照片| 亚洲不卡免费看| 日本黄大片高清| 欧美日韩中文字幕国产精品一区二区三区| 成年女人毛片免费观看观看9| 91狼人影院| 免费不卡的大黄色大毛片视频在线观看 | 国产精品爽爽va在线观看网站| 中文字幕免费在线视频6| 人人妻,人人澡人人爽秒播| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 乱人视频在线观看| 日韩欧美在线乱码| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 波多野结衣高清作品| 亚洲国产精品久久男人天堂| 搞女人的毛片| 欧美日本亚洲视频在线播放| 搞女人的毛片| 国产 一区精品| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 精品福利观看| 天堂√8在线中文| 成人性生交大片免费视频hd| 国产av麻豆久久久久久久| 日韩欧美精品v在线| 午夜激情福利司机影院| 亚洲av中文av极速乱 | 中文字幕久久专区| 国产亚洲91精品色在线| 国产高清三级在线| 色吧在线观看| 十八禁国产超污无遮挡网站| 亚洲av美国av| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 久99久视频精品免费| 精品一区二区三区视频在线观看免费| 性插视频无遮挡在线免费观看| 欧美高清成人免费视频www| 日本精品一区二区三区蜜桃| 老师上课跳d突然被开到最大视频| 91av网一区二区| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 久99久视频精品免费| av在线蜜桃| 毛片女人毛片| 亚洲国产高清在线一区二区三| 成人午夜高清在线视频| 成人国产综合亚洲| 亚洲av第一区精品v没综合| 日日摸夜夜添夜夜添av毛片 | 免费人成在线观看视频色| 高清日韩中文字幕在线| 级片在线观看| 免费不卡的大黄色大毛片视频在线观看 | 91精品国产九色| 在线观看av片永久免费下载| 美女 人体艺术 gogo| 成人毛片a级毛片在线播放| 美女cb高潮喷水在线观看| 99视频精品全部免费 在线| 国产精品久久视频播放| 熟女电影av网| 桃红色精品国产亚洲av| 国产一区二区三区视频了| 免费观看的影片在线观看| 熟女电影av网| 午夜精品在线福利| 欧美一区二区精品小视频在线| 久久久久九九精品影院| 999久久久精品免费观看国产| 两个人视频免费观看高清| 日韩国内少妇激情av| 午夜亚洲福利在线播放| 成人高潮视频无遮挡免费网站| 天堂动漫精品| 免费在线观看日本一区| 亚洲av日韩精品久久久久久密| 日本一本二区三区精品| 深夜精品福利| 美女被艹到高潮喷水动态| 哪里可以看免费的av片| 日日摸夜夜添夜夜添小说| 国产精品一区二区三区四区久久| 国内久久婷婷六月综合欲色啪| 精品久久久久久,| 夜夜夜夜夜久久久久| 波多野结衣高清无吗| 亚洲人成伊人成综合网2020| 亚洲精品一区av在线观看| 在线看三级毛片| 自拍偷自拍亚洲精品老妇| 最新在线观看一区二区三区| 国产真实伦视频高清在线观看 | 婷婷色综合大香蕉| 欧美极品一区二区三区四区| 久久亚洲精品不卡| 午夜久久久久精精品| 久久中文看片网| 国产色爽女视频免费观看| 免费不卡的大黄色大毛片视频在线观看 | 极品教师在线视频| 国产精品av视频在线免费观看| 特大巨黑吊av在线直播| 日韩精品有码人妻一区| 久久久久久久久久久丰满 | 非洲黑人性xxxx精品又粗又长| 两个人的视频大全免费| 中文字幕精品亚洲无线码一区| 夜夜看夜夜爽夜夜摸| 欧美黑人巨大hd| 亚洲精华国产精华液的使用体验 | av在线亚洲专区| 在线观看美女被高潮喷水网站| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 美女被艹到高潮喷水动态| 久久久久久久午夜电影| 日日摸夜夜添夜夜添av毛片 | av在线亚洲专区| 99久久中文字幕三级久久日本| 成人鲁丝片一二三区免费| 97人妻精品一区二区三区麻豆| 51国产日韩欧美| 69av精品久久久久久| 欧美不卡视频在线免费观看| 中国美女看黄片| 熟女电影av网| 能在线免费观看的黄片| 久久人人精品亚洲av| 国产69精品久久久久777片| 国产午夜精品久久久久久一区二区三区 | 狠狠狠狠99中文字幕| 日韩,欧美,国产一区二区三区 | 成人国产综合亚洲| 成人无遮挡网站| 99视频精品全部免费 在线| 国产精华一区二区三区| 1000部很黄的大片| 99热这里只有是精品在线观看| 久久热精品热| 又爽又黄无遮挡网站| 久久精品久久久久久噜噜老黄 | 2021天堂中文幕一二区在线观| 无人区码免费观看不卡| 日本与韩国留学比较| 老司机福利观看| 亚洲综合色惰| 亚洲三级黄色毛片| 欧美性猛交╳xxx乱大交人| 成人永久免费在线观看视频| 亚洲精品影视一区二区三区av| 亚洲av五月六月丁香网| 给我免费播放毛片高清在线观看| 悠悠久久av| 亚洲人与动物交配视频| 18禁黄网站禁片免费观看直播| 免费在线观看成人毛片| 一a级毛片在线观看| 精品人妻一区二区三区麻豆 | 精品午夜福利在线看| 看片在线看免费视频| or卡值多少钱| 亚洲av五月六月丁香网| 一本一本综合久久| av在线观看视频网站免费| 人妻久久中文字幕网| 99久久精品热视频| 国产精品电影一区二区三区| 特级一级黄色大片| 偷拍熟女少妇极品色| 欧美不卡视频在线免费观看| 亚洲最大成人手机在线| 日本五十路高清| 91精品国产九色| 亚洲综合色惰| 99精品在免费线老司机午夜| 淫秽高清视频在线观看| 欧美不卡视频在线免费观看| 久久久色成人| 村上凉子中文字幕在线| 亚洲国产精品久久男人天堂| АⅤ资源中文在线天堂| 内射极品少妇av片p| 久久天躁狠狠躁夜夜2o2o| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 欧美日韩黄片免| 久久久久九九精品影院| 亚洲av五月六月丁香网| 日本免费a在线| 亚洲中文日韩欧美视频| 国国产精品蜜臀av免费| 精品国内亚洲2022精品成人| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| 美女免费视频网站| 国产精品一区www在线观看 | 欧美激情在线99| av在线亚洲专区| 国国产精品蜜臀av免费| 日韩欧美在线乱码| 亚洲黑人精品在线| 中文亚洲av片在线观看爽| 嫁个100分男人电影在线观看| 国产探花极品一区二区| 男插女下体视频免费在线播放| 又紧又爽又黄一区二区| 最后的刺客免费高清国语| 国产乱人视频| 国产精品伦人一区二区| 婷婷色综合大香蕉| 九色国产91popny在线| 久久精品综合一区二区三区| 色综合色国产| 成人av在线播放网站| 午夜福利在线在线| 欧美极品一区二区三区四区| 久久人人爽人人爽人人片va| 国产探花极品一区二区| 久久九九热精品免费| 国产三级在线视频| 国产aⅴ精品一区二区三区波| 最新中文字幕久久久久| 国产欧美日韩精品一区二区| 女人十人毛片免费观看3o分钟| 伦理电影大哥的女人| 精品一区二区三区视频在线| 精品久久久久久,| 老师上课跳d突然被开到最大视频| 日本与韩国留学比较| 在线播放国产精品三级| 国国产精品蜜臀av免费| 国产白丝娇喘喷水9色精品| 麻豆成人午夜福利视频| 亚洲国产精品sss在线观看| 十八禁国产超污无遮挡网站| 97超视频在线观看视频| 亚洲,欧美,日韩| 日韩欧美精品免费久久| 神马国产精品三级电影在线观看| 午夜爱爱视频在线播放| 最新中文字幕久久久久| a级一级毛片免费在线观看| 日韩欧美免费精品| 人妻制服诱惑在线中文字幕| 国产成年人精品一区二区| 欧美zozozo另类| 国产伦在线观看视频一区| a级毛片a级免费在线| 欧美+日韩+精品| 免费大片18禁| 欧美精品啪啪一区二区三区| 国产一区二区在线观看日韩| 久99久视频精品免费| 91午夜精品亚洲一区二区三区 | 看免费成人av毛片| 国产精品久久视频播放| 国产老妇女一区| 亚洲四区av| av在线蜜桃| 午夜激情欧美在线| 精品久久久久久,| 亚洲男人的天堂狠狠| 高清日韩中文字幕在线| 日韩大尺度精品在线看网址| 琪琪午夜伦伦电影理论片6080| 国产精品电影一区二区三区| 丰满乱子伦码专区| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 亚洲av一区综合| 黄色配什么色好看| 日韩欧美精品免费久久| 国产免费一级a男人的天堂| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 欧美高清性xxxxhd video| 国产精品,欧美在线| 午夜精品一区二区三区免费看| 五月玫瑰六月丁香| 精品久久久久久久人妻蜜臀av| 久久精品综合一区二区三区| 国产伦在线观看视频一区| 免费人成视频x8x8入口观看| 国产视频一区二区在线看| av女优亚洲男人天堂| 日本与韩国留学比较| 男人狂女人下面高潮的视频| 久久精品国产亚洲av香蕉五月| 欧洲精品卡2卡3卡4卡5卡区| 99热网站在线观看| 99精品在免费线老司机午夜| 俄罗斯特黄特色一大片| 成年版毛片免费区| 国产真实乱freesex| 九九久久精品国产亚洲av麻豆| 亚洲av中文字字幕乱码综合| 亚洲精品影视一区二区三区av| 在线观看av片永久免费下载| 尤物成人国产欧美一区二区三区| 久久久久久大精品| 亚洲电影在线观看av| 成熟少妇高潮喷水视频| 欧美日韩中文字幕国产精品一区二区三区| 国产三级中文精品| av黄色大香蕉| 国产黄色小视频在线观看| 国产精品一区二区性色av| 免费黄网站久久成人精品| 舔av片在线| 三级国产精品欧美在线观看| 成人二区视频| 久久久久久国产a免费观看| 女同久久另类99精品国产91| 久久久国产成人精品二区| 色哟哟·www| 伦理电影大哥的女人| xxxwww97欧美| 99热只有精品国产| 伦理电影大哥的女人| 麻豆精品久久久久久蜜桃| 国产亚洲91精品色在线| av福利片在线观看| 国产精品伦人一区二区| 欧美一区二区国产精品久久精品| 黄色配什么色好看| 18禁黄网站禁片午夜丰满| 国产精品免费一区二区三区在线| 国产精品不卡视频一区二区| 好男人在线观看高清免费视频| 精品久久久久久久久亚洲 | 国产精品精品国产色婷婷|