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

    低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液的代謝物變化

    2022-05-16 05:11:20李青林張文濤徐慧孫京京
    中國農(nóng)業(yè)科學(xué) 2022年8期
    關(guān)鍵詞:韌皮部木質(zhì)部汁液

    李青林,張文濤,徐慧,孫京京

    低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液的代謝物變化

    李青林*,張文濤,徐慧,孫京京

    江蘇大學(xué)農(nóng)業(yè)工程學(xué)院,江蘇鎮(zhèn)江 212013

    【背景】木質(zhì)部汁液負(fù)責(zé)為地上部分運輸水和營養(yǎng)物質(zhì),韌皮部汁液則主要負(fù)責(zé)光合產(chǎn)物的運輸,在植株受到非生物脅迫時,它們會發(fā)生不同的變化?!灸康摹垦芯康土酌{迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)各種代謝物含量變化,探究黃瓜在低磷脅迫下的代謝變化規(guī)律,為黃瓜磷匱乏狀態(tài)診斷提供新思路與依據(jù)?!痉椒ā坎杉A姿剑?.25 mmol?L-1)與低磷脅迫(0.3125 mmol?L-1)下黃瓜植株的木質(zhì)部與韌皮部汁液,經(jīng)三甲基硅烷化衍生化處理后用氣相色譜質(zhì)譜聯(lián)用技術(shù)(GC-MS)進(jìn)行檢測和定性定量分析。【結(jié)果】低磷脅迫下木質(zhì)部汁液內(nèi)糖類沒有發(fā)生顯著變化,多數(shù)氨基酸的含量與正常情況相比明顯上升,有機(jī)酸中的蘋果酸呈顯著下降趨勢,軟脂酸、油酸等則顯著高于正常水平;低磷脅迫下的韌皮部汁液中,糖類與氨基酸多數(shù)低于正常水平,草酸、琥珀酸等有機(jī)酸的含量也明顯降低。此外,為篩選黃瓜木質(zhì)部與韌皮部汁液內(nèi)潛在的生物標(biāo)志物,使用OPLS-DA與獨立樣本T檢測方法對GC-MS試驗數(shù)據(jù)進(jìn)行進(jìn)一步分析,結(jié)果表明,在黃瓜木質(zhì)部汁液檢出8種低磷脅迫下的潛在生物標(biāo)志物(VIP>1,<0.05),分別為尸胺、蘋果酸、四氨基丁酸、軟脂酸、肌醇、纖維二糖、N, N-二甲基十二酰胺與單硬脂酸甘油酯;在黃瓜韌皮部汁液中檢出7種低磷脅迫下的潛在生物標(biāo)志物(VIP>1,<0.05),分別為絲氨酸、亞油酸、10, 13-環(huán)氧-11-甲基十八碳-10, 12-二烯酸、脫氫樅酸、N, N-二甲基十二酰胺、2, 2'-亞甲基雙(6-叔丁基-4-甲基苯酚)、L-2-哌啶酸?!窘Y(jié)論】低磷脅迫下,黃瓜的木質(zhì)部與韌皮部汁液內(nèi)的糖類、氨基酸和有機(jī)酸均發(fā)生了變化,但是變化趨勢并不相同。分別在木質(zhì)部和韌皮部汁液中篩選出8種和7種生物標(biāo)志物,可為黃瓜磷匱乏狀態(tài)診斷提供新思路與依據(jù)。生物標(biāo)志物中,蘋果酸、四氨基丁酸、軟脂酸、絲氨酸與亞油酸的檢測手段已經(jīng)較為成熟,以其為檢測對象的黃瓜磷匱乏狀態(tài)診斷方法具有較好的實用性。

    黃瓜;木質(zhì)部;韌皮部;磷脅迫;代謝組學(xué);生物標(biāo)志物

    0 引言

    【研究意義】木質(zhì)部與韌皮部汁液在植物養(yǎng)分和信號的傳輸與分配中起著關(guān)鍵作用[1],木質(zhì)部汁液主要將根系吸收的水分和礦質(zhì)養(yǎng)分從土壤輸送到植物的地上部分,其中還含有許多其他有機(jī)化合物,包括羧酸鹽、激素、氨基酸、肽和蛋白質(zhì)[2];韌皮部汁液的主要作用是將光合產(chǎn)物從進(jìn)行光合作用的源器官運輸?shù)綆炱鞴?,其主要成分是糖,還有蛋白質(zhì)、氨基酸、酰氨、維生素、無機(jī)鹽和植物激素等[3]。研究低磷脅迫下植株木質(zhì)部與韌皮部汁液內(nèi)代謝物含量的變化規(guī)律,對加深植物在低磷脅迫下代謝變化規(guī)律的認(rèn)識,探究新的植物磷匱乏檢測手段具有重要意義?!厩叭搜芯窟M(jìn)展】在受到非生物脅迫時,植株會調(diào)整代謝網(wǎng)絡(luò)以保證植株正常生長發(fā)育[4],從而影響體內(nèi)各種代謝物的含量。代謝組學(xué)可以通過對植株體內(nèi)的各種小分子化合物進(jìn)行定性定量分析,研究代謝物的變化規(guī)律[5],進(jìn)而反映植物在非生物脅迫下的生理變化。同時,代謝組學(xué)分析能夠發(fā)現(xiàn)在不同組間具有顯著性差異的某種或若干種代謝物,即為生物標(biāo)志物,其在生物體代謝過程和生理變化機(jī)制的研究中具有重要作用[6]。Sung等[7]通過對番茄葉片和根系在氮、磷或鉀脅迫下的新陳代謝變化進(jìn)行分析,發(fā)現(xiàn)這3種礦物質(zhì)脅迫都會影響能量生產(chǎn)和氨基酸代謝。zhao等[8]通過高效液相色譜質(zhì)譜分析低溫下草菇的代謝變化,確定了216種潛在生物標(biāo)志物,揭示了耐低溫草菇菌株的耐低溫機(jī)理。Zhang等[9]使用核磁共振波譜分析鹽脅迫下煙草植株的代謝響應(yīng)。對植株進(jìn)行代謝組學(xué)研究,能探究植株在不同環(huán)境影響下的代謝變化規(guī)律。對生物標(biāo)志物的含量進(jìn)行測量,能在一定程度上反映植株的生理狀態(tài),為植株生長環(huán)境的調(diào)控提供指導(dǎo)?!颈狙芯壳腥朦c】近年來傳感器工藝不斷發(fā)展,其在農(nóng)業(yè)上也有所運用[10],借助微針傳感器可以在低損傷的情況下對木質(zhì)部或韌皮部汁液的某些生物特性參數(shù)進(jìn)行測量,以黃瓜木質(zhì)部或韌皮部汁液為檢測對象的低損傷檢測方法成為可能。但是現(xiàn)階段針對黃瓜的代謝組學(xué)研究主要對象多為黃瓜的根部或葉片[11-13],針對黃瓜木質(zhì)部與韌皮部汁液在非生物脅迫下的代謝變化規(guī)律和生物標(biāo)志物的研究較少?!緮M解決的關(guān)鍵問題】采用GC-MS技術(shù)對正常與低磷脅迫下黃瓜木質(zhì)部和韌皮部汁液內(nèi)的代謝物進(jìn)行定性定量,分析低磷脅迫下汁液內(nèi)糖類、氨基酸和有機(jī)酸含量的變化,并篩選出潛在生物標(biāo)志物,旨在使用代謝組學(xué)技術(shù)加深對低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)代謝物組分變化規(guī)律的認(rèn)識,為通過木質(zhì)部或韌皮部汁液進(jìn)行黃瓜磷養(yǎng)分狀態(tài)的診斷提供依據(jù),為新的黃瓜磷匱乏狀態(tài)檢測方法提供方向,為黃瓜的精細(xì)化種植提供參考。

    1 材料與方法

    試驗于2020年在江蘇大學(xué)現(xiàn)代農(nóng)業(yè)裝備與技術(shù)教育部重點試驗室玻璃溫室進(jìn)行。

    1.1 試驗材料與儀器

    試驗材料:黃瓜品種為‘津研四號’,種植期間溫室溫度范圍為18—35℃,光照采用自然光照,種植期間標(biāo)準(zhǔn)光照下平均日照時間為5.13 h,平均日輻照量13.2 MJ?m-2。

    試驗儀器與設(shè)備:GCMS-QP2010氣相色譜質(zhì)譜分析儀(日本島津公司);HP-5MS毛細(xì)管柱(30 m×0.25 mm×0.25 μm,美國安捷倫公司);超聲洗浴儀(上??茖?dǎo)公司);臺式高速冷凍離心機(jī)(Neofuge 18R,力康生物醫(yī)療科技控股有限公司);FreeZone12Plus冷凍干燥機(jī)(北京中科科爾儀器有限公司);DH-101-2BY電熱鼓風(fēng)干燥箱(天津中環(huán)電爐股份有限公司)。

    1.2 樣本培育

    采用鎮(zhèn)江培蕾基質(zhì)科技發(fā)展有限公司的植物培養(yǎng)基進(jìn)行育苗,待黃瓜苗生長到3片真葉時進(jìn)行移栽工作,移栽后采用無土栽培方式進(jìn)行樣本培育,基質(zhì)為2—4 mm珍珠巖,使用霍格蘭配方營養(yǎng)液進(jìn)行澆灌,營養(yǎng)液的組成成分:A溶液(Ca(NO3)2·4H2O 945 mg?L-1、KNO379 mg?L-1、NH4NO340 mg?L-1、Fe-EDTA (13.1%) 7 mg?L-1);B溶液(KNO3604 mg?L-1、KH2PO4170 mg?L-1、MgSO4·7H2O 339 mg?L-1、MnSO4?4H2O 1.7 mg?L-1、ZnSO4·7H2O 1.45 mg?L-1、Na2B4O7·10H2O 2.45 mg?L-1、CuSO4·5H2O 0.19 mg?L-1、Na2MoO4·2H2O 0.12 mg?L-1),將A、B兩溶液等量混合進(jìn)行澆灌,營養(yǎng)液的補(bǔ)充時間為每天早上7:00—8:00??傸S瓜培育數(shù)為60株,當(dāng)黃瓜生長到七葉一心時,選取長勢基本一致的12株黃瓜幼苗(株高、莖粗上下偏差在20%以內(nèi)),并將其分為試驗組(6株)與對照組(6株)。對試驗組進(jìn)行低磷脅迫處理,缺磷處理下黃瓜試驗組的營養(yǎng)液中磷含量為正常水平的25%,實施方式是將B溶液中的KH2PO4降低至42.5mg?L-1,并使用63.75mg?L-1的K2SO4保證K+濃度正常,對照組不做任何處理,按照原營養(yǎng)液配方進(jìn)行培育。

    1.3 黃瓜木質(zhì)部與韌皮部汁液取樣

    每周對黃瓜植株的株高和離地3 cm處的莖粗進(jìn)行測量,確定黃瓜植株每周的株高和莖粗增長量,測量結(jié)果見表1,可以看出在低磷脅迫的第一周,植株株高和莖粗的增長就已有所放緩,但是并沒有顯著性差異。在第二周內(nèi),植株的株高和莖粗增長量均與正常處理植株有顯著區(qū)別,證實了兩周的低磷脅迫會導(dǎo)致黃瓜植株的性狀與正常狀態(tài)產(chǎn)生差異[14],表明黃瓜的生理代謝已經(jīng)產(chǎn)生了明顯變化,此時的檢測結(jié)果更為可靠,更能反映黃瓜木質(zhì)部與韌皮部汁液內(nèi)代謝物在低磷脅迫下的含量變化規(guī)律,因此,本研究在低磷脅迫的14 d后進(jìn)行取樣。取樣時,將植株莖稈切斷,從近地段切口取得汁液為木質(zhì)部汁液[15];從另一端切口取得的為韌皮部汁液[16]。

    表1 不同磷處理下黃瓜植株的表型參數(shù)每周變化量

    同列不同小寫字母表示差異顯著(<0.05)Different lowercase letters in the same column indicate significant differences (<0.05)

    使用消毒后的鋒利刀片,在離地3—5 cm處切斷黃瓜莖稈,使用去離子水清洗兩端創(chuàng)傷面,再使用消毒后的脫脂棉吸干水分,棄去第一滴后開始分別從兩端切口收集木質(zhì)部與韌皮部汁液。收集到的汁液放于樣本管中,置于液氮暫存。每株黃瓜取得木質(zhì)部汁液樣本與韌皮部汁液樣本各1個,取樣量為100—300 μL,共取得樣本24個(表2)。

    1.4 樣本三甲基化衍生化處理

    每例樣本30 μL,加入170 μL去離子水,800 μL甲醇乙腈(1﹕1,v/v),53 KHz冰浴超聲1 h, -20℃孵育1 h沉淀蛋白,14 000×4℃離心20 min,取上清液400 μL,加入60 μL核糖醇(0.2 mg?mL-1)作為內(nèi)參代謝物,渦旋10 s,然后進(jìn)行真空冷凍干燥,得到干燥提取物;往干燥后的樣本中加入80 μL的15 mg?mL-1甲氧胺吡啶溶液復(fù)溶,渦旋振蕩30 s。使用恒溫箱37℃反應(yīng)90 min然后快速加入80 μL的BSTFA試劑(含1% TMCS);使用恒溫箱70℃反應(yīng)60 min,取出后等待測定。

    表2 黃瓜分組表

    組別名中,第一個英文字母“P”代表試驗方法為磷脅迫,中間的數(shù)字“25”“100”分別代表磷含量為25%與100%,最后的字母“X”與“P”分別為木質(zhì)部和韌皮部的英文首字母(Xylem,Phloem)。下同

    n the group aliases, the first English letter “P” represents that the test method is phosphorus stress, the numbers “25” and “100” in the middle represent the phosphorus content of 25% and 100%, respectively, and the last letters “X” and “P” are the English initials of xylem and phloem, respectively.The same as below

    1.5 GC-MS分析

    使用島津GCMS-QP2010氣相色譜質(zhì)譜分析儀對前處理后的樣本進(jìn)行檢測。

    色譜條件:進(jìn)樣量:1 μL;載氣:氦氣;前進(jìn)樣口吹掃流速:3 mL·min-1;流速:1 mL·min-1;柱溫:初始溫度90℃保持15 s,以10℃·min-1升至180℃,再以5℃·min-1升至240℃,再以20℃·min-1繼續(xù)升至285℃保持11.5 min。

    質(zhì)譜條件:前進(jìn)樣口溫度:280℃;傳輸線溫度:245℃;離子源溫度:220℃;電離電壓:-70 eV;掃描方式:20—600 m/z(全掃描);掃描速度:100 spectra/s;溶劑延遲:480 s。

    1.6 數(shù)據(jù)處理與分析

    1.6.1 代謝物的定性定量 使用島津儀器的GCMS Real Time Analysis軟件,結(jié)合NIST數(shù)據(jù)庫確定檢出代謝物名稱,根據(jù)核糖醇內(nèi)標(biāo)物的濃度,使用內(nèi)標(biāo)標(biāo)準(zhǔn)化法確定每種代謝物的濃度。

    1.6.2 代謝物的PCA分析 PCA是多元統(tǒng)計分析中最常見的無監(jiān)督模式的數(shù)據(jù)分析方法,能將大量相關(guān)變量通過線性變換轉(zhuǎn)化為一組最能代表數(shù)據(jù)特征的不相關(guān)的變量,展示數(shù)據(jù)組內(nèi)的重復(fù)性和組間的差異性,評估數(shù)據(jù)的可重復(fù)性,發(fā)現(xiàn)異常值,提高數(shù)據(jù)精準(zhǔn)度[17]。本研究使用PCA判斷低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)代謝物含量是否發(fā)生了明顯變化。

    1.6.3 糖、氨基酸和有機(jī)酸含量分析 糖、氨基酸和有機(jī)酸都是植物正常生長發(fā)育必不可缺的營養(yǎng)成分,糖為植物的生長發(fā)育儲備和提供能量[18];氨基酸參與植物體內(nèi)各種生理代謝活動,增強(qiáng)植株抗逆性[19];有機(jī)酸同樣涉及到抗逆性等生理代謝反應(yīng)[20]。為了解植株在磷匱乏時的生理反應(yīng),將黃瓜木質(zhì)部與韌皮部汁液在正常與低磷脅迫下的糖、氨基酸和有機(jī)酸的含量進(jìn)行對比,分析是否有顯著性差異。

    1.6.4 潛在生物標(biāo)志物篩選 使用正交偏最小二乘法(OPLS-DA)和獨立樣本t檢測對數(shù)據(jù)進(jìn)行處理并篩選潛在生物標(biāo)志物。OPLS-DA能選擇區(qū)分各組的特征變量,確定樣本間的關(guān)系,建立代謝物表達(dá)量與分組關(guān)系之間的模型,獲取組間差異信息,預(yù)測樣品分組[21]。獨立樣本t檢測用于分析定性數(shù)據(jù)與定量數(shù)據(jù)之間的關(guān)系,可檢驗數(shù)據(jù)間的差異[22]。將OPLS-DA與獨立樣本T檢測相結(jié)合,是常見的生物標(biāo)志物篩選方式[23]。

    2 結(jié)果

    2.1 木質(zhì)部與韌皮部汁液的主成分分析

    圖1分別為使用SIMCA-P軟件(14.1版,瑞典)得到的低磷脅迫與正常處理黃瓜木質(zhì)部汁液和韌皮部汁液內(nèi)代謝物的PCA得分圖,圖1-A是木質(zhì)部汁液的PCA分析,此時第一和第二主成分的解釋率分別為33.0%和19.7%;圖1-B是韌皮部汁液的PCA分析,第一和第二主成分的解釋率分別為37.0%和17.2%。根據(jù)樣本在PCA得分圖上的分布來看,木質(zhì)部和韌皮部汁液中的低磷組和正常磷組被明顯區(qū)分,說明磷脅迫下,黃瓜木質(zhì)部汁液和韌皮部汁液內(nèi)各種代謝物的含量變化明顯。

    2.2 磷脅迫對黃瓜木質(zhì)部與韌皮部汁液內(nèi)糖、氨基酸和有機(jī)酸含量的影響

    2.2.1 對糖組分含量的影響 圖2為低磷脅迫下木質(zhì)部與韌皮部汁液內(nèi)糖分的變化。其中,木質(zhì)部汁液內(nèi)檢出蔗糖、果糖、葡萄糖、山梨醇共4種糖類;韌皮部汁液中,除了同樣能檢出木質(zhì)部中的4種糖類外,還檢出海藻糖。在韌皮部汁液中,糖組分的含量普遍高于木質(zhì)部汁液,這符合木質(zhì)部與韌皮部汁液分別運輸根系合成物與葉片合成物的功能特性。低磷脅迫下,木質(zhì)部汁液內(nèi)的各種糖類均高于正常磷供應(yīng)情況,平均提升15%以上,但在統(tǒng)計學(xué)上,木質(zhì)部汁液內(nèi)的各種糖類含量均未發(fā)生顯著變化。低磷脅迫下韌皮部汁液中不同糖組分的變化趨勢則略有不同,其中蔗糖、海藻糖和山梨醇的含量顯著降低(< 0.05),分別是正常水平的72.3%、52.1%和50.6%,葡萄糖和果糖的含量雖然有所升高,但未達(dá)到顯著性差異水平。

    圖1 木質(zhì)部(A)與韌皮部(B)汁液的PCA得分圖

    不同小寫字母表示差異顯著(P<0.05)。下同Different lowercase letters mean significant difference (P<0.05).The same as below

    2.2.2 對氨基酸含量的影響 木質(zhì)部汁液中主要檢出谷氨酸、甘氨酸等9種氨基酸,韌皮部共檢出14種氨基酸,比木質(zhì)部汁液多檢出鳥氨酸、脯氨酸等(圖3)。低磷脅迫的木質(zhì)部汁液中,纈氨酸的含量顯著低于正常水平(<0.05),僅為正常狀態(tài)的43.5%;而半胱氨酸、谷氨酸、蘇氨酸、天冬氨酸、賴氨酸則顯著高于正常狀態(tài)(<0.05),分別是正常水平的165.3%、146.4%、200.5%、122.6%和137.4%。低磷處理下韌皮部汁液檢出的氨基酸中,僅鳥氨酸的含量顯著高于正常水平(<0.05),而蘇氨酸、甘氨酸、絲氨酸、丙氨酸、脯氨酸、纈氨酸、谷氨酸和瓜氨酸均顯著低于正常狀態(tài)(<0.05),分別是正常水平的83.0%、80.7%、59.3%、61.8%、83.1%、79.0%、59.0%和73.9%。

    圖3 低磷脅迫下木質(zhì)部與韌皮部汁液內(nèi)氨基酸含量的變化

    2.2.3 對常見有機(jī)酸含量的影響 木質(zhì)部與韌皮部汁液內(nèi)的有機(jī)酸種類有一定差別(圖4)。低磷脅迫下的木質(zhì)部汁液中,蘋果酸的含量顯著低于正常水平(<0.05),僅為正常狀態(tài)的38.2%,而乳酸、軟脂酸、油酸、肉豆蔻酸、四氨基丁酸、蘇糖酸的含量則顯著高于正常水平(<0.05),分別為正常水平的128.8%、108.0%、143.5%、111.8%、159.3%和172.5%。低磷脅迫的韌皮部汁液中,草酸、琥珀酸、亞油酸、肉豆蔻酸、四氨基丁酸、蘇糖酸和丙二酸的含量與正常水平相比有顯著降低(<0.05),分別為正常水平的79.4%、72.1%、55.2%、87.3%、64.3%、67.3%和57.8%,未發(fā)現(xiàn)顯著高于正常水平的有機(jī)酸。

    2.3 木質(zhì)部與韌皮部汁液內(nèi)代謝物的潛在生物標(biāo)志物篩選

    2.3.1 汁液內(nèi)代謝物的OPLS-DA分析 對木質(zhì)部與韌皮部汁液內(nèi)代謝物進(jìn)行OPLS-DA分析,得到OPLS- DA得分圖(圖5-A、B)。在OPLS-DA中,R2Y和Q2為模型的評價指標(biāo),R2Y和Q2均大于0.5,說明模型效果較好。對于木質(zhì)部汁液,R2Y=0.992,Q2= 83.2%;對于韌皮部汁液,R2Y=1,Q2=97.8%,說明OPLS-DA得到的模型良好。

    圖4 低磷脅迫下木質(zhì)部與韌皮部汁液內(nèi)有機(jī)酸含量的變化

    使用OPLS-DA模型生成木質(zhì)部與韌皮部汁液代謝物的S-plot圖(圖5-C、D)。S-plot圖能同時體現(xiàn)每個變量與模型某一成分之間的協(xié)變性和相關(guān)性,它們分別代表變量對模型的貢獻(xiàn)大小和變量的權(quán)重系數(shù),可以用于篩選有意義的變量。S-plot圖中S型曲線上的每一個點代表一個代謝物,點與原點的距離越大,對分類的貢獻(xiàn)就越大,越靠近左下角和右上角的代謝物對樣品的分類貢獻(xiàn)越大,可以用變量權(quán)重重要性排序值(VIP)評價變量的貢獻(xiàn)大小,VIP>1即認(rèn)為此代謝物有較大貢獻(xiàn)。

    為了防止模型過度擬合,使用200次置換檢驗驗證模型(圖5-E、F),置換檢驗要求Q2的回歸截距小于0.05,圖中木質(zhì)部汁液中Q2的回歸截距為-0.77,韌皮部汁液中為-1.16,均低于0.05,說明兩個模型均可靠,不存在過擬合情況。

    2.3.2 木質(zhì)部與韌皮部汁液的獨立樣本t檢測和潛在生物標(biāo)志物篩選 結(jié)合OPLS-DA的VIP值與獨立樣本t檢測的值,按照VIP>1且<0.05的要求篩選磷脅迫下木質(zhì)部與韌皮部汁液的代謝物,最終分別篩選出8種和7種代謝物。木質(zhì)部汁液篩選出的8種代謝物(VIP>1,<0.05)分別為尸胺、蘋果酸、四氨基丁酸、軟脂酸、肌醇、纖維二糖、N, N-二甲基十二酰胺與單硬脂酸甘油酯(圖6)。韌皮部汁液篩選出的7種代謝物(VIP>1,<0.05)分別為絲氨酸、亞油酸、10, 13-環(huán)氧-11-甲基十八碳-10, 12-二烯酸、脫氫樅酸、N, N-二甲基十二酰胺、2, 2'-亞甲基雙(6-叔丁基-4-甲基苯酚)、L-2-哌啶酸(圖7)。這些符合條件的代謝物均在低磷與正常處理兩組樣本的分離中發(fā)揮了重要作用,被認(rèn)為是低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)的潛在生物標(biāo)志物。

    圖5 木質(zhì)部與韌皮部汁液的OPLS-DA得分圖(木質(zhì)部A、韌皮部B),S-plot圖(木質(zhì)部C、韌皮部D),200次置換檢驗結(jié)果(木質(zhì)部E、韌皮部F)

    3 討論

    3.1 磷脅迫影響黃瓜中糖、氨基酸和有機(jī)酸代謝

    黃瓜在低磷情況下會對代謝進(jìn)行適當(dāng)調(diào)整,以保證黃瓜正常生長。根據(jù)低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)糖分變化,可以看出低磷脅迫阻礙了植株地上組織內(nèi)蔗糖、海藻糖和山梨醇的合成,導(dǎo)致向根部供應(yīng)的養(yǎng)分大量減少。根據(jù)低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)氨基酸的變化,可以看出黃瓜韌皮部汁液中多種氨基酸的含量均減小,說明低磷脅迫同樣抑制植株地上組織的氨基酸合成。與此同時,植株根系增加了多種氨基酸的合成以抵抗外界脅迫對正常生長代謝的干擾,并對地上部分氨基酸合成量的減少進(jìn)行補(bǔ)償,比如大量生產(chǎn)谷氨酸以增加對磷的吸收效率,大量生產(chǎn)蘇氨酸增加對不利條件的耐受性。根據(jù)低磷脅迫下黃瓜木質(zhì)部與韌皮部汁液內(nèi)有機(jī)酸的變化,可以看出低磷脅迫影響各種有機(jī)酸的合成,特別是四氨基丁酸、蘇糖酸、肉豆蔻酸,它們在低磷脅迫下的木質(zhì)部和韌皮部汁液中的含量均與正常水平有顯著性差異,但是變化趨勢并不相同,這可能與植株不同部位應(yīng)對磷脅迫時的反應(yīng)機(jī)理不同有關(guān)。

    圖6 木質(zhì)部汁液潛在生物標(biāo)志物的箱型圖

    Fig.6 Boxplots of potential biomarkers of xylem sap

    圖7 韌皮部汁液潛在生物標(biāo)志物的箱型圖

    Fig.7 Boxplots of potential biomarkers of phloem sap

    3.2 生物標(biāo)志物可用于估測黃瓜磷營養(yǎng)供應(yīng)

    低磷脅迫下,植株會調(diào)節(jié)糖、氨基酸和有機(jī)酸的代謝,以減小磷供應(yīng)不足造成的不利影響,維持植株的正常生理代謝[24-25]。且這種變化會先于植株外部特征的變化,因此檢測植株的代謝物變化能夠較早反映植株的養(yǎng)分狀態(tài)。但相較于對所有代謝物進(jìn)行檢測,僅檢測生物標(biāo)志物即可反映植株的磷養(yǎng)分狀態(tài),且檢測流程更短,檢測效率更高。如Gao等[26]通過比較低氮、低磷脅迫下萵苣葉片各種代謝產(chǎn)物的含量變化并篩選出代謝物標(biāo)志物,為萵苣早期營養(yǎng)脅迫的診斷提供依據(jù)。但針對植物木質(zhì)部或韌皮部汁液進(jìn)行代謝組學(xué)分析的案例較少。

    本研究在低磷脅迫下的木質(zhì)部汁液中發(fā)現(xiàn)了8種具有明顯差異的代謝物,在韌皮部汁液中發(fā)現(xiàn)了7種具有明顯差異的代謝物,這些潛在生物標(biāo)志物能在一定程度上反映出黃瓜的磷素代謝狀況。在木質(zhì)部汁液中,參與眾多代謝途徑的關(guān)鍵代謝物蘋果酸的含量大幅下降,說明低磷脅迫影響了植物體內(nèi)的正常生長代謝,限制了植物生長發(fā)育;能增強(qiáng)植物抗逆性的尸胺與四氨基丁酸的含量顯著上升,說明植株正加強(qiáng)自身抗逆性以適應(yīng)低磷脅迫;作為重要生長因子的軟脂酸和肌醇的含量顯著上升,說明植株增加了生長因子的產(chǎn)出以維持生長發(fā)育。在韌皮部汁液中,增強(qiáng)植物抗逆性的L-2-哌啶酸顯著增加[27],也是植株抗逆性增強(qiáng)的體現(xiàn);參與植株代謝的絲氨酸和亞油酸顯著降低,顯示出低磷脅迫對于植株生理代謝反應(yīng)的抑制;誘導(dǎo)植株系統(tǒng)性獲得抗性的脫氫樅酸顯著降低[28],說明低磷脅迫下的植株對外界病菌局部侵染的抵抗力顯著下降。

    同時,與正常磷供應(yīng)狀態(tài)相比,低磷脅迫下木質(zhì)部汁液的8種潛在生物標(biāo)志物中有7種高于正常水平,而韌皮部汁液中的7種生物標(biāo)志物含量則有6種低于正常水平,結(jié)合磷脅迫對木質(zhì)部與韌皮部汁液內(nèi)糖、氨基酸和無機(jī)酸含量的不同影響,發(fā)現(xiàn)在同樣的低磷脅迫下,黃瓜木質(zhì)部與韌皮部汁液的響應(yīng)機(jī)制有明顯差異,其他植物在同種非生物脅迫下也出現(xiàn)不同部位響應(yīng)機(jī)制有明顯差異的情況[29-30]。

    3.3 易于實際檢測的生物標(biāo)志物

    目前,對于植物體內(nèi)代謝物的檢測方法有核磁共振(NMR)、GC-MS、LC-MS等,這些檢測方法都需要離體檢測,且檢測設(shè)備體型巨大、成本高昂。相對而言,借助敏感電極進(jìn)行檢測的傳感器雖然精度有所不足,但是更方便、廉價,具有較高的普適性?,F(xiàn)階段有大量用于氨基酸、有機(jī)酸等有機(jī)物檢測的傳感器,也有直接植入植物體內(nèi)、對植物某一生長指標(biāo)進(jìn)行長期檢測的傳感器。因此,相較于色譜、質(zhì)譜等方法,通過體積小、成本低的傳感器檢測黃瓜磷脅迫下的生物標(biāo)志物,檢測成本更低,檢測效率更高。目前,在木質(zhì)部與韌皮部汁液發(fā)現(xiàn)的生物標(biāo)志物中,蘋果酸、四氨基丁酸、軟脂酸、絲氨酸與亞油酸的傳感器檢測較為完善,相對于其他生物標(biāo)志物,它們的實用性更強(qiáng)。

    蘋果酸、四氨基丁酸、軟脂酸為木質(zhì)部汁液中的生物標(biāo)志物。DALAL等[31]基于蘋果酸濃度這一果實成熟的重要指標(biāo),通過將蘋果酸酶固定在工作電極表面,制備了羧化多壁碳納米管(C-MWCNT)修飾的絲網(wǎng)印刷電化學(xué)傳感器,可用于蘋果酸檢測。SANGUBOTLA等[32]以玉米汁為原料合成碳點(CDs),并用3-氨基苯基硼酸(APBA)和煙酰胺腺嘌呤二核苷酸磷酸(NADP+)進(jìn)行功能化,得到了CANs;制備的CANs和酶GABase(四氨基丁酸氨基轉(zhuǎn)移酶和琥珀酸半醛脫氫酶的混合物)能以較高的靈敏度對四氨基丁酸進(jìn)行檢測。DUVANOVA等[33]用PM聚酰亞胺(1, 2, 4, 5-苯四甲酸二酐和4, 4′-二氨基二苯氧基)修飾電極,制備出可以檢測植物油中軟脂酸濃度的分子印跡聚合物(MIPs)壓電傳感器。

    絲氨酸與亞油酸為韌皮部汁液中的生物標(biāo)志物。SINGH基于分子印跡技術(shù),制備功能化金納米粒子(AuNPs),并用其修飾EQCM金電極,使用4-氨基硫酚作為電聚合組分,使用半胱氨酸單元通過兩性離子和氫鍵間的相互作用為識別氨基酸提供配體,制備出一種能夠檢測絲氨酸的壓電傳感器[34]。LU等[35]同時基于膳食脂肪的疏水性,采用氧化石墨烯和氣味結(jié)合蛋白(OBPs)修飾絲網(wǎng)印刷電極,研制了一種用于脂肪味覺檢測的味覺生物傳感器,能夠用于亞油酸的檢測。

    綜上,現(xiàn)階段已經(jīng)有可用于蘋果酸、氨酪酸、軟脂酸、絲氨酸與亞油酸檢測的傳感器。相較于其他生物標(biāo)志物,在通過傳感器確定其在木質(zhì)部或韌皮部汁液中的含量,并根據(jù)生物標(biāo)志物的含量進(jìn)行黃瓜磷匱乏診斷,更易于生產(chǎn)實際應(yīng)用。

    4 結(jié)論

    低磷脅迫下,植物汁液內(nèi)糖類、四氨基丁酸與有機(jī)酸的含量受到明顯影響,且木質(zhì)部與韌皮部汁液對低磷脅迫的響應(yīng)機(jī)制有很大差異。其中,低磷脅迫下,木質(zhì)部汁液中的糖無顯著變化,氨基酸和有機(jī)酸多呈顯著增加的趨勢,而韌皮部汁液中糖、氨基酸和有機(jī)酸則以顯著降低為主。同時,在黃瓜韌皮部汁液中篩選出蘋果酸、四氨基丁酸等8種潛在生物標(biāo)志物,在韌皮部汁液中篩選出L-2-哌啶酸、絲氨酸等7種潛在生物標(biāo)志物,這些潛在生物標(biāo)志物可用于黃瓜磷養(yǎng)分供應(yīng)狀況的評估,為植物的養(yǎng)分檢測以及磷肥的精準(zhǔn)供應(yīng)提供指導(dǎo)。

    References

    [1] RODRíGUEZ-CELMA J, CEBALLOS-LAITA L, GRUSAK M A, ABADíA J, LóPEZ-MILLáN A F.Plant fluid proteomics: Delving into the xylem sap, phloem sap and apoplastic fluid proteomes.Biochimica et Biophysica Acta, 2016, 1864(8): 991-1002.doi: 10.1016/j.bbapap.2016.03.014.

    [2] LUCAS W J, GROOVER A, LICHTENBERGER R, FURUTA K, YADAV S R, HELARIUTTA Y, HE X Q, FUKUDA H, KANG J, BRADY S M, PATRICK J W, SPERRY J, YOSHIDA A, LóPEZ-MILLáN A F, GRUSAK M A, KACHROO P.The plant vascular system: Evolution, development and functions.Journal of Integrative Plant Biology, 2013, 55(4): 294-388.

    [3] 楊春, 譚太龍, 余佳玲, 廖瓊, 張曉龍, 張振華, 宋海星, 官春云.大氣CO2濃度倍增對油菜韌皮部汁液成分及根部氮素積累的影響.植物生態(tài)學(xué)報, 2014, 38(7): 776-784.doi:10.3724/SP.J.1258.2014.00073.

    YANG C, TAN T L, YU J L, LIAO Q, ZHANG X L, ZHANG Z H, SONG H X, GUAN C Y.Effects of atmospheric CO2enrichment on phloem sap composition and root nitrogen accumulation in oilseed rape.Chinese Journal of Plant Ecology, 2014, 38(7): 776-784.doi: 10.3724/SP.J.1258.2014.00073.(in Chinese)

    [4] ZHANG Y, MA X M, WANG X C, LIU J H, HUANG B Y, GUO X Y, XIONG S P, LA G X.UPLC-QTOF analysis reveals metabolomic changes in the flag leaf of wheat (L.) under low-nitrogen stress.Plant Physiology and Biochemistry, 2017, 111: 30-38.doi: 10.1016/j.plaphy.2016.11.009.

    [5] 許彥陽, 姚桂曉, 劉平香, 趙潔, 王昕璐, 孫君茂, 錢永忠.代謝組學(xué)在農(nóng)產(chǎn)品營養(yǎng)品質(zhì)檢測分析中的應(yīng)用.中國農(nóng)業(yè)科學(xué), 2019, 52(18): 3163-3176.doi: 10.3864/j.issn.0578-1752.2019.18.009.

    XU Y Y, YAO G X, LIU P X, ZHAO J, WANG X L, SUN J M, QIAN Y Z.Review on the application of metabolomic approaches to investigate and analysis the nutrition and quality of agro-products.Scientia Agricultura Sinica, 2019, 52(18): 3163-3176.doi: 10.3864/ j.issn.0578-1752.2019.18.009.(in Chinese)

    [6] 吳昊, 于小紅, 馬光朝, 馬致潔, 章從恩.基于LC-MS的甘草炮制雷公藤降低肝毒性的代謝組學(xué)研究.中草藥, 2020, 51(21): 5501-5508.doi: 10.7501/j.issn.0253-2670.2020.21.014.

    WU H, YU X H, MA G C, MA Z J, ZHANG C E.Metabolomics study on reduction of hepatotoxic ofprocessed by liquorice based on LC-MS.Acupuncture Research, 2020, 51(21): 5501-5508.doi: 10.7501/j.issn.0253-2670.2020.21.014.(in Chinese)

    [7] SUNG J, LEE S, LEE Y J, HA S, SONG B, KIM T, WATERS B M, KRISHNAN H B.Metabolomic profiling from leaves and roots of tomato (L.) plants grown under nitrogen, phosphorus or potassium-deficient condition.Plant Science, 2015, 241: 55-64.doi: 10.1016/j.plantsci.2015.09.027.

    [8] ZHAO X, CHEN M J, LI Z P, ZHAO Y, YANG H L, ZHA L, YU C X, WU Y J, SONG X X.The response ofto low-temperature stress based on metabonomics.Frontiers in Microbiology, 2020, 11: 1787.doi: 10.3389/fmicb.2020.01787.

    [9] ZHANG J T, ZHANG Y, DU Y Y, CHEN S Y, TANG H R.Dynamic metabonomic responses of tobacco () plants to salt stress.Journal of Proteome Research, 2011, 10(4): 1904-1914.doi: 10.1021/pr101140n.

    [10] JEON E, CHOI S, YEO K H, PARK K S, RATHOD M L, LEE J.Development of electrical conductivity measurement technology for key plant physiological information using microneedle sensor.Journal of Micromechanics and Microengineering, 2017, 27(8): 85009.doi: 10.1088/1361-6439/aa7362.

    [11] LU L, HUANG M, HUANG Y X, CORVINI P F X, JI R, ZHAO L J.Mn3O4nanozymes boost endogenous antioxidant metabolites in cucumber () plant and enhance resistance to salinity stress.Environmental Science: Nano, 2020, 7(6): 1692-1703.doi: 10.1039/D0EN00214C.

    [12] HE J, BOUWMEESTER H J, DICKE M, KAPPERS I F.Correction to: Transcriptional and metabolite analysis reveal a shift in direct and indirect defences in response to spider?mite infestation in cucumber ().Plant Molecular Biology, 2020, 103(4/5): 507-509.doi: 10.1007/s11103-020-01009-8.

    [13] LI T F, LU M, GAO Y H, HUANG X D, LIU G Y, XU D H.Double layer MOFs M-ZIF-8@ZIF-67: The adsorption capacity and removal mechanism of fipronil and its metabolites from environmental water and cucumber samples.Journal of Advanced Research, 2020, 24: 159-166.doi: 10.1016/j.jare.2020.03.013.

    [14] 趙靜, 郭勇祥, 林志豪, 陳曉鑫, 黃三文, 廖紅.黃瓜磷效率評價及其基因型差異.廣東農(nóng)業(yè)科學(xué), 2015, 42(19): 27-35, 2.doi: 10.16768/j.issn.1004-874x.2015.19.019.

    ZHAO J, GUO Y X, LIN Z H, CHEN X X, HUANG S W, LIAO H.Evaluation of phosphorus efficiency and its genotypic variations in cucumber.Guangdong Agricultural Sciences, 2015, 42(19): 27-35, 2.doi: 10.16768/j.issn.1004-874x.2015.19.019.(in Chinese)

    [15] 葛體達(dá), 姜武, 宋世威, 黃丹楓.無機(jī)氮和氨基酸態(tài)氮處理對番茄幼苗木質(zhì)部和韌皮部汁液中礦質(zhì)養(yǎng)分的影響.園藝學(xué)報, 2009, 36(3): 347-354.doi: 10.16420/j.issn.0513-353x.2009.03.027.

    GE T D, JIANG W, SONG S W, HUANG D F.Influence of inorganic and amino acid nitrogen on mineral nutrient contents in xylem and phloem sap of different tomato () cultivars.Acta Horticulturae Sinica, 2009, 36(3): 347-354.doi: 10.16420/j.issn.0513-353x.2009.03.027.(in Chinese)

    [16] FAN H F, XU Y L, DU C X, WU X.Phloem sap proteome studied by iTRAQ provides integrated insight into salinity response mechanisms in cucumber plants.Journal of Proteomics, 2015, 125: 54-67.doi: 10.1016/j.jprot.2015.05.001.

    [17] 甄麗莎, 谷潔, 胡婷, 劉晨, 賈鳳安, 呂睿.石油烴類污染物降解動力學(xué)和微生物群落多樣性分析.農(nóng)業(yè)工程學(xué)報, 2015, 31(15): 231-238.

    ZHEN L S, GU J, HU T, LIU C, JIA F A, Lü R.Kinetics of petroleum hydrocarbon degradation in soil and diversity of microbial community during composting.Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(15): 231-238.(in Chinese)

    [18] AGUIRRE M, KIEGLE E, LEO G, EZQUER I.Carbohydrate reserves and seed development: An overview.Plant Reproduction, 2018, 31(3): 263-290.doi: 10.1007/s00497-018-0336-3.

    [19] ZHAO D K, YU Y, SHEN Y, LIU Q, ZHAO Z W, SHARMA R, REITER R J.Melatonin synthesis and function: evolutionary history in animals and plants.Frontiers in Endocrinology, 2019, 10: 249.doi: 10.3389/fendo.2019.00249.

    [20] LóPEZ-BUCIO J, NIETO-JACOBO M F, RAMíREZ-RODRíGUEZ V, HERRERA-ESTRELLA L.Organic acid metabolism in plants: from adaptive physiology to transgenic varieties for cultivation in extreme soils.Plant Science, 2000, 160(1): 1-13.doi: 10.1016/s0168- 9452(00)00347-2.

    [21] 程少禹, 宣鈴娟, 董彬, 顧翠花, 申亞梅, 張明如, 戴夢怡, 王卓為, 章穎佳, 陸丹迎.‘紅元寶’紫玉蘭兩次花芽分化差異代謝通路及關(guān)鍵調(diào)控基因篩選.園藝學(xué)報, 2020, 47(8): 1490-1504.doi: 10.16420/j.issn.0513-353x.2019-0963.

    CHENG S Y, XUAN L J, DONG B, GU C H, SHEN Y M, ZHANG M R, DAI M Y, WANG Z W, ZHANG Y J, LU D Y.Identification of differential metabolic pathways and key regulatory genes in the two flower bud differentiation processes of.Acta Horticulturae Sinica, 2020, 47(8): 1490-1504.doi: 10.16420/j.issn.0513-353x.2019-0963.(in Chinese)

    [22] 滑金杰, 王華杰, 王近近, 李佳, 江用文, 王岳梁, 袁海波.采用PLS-DA分析毛火方式對工夫紅茶品質(zhì)的影響.農(nóng)業(yè)工程學(xué)報, 2020, 36(8): 260-270.doi: 10.11975/j.issn.1002-6819.2020.08.032.

    HUA J J, WANG H J, WANG J J, LI J, JIANG Y W, WANG Y L, YUAN H B.Influences of first-drying methods on the quality of Congou black tea using partial least squares-discrimination analysis.Transactions of the Chinese Society of Agricultural Engineering, 2020, 36(8): 260-270.doi: 10.11975/j.issn.1002-6819.2020.08.032.(in Chinese)

    [23] 劉曉燕, 劉艷秋, 程孟春, 肖紅斌.超高效液相色譜-質(zhì)譜聯(lián)用技術(shù)在藥物肝損傷代謝組學(xué)研究中的應(yīng)用.色譜, 2015, 33(7): 683-690.doi: 10.3724/SP.J.1123.2015.04007.

    LIU X Y, LIU Y Q, CHENG M C, XIAO H B.Application of ultra high performance liquid chromatography-mass spectrometry to metabolomics study of drug-induced hepatotoxicity.Chinese Journal of Chromatography, 2015, 33(7): 683-690.doi: 10.3724/SP.J.1123.2015.04007.(in Chinese)

    [24] LIU Y L, HOU W P, JIN J, CHRISTENSEN M J, GU L J, CHENG C, WANG J F.Epichlo? gansuensis increases the tolerance of achnatherum inebrians to low-P stress by modulating amino acids metabolism and phosphorus utilization efficiency.Journal of Fungi (Basel, Switzerland), 2021, 7(5): 390.

    [25] 于新超, 王晶, 朱美玉, 姜晶.碳水化合物代謝參與番茄響應(yīng)低磷脅迫的分子機(jī)制.分子植物育種, 2015, 13(12): 2833-2842.doi: 10.13271/j.mpb.013.002833.

    YU X C, WANG J, ZHU M Y, JIANG J.Molecular mechanism of carbohydrate metabolism participation in tomato response to low phosphorus stress.Molecular Plant Breeding, 2015, 13(12): 2833-2842.doi: 10.13271/j.mpb.013.002833.(in Chinese)

    [26] GAO H Y, MAO H P, ULLAH I.Analysis of metabolomic changes in lettuce leaves under low nitrogen and phosphorus deficiencies stresses.Agriculture (Basel), 2020, 10(9): 406.

    [27] 潘航, 馮纓, 王喜勇, 李巖.荒漠環(huán)境下10種沙拐棗的生理特征比較研究.草業(yè)學(xué)報, 2017, 26(6): 68-75.

    PAN H, FENG Y, WANG X Y, LI Y.Examination and comparison of the physiological characteristics of tenspecies in a desert environment.Acta Prataculturae Sinica, 2017, 26(6): 68-75.(in Chinese)

    [28] 黃子凌, 李大勇, 宋鳳鳴.植物系統(tǒng)獲得抗性中的系統(tǒng)信號及其作用機(jī)制.植物生理學(xué)報, 2020, 56(7): 1346-1360.doi: 10.13592/j.cnki.ppj.2020.0071.

    HUANG Z L, LI D Y, SONG F M.Systemic signals and their mechanisms in plant systemic acquired resistance.Plant Physiology Journal, 2020, 56(7): 1346-1360.doi: 10.13592/j.cnki.ppj.2020.0071.(in Chinese)

    [29] 郭鵬飛, 雷健, 羅佳佳, 劉攀道, 虞道耿, 羅麗娟.柱花草苯丙氨酸解氨酶(SgPALs)對生物脅迫與非生物脅迫的響應(yīng).熱帶作物學(xué)報, 2019, 40(9): 1742-1751.doi: 10.3969/j.issn.1000-2561.2019.09.011.

    GUO P F, LEI J, LUO J J, LIU P D, YU D G, LUO L J.Response of phenylpropane ammonia-lyase on biotic and abiotic stress in.Chinese Journal of Tropical Crops, 2019, 40(9): 1742-1751.doi: 10.3969/j.issn.1000-2561.2019.09.011.(in Chinese)

    [30] DU Y L, ZHAO Q, CHEN L R, YAO X D, ZHANG W, ZHANG B, XIE F T.Effect of drought stress on sugar metabolism in leaves and roots of soybean seedlings.Plant Physiology and Biochemistry, 2020, 146: 1-12.doi: 10.1016/j.plaphy.2019.11.003.

    [31] DALAL A, RANA J S, KUMAR A.Ultrasensitive nanosensor for detection of malic acid in tomato as fruit ripening indicator.Food Analytical Methods, 2017, 10(11): 3680-3686.doi: 10.1007/s12161- 017-0919-x.

    [32] SANGUBOTLA R, KIM J S.A facile enzymatic approach for selective detection of γ-aminobutyric acid using corn-derived fluorescent carbon dots.Applied Surface Science, 2019, 490: 61-69.

    [33] DUVANOVA O V, KRIVONOSOVA I A, ZYABLOV A N, FALALEEV A V, SELEMENEV V F, SOKOLOVA S A.Use of piezoelectric sensors for the determination of oleic and palmitic acids in vegetable oils.Inorganic Materials, 2018, 54(14): 1387-1391.doi: 10.1134/S002016851814008X.

    [34] SINGH A K, SINGH M.Designing L-serine targeted molecularly imprinted polymer via theoretical investigation.Journal of Theoretical and Computational Chemistry, 2016, 15(5): 1650041.doi: 10.1142/ S0219633616500413.

    [35] LU Y L, HUANG Y X, LI S, ZHANG Q, WU J J, XIONG Z Y, XIONG L X, WAN Q Q, LIU Q J.Fat taste detection with odorant- binding proteins (OBPs) on screen-printed electrodes modified by reduced graphene oxide.Sensors & Actuators: B Chemical, 2017, 252: 973-982.doi: 10.1016/j.snb.2017.06.100.

    Metabolites Changes of Cucumber Xylem and Phloem Sap Under Low Phosphorus Stress

    LI QingLin*, ZHANG WenTao, XU Hui, SUN JingJing

    College of Agricultural Engineering, Jiangsu University, Zhenjiang 212013, Jiangsu

    【Background】 Xylem sap is responsible for the transport of water and nutrients for the aerial part, while phloem sap is mainly responsible for the transport of photosynthetic products.When plants are subjected to abiotic stress, the sap of xylem and phloem undergoes different changes.【Objective】 The aim of this study was to investigate the content changes of various metabolites in cucumber xylem and phloem sap under low phosphorus stress, and to explore the metabolic changes of cucumber under low phosphorus stress, so as to provide new ideas and basis for the diagnosis of cucumber phosphorus deficiency.【Method】 The xylem and phloem sap were collected from cucumber plants under normal phosphorus level (1.25 mmol?L-1) and low phosphorus stress (0.3125 mmol?L-1).These sap samples were analyzed qualitatively and quantitatively by GC-MS after trimethylsilyl derivatization.【Result】In xylem sap under low phosphorus stress, there was no significant change in carbohydrates, and the content of most amino acids increased significantly.Among organic acids, the malic acid showed a significant downward trend, while others, such as palmitic acid and oleic acid, were significantly higher than the normal levels.In phloem sap under low phosphorus stress, most carbohydrates and amino acids were lower than the normal level, and the contents of organic acids, such as oxalic acid and succinic acid, were also significantly reduced.Further, in order to screen potential biomarkers in cucumber xylem and phloem sap, the detection data were further analyzed by OPLS-DA and Student's t test.Eight potential biomarkers (VIP>1,<0.05) under low phosphorus stress were detected in cucumber xylem sap, which were cadaverine, malic acid, 4-aminobutyric acid, palmitic acid, myo-inositol, cellobiose acid, N, N-dimethyldodecamide, and glycerol monostearate.Seven potential biomarkers (VIP>1,<0.05) under low phosphorus stress were detected in cucumber phloem sap, which were serine, linoleic acid, 9-(3-methyl-5-pentylfuran-2-yl) nonanoic acid, dehydroabietic acid, N, N-dimethyldodecamide, 2, 2'-methylenebis (6-tert-butyl-4-methylphenol), and L-2-piperic acid.【Conclusion】 Under low phosphorus stress, the carbohydrates, amino acids and organic acids in the xylem and phloem sap of cucumber changed, but the change trends were not the same.Eight and seven biomarkers were screened from xylem and phloem sap, which could provide new ideas and basis for the diagnosis of cucumber phosphorus deficiency.Among the biomarkers, the detection methods of malic acid, 4-aminobutyric acid, palmitic acid, serine and linoleic acid have been relatively mature, and the cucumber phosphorus deficiency diagnosis method based on them was more practical.

    cucumber; xylem; phloem; phosphorus stress; metabolomics; biomarkers

    2021-08-13;

    2021-10-13

    國家自然科學(xué)基金(62101217)、江蘇大學(xué)農(nóng)業(yè)裝備學(xué)部項目(NZXB20200203)、江蘇省自然科學(xué)基金(BK20190858)、中國博士后基金(2019M661754)

    通信作者李青林,E-mail:776767151@qq.com

    (責(zé)任編輯 趙伶俐)

    猜你喜歡
    韌皮部木質(zhì)部汁液
    為什么蟬會撒尿噴人?
    不同品種吊蘭根初生木質(zhì)部原型的觀察與比較
    植物研究(2021年2期)2021-02-26 08:40:10
    霜花:冰凍的植物汁液
    給麒麟掌給麒麟掌“打針”
    鎘污染來源對蘿卜鎘積累特性的影響
    蘋果樹木質(zhì)部及韌皮部組織基因組DNA的提取及質(zhì)量檢測
    木薯塊根膨大期韌皮部和木質(zhì)部比較蛋白組學(xué)初步研究
    鹽分脅迫對2種楊樹次生木質(zhì)部導(dǎo)管特征的影響
    提取庫爾勒香梨葉片質(zhì)外體汁液離心力的確認(rèn)
    秋季斷根有利于渭北黃土高原蘋果安全越冬
    亚洲欧洲日产国产| 狠狠狠狠99中文字幕| 国产黄片视频在线免费观看| 精品国产露脸久久av麻豆 | 日韩av不卡免费在线播放| 国产精品,欧美在线| 日本爱情动作片www.在线观看| 99在线视频只有这里精品首页| 99久久成人亚洲精品观看| 亚洲最大成人中文| 亚洲精品日韩av片在线观看| 青青草视频在线视频观看| 午夜爱爱视频在线播放| 免费人成在线观看视频色| 亚洲电影在线观看av| 九九久久精品国产亚洲av麻豆| 久久综合国产亚洲精品| 色视频www国产| 麻豆一二三区av精品| 亚洲av日韩在线播放| 色吧在线观看| 亚洲精品国产av成人精品| 久久久久久久国产电影| 精品酒店卫生间| 国产精品福利在线免费观看| 男女国产视频网站| 少妇熟女欧美另类| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在| 久久精品国产自在天天线| 少妇裸体淫交视频免费看高清| 少妇人妻精品综合一区二区| 18+在线观看网站| 中文天堂在线官网| av又黄又爽大尺度在线免费看 | 亚洲欧美日韩东京热| 人人妻人人澡欧美一区二区| 国产精品女同一区二区软件| 久久久国产成人精品二区| 18+在线观看网站| 亚洲av中文av极速乱| 国产一区二区在线av高清观看| 国产精品熟女久久久久浪| 国产又黄又爽又无遮挡在线| 欧美三级亚洲精品| 国产高清视频在线观看网站| 亚洲国产精品成人久久小说| 亚洲av中文字字幕乱码综合| 欧美一区二区精品小视频在线| 亚洲国产欧洲综合997久久,| 国产精品一区www在线观看| 午夜精品一区二区三区免费看| 国产精品久久久久久久电影| 亚洲18禁久久av| 免费观看的影片在线观看| 国产又黄又爽又无遮挡在线| av在线蜜桃| 亚洲av成人av| 国产黄色视频一区二区在线观看 | 日日干狠狠操夜夜爽| 一级av片app| 亚洲在久久综合| 亚洲欧美精品专区久久| 欧美三级亚洲精品| 久久这里只有精品中国| 日韩高清综合在线| 国产成人精品一,二区| 精品熟女少妇av免费看| 亚洲精品乱久久久久久| 欧美色视频一区免费| 成人美女网站在线观看视频| 欧美精品一区二区大全| 亚洲最大成人av| 久久人人爽人人片av| 国产又色又爽无遮挡免| 欧美又色又爽又黄视频| 观看美女的网站| 国产精品伦人一区二区| 三级经典国产精品| 亚洲怡红院男人天堂| 嫩草影院入口| 十八禁国产超污无遮挡网站| 一区二区三区高清视频在线| 精品久久久久久久久av| 国产乱来视频区| 亚洲国产日韩欧美精品在线观看| 国产 一区精品| 最近中文字幕高清免费大全6| 婷婷色av中文字幕| 亚洲色图av天堂| 搡女人真爽免费视频火全软件| 一区二区三区免费毛片| 久久久久久久亚洲中文字幕| 欧美精品一区二区大全| 偷拍熟女少妇极品色| 干丝袜人妻中文字幕| 国产免费男女视频| 亚洲国产精品专区欧美| 好男人在线观看高清免费视频| 欧美97在线视频| 国产爱豆传媒在线观看| 美女xxoo啪啪120秒动态图| 亚洲一区高清亚洲精品| 国产精品一区二区在线观看99 | 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 国产老妇女一区| 亚洲伊人久久精品综合 | 男人的好看免费观看在线视频| 日韩三级伦理在线观看| 亚洲精品国产成人久久av| 国产精品嫩草影院av在线观看| 极品教师在线视频| 中文字幕久久专区| 成人毛片60女人毛片免费| 久久久久久久午夜电影| 建设人人有责人人尽责人人享有的 | 国产伦在线观看视频一区| 日本午夜av视频| 男女视频在线观看网站免费| 精品人妻一区二区三区麻豆| 纵有疾风起免费观看全集完整版 | 少妇猛男粗大的猛烈进出视频 | 成年女人永久免费观看视频| 亚洲人与动物交配视频| 特级一级黄色大片| 国产精品一及| 成人性生交大片免费视频hd| 日韩国内少妇激情av| 亚洲欧洲国产日韩| 国产淫片久久久久久久久| 国产午夜精品一二区理论片| 亚洲欧美精品专区久久| 插阴视频在线观看视频| 18+在线观看网站| 婷婷色综合大香蕉| 国产精华一区二区三区| 中文在线观看免费www的网站| 国产淫语在线视频| 午夜激情福利司机影院| 晚上一个人看的免费电影| 国产精品蜜桃在线观看| 午夜免费男女啪啪视频观看| 国产午夜福利久久久久久| 国产精品久久久久久精品电影| 一夜夜www| 中文字幕精品亚洲无线码一区| 日韩视频在线欧美| 国产熟女欧美一区二区| 国产精品久久久久久精品电影| 最近中文字幕高清免费大全6| 欧美97在线视频| 秋霞伦理黄片| 亚洲av二区三区四区| av又黄又爽大尺度在线免费看 | 久久韩国三级中文字幕| 在线免费观看不下载黄p国产| 国产在线一区二区三区精 | 欧美日韩在线观看h| 国产av不卡久久| 免费观看性生交大片5| 一本一本综合久久| 国产高清三级在线| 国产精品一区二区性色av| 淫秽高清视频在线观看| 在线播放国产精品三级| 日韩精品青青久久久久久| 国产黄片视频在线免费观看| 国产探花在线观看一区二区| 亚洲自拍偷在线| 美女内射精品一级片tv| 波多野结衣巨乳人妻| 26uuu在线亚洲综合色| 男插女下体视频免费在线播放| 中文字幕av在线有码专区| 日本午夜av视频| 国产精华一区二区三区| 日韩国内少妇激情av| 最近最新中文字幕免费大全7| 国产精品一区www在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 又爽又黄无遮挡网站| 精品国产露脸久久av麻豆 | 一个人观看的视频www高清免费观看| 亚洲av成人精品一区久久| 只有这里有精品99| 男人和女人高潮做爰伦理| 级片在线观看| 亚洲乱码一区二区免费版| 好男人视频免费观看在线| 一级毛片aaaaaa免费看小| 91av网一区二区| 插阴视频在线观看视频| 亚洲在久久综合| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美精品综合久久99| 欧美日韩一区二区视频在线观看视频在线 | 日韩在线高清观看一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩强制内射视频| 亚洲人成网站高清观看| 少妇人妻一区二区三区视频| 国产老妇伦熟女老妇高清| 麻豆成人av视频| 日韩欧美三级三区| 午夜精品国产一区二区电影 | 看片在线看免费视频| 日韩中字成人| 日韩 亚洲 欧美在线| 黄色一级大片看看| 欧美性猛交黑人性爽| 久久久久久久国产电影| 2022亚洲国产成人精品| 男人和女人高潮做爰伦理| 国产探花极品一区二区| 精品久久久久久久久av| 亚洲天堂国产精品一区在线| 村上凉子中文字幕在线| 精品人妻熟女av久视频| 熟女电影av网| 69av精品久久久久久| 免费人成在线观看视频色| 久久精品国产自在天天线| 成人性生交大片免费视频hd| 国产白丝娇喘喷水9色精品| 久久久精品94久久精品| 国产成人a区在线观看| 亚洲成人中文字幕在线播放| 欧美xxxx性猛交bbbb| 在线观看66精品国产| 亚洲av电影在线观看一区二区三区 | 男女视频在线观看网站免费| 国产伦在线观看视频一区| 亚洲精品色激情综合| 久久草成人影院| 国产高潮美女av| 国产一区二区三区av在线| 日本欧美国产在线视频| 女人久久www免费人成看片 | 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 日韩在线高清观看一区二区三区| 国产av不卡久久| 欧美高清性xxxxhd video| 91午夜精品亚洲一区二区三区| 成人无遮挡网站| 国产亚洲精品av在线| 干丝袜人妻中文字幕| 亚洲av日韩在线播放| 亚洲aⅴ乱码一区二区在线播放| 看片在线看免费视频| 色综合亚洲欧美另类图片| 国产乱来视频区| 亚洲中文字幕一区二区三区有码在线看| 两个人视频免费观看高清| 国产激情偷乱视频一区二区| 亚洲18禁久久av| 99热这里只有精品一区| 男人舔奶头视频| 插阴视频在线观看视频| 亚洲经典国产精华液单| 91在线精品国自产拍蜜月| 一级毛片aaaaaa免费看小| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类| 成年女人看的毛片在线观看| 亚洲精品影视一区二区三区av| 亚洲成av人片在线播放无| 日韩大片免费观看网站 | 只有这里有精品99| 午夜激情福利司机影院| 国产爱豆传媒在线观看| 久久久久网色| 亚洲国产高清在线一区二区三| 久久99热这里只频精品6学生 | 亚洲精品乱久久久久久| 高清在线视频一区二区三区 | 免费黄色在线免费观看| 1000部很黄的大片| 亚洲最大成人手机在线| 欧美日韩一区二区视频在线观看视频在线 | 国产av码专区亚洲av| 国产高清有码在线观看视频| 欧美+日韩+精品| 国产黄a三级三级三级人| 中文字幕人妻熟人妻熟丝袜美| 美女大奶头视频| 最近2019中文字幕mv第一页| 午夜福利在线观看免费完整高清在| 91精品伊人久久大香线蕉| 久久亚洲国产成人精品v| 日本猛色少妇xxxxx猛交久久| 直男gayav资源| 欧美日韩在线观看h| 最近2019中文字幕mv第一页| 成人亚洲欧美一区二区av| 丰满乱子伦码专区| 99国产精品一区二区蜜桃av| 天美传媒精品一区二区| 成人鲁丝片一二三区免费| 欧美激情国产日韩精品一区| 国产一区二区三区av在线| 91久久精品国产一区二区成人| 晚上一个人看的免费电影| 精品不卡国产一区二区三区| 黑人高潮一二区| 亚洲国产精品成人久久小说| 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 观看美女的网站| 蜜桃亚洲精品一区二区三区| 国产黄片美女视频| 久久久成人免费电影| 99久久精品国产国产毛片| 啦啦啦韩国在线观看视频| 中文字幕av成人在线电影| 国产黄片视频在线免费观看| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 熟女人妻精品中文字幕| 天堂中文最新版在线下载 | 亚洲精品日韩在线中文字幕| 成人性生交大片免费视频hd| 国产高潮美女av| 春色校园在线视频观看| 欧美区成人在线视频| 一二三四中文在线观看免费高清| 亚洲最大成人中文| 狠狠狠狠99中文字幕| 欧美人与善性xxx| 色综合站精品国产| kizo精华| 综合色av麻豆| 99久久精品一区二区三区| 黑人高潮一二区| 日韩欧美国产在线观看| 亚洲精品乱久久久久久| 夜夜爽夜夜爽视频| 国产男人的电影天堂91| 成人欧美大片| 淫秽高清视频在线观看| 久久久久免费精品人妻一区二区| 少妇裸体淫交视频免费看高清| 日韩成人伦理影院| 国产精品精品国产色婷婷| 欧美变态另类bdsm刘玥| 男女国产视频网站| 尾随美女入室| 日本爱情动作片www.在线观看| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 免费一级毛片在线播放高清视频| 男女边吃奶边做爰视频| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 赤兔流量卡办理| 亚洲内射少妇av| 日本一本二区三区精品| 高清av免费在线| 青春草视频在线免费观看| 黑人高潮一二区| 成人美女网站在线观看视频| 亚洲内射少妇av| 欧美97在线视频| 精品久久久久久久久亚洲| 99久久精品国产国产毛片| 国产精品久久久久久久电影| 女人十人毛片免费观看3o分钟| 久久久精品大字幕| 日韩欧美三级三区| 精品国内亚洲2022精品成人| 亚洲精品色激情综合| 欧美成人精品欧美一级黄| 久久热精品热| 亚洲av一区综合| 高清毛片免费看| 色网站视频免费| 18禁在线播放成人免费| 国产色爽女视频免费观看| 日韩一本色道免费dvd| 日本免费a在线| 观看免费一级毛片| 2021天堂中文幕一二区在线观| 性色avwww在线观看| 91狼人影院| 色吧在线观看| 久久精品国产99精品国产亚洲性色| 永久网站在线| 日日干狠狠操夜夜爽| 美女被艹到高潮喷水动态| 国产精品av视频在线免费观看| 久久人人爽人人片av| 欧美三级亚洲精品| av在线播放精品| 波多野结衣巨乳人妻| 色吧在线观看| 女的被弄到高潮叫床怎么办| 久久久午夜欧美精品| 麻豆乱淫一区二区| 国产 一区精品| 成人亚洲精品av一区二区| 最新中文字幕久久久久| 日韩欧美国产在线观看| 国产av不卡久久| 少妇的逼水好多| 嫩草影院精品99| 熟妇人妻久久中文字幕3abv| 免费搜索国产男女视频| 精品国产一区二区三区久久久樱花 | 成人亚洲精品av一区二区| 精品国产三级普通话版| 国产 一区 欧美 日韩| 国产精品熟女久久久久浪| 国产成人精品久久久久久| 中文乱码字字幕精品一区二区三区 | 免费看a级黄色片| 亚洲精品456在线播放app| 在线观看一区二区三区| 国产大屁股一区二区在线视频| 久久久久久久久久黄片| 成人一区二区视频在线观看| 中文字幕制服av| 天堂影院成人在线观看| 国产极品天堂在线| 国产精华一区二区三区| 国内精品一区二区在线观看| 免费看美女性在线毛片视频| 久久久久久久久久久免费av| 最近手机中文字幕大全| 欧美精品一区二区大全| 干丝袜人妻中文字幕| 久久久久精品久久久久真实原创| 国产大屁股一区二区在线视频| 国产男人的电影天堂91| 久久亚洲国产成人精品v| 一个人免费在线观看电影| 亚洲最大成人中文| 最后的刺客免费高清国语| 三级国产精品欧美在线观看| 免费在线观看成人毛片| 国产一区二区在线观看日韩| videos熟女内射| 长腿黑丝高跟| 亚洲成av人片在线播放无| 69人妻影院| 在线a可以看的网站| 99久国产av精品国产电影| 午夜福利高清视频| 亚洲一区高清亚洲精品| 国产不卡一卡二| 日本三级黄在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人漫画全彩无遮挡| 搞女人的毛片| 日韩高清综合在线| 91午夜精品亚洲一区二区三区| 亚洲一区高清亚洲精品| 国产精品久久久久久精品电影| 成人综合一区亚洲| 中国国产av一级| 色综合色国产| 精品欧美国产一区二区三| 久久久午夜欧美精品| 亚洲av电影不卡..在线观看| 在线免费观看的www视频| 成人欧美大片| 观看美女的网站| 色网站视频免费| 国产精品一区www在线观看| 亚洲自拍偷在线| 别揉我奶头 嗯啊视频| 国产男人的电影天堂91| 国产伦一二天堂av在线观看| 国产精品久久视频播放| 欧美精品一区二区大全| 亚洲欧美精品综合久久99| 国产成人免费观看mmmm| 久久久久久久久中文| av视频在线观看入口| av.在线天堂| 午夜激情欧美在线| 久久亚洲精品不卡| 成人一区二区视频在线观看| 狠狠狠狠99中文字幕| 亚洲国产高清在线一区二区三| 最近的中文字幕免费完整| 精品久久久久久久久亚洲| 男女那种视频在线观看| 嫩草影院精品99| 亚洲av二区三区四区| 简卡轻食公司| 美女大奶头视频| av在线老鸭窝| 欧美一区二区亚洲| 国产精品久久久久久久电影| 国产乱来视频区| www日本黄色视频网| 欧美另类亚洲清纯唯美| 少妇丰满av| 国产精品国产三级国产av玫瑰| 少妇人妻精品综合一区二区| 永久免费av网站大全| 国产精品久久久久久av不卡| 人妻制服诱惑在线中文字幕| 少妇裸体淫交视频免费看高清| 内射极品少妇av片p| av在线天堂中文字幕| 日韩欧美在线乱码| 国产亚洲午夜精品一区二区久久 | 成人毛片60女人毛片免费| 天天躁日日操中文字幕| 精品国内亚洲2022精品成人| 日韩制服骚丝袜av| 内地一区二区视频在线| 亚洲国产精品sss在线观看| 日韩一区二区三区影片| 久久热精品热| 麻豆成人午夜福利视频| 国内精品美女久久久久久| 美女cb高潮喷水在线观看| 日本免费在线观看一区| 能在线免费看毛片的网站| 人体艺术视频欧美日本| 六月丁香七月| 三级经典国产精品| 日日干狠狠操夜夜爽| 免费av观看视频| 日韩三级伦理在线观看| 中文在线观看免费www的网站| 日韩欧美精品v在线| 国产黄色视频一区二区在线观看 | 乱人视频在线观看| av黄色大香蕉| 亚洲av成人精品一二三区| 国产精品一二三区在线看| 日本五十路高清| 亚洲一区高清亚洲精品| 日本色播在线视频| 成年av动漫网址| 搡女人真爽免费视频火全软件| 亚洲美女搞黄在线观看| 国产精品乱码一区二三区的特点| 亚洲自偷自拍三级| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 国产激情偷乱视频一区二区| 国产精品人妻久久久影院| 国产在线一区二区三区精 | 国产精品爽爽va在线观看网站| 久久这里有精品视频免费| 亚洲国产欧洲综合997久久,| 天堂影院成人在线观看| 久久婷婷人人爽人人干人人爱| 男人狂女人下面高潮的视频| 男女啪啪激烈高潮av片| 国产精品一区www在线观看| 亚洲欧美精品自产自拍| 18禁在线无遮挡免费观看视频| 美女脱内裤让男人舔精品视频| 欧美xxxx黑人xx丫x性爽| 日本免费在线观看一区| 好男人在线观看高清免费视频| 久久久亚洲精品成人影院| 久久精品国产亚洲网站| 日本五十路高清| 亚洲av不卡在线观看| 国产精品久久久久久精品电影| 国产高清三级在线| 免费电影在线观看免费观看| 一级av片app| 久久99精品国语久久久| 在线免费观看不下载黄p国产| 婷婷色麻豆天堂久久 | 免费无遮挡裸体视频| 国产免费视频播放在线视频 | 国产成人一区二区在线| 久久99热这里只频精品6学生 | 免费av毛片视频| 精品久久久噜噜| 噜噜噜噜噜久久久久久91| 久久久久久久久久成人| 欧美激情在线99| 三级经典国产精品| 亚洲综合精品二区| 久久精品91蜜桃| 毛片女人毛片| 欧美另类亚洲清纯唯美| 午夜日本视频在线| 久久人人爽人人片av| 黄片无遮挡物在线观看| 国产乱来视频区| 亚洲欧美精品专区久久| 桃色一区二区三区在线观看| 日韩一本色道免费dvd| 国产人妻一区二区三区在| 国内少妇人妻偷人精品xxx网站| 久久精品夜色国产| 国产精品无大码| 别揉我奶头 嗯啊视频| 亚洲熟妇中文字幕五十中出| 国产伦精品一区二区三区四那| 亚洲最大成人av| 熟妇人妻久久中文字幕3abv| 综合色丁香网| 亚洲在久久综合| 日日摸夜夜添夜夜爱| 久久久久久久久久久免费av| 一本久久精品| 大话2 男鬼变身卡| 人妻少妇偷人精品九色| 午夜福利在线在线| 成人鲁丝片一二三区免费| 3wmmmm亚洲av在线观看| 老师上课跳d突然被开到最大视频|