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

    不同耐旱性青稞葉片差異蛋白分析

    2021-04-26 12:08:56姚曉華吳昆侖
    作物學報 2021年7期
    關鍵詞:旱地大麻青稞

    李 潔 付 惠 姚曉華 吳昆侖

    不同耐旱性青稞葉片差異蛋白分析

    李 潔**付 惠**姚曉華 吳昆侖*

    青海大學農林科學院 / 青海省農林科學院 / 青海省青稞遺傳育種重點實驗室 / 國家麥類改良中心青海青稞分中心, 青海西寧 810016

    為從蛋白質水平揭示不同青稞品種響應干旱脅迫的差異, 分析抗旱蛋白質分子機制, 本研究以干旱脅迫不敏感的旱地紫青稞(HDZ)和干旱脅迫敏感的大麻青稞(DM)為研究材料, 以盆栽限量供水種植方法, 對干旱處理不同梯度的青稞葉片進行葉綠素、可溶性蛋白、丙二醛含量及相對電導率4項生理指標的測定, 同時利用iTRAQ技術對深度干旱脅迫青稞葉片全蛋白組進行差異蛋白分析。結果表明: 隨著干旱處理時間的延長, 兩種青稞的葉綠素和可溶性蛋白含量逐漸下降, 電導率及丙二醛的含量逐漸升高, 在相同處理條件下, 大麻青稞的葉綠素和可溶性蛋白含量降低幅度、電導率及丙二醛含量的增高幅度明顯高于旱地紫青稞; 對兩個青稞品種的干旱脅迫和正常培養(yǎng)的比較組進行iTRAQ分析, 共定量出4163個蛋白(多肽), 其中旱地紫青稞比較組中對比正常培養(yǎng), 篩選到表達上調的蛋白68個, 下調的蛋白63個, 在大麻青稞的比較組中篩選出表達上調蛋白21個, 下調蛋白32個。KEGG通路分析表明, 富集程度位于前3位的通路是代謝、氨基酸的生物合成以及次級代謝產物的生物合成, 主要涉及到檸檬酸循環(huán)、碳循環(huán)等代謝通路, 丙氨酸、精氨酸等氨基酸的合成降解以及花生四烯酸、亞麻酸等次級代謝產物的合成。本研究從蛋白組水平篩選了青稞干旱脅迫響應相關的代謝通路和其他相關功能的蛋白, 為揭示青稞干旱脅迫的應答分子調控機制提供了一定的理論基礎。

    青稞; 干旱脅迫; 葉片; 差異蛋白; iTRAQ

    干旱是影響面積最廣、危害程度最深的非生物脅迫, 每年因干旱導致作物減產高達50%以上[1]。植物為適應干旱環(huán)境, 其葉片、根系等組織感受到脅迫后, 會通過信號轉導級聯反應, 轉錄激活、調節(jié)大量基因及蛋白質的轉錄與表達, 從而協(xié)助植物得以在逆境中存活[2-4]。因此, 絕大部分基因及其功能有待于在蛋白質層面予以揭示和闡述。

    青稞(L. var.Hook. f.)屬于禾本科小麥族大麥屬普通大麥種, 因其內外穎殼分離, 籽粒裸露, 又稱為裸大麥[5], 是青藏高原等高寒地區(qū)的主要作物。青稞主要種植在高寒缺氧、環(huán)境惡劣的地區(qū), 在長期的適應性進化過程中產生了大量的優(yōu)異抗逆基因, 具有耐寒、耐旱、耐瘠薄等優(yōu)異種性[6]。在青海, 青稞種植區(qū)域多在海拔2700 m以上灌溉設施不能滿足的高寒農區(qū)和農牧交錯區(qū), 其受旱時期主要發(fā)生在春季及夏季, 嚴重影響青稞幼苗分蘗的形成與籽粒灌漿, 最終導致籽粒不飽滿以及減產等問題。因此選擇耐旱青稞品種, 可以使廣闊、干旱、貧瘠的土地最大資源化[5,7]。

    目前對青稞抗逆基因[8-10]的研究已經有了一些報道, 但利用差異蛋白質研究方法篩選、鑒定相關耐旱基因在青稞抗逆研究中應用鮮見報道。本研究采用同位素標記相對和絕對定量(isobaric tags for relative and absolute quantification, iTRAQ)技術, 比較兩種對水分脅迫敏感程度不同的青稞在干旱脅迫下葉片蛋白表達的差異, 進行定性定量分析, 旨在從蛋白水平揭示青稞品種間的耐旱差異, 為青稞耐旱研究提供一些方法和手段。

    1 材料與方法

    1.1 試驗材料

    根據之前28個青稞品種干旱脅迫篩選結果[5], 從中選用敏感與不敏感的2個青稞品種為供試青稞品種, 材料分別為旱地紫青稞(不敏感, HDZ)、大麻青稞(敏感, DM)。

    1.2 試驗處理

    材料采用玻璃溫室盆栽限量供水種植, 每個花盆中放風干土6 kg, 種植前兩天用水將花盆中的土壤澆透。選擇籽粒光滑飽滿的種子播種, 每盆播25粒, 出芽后每盆定苗20株, 溫室溫度22~25℃。待青稞植株長至三葉一心時處理, 設干旱脅迫與正常供水兩種處理。干旱梯度形成前確保各盆土壤相對含水量基本達到飽和, 然后停止?jié)菜屍渥匀桓珊抵猎O定標準, 每處理重復3次。在干旱處理第7天(土壤相對含水量為對照的65%)及第10天(土壤相對含水量為對照的45%)上午10:00, 收取長勢一致的植株倒二葉同一部位葉片約0.5 g, 放置于–80℃超低溫冰箱保存。

    1.3 生理指標的測量

    在干旱處理第7天和第10天, 采植株葉片進行葉綠素含量、蛋白質含量、相對電導率及丙二醛含量等生理生化指標進行測定。葉綠素含量采用丙酮乙醇混合液法測定[11]; 蛋白質含量采用Bradford法測定[12]; 電導率采用電導率儀法測定, 以相對電導率(%)表示, 相對電導率= (浸泡電導率值/煮沸后電導率值)×100%, 丙二醛(MDA)含量采用硫代巴比妥酸顯色法測定[13]。

    1.4 葉片總蛋白質提取

    采用TCA-丙酮沉淀-酚/SDS-聯合抽提法提取青稞葉片總蛋白[14]。提取樣品為干旱處理第10天2個青稞的葉片及正常培養(yǎng)對照葉片??偟鞍滋崛r, 預冷研磨器具及其他需要接觸到樣品的工具; 樣品加入液氮研磨, 加入含0.2% DTT的10% TCA-丙酮(w/v) (–20℃預冷), 充分混合, –20℃沉淀2 h, 20,000×, 4℃, 離心30 min; 棄上清, 取沉淀; 沉淀用含100 mmol L–1醋酸銨的80%甲醇溶液洗1次, 80%丙酮洗1次, 通風櫥內室溫干燥10 min; 加入10 mL酚/SDS抽提液振蕩混勻溫育5 min, 4℃, 20,000×離心10 min, 轉移上層酚相至干凈新管中, 加入5倍體積的預冷的含100 mmol L–1醋酸銨的甲醇溶液, –20℃沉淀4 h或過夜, 20,000×離心10 min, 棄上清液, 沉淀中再加入預冷的80%甲醇和80%丙酮各一次, 懸起沉淀并充分清洗后將溶液吸出, 通風櫥內自然風干, –80℃冰箱保存?zhèn)溆谩?0 mL酚/SDS抽提液: pH 8.0 Tris飽和酚與SDS緩沖液(30%蔗糖、2% SDS、5% β-巰基乙醇、0.1 mol L–1pH 8.0 Tris-HCl)等量混合。

    1.5 蛋白質裂解與定量

    采用Bradford[12]法對蛋白進行定量。稱取20 mg蛋白質干粉加入200 μL裂解緩沖液, 30℃水浴1 h, 14℃、20,000×離心30 min, 取上清液即為葉片總蛋白溶液, 測定蛋白濃度, 每個樣品測量3個生物學重復。

    1.6 蛋白組學測定

    按照美國AB Sciex公司生產的iTRAQ Reagent 8 Plex Multi-plex試劑盒說明書, 對2組共計12個蛋白質樣品進行酶解與標記。混合標記后的樣品用gemini NX (Phenomemex, 00F-4453-E0)色譜柱進行分級, 隨后用1100系列高效液相色譜儀(Agilent公司, 美國)和Orbi-trap Elite質譜儀(Thero Scientific公司, 美國)進行液相分離和質譜分析[15-16], 此部分工作由上海維基生物科技有限公司承擔。

    1.7 數據處理與分析

    1.7.1 原始數據 根據鑒定蛋白質的FASTA格式序列文本, 利用Compute pI/Mw工具計算每一個蛋白質的分子量和等電點, 將定量及鑒定結果進行合并處理, 采用GO數據庫的注釋對鑒定蛋白質的功能進行分類。MS/MS數據通過搜索NCBI和EST數據庫。搜索條件: 胰酶消化, 最大遺漏酶切位點1個, 肽段質量精確度1×10–3%, MS/MS 質量精確度±0.50 Da, 固定修飾Carbamidomethyl (C)。根據蛋白質譜搜庫檢索和差異蛋白篩選結果, 進行Blast分析。

    1.7.2 GO分析及KEGG代謝通路分析 所有鑒定差異表達蛋白通過ID號到Unipro數據庫中批量獲取對應的注釋文件(*.dat) FASTA文件, 在Customize選項中選Geneontology鏈接到GO數據庫, 得到每個蛋白的GO注釋信息。將差異表達蛋白的FASTA格式文件上傳到KOBAS2.0中, 選擇map到KO中, 得到差異表達蛋白所對應的KO號, 根據所得到的KO號在KEGG ORTHOLOGY數據庫中做代謝通路分析, 得到差異表達蛋白質作圖到各個通路上的結果。

    2 結果與分析

    2.1 葉綠素含量

    干旱脅迫處理對2個青稞品種的葉綠素含量產生了不同的影響(圖1)。干旱脅迫后的2個品種的葉片葉綠素、葉綠素、總葉綠素含量與對照相比呈逐漸降低的趨勢, 干旱脅迫時間越久, 葉綠素含量對比對照降低越多, 其中大麻青稞葉片葉綠素、葉綠素、總葉綠素含量的降低趨勢更為明顯, 說明在重度干旱脅迫時大麻青稞葉片的光合能力受到的抑制明顯高于旱地紫青稞。

    A: 葉綠素; B: 葉綠素; C: 總葉綠素。大寫字母表示0.01極顯著水平, 小寫字母表示0.05顯著水平。

    A: chlorophyll; B: chlorophyll; C: total chlorophyll. Uppercase letters indicate extremely significant differences at the 0.01 proba-bility level, lowercase letters indicate significant differences at the 0.05 probability level.

    2.2 可溶性蛋白質含量

    干旱脅迫條件下, 旱地紫青稞和大麻青稞葉片可溶性蛋白質含量隨干旱脅迫時間的延長呈逐漸降低的趨勢(表1), 干旱脅迫第7天旱地紫青稞葉片可溶性蛋白質含量降低為對照的94.38%, 干旱脅迫第10天, 降低為對照的92.47%。大麻青稞葉片在干旱脅迫第7天可溶性蛋白含量降低為對照的78.70%, 干旱脅迫第10天, 降低為對照的76.50%, 降低趨勢明顯大于旱地紫青稞。推測可溶性蛋白與植物細胞的滲透調節(jié)有關, 高含量的可溶性蛋白可使細胞維持較低的滲透勢, 抵抗水分脅迫帶來的傷害, 但是當水分虧缺嚴重到影響了植物蛋白合成代謝時, 可溶性蛋白含量就會下降。由此可以看出, 旱地紫青稞比大麻青稞具有更強的抗?jié)B透脅迫能力。

    2.3 相對電導率和丙二醛含量

    干旱脅迫處理后, 旱地紫青稞葉片及大麻青稞葉片相對電導率和丙二醛含量均呈現增高趨勢(表1)。但在相同干旱脅迫下大麻青稞葉片的相對電導率和丙二醛含量增加趨勢更加明顯, 證明其細胞膜遭到破壞更大, 導致膜透性增大, 電解質外滲更多, 同時膜脂過氧化程度更高, 釋放出的丙二醛還會抑制蛋白合成, 從而增加細胞的滲透勢。這些結果表明, 大麻青稞膜系統(tǒng)對干旱脅迫的耐受程度明顯低于旱地紫青稞。

    2.4 蛋白鑒定

    分別對干旱處理10 d (HDZD)和正常水分培養(yǎng)(HDZCK)的旱地紫青稞葉片及干旱處理10 d (DMD)和正常水分培養(yǎng)(DMCK)的大麻青稞葉片, 采用iTRAQ技術進行蛋白質的定性定量工作, 共鑒定出蛋白質及多肽4163個。由表2可知, 肽段覆蓋率在10%以下的蛋白占鑒定總蛋白的30.55%, 肽段覆蓋率在10%~20%的蛋白占鑒定總蛋白的37.45%, 肽段覆蓋率在20%~30%以上的蛋白占鑒定總蛋白的19.19%, 肽段覆蓋率在30%~40%以上的蛋白占鑒定總蛋白的8.82%, 肽段覆蓋率在40%以上的蛋白占鑒定總蛋白的4.01%, HDZD- HDZCK比較組葉片經過質譜分析共測出蛋白質及多肽共2150個, DMD- DMCK比較共測出蛋白質及多肽2013個。

    表1 干旱處理對青稞葉片蛋白質、相對電導率及丙二醛含量的影響

    大寫字母表示0.01極顯著水平, 小寫字母表示0.05顯著水平。

    Uppercase letters indicate extremely significant differences at the 0.01 probability level; lowercase letters indicate significant differences at the 0.05 probability level.

    表2 鑒定蛋白的覆蓋率分布

    2.5 差異表達蛋白

    根據實驗條件篩選并注釋了4163個差異蛋白, 對每個蛋白質的差異倍數以2為底數取其對數, 以此作為橫坐標, 用值以10為底取其對數的絕對值為縱坐標, 做出差異蛋白的火山圖(圖2)。結果顯示, 2個品種青稞的差異表達蛋白總數是184個, 其中, HDZD-HDZCK比較組中表達上調的蛋白68個, 表達下調的蛋白63個, 共131個。蛋白含量相差4.0倍以上的有4個, 其中最大為7.495倍。含量相差3.0倍以上的有11個, 相差3.0~2.0倍的有17個, 相差1.9~1.0倍有31個。在DMD-DMCK比較組中表達上調蛋白21個, 表達下調蛋白32個, 共53個。蛋白含量相差3.0倍以上的有1個, 含量相差2.9~2.0倍2個, 含量相差1.9~1.0倍有13個。

    2.6 差異表達蛋白的生物信息功能分析

    2.6.1 青稞差異蛋白GO富集分析 對差異蛋白主要功能進行GO功能注釋, 分析結果如圖3, 差異蛋白的生物學功能覆蓋范圍分布較廣。HDZD-HDZCK對照組中, 分子功能類占比66%, 主要涉及催化活性及結合等功能; 生物過程類占比17%, 主要涉及生物合成過程、細胞過程以及代謝過程; 細胞組成類占比17%, 主要涉及細胞膜、細胞器以及大分子復合物等細胞組分。DMD-DMCK比較組中, 分子功能類占比71%, 主要涉及催化活性及結合功能; 生物過程類占比14%, 主要涉及單一生物合成過程、代謝過程以及脂質轉運過程等; 細胞組成類占比15%, 主要由細胞器、細胞膜、大分子復合物以及胞質溶膠四部分組成。HDZD-HDZCK和DMD-DMCK的差異蛋白共有184個, 對其行使的功能進行分類, 可以分為12類(表3), 分屬于碳水化合物和能量代謝類、蛋白質生物合成類、氧化還原類、氨基酸轉運與代謝類、蛋白質代謝類、應激與防御類、電子載體類等。

    A: 屬于生物過程(BP); B: 屬于細胞組分(CC); C、D、E、F、G屬于分子功能(MF); A: 防御響應; B: 胞外部分; C: 催化活性; D: ATP結合; E: 結構分子活性; F: 轉移活性; G: 運輸活性; H: 其他活性。

    A: biological process (BP), B: cell component (cc), and C–G: molecular functions (MF). A: defense response; B: extracellular part; C: catalytic activity; D: ATP binding; E: structural molecule activity; F: transfer activity; G: transport activity; H: other activities.

    表3 HDZD-HDZCKA和DMD-DMCK差異表達蛋白功能分類

    2.6.2 差異表達蛋白的KEGG分析 通過對差異蛋白進行KEGG分析, 兩組對照注釋蛋白通路見圖4。其中占比最大的3條通路分別是代謝通路、氨基酸的生物合成以及次級代謝產物的生物合成。差異表達蛋白主要涉及13條代謝通路, 包括檸檬酸循環(huán)、光合作用、半胱氨酸和蛋氨酸等通路。其中表達上調和下調蛋白共同參與的途徑有6條, 半胱氨酸和蛋氨酸代謝通路在干旱脅迫的2個品種中都表現為表達上調。在HDZD-HDZCK中, 與呼吸和光合作用相關的通路均顯著上調, 其中包括三羧酸循環(huán)(TCA)、碳代謝、丙酮酸代謝、乙醛酸和二羧酸的代謝以及光合作用中的碳固定通路; 在DMD-DMCK中, 關于蛋白質的合成通路均顯著上調, 其中包括核糖體組裝、RNA轉運、內質網中的蛋白質加工以及可提供能量的氧化磷酸化通路。

    表4 顯著差異蛋白的KEGG通路

    3 討論

    本研究以2份干旱脅迫耐受性不同的青稞品種(旱地紫青稞和大麻青稞)為研究材料, 對其葉綠素含量、可溶性蛋白含量等進行了生理指標的測定, 其次采用同位素標記相對和絕對定量(iTRAQ)技術, 分析了干旱脅迫后的葉片全蛋白組的差異表達。

    3.1 干旱處理對青稞葉綠素含量的影響

    在干旱處理7 d和10 d后, 旱地紫青稞和大麻青稞葉綠素、及葉綠素總量對比對照均呈現為持續(xù)下降的趨勢, 與之前對馬鈴薯[17]、文冠果[18]等研究結果不同, 沒有出現抗旱性強的植物葉綠素含量隨著脅迫時間的延長呈現先上升再下降的趨勢, 可能是由于不同植物內葉綠素響應干旱脅迫的機制存在差異, 但大麻青稞葉片葉綠素、葉綠素、總葉綠素含量的降低趨勢更為明顯, 證明在重度干旱脅迫時大麻青稞葉片的光合作用能力受到的抑制高于旱地紫青稞。干旱組的2個品種的葉綠素含量較對照組比葉綠素含量下降趨勢明顯, 說明重度干旱脅迫會造成葉綠素分解, 葉綠素含量降低, 并且葉綠素相比較葉綠素對干旱脅迫更加敏感, 這與對徐銀萍等對大麥干旱脅迫后葉片的葉綠素研究結果一致[19]。

    3.2 干旱處理對青稞滲透調節(jié)物質的影響

    植物受到干旱脅迫時, 會產生與抗逆活動相關的蛋白質并會水解一些結構蛋白以維持正常的生理需要, 導致細胞內可溶性蛋白含量增加[20]。在正常生理環(huán)境下, 旱地紫青稞體內的可溶性蛋白質含量低于大麻青稞。隨著干旱脅迫時間的延長, 2個品種青稞內的可溶性蛋白含量均表現出先快后慢的下降趨勢, 大麻青稞蛋白質含量的減少值和降低比例始終大于旱地紫青稞, 這與對甘薯研究結果相一致[21]。

    相對電導率在衡量植物細胞膜穩(wěn)定性中具有重要作用。在正常生理環(huán)境下, 旱地紫青稞的相對電導率低于大麻青稞的相對電導率。隨著干旱脅迫時間的延長, 兩個品種青稞的相對電導率均表現為先慢后快的增加趨勢, 這說明干旱脅迫破壞了細胞膜的穩(wěn)定性, 大量溶質滲出, 導致相對電導率增大, 這種現象在易家寧等[22]對紫蘇以及鄧輝茗等[23]對綿毛水蘇幼苗的水分脅迫的研究中都得到了印證。

    丙二醛含量主要反映細胞膜脂過氧化的程度, 含量越高細胞膜脂過氧化程度越高, 細胞膜穩(wěn)定性越差[24]。在正常生理環(huán)境下, 旱地紫青稞的丙二醛含量高于大麻青稞。隨著干旱脅迫時間的延長, 前者的丙二醛含量持續(xù)上升, 而大麻青稞的丙二醛含量的增加呈現出先慢后快的趨勢, 并且在干旱處理10 d的時候, 其丙二醛含量超過了旱地紫青稞。由此暗示, 旱地紫青稞中的細胞膜能維持相對的穩(wěn)定性, 受干旱影響較小。

    3.3 干旱處理對青稞蛋白表達的影響

    本研究共篩選到184個差異表達蛋白, 其中上調表達蛋白89個, 下調表達蛋白95個。旱地紫青稞和大麻青稞表達差異量最高的蛋白質分別是異檸檬酸裂解酶和β-1,3-葡聚糖酶, 上調倍數分別是16.64和24.03; 最顯著下調的蛋白質是DUF3700 domain-containing protein以及DUF3479 domain-containing protein, 下調倍數分別為7.69和2.25。D1-吡咯啉-5-羧酸合成酶(P5CS), 脫水蛋白及BURP domain-containing protein等蛋白質在2個品種的青稞中的表達量均上調, 同時脂氧合酶、類萌發(fā)蛋白等蛋白質也出現不同程度的差異表達。

    異檸檬酸裂解酶在旱地紫青稞中上調表達倍數最大, 為16.64倍。異檸檬酸裂解酶是乙醛酸循環(huán)體中特有的酶, 可將異檸檬酸裂解為琥珀酸和乙醛酸。后期生成的草酰乙酸可參加機體內的能量代謝過程或生成其他代謝物質, 乙醛酸循環(huán)可將油脂轉化為糖類[25]。目前有關異檸檬酸裂解酶在非生物脅迫下的作用研究還較少, 而本實驗中在旱地紫青稞中異檸檬酸酸裂解酶含量顯著上調, 其作用可能是在干旱脅迫下為機體供能以維持正常的生理代謝, 防止植物受到損害, 同時可能是利用乙醛酸循環(huán)將油脂轉化為糖類, 可增強細胞的滲透調節(jié)。

    在大麻青稞中受干旱脅迫上調表達倍數最高的蛋白為β-1,3-葡聚糖酶, 上調24.03倍, 其隸屬于糖基水解酶, 參與植物體的器官發(fā)育、果實成熟及衰老等重要生命活動[26]。之前的研究結果顯示, β-1,3-葡聚糖酶可提高植物的抗病性[27-28]、還可受到低溫脅迫[29]、水楊酸、乙烯等外源植物激素的誘導, 增強植物的抗性。徐小萍等[29]對三明野生蕉低溫處理發(fā)現該酶對低溫的響應活性不僅存在時間的臨界點, 還存在溫度臨界點, 在短時間內對低溫響應強烈, 且僅在某一臨界溫度下響應低溫脅迫, 是一種低溫保護酶[30]。在其他非生物脅迫中, β-1,3-葡聚糖酶的研究還知之甚少, 尤其是在干旱脅迫下, 發(fā)揮怎樣的作用及其作用機理, 也是需要后期著重研究的內容。

    D1-吡咯啉-5-羧酸合成酶(P5CS), 在旱地紫青稞和大麻青稞的干旱組中表現為表達量上調, 上調倍數分別是2.40和2.15, 在旱地紫青稞的表達量略高。P5CS是合成脯氨酸的關鍵酶, 而脯氨酸是植物內一種重要的滲透調節(jié)物質, 在非生物脅迫下發(fā)揮重要的作用[31]。前人已對煙草[32]、菜豆[33]、枸杞[34]、堿蓬[35]中的P5CS基因進行研究, 發(fā)現其在干旱脅迫、鹽脅迫、低溫脅迫以及ABA脅迫處理后表達上調、抗性增強。堿蓬的基因導入擬南芥并過量表達, 干旱脅迫后, 發(fā)現游離脯氨酸含量增加, 丙二醛的含量顯著降低, 擬南芥抗旱性增加, 說明P5CS是一種受干旱脅迫的正調控的蛋白[35]。

    DUFs雖然是一類功能還未研究透徹的蛋白家族, 但經研究發(fā)現DUFs在植物正常生長發(fā)育、生物脅迫防御以及響應環(huán)境逆境危害等方面發(fā)揮著至關重要的作用[36]。在不同作物中, DUFs受到不同非生物脅迫的調控。擬南芥的DUF6對鹽、干旱脅迫敏感[37]; 玉米的DUF642受酸鋁脅迫調控[38]; 小麥的DUF860響應熱脅迫[39]; 沙冬青的DUF1517具有一定的抗寒作用[40]; 水稻中的DUF則受到干旱脅迫調控[41]。大量實驗證據表明, DUFs受到非生物脅迫的正調控, 但在本次實驗中, 受到干旱脅迫的旱地紫青稞和大麻青稞的DUF3700和DUF3479蛋白均呈下調表達, 下調倍數分別是7.69倍和2.22倍, 這與現有的研究報道結果不相一致, 可能與在青稞中的其他調控網絡有關, 還需要后續(xù)在其他組學中進行比對分析其原因。

    DHNs是LEA蛋白第二家族成員, 因其C端有一段保守的K片段, 在植物失水時可以形成a-螺旋, 以防止細胞膜受損以及細胞內蛋白質變形[42]。DHN5在‘旱地紫’青稞和DHN8在‘大麻’青稞的干旱組中均表現為表達量上調, 分別上調2.89倍和15.65倍。DHNs雖均屬于同一個家族, 但其功能及作用機理不盡相同[43], 有主要受干旱和ABA誘導表達的, 有被低溫和寒冷誘導型, 也有與植物抗旱和抗旱相關型, 在植物中多表現為組成型行使其抗逆功能。在小麥中, 表現為對干旱和鹽脅迫有響應[44-45], 對冷脅迫下的大麥也表現出有增強其抗逆性的作用[46], 在擬南芥中過表達的DHN5還可參與到JA信號轉導中, 下調3個JA信號轉導負調控因子JAZ家族成員的蛋白表達, 從而增強其對植物病原菌的防御作用[47]。在早稻中的DHN8可受干旱和高鹽誘導表達, 但不受低溫誘導表達。DHN5及DHN8在這兩種青稞種質中抵御脅迫時的作用機理可能存在差異, 且其表達量上調倍數相差懸殊, 但具體機理不明, 這是后續(xù)需要著重研究的部分。

    脂氧合酶(LOX)催化以亞油酸、亞麻酸等多元不飽和脂肪酸為底物的膜脂過氧化反應, 生成氫過氧化物, 破壞膜的完整性[48]。馬明杰等[49]通過對冷敏性蔬菜青椒低溫處理, 發(fā)現LOX活性增加, 說明LOX對低溫敏感, 活性增強后破壞膜脂從而破壞膜的完整性, 導致青椒的冷敏性。王成慧[50]通過轉基因手段將甜瓜轉入擬南芥, 并構建沉默植株以及過表達植株, 進行干旱處理后發(fā)現過表達植株LOX活性升高, 過氧化氫含量降低, 減少了氧化傷害, 增強了植株的抗旱性; 而沉默植株LOX活性降低, 過氧化氫積累, 對干旱脅迫敏感。孫婷婷等[51]將甘蔗的基因導入本氏煙草中, 在氯化鈉和PEG處理后, 該基因表達量顯著上升, 增強了植物的抗旱性。本次實驗中, LOX的活性在旱地紫青稞中表達下調, 在大麻青稞中相反, 說明在大麻青稞中LOX的活性更容易受到干旱脅迫的調控, 推測與大麻青稞是干旱敏感型有關。通過比較二者干旱第7天與第10天的電導率結果發(fā)現, 這與干旱脅迫下大麻青稞的電導率高于同等處理時間旱地紫青稞變化的生理狀態(tài)一致。

    類萌發(fā)素蛋白(GLPs)是一種植物體內普遍存在的可溶性糖蛋白, 既可以作為酶和結構蛋白催化調控植物體內生理反應, 也可作為受體參加信號轉導[52]。GLPs在兩個品種青稞中的表達結果相反, 在旱地紫青稞中下調表達, 下調1.92倍; 在大麻青稞中表達上調, 上調倍數是2.19。之前的研究發(fā)現, GLPs既可以通過介導細胞壁增厚來阻止病原菌入侵響應生物脅迫, 也可以對鹽脅迫、干旱脅迫、低溫脅迫[53]以及重金屬[54]等非生物脅迫做出具有組織特異性的應答。Hurkman等[55]在研究鹽脅迫大麥幼苗時發(fā)現, 在大麥幼苗受到鹽脅迫時, 僅在根部特異性上調GLPs表達。本次實驗GLPs的差異表達說明旱地紫青稞中類萌發(fā)蛋白對干旱脅迫的敏感程度低于大麻青稞, 這也可能與旱地紫青稞的耐旱性強于大麻青稞有關。

    4 結論

    干旱脅迫對青稞的生長和代謝具有一定的破壞, 從而導致青稞減產減量。干旱脅迫致使青稞葉片細胞膜穩(wěn)定性下降, 膜脂過氧化程度升高, 滲透調節(jié)物質大量積累, 葉綠素含量降低, 影響了青稞的正常生理代謝過程。異檸檬酸裂解酶、β-1,3-葡聚糖酶、D1-吡咯啉-5-羧酸合成酶、類萌發(fā)素蛋白等代謝物也出現顯著差異表達的現象, 通過表達的上調或下調, 增加青稞的抗逆性, 使其在干旱環(huán)境中足以存活且保證一定的產量。

    [1] 唐益苗, 趙昌平, 高世慶, 田立平, 單福華, 吳敬新. 植物抗旱相關基因研究進展. 麥類作物學報, 2009, 29: 166–173. Tang Y M, Zhao C P, Gao S Q, Tian L P, Shan F H, Wu J X. Advance in genes related to plant drought tolerance.,2009, 29: 166–173 (in Chinese with English abstract).

    [2] 蔡昆爭, 吳學祝, 駱世明. 不同生育期水分脅迫對水稻根葉滲透調節(jié)物質變化的影響. 植物生態(tài)學報, 2008, 32: 491–500. Cai K Z, Wu X Z, Luo S M. Effects of water stress on osmolytes at different growth stages in rice leaves and roots., 2008, 32: 491–500 (in Chinese with English abstract).

    [3] 王一鳴, 汪湖, 龍勝舉, 趙英鵬, 陳延, 賀忠群. 干旱脅迫對蒲公英滲透調節(jié)物質?酶保護系統(tǒng)及質膜水孔蛋白基因表達的影響. 華北農學報, 2017, 32(4): 85–90. Wang Y M, Wang H, Long S J, Zhao Y P, Chen Y, He Z Q. Effects of drought stress on osmoregulation substances, enzyme protection system and expression of plasma membrane hole protein gene PIP2-3 inhand-mazz., 2017, 32(4): 85–90 (in Chinese with English abstract).

    [4] 李春香, 王瑋, 李德全. 長期水分脅迫對小麥生育中后期根葉滲透調節(jié)能力?滲透調節(jié)物質的影響. 西北植物學報, 21: 924–930. Li C X, Wang W, Li D Q. Effects of long-term water stress on osmotic adjustment and osmolytes in wheat roots and leaves., 2001, 21: 924–930 (in Chinese with English abstract).

    [5] 姚曉華, 吳昆侖. PEG預處理對青稞種子萌發(fā)和幼苗生理特性的影響. 西北植物學報, 2012, 32: 1403–1411. Yao X H, Wu K L. Effect of PEG pretreatment on germination and growth physiology of hulless barley., 2012, 32: 1403–1411 (in Chinese with English abstract).

    [6] 王玉林, 徐齊君, 原紅軍, 曾興權, 尼瑪扎西. PEG模擬干旱脅迫處理對青稞幼苗生長和生理特性的影響. 大麥與谷類科學, 2018, 35(1): 6–12. Wang Y L, Xu Q J, Yuan H J, Zeng X Q, Nima T S. Impact of drought stress induced by polyethylene glycol on the growth and physiological characteristics of hulless barley seedlings., 2018, 35(1): 6–12 (in Chinese with English abstract).

    [7] 徐齊君, 原紅軍, 曾興權, 王玉林, 扎桑, 于明寨, 頓珠加布, 尼瑪扎西. 西藏青稞抗旱研究進展. 西藏農業(yè)科技, 2019, 41(1): 64–67. Xu Q J, Yuan H J, Zeng X Q, Wang Y L, Zha S, Yu M Z, Dunzhu J B, Nima Z X. Research progress on drought resistance of Tibetan highland barley., 2019, 41(1): 64–67 (in Chinese with English abstract).

    [8] 白羿雄, 姚曉華, 姚有華, 吳昆侖. 青稞抗倒伏性狀的基因型差異. 中國農業(yè)科學, 2019, 52: 228–238. Bai Y X, Yao X H, Yao Y H, Wu K L. Difference of traits relating to lodging resistance in hulless barley genotypes., 2019, 52: 228–238 (in Chinese with English abstract).

    [9] 蔡靜. 溫度對青稞非生物脅迫抗性基因節(jié)律性表達的影響研究. 西北大學碩士學位論文, 陜西西安, 2018. Cai J. Study on the Effect of Environmental Temperature on Rhythmic Expression of Abiotic Stress Resistance Genes in Tibetan Hulless Barley. MS Thesis of Northwest University, Xi’an, Shaanxi, China, 2018 (in Chinese with English abstract).

    [10] 姚曉華, 吳昆侖. 青稞脂質轉運蛋白基因的克隆及其對非生物脅迫的響應. 作物學報, 2016, 42: 399–406. Yao X H, Wu K L. Isolation ofgene encoding LTP protein in hulless barley and its response to abiotic stresses., 2016, 42: 399–406 (in Chinese with English abstract).

    [11] 賴海濤, 呂禹澤, 蘇國成, 林麗娟. 螺旋藻葉綠素的提取工藝研究. 廣東化工, 2020, 47(1): 11–12. Lai H T, Lyu Y Z, Su G C, Lin L J. Study on the extraction technology of chlorophyll from spirulina., 2020, 47 (1): 11–12 (in Chinese with English abstract).

    [12] Bradford M M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding., 1976, 72: 248–254.

    [13] 趙世杰, 許長成, 鄒琦, 孟慶偉. 植物組織中丙二醛測定方法的改進. 植物生理學通訊, 1994, 30: 207–210. Zhao S J, Xu C C, Zou Q, Meng Q W. Improvement of measurement of malondialdehvde in plant tissues., 1994, 30: 207–210 (in Chinese with English abstract).

    [14] 李紅兵, 康振生. 適于小麥葉片蛋白質組分析的樣品提取方法研究. 西北植物學報, 2011, 31: 632–638. Li H B, Kang Z S. Sample preparation methods suitable for wheat leaf proteome analysis., 2011, 31: 632–638 (in Chinese with English abstract).

    [15] Evans C, Noirel J, Ow S Y, Salim M, Pereira M A, Couto N, Pandhal J, Smith D, Pham T K, Karunakaran E, Zou X, Biggs C A, Phillip C. An insight into iTRAQ: where do we stand now?, 2012, 404: 1011–1027.

    [16] Kocher T, Pichler P, Schutzbier M, Stingl C, Kaul A, Teucher N, Hasenfuss G, Penninger J M, Mechtler K. High precision quantitative proteomics using iTRAQ on an LTQ orbitrap: a new mass spectrometric method combining the benefits of all., 2009, 8: 4743–4752.

    [17] 楊喜珍, 楊利, 覃亞, 劉磊, 楊歡, 謝婉. PEG-8000模擬干旱脅迫對馬鈴薯組培苗葉綠素和類胡蘿卜素含量的影響. 中國馬鈴薯, 2019, 33(4): 193–202. Yang X Z, Yang L, Qin Y, Liu L, Yang H, Xie W. Effects of PEG-8000 stress on contents of chlorophyll and carotenoid of potato plantlets.,2019, 33(4): 193–202 (in Chinese with English abstract).

    [18] 權伍榮. 干旱脅迫下不同文冠果品種抗氧化酶活性及葉綠素含量變化. 延邊大學農學學報, 2020, 42(1): 15–21. Quan W R. Changes of antioxidant enzyme activities and chlorophyll content in differentvarieties under drought stress., 2020, 42(1): 15–21 (in Chinese with English abstract).

    [19] 徐銀萍, 潘永東, 劉強德, 任誠, 姚元虎, 賈延春, 陳文慶, 火克倉, 包奇軍, 趙鋒, 張華瑜. 生長后期干旱復水對飼草大麥產量?品質及葉綠素含量的影響. 中國土壤與肥料, 2020, (2): 192–197. Xun Y P, Pan Y D, Liu Q D, Ren C, Yao Y H, Jia Y C, Chen W Q, Huo K C, Bao Q J, Zhao F, Zhang H Y. Effects of rewatering after drought stress on yield quality and chlorophyll of forage barley in late growth period., 2020, (2): 192–197 (in Chinese with English abstract).

    [20] 尹啟琳, 郭丁預, 姜倩倩, 張立培, 陳磊, 趙婧, 宋建成, 趙吉強. 干旱脅迫對不同小麥品種苗期抗旱生理指標的影響. 煙臺大學學報(自然科學與工程版), 2020, (3): 289–297. Yin Q L, Guo D Y, Jiang Q Q, Zhang L P, Chen L, Zhao J, Song J C, Zhao Q. Effects of drought stress on drought tolerant physiological indexes of different wheat varieties at seedling stage.(Nat Sci Eng Edn), 2020, (3): 289–297 (in Chinese with English abstract).

    [21] 張海燕, 汪寶卿, 馮向陽, 李廣亮, 解備濤, 董順旭, 段文學, 張立明. 不同時期干旱脅迫對甘薯生長和滲透調節(jié)能力的影響. 作物學報, 2020, 46: 1760–1770. Zhang H Y, Wang B Q, Feng X Y, Li G L, Xie B T, Dong S X, Duan W X, Zhang L M. Effects of drought treatments at different growth stages on growth and the activity of osmotic adjustment in sweet potato [(L.) Lam]., 2020, 46: 1760–1770 (in Chinese with English abstract).

    [22] 易家寧, 王康才, 張琪綺, 董雨青, 毛曉敏, 鄧艷婷. 干旱脅迫對紫蘇生長及品質的影響. 核農學報, 2020, 34: 1320–1326. Yi J N, Wang K C, Zhang Q Q, Dong Y Q, Mao X M, Deng Y T. Effects of drought stress on growth and quality of., 2020, 34: 1320–1326 (in Chinese with English abstract).

    [23] 鄧輝茗, 龍聰穎, 蔡仕珍, 宋宇, 鄢如霞, 車亦然, 王長見, 肖瑤. 不同水分脅迫對綿毛水蘇幼苗形態(tài)和生理特性的影響. 西北植物學報, 2018, 38: 1099–1108. Deng H M, Long C Y, CAI S Z, Song Y, Yan R X, Che Y R, Wang C J, Xiao Y. Morphology and physiological characteristics ofseeding under water stress.,2018, 38: 1099–1108 (in Chinese with English abstract).

    [24] 韓志順, 鄭敏娜, 梁秀芝, 康佳惠, 陳燕妮. 干旱脅迫對不同紫花苜蓿品種形態(tài)特征和生理特性的影響. 中國草地學報, 2020, 42(3): 37–43. Han Z S, Zheng M N, Liang X Z, Kang J H, Chen Y N. Effects of drought stress on morpho logical and physioloical characteristics of different alfalfa cultivars., 2020, 42(3): 37–43 (in Chinese with English abstract).

    [25] 張潔瓊. 異檸檬酸裂解酶基因克隆與表達特征研究. 長春工業(yè)大學碩士學位論文, 吉林長春, 2010. Zhang J Q. Cloning and Expression of Isocitrate Lyase Gene. MS Thesis of Changchun University of Technology, Changchun, Jilin, China, 2010 (in Chinese with English abstract).

    [26] 歐陽波, 李漢霞, 葉志彪. 植物β-1,3-葡聚糖酶及其基因. 中國生物工程雜志, 2002, (6): 18–23. Ou-Yang B, Li H X, Ye Z B. Plant β-1,3-glucanase and its related genes., 2002, (6): 18–23 (in Chinese with English abstract).

    [27] Lusso M, Kuc J. The effect of sense and antisense expression of PRN gene for β-1,3-glucanaseon disease resistance of tobacco to fungi and virus., 1996, 49: 267–283.

    [28] Nakamura Y, Sawada H, Kobayashi S, Yoshikawa M. Expression of soybean β-1,3-endoglucanase cDNA and effect on disease tolerance in kiwi fruit plants., 1999, 18: 527–532.

    [29] 徐小萍, 謝燕萍, 陳芳蘭, 陳曉慧, 陳裕坤, 張梓浩, 程春振, 林玉玲, 賴鐘雄. 三明野生蕉β-1,3-葡聚糖酶基因克隆及其在低溫處理下的表達分析. 熱帶作物學報, 2020, 41(2): 292–299. Xu X P, Xie Y P, Chen F L, Chen X H, Chen Y K, Zhang Z H, Chen C Z, Lin Y L, Lai Z X. Cloning and characterization of, a β-1,3-glucanase gene of wild banana germplasm from Sanming, Fujian, China, under cold treatment., 2020, 41: 292–299 (in Chinese with English abstract).

    [30] 陳芳蘭. 野生蕉β-1,3葡聚糖酶基因克隆及抗寒相關功能分析. 福建農林大學碩士學位論文, 福建福州, 2016. Chen F L. Cloning and Cold Resistance Analysis of β-1,3 Glucanase Gene () from Wild Banana. MS Thesis of Fujian Agriculture and Forestry University, Fujian, Fuzhou, China, 2016 (in Chinese with English abstract).

    [31] 李健, 杜成忠, 王露蓉, 邢永秀, 李楊瑞, 楊麗濤. 甘蔗D1-吡咯啉-5-羧酸合成酶基因轉化煙草的抗旱性分析. 華中農業(yè)大學學報, 2018, 37(2): 34–42. Li J, Do T R, Wang L R, Xing Y X, Li Y R, Yang L T. Drought resistance of transgenic tobacco with sugarcaneD1-pyrroline-5- carboxylate synthase (SoP5CS) gene., 2018, 37(2): 34–42 (in Chinese with English abstract).

    [32] 焦蓉, 劉貫山, 劉好寶, 王樹林, 侯娜, 王全貞, 靳義榮, 白巖, 馮廣林, 劉朝科, 馮祥國, 胡曉明. 普通煙草抗?jié)B透脅迫基因的克隆與表達分析. 中國煙草學報, 2012, 18(2): 49–57. Jiao R, Liu G S, Liu H B, Wang S L, Hou N, Wang Q Z, Jin Y R, Bai Y, Feng G L, Liu C K, Feng X G, Hu X M. Cloning and expression analysis of osmotic stress tolerance genein., 2012, 18(2): 49–57 (in Chinese with English abstract).

    [33] 陳吉寶. 普通菜豆基因的克隆?功能驗證及單核苷酸多態(tài)性. 中國農業(yè)科學院博士學位論文, 北京, 2008. Chen J B. Cloning, Function Analysis and Single Nucleotide Polymorphism of Common Bean (L.)Gene. PhD Dissertation of Chinese Academy of Agricultural Science, Beijing, China, 2008 (in Chinese with English abstract).

    [34] 馮遠航, 王罡, 季靜, 關春峰, 金超. 枸杞基因的克隆及表達分析. 中國生物工程雜志, 2013, 33(1): 33–40. Feng Y H, Wang G, Ji J, Guan C F, Jin C. Cloning and expression analysis ofgene fromMiller., 2013, 31(4): 572–578 (in Chinese with English abstract).

    [35] 郭丹丹, 楊清華, 朱丹華, 金杭霞. 堿蓬基因過表達提高擬南芥耐旱性. 浙江農業(yè)學報, 2019, 31: 572–578. Guo D D, Yang Q H, Zhu D H, Jin H X. Overexpression ofgene fromenhances drought tolerance of transgenic., 2019, 31: 572–578 (in Chinese with English abstract).

    [36] 楊杞, 牛肖翠, 王瑞剛, 李國婧, 張國盛. 蒺藜苜蓿DUF221基因家族全基因組鑒定及鹽響應相關基因篩選. 分子植物育種, 2019, 17: 5255–5262. Yang Q, Niu X C, Wang R G, Li G J, Zhang G S. Genome-wide characterization of DUF221 gene family inand screening for salt response genes., 2019, 17: 5255–5262 (in Chinese with English abstract).

    [37] 孫書琦, 張銳, 郭三堆. 擬南芥Nodulin Mt N21家族基因功能的初步研究. 核農學報, 2009, 23: 429–434. Sun S Q, Zhang R, Guo S D. Primary analysis of the gunction of Nodulin Mt N21 genein., 2009, 23: 429–434 (in Chinese with English abstract).

    [38] 劉天宇. 玉米DUF642基因家族的鑒定和分析. 四川農業(yè)大學碩士學位論文, 四川成都, 2018. Liu T Y. Identification and Characterization of the DUF642 Family in Maize. MS Thesis of Sichuan Agricultural University, Chengdu, Sichuan, China, 2018 (in Chinese with English abstract).

    [39] 秦丹丹, 謝頌朝, 劉剛, 倪中福, 姚穎垠, 孫其信, 彭惠茹. 小麥中編碼未知蛋白的熱脅迫響應基因的克隆和功能分析. 植物學報, 2013, 48: 34–41. Qin D D, Xie S C, Liu G, Ni Z F, Yao Y Y, Sun Q X, Peng H R. Isolation and functional characterization of heat-stress responsive genefrom wheat., 2013, 48: 34–41 (in Chinese with English abstract).

    [40] Gu L, Cheng H. Isolation, molecular cloning and characterization of a cold-responsive gene, Am DUF1517, from., 2014, 117: 201–211.

    [41] 周炎, 樊帆, 雷東陽, 盧學丹. 水稻逆境響應蛋白OsSGL的生物信息學分析及原核表達條件優(yōu)化. 分子植物育種, 2021, 19: 1531–1540. Zhou Y, Fan F, Lei D Y, Lu X D. Bioinformatics analysis and opitmization of protein expression in prokaryotic system of rice stress responsive protein OsSGL., 2021, 19: 1531–1540.

    [42] 齊玉紅. 干旱脅迫下小麥脫水素WZY1-2?WZY2互作蛋白的篩選. 西北農林科技大學碩士學位論文, 陜西楊凌, 2018. Qi Y H. Screening of Wheat Dehydrin WZY1-2 and WZY2 Interacting Proteins Under Drought Stress. MS Thesis of Northwest Agricultural University, Yangling, Shaanxi, China, 2018 (in Chinese with English abstract).

    [43] Rorat T. Plant dehydrins-tissue location, structure and function., 2006, 11: 536–556.

    [44] Brini F, Saibi W, Amara I, GargouriI A, Masmoudi K, Hanin M. Wheat dehydrin DHN5 exerts a heat-protective effect on β-glucosidase and glucose oxidase activities., 2010,74: 1050–1054.

    [45] Brini F, Hanin M, Lumbreras V, Amara I, Khoudi H, Hassairi A, Pagés M, Masmoudi K. Overexpression of wheat dehydrin DHN-5 enhances tolerance to salt and osmotic stress in., 2007, 26: 2017–2026.

    [46] Bravoa L A, Gallardoa J, Navarretea A, Olavea N, Martínezb J, Alberdic M, Closed T J, Corcuera L J. Cryoprotective activity of a cold-induced dehydrin purified from barley., 2003, 118: 262–269.

    [47] Brini F, Yamamoto A, Jlaiel L, Takeda S, Hobo T, Dinh H, Hattori T, Masmoudi K, Hanin M. Pleiotropic effects of the wheat dehydrin DHN-5 on stress responses in., 2011, 52: 676–688.

    [48] 吳桂玲, 馮定坤. 植物脂氧合酶的研究進展. 廣州化工, 2019, 47(17): 37–39. Wu G L, Feng D K. Research progress on plant lipoxygenase., 2019, 47(17): 37–39 (in Chinese with English abstract).

    [49] 馬明杰, 程順昌, 紀淑娟, 胡美斯, 陰曉晨, 賈文韜, 魏寶東. 低溫脅迫對青椒膜脂代謝的影響. 包裝工程, 2020, 41(3): 21–27.Ma M J, Cheng S C, Ji S J, Hu M S, Yin X C, Jia W T, Wei B D. Effects of low temperature stress on membrane lipid metabolism of., 2020, 41(3): 21–27 (in Chinese with English abstract).

    [50] 王成慧. 薄皮甜瓜基因的克隆及其在干旱和鹽脅迫中的功能分析. 沈陽農業(yè)大學博士學位論文, 遼寧沈陽, 2019. Wang C H. Cloning and Functional Analysis ofin Drought and Salt Stresses from Oriental Melon (var.MakuwaMakino). PhD Dissertation of Shenyang Agricultural University, Shenyang, Liaoning, China, 2019 (in Chinese with English abstract).

    [51] 孫婷婷, 王文舉, 婁文月, 劉峰, 張旭, 王玲, 陳玉鳳, 闕友雄, 許莉萍, 李大妹, 蘇亞. 甘蔗脂氧合酶基因的克隆與表達分析. 作物學報, 2019, 45: 1002–1016. Sun T T, Wang W J, Lou W Y, Liu F, Zhang X, Wang L, Chen Y F, Que Y X, Xu L P, Li D M, Su Y C. Cloning and expression analysis of sugarcane lipoxygenase gene., 2019, 45: 1002–1016 (in Chinese with English abstract).

    [52] 李瑩, 柳參奎. 植物類萌發(fā)素蛋白研究進展. 中國農學通報, 2014, 30(30): 246–254. Li Y, Liu C K. Research progresson plant germin-like proteins., 2014, 30(30): 246–254 (in Chinese with English abstract).

    [53] Li H Y, Jiang J, Wang S, Liu F F. Expression analysis of ThGLP, a new germin-like protein gene, in., 2010, 21: 323–330.

    [54] Berna A, Bernier F. Regulation by biotic and abiotic stress of a wheat germin gene encoding oxalate oxidase, a H2O2-producing enzyme., 1999, 39: 539–549.

    [55] Hurkman W J, Tao H P, Tanaka C K. Germin-like polypeptides increase in barley roots during salt stress., 1991, 97: 366–374.

    Differentially expressed protein analysis of different drought tolerance hulless barley leaves

    LI Jie**, FU Hui**, YAO Xiao-Hua, and WU Kun-Lun*

    Academy of Agricultural and Forestry Sciences, Qinghai University / Qinghai Academy of Agricultural and Forestry Sciences / Qinghai Key Laboratory of Hulless Barley Genetics and Breeding / Qinghai Subcenter of National Hulless Barley Improvement, Xining 810016, Qinghai, China

    In order to reveal the differences in response to drought stress among different hulless barley varieties from the protein level and to analyze the protein molecular mechanism of drought tolerance, Handizi Barley (HDZ) resistant to drought stress and Dama Barley (DM) sensitive to drought stress were used as research materials in this study. Drought treatment were determined by potted-planting method with limited water supply, four physiological indexes of hulless barley leaves with different drought gradients, including chlorophyll, soluble protein, malondialdehyde content, and relative conductivity were investigated. iTRAQ technology was used to conduct differential protein analysis on the whole protein group of barley leaves under deep drought stress. The results showed that with the extension of the drought treatment, the chlorophyll and soluble protein content of two hulless barleys under drought stress gradually decreased, the electrical conductivity and malondialdehyde content gradually increased, and the decrease in chlorophyll and soluble protein content, the increase in electrical conductivity and the content of malondialdehyde in Dama were greater than that of Handizi; 4163 proteins (polypeptides) were quantified, among them, compared with normal culture in the Handizi comparison group, 68 up-regulated proteins and 63 down-regulated proteins were screened by iTRAQ; in the comparative group of Dama, 21 up-regulated proteins and 32 down-regulated proteins were screened. KEGG pathway showed that the top three enrichment pathways were metabolic, amino acid biosynthesis, and secondary metabolite biosynthesis. The first one mainly related to citric acid cycle, carbon cycle, and other metabolic pathways. The synthesis and degradation of amino acids were mainly involved arginine and alanine. The synthesis of secondary metabolites were about arachidonic acid and linolenic acid. This study screened the proteins related to the metabolic pathways and other related functions in response to drought stress on proteome level in hulless barley, providing a theoretical basis for revealing the molecular regulation mechanism in response to drought stress.

    hulless barley; drought stress; leaf; differentially expressed protein; iTRAQ

    10.3724/SP.J.1006.2021.01062

    本研究由國家自然科學基金項目(32060480, 31660388, 31960427), 青海省農林科學院創(chuàng)新基金(2018-NKY-12), 國家現代農業(yè)產業(yè)技術體系建設專項(CARS-05)和青海省科技廳應用基礎研究項目(2019-ZJ-7075)資助。

    This study was supported by the National Natural Science Foundation of China (32060480, 31660388, 31960427), the Qinghai Provincial Academy of Agriculture and Forestry Innovation Fund (2018-NKY-12), the China Agriculture Research System (CARS-05), and the Qinghai Science and Technology Support Project (2019-ZJ-7075).

    吳昆侖, E-mail: wklqaaf@163.com

    **同等貢獻(Contributed equally to this work)

    E-mail: lij_28257@163.com

    2020-08-05;

    2020-12-01;

    2021-01-04.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20210104.1355.010.html

    猜你喜歡
    旱地大麻青稞
    半碗青稞面
    旱地麥田夏閑期復種綠肥 保持土壤的可持續(xù)生產力
    工業(yè)大麻中大麻二酚的研究進展
    中成藥(2021年5期)2021-07-21 08:38:48
    大麻二酚在醫(yī)學上的應用前景
    旱地冰球運動開展價值的研究
    冰雪運動(2020年6期)2020-07-21 03:28:16
    這一次我又寫到青稞
    翠苑(2019年4期)2019-11-11 04:27:26
    旱地冰球運動推廣及發(fā)展策略
    冰雪運動(2019年5期)2019-08-24 08:04:48
    大麻是個啥?
    華聲(2016年20期)2016-11-19 12:38:03
    “藏青2000”青稞新品種示范推廣的成效與做法
    西藏科技(2016年8期)2016-09-26 09:00:24
    長期施肥下紅壤旱地土壤CO2排放及碳平衡特征
    乱系列少妇在线播放| 久久精品影院6| 色综合亚洲欧美另类图片| 国产私拍福利视频在线观看| 国产精品,欧美在线| 亚洲欧美清纯卡通| 91久久精品国产一区二区成人| 卡戴珊不雅视频在线播放| 免费观看的影片在线观看| 99久国产av精品国产电影| 日韩一区二区三区影片| 国产91av在线免费观看| 丝袜喷水一区| 国产伦一二天堂av在线观看| 成年版毛片免费区| 天天一区二区日本电影三级| 欧美潮喷喷水| 国产精品99久久久久久久久| 网址你懂的国产日韩在线| 永久网站在线| 亚洲成人精品中文字幕电影| www.色视频.com| 变态另类丝袜制服| 91午夜精品亚洲一区二区三区| 国产激情偷乱视频一区二区| 黄色一级大片看看| 久久久久久九九精品二区国产| 久久久久久久久久久免费av| 成人午夜精彩视频在线观看| 亚洲,欧美,日韩| 人人妻人人澡人人爽人人夜夜 | 国产一区亚洲一区在线观看| 国产精品日韩av在线免费观看| 观看美女的网站| 日韩av在线大香蕉| 99热这里只有是精品在线观看| 亚洲人成网站高清观看| 精品酒店卫生间| 国产精品一区二区在线观看99 | 日本wwww免费看| 亚洲电影在线观看av| 身体一侧抽搐| 欧美日韩一区二区视频在线观看视频在线 | av专区在线播放| 国产欧美日韩精品一区二区| 亚洲人成网站在线观看播放| 国产白丝娇喘喷水9色精品| 天天躁日日操中文字幕| 精品久久国产蜜桃| 乱人视频在线观看| 国产精品久久视频播放| 亚洲国产最新在线播放| 一本久久精品| 人妻少妇偷人精品九色| 男人狂女人下面高潮的视频| 国产色爽女视频免费观看| 欧美日本视频| 插逼视频在线观看| 成人三级黄色视频| 日本猛色少妇xxxxx猛交久久| 嫩草影院精品99| 美女cb高潮喷水在线观看| 国产在视频线在精品| 一级黄色大片毛片| 国产伦精品一区二区三区视频9| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 最近手机中文字幕大全| 天堂网av新在线| 国产免费福利视频在线观看| 国产视频首页在线观看| 美女高潮的动态| 在线播放无遮挡| 人人妻人人澡欧美一区二区| 欧美zozozo另类| 中文字幕亚洲精品专区| 日本-黄色视频高清免费观看| 美女黄网站色视频| 国产麻豆成人av免费视频| 亚洲av中文字字幕乱码综合| 美女大奶头视频| 狠狠狠狠99中文字幕| 久热久热在线精品观看| 成人高潮视频无遮挡免费网站| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影| 又爽又黄无遮挡网站| 久久鲁丝午夜福利片| 3wmmmm亚洲av在线观看| 男人的好看免费观看在线视频| 老女人水多毛片| 高清午夜精品一区二区三区| 久久久久久久午夜电影| 国产精品蜜桃在线观看| 国产av码专区亚洲av| www.色视频.com| 丰满少妇做爰视频| 国产精品电影一区二区三区| 免费看a级黄色片| 午夜亚洲福利在线播放| 乱码一卡2卡4卡精品| 国产极品精品免费视频能看的| 亚洲18禁久久av| 亚洲av免费高清在线观看| 青春草视频在线免费观看| 22中文网久久字幕| 成人午夜高清在线视频| 国产色婷婷99| 少妇人妻一区二区三区视频| 禁无遮挡网站| 日韩在线高清观看一区二区三区| 久久久久精品久久久久真实原创| 午夜精品在线福利| 欧美高清成人免费视频www| 国产极品天堂在线| 男女视频在线观看网站免费| 国产亚洲av片在线观看秒播厂 | 三级经典国产精品| 免费看光身美女| 男女啪啪激烈高潮av片| 午夜亚洲福利在线播放| 国产一区有黄有色的免费视频 | 国产免费一级a男人的天堂| 亚洲内射少妇av| 免费在线观看成人毛片| 亚洲欧美日韩东京热| 热99在线观看视频| 午夜福利视频1000在线观看| 婷婷色av中文字幕| 国产激情偷乱视频一区二区| 亚洲最大成人中文| 国内精品美女久久久久久| 成人毛片60女人毛片免费| 97人妻精品一区二区三区麻豆| 99在线人妻在线中文字幕| 日日摸夜夜添夜夜爱| 三级毛片av免费| 国产在视频线在精品| 一个人看视频在线观看www免费| 久久热精品热| 国产成年人精品一区二区| 亚洲五月天丁香| 99热6这里只有精品| 天堂网av新在线| 全区人妻精品视频| 国产伦精品一区二区三区视频9| 一级二级三级毛片免费看| 亚洲国产成人一精品久久久| 三级国产精品欧美在线观看| 免费大片18禁| 国产精品麻豆人妻色哟哟久久 | 国产一区二区亚洲精品在线观看| 久久久久久久国产电影| 在线观看av片永久免费下载| 免费搜索国产男女视频| 亚洲av中文字字幕乱码综合| 久久久成人免费电影| 国产精品久久久久久精品电影小说 | 蜜臀久久99精品久久宅男| 日日干狠狠操夜夜爽| 国产 一区 欧美 日韩| 国内精品美女久久久久久| 国产一区二区三区av在线| 在线播放无遮挡| 网址你懂的国产日韩在线| 亚洲久久久久久中文字幕| 国产成人精品一,二区| av在线亚洲专区| 亚洲精品乱码久久久久久按摩| 亚洲av福利一区| 午夜福利在线观看吧| 一卡2卡三卡四卡精品乱码亚洲| 高清视频免费观看一区二区 | 国产免费福利视频在线观看| 99热网站在线观看| 日韩av在线大香蕉| 干丝袜人妻中文字幕| 国产亚洲精品av在线| 亚洲精品日韩在线中文字幕| a级毛片免费高清观看在线播放| 国产精品嫩草影院av在线观看| 国产综合懂色| av天堂中文字幕网| 午夜福利在线观看免费完整高清在| 国产探花在线观看一区二区| 国产爱豆传媒在线观看| 精品人妻偷拍中文字幕| 国产一区二区在线av高清观看| 午夜福利视频1000在线观看| 色尼玛亚洲综合影院| 不卡视频在线观看欧美| 插逼视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 天堂av国产一区二区熟女人妻| 亚洲乱码一区二区免费版| 天天躁日日操中文字幕| 国产亚洲精品av在线| 免费av毛片视频| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 国产成人福利小说| 天堂影院成人在线观看| 国产 一区精品| 美女cb高潮喷水在线观看| 91精品一卡2卡3卡4卡| 日本wwww免费看| 午夜激情欧美在线| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 看片在线看免费视频| 免费播放大片免费观看视频在线观看 | videos熟女内射| 亚洲精品亚洲一区二区| 免费观看的影片在线观看| 伦精品一区二区三区| 黑人高潮一二区| 欧美丝袜亚洲另类| 在线a可以看的网站| a级毛片免费高清观看在线播放| 亚洲不卡免费看| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 91狼人影院| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 国内少妇人妻偷人精品xxx网站| 男女视频在线观看网站免费| 国产精品不卡视频一区二区| 少妇的逼好多水| 亚洲不卡免费看| 国产极品天堂在线| 国产真实乱freesex| 国产精品.久久久| 国产一区二区亚洲精品在线观看| 国产日韩欧美在线精品| 国产精品av视频在线免费观看| 精品无人区乱码1区二区| 欧美日韩综合久久久久久| 中文字幕熟女人妻在线| 中文字幕久久专区| 日韩欧美三级三区| 在线播放国产精品三级| 国产精品无大码| 亚洲精品影视一区二区三区av| 少妇熟女欧美另类| 成年女人永久免费观看视频| 午夜免费激情av| 乱码一卡2卡4卡精品| 人妻制服诱惑在线中文字幕| 黄色配什么色好看| 国产成人免费观看mmmm| 性插视频无遮挡在线免费观看| 如何舔出高潮| 日韩av不卡免费在线播放| 欧美成人一区二区免费高清观看| 99久久九九国产精品国产免费| 亚洲最大成人av| 精品人妻一区二区三区麻豆| 男人舔女人下体高潮全视频| 三级国产精品片| 日本黄大片高清| 亚洲欧美日韩卡通动漫| 国产精品人妻久久久久久| 亚洲自偷自拍三级| 精品人妻一区二区三区麻豆| 国产伦精品一区二区三区视频9| 美女高潮的动态| 色综合色国产| 日日啪夜夜撸| 亚洲成人av在线免费| 男人舔女人下体高潮全视频| 国产精品嫩草影院av在线观看| 小蜜桃在线观看免费完整版高清| 亚洲av中文av极速乱| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 大话2 男鬼变身卡| 欧美激情在线99| 啦啦啦韩国在线观看视频| 亚洲国产精品sss在线观看| 久久久久九九精品影院| 可以在线观看毛片的网站| 91aial.com中文字幕在线观看| 男插女下体视频免费在线播放| 18+在线观看网站| 久久久久久伊人网av| 日韩亚洲欧美综合| 中国美白少妇内射xxxbb| av免费观看日本| 嫩草影院新地址| 99热这里只有精品一区| 欧美日韩一区二区视频在线观看视频在线 | 亚州av有码| 亚洲欧美一区二区三区国产| 亚洲无线观看免费| 九九久久精品国产亚洲av麻豆| 大香蕉97超碰在线| a级一级毛片免费在线观看| 精品人妻视频免费看| 91久久精品国产一区二区成人| 成人美女网站在线观看视频| 波多野结衣高清无吗| 你懂的网址亚洲精品在线观看 | 啦啦啦韩国在线观看视频| 男人舔奶头视频| 我要看日韩黄色一级片| 最近2019中文字幕mv第一页| 69av精品久久久久久| 亚洲成人久久爱视频| 国产精品人妻久久久久久| 久久99热这里只有精品18| 国产精品久久久久久久久免| 精华霜和精华液先用哪个| 三级国产精品片| 日韩av在线免费看完整版不卡| 久久久国产成人免费| 日本免费在线观看一区| 99久久精品国产国产毛片| 三级经典国产精品| 天美传媒精品一区二区| 国产免费男女视频| 国产高清有码在线观看视频| 日本熟妇午夜| 一区二区三区四区激情视频| 亚洲最大成人中文| 日日摸夜夜添夜夜爱| 亚洲欧美日韩东京热| 在线观看美女被高潮喷水网站| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 午夜视频国产福利| 亚洲最大成人av| 国产在线男女| 国产高清国产精品国产三级 | 亚州av有码| 中文字幕人妻熟人妻熟丝袜美| 少妇丰满av| 久久精品久久久久久久性| 亚洲精品亚洲一区二区| 麻豆精品久久久久久蜜桃| 神马国产精品三级电影在线观看| a级毛色黄片| 午夜福利在线在线| 色视频www国产| av在线观看视频网站免费| 非洲黑人性xxxx精品又粗又长| 国产乱人视频| 少妇熟女欧美另类| 91aial.com中文字幕在线观看| 日本黄色片子视频| 久久99蜜桃精品久久| 99热全是精品| av播播在线观看一区| 韩国高清视频一区二区三区| 久久99蜜桃精品久久| 色视频www国产| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 国产真实伦视频高清在线观看| 成人鲁丝片一二三区免费| 黄片无遮挡物在线观看| 国产69精品久久久久777片| 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区| 久久精品人妻少妇| 麻豆av噜噜一区二区三区| 亚洲自拍偷在线| 久久久国产成人精品二区| 精品人妻视频免费看| 免费看日本二区| 婷婷色综合大香蕉| av线在线观看网站| 免费电影在线观看免费观看| 亚洲av成人精品一区久久| 亚洲av免费在线观看| 老女人水多毛片| 亚州av有码| 国产免费又黄又爽又色| 97在线视频观看| 欧美高清成人免费视频www| 我要看日韩黄色一级片| 国产成人freesex在线| 国产黄片视频在线免费观看| 国产精品1区2区在线观看.| 97热精品久久久久久| 插逼视频在线观看| av女优亚洲男人天堂| 男的添女的下面高潮视频| 最近2019中文字幕mv第一页| av播播在线观看一区| 国产单亲对白刺激| 国产精品av视频在线免费观看| 久久精品国产鲁丝片午夜精品| 日日撸夜夜添| 联通29元200g的流量卡| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜添av毛片| 免费观看性生交大片5| 精品酒店卫生间| 日日摸夜夜添夜夜爱| 免费av不卡在线播放| 成人无遮挡网站| 久久精品91蜜桃| 国产片特级美女逼逼视频| 中文字幕av成人在线电影| 日韩一区二区三区影片| 黄片wwwwww| 桃色一区二区三区在线观看| 熟女电影av网| 亚洲精品久久久久久婷婷小说 | 欧美精品国产亚洲| 免费看日本二区| 亚洲第一区二区三区不卡| 嫩草影院新地址| 国产精品人妻久久久久久| 在线免费观看的www视频| 国产精品久久久久久久久免| 男女啪啪激烈高潮av片| 免费观看的影片在线观看| 大香蕉久久网| 美女xxoo啪啪120秒动态图| 成人欧美大片| 亚洲综合色惰| 成年av动漫网址| 亚洲精品日韩av片在线观看| 三级毛片av免费| 大话2 男鬼变身卡| videossex国产| 国产精品乱码一区二三区的特点| 亚洲精品国产成人久久av| 国产一区有黄有色的免费视频 | 日本五十路高清| 国产69精品久久久久777片| 久久精品国产自在天天线| av黄色大香蕉| 男女边吃奶边做爰视频| 熟妇人妻久久中文字幕3abv| 嫩草影院入口| 国产精品一区二区性色av| 午夜爱爱视频在线播放| 深夜a级毛片| 精品久久久久久成人av| 亚洲18禁久久av| 人人妻人人看人人澡| 我要搜黄色片| 亚洲成人av在线免费| 久久久久网色| 亚洲自拍偷在线| 内地一区二区视频在线| 国产亚洲91精品色在线| 91狼人影院| av卡一久久| 亚洲av日韩在线播放| 欧美xxxx性猛交bbbb| 精品人妻视频免费看| 久久精品国产鲁丝片午夜精品| 国产高清视频在线观看网站| 中文字幕av在线有码专区| 日本爱情动作片www.在线观看| 成人亚洲欧美一区二区av| 午夜精品在线福利| 国产精品福利在线免费观看| 2021天堂中文幕一二区在线观| 在线观看66精品国产| 我要看日韩黄色一级片| 波多野结衣巨乳人妻| 亚洲精品色激情综合| 亚洲精品456在线播放app| 日韩欧美 国产精品| 春色校园在线视频观看| 久久精品熟女亚洲av麻豆精品 | 熟妇人妻久久中文字幕3abv| 国产亚洲最大av| 91精品国产九色| 国产精品无大码| 能在线免费观看的黄片| 免费观看在线日韩| 丝袜美腿在线中文| av专区在线播放| 国产精品人妻久久久影院| 99久久成人亚洲精品观看| 国产 一区精品| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 好男人在线观看高清免费视频| 伊人久久精品亚洲午夜| 日韩强制内射视频| 赤兔流量卡办理| av又黄又爽大尺度在线免费看 | 亚洲欧美精品自产自拍| 1024手机看黄色片| 日韩av不卡免费在线播放| 午夜激情福利司机影院| 波多野结衣高清无吗| 韩国av在线不卡| 日产精品乱码卡一卡2卡三| 亚洲人与动物交配视频| 寂寞人妻少妇视频99o| 岛国在线免费视频观看| 国产探花极品一区二区| 激情 狠狠 欧美| 国产伦在线观看视频一区| 一级二级三级毛片免费看| 国产精品99久久久久久久久| 97在线视频观看| 中国美白少妇内射xxxbb| 国产黄片视频在线免费观看| 看十八女毛片水多多多| 少妇猛男粗大的猛烈进出视频 | 国产69精品久久久久777片| 国产真实伦视频高清在线观看| 女的被弄到高潮叫床怎么办| 亚洲av成人精品一二三区| 久久精品国产鲁丝片午夜精品| 七月丁香在线播放| 亚洲在线自拍视频| 国产黄a三级三级三级人| 国产在线男女| 成人漫画全彩无遮挡| 国产免费男女视频| 两个人视频免费观看高清| 欧美激情久久久久久爽电影| 欧美成人a在线观看| 九九在线视频观看精品| 精品欧美国产一区二区三| 两个人的视频大全免费| 久久久久网色| 黄色日韩在线| 色哟哟·www| 国产亚洲av片在线观看秒播厂 | 精品久久久噜噜| 长腿黑丝高跟| 久久久欧美国产精品| 一级二级三级毛片免费看| 两性午夜刺激爽爽歪歪视频在线观看| av在线播放精品| 白带黄色成豆腐渣| 国产高清有码在线观看视频| 看免费成人av毛片| 一二三四中文在线观看免费高清| av在线亚洲专区| 中国国产av一级| 国产精品99久久久久久久久| 偷拍熟女少妇极品色| 黄色配什么色好看| 一个人观看的视频www高清免费观看| 欧美日韩综合久久久久久| 亚洲婷婷狠狠爱综合网| 午夜福利在线观看免费完整高清在| 1024手机看黄色片| 亚洲精品日韩在线中文字幕| 一级av片app| 亚洲成人中文字幕在线播放| 老司机福利观看| 久久国产乱子免费精品| 如何舔出高潮| 亚洲伊人久久精品综合 | 国产免费福利视频在线观看| 国产美女午夜福利| av在线观看视频网站免费| 人妻系列 视频| 黄片wwwwww| 最近最新中文字幕大全电影3| 黑人高潮一二区| 国产精品野战在线观看| 青春草亚洲视频在线观看| 看免费成人av毛片| 午夜免费激情av| 边亲边吃奶的免费视频| 久久久久性生活片| 日本午夜av视频| a级毛色黄片| 国产一区二区在线观看日韩| 久久久久久久久久久丰满| 97超视频在线观看视频| 国产精品电影一区二区三区| 青春草亚洲视频在线观看| 亚洲欧美日韩东京热| 久久午夜福利片| 久久精品91蜜桃| 在线免费十八禁| 精品久久久久久电影网 | 毛片一级片免费看久久久久| 久久久午夜欧美精品| 日本一本二区三区精品| 久久6这里有精品| 国产乱人视频| 久久人妻av系列| 亚洲av熟女| 久久精品人妻少妇| 国产色爽女视频免费观看| 一区二区三区四区激情视频| 亚洲精品乱码久久久久久按摩| 最近最新中文字幕大全电影3| 久久精品久久精品一区二区三区| 久久久亚洲精品成人影院| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站| 欧美精品一区二区大全| 级片在线观看| 91在线精品国自产拍蜜月| 男人舔奶头视频| 亚洲在线自拍视频| 天天躁日日操中文字幕| av播播在线观看一区| 欧美潮喷喷水| av国产久精品久网站免费入址| 深夜a级毛片| 波多野结衣巨乳人妻| 久久国内精品自在自线图片| 精品久久久久久久久亚洲| 看黄色毛片网站| 国产亚洲av嫩草精品影院| 亚洲欧洲国产日韩| 中国美白少妇内射xxxbb|