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

    青海高原耐旱蠶豆品種青海13號響應(yīng)干旱脅迫蛋白質(zhì)組學(xué)分析

    2019-01-17 01:35:52侯萬偉劉玉皎
    作物學(xué)報(bào) 2019年2期
    關(guān)鍵詞:抗旱性蠶豆組學(xué)

    李 萍 侯萬偉 劉玉皎,*

    ?

    青海高原耐旱蠶豆品種青海13號響應(yīng)干旱脅迫蛋白質(zhì)組學(xué)分析

    李 萍1,2侯萬偉1,2劉玉皎1,2,*

    1青海大學(xué)農(nóng)林科學(xué)院, 青海西寧 810000;2青海大學(xué)省部共建三江源生態(tài)與高原農(nóng)牧業(yè)國家重點(diǎn)實(shí)驗(yàn)室, 青海西寧 810000

    蛋白質(zhì)組學(xué)研究在功能基因組時代發(fā)揮著越來越重要的作用, 利用雙向電泳技術(shù)和質(zhì)譜鑒定技術(shù), 可大量研究作物逆境脅迫后蛋白質(zhì)組的變化, 增加作物響應(yīng)干旱脅迫機(jī)制的認(rèn)識和理解。為探討一種抗旱性蠶豆青海13號品種的耐旱機(jī)制, 本研究對其幼苗期進(jìn)行干旱脅迫處理, 應(yīng)用上述技術(shù)進(jìn)行差異蛋白質(zhì)組分析, 經(jīng)檢驗(yàn)發(fā)現(xiàn)32個差異表達(dá)蛋白點(diǎn), 部分呈現(xiàn)上調(diào)表達(dá), 部分呈現(xiàn)下調(diào)表達(dá), 還有7個消失蛋白點(diǎn)和1個新增蛋白點(diǎn)。采用MALDI-TOF/TOF鑒定和生物信息學(xué)分析發(fā)現(xiàn), 成功鑒定的21個蛋白點(diǎn)按其所參與的代謝途徑和生化功能可分為七大類, 參與信號轉(zhuǎn)導(dǎo)的2個, 參與自由基清除的1個, 參與防衛(wèi)反應(yīng)的1個, 參與代謝的8個, 參與蛋白加工的1個, 參與光合的5個, 未知功能蛋白3個。22 kD干旱誘導(dǎo)蛋白、應(yīng)激誘導(dǎo)蛋白、17.5 kD一級熱激蛋白、超氧化物歧化酶是與抗旱性有直接關(guān)聯(lián)的蛋白點(diǎn), 其相對表達(dá)量的上調(diào)可能是青海13號蠶豆具有較強(qiáng)抗旱性的重要原因。

    蠶豆; 干旱脅迫; 蛋白質(zhì)組; 差異表達(dá); 功能分類

    干旱是脅迫因素中對植物的生長發(fā)育影響最大的因子之一, 而對農(nóng)作物造成的損更是嚴(yán)重[1], 在眾多非生物脅迫中占首位。有學(xué)者在關(guān)于植物的抗性研究中指出在水分脅迫下植物能夠產(chǎn)生一系列包括形態(tài)、生理生化及分子生物學(xué)等多方面的變化而表現(xiàn)出抗旱能力, 這些適應(yīng)性調(diào)節(jié)反應(yīng)是植物自身的保護(hù)性選擇[2-3]。蛋白質(zhì)組學(xué)是研究生物體在特定時間或空間內(nèi)所表達(dá)的全部蛋白質(zhì)的特征, 是揭示植物與脅迫之間相互作用分子機(jī)制的一種新手段。隨著蛋白質(zhì)組學(xué)的發(fā)展和植物蛋白質(zhì)組學(xué)研究的日趨深入, 以通過研究植物在非生物逆境脅迫過程中抗逆相關(guān)蛋白質(zhì)的表達(dá)而進(jìn)一步了解逆境脅迫的傷害機(jī)制及植物的逆境適應(yīng)機(jī)制的逆境蛋白質(zhì)組學(xué)研究成為熱點(diǎn)[4–6]。雙向凝膠電泳是蛋白質(zhì)組學(xué)研究的一個重要支撐技術(shù), 結(jié)合質(zhì)譜技術(shù),仍是目前最常用和可靠的蛋白質(zhì)分離鑒定技術(shù)平臺, 涉及豆科作物逆境脅迫后蛋白質(zhì)組的變化, 屢有對獲得的差異蛋白進(jìn)行功能研究的報(bào)道。Mohammadi等[7]對干旱脅迫和PEG滲透脅迫處理的大豆幼苗研究發(fā)現(xiàn), 脅迫處理下葉片中代謝相關(guān)蛋白表達(dá)上調(diào), 能量和蛋白質(zhì)合成相關(guān)蛋白下調(diào)表達(dá), 葉片、根、下胚軸這3個器官中蛋氨酸合成酶的mRNA與蛋白水平均下調(diào), 蛋氨酸合成酶是一種干旱反應(yīng)蛋白, 其下調(diào)表達(dá)抑制了大豆幼苗的正常生長。蛋白質(zhì)組學(xué)在菜豆干旱脅迫[8]、冷脅迫[9]和紫花苜蓿[10]鹽脅迫、箭舌豌豆[11]化學(xué)脅迫等方面得以廣泛應(yīng)用, 為深入開展豆科作物抗逆蛋白質(zhì)組學(xué)研究提供了參考和借鑒。本課題組前期對水分脅迫下蠶豆葉片蛋白質(zhì)組變化初探發(fā)現(xiàn), 脅迫后的葉片差異蛋白中下調(diào)表達(dá)的占大多數(shù)且主要與脅迫防御、代謝和能量、細(xì)胞骨架以及氧化平衡有關(guān), 上調(diào)表達(dá)的相對較少, 主要參與蛋白折疊與聚集以及光合系統(tǒng)[12]。

    蠶豆(L.)屬蝶形花科野豌豆屬, 是世界上重要的豆科作物, 也是青海高原重要的糧食、蔬菜、副食、飼料和養(yǎng)地兼有作物, 在該省農(nóng)業(yè)經(jīng)濟(jì)生產(chǎn)和生態(tài)保護(hù)中占有重要的地位。本研究采用人工控水模擬干旱脅迫, 研究葉片蛋白質(zhì)水平的響應(yīng), 發(fā)掘與蠶豆抵御干旱脅迫相關(guān)蛋白, 為后續(xù)的抗旱基因克隆和抗逆分子育種工作提供有力技術(shù)支持和理論基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料及處理

    青海13號是農(nóng)業(yè)生產(chǎn)中表現(xiàn)強(qiáng)抗旱性的蠶豆品種, 由青海省農(nóng)林科學(xué)院提供。挑選色澤正常無病斑的種子, 經(jīng)10%次氯酸溶液消毒及蒸餾水漂洗, 置玻璃培養(yǎng)皿中25℃恒溫光照培養(yǎng)催芽, 24 h后挑選發(fā)芽一致的蠶豆種子種植在抗旱棚中, 播種4周后待蠶豆幼苗長至二葉一心時進(jìn)行不同梯度人工控水模擬干旱脅迫處理。第1天, 對第10排蠶豆(DS9)停止供水, 由其蒸騰作用開始進(jìn)行自然干旱; 第2天, 對第9排蠶豆(DS8)停止供水開始自然干旱; 以此類推, 直到第2排蠶豆(DS1)開始干旱脅迫, 以正常供水為對照(CK)。

    1.2 IEF/SDS-PAGE雙向電泳技術(shù)及質(zhì)譜分析

    參照改良后的TCA/丙酮法[13], 略有改進(jìn), 提取葉片總蛋白。第一向IEF采用17 cm, pH 3~10 (NL)固相膠條, 上樣量900 μg; 第二向采用12%的聚丙烯酰胺分離膠進(jìn)行SDS-PAGE電泳。電泳結(jié)束后以考馬斯亮藍(lán)染色經(jīng)UMAX Power look 2100 XL型光密度掃描儀掃描, 以300 dpi分辨率采集圖像。采用PD Quest 8.0.1軟件比較處理組和對照組的葉片總蛋白2-DE圖譜, 通過差異蛋白點(diǎn)檢測, 凝膠圖譜標(biāo)準(zhǔn)化處理, 蛋白質(zhì)點(diǎn)匹配和生物統(tǒng)計(jì), 確定差異表達(dá)的蛋白質(zhì)點(diǎn), 實(shí)驗(yàn)重復(fù)3次數(shù)據(jù)進(jìn)行-test檢驗(yàn), 倍數(shù)變化>2或<0.5為差異蛋白。挖取差異表達(dá)蛋白質(zhì)點(diǎn)送上海鹿明生物有限公司進(jìn)行MALDI- TOF/TOF串聯(lián)質(zhì)譜分析, 采用ABI5800串聯(lián)飛行時間質(zhì)譜儀進(jìn)行質(zhì)譜點(diǎn)靶鑒定, 獲得的PMF的質(zhì)譜掃描范圍為800~3500 Da, 對選擇強(qiáng)度最大的10個峰進(jìn)行二級質(zhì)譜, 將一級和二級質(zhì)譜數(shù)據(jù)整合并使用GPS 3.6 (Applied Biosystems)和Mascot2.3(Matrix Science)對質(zhì)譜數(shù)據(jù)進(jìn)行分析和蛋白鑒定。進(jìn)行NCBInr-Other green plants數(shù)據(jù)庫檢索, 消化多肽的酶為胰蛋白酶, 允許最大漏切位點(diǎn)為1, 固定修飾為Carbamidomethyl (C), 可變修飾為Acetyl (Protein N-term)、Deamidated (NQ)、Dioxidation (W)和Oxidation (M), MS tolerance為0.0001 Da, MS/MS tolerance為0.3 Da, Protein score C.I.%大于95%為鑒定成功。

    2 結(jié)果與分析

    2.1 干旱脅迫最適時期確定

    土壤含水量(SWC)是表征干旱脅迫程度的一個重要的正相關(guān)指標(biāo)。通過圖1可以看出, SWC呈顯著下降趨勢, 從42.67% (SWCCK)下降到6.7% (SWCDS9),與對照相比各處理間差異極顯著, 說明本實(shí)驗(yàn)采用的不同干旱脅迫梯度設(shè)置成功, 從DS1到DS9干旱脅迫程度逐漸加強(qiáng)。其中, DS3的SWC下降迅速, 出現(xiàn)一個“折點(diǎn)”, 脅迫繼續(xù)增加, DS4和DS5與DS6至DS9脅迫間差異不顯著。由此可見, 可在后續(xù)蠶豆抗旱性研究中直接選擇DS3作為干旱脅迫的最適處理期。植物葉片相對含水量(RWC)是衡量植物持水能力強(qiáng)弱的重要指標(biāo), 在一定程度上體現(xiàn)植物抗旱能力。由圖1可以看出, 隨著脅迫時間的延長, 蠶豆葉片RWC呈下降趨勢, 但與正常供水對照相比, 直到干旱脅迫DS4時其RWC含量才與CK差異顯著。值得注意的是, 在干旱脅迫處理下, SWC在DS3處理期極顯著減小(<0.01)的情況下蠶豆葉片的RWC與對照相比在0.05水平上無顯著差異, 這也是實(shí)驗(yàn)材料青海13號具有抗旱性的一個生理表現(xiàn), 選擇此材料非常具有代表性。

    圖1 不同梯度干旱脅迫下土壤含水量(SWC)和蠶豆葉片的相對含水量(RWC)變化

    大寫字母表示不同脅迫處理間0.01水平上顯著性差異, 小寫字母表示0.05水平上差異顯著。DS1~DS9為不同梯度干旱脅迫處理, 從DS1到DS9干旱脅迫程度逐漸加強(qiáng)。

    The letter denotes differences between stress degree: the same letters indicate significant difference at the 0.05 probability level; Capital letters indicate significant difference at the 0.01 probability level. DS1–DS9 are drought stress treatments and the drought stress degrees gradually strengthened from DS1 to DS9.

    2.2 雙向電泳圖譜構(gòu)建和質(zhì)譜分析

    DS3脅迫處理是適宜的脅迫條件, 故取此處理時期的蠶豆葉片進(jìn)行下一步分析。利用本課題組前期建立的改良后的TCA/丙酮法分別提取對照組和處理組葉片總蛋白進(jìn)行雙向電泳, 經(jīng)考馬斯亮藍(lán)染色后獲得最終電泳圖譜(圖2)。

    圖2 干旱脅迫下蠶豆葉片蛋白質(zhì)組雙向電泳圖譜

    應(yīng)用PD Quest8.0.1軟件對凝膠圖譜進(jìn)行分析, 結(jié)合人工篩選,-檢測, 共有38個蛋白點(diǎn)的表達(dá)發(fā)生顯著變化(Fold>2.0或<0.5,<0.05)。處理組相對于正常對照組有10 個差異蛋白點(diǎn)的相對表達(dá)量上調(diào),占總差異蛋白點(diǎn)的26.3%, 其中特異表達(dá)的差異蛋白1個(3809)。而28個差異蛋白點(diǎn)的相對表達(dá)量下調(diào), 占總差異蛋白點(diǎn)的73.7%, 其中消失的蛋白點(diǎn)7個。挖取蠶豆葉片2-DE凝膠對差異蛋白點(diǎn)進(jìn)行MALDI-TOF-TOF質(zhì)譜分析, 通過MASCOT數(shù)據(jù)庫對獲得的數(shù)據(jù)信息進(jìn)行搜索比對, Mowse Score大于65分的蛋白即為鑒定的蛋白。通過分析成功鑒定了38個差異蛋白點(diǎn)中的27個, 鑒定率為71.1%, 這些成功鑒定的蛋白點(diǎn)的具體信息列于表1。鑒定結(jié)果表明,有多個蛋白點(diǎn)被鑒定為一個相同蛋白, 如葉綠素/結(jié)合蛋白AB80 (gi|115788: 蛋白點(diǎn)2036和2309), ATP合酶CF1 (gi|528749836: 蛋白點(diǎn)3815和3816), 質(zhì)體轉(zhuǎn)酮醇酶(gi|357445031: 蛋白點(diǎn)5808和5809), 碳酸酐酶葉綠體X1 (gi|502090577: 蛋白點(diǎn)6303、7307和7311), 核酮糖-二磷酸羧化小鏈3C (gi|132097: 蛋白點(diǎn)7005和8004), 最后27個成功鑒定的蛋白代表了21個非冗余的蛋白質(zhì)。

    2.3 差異蛋白的功能及生物信息學(xué)分析

    根據(jù)鑒定出來的21個差異蛋白所參與的代謝途徑和生化功能將其分為7類(圖3)。(1)參與信號轉(zhuǎn)導(dǎo)的2個蛋白, 占所有蛋白的10%, 分別為22 kD干旱誘導(dǎo)蛋白(drought inducible 22 kD protein, spot 2705)和應(yīng)激誘導(dǎo)蛋白(stress-inducible protein, spot 5807)。(2)參與自由基清除的1個蛋白, 超氧化物歧化酶(superoxide dismutase, spot 3106)。(3)參與防衛(wèi)反應(yīng)的1個蛋白, 17.5 kD一級熱激蛋白(17.5 kD classⅠHSPza, partial, spot 3208)。(4)參與代謝的8個蛋白占到所有蛋白的38%, 屬于最大一類, 分別為莽草酸酯脫氫酶(Shikimate dehydrogenase substrate binding, N-terminal, partial (ISS), spot 3809)、乙酰絲氨酸裂解酶(0-acetylserin(thiol)lyase, spot 5605)、胰蛋白酶抑制劑(Kunitz-type trypsin inhibitor-like 2 protein, spot 7206)、1,4-前體D-麥芽糖水解酶(1,4-alpha- D-glucan maltohydrolase, spot 1815)、ATP合酶CF1 (ATP synthase CF1 alpha subunit, spot 3815)、S-腺苷甲硫氨酸合酶(S-adenosylmethionine synthase 2, spot 4707)、質(zhì)體轉(zhuǎn)酮醇酶(plastid transketolase, spot 5808)和脫鹵素酶水解酶(haloacid dehalogenase-like hydrolase domain-containing protein At3g48420, spot 1404)。(5)參與蛋白加工的1個蛋白, nascent肽-關(guān)聯(lián)復(fù)合體亞基蛋白I (nascent polypeptide-associated complex subunit alpha-like protein 1, spot 408)。(6)參與光合的5個蛋白, 分別是葉綠體/結(jié)合蛋白AB80 (chlorophyll/binding protein AB80, spot 2306)、葉綠體/結(jié)合蛋白8 (chlorophyll a/b binding protein 8, chloroplastic-like, spot 5314)、RuBisCO小亞基3C (Ribulose bisphosphate carboxylase small chain 3C, spot 7005)、碳酸酐酶、葉綠體同型X1 (carbonic anhydrase, chloroplastic isoform X1, spot 7307)和光合系統(tǒng)II修補(bǔ)蛋白PSB27-H1 (photosystem II repair protein PSB27-H1, spot 8010)。(7)未知功能蛋白3個, 假定蛋白TSUD-183880 (hypothetical protein TSUD_183880, spot 2813)、未知蛋白(gi|217071344, spot 5202)和未命名的蛋白(unnamed protein product, spot 7608)。

    基于主流的數(shù)據(jù)庫David 6.7和QuickGO (http://www.ebi.ac.uk/QuickGO/)將21個差異表達(dá)蛋白點(diǎn)進(jìn)行質(zhì)譜搜庫檢索得到ID, 同源映射到模式植物擬南芥上進(jìn)行GO (gene ontology)富集差異表達(dá)蛋白生物信息學(xué)分析(圖4)。從圖中可以看出, 這21個差異蛋白在生物過程富集分析中參與了20個生物學(xué)過程, 主要參與光合系統(tǒng)、有機(jī)化合物代謝、有機(jī)物生物合成、響應(yīng)非生物刺激等生物過程; 在細(xì)胞定位富集中可以定位到類囊體、細(xì)胞質(zhì)、細(xì)胞器部分、胞內(nèi)細(xì)胞器、核糖核蛋白復(fù)合物、外部包裝結(jié)構(gòu)、蛋白酶體核心復(fù)合物7個細(xì)胞中, 其中在細(xì)胞質(zhì)、細(xì)胞器部分及胞內(nèi)細(xì)胞器富集蛋白百分比均超過50%; 在分子功能富集分析中, 21個差異蛋白在未折疊蛋白、前蛋白、核苷磷酸、核苷酸、核苷、核糖核苷酸、同類蛋白、幾丁質(zhì)、熱休克蛋白等的結(jié)合及轉(zhuǎn)移酶活性、陰離子配位和異構(gòu)酶活性中發(fā)揮功能, 其中參與核苷磷酸結(jié)合和核苷酸結(jié)合的蛋白富集百分比均為38.2%, 參與核苷結(jié)合和核糖核苷酸結(jié)合的蛋白富集百分比均為31.1%, 而參與陰離子配位的蛋白富集百分比為34.6%。

    3 討論

    植物在遭受干旱等逆境脅迫后基因表達(dá)發(fā)生改變,表現(xiàn)出原來一些蛋白的合成受到抑制, 也會新合成一些蛋白[14-15], 它們是植物抗旱一系列生理生化反應(yīng)的基礎(chǔ), 因此分析脅迫前后蛋白組的變化對從蛋白質(zhì)水平了解植物的抗逆機(jī)制是很有幫助的[16-18]。比起蛋白的增加或新蛋白的出現(xiàn), 干旱更容易引起蛋白含量的減少或丟失[19], 本研究中雙向電泳分析結(jié)果也得出了相同的結(jié)論, 在分離到的38個差異蛋白點(diǎn)中, 下調(diào)表達(dá)蛋白數(shù)量高于上調(diào)表達(dá)蛋白, 其中消失的蛋白有7個(spot 1504、spot 5202、spot 5314、spot 5814、spot 7311、spot 7808、spot 8102),說明干旱脅迫抑制了蠶豆植物體某些蛋白的合成來提高對干旱環(huán)境的適應(yīng)性, 亦或提高一些蛋白的表達(dá)量(8個上調(diào)表達(dá)蛋白), 同時又誘導(dǎo)產(chǎn)生一些新的蛋白如莽草酸之脫氫酶(Shikimate dehydrogenase substrate binding, N-terminal, partial (ISS), spot 3809), 從而提高蠶豆植株對干旱的耐受能力。利用MALDI-TOF/TOF串聯(lián)質(zhì)譜分析技術(shù)成功鑒定了21個非冗余的干旱響應(yīng)相關(guān)差異蛋白, 涉及信號轉(zhuǎn)導(dǎo)、防衛(wèi)反應(yīng)、光合作用、能量代謝、自由基清除、蛋白加工等過程, 還有部分獲得的蛋白為未知功能蛋白。

    圖4 差異蛋白GO功能分類

    本研究中參與信號轉(zhuǎn)導(dǎo)的Drought inducible 22 kD protein (spot 2705)和stress-inducible protein (spot 5807)等蛋白表達(dá)豐度均上調(diào), 推測蛋白酶體的增加有利于清除在干旱脅迫條件下產(chǎn)生的非正常蛋白。Desclos等[20]發(fā)現(xiàn)油菜干旱誘導(dǎo)的22 kD蛋白同時具有水溶性葉綠素結(jié)合蛋白活性和胰蛋白酶抑制劑活性, 它能通過保持蛋白的完整性和正常的光合作用保護(hù)油菜幼嫩組織不受逆境的影響。這說明Drought inducible 22 kD protein和Stress-inducible protein表達(dá)量的上升對青海13號蠶豆的抗旱性發(fā)揮了重要的作用。在植物應(yīng)對非生物脅迫中活性氧(ROS)傷害植物的膜系統(tǒng)和體內(nèi)大分子, 是新陳代謝的有毒副產(chǎn)物, 而超氧化物歧化酶(SOD)是可解除脅迫條件下植物體內(nèi)產(chǎn)生的多余活性氧分子的一類重要的活性氧清除酶, 此類表達(dá)量的增加使植物減少或免除脅迫對植物的傷害, 提高其抗逆性[21-22]。我們發(fā)現(xiàn), 干旱脅迫下超氧化物歧化酶(superoxide dismutase, spot 3106)被誘導(dǎo)表達(dá), 這與豌豆[23]、水稻[24]等作物中的研究結(jié)果一致, 是青海13號蠶豆品種對干旱脅迫的一種應(yīng)激性生理響應(yīng), 推測這也是青海13號具有強(qiáng)抗旱性的原因所在。熱激蛋白(heat shock protein, HSP)是干旱脅迫物質(zhì)形成的必需因子, 表達(dá)量的升高是植物適應(yīng)脅迫的一種機(jī)制,有利于維持脅迫調(diào)節(jié)下蛋白的穩(wěn)定性和維持植物生存與生長, 可以提高植株的抗逆能力, 在脅迫中的作用己經(jīng)被廣泛地研究[25-26]。研究者Sato和Yokoya[27]發(fā)現(xiàn)超表達(dá)sHSP17.7基因可提高干旱處理的水稻植株抗旱性,在水稻葉片蛋白質(zhì)組分析中也發(fā)現(xiàn)干旱脅迫調(diào)節(jié)多個熱激蛋白上調(diào)表達(dá)[28-29]。本研究中發(fā)現(xiàn)干旱脅迫使17.5 kD一級熱激蛋白(spot 3208)表達(dá)上調(diào), 這說明青海13號蠶豆在干旱脅迫下通過誘導(dǎo)熱激蛋白的增加來對逆境適應(yīng)起保護(hù)作用這也是該品種具有抗旱性的原因所在。氧乙酰絲氨酸裂解酶(0-acetylserin(thiol)lyase)是半胱氨酸合成最后一步的關(guān)鍵酶, 在半胱氨酸的合成調(diào)節(jié)中占據(jù)核心地位, 植物體內(nèi)的一些含硫化合物(如谷胱甘肽)可通過一些生化反應(yīng)途徑淬滅這些游離基團(tuán), 從而提高植物體的抗逆性, 我們在干旱脅迫處理下發(fā)現(xiàn)了差異上調(diào)表達(dá)蛋白中包括此蛋白(spot 5605)。干旱誘導(dǎo)蛋白是指植物在受到干旱脅迫時新合成或合成增多的一類蛋白, 干旱誘導(dǎo)蛋白的形成是植物抵御干旱脅迫的主動保護(hù)機(jī)制[30]。本研究發(fā)現(xiàn)在分離鑒定的眾多差異蛋白中, 莽草酸酯脫氫酶(shikimate dehydrogenase substrate binding, N-terminal, partial (ISS), spot 3809)是新產(chǎn)生的一類干旱誘導(dǎo)蛋白。莽草酸途徑是存在于植物中的一條重要的代謝途徑, 與植物應(yīng)對環(huán)境刺激相關(guān)聯(lián), 而莽草酸酯脫氫酶是一類促進(jìn)莽草酸途徑的關(guān)鍵性酶之一[31–32], 我們推測此干旱誘導(dǎo)蛋白對青海13號蠶豆的抗旱性有貢獻(xiàn)。植物抵抗環(huán)境脅迫的一個重要機(jī)制是調(diào)節(jié)基因的表達(dá), 嚴(yán)順平[24]在水稻抗逆性研究中推測新生肽結(jié)合蛋白復(fù)合物的α亞基(α-NAC)可能參與了水稻抗逆的基因表達(dá)的調(diào)節(jié)。也有證據(jù)表明新生肽結(jié)合復(fù)合物的α亞基可以作為轉(zhuǎn)錄的激活子[33], 在我們的研究中α-NAC (spot 408)在干旱脅迫下表達(dá)量上調(diào), 這可能影響了整個NAC的功能, 維持脅迫下基因轉(zhuǎn)錄、蛋白質(zhì)翻譯和定位的正常進(jìn)行, 保持了蠶豆正常的生理功能, 表現(xiàn)出了抗旱性。

    光合作用是綠色植物產(chǎn)生糖類等各種有機(jī)物的起始反應(yīng),是蠶豆產(chǎn)量構(gòu)成的唯一來源, 對于干旱、鹽漬等很多非生物脅迫來說都是一個非常敏感的生物學(xué)過程[34-35], 所以在本研究中發(fā)現(xiàn)干旱脅迫處理引起許多參與光合作用的蛋白表達(dá)發(fā)生變化并不奇怪。我們通過蛋白質(zhì)組學(xué)的分析發(fā)現(xiàn),很多參與光合作用的蛋白在干旱脅迫下被降解或者差異表達(dá), 這也表明了光合作用對干旱脅迫信號也是非常敏感的。Chlorophyll/binding protein AB80 (spot 2306)、chlorophyll/binding protein 8 (spot 5314)和photosystem II repair protein PSB27-H1 (spot 8010)都是與光合系統(tǒng)密切相關(guān)的蛋白, 蛋白表達(dá)的下降也驗(yàn)證了在干旱脅迫下, 植物葉片的光合作用受到抑制。ATP合酶CFI亞基(spot 3815)是與能量代謝相關(guān)的酶, 這些ATP合酶亞基在干旱脅迫下被降解將不可避免地影響光合磷酸化產(chǎn)生ATP進(jìn)一步影響光合作用的卡爾文循環(huán)。在干旱脅迫條件下, 植物通過改變體內(nèi)的初級代謝如碳的能量代謝使其達(dá)到新的平衡狀態(tài)[36], 本研究發(fā)現(xiàn)參與糖代謝途徑的磷酸戊糖途徑的關(guān)鍵酶, 轉(zhuǎn)酮醇酶(transketolase, TK, spot 5808)在干旱脅迫下表達(dá)量下調(diào), 說明干旱脅迫下磷酸戊糖途徑受到抑制。S-腺苷甲硫氨酸合酶2 (SAMS2)是植物合成木質(zhì)素的關(guān)鍵酶, 其下降表達(dá)表明木質(zhì)素以及其他相關(guān)的次生代謝產(chǎn)物的合成下降, 以利于減少不必要的能量和物質(zhì)的消耗。本研究中干旱脅迫下SAMS2 (spot 4707)蛋白表達(dá)下調(diào), 前人在研究水稻響應(yīng)鹽脅迫和干旱脅迫的蛋白質(zhì)組變化時也發(fā)現(xiàn)SAMS2的蛋白水平被下調(diào)[37], 推測此蛋白也與青海13號蠶豆抗旱性有關(guān)聯(lián)。

    4 結(jié)論

    在干旱脅迫下, 青海 13 號蠶豆葉片與信號轉(zhuǎn)導(dǎo)、自由基清除、防衛(wèi)反應(yīng)、能量代謝及蛋白加工等相關(guān)的蛋白主要表現(xiàn)為上調(diào)表達(dá)中發(fā)現(xiàn)了一個新的干旱脅迫響應(yīng)蛋白莽草酸酯脫氫酶(shikimate dehydrogenase substrate binding, N-terminal, partial (ISS), spot 3809), 推測是青海13號蠶豆具有抗旱性機(jī)制所在。這些為我們深入了解蠶豆對干旱脅迫的響應(yīng)提供了新的線索, 也是利用遺傳等手段進(jìn)一步研究這些基因功能的開始。

    [1] Pandey A, Chakraborty S, Datta A, Chakraborty N. Proteomics approach to identify dehydration responsive nuclear proteins from chickpea (L.)., 2008, 7: 88–107.

    [2] Marie P R, Vartavian N, Sallantin M. Characterization of a novel protein induced by progressive or rapid drought and salinity inleaves., 1992, 100: 1486–1493.

    [3] Ingram J, Bartels D. The molecular basis of dehydration tolerance in plant., 1996, 47: 377–403.

    [4] Levenfors J P, Wikstr M M, Persson L, Gerhardson B. Pathogenicity of aphanomyces spp from different leguminous crops in Sweden., 2003, 10: 535–543.

    [5] Agrawal G K, Pedreschi R, Barkla B J, Bindschedler L V, Cramer R, Sarkar A, Renaut J, Job D, Rakwal R. Translational plant proteomics: A perspective., 2012, 75: 4588–4601.

    [6] DO Thanh-Trung, 李健, 張風(fēng)娟, 楊麗濤, 李楊瑞, 邢永秀. 甘蔗與抗旱性相關(guān)差異蛋白質(zhì)組分析. 作物學(xué)報(bào), 2017, 43: 1337–1346. DO T T, Li J, Zhang F J, Yang L T, Li Y R, Xing Y X. Analysis of differential proteome in relation to drought resistance in Sugarcane., 2017, 43: 1337–1346 (in Chinese with English abstract).

    [7] Mohammadi P P, Moieni A, Hiraga S, Komatsu S. Organ-specific proteomic analysis of drought-stressed soybean seedlings., 2012, 75: 1906–1923.

    [8] Zadra?nik T, Hollung K, Egge-Jacobsen W. Differential proteomic analysis of drought stress response in leaves of common bean (L.)., 2013, 78: 254–272.

    [9] Badowiec A, Weidner S. Proteomic changes in the roots of germinatingseeds in response to chilling stress and post-stress recovery., 2014, 2: 389–398.

    [10] 熊軍波, 楊青川, 蔡化, 田宏, 張賀山, 劉洋. 紫花苜蓿根響應(yīng)鹽脅迫的比較蛋白質(zhì)組學(xué)分析. 湖北農(nóng)業(yè)科學(xué), 2015, 54: 5422–5428. Xiong J B, Yang Q H, Cai H, Tian H, Zhang H S, Li Y. Comparative proteomic analysis of salt-stress response of alfalfa proteins in Root., 2015, 54: 5422–5428 (in Chinese with English abstract).

    [11] 芮海云, 莊凱, 沈振國, 張芬琴. 兩個箭舌豌豆品種根響應(yīng)鎘脅迫的蛋白質(zhì)組學(xué)分析. 植物生理學(xué)報(bào), 2016, 52: 1089–1098. Rui H Y, Zhuang K, Sheng Z G, Zhang F Q. Proteomics analysis of cadmium stress responses in the roots of twovarieties differing in Cd tolerance., 2016, 52: 1089–1098 (in Chinese with English abstract).

    [12] 李萍, 張雁霞, 劉玉皎. 干旱脅迫對抗旱性蠶豆幼苗生長特性影響及葉片差異蛋白質(zhì)組學(xué)研究. 基因組學(xué)與應(yīng)用生物學(xué), 2017, 36: 340–351. Li P, Zhang Y X, Liu Y J. Effects of drought stress on seeding growth characteristics of drought resistance faba bean (L.) and study on differential proteomics of leaves., 2017, 36: 340–351 (in Chinese with English abstract).

    [13] 張小靜. 馬鈴薯塊莖蛋白質(zhì)雙向電泳技術(shù)體系建立及發(fā)育相關(guān)蛋白質(zhì)的分析. 甘肅農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 甘肅蘭州, 2008.Zhang X J. Establishment of Two-dimensional Electrophoresis System for Analysis of Proteins during Potato (L.) Tuber Development. MS Thesis of Gansu Agricultural University, Lanzhou, Gansu, China, 2008 (in Chinese with English abstract).

    [14] Demirevska K, Zasheva D, Dimitrov R, Simova-Stoilova L, Stamenova M, Feller U. Drought stress effects on rubisco in wheat: changes in the rubisco large subunit., 2009, 31: 1129–1138.

    [15] Salekdeh G H, Siopongco J, Wade L J, Ghareyazie B, Benett J. Proteomic analysis of rice leaves during drought stress and recovery., 2002, 2: 1131–1145.

    [16] Mohammadi P P, Moieni A, Hiraga S. Organ-specific proleomic analysis of drought-stressed soybean seedling., 2012, 75: 1960–1923.

    [17] Mustafa G, Komatsu S. Quantitative proteomics reveals the effect of protein glycosylation in soybean root under flooding stress., 2014, 5: 627.

    [18] Ma H Y, Song L R, Huang Z G. Comparative proteomic analysis reveals molecular mechanism of seedling roots of different salt tolerant soybean genotypes in responses to salinity stress., 2014, 4: 40–57.

    [19] Celis J E, Gromov P. 2D Protein electrophoresis: can it be perfected?, 1999, 10: 16–21.

    [20] Desclos M, Dubousset L, Etienne P, Le Caherec F, Satoh H, Bonnefoy J, Ourry A, Avice J C. A proteomic profiling approach to reveal a novel role ofdrought 22 kD/water- soluble chlorophyll-binding protein in young leaves during nitrogen remobilization induced by stressful conditons., 2008, 147: 1830–1844.

    [21] Balestrasse K B, Gardey L, Gallego S M. Response of antioxidant defence systerm in soybean nodules and roots subjected to cadmium stress., 2001, 28: 497–453.

    [22] Malgorzata G, Waldemar B. Effects of a short-term hypoxic treatment followed by re-aeration on free radicals level and anti-oxidative enzymes in lupine roots., 2004, 42: 233–240.

    [23] Hernandez J A, Jimenez A, Mullineaux P, Sevilla F. Tolerance of pea (L.) to long-term salt stress is associated with induction of antioxidant defences., 2000, 23: 853–862.

    [24] 嚴(yán)順平. 水稻響應(yīng)鹽脅迫和低溫脅迫的蛋白質(zhì)組研究. 中國科學(xué)院博士學(xué)位論文, 上海, 2006. Yan S P. Proteomic Analysis of Salt Stress Responsive and Chilling Stress Response in Rice. PhD Dissertation of Graduate University of the Chinese Academy of Sciences, Shanghai, China, 2006 (in Chinese with English abstract).

    [25] Hajheidari M, Abdollahian-Noghabi M, Askari H, Heidari M, Sadeghian S Y, Ober E S, Salekdeh G H. Proteome analysis of sugar beet leaves under drought stress., 2005, 5: 950–960.

    [26] Wang W, Vinocur B, Shoseyov O, Altman A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response., 2004, 9: 244–252.

    [27] Sato Y, Yokoya S. Enhanced to tolerance to drought stress in transgenic rice plants overexpressing a small heat-shock protein, sHSP17.7., 2008, 27: 329–334.

    [28] 丁偉. 水稻干旱脅迫蛋白質(zhì)組相關(guān)數(shù)據(jù)和生物信息分析研究. 華中農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 湖北武漢, 2009. Ding W. Analysis of Bioinformation and Data Related to Proteome of Rice under Drought Stress. MS Thesis of Huazhong Agricultural University, Wuhan, Hubei, China, 2009 (in Chinese with English abstract).

    [29] 舒烈波. 水稻葉片響應(yīng)干旱和滲透脅迫的蛋白質(zhì)組學(xué)研究. 華中農(nóng)業(yè)大學(xué)博士學(xué)位論文, 湖北武漢, 2010. Shu L B. Proteomic Analysis of Rice Leaves in Response to Drought and Osmotic Stress. PhD Dissertation of Huazhong Agricultural University, Wuhan, Hubei, China, 2010 (in Chinese with English abstract).

    [30] 李麗芳, 羅曉芳, 王華芳. 植物抗旱基因工程研究進(jìn)展. 西北林學(xué)院學(xué)報(bào), 2004, 19: 53–57. Li L F, Luo X F, Wang H F. Advances in the studies of gene engineering on plant drought-resistance., 2004, 19: 53–57 (in Chinese with English abstract).

    [31] Weaver L M, Herrmann K M. Dynamics of the shikimate pathway in plants., 1997, 2: 346–351.

    [32] Diaz J, Bernal A, Pomar F, Merino F. Induction of shikimate dehydrogenase and peroxidase in pepper (L.) seedlings in response to copper stress and its relation to lignification., 2001, 161: 179–188.

    [33] Rospert S, Dubaquie Y, Gautschi M. Nascentpolypeptide- associated complex., 2002, 59: 1632–1639.

    [34] Chaves M M, Flexas J, Pinheiro C. Photosynthesis under drought and salt stress: Regulation mechanisms from whole plant to cell., 2009, 103: 551–560.

    [35] Reddy A R, Chaitanya K V, Vivekanandan M. Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants., 2004, 161: 1189–1202.

    [36] Thomashow M F. Plant cold acclimation: freezing tolerance genes and regulatory mechanisms., 1999, 50: 571–599.

    [37] Li X J, Yang M F, Chen H, Qu L Q, Chen F, Shen S H. Abscisic acid pretreatment enhances salt tolerance of rice seedlings: proteomic evidence., 2010, 1804: 929–940.

    Proteomic analysis of drought stress response on drought resistance forL. variety ‘Qinghai 13’ in Qinghai Plateau of China

    LI Ping1,2, HOU Wan-Wei1,2, and LIU Yu-Jiao1,2,*

    1Academy of Agriculture and Forestry Science of Qinghai University, Xining 810016, Qinghai, China;2State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining 810016, Qinghai, China

    Proteomics is playing an increasingly important role in the functional genomics era. Two-dimensional electrophoresis and mass spectrometry can be used to study the changes of proteome of crop under stress and increase the recognization and comprehension of crop response to drought stress. In order to explore the mechanism of drought resistance ofL. variety ‘Qinghai 13’, we treated the seedlings with three days for water stress and analyzed by two-dimensional gel electrophoresis combined with mass spectrometry analysis. Bytest, 32 differentially expressed proteins spots were detected between normal and drought-stress treatments, respectively, including up- and down-regulated proteins, seven disappear protein spots and a new protein spot. Twenty-one differentially expressed proteins in seven function categories were identified and confirmed by MALDI-TOF/TOF. Among them, two participated in information transfer, one in oxygen radical scavenging, one in protective response, eight in energy metabolism, one in protein processing, five in photosynthesis, and three unknown in function. These results indicate that drought-inducible 22 kD protein, stress-inducible protein, superoxide dismutase and 17.5 kD class HSP are directly related with drought resistance, which may be the important reasons for strong drought resistance of ‘Qinghai 13’.

    faba bean; drought stress; proteomic; differential expression; functional classification

    ): 2018-05-29;

    2018-10-08;

    2018-11-03.

    10.3724/SP.J.1006.2019.84075

    劉玉皎, E-mail: 13997058356@163.com

    E-mail: lping2008@ sohu.com

    本研究由國家自然科學(xué)基金項(xiàng)目(31460377)和國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-09)資助。

    This study was supported by the National Natural Science Foundation of China (31460377) and the China Agriculture Research System (CARS-09).

    URL: http://kns.cnki.net/kcms/detail/11.1809.S.20181031.1457.004.html

    猜你喜歡
    抗旱性蠶豆組學(xué)
    蠶豆花開
    小讀者(2021年2期)2021-03-29 05:03:30
    且將蠶豆伴青梅
    口腔代謝組學(xué)研究
    蠶豆大嘴巴
    基于UHPLC-Q-TOF/MS的歸身和歸尾補(bǔ)血機(jī)制的代謝組學(xué)初步研究
    基于蠶豆產(chǎn)業(yè)鏈的研究與開發(fā)
    不同光質(zhì)對黃瓜幼苗抗旱性的影響
    基于稱重法的煙草伸根期抗旱性鑒定
    樹木抗旱性及抗旱造林技術(shù)研究綜述
    代謝組學(xué)在多囊卵巢綜合征中的應(yīng)用
    亚洲最大成人中文| 搞女人的毛片| 国产黄片美女视频| 人人妻人人看人人澡| 欧美人与善性xxx| 国产精品国产高清国产av| 看十八女毛片水多多多| 啦啦啦韩国在线观看视频| 欧美激情久久久久久爽电影| 亚洲中文字幕一区二区三区有码在线看| 大型黄色视频在线免费观看| 久久精品国产清高在天天线| 欧美成人性av电影在线观看| 18禁黄网站禁片午夜丰满| 九色成人免费人妻av| 国产亚洲欧美98| 亚洲av日韩精品久久久久久密| 日韩欧美精品v在线| 日韩精品青青久久久久久| 老熟妇仑乱视频hdxx| ponron亚洲| 久久热精品热| 三级男女做爰猛烈吃奶摸视频| 日日夜夜操网爽| 91精品国产九色| 成人av在线播放网站| 午夜福利视频1000在线观看| 亚洲av熟女| 亚洲av一区综合| 美女高潮喷水抽搐中文字幕| 少妇人妻一区二区三区视频| 亚洲av中文av极速乱 | 国产av在哪里看| ponron亚洲| 3wmmmm亚洲av在线观看| 国产av不卡久久| 亚洲av熟女| 国产v大片淫在线免费观看| 啦啦啦韩国在线观看视频| 白带黄色成豆腐渣| 国产av在哪里看| 国产精品一区二区三区四区久久| 日韩欧美精品免费久久| 一卡2卡三卡四卡精品乱码亚洲| 中文在线观看免费www的网站| 搞女人的毛片| 搡老岳熟女国产| 哪里可以看免费的av片| 欧美色视频一区免费| 亚洲国产精品久久男人天堂| 欧美日韩瑟瑟在线播放| 欧美+亚洲+日韩+国产| 久久久久久伊人网av| 制服丝袜大香蕉在线| 简卡轻食公司| 日日摸夜夜添夜夜添小说| 在线播放无遮挡| 熟女电影av网| 久久精品国产亚洲av香蕉五月| 听说在线观看完整版免费高清| 日本免费a在线| 国产欧美日韩一区二区精品| 又爽又黄a免费视频| 欧美+日韩+精品| 一区二区三区高清视频在线| 美女高潮喷水抽搐中文字幕| 国产成人影院久久av| 婷婷丁香在线五月| 国产伦一二天堂av在线观看| 午夜福利在线在线| 国产一区二区激情短视频| 午夜免费激情av| av黄色大香蕉| 国产免费一级a男人的天堂| 精品午夜福利在线看| 春色校园在线视频观看| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 丰满乱子伦码专区| 在线看三级毛片| 日日夜夜操网爽| 成人特级黄色片久久久久久久| 少妇被粗大猛烈的视频| 成人午夜高清在线视频| 亚洲精品456在线播放app | 国产精品爽爽va在线观看网站| 免费观看的影片在线观看| 中文字幕免费在线视频6| 日日摸夜夜添夜夜添av毛片 | 国产淫片久久久久久久久| 国产精品久久久久久久久免| 亚洲午夜理论影院| 精品一区二区三区人妻视频| 天堂动漫精品| 97热精品久久久久久| 欧美丝袜亚洲另类 | 香蕉av资源在线| 最近中文字幕高清免费大全6 | 免费观看精品视频网站| 在线观看66精品国产| 在线观看一区二区三区| 国产精品久久久久久av不卡| 老司机福利观看| 中国美白少妇内射xxxbb| 久久精品人妻少妇| 日韩强制内射视频| 韩国av在线不卡| 午夜精品一区二区三区免费看| 国产视频一区二区在线看| 成人国产一区最新在线观看| 亚洲熟妇熟女久久| 熟女人妻精品中文字幕| 国产亚洲精品综合一区在线观看| 午夜免费成人在线视频| 欧美性猛交╳xxx乱大交人| 亚洲人与动物交配视频| 啦啦啦韩国在线观看视频| 亚洲精品亚洲一区二区| a级毛片免费高清观看在线播放| a在线观看视频网站| 日本欧美国产在线视频| 一进一出抽搐gif免费好疼| 精品福利观看| 成人毛片a级毛片在线播放| 美女大奶头视频| 丰满的人妻完整版| 久久久久久久午夜电影| 校园人妻丝袜中文字幕| 岛国在线免费视频观看| 18禁裸乳无遮挡免费网站照片| 国产精品一区www在线观看 | 狠狠狠狠99中文字幕| 亚洲精华国产精华精| 一夜夜www| 国产精品久久视频播放| 久久精品国产鲁丝片午夜精品 | 亚洲欧美日韩卡通动漫| 久久99热6这里只有精品| 色综合婷婷激情| 成人国产一区最新在线观看| or卡值多少钱| 欧美中文日本在线观看视频| 亚洲国产精品合色在线| 日韩欧美一区二区三区在线观看| 国产一区二区在线av高清观看| 久久九九热精品免费| 亚洲天堂国产精品一区在线| 成年女人看的毛片在线观看| 99久久精品热视频| 久久国内精品自在自线图片| 91在线观看av| 免费电影在线观看免费观看| 熟女人妻精品中文字幕| h日本视频在线播放| 五月玫瑰六月丁香| 91麻豆av在线| 久久热精品热| 午夜福利在线观看吧| 国产探花在线观看一区二区| 婷婷色综合大香蕉| 国产精品嫩草影院av在线观看 | 男人舔女人下体高潮全视频| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播| 伊人久久精品亚洲午夜| 久久婷婷人人爽人人干人人爱| 日本a在线网址| 午夜福利成人在线免费观看| 精品福利观看| 日日撸夜夜添| 人妻少妇偷人精品九色| 国产伦在线观看视频一区| 午夜福利18| 一个人免费在线观看电影| 久久精品国产99精品国产亚洲性色| 国产精品99久久久久久久久| 少妇丰满av| 毛片一级片免费看久久久久 | 美女cb高潮喷水在线观看| 国产av一区在线观看免费| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 欧美xxxx性猛交bbbb| 色尼玛亚洲综合影院| 午夜久久久久精精品| 亚洲狠狠婷婷综合久久图片| 亚洲精品粉嫩美女一区| 欧美一区二区亚洲| 国模一区二区三区四区视频| 国产人妻一区二区三区在| 如何舔出高潮| 啦啦啦韩国在线观看视频| 亚洲精品亚洲一区二区| 女的被弄到高潮叫床怎么办 | 中国美白少妇内射xxxbb| 伦理电影大哥的女人| 12—13女人毛片做爰片一| 午夜福利欧美成人| 乱人视频在线观看| 一本久久中文字幕| 久久精品国产亚洲av天美| 色5月婷婷丁香| 人妻夜夜爽99麻豆av| 亚洲欧美日韩东京热| 偷拍熟女少妇极品色| 日本与韩国留学比较| 亚洲精品色激情综合| 国产三级在线视频| 男女边吃奶边做爰视频| 丝袜美腿在线中文| 精品欧美国产一区二区三| av在线天堂中文字幕| av国产免费在线观看| 人妻夜夜爽99麻豆av| 亚洲va日本ⅴa欧美va伊人久久| 尤物成人国产欧美一区二区三区| 欧美高清性xxxxhd video| 亚洲一区高清亚洲精品| 91精品国产九色| 欧美极品一区二区三区四区| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 一级a爱片免费观看的视频| 久久6这里有精品| 欧美另类亚洲清纯唯美| 亚洲av五月六月丁香网| 十八禁网站免费在线| 男女视频在线观看网站免费| 成人二区视频| 成人国产综合亚洲| 久久精品国产亚洲av香蕉五月| 久久久久国内视频| 日韩一区二区视频免费看| 久久欧美精品欧美久久欧美| 岛国在线免费视频观看| 特大巨黑吊av在线直播| 特级一级黄色大片| 久久精品国产清高在天天线| 18禁黄网站禁片午夜丰满| 热99re8久久精品国产| 桃红色精品国产亚洲av| 成人一区二区视频在线观看| 91av网一区二区| 少妇的逼好多水| 精品一区二区三区视频在线| 成人精品一区二区免费| 91午夜精品亚洲一区二区三区 | 黄色女人牲交| 乱系列少妇在线播放| 国产v大片淫在线免费观看| 国产精品人妻久久久影院| 亚洲综合色惰| 免费av毛片视频| 亚洲最大成人av| 亚洲四区av| 亚洲人成网站在线播| 成人国产麻豆网| 国产一区二区在线av高清观看| 免费av观看视频| 草草在线视频免费看| 啪啪无遮挡十八禁网站| 啦啦啦观看免费观看视频高清| 99久久久亚洲精品蜜臀av| 国产精品爽爽va在线观看网站| 亚洲男人的天堂狠狠| 一区二区三区激情视频| 极品教师在线免费播放| 精品久久久噜噜| 2021天堂中文幕一二区在线观| 18禁裸乳无遮挡免费网站照片| 中文字幕av在线有码专区| 欧美日韩精品成人综合77777| 一本久久中文字幕| 亚洲,欧美,日韩| 国产精品国产三级国产av玫瑰| 国内精品久久久久久久电影| 亚洲一区二区三区色噜噜| 亚洲av一区综合| 真实男女啪啪啪动态图| 日韩中字成人| 中文字幕精品亚洲无线码一区| 网址你懂的国产日韩在线| 国产一区二区在线av高清观看| 色尼玛亚洲综合影院| av天堂在线播放| 午夜免费激情av| 熟妇人妻久久中文字幕3abv| 国产成人福利小说| 国产又黄又爽又无遮挡在线| avwww免费| 很黄的视频免费| 日本黄色视频三级网站网址| 99国产精品一区二区蜜桃av| 尾随美女入室| a级毛片免费高清观看在线播放| av黄色大香蕉| 国产精品av视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 偷拍熟女少妇极品色| 最近最新中文字幕大全电影3| 国产精品无大码| 99在线视频只有这里精品首页| 黄色一级大片看看| 欧美性猛交╳xxx乱大交人| 亚洲人成网站高清观看| 精品一区二区免费观看| 欧美成人免费av一区二区三区| 最后的刺客免费高清国语| 91狼人影院| 国产av不卡久久| 日韩欧美在线二视频| av中文乱码字幕在线| 日韩亚洲欧美综合| 国产亚洲精品久久久久久毛片| 国产精品女同一区二区软件 | 欧美又色又爽又黄视频| 国产淫片久久久久久久久| 免费无遮挡裸体视频| 国产精品自产拍在线观看55亚洲| 国内精品久久久久精免费| 久久精品国产自在天天线| 精品无人区乱码1区二区| 亚洲第一区二区三区不卡| 岛国在线免费视频观看| 一个人看视频在线观看www免费| 亚洲av日韩精品久久久久久密| 中文字幕人妻熟人妻熟丝袜美| 欧美日本亚洲视频在线播放| 乱系列少妇在线播放| 一级黄色大片毛片| 久久6这里有精品| 在线国产一区二区在线| 亚洲国产欧美人成| 天堂网av新在线| 婷婷六月久久综合丁香| 亚洲欧美精品综合久久99| 神马国产精品三级电影在线观看| 国产成人av教育| 动漫黄色视频在线观看| 97热精品久久久久久| 日本熟妇午夜| 国产精品女同一区二区软件 | 俄罗斯特黄特色一大片| 中文字幕高清在线视频| 色哟哟·www| 国产大屁股一区二区在线视频| 18禁黄网站禁片免费观看直播| 一级黄片播放器| 成人毛片a级毛片在线播放| 国产精华一区二区三区| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 日韩,欧美,国产一区二区三区 | 午夜免费男女啪啪视频观看 | 免费观看的影片在线观看| 日韩精品中文字幕看吧| 69人妻影院| av在线亚洲专区| 国产精品三级大全| 美女被艹到高潮喷水动态| 免费不卡的大黄色大毛片视频在线观看 | 男人和女人高潮做爰伦理| 国产白丝娇喘喷水9色精品| 午夜亚洲福利在线播放| 国产精品亚洲美女久久久| 亚洲狠狠婷婷综合久久图片| 在线观看舔阴道视频| 最近视频中文字幕2019在线8| 成人一区二区视频在线观看| 九色成人免费人妻av| 国内精品久久久久精免费| 久久久久久久久大av| 色av中文字幕| 精品久久久久久久久亚洲 | 日韩欧美三级三区| 又黄又爽又刺激的免费视频.| 国产精品久久久久久精品电影| 成熟少妇高潮喷水视频| 国内揄拍国产精品人妻在线| 亚洲图色成人| 丰满的人妻完整版| 真实男女啪啪啪动态图| 天堂动漫精品| 联通29元200g的流量卡| 亚洲18禁久久av| 久久久久精品国产欧美久久久| 男女下面进入的视频免费午夜| 久久精品国产99精品国产亚洲性色| 午夜免费成人在线视频| 欧美色欧美亚洲另类二区| 大型黄色视频在线免费观看| 午夜爱爱视频在线播放| 欧美+亚洲+日韩+国产| 51国产日韩欧美| 91久久精品电影网| 成人av在线播放网站| 午夜福利高清视频| 简卡轻食公司| 亚洲成人免费电影在线观看| 一本一本综合久久| 中文亚洲av片在线观看爽| 国产成人影院久久av| 欧美三级亚洲精品| 国产欧美日韩一区二区精品| 日韩欧美精品免费久久| 国产高清不卡午夜福利| 日韩欧美精品免费久久| av.在线天堂| 国产精品嫩草影院av在线观看 | 午夜福利在线观看吧| 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 欧美日韩综合久久久久久 | 国产探花极品一区二区| 狠狠狠狠99中文字幕| 黄色一级大片看看| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 欧美色视频一区免费| 国产乱人视频| 97超级碰碰碰精品色视频在线观看| 欧美又色又爽又黄视频| 最新在线观看一区二区三区| 搡老岳熟女国产| 99久国产av精品| 亚洲成人中文字幕在线播放| 99热这里只有精品一区| 丰满的人妻完整版| 午夜精品一区二区三区免费看| 国产av不卡久久| 最近最新中文字幕大全电影3| 我的女老师完整版在线观看| 老女人水多毛片| 日韩欧美三级三区| 大又大粗又爽又黄少妇毛片口| 变态另类丝袜制服| 中文字幕av在线有码专区| 欧美性猛交╳xxx乱大交人| 特级一级黄色大片| 久久6这里有精品| 日本成人三级电影网站| 桃红色精品国产亚洲av| ponron亚洲| 免费在线观看成人毛片| 亚洲欧美日韩无卡精品| 精品99又大又爽又粗少妇毛片 | 国产欧美日韩精品一区二区| 丰满乱子伦码专区| av在线亚洲专区| 国产伦人伦偷精品视频| 国产精品一区二区三区四区免费观看 | 亚洲五月天丁香| 麻豆成人av在线观看| 嫁个100分男人电影在线观看| 国产极品精品免费视频能看的| 欧美精品国产亚洲| 99精品在免费线老司机午夜| 国产精品自产拍在线观看55亚洲| 国产免费男女视频| 国产精品久久视频播放| 国产一级毛片七仙女欲春2| 午夜精品在线福利| 日本免费一区二区三区高清不卡| 少妇人妻精品综合一区二区 | 免费av观看视频| 免费大片18禁| 精品福利观看| 国产亚洲av嫩草精品影院| 午夜福利在线观看免费完整高清在 | 日韩中字成人| 九九爱精品视频在线观看| 色综合亚洲欧美另类图片| av在线老鸭窝| 亚洲精品粉嫩美女一区| 中文资源天堂在线| 一个人观看的视频www高清免费观看| 国产精品98久久久久久宅男小说| 一本久久中文字幕| 久久久久久九九精品二区国产| 非洲黑人性xxxx精品又粗又长| 国产 一区精品| 别揉我奶头 嗯啊视频| 日本精品一区二区三区蜜桃| 小蜜桃在线观看免费完整版高清| av在线天堂中文字幕| 欧美绝顶高潮抽搐喷水| 性插视频无遮挡在线免费观看| 国产高清不卡午夜福利| 久久久国产成人免费| 波多野结衣高清无吗| 男人狂女人下面高潮的视频| 热99在线观看视频| 国产成人a区在线观看| 久久精品综合一区二区三区| 熟妇人妻久久中文字幕3abv| 色精品久久人妻99蜜桃| 国产一区二区三区视频了| 51国产日韩欧美| 亚洲欧美日韩无卡精品| 级片在线观看| 欧美性猛交黑人性爽| 亚洲欧美日韩东京热| 丝袜美腿在线中文| 精品一区二区三区视频在线观看免费| 可以在线观看的亚洲视频| 欧美日本亚洲视频在线播放| 高清在线国产一区| 精品久久国产蜜桃| 成人永久免费在线观看视频| 免费搜索国产男女视频| 精品一区二区三区视频在线观看免费| 午夜激情福利司机影院| 最新中文字幕久久久久| 亚洲欧美日韩高清在线视频| 亚洲精品国产成人久久av| 搡老熟女国产l中国老女人| 女人十人毛片免费观看3o分钟| 国产精品98久久久久久宅男小说| 国产真实乱freesex| 自拍偷自拍亚洲精品老妇| 午夜福利成人在线免费观看| 88av欧美| 亚洲精华国产精华液的使用体验 | 亚洲综合色惰| 黄色欧美视频在线观看| 最近视频中文字幕2019在线8| 欧美绝顶高潮抽搐喷水| 性色avwww在线观看| 久久热精品热| 精品日产1卡2卡| 亚洲精品色激情综合| 国产精品不卡视频一区二区| 丝袜美腿在线中文| 亚洲一区二区三区色噜噜| 干丝袜人妻中文字幕| 美女免费视频网站| 黄色丝袜av网址大全| 国产精品三级大全| 老司机深夜福利视频在线观看| 日本 欧美在线| 国内少妇人妻偷人精品xxx网站| 一卡2卡三卡四卡精品乱码亚洲| 哪里可以看免费的av片| 伦理电影大哥的女人| 欧美色视频一区免费| 精品久久久久久久久亚洲 | 我要看日韩黄色一级片| h日本视频在线播放| 18禁黄网站禁片午夜丰满| 草草在线视频免费看| 亚洲中文日韩欧美视频| 99久久无色码亚洲精品果冻| 国产三级中文精品| 黄色丝袜av网址大全| 我要看日韩黄色一级片| av.在线天堂| 久久草成人影院| 亚洲国产欧洲综合997久久,| 黄色一级大片看看| 欧美xxxx性猛交bbbb| 国产免费av片在线观看野外av| 国产亚洲精品综合一区在线观看| 亚洲va在线va天堂va国产| 午夜影院日韩av| 久久人妻av系列| 天美传媒精品一区二区| 1000部很黄的大片| 九九爱精品视频在线观看| 淫秽高清视频在线观看| 日韩中文字幕欧美一区二区| 中出人妻视频一区二区| 国产一区二区激情短视频| .国产精品久久| 99热6这里只有精品| 国产精品,欧美在线| 久久久久久久久久成人| 久久久国产成人免费| 热99在线观看视频| 亚洲国产日韩欧美精品在线观看| 春色校园在线视频观看| 美女免费视频网站| a级一级毛片免费在线观看| 毛片一级片免费看久久久久 | 亚洲国产色片| 日本 av在线| 在线观看舔阴道视频| 午夜精品在线福利| 日韩精品有码人妻一区| 赤兔流量卡办理| 亚洲欧美清纯卡通| 久久国产乱子免费精品| 窝窝影院91人妻| 97热精品久久久久久| 亚洲国产精品sss在线观看| 国内精品久久久久精免费| 男女视频在线观看网站免费| 国产男靠女视频免费网站| 精品久久久噜噜| 天天躁日日操中文字幕| 日本与韩国留学比较| 免费观看精品视频网站| 最近最新免费中文字幕在线| 国产男靠女视频免费网站| 淫妇啪啪啪对白视频| 91久久精品国产一区二区成人| 婷婷精品国产亚洲av在线| 亚洲人成网站高清观看| 免费搜索国产男女视频| 2021天堂中文幕一二区在线观| 大又大粗又爽又黄少妇毛片口| 中亚洲国语对白在线视频| 99热这里只有精品一区| 亚洲最大成人av| 动漫黄色视频在线观看|