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

    灌漿期夜間高溫脅迫下耐熱和熱敏感水稻籽粒的比較蛋白質(zhì)組分析

    2017-02-20 07:49:02黃小平張宏玉雷剛王志美章智賀超廖江林黃英金
    中國(guó)水稻科學(xué) 2017年1期
    關(guān)鍵詞:鋅指籽粒灌漿

    黃小平張宏玉雷剛王志美章智賀超廖江林,,?黃英金,,?

    (1江西農(nóng)業(yè)大學(xué)作物生理生態(tài)與遺傳育種教育部重點(diǎn)實(shí)驗(yàn)室,南昌330045;2江西農(nóng)業(yè)大學(xué)農(nóng)業(yè)應(yīng)對(duì)氣候變化南昌市重點(diǎn)實(shí)驗(yàn)室,南昌330045;?通訊聯(lián)系人,E-mail:jlliao514815@163.com,yjhuang_cn@126.com)

    灌漿期夜間高溫脅迫下耐熱和熱敏感水稻籽粒的比較蛋白質(zhì)組分析

    黃小平1張宏玉2雷剛1王志美1章智1賀超1廖江林1,2,?黃英金1,2,?

    (1江西農(nóng)業(yè)大學(xué)作物生理生態(tài)與遺傳育種教育部重點(diǎn)實(shí)驗(yàn)室,南昌330045;2江西農(nóng)業(yè)大學(xué)農(nóng)業(yè)應(yīng)對(duì)氣候變化南昌市重點(diǎn)實(shí)驗(yàn)室,南昌330045;?通訊聯(lián)系人,E-mail:jlliao514815@163.com,yjhuang_cn@126.com)

    【目的】為了闡明水稻灌漿期耐夜間高溫的分子機(jī)制,提高水稻耐熱性,鑒定了水稻近等基因系耐熱純系XN0437T與熱敏感純系XN0437S在灌漿期夜間高溫脅迫下的差異表達(dá)蛋白質(zhì)?!痉椒ā坎捎猛霸苑椒ㄅ嘤?于開花期標(biāo)記同一天揚(yáng)花的穎花以保障所取穎花樣品的生育進(jìn)程一致;于水稻灌漿初期(花后第8天)進(jìn)行夜間高溫處理。高溫處理結(jié)束后剪取帶標(biāo)記的穎花提取水稻籽??偟鞍踪|(zhì),采用8-plex i TRAQ試劑盒進(jìn)行蛋白質(zhì)樣品的差異顯色標(biāo)記,標(biāo)記樣品采用高效液相系統(tǒng)進(jìn)行質(zhì)譜鑒定及定量分析。【結(jié)果】鑒定并定量了3130個(gè)蛋白質(zhì),耐熱與熱敏感水稻純系間存在36個(gè)差異表達(dá)蛋白質(zhì)。蛋白質(zhì)功能注釋結(jié)果顯示,鑒定的36個(gè)差異表達(dá)蛋白質(zhì)中僅14個(gè)蛋白質(zhì)(占38.9%)注釋了功能,12個(gè)蛋白質(zhì)(占33.3%)為推測(cè)性功能注釋,10個(gè)(占27.8%)為功能未知的蛋白質(zhì)。已注釋功能的14個(gè)蛋白質(zhì)中,5個(gè)蛋白質(zhì)參與能量代謝,3個(gè)蛋白質(zhì)參與物質(zhì)轉(zhuǎn)運(yùn)與代謝,2個(gè)蛋白質(zhì)參與光合作用,3個(gè)蛋白質(zhì)為響應(yīng)逆境的鋅指蛋白質(zhì),1個(gè)蛋白質(zhì)為響應(yīng)逆境的熱激蛋白質(zhì)。【結(jié)論】灌漿期夜間高溫影響水稻籽粒細(xì)胞內(nèi)參與能量代謝、物質(zhì)轉(zhuǎn)運(yùn)與代謝、光合作用等相關(guān)蛋白質(zhì)的表達(dá)模式。水稻籽粒細(xì)胞中鋅指蛋白質(zhì)Q67 TK9、Q10N88上調(diào)表達(dá),鋅指蛋白質(zhì)Q5YLY5下調(diào)表達(dá),有利于提高水稻灌漿期對(duì)夜間高溫的耐熱性。

    水稻;灌漿期;夜間高溫;蛋白質(zhì)組學(xué);iTRAQ技術(shù)

    全球氣候變暖,尤其是夏季頻繁高溫對(duì)水稻等作物產(chǎn)量和品質(zhì)的不利影響,已受到廣泛關(guān)注[1,2]。在全球氣溫升高的大趨勢(shì)下,夜間氣溫的升高幅度顯著大于白天。而且夜間平均氣溫每升高1℃可導(dǎo)致水稻減產(chǎn)10%,而白天氣溫的升高與水稻產(chǎn)量間無(wú)顯著相關(guān)性[3],表明夜間高溫是直接導(dǎo)致水稻減產(chǎn)的主要原因[4]。在我國(guó),夏季高溫?zé)岷χ饕l(fā)生在長(zhǎng)江中下游地區(qū),尤其在7月份,該區(qū)域受西太平洋副熱帶高壓控制,天氣晴朗、太陽(yáng)輻射強(qiáng)烈,極易出現(xiàn)高溫酷熱天氣[5,6]。該區(qū)域是我國(guó)雙季稻的主產(chǎn)區(qū),而且該區(qū)域雙季早稻的灌漿初期正值7月上、中旬。水稻灌漿初期(揚(yáng)花后第8~14天)是籽粒胚乳細(xì)胞形成和充實(shí)的關(guān)鍵時(shí)期,期間遇高溫?zé)岷?huì)加快胚乳細(xì)胞的分化而增加畸形胚乳細(xì)胞的比例,結(jié)果不利于同化物向胚乳細(xì)胞的運(yùn)輸和積累從而導(dǎo)致水稻籽粒的充實(shí)度下降、粒重降低[7]。同時(shí),高溫影響水稻籽粒細(xì)胞內(nèi)的酶活性,不利于籽粒儲(chǔ)藏物質(zhì)的轉(zhuǎn)化和積累進(jìn)而導(dǎo)致稻米品質(zhì)變劣[8,9]。長(zhǎng)江中下游地區(qū)的雙季早稻是其他作物難以替代的優(yōu)勢(shì)糧食作物,我國(guó)稻谷總產(chǎn)量的幾次大波動(dòng)都與雙季早稻的減產(chǎn)密切相關(guān)。因此,研究水稻灌漿期耐夜間高溫的分子機(jī)制,加快選育對(duì)夜溫升高鈍感的優(yōu)良水稻新品種以保證我國(guó)稻谷產(chǎn)量和品質(zhì)的穩(wěn)定,進(jìn)而對(duì)保障國(guó)家的口糧安全具有重要意義。水稻灌漿期的耐熱性受多基因調(diào)控,其分子調(diào)控機(jī)制相對(duì)復(fù)雜[10,11]。明確哪些基因及其編碼的蛋白質(zhì)參與調(diào)控水稻灌漿期耐熱性是解析其分子機(jī)制的先決內(nèi)容。已有研究表明,灌漿期高溫抑制淀粉支解酶Ⅱb(BEⅡb)活性、增加直鏈淀粉支鏈長(zhǎng)度而減少直鏈淀粉的長(zhǎng)鏈數(shù)量[12];灌漿期高溫提高水稻胚乳細(xì)胞中焦磷酸化酶(AGPase)、可溶性淀粉合成酶(SSS)和淀粉分支酶(SBE)等與淀粉合成有關(guān)的酶活性[13]。在基因表達(dá)方面,灌漿期高溫抑制顆粒淀粉合成酶Ⅰ(GBSSⅠ)、顆粒淀粉合成酶Ⅱb(BEⅡb)、ADP-葡萄糖焦磷酸化酶(包括AGPS2b, AGPS1和AGPL2)、ADP-葡萄糖轉(zhuǎn)運(yùn)體等與淀粉合成相關(guān)基因以及丙酮酸磷酸雙激酶(PPDK)基因的表達(dá)[14,15];灌漿期高溫誘導(dǎo)α-淀粉酶(amylase)基因(包括Amy1A,Amy1C,Amy3 A,Amy3D和Amy3E)表達(dá),加劇胚乳細(xì)胞內(nèi)淀粉的水解而增大籽粒堊白面積[16]。在蛋白質(zhì)表達(dá)方面,灌漿期高溫誘導(dǎo)18.1 k D和17.9 k D熱激蛋白質(zhì)(HSP)、醇溶谷蛋白質(zhì)(Prolamin)和谷蛋白質(zhì)(Glutelin)表達(dá),抑制顆粒淀粉合成酶(granule-boundstarch synthase)、類過(guò)敏原蛋白質(zhì)(allergen-like protein)和延長(zhǎng)因子1β(elongation factor 1β)的表達(dá)[14,17]。在調(diào)控途徑方面,灌漿期夜間高溫誘導(dǎo)線粒體跨膜電子傳遞相關(guān)基因上調(diào)表達(dá)[18],影響三羧酸循環(huán)(TCA)、氨基酸和多胺的生物合成途徑[19]。這些基因和蛋白質(zhì)的鑒定及其表達(dá)模式的明確為進(jìn)一步闡明水稻灌漿期耐熱的分子機(jī)制奠定了研究基礎(chǔ)。然而,水稻在遭遇灌漿期夜間高溫脅迫時(shí),到底哪些基因和蛋白質(zhì)的表達(dá)及其呈現(xiàn)何種表達(dá)模式,其有利于提高水稻對(duì)灌漿期夜間高溫的耐性,尚鮮見報(bào)道[18]。本研究以近等基因系的耐熱水稻純系XN0437T和熱敏感純系XN0437S為供試材料,利用i TRAQ(同位素標(biāo)記相對(duì)和絕對(duì)定量,isobaric tags for relative and absolute quantitation)技術(shù)和高效液相質(zhì)譜檢測(cè)技術(shù)鑒定耐熱和熱敏感水稻純系在應(yīng)答灌漿期夜間高溫過(guò)程中的差異表達(dá)蛋白質(zhì),明確哪些蛋白質(zhì)的上調(diào)表達(dá)或下調(diào)表達(dá)有利于提高水稻灌漿期對(duì)夜間高溫的耐熱性,為進(jìn)一步闡明水稻灌漿期耐夜間高溫的分子機(jī)制奠定研究基礎(chǔ)。

    1 材料與方法

    1.1 供試水稻材料

    供試水稻材料為江西農(nóng)業(yè)大學(xué)早稻灌漿期耐熱抗早衰育種技術(shù)研究課題組選育的一對(duì)近等基因系,即耐熱純系XN0437T和熱敏感純系XN0437S。此對(duì)近等基因系源于協(xié)青早B/N22//協(xié)青早B的回交重組自交系,在江西南昌作為雙季早稻在常規(guī)栽培條件下種植,耐熱純系XN0437T和熱敏感純系XN0437S的株高分別為81.3 cm和81.0 cm,穗長(zhǎng)分別為18.4 cm和18.1 cm,每穗總粒數(shù)分別為94.3粒和95.0粒,全生育期均為106 d。SSR分子標(biāo)記檢測(cè)結(jié)果顯示此對(duì)近等基因系的基因組多態(tài)性僅為1.8%。兩水稻純系于水稻開花后第8天開始進(jìn)行48 h(38.0±0.5)℃高溫處理后,在常溫條件(25.0±0.5)℃下恢復(fù)生長(zhǎng)至成熟并調(diào)查其籽粒充實(shí)度。結(jié)果顯示耐熱水稻純系XN0437T(原名703T)的籽粒充實(shí)度為92.4%,熱敏感純系XN0437S(原名704S)的籽粒充實(shí)度為43.6%,籽粒充實(shí)度差異達(dá)到極顯著水平[20]。

    1.2 供試水稻材料的培育、高溫處理及取樣

    供試水稻材料的培育、高溫處理及取樣參照Liao等的方法進(jìn)行[18]。水稻的培育采用桶栽的方法,水稻抽穗期、揚(yáng)花期標(biāo)記同一天抽穗且同一天開花的穎花以保證所取樣品的生長(zhǎng)發(fā)育進(jìn)程一致。水稻揚(yáng)花后第8天開始在人工氣候室內(nèi)進(jìn)行高溫處理,夜間溫度設(shè)為38.0℃±0.5℃(高溫處理)和25.0℃±0.5℃(對(duì)照),處理和對(duì)照的白天溫度均設(shè)為(26.0±0.5)℃。暗期時(shí)長(zhǎng)10 h、光期時(shí)長(zhǎng)14 h,連續(xù)處理48 h。設(shè)置3次生物學(xué)重復(fù)。高溫處理結(jié)束時(shí),取處理與對(duì)照的帶標(biāo)記穎花各約5 g,采用錫箔紙包裹,置于液氮中速凍,并于-80℃下保存?zhèn)溆谩?/p>

    1.3 水稻籽??偟鞍踪|(zhì)提取及質(zhì)量檢測(cè)

    取水稻籽粒樣品3 g,經(jīng)液氮研磨成粉末后,轉(zhuǎn)移到15 m L離心管中,加入蛋白質(zhì)提取緩沖液(8 mol/L尿素,0.1%SDS,1 mmol/L PMSF,1 mmol/L DTT)于自動(dòng)渦旋混合器上室溫振蕩提取3 h后,4℃、14 000 r/min下離心15 min,取上清液轉(zhuǎn)移至新的15 m L離心管中,加入6倍體積的預(yù)冷丙酮,-20℃下沉淀過(guò)夜。沉淀過(guò)夜樣品于4℃、14 000 r/min下離心15 min,去上清液,將沉淀于4℃真空干燥后獲得總蛋白質(zhì)干粉,-80℃下保存?zhèn)溆谩?/p>

    取蛋白質(zhì)干粉50 mg,經(jīng)裂解液(含8 mol/L尿素,2 mol/L硫脲,質(zhì)量分?jǐn)?shù)4%CHAPS,p H3~10體積分?jǐn)?shù)0.5%兩性電解質(zhì),50 mmol/L DTT和1.0 mmol/L PMSF)超聲溶解后,于4℃、14 000 r/ min下離心15 min后,上清液(總蛋白質(zhì)溶液)轉(zhuǎn)移至1.5 m L離心管。采用SDS-PAGE電泳和考馬斯亮藍(lán)染色法檢測(cè)總蛋白質(zhì)的完整性,SDS-PAGE電泳參數(shù)與方法參照Liao等[21]。以牛血清白蛋白質(zhì)為標(biāo)準(zhǔn)品,采用Bradford法測(cè)定蛋白質(zhì)樣品濃度。

    1.4 蛋白質(zhì)酶解及iTRAQ標(biāo)記

    蛋白質(zhì)樣品先經(jīng)還原和烷基化后,采用胰蛋白質(zhì)酶酶解。每個(gè)樣品取100μg蛋白質(zhì),經(jīng)10 mmol/L二硫蘇糖醇(DTT)于37℃溫浴1 h、55 mmol/L碘乙酰胺于室溫烷化1 h后,加入3.3μg胰蛋白質(zhì)酶于37℃下消化12 h,消化結(jié)束時(shí)加入體積分?jǐn)?shù)1%的蟻酸100μL終止酶解反應(yīng),真空抽干獲得蛋白質(zhì)肽段干粉。蛋白質(zhì)肽段干粉復(fù)溶于8 mol/L尿素(含0.1%SDS)和500 mmol/L三乙基碳酸氫銨的水溶液后,采用AB Sciex公司的8-plex i TRAQ試劑盒進(jìn)行差異顯色標(biāo)記。將試劑盒中8管標(biāo)記試劑(113、114、115、116、117、118、119和121)經(jīng)50μL異丙醇稀釋后分別與相應(yīng)的蛋白質(zhì)肽段干粉復(fù)溶樣品混合,室溫放置1 h,各管中加入100μL去離子水使標(biāo)記物失活,混合8組樣品,凍干備用。

    隨機(jī)取3次高溫處理重復(fù)中的2次重復(fù)進(jìn)行i TRAQ分析。耐熱水稻純系的2個(gè)對(duì)照重復(fù)分別命名為TC1和TC2,2個(gè)夜間高溫脅迫重復(fù)分別命名為TT1和TT2,并用iTRAQ試劑113、114、115、116依次進(jìn)行標(biāo)記;熱敏感水稻純系的2個(gè)對(duì)照重復(fù)分別命名為SC1和SC2,2個(gè)夜間高溫脅迫重復(fù)分別命名為ST1和ST2,并用iTRAQ試劑117、118、119、121依次進(jìn)行標(biāo)記。

    1.5 iTRAQ標(biāo)記樣品的強(qiáng)陽(yáng)離子交換柱分級(jí)及質(zhì)譜檢測(cè)

    iTRAQ標(biāo)記的8個(gè)樣品等量混合后采用C18強(qiáng)陽(yáng)離子柱進(jìn)行樣品分級(jí)。經(jīng)強(qiáng)陽(yáng)離子色譜柱平衡緩沖液A(含p H 2.55、5 mmol/L KH2PO4,體積比為20%的乙腈、H3PO4)平衡25 min后,取標(biāo)記樣品混合物300μL,用強(qiáng)陽(yáng)離子色譜柱平衡緩沖液B (含p H 2.75、5 mmol/L KH2PO4,體積比為20%的乙腈,600 mmol/L KCl,H3PO4)稀釋7倍,正磷酸調(diào)p H值至2.5,離心取上清液進(jìn)行梯度洗脫,每1 min收集1管樣品洗脫液,洗脫液流速為0.2 ml/ min。線性洗脫程序如下:0-3 min,緩沖液B的體積分?jǐn)?shù)由0%升至5%;3-21 min,緩沖液B的體積分?jǐn)?shù)由5%升至20%;21-27 min,緩沖液B的體積分?jǐn)?shù)由20%升至30%;27-33min,緩沖液B的體積分?jǐn)?shù)由30%升至100%。梯度洗脫結(jié)束后,根據(jù)色譜圖的峰形和時(shí)間合并連續(xù)的多肽含量較少的樣品,共形成20個(gè)組分用于肽段的質(zhì)譜檢測(cè)。

    質(zhì)譜檢測(cè)采用Thermo Scientific公司生產(chǎn)的Q Exactive質(zhì)譜儀,高效液相系統(tǒng)為戴安NCS3500系統(tǒng)。液相色譜A相為體積分?jǐn)?shù)為99.9%的dd H2O和0.1%的甲酸,B相為體積分?jǐn)?shù)為99.9%的乙腈、0.1%的甲酸。液相色譜洗脫液流速為300 n L/min,液相色譜線性梯度洗脫程序如下:0-5 min,B相體積分?jǐn)?shù)保持3%不變;5-6 min,B相體積分?jǐn)?shù)由3%升至5%;6-36 min,B相體積分?jǐn)?shù)由5%升至20%;36-46 min,B相體積分?jǐn)?shù)由20%升至30%;46-50 min,B相體積分?jǐn)?shù)由30%升至80%,并在之后5 min內(nèi)B相體積分?jǐn)?shù)保持不變;55 -55.5 min,B相體積分?jǐn)?shù)由80%降至3%,并在之后19.5 min內(nèi)B相體積百分?jǐn)?shù)保持3%不變。檢測(cè)質(zhì)譜全掃描范圍m/z=350~1800,全掃描中選擇離子強(qiáng)度前20位的母離子用標(biāo)準(zhǔn)碰撞能為30 eV的HCD模式碎裂后,進(jìn)行二級(jí)質(zhì)譜序列測(cè)定,以報(bào)告離子的比例進(jìn)行定量。

    1.6 蛋白質(zhì)定量及功能注釋

    將質(zhì)譜鑒定結(jié)果導(dǎo)入Proteome Discoverer軟件,篩選選擇肽段可信度高(Peptide Confidence=high)的肽段,并由軟件根據(jù)同一特異肽段(Unique Peptide)在各樣本中的報(bào)告離子峰的面積比值計(jì)算蛋白質(zhì)的表達(dá)量。利用Proteome Discoverer軟件搜索Uniport數(shù)據(jù)庫(kù)(http://ww w.uniprot.org/),注釋蛋白質(zhì)的功能。

    1.7 蛋白質(zhì)編碼基因的GO功能和KEGG通路的顯著富集分析

    質(zhì)譜數(shù)據(jù)經(jīng)Proteome Discoverer軟件搜索Uniport數(shù)據(jù)庫(kù)鑒定蛋白質(zhì)后,按照蛋白質(zhì)編碼基因ID、利用KOBAS軟件搜索在線GO(Gene ontology)數(shù)據(jù)庫(kù)(http://geneontology.org/),采用Qvalue方法進(jìn)行FDR校正,從而注釋蛋白質(zhì)編碼基因產(chǎn)物屬性和蛋白質(zhì)編碼基因參與的通路并統(tǒng)計(jì)顯著富集的通路。

    1.8 差異表達(dá)蛋白質(zhì)分析

    根據(jù)軟件Proteome Discoverer計(jì)數(shù)獲得的蛋白質(zhì)表達(dá)量,分2步進(jìn)行差異表達(dá)蛋白質(zhì)分析。首先,依據(jù)蛋白質(zhì)表達(dá)量分析耐熱純系與熱敏感純系各自的夜間高溫脅迫與對(duì)照的差異表達(dá)蛋白質(zhì)(組內(nèi)分析)。若蛋白質(zhì)在夜間高溫脅迫中的表達(dá)量與其在對(duì)照中的表達(dá)量的比值大于1.5倍,則定義為組內(nèi)上調(diào)1.5倍的表達(dá)蛋白質(zhì);上述比值小于0.667(1/1.5),則定義為組內(nèi)下調(diào)1.5倍的表達(dá)蛋白質(zhì)。然后,利用組內(nèi)分析結(jié)果比較耐熱和熱敏感純系之間差異表達(dá)蛋白質(zhì)(組間分析)。以下兩種情況均定義為耐熱純系與熱敏感純系間的差異表達(dá)蛋白質(zhì):1)相同表達(dá)模式的差異蛋白質(zhì):在2個(gè)水稻純系中均呈現(xiàn)相同的上調(diào)或下調(diào)表達(dá)模式,且耐熱純系與熱敏感純系之間的上調(diào)或下調(diào)表達(dá)量的比值大于1.5倍或小于0.667(1/1.5)倍的蛋白質(zhì);2)不同表達(dá)模式的差異蛋白質(zhì),在2個(gè)水稻純系中呈現(xiàn)不同的表達(dá)模式,即在耐熱純系中上調(diào)表達(dá)而在熱敏感純系中下調(diào)表達(dá)的蛋白質(zhì),或在耐熱純系中下調(diào)表達(dá)而在熱敏感純系中上調(diào)表達(dá)的蛋白質(zhì)。

    1.9 蛋白質(zhì)編碼基因的mRNA表達(dá)量驗(yàn)證

    為了驗(yàn)證iTRAQ數(shù)據(jù)的可靠性,在不同功能分類中選擇了7個(gè)有代表性的差異表達(dá)蛋白質(zhì),根據(jù)其編碼基因在Gen Bank的EST數(shù)據(jù)庫(kù)(w ww.ncbi.nlm.nih.gov/dbEST)中mRNA序列,利用NCBI在線設(shè)計(jì)軟件primer-BLAST設(shè)計(jì)引物,采用熒光定量PCR方法檢測(cè)差異表達(dá)蛋白質(zhì)編碼基因在水稻應(yīng)答灌漿期夜間高溫的表達(dá)量變化。

    熒光定量PCR的水稻籽粒樣品來(lái)自“1.2供試水稻材料的培育、高溫處理及取樣”的3次夜間高溫處理重復(fù),樣品的總RNA提取和RNA反轉(zhuǎn)錄合成cDNA第一鏈、qPCR和表達(dá)量計(jì)算參照文獻(xiàn)[18]和[22]的方法進(jìn)行。

    2 結(jié)果與分析

    2.1 蛋白質(zhì)的鑒定及定量

    利用i TRAQ標(biāo)記技術(shù)和液相色譜串聯(lián)質(zhì)譜檢測(cè)技術(shù),供試的8個(gè)樣品共檢測(cè)到肽段次數(shù)為83 254次。通過(guò)軟件Proteome Discoverer過(guò)濾,得到了16 437條高可信度(Peptide Confidence=high)的肽段,其中可用于蛋白質(zhì)定量的特異肽段為12 911條??尚烹亩魏吞禺愲亩谓?jīng)Uniport數(shù)據(jù)庫(kù)搜索,最終鑒定并定量了3130個(gè)蛋白質(zhì)。

    蛋白質(zhì)的表達(dá)量的差異分析結(jié)果顯示,耐熱水稻純系中上調(diào)1.5倍的蛋白質(zhì)有61個(gè),下調(diào)1.5倍的蛋白質(zhì)有20個(gè)(圖1)。熱敏感水稻純系中上調(diào)1.5倍的蛋白質(zhì)有38個(gè),下調(diào)1.5倍的蛋白質(zhì)有18個(gè)。其中,在耐熱純系和熱敏感純系中均上調(diào)表達(dá)的蛋白質(zhì)有17個(gè),均下調(diào)表達(dá)的蛋白質(zhì)有8個(gè)。

    2.2 GO功能和KEGG通路的顯著富集結(jié)果

    圖1 耐熱純系與熱敏感純系的應(yīng)答夜間高溫的蛋白質(zhì)數(shù)量Fig.1.Protein amount from the heat-tolerant line and heat-sensitive line responding to high night temperature.

    蛋白質(zhì)編碼基因的GO功能(A)和KEGG通路(B)的顯著富集結(jié)果如圖2所示。蛋白質(zhì)編碼基因的GO顯著富集結(jié)果顯示(圖2-A),質(zhì)譜鑒定蛋白質(zhì)的GO功能產(chǎn)物屬性顯著富集(-logP值≥1.5)在生物學(xué)過(guò)程(Biological process)的萜類化合物生物合成過(guò)程(GO ID=0016114:Terpenoid biosynthetic process)、類異戊二烯生物合成過(guò)程(GO ID=0008299:Isoprenoid biosynthetic process)、單羧酸生物合成過(guò)程(GO ID=0072330和0032787: Monocarboxylic acid biosynthetic process)、有機(jī)含羥基化合物生物合成過(guò)程(GO ID=1901617和1901615:Organic hydroxy compound biosynthetic process)、有機(jī)酸生物合成過(guò)程(GO ID=0016053: Organic acid biosynthetic process)、羧酸生物合成過(guò)程(GO ID=0046394:Carboxylic acid biosyntheticprocess)、單一有機(jī)體生物合成過(guò)程(GOID=0044711:Single-organism biosynthetic process)、脂質(zhì)生物合成過(guò)程(GO ID=0008610:Lipid biosynthetic process)和小分子生物合成過(guò)程(GO ID=0044283:Small molecule biosynthetic process);其次,GO功能產(chǎn)物屬性顯著富集(-logP值≥1.5)在細(xì)胞組分(Cellular component)的酶抑制劑活性(GO ID=0004857:Enzyme inhibitor activity)、肽鏈內(nèi)切酶抑制劑活性(GO ID=004866:Endopeptidase inhibitor activity)、肽酶調(diào)節(jié)因子活性(GO ID=0061134:Peptidase regulator activity)、肽酶抑制因子活性(GO ID=0030414:Peptidase inhibitor activity)、肽鏈內(nèi)切酶調(diào)節(jié)因子活性(GO ID=0061135:Endopeptidase regulator activity)、絲氨酸蛋白質(zhì)酶抑制劑活性(GO ID=0004867:Serinetype endopeptidase inhibitor activity)、分子功能調(diào)節(jié)因子(GO ID=0098772:Molecular function regulator)和酶調(diào)節(jié)因子活性(GO ID=0030234:Enzyme regulator activity)。蛋白質(zhì)編碼基因的KEGG通路顯著富集結(jié)果顯示(圖2-B),質(zhì)譜鑒定蛋白質(zhì)的KEGG通路顯著富集(-logP值≥4)在內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工通路上(KEGG ID=Osa04141: Protein processing in endoplasmic reticulum)。

    2.3 差異表達(dá)蛋白質(zhì)的功能注釋及其分類

    依據(jù)1.7的差異表達(dá)蛋白質(zhì)判定方法,耐熱水稻純系與熱敏感水稻純系之間存在36個(gè)差異表達(dá)蛋白質(zhì)。這些差異表達(dá)蛋白質(zhì)經(jīng)Uniport數(shù)據(jù)庫(kù)搜索,發(fā)現(xiàn)只有14個(gè)蛋白質(zhì)(占38.9%)有功能注釋、還有12個(gè)蛋白質(zhì)(占33.3%)具有推測(cè)性功能注釋,另有10個(gè)(占27.8%)為功能未知蛋白質(zhì)。14個(gè)有功能注釋的蛋白質(zhì)中(表1),有5個(gè)蛋白質(zhì)(Q9FU69、Q7Y179、Q0IPF8、Q2R1N0和Q84T00)參與能量代謝;有3個(gè)蛋白質(zhì)(Q8GTK2、Q7Y007和Q10BU2)參與物質(zhì)轉(zhuǎn)運(yùn)與代謝;有2個(gè)蛋白質(zhì)(Q2R237和Q658I1)參與光合作用;另外4個(gè)蛋白質(zhì)(Q67TK9、Q10N88、Q84Q72和Q5YLY5)為應(yīng)答逆境的18.1 k D熱激蛋白質(zhì)和鋅指蛋白質(zhì)。

    僅有推測(cè)性功能注釋及功能未知的22個(gè)差異表達(dá)蛋白質(zhì)中,在2水稻純系中均呈現(xiàn)上調(diào)表達(dá)1.2倍以上的蛋白質(zhì)有10個(gè)(Q6I587、Q7XI37、Q69MT6、Q338P8、Q7XI22、Q8 H8U1、Q94GQ6、Q0E3F2、Q6YUH8、Q6ES31),均呈現(xiàn)下調(diào)表達(dá)1.2倍以上的蛋白質(zhì)有2個(gè)(Q7XQ85、Q0E225);在耐熱水稻純系中上調(diào)表達(dá)1.2倍以上而在熱敏感水稻純系中下調(diào)表達(dá)1.2倍以上的蛋白質(zhì)有8個(gè)(Q6ZJ47、Q6ZCC9、Q6ZCP8、Q53M16、Q9LWS6、C7JA46、Q0E0B4、Q7XVN6),在耐熱水稻純系中下調(diào)表達(dá)1.2倍以上而在熱敏感水稻純系中上調(diào)表達(dá)1.2倍以上的蛋白質(zhì)有2個(gè)(Q0DTK3、Q5N9C8)。

    2.4 差異表達(dá)蛋白質(zhì)的mRNA表達(dá)模式驗(yàn)證

    有代表性的7個(gè)差異表達(dá)蛋白質(zhì)及其編碼基因的m RNA表達(dá)模式如圖3所示。通過(guò)iTRAQ差異標(biāo)記結(jié)合蛋白質(zhì)譜分析獲得的蛋白質(zhì)表達(dá)模式與熒光定量PCR中的m RNA表達(dá)模式的趨勢(shì)基本一致,表明通過(guò)i TRAQ差異標(biāo)記結(jié)合高效液相色譜技術(shù)鑒定并定量的蛋白質(zhì)組其結(jié)果可靠。

    3 討論

    圖2 差異表達(dá)蛋白質(zhì)編碼基因的GO功能(A)和KEGG通路(B)顯著富集結(jié)果Fig.2.Significant enrichment results of Gene Ontology(A)and Kyoto Encyclopedia of Genes(B).

    iTRAQ差異標(biāo)記具有在同一個(gè)體系中同時(shí)標(biāo)記多個(gè)(4或6或8個(gè))樣品的優(yōu)點(diǎn),不僅可以增加后續(xù)質(zhì)譜分析的樣品通量,而且可以減少蛋白質(zhì)色譜分析過(guò)程中的定量誤差[23]。本研究通過(guò)iTRAQ差異標(biāo)記和液相質(zhì)譜技術(shù)鑒定并定量了耐熱與熱敏感水稻近等基因系籽粒中的3130個(gè)蛋白質(zhì)。這些蛋白質(zhì)中,在耐熱與熱敏感水稻純系中相對(duì)表達(dá)量超過(guò)1.5倍的蛋白質(zhì)分別有81個(gè)和56個(gè);進(jìn)一步比較這些蛋白質(zhì)在2個(gè)水稻純系之間的表達(dá)模式差異,結(jié)果發(fā)現(xiàn)耐熱與熱敏感水稻純系之間存在36個(gè)差異表達(dá)蛋白質(zhì)。這些差異表達(dá)蛋白質(zhì)經(jīng)數(shù)據(jù)庫(kù)搜索,發(fā)現(xiàn)僅14個(gè)蛋白質(zhì)已注釋功能,其余蛋白質(zhì)只注釋了推測(cè)性功能或功能未知。這14個(gè)已注釋功能的蛋白質(zhì)中(表1),有5個(gè)蛋白質(zhì)參與能量代謝, 3個(gè)蛋白質(zhì)參與物質(zhì)轉(zhuǎn)運(yùn)與代謝,2個(gè)蛋白質(zhì)參與光合作用,4個(gè)為逆境應(yīng)答蛋白質(zhì)(18.1 k D熱激蛋白質(zhì)和鋅指蛋白質(zhì))。

    鋅指蛋白質(zhì)是一類轉(zhuǎn)錄調(diào)控因子蛋白質(zhì),通過(guò)蛋白質(zhì)的鋅、鐵中心與DNA結(jié)合從而調(diào)控基因的表達(dá)[24]。它們不僅在植物生長(zhǎng)發(fā)育中起著重要的基因表達(dá)調(diào)控功能[25,26],而且在植物應(yīng)答生物脅迫和非生物脅迫過(guò)程中扮演著非常重要的調(diào)控基因表達(dá)的角色[27-31]。Chen等[25]的研究報(bào)道鋅指蛋白質(zhì)OsGZF1在水稻種胚發(fā)育過(guò)程中具有調(diào)控籽粒中谷蛋白質(zhì)形成功能。Cao等分析了水稻的Cys2-His2(C2H2)鋅指蛋白質(zhì)家族,發(fā)現(xiàn)其中的20個(gè)鋅指蛋白質(zhì)參與調(diào)控水稻的生長(zhǎng)發(fā)育,28個(gè)參與調(diào)控水稻的授粉受精,22個(gè)參與調(diào)控水稻的抗病性,而且鋅指蛋白質(zhì)VRF1是水稻抗病的不可或缺的調(diào)控因子[32]。Huang等發(fā)現(xiàn)CHY鋅指蛋白質(zhì)具有通過(guò)調(diào)控水稻葉片氣孔的開閉從而提高水稻抗旱性和耐鹽性的功能[27,32]。本研究通過(guò)比較耐熱與熱敏感水稻近等基因系的蛋白質(zhì)組的表達(dá)模式,發(fā)現(xiàn)鋅指蛋白質(zhì)Q67TK9(Zinc knuckle containing proteinlike)和Q10N88(ARF GAP-like zinc finger-containing protein ZIGA2)在耐熱水稻純系中上調(diào)表達(dá)2.3倍以上,而在熱敏感水稻純系中無(wú)明顯表達(dá)變化(<1.5倍);鋅指蛋白質(zhì)Q5YLY5在耐熱與熱敏感水稻純系中均呈現(xiàn)下調(diào)表達(dá),但是該蛋白質(zhì)在耐熱水稻純系中下調(diào)表達(dá)了4.167倍,而在熱敏感水稻純系中僅下調(diào)表達(dá)1.56倍。表明鋅指蛋白質(zhì)Q67TK9、Q10N88的上調(diào)表達(dá),以及鋅指蛋白質(zhì)Q 5 YLY 5的下調(diào)表達(dá)有利于提高水稻灌漿期對(duì)夜間高溫的耐熱性。

    表1 耐熱純系與熱敏感純系的差異表達(dá)蛋白質(zhì)的表達(dá)模式及其功能注釋Table 1.Expression model and functional illustration of the differentially expressed proteins between the heat-tolerant and the heat-sensitive lines.

    圖3 差異表達(dá)蛋白質(zhì)及其編碼基因的表達(dá)模式Fig.3.Expression pattern of the differentially expressed protein and its coding gene.

    研究表明,籽粒ATP酶促進(jìn)外源營(yíng)養(yǎng)物質(zhì)經(jīng)質(zhì)外體和共質(zhì)體的協(xié)同作用向淀粉胚乳的運(yùn)輸,同時(shí)ATP酶還為胚乳細(xì)胞吸收并積累蛋白質(zhì),進(jìn)一步發(fā)育成蛋白體提供動(dòng)力[33,34]。本研究發(fā)現(xiàn),在耐熱和熱敏感水稻應(yīng)答夜間高溫過(guò)程中,籽粒中與ATP合成過(guò)程相關(guān)的ATP合成酶亞基β(Q0IPF8、LOC_Os12g10570)和ATP酶家族蛋白質(zhì)(Q84T00、LOC_Os03g58800)以及與物質(zhì)運(yùn)轉(zhuǎn)過(guò)程相關(guān)的羧肽酶(Q8GTK2、LOC_Os07g46350)、GDSL脂肪酶或酰基水解酶家族蛋白質(zhì)(Q7Y007、LOC_Os03g47940)和種子萌發(fā)蛋白質(zhì)3-7 (Q10BU2、LOC_Os03g58980)均呈現(xiàn)差異表達(dá)模式。表明夜間高溫影響了ATP酶活性,進(jìn)而阻礙了營(yíng)養(yǎng)物質(zhì)向淀粉胚乳的運(yùn)轉(zhuǎn)和積累。

    熱激蛋白質(zhì)普遍存在于動(dòng)植物細(xì)胞內(nèi),它們通過(guò)控制蛋白質(zhì)的折疊和積累從而調(diào)控蛋白質(zhì)表達(dá)的分子或分子伴侶[35]。據(jù)報(bào)道,干旱、高鹽和高溫等環(huán)境條件均可誘導(dǎo)熱激蛋白質(zhì)表達(dá)[36-40]。前期研究中,我們通過(guò)蛋白質(zhì)雙向電泳技術(shù)檢測(cè)了水稻應(yīng)答灌漿期(白天+夜間)高溫的蛋白質(zhì)組[17],發(fā)現(xiàn)18.1 kD熱激蛋白質(zhì)和17.9 kD熱激蛋白質(zhì)在耐熱與熱敏感水稻純系中均呈現(xiàn)上調(diào)表達(dá),而且在熱敏感水稻純系中的上調(diào)表達(dá)倍數(shù)高于其在耐熱水稻純系中的上調(diào)表達(dá)倍數(shù)。本研究利用i TRAQ差異標(biāo)記和液相質(zhì)譜技術(shù)比較耐熱與熱敏感水稻近等基因系應(yīng)答夜間高溫脅迫的蛋白質(zhì)組差異,也檢測(cè)到了18.1 kD熱激蛋白質(zhì)在夜間高溫脅迫與對(duì)照、耐熱與熱敏水稻純系中呈現(xiàn)差異表達(dá),而且該蛋白質(zhì)在熱敏感水稻純系中的上調(diào)表達(dá)倍數(shù)(2.975倍)高于其在耐熱純系中的上調(diào)表達(dá)倍數(shù)(1.906倍)。此外,本研究中還檢測(cè)到17.9 k D熱激蛋白質(zhì)(Q84Q77)在耐熱水稻純系及熱敏感水稻純系中均上調(diào)表達(dá),在耐熱水稻純系中上調(diào)了1.829倍、在熱敏感水稻純系中上調(diào)了2.292倍,但該蛋白質(zhì)在2個(gè)水稻純系間的表達(dá)倍數(shù)差異小于1.5倍。據(jù)此認(rèn)為18.1 kD熱激蛋白質(zhì)的表達(dá)豐度似可以作為水稻對(duì)灌漿期夜間高溫耐熱性的評(píng)判指標(biāo)。

    4 結(jié)論

    灌漿期夜間高溫影響水稻籽粒細(xì)胞內(nèi)參與能量代謝、物質(zhì)轉(zhuǎn)運(yùn)與代謝、光合作用等相關(guān)蛋白質(zhì)的表達(dá)模式。水稻籽粒細(xì)胞中鋅指蛋白質(zhì)Q67TK9、Q10N88的上調(diào)表達(dá),以及鋅指蛋白質(zhì)Q5YLY5的下調(diào)表達(dá),有利于提高水稻灌漿期對(duì)夜間高溫的耐熱性。

    [1] Wheeler T,von Braun J.Climate change impacts on global food security.Science,2013,341(6145):508-513.

    [2] Ray D K,Gerber J S,Macdonald G K,West P C.Climate variation explains a third of global crop yield variability.Nat Com mun,2015,6:5989.

    [3] Peng S,Huang J,Sheehy J E,Laza R C,Visperas R M, Zhong X,Centeno G S,Khush G S,Cassman K G.Rice yields decline with higher night temperature from global war-ming.Proc Natl Acad Sci USA,2004,101(27):9971-9975.

    [4] Shi W J,Yin X Y,Struik P C,Xie F M,Schmidt R C,Jagadish K S V.Grain yield and quality responses of tropical hybrid rice to high night-time temperature.Field Crop Res,2016, 190:18-25.

    [5] 謝曉金,李秉柏,王琳,戴秦如,申雙和.長(zhǎng)江中下游地區(qū)高溫時(shí)空分布及水稻花期的避害對(duì)策.中國(guó)農(nóng)業(yè)氣象,2010,31 (1):144-150.Xie X J,Li B B,Wang L,Dai Q R,Shen S H.Spatial and temporal distribution of high temperatureand strategies to rice florescence harm in the lower-middle reaches of Yangtze River.Chi J Agrometeorol,2010,31(1):144-150.(in Chinese with English abstract)

    [6] 張校瑋.我國(guó)極端氣候時(shí)空特征及風(fēng)險(xiǎn)分析.上海:上海師范大學(xué),2012.Zhang J W.Spatial and temporal characteristics of extreme climate and risks analysis in China.Shanghai:Shanghai Normal University,2012.(in Chinese with English abstract)

    [7] Morita S,Yonemaru J,Takanashi J.Grain growth and endosperm cell size under high night temperatures in rice(Oryza sativa L.).Ann Bot,2005,95(4):695-701.

    [8] Lin C J,Li C Y,Lin S K,Yang F H,Huang J J,Liu Y H, Lur H S.Influence of high temperature during grain filling on the accumulation of storage proteins and grain quality in rice (Oryza sativa L.).J Agric Food Chem,2010,58(19): 10545-10552.

    [9] Ahmed N,Tetlow I J,Nawaz S,Iqbal A,Mubin M,Nawaz U R M S,Butt A,Lightfoot D A,Maekawa M.Effect of high temperature on grain filling period,yield,amylose content and activity of starch biosynthesis enzymes in endosperm of basmati rice.J Sci Food Agric,2015,95(11):2237-2243.

    [10]Jagadish S V K,Murty M V R,Quick W P.Rice responses to rising temperatures-challenges,perspectives and future directions.Plant,Cell Environ,2015,38:1686-1698.

    [11]Liu J,Feng L,Li J,He Z.Genetic and epigenetic control of plant heat responses.Front Plant Sci,2015,6:267.

    [12]Jiang H,Dian W,Wu P.Effect of high temperature on fine structure of amylopectin in rice endosperm by reducing the activity of the starch branching enzyme.Phytochemistry,2003, 63(1):53-59.

    [13]張桂蓮,廖斌,武小金,肖應(yīng)輝,肖浪濤,陳立云.高溫對(duì)水稻胚乳淀粉合成關(guān)鍵酶活性及內(nèi)源激素含量的影響.植物生理學(xué)報(bào),2014,50(12):1840-1844.Zhang G L,Liao B,Wu X J,Xiao Y H,Xiao L T,Chen L Y.Effect of high temperature on activities of enzymes associated with starch synthesis and hormones contents in endosperm of rice.J Plant Physiol,2014,50(12):1840-1844.(in Chinese with English abstract)

    [14]Yamakawa H,Hirose T,Kuroda M,Yamaguchi T.Comprehensive expression profiling of rice grain filling-related genes under high temperature using DNA microarray.Plant Physiol,2007,144(1):258-277.

    [15]Mitsui T,Shiraya T,Kaneko K,Wada K.Proteomics of rice grain under high temperature stress.Front Plant Sci,2013, 4:36.

    [16]Hakata M,Kuroda M,Miyashita T,Yamaguchi T,Kojima M,Sakakibara H,Mitsui T,Yamakawa H.Suppression of alpha-amylase genes improves quality of rice grain ripened under high temperature.Plant Biotechnol J,2012,10(9): 1110-1117.

    [17]廖江林,宋宇,鐘平安,周會(huì)汶,張宏玉,黃英金.耐熱和熱敏感水稻應(yīng)答灌漿初期高溫脅迫過(guò)程中的差異表達(dá)蛋白質(zhì)鑒定.中國(guó)農(nóng)業(yè)科學(xué),2014,47(16):3121-3131.Liao J L,Song Y,Zhong P A,Zhou H W,Zhang H Y, Huang Y J.Identification of the differentially expressed proteins between heat-tolerant and heat-sensitive rice responding to high-temperature stress at the early milky stage.Sci Agric Sin,2014,47(16):3121-3131.(in Chinese with English abstract)

    [18]Liao J L,Zhou H W,Peng Q,Zhong P A,Zhang H Y,He C,Huang Y J.Transcriptome changes in rice(Oryza sativa L.)in response to high night temperature stress at the early milky stage.BMC Genom,2015,16:18.

    [19]Yamakawa H,Hakata M.Atlas of rice grain filling-related metabolism under high temperature:Joint analysis of metabolome and transcriptome demonstrated inhibition of starch accumulation and induction of amino acid accumulation.Plant Cell Physiol,2010,51(5):795-809.

    [20]Liao J L,Zhang H Y,Shao X L,Zhong P A,Huang Y J.I-dentification for heat tolerance in backcross recombinant lines and screening of backcross introgression lines with heat tolerance at milky stage in rice.Rice Sci,2011,18(4):279-286.

    [21]Liao J L,Huang Y J.Evaluation of protocols used in 2-D electrophoresis for proteome analysis of young rice caryopsis.GPB,2011,9(6):229-237.

    [22]Liao J L,Zhang H Y,Liu J B,Zhong P A,Huang Y J.Identification of candidate genes related to rice grain weight under high-temperature stress.Plant Sci,2012,196:32-43.

    [23]Frost D C,Greer T,Xiang F,Liang Z,Li L.Development and characterization of novel 8-plex Di Leu isobaric labels for quantitative proteomics and peptidomics.Rapid Com mun Mass Spectrom,2015,29(12):1115-1124.

    [24]Klug A.The discovery of zinc fingers and their applications in gene regulation and genome manipulation.Annu Rev Biochem,2010,79:213-231.

    [25]Chen Y,Sun A,Wang M,Zhu Z,Ouwerkerk P B.Functions of the CCCH type zinc finger protein OsGZF1 in regulation of the seed storage protein GluB-1 from rice.Plant Mol Biol, 2014,84(6):621-634.

    [26]Sekhar S,Gharat S A,Panda B B,Mohaptra T,Das K,Kariali E,Mohapatra P K,Shaw B P.Identification and characterization of differentially expressed genes in inferior and superiorspikelets of rice cultivars with contrasting panicle-compactness and grain-filling properties.PLOS One,2015,10(12): e145749.

    [27]Huang X Y,Chao D Y,Gao J P,Zhu M Z,Shi M,Lin H X.A previously unknown zinc finger protein,DST,regulates drought and salt tolerance in rice via stomatal aperture control.Gene Dev,2009,23(15):1805-1817.

    [28]Tyagi H,Jha S,Sharma M,Giri J,Tyagi A K.Rice saps are responsive to multiple biotic stresses and overexpression of OsSAP1,an A20/AN1 zinc-finger protein,enhances the basal resistance against pathogen infection in tobacco.Plant Sci, 2014,225:68-76.

    [29]Wang F,Tong W,Zhu H,Kong W,Peng R,Liu Q,Yao Q.A novel Cys2/His2 zinc finger protein gene from sweetpotato, IbZFP1,is involved in salt and drought tolerance in transgenic Arabidopsis.Planta,2016,243(3):783-797.

    [30]Yue X,Que Y,Xu L,Deng S,Peng Y,Talbot N J,Wang Z.ZNF1 encodes a putative C2 H2 zinc-finger protein essential for appressorium differentiation by the rice blast fungus Magnaporthe oryzae.Mol Plant Microbe Interact,2016,29(1):22-35.

    [31]Zhang Y,Lan H,Shao Q,Wang R,Chen H,Tang H,Zhang H,Huang J.An A20/AN1-type zinc finger protein modulates gibberellins and abscisic acid contents and increases sensitivity to abiotic stress in rice(Oryza sativa L.).J Exp Bot,2016, 67(1):315-326.

    [32]Cao H,Huang P,Zhang L,Shi Y,Sun D,Yan Y,Liu X, Dong B,Chen G,Snyder J H,Lin F,Lu J.Characterization of 47 Cys2-His2 zinc finger proteins required for the development and pathogenicity of the rice blast fungus Magnaporthe oryzae.New Phytol,2016,doi:10.1111/nph.13948.

    [33]Wei C X,Lan S Y,Xu Z X.Ultracytochemical localization and functional analysis of ATPase during the endosperm development in rice(Oryza sativa L.).Sci Agric Sin,2003,36 (3):259-262.

    [34]Zhou Z Q,Lan S Y,Xu Z X,Yang Z M.Dynamic change of ATPase activity on amyloplasts and protein bodies during the endosperm development in rice(Oryza sativa L.).Acta Biol Ex p Sin,2005,38(1):7-15.

    [35]Al-Whaibi M H.Plant heat-shock proteins:A mini review.J King Saud Univ-Sci,2011,23(2):139-150.

    [36]Gonzalez-Schain N,Dreni L,Lawas L M,Galbiati M,Colombo L,Heuer S,Jagadish K S,Kater M M.Genome-wide transcriptome analysis during anthesis reveals new insights into the molecular basis of heat stress responses in tolerant and sensitive rice varieties.Plant Cell Physiol,2016,57(1):57-68.

    [37]Das S,Krishnan P,Mishra V,Kumar R,Ramakrishnan B, Sinqh N K.Proteomic changes in rice leaves grown under open field high temperature stress conditions.Mol Biol Rep,2015, 42(11):1545-1558.

    [38]Xiang J H,Ran J,Zou J,Zhou X,Liu A,Zhang X,Peng Y, Tang N,Luo G,Chen X.Heat shock factor Os Hsf B2b negatively regulates drought and salt tolerance in rice.Plant Cell Rep,2013,32(11):1795-1806.

    [39]Liu A L,Zou J,Liu C F,Zhou X Y,Zhang X W,Luo G Y, Chen X B.Over-expression of Os Hsf A7 enhanced salt and drought tolerance in transgenic rice.BMB Rep,2013,46(1): 31-36.

    [40]Zou J,Liu C,Liu A,Zou D,Chen X.Overexpression of Os Hsp17.0 and Os Hsp23.7 enhances drought and salt tolerance in rice.J Plant Physiol,2012,169(6):628-635.

    Analysis on Comparative Proteomics of Rice Grain Between Heat-tolerant and Heat-sensitive Lines Under High Night Temperature Stress at Filling Stage

    HUANG Xiaoping1,ZHANG Hongyu2,LEI Gang1,WANG Zhimei1,ZHANG Zhi1,HE Chao1, LIAO Jianglin1,2,?,HUANG Yingjin1,2,?
    (1Key Laboratory of Crop Physiology,Ecology and Genetic Breeding,Ministry of Education,Jiangxi Agricultural University, Nanchang 330045;2 Key Laboratory of Agriculture Responding to Climate Change,Jiangxi Agricultural University,Nanchang 330045;?Corresponding author,E-mail:jlliao514815@163.com;yjhuang_cn@126.com)

    【Objective】To understand molecular mechanism of rice tolerance to high night temperature at filling stage, we identify the differentially expressed proteins(DEPs)and screen the proteins involving in regulating the heat-tolerance in rice.【Method】The heat-tolerant rice line XN0437T and the heat-sensitive rice line XN0437S originated from a near-genetic lines were used as plant materials.Rice was pot-cultivated under traditional management.To ensure that only grain samples with uniform growth were used for proteomes analysis,rice panicle bloomed on the same date were labeled with same markers and spikelets flowering on the same date on the labeled panicle were further labeled.On the 8th day after flowering,rice plants with same label were transferred to chambers for high night temperature treatment.Rice grains with same label were harvested and total proteins of the grain samples were extracted after high temperature treatment.The 8-plex iTRAQ kits combined with the LC-MS/MS technology were used to analyze the DEPs between the heat-tolerant and heat-sensitive lines.【Result】After protein’s databases searching and differentially expressed proteins analysis,3130 proteins were finally identified and 36 proteins showed differential expression levels between the heat-tolerant and heat-sensitive line.For the 36 DEPs,14 proteins(38.9%) had functional annotation,12(33.3%)had putative function,and 10(33.3%)were functional unknown proteins.The 14 function-known proteins were involved in energy metabolism(five proteins),transport and metabolism(three proteins),photosynthesis(two proteins)and defense response(four proteins).【Conclusion】High night temperature impacts the expressed patterns of the proteins involved in energy metabolism,matter transport and metabolism, photosynthesis and defense response in rice at filling stage.We suggested that up-regulated expression of the zinc finger proteins(Q67TK9 and Q10N88),and down-regulated expression of the zinc finger protein(Q5YLY5)in rice grain could enhance the heat-tolerance in rice.

    rice;filling stage;high night temperature;proteomics;iTRAQ technique

    Q948.112+.2;S511.01

    A

    1001-7216(2017)01-0013-10

    2016-06-24;修改稿收到日期:2016-09-17。

    江西省青年科學(xué)基金重點(diǎn)資助項(xiàng)目(20133 ACB21004);國(guó)家自然科學(xué)基金資助項(xiàng)目(31260315,31471467)。

    猜你喜歡
    鋅指籽粒灌漿
    鋅指蛋白與肝細(xì)胞癌的研究進(jìn)展
    籽粒莧的飼用價(jià)值和高產(chǎn)栽培技術(shù)
    籽粒莧的特性和種植技術(shù)
    谷子灌漿期噴施硫酸鋅增產(chǎn)
    C2H2型鋅指蛋白在腫瘤基因調(diào)控中的研究進(jìn)展
    玉米機(jī)械脫粒籽粒含水量與破碎率的相關(guān)研究
    無(wú)蓋重固結(jié)灌漿在烏弄龍水電站的應(yīng)用
    江西建材(2018年4期)2018-04-10 12:37:06
    自流可控灌漿在堰塞體防滲加固處理中的應(yīng)用
    商麥1619 籽粒灌漿的特性
    Myc結(jié)合的鋅指蛋白1評(píng)估實(shí)驗(yàn)性急性胰腺炎疾病嚴(yán)重程度的價(jià)值
    欧美成人a在线观看| 99视频精品全部免费 在线| 三级毛片av免费| 色在线成人网| 欧美高清成人免费视频www| 国产精华一区二区三区| 午夜福利欧美成人| 成人高潮视频无遮挡免费网站| 精品日产1卡2卡| 免费观看人在逋| avwww免费| 色5月婷婷丁香| 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 欧美最黄视频在线播放免费| 久久午夜亚洲精品久久| 美女高潮喷水抽搐中文字幕| 日本熟妇午夜| 美女cb高潮喷水在线观看| 国产成年人精品一区二区| 日韩免费av在线播放| 中文字幕久久专区| 亚洲欧美日韩高清专用| 精品久久久久久久久av| 中文资源天堂在线| 色播亚洲综合网| 在线天堂最新版资源| a级一级毛片免费在线观看| 亚洲无线在线观看| 日本五十路高清| 内射极品少妇av片p| 美女 人体艺术 gogo| 五月伊人婷婷丁香| 99国产精品一区二区三区| 日本在线视频免费播放| 亚洲欧美日韩高清专用| 九九热线精品视视频播放| 亚洲,欧美,日韩| 又爽又黄a免费视频| 欧美黑人巨大hd| 美女黄网站色视频| 熟女电影av网| 亚洲无线观看免费| 精品久久久久久成人av| 男人狂女人下面高潮的视频| av视频在线观看入口| 草草在线视频免费看| 国内少妇人妻偷人精品xxx网站| 美女 人体艺术 gogo| 亚洲av一区综合| 国产蜜桃级精品一区二区三区| 日本免费一区二区三区高清不卡| 久久久久久久亚洲中文字幕 | 免费无遮挡裸体视频| 免费在线观看成人毛片| 一个人观看的视频www高清免费观看| 波多野结衣高清作品| 国产黄片美女视频| 久久久精品大字幕| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 丰满人妻熟妇乱又伦精品不卡| 青草久久国产| 亚洲av电影在线进入| av欧美777| 九九在线视频观看精品| 黄色一级大片看看| xxxwww97欧美| 久久午夜福利片| 88av欧美| 亚洲自拍偷在线| 亚洲色图av天堂| 看黄色毛片网站| 国产又黄又爽又无遮挡在线| 99久久精品国产亚洲精品| 91av网一区二区| 成人国产综合亚洲| 极品教师在线视频| www.www免费av| 免费电影在线观看免费观看| 悠悠久久av| 在线国产一区二区在线| 精品日产1卡2卡| 免费观看的影片在线观看| 免费av毛片视频| 欧美又色又爽又黄视频| 国产乱人伦免费视频| 欧美一级a爱片免费观看看| 夜夜夜夜夜久久久久| av视频在线观看入口| 一二三四社区在线视频社区8| av国产免费在线观看| 欧美日本亚洲视频在线播放| 91久久精品电影网| 在线国产一区二区在线| 老女人水多毛片| 亚洲欧美激情综合另类| 国内毛片毛片毛片毛片毛片| 国产欧美日韩精品亚洲av| 美女免费视频网站| 日本黄大片高清| 一个人看视频在线观看www免费| 中文在线观看免费www的网站| 99在线视频只有这里精品首页| 亚洲av一区综合| 美女高潮喷水抽搐中文字幕| 久久国产乱子免费精品| 精品一区二区三区人妻视频| 国产白丝娇喘喷水9色精品| 欧美另类亚洲清纯唯美| 日韩有码中文字幕| 看片在线看免费视频| 色综合站精品国产| 亚洲av免费在线观看| 亚洲国产精品sss在线观看| av福利片在线观看| 18+在线观看网站| 国内精品久久久久久久电影| 亚洲av电影在线进入| 88av欧美| 国产主播在线观看一区二区| 精品人妻偷拍中文字幕| 激情在线观看视频在线高清| 麻豆一二三区av精品| 亚洲av电影不卡..在线观看| 亚洲五月婷婷丁香| 久久精品国产亚洲av香蕉五月| 国产精品亚洲av一区麻豆| 久久久成人免费电影| 国产一级毛片七仙女欲春2| 天堂网av新在线| 日本三级黄在线观看| 俺也久久电影网| 亚洲精品粉嫩美女一区| 丰满的人妻完整版| 搡女人真爽免费视频火全软件 | 精品久久久久久成人av| 精品熟女少妇八av免费久了| 十八禁人妻一区二区| 中文字幕熟女人妻在线| 精品一区二区免费观看| 精品人妻1区二区| 两个人的视频大全免费| 婷婷精品国产亚洲av| 白带黄色成豆腐渣| netflix在线观看网站| 国产在线男女| 又爽又黄a免费视频| 亚洲精品乱码久久久v下载方式| 在线看三级毛片| www.www免费av| 一个人免费在线观看电影| 亚洲性夜色夜夜综合| 色哟哟哟哟哟哟| 又紧又爽又黄一区二区| 搡老熟女国产l中国老女人| 国产精品电影一区二区三区| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久久久免 | 久久人人爽人人爽人人片va | 真人一进一出gif抽搐免费| 亚洲精华国产精华精| 搡老妇女老女人老熟妇| 国产精品女同一区二区软件 | 欧美日本视频| 免费人成在线观看视频色| 国产色爽女视频免费观看| 又爽又黄无遮挡网站| 日本撒尿小便嘘嘘汇集6| 男人的好看免费观看在线视频| 国产精品久久视频播放| 午夜亚洲福利在线播放| 中文字幕熟女人妻在线| 丰满乱子伦码专区| 国产精品三级大全| 亚洲av美国av| 免费高清视频大片| 黄色女人牲交| 嫩草影视91久久| 人妻制服诱惑在线中文字幕| 中文资源天堂在线| 午夜视频国产福利| 欧美日韩国产亚洲二区| 久久精品国产亚洲av天美| 亚洲成a人片在线一区二区| 99精品久久久久人妻精品| 亚洲人与动物交配视频| 最后的刺客免费高清国语| а√天堂www在线а√下载| 最近视频中文字幕2019在线8| 亚洲综合色惰| 99国产精品一区二区三区| 91在线观看av| av专区在线播放| 老女人水多毛片| 久久九九热精品免费| 俄罗斯特黄特色一大片| 给我免费播放毛片高清在线观看| 能在线免费观看的黄片| 久久伊人香网站| 国产av麻豆久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 动漫黄色视频在线观看| 欧美日本亚洲视频在线播放| 日韩欧美三级三区| 我要搜黄色片| 一个人观看的视频www高清免费观看| 久久午夜亚洲精品久久| or卡值多少钱| 欧美激情在线99| 色播亚洲综合网| 老鸭窝网址在线观看| 毛片一级片免费看久久久久 | 制服丝袜大香蕉在线| 欧美黑人欧美精品刺激| 久久这里只有精品中国| 久久久成人免费电影| 久久99热这里只有精品18| 国产蜜桃级精品一区二区三区| 亚洲国产精品sss在线观看| 18禁黄网站禁片免费观看直播| 亚洲成人久久性| 欧美日韩黄片免| 国产视频一区二区在线看| 亚洲三级黄色毛片| 欧美高清成人免费视频www| 人人妻人人看人人澡| 18+在线观看网站| 国产熟女xx| 欧美日韩瑟瑟在线播放| av福利片在线观看| 久久精品91蜜桃| 男人舔奶头视频| 人妻夜夜爽99麻豆av| 小说图片视频综合网站| 亚洲经典国产精华液单 | 亚洲在线观看片| 精品99又大又爽又粗少妇毛片 | 老熟妇乱子伦视频在线观看| 中文字幕人妻熟人妻熟丝袜美| xxxwww97欧美| 老鸭窝网址在线观看| 十八禁人妻一区二区| 岛国在线免费视频观看| 岛国在线免费视频观看| 色精品久久人妻99蜜桃| www.www免费av| 亚洲最大成人手机在线| 在现免费观看毛片| 国产欧美日韩一区二区精品| 午夜精品久久久久久毛片777| 成人av在线播放网站| 日韩有码中文字幕| 简卡轻食公司| av女优亚洲男人天堂| 一级作爱视频免费观看| 成年免费大片在线观看| 久久久国产成人精品二区| 特大巨黑吊av在线直播| 天堂影院成人在线观看| 久久久久久久亚洲中文字幕 | 有码 亚洲区| 亚洲精品亚洲一区二区| 久久人人精品亚洲av| 男女做爰动态图高潮gif福利片| 男女床上黄色一级片免费看| 国产av不卡久久| 偷拍熟女少妇极品色| 啪啪无遮挡十八禁网站| 少妇的逼好多水| 亚洲无线在线观看| 欧美性猛交╳xxx乱大交人| 亚洲性夜色夜夜综合| 深爱激情五月婷婷| 中国美女看黄片| 亚洲成av人片免费观看| 亚洲人与动物交配视频| 欧美中文日本在线观看视频| 如何舔出高潮| 亚洲av成人精品一区久久| 搡老熟女国产l中国老女人| 亚洲色图av天堂| 一级黄色大片毛片| 麻豆久久精品国产亚洲av| 日韩欧美国产一区二区入口| 精品久久久久久,| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 91av网一区二区| 男女床上黄色一级片免费看| 久久婷婷人人爽人人干人人爱| a在线观看视频网站| 亚洲精品影视一区二区三区av| 午夜精品久久久久久毛片777| 成熟少妇高潮喷水视频| 欧美另类亚洲清纯唯美| 又黄又爽又刺激的免费视频.| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美性猛交╳xxx乱大交人| 久久久久九九精品影院| 在线看三级毛片| 欧美成人一区二区免费高清观看| 国产真实乱freesex| 亚洲人与动物交配视频| 内地一区二区视频在线| av在线老鸭窝| 观看免费一级毛片| 久久久久久久久中文| 亚洲成人久久性| 国产高潮美女av| 一a级毛片在线观看| 国产黄a三级三级三级人| 免费看光身美女| 中文字幕高清在线视频| 欧美又色又爽又黄视频| 九九热线精品视视频播放| 美女 人体艺术 gogo| 在线国产一区二区在线| 婷婷精品国产亚洲av在线| 听说在线观看完整版免费高清| 97超级碰碰碰精品色视频在线观看| 免费av观看视频| 日本五十路高清| 免费看美女性在线毛片视频| 性欧美人与动物交配| 亚洲,欧美,日韩| 欧美又色又爽又黄视频| 亚洲av日韩精品久久久久久密| 亚洲人成伊人成综合网2020| 少妇的逼水好多| 自拍偷自拍亚洲精品老妇| a级毛片a级免费在线| 99热6这里只有精品| 午夜视频国产福利| 欧美一区二区精品小视频在线| 久久人人爽人人爽人人片va | 亚洲人成伊人成综合网2020| 国内精品一区二区在线观看| 18禁黄网站禁片免费观看直播| 嫩草影院入口| 亚洲人成网站高清观看| 亚洲 欧美 日韩 在线 免费| 久久久成人免费电影| 欧美成人a在线观看| 老司机福利观看| 欧美+日韩+精品| 午夜福利18| 欧美激情在线99| 国产又黄又爽又无遮挡在线| 日本成人三级电影网站| 国产精品一及| 一本久久中文字幕| 欧美最黄视频在线播放免费| 国产男靠女视频免费网站| 最近中文字幕高清免费大全6 | 亚洲无线在线观看| 久久久久九九精品影院| 级片在线观看| 日韩中文字幕欧美一区二区| 国产高清视频在线观看网站| 成人av在线播放网站| 好看av亚洲va欧美ⅴa在| av国产免费在线观看| 国产成年人精品一区二区| 国产真实乱freesex| 韩国av一区二区三区四区| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇熟女久久| 亚洲人与动物交配视频| a级一级毛片免费在线观看| 精品人妻1区二区| 99国产综合亚洲精品| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 非洲黑人性xxxx精品又粗又长| 久久久久免费精品人妻一区二区| 国产野战对白在线观看| 69av精品久久久久久| 亚洲在线自拍视频| 亚洲av成人不卡在线观看播放网| 日本免费a在线| 极品教师在线视频| 网址你懂的国产日韩在线| 88av欧美| 三级毛片av免费| 特级一级黄色大片| 久久久久久久久大av| av女优亚洲男人天堂| 特级一级黄色大片| 亚洲av.av天堂| 三级毛片av免费| 激情在线观看视频在线高清| 欧美午夜高清在线| 99精品久久久久人妻精品| 国产精品99久久久久久久久| av欧美777| 亚洲最大成人手机在线| 国产成+人综合+亚洲专区| 久久精品久久久久久噜噜老黄 | 亚洲av二区三区四区| 床上黄色一级片| 欧美3d第一页| 又紧又爽又黄一区二区| 亚洲国产欧美人成| 国产精品亚洲av一区麻豆| 亚洲中文字幕日韩| 成人无遮挡网站| www.熟女人妻精品国产| 久久热精品热| 两个人视频免费观看高清| 国产高清三级在线| 国产野战对白在线观看| 最近最新免费中文字幕在线| 久久久久久久亚洲中文字幕 | 欧美成人一区二区免费高清观看| 国产精品一区二区免费欧美| 久久精品人妻少妇| 成人午夜高清在线视频| 国产精品永久免费网站| 亚洲美女黄片视频| 岛国在线免费视频观看| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 欧美成人a在线观看| 日日夜夜操网爽| 国产精品久久久久久亚洲av鲁大| 亚洲五月天丁香| 成人美女网站在线观看视频| 亚洲成人精品中文字幕电影| 日韩欧美精品免费久久 | 88av欧美| 极品教师在线视频| 人妻夜夜爽99麻豆av| 亚洲av熟女| 高潮久久久久久久久久久不卡| 午夜福利18| 亚洲欧美清纯卡通| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av免费在线观看| 两人在一起打扑克的视频| 午夜福利在线观看免费完整高清在 | 乱码一卡2卡4卡精品| 99国产极品粉嫩在线观看| 中国美女看黄片| 18+在线观看网站| 在线免费观看的www视频| 首页视频小说图片口味搜索| 亚洲第一区二区三区不卡| 超碰av人人做人人爽久久| 成人一区二区视频在线观看| 国内揄拍国产精品人妻在线| 久久午夜亚洲精品久久| 欧美绝顶高潮抽搐喷水| 国产欧美日韩精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 不卡一级毛片| 一区二区三区四区激情视频 | 成人亚洲精品av一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久久久人妻精品电影| 丁香六月欧美| 亚洲成人中文字幕在线播放| 久久天躁狠狠躁夜夜2o2o| 免费观看精品视频网站| 日本五十路高清| 国模一区二区三区四区视频| 九九热线精品视视频播放| 99在线视频只有这里精品首页| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 亚洲第一欧美日韩一区二区三区| 日韩中文字幕欧美一区二区| 最近中文字幕高清免费大全6 | 亚洲最大成人av| 久久草成人影院| 欧美日韩黄片免| 天美传媒精品一区二区| 最近最新中文字幕大全电影3| 精品人妻偷拍中文字幕| 免费观看人在逋| 久久久精品大字幕| 成人三级黄色视频| 久久九九热精品免费| 91狼人影院| 亚洲熟妇中文字幕五十中出| 中出人妻视频一区二区| 此物有八面人人有两片| 最后的刺客免费高清国语| a级一级毛片免费在线观看| 日本 av在线| 精品无人区乱码1区二区| 免费av不卡在线播放| 嫩草影院精品99| 蜜桃久久精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 性插视频无遮挡在线免费观看| 日韩精品中文字幕看吧| 90打野战视频偷拍视频| 天堂影院成人在线观看| 午夜福利18| 国产淫片久久久久久久久 | 嫩草影视91久久| 18禁裸乳无遮挡免费网站照片| 日日摸夜夜添夜夜添av毛片 | 欧美三级亚洲精品| 欧美成人性av电影在线观看| 又黄又爽又刺激的免费视频.| 无遮挡黄片免费观看| 国产欧美日韩一区二区精品| 99久久九九国产精品国产免费| 我要看日韩黄色一级片| 日韩欧美国产在线观看| 国产男靠女视频免费网站| 在线观看午夜福利视频| 免费av观看视频| 女生性感内裤真人,穿戴方法视频| 国内久久婷婷六月综合欲色啪| av中文乱码字幕在线| 国产淫片久久久久久久久 | 亚洲精品一区av在线观看| 精品欧美国产一区二区三| 亚洲av五月六月丁香网| 精品免费久久久久久久清纯| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 99在线视频只有这里精品首页| 国产午夜精品论理片| 国产一区二区在线av高清观看| 校园春色视频在线观看| 少妇丰满av| 成年人黄色毛片网站| 亚洲av二区三区四区| 亚洲最大成人av| 亚洲人成网站高清观看| 久久精品国产清高在天天线| 久久久国产成人精品二区| 亚洲精品在线观看二区| 国产人妻一区二区三区在| 国产精品一区二区性色av| 最后的刺客免费高清国语| 国产乱人伦免费视频| 日本免费a在线| 国内精品久久久久精免费| av中文乱码字幕在线| 国产欧美日韩一区二区精品| 国产精品99久久久久久久久| 麻豆国产97在线/欧美| 一本精品99久久精品77| 一区福利在线观看| 蜜桃亚洲精品一区二区三区| 人人妻人人澡欧美一区二区| 国产欧美日韩一区二区精品| 女同久久另类99精品国产91| 久久久久国内视频| 婷婷精品国产亚洲av| 日韩av在线大香蕉| 国产伦在线观看视频一区| 免费看美女性在线毛片视频| 真人做人爱边吃奶动态| 淫秽高清视频在线观看| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 国产色爽女视频免费观看| 我的女老师完整版在线观看| 中文字幕人妻熟人妻熟丝袜美| 最新中文字幕久久久久| xxxwww97欧美| 亚洲美女黄片视频| 男人狂女人下面高潮的视频| 国产av在哪里看| 欧美日韩福利视频一区二区| 成年女人看的毛片在线观看| 久久6这里有精品| 亚洲最大成人中文| 我要搜黄色片| 看黄色毛片网站| 国产亚洲欧美在线一区二区| 精品国内亚洲2022精品成人| 亚洲18禁久久av| 淫妇啪啪啪对白视频| 精品午夜福利在线看| 成人一区二区视频在线观看| 国产成年人精品一区二区| 九色国产91popny在线| 五月玫瑰六月丁香| av专区在线播放| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 久久久久久国产a免费观看| 在线观看舔阴道视频| 国产精品自产拍在线观看55亚洲| 在线观看舔阴道视频| 日本免费a在线| 少妇熟女aⅴ在线视频| 一区二区三区免费毛片| 国产精品1区2区在线观看.| 男女视频在线观看网站免费| 一个人看的www免费观看视频| 欧美成人一区二区免费高清观看| 在线观看免费视频日本深夜| 美女cb高潮喷水在线观看| 人妻夜夜爽99麻豆av| 欧美日韩福利视频一区二区| 成人国产一区最新在线观看| aaaaa片日本免费| 婷婷色综合大香蕉| 久久精品国产亚洲av天美| 一进一出好大好爽视频| 成人美女网站在线观看视频| 九色成人免费人妻av| 在线观看舔阴道视频| а√天堂www在线а√下载| 亚洲色图av天堂| 97热精品久久久久久|