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

    水稻穗頂端退化突變體paa21的表型分析及基因克隆

    2022-02-02 03:34:14赫磊路凱趙春芳姚姝周麗慧趙凌陳濤朱鎮(zhèn)趙慶勇梁文化王才林朱麗張亞東
    中國農(nóng)業(yè)科學(xué) 2022年24期
    關(guān)鍵詞:小穗枝梗幼穗

    赫磊,路凱,趙春芳,姚姝,周麗慧,趙凌,陳濤,朱鎮(zhèn),趙慶勇,梁文化,王才林,朱麗,張亞東

    水稻穗頂端退化突變體的表型分析及基因克隆

    1江蘇省農(nóng)業(yè)科學(xué)院糧食作物研究所/國家耐鹽堿水稻技術(shù)創(chuàng)新中心華東中心/江蘇省優(yōu)質(zhì)水稻工程技術(shù)研究中心/國家水稻改良中心南京分中心/江蘇省農(nóng)業(yè)生物學(xué)重點(diǎn)實(shí)驗(yàn)室,南京 210014;2中國水稻研究所/水稻生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,杭州 310006

    【】水稻穗頂端退化嚴(yán)重影響產(chǎn)量,鑒定與克隆水稻穗頂端退化相關(guān)基因,可以豐富水稻穗發(fā)育調(diào)控的分子機(jī)理,為水稻高產(chǎn)分子設(shè)計(jì)育種提供理論基礎(chǔ)和基因資源?!尽繌木酒贩N武運(yùn)粳30號(hào)EMS突變體庫篩選到一份穩(wěn)定遺傳的穗頂端退化突變體()。對(duì)退化一次枝梗比例、每穗退化粒數(shù)占比、每穗粒數(shù)、株高、穗長、單株產(chǎn)量等農(nóng)藝性狀進(jìn)行統(tǒng)計(jì)。使用臺(tái)盼藍(lán)和伊文思藍(lán)染色檢測(cè)頂端小穗是否發(fā)生程序性細(xì)胞死亡。測(cè)定WT和不同發(fā)育時(shí)期幼穗和不同穗部位的H2O2含量。分別與秈稻II-32B、9311正反交進(jìn)行遺傳分析。利用與秈稻II-32B雜交構(gòu)建的F2群體進(jìn)行基因定位和克隆。使用SWISS-MODEL網(wǎng)站預(yù)測(cè)野生型和突變體蛋白的三維結(jié)構(gòu)。利用RT-qPCR分析ROS響應(yīng)標(biāo)志基因、程序性細(xì)胞死亡相關(guān)基因、過氧化氫酶相關(guān)基因的表達(dá)量。突變體發(fā)生嚴(yán)重的穗頂端退化,統(tǒng)計(jì)所有一次枝梗退化情況,發(fā)現(xiàn)退化小穗主要位于頂端的一次枝梗上。與WT相比,的株高、每穗粒數(shù)、穗長和單株產(chǎn)量均降低。通過觀察不同發(fā)育時(shí)期的幼穗,發(fā)現(xiàn)在突變體幼穗發(fā)育至12 cm時(shí),可見穗頂端退化表型。臺(tái)盼藍(lán)和伊文思藍(lán)染色結(jié)果表明突變體頂端小穗發(fā)生程序性細(xì)胞死亡。在退化的頂端小穗中觀察到更強(qiáng)烈的DAB染色;H2O2含量測(cè)定結(jié)果表明,與WT相比,穗中積累更高水平的ROS。遺傳分析表明突變表型受一對(duì)隱性核基因控制。圖位克隆結(jié)果發(fā)現(xiàn)中第二外顯子發(fā)生一個(gè)C到T的突變,導(dǎo)致丙氨酸突變?yōu)槔i氨酸。該基因編碼一個(gè)鋁激活的蘋果酸轉(zhuǎn)運(yùn)蛋白ALMT7。突變位點(diǎn)位于第4個(gè)跨膜螺旋上。SWISS-MODEL預(yù)測(cè)結(jié)果表明,該突變位點(diǎn)并未對(duì)突變體蛋白三維結(jié)構(gòu)造成明顯影響。RT-qPCR結(jié)果表明,在幼穗發(fā)育至10 cm時(shí),中ROS響應(yīng)標(biāo)志基因、和,程序性細(xì)胞死亡相關(guān)基因和,過氧化氫酶編碼基因、、的表達(dá)量較WT大幅升高。此外,10 cm幼穗中過氧化氫酶(CAT)的活性較WT明顯下降。幼穗在發(fā)育后期頂端小穗中積累過量的ROS,產(chǎn)生程序性細(xì)胞死亡,最終導(dǎo)致頂端小穗發(fā)生退化。

    水稻;穗頂端退化;活性氧;程序性細(xì)胞死亡

    0 引言

    【研究意義】水稻是中國重要的糧食作物,在保障國家糧食安全方面具有重要作用。水稻產(chǎn)量由千粒重、每穗實(shí)粒數(shù)和單位面積有效穗數(shù)決定[1]。水稻穗部發(fā)育關(guān)系到產(chǎn)量,水稻穗的形成涉及復(fù)雜的生理生化過程,包括分生組織發(fā)育、花序建立和籽粒發(fā)育等。深入研究水稻穗發(fā)育相關(guān)基因?qū)μ岣咚井a(chǎn)量具有重要意義?!厩叭搜芯窟M(jìn)展】目前,已經(jīng)報(bào)道了許多與穗發(fā)育相關(guān)的重要基因,涉及幼穗發(fā)育的不同階段,如腋生分生組織發(fā)育、分生組織之間的轉(zhuǎn)化、分枝伸長和末端小穗等的調(diào)控[2]。如,突變體()的枝梗數(shù)和每穗粒數(shù)的數(shù)量嚴(yán)重減少;相反,過量表達(dá)增加了枝梗和小穗的數(shù)目[3-5]。突變體()穗子顯著變大,分枝數(shù)和每穗小穗數(shù)也顯著增加[6]。()突變后分生組織活性增強(qiáng),穗軸與分枝長度縮短,一次、二次枝梗數(shù)目增加[7]。()編碼一個(gè)ERF/AP2轉(zhuǎn)錄因子,在禾本科植物中高度保守。強(qiáng)突變體的穗子嚴(yán)重卷曲、末端小穗被更高階的輪狀分枝取代[8]。在水稻、玉米、谷子等作物的產(chǎn)量形成過程中,穗頂端退化對(duì)最終的產(chǎn)量形成有重要影響[9-11]。頂端小穗退化主要由遺傳因素調(diào)控,同時(shí)也受溫度和肥水等非生物脅迫的影響,是水稻生產(chǎn)中的重大問題之一[12-14]。目前,已經(jīng)報(bào)道了一些與穗頂端退化相關(guān)的QTL,如—、、、等[14-17]。此外,利用穗頂端退化突變體已克隆數(shù)個(gè)涉及穗頂端退化的基因。()是第一個(gè)被克隆的水稻穗頂端退化的基因,編碼一個(gè)定位于細(xì)胞質(zhì)的包含SHD和VCA結(jié)構(gòu)域的SCAR/WAVE蛋白。與野生型相比,表現(xiàn)出株高變矮、葉尖枯萎、穗頂端嚴(yán)重退化、花藥花粉發(fā)育異常等多重表型[11]。Heng等[18]利用一個(gè)穗頂端退化突變體()克隆到一個(gè)控制穗頂端退化的顯性基因。該基因編碼一個(gè)具有轉(zhuǎn)運(yùn)蘋果酸功能的質(zhì)膜蛋白。細(xì)胞學(xué)分析發(fā)現(xiàn),穗頂端退化是由發(fā)生在幼穗發(fā)育后期的程序性細(xì)胞死亡造成的。()編碼一個(gè)定位于線粒體的含有胱硫醚β-合酶結(jié)構(gòu)域和DUF21結(jié)構(gòu)域的蛋白,突變體表現(xiàn)為穗頂端小穗退化和穗中間小穗的育性降低。研究表明,DPS1能夠與線粒體硫氧還蛋白Trx1和Trx20互作,并參與ROS清除。突變以后,突變體中脂肪酸代謝和ROS穩(wěn)態(tài)相關(guān)的生物過程受到影響,從而導(dǎo)致穗頂端退化[19]。Peng等[20]利用一個(gè)退化率高達(dá)60%的穗頂端退化突變體()克隆到另外一個(gè)控制水稻穗頂端退化的基因。編碼一個(gè)定位于細(xì)胞質(zhì)的類鈣調(diào)神經(jīng)磷酸酶B亞基互作蛋白激酶。的突變導(dǎo)致ROS在幼穗中積累,最終導(dǎo)致水稻穗細(xì)胞死亡產(chǎn)生頂端退化表型。編碼一個(gè)鈣氫離子交換蛋白,研究發(fā)現(xiàn)OsCAX1a參與Ca、Mn離子的轉(zhuǎn)運(yùn),突變后導(dǎo)致突變體中Ca、Mn等元素含量下降,產(chǎn)生葉尖枯和穗部頂端退化表型[21]?!颈狙芯壳腥朦c(diǎn)】目前,已克隆的水稻穗頂端退化相關(guān)基因涉及多種調(diào)控途徑,但仍不能完全解釋水稻穗頂端退化相關(guān)的分子機(jī)理。因此,穗頂端退化的分子機(jī)理仍需深入研究?!緮M解決的關(guān)鍵問題】本研究以粳稻品種武運(yùn)粳30號(hào)EMS誘變得到的穗頂端退化突變體為研究對(duì)象,通過對(duì)其進(jìn)行表型考察、組織化學(xué)分析、遺傳分析、基因定位以及基因表達(dá)分析,以期解析調(diào)控水稻穗頂端退化的分子機(jī)制,為提高水稻產(chǎn)量奠定基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料

    利用EMS(甲基磺酸乙酯)處理粳稻品種武運(yùn)粳30號(hào),田間篩選得到一份穗頂端退化的材料,經(jīng)多年種植穗頂端退化表型能穩(wěn)定遺傳,命名為(,)。將突變體分別與秈稻9311、II-32B雜交構(gòu)建F2分離群體,用于遺傳分析與基因定位。

    在水稻抽穗后,WT和分別取10株主穗,統(tǒng)計(jì)每穗退化粒數(shù)等表型。于成熟期統(tǒng)計(jì)株高、每穗粒數(shù)、穗長等表型。利用Excel和GraphPad Prism 8軟件統(tǒng)計(jì)、分析數(shù)據(jù)。

    1.2 臺(tái)盼藍(lán)、伊文思藍(lán)和DAB染色

    將WT和的小穗放入0.4%臺(tái)盼藍(lán)染液中,100℃染色10 min,然后,用2.5 mg·mL-1水合氯醛溶液脫色。將WT和的小穗放入0.25%伊文思藍(lán)染液中,室溫染色過夜。之后用無水乙醇脫色,掃描拍照保存。將WT和突變體小穗置于DAB染液(1 mg·mL-1)中,室溫染色8 h,染色完成后使用無水乙醇脫色,掃描拍照保存。

    1.3 H2O2含量、CAT酶活性測(cè)定

    取不同發(fā)育時(shí)期的WT和的小穗0.1—0.2 g,按1﹕1 000比例分別加入H2O2、CAT提取液冰浴勻漿(每個(gè)樣品3次重復(fù))。8 500 r/min離心10 min,取上清放置冰上待測(cè)。分別使用蘇州科銘公司的H2O2含量測(cè)定試劑盒和CAT酶活性測(cè)定試劑盒測(cè)定H2O2含量、CAT酶活性。

    1.4 遺傳分析及基因定位

    將與秈稻II-32B、9311進(jìn)行正反交,獲得F1株系和F2群體,分別統(tǒng)計(jì)F1、F2中穗發(fā)育正常和穗頂端退化表型分離比進(jìn)行遺傳分析。采用CTAB法提取葉片DNA。利用水稻生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室保存的均勻分布于12條染色體上的標(biāo)記引物分析親本多態(tài)性。在與秈稻II-32B構(gòu)建的F2群體中選取與表型一致的植株,利用多態(tài)性好的標(biāo)記對(duì)目的基因進(jìn)行連鎖分析。將目的基因初步定位于染色體某個(gè)區(qū)間。

    利用RAP-DB(https://rapdb.dna.affrc.go.jp/)網(wǎng)站分析定位區(qū)間內(nèi)的開放閱讀框(open reading frame,ORF),在Primer-Blast網(wǎng)站(https://www.ncbi.nlm.nih. gov/tools/primer-blast/)設(shè)計(jì)引物對(duì)WT和的候選基因測(cè)序。利用seqman軟件分析測(cè)序結(jié)果。

    1.5 蛋白三維結(jié)構(gòu)分析

    將WT和的氨基酸序列在SWISS-MODEL網(wǎng)站(https://swissmodel.expasy.org/)進(jìn)行同源建模,分析WT和的蛋白結(jié)構(gòu)。

    1.6 實(shí)時(shí)熒光定量PCR

    水稻組織提取RNA后,使用ReverTra Ace qPCR PCR試劑盒反轉(zhuǎn)獲得cDNA備用。使用AceQ qPCR SYBR Green Master Mix和熒光定量PCR儀(Applied Biosystems 7900HT)進(jìn)行RT-qPCR檢測(cè),以作為內(nèi)參基因,利用2–ΔΔCt計(jì)算基因的表達(dá)量。

    2 結(jié)果

    2.1 paa21表型分析

    2.1.1表現(xiàn)為嚴(yán)重的穗部頂端退化(,)多年種植,穗部頂端退化性狀穩(wěn)定,抽穗后,的頂端小穗發(fā)白、不能夠正常發(fā)育、開花,在成熟期后,頂端小穗風(fēng)干、脫落(圖1-A—C)。統(tǒng)計(jì)所有一次枝梗退化情況,發(fā)現(xiàn)絕大多數(shù)退化小穗位于頂端一次枝梗,約35%的一次枝梗發(fā)生退化,約20%的小穗發(fā)生退化(圖1-D—G)。此外,的株高低于WT,每穗粒數(shù)、穗長、單株產(chǎn)量均下降(圖1-H—K)。

    2.1.2 頂端小穗的退化發(fā)生在穗發(fā)育的后期 在種植過程中,發(fā)現(xiàn)其頂端小穗是在幼穗發(fā)育后期產(chǎn)生的退化現(xiàn)象。為了明確頂端小穗發(fā)生退化的時(shí)期,從幼穗發(fā)育1 cm時(shí)開始觀察,在幼穗生長至12 cm之前,的頂端小穗與WT沒有區(qū)別(圖2-A—D)。當(dāng)幼穗生長至12 cm時(shí),能夠觀察到明顯的頂端小穗退化現(xiàn)象(圖2-E)。隨著幼穗后續(xù)不斷發(fā)育至成熟,這些小穗停止發(fā)育,逐漸萎縮干癟,在穗子成熟以后,十分容易脫落(圖2-F)。

    A:WT(左)和paa21突變體(右)的株型,bars=10 cm;B:WT(左)和paa21突變體(右)抽穗期的穗子,bars=2 cm;C:WT(左)和paa21突變體(右)成熟的穗子,bars=2 cm;D:來自一個(gè)穗子的一次枝梗,1—16表示paa21中從穗頂部到底部的一次枝梗。E:D圖中對(duì)應(yīng)枝梗的正常小穗和退化小穗統(tǒng)計(jì),F(xiàn):WT與paa21突變體退化一次枝梗數(shù)目所占比例;G:每穗退化粒數(shù)所占比例;H:株高;I:每穗粒數(shù);J:穗長;K:單株產(chǎn)量。數(shù)據(jù)是均值±標(biāo)準(zhǔn)誤;**表示t檢驗(yàn)下顯著差異(P≤0.01)。下同

    A—F:WT(左)和paa21(右)穗發(fā)育不同時(shí)期的表型。依次為1、2、5、7、12和17 cm。bars=1 cm

    2.2 突變體小穗發(fā)生程序性細(xì)胞死亡

    由上述結(jié)果可知,的穗頂端正常發(fā)育出小穗,但是頂端小穗生長到后期的時(shí)候停止生長,種子成熟后風(fēng)干脫落。據(jù)此,推斷頂端小穗的細(xì)胞發(fā)生了程序性細(xì)胞死亡(programmed cell death,PCD)。為了驗(yàn)證這個(gè)猜測(cè),利用臺(tái)盼藍(lán)和伊文思藍(lán)對(duì)WT和10—11 cm穗的頂端小穗進(jìn)行染色。結(jié)果發(fā)現(xiàn),與WT相比,突變體頂端小穗著色變深(圖3)。表明突變體的頂端小穗發(fā)生了程序性細(xì)胞死亡。

    活性氧(reactive oxygen species,ROS)被認(rèn)為是PCD的重要觸發(fā)因素,較高濃度的ROS能夠誘導(dǎo)植物和動(dòng)物組織中的細(xì)胞死亡[22-23]。H2O2是ROS的一種,據(jù)報(bào)道,水稻穗頂部退化多數(shù)是由H2O2積累導(dǎo)致[18-19]。使用DAB染色檢測(cè)頂端小穗中H2O2的含量。與WT相比,在退化的頂端小穗中觀察到更強(qiáng)烈的DAB染色(圖4-A),表明穗中的H2O2含量更高(圖4-A)。

    A:臺(tái)盼藍(lán)染色;B:伊文思藍(lán)染色

    對(duì)幼穗不同發(fā)育階段的H2O2含量進(jìn)行測(cè)定,結(jié)果表明,5 cm長的幼穗中,WT和的H2O2含量沒有明顯差異,但是,隨著幼穗的發(fā)育,幼穗中的H2O2含量較WT逐漸升高。特別是在10和15 cm長的幼穗中(圖4-B)。該結(jié)果與之前的觀察結(jié)果一致,即當(dāng)穗長達(dá)到12 cm時(shí),穗頂端退化開始(圖2),表明ROS水平與穗退化之間可能存在相關(guān)性。

    對(duì)的10 cm幼穗不同穗部的H2O2進(jìn)一步測(cè)定,與基部小穗相比,退化的頂端小穗和中間小穗的H2O2含量更高(圖4-C)。表明頂端小穗中積累的過量ROS是導(dǎo)致其發(fā)生PCD的原因。

    A:WT和paa21頂端小穗DAB染色;B:不同發(fā)育時(shí)期WT和paa21頂端小穗中過氧化氫含量的測(cè)定;C:水稻穗不同部位小穗中過氧化氫含量的測(cè)定,不同字母表示通過方差分析和LSD測(cè)試進(jìn)行平均值比較的結(jié)果存在顯著差異(p<0.05),ap:穗頂端部分,mp:穗中間部分,bp:穗底部

    2.3 穗頂端退化性狀的遺傳分析及PAA21的基因定位

    為了確定的顯隱性以及是否受單基因控制,利用秈稻品種II-32B、9311與突變體進(jìn)行雜交,所有F1植株均表現(xiàn)出野生型表型。F1自交產(chǎn)生的F2群體中出現(xiàn)了野生型表型和穗頂端退化的表型。通過對(duì)F2群體中野生型和突變型的個(gè)體數(shù)目進(jìn)行統(tǒng)計(jì),4個(gè)群體中χ2值均小于3.84,符合3﹕1的分離比(表1),表明穗頂端退化的性狀由一對(duì)隱性核基因控制。

    表1 paa21遺傳分析

    為了克隆,利用和II-32B構(gòu)建的F2群體對(duì)進(jìn)行定位。從F2群體中選取21個(gè)表型個(gè)體,利用本實(shí)驗(yàn)室保存的232對(duì)SSR、InDel標(biāo)記將其初步定位于第2染色體上的標(biāo)記B2-19和B2-21之間。

    為進(jìn)一步精細(xì)定位,利用418株單株,運(yùn)用新開發(fā)的分子標(biāo)記,最終將其定位于H-2-6和H-2-13之間102 kb區(qū)間內(nèi),包含11個(gè)開放閱讀框(open reading frame,ORF),對(duì)其測(cè)序比對(duì),發(fā)現(xiàn)基因第二外顯子發(fā)生一個(gè)C到T的突變,導(dǎo)致丙氨酸突變?yōu)槔i氨酸(圖5)。RAP-DB網(wǎng)站注釋編碼一個(gè)鋁激活蘋果酸轉(zhuǎn)運(yùn)蛋白(aluminum-activated malate transporter,ALMT),是的等位基因[18]。

    圖5 PAA21的圖位克隆

    ALMT7蛋白有7個(gè)跨膜螺旋[18],通過分析WT和蛋白氨基酸序列,結(jié)果表明,突變位點(diǎn)在第4個(gè)跨膜螺旋上(圖6-A)。此外,利用SWISS-MODEL網(wǎng)站對(duì)蛋白三維結(jié)構(gòu)進(jìn)行分析。結(jié)果表明,第二外顯子C到T的突變并未對(duì)蛋白三維結(jié)構(gòu)造成明顯影響(圖6-B—C)。

    2.4 paa21中ROS相關(guān)基因的表達(dá)受到影響

    由于的穗中積累大量的ROS,為了研究ROS產(chǎn)生的分子效應(yīng),對(duì)3個(gè)ROS響應(yīng)標(biāo)志基因、和的表達(dá)進(jìn)行分析。這些基因?qū)儆赪RKY轉(zhuǎn)錄因子家族,在水稻、擬南芥中受ROS的強(qiáng)烈誘導(dǎo)表達(dá)[24-25]。在幼穗發(fā)育至10 cm時(shí),中3個(gè)基因的表達(dá)量較WT大幅升高(圖7-A—C)。

    、編碼液泡加工酶(vacuolar processing enzymes,VPE),是植物程序性細(xì)胞死亡的關(guān)鍵調(diào)控因子[26-27]。RT-qPCR結(jié)果表明,幼穗發(fā)育至10 cm時(shí),中和的表達(dá)量較WT顯著升高(圖7-D—E)。

    、、編碼過氧化氫酶,過氧化氫酶能夠催化H2O2形成H2O和O2,在ROS清除系統(tǒng)中發(fā)揮重要作用。檢測(cè)發(fā)育至5和10 cm的幼穗中這3個(gè)基因的表達(dá)情況。結(jié)果表明,5 cm幼穗中,3個(gè)基因在突變體和WT之間的表達(dá)量相似(圖7-F—H)。而在10 cm幼穗中,3個(gè)基因的表達(dá)量較WT明顯上升(圖7-F—H)。

    此外,還對(duì)10 cm幼穗中過氧化氫酶(CAT)的活性進(jìn)行了測(cè)定。發(fā)現(xiàn)中CAT的活性明顯下降。說明突變體中ROS的清除能力下降(圖7-I)。

    M1—M6表示PAA21蛋白的跨膜螺旋結(jié)構(gòu),紅框指示突變位點(diǎn)

    綜上所述,幼穗在發(fā)育后期頂端小穗中積累過量的ROS,產(chǎn)生程序性細(xì)胞死亡,最終導(dǎo)致頂端小穗發(fā)生退化。

    3 討論

    3.1 paa21是paab1-1新的等位突變體

    水稻穗頂端退化會(huì)造成產(chǎn)量上的損失。迄今為止,已經(jīng)報(bào)道了多個(gè)導(dǎo)致水稻穗部頂端退化的基因,如、、、、、和[11, 18-20, 28-31]。這些基因涉及內(nèi)質(zhì)網(wǎng)應(yīng)激、營養(yǎng)物質(zhì)運(yùn)輸、ROS穩(wěn)態(tài)調(diào)節(jié)等多個(gè)生理過程,穎花退化的分子機(jī)制和生理機(jī)理還有待進(jìn)一步豐富[32]。

    本研究在EMS突變體庫中篩選獲得一個(gè)以武運(yùn)粳30號(hào)為背景的穗頂端退化突變體。的穗頂端會(huì)發(fā)生嚴(yán)重的退化現(xiàn)象,導(dǎo)致其穗子變短、每穗粒數(shù)減少、單株產(chǎn)量下降(圖1)。退化小穗數(shù)占每穗總穎花數(shù)的比例約為20%,退化小穗主要分布在穗頂端,并且處于穗軸上部的一次枝梗退化的穎花占比大(圖1)。發(fā)現(xiàn)的幼穗發(fā)育到12 cm時(shí),頂端小穗發(fā)生退化,表明參與維持水稻穗發(fā)育(圖2)。遺傳分析結(jié)果表明,表型受一對(duì)隱性核基因控制(表1)。基因定位結(jié)果表明,突變體在的第二個(gè)外顯子中具有單核苷酸替換(C突變A),導(dǎo)致在PAA21第129個(gè)殘基處丙氨酸(Ala)被纈氨酸(Val)替換,但并沒有造成蛋白三維結(jié)構(gòu)發(fā)生明顯改變(圖5和圖6)。編碼一個(gè)鋁激活的蘋果酸轉(zhuǎn)運(yùn)蛋白,預(yù)測(cè)有7個(gè)跨膜螺旋。的突變位點(diǎn)在第4個(gè)跨膜螺旋上(圖6-A),推斷該突變可能導(dǎo)致跨膜結(jié)構(gòu)改變,無法正常轉(zhuǎn)運(yùn)蘋果酸,產(chǎn)生了穗頂端退化的表型。

    圖7 活性氧、PCD相關(guān)基因的表達(dá)量及過氧化氫酶活性

    RA-PDB網(wǎng)站注釋編碼一個(gè)鋁激活的蘋果酸轉(zhuǎn)運(yùn)蛋白OsALMT。Heng等[18]利用以Kittake為背景的穗頂端退化突變體研究了該基因通過參與蘋果酸的轉(zhuǎn)運(yùn)維持穗發(fā)育的功能,并將其命名為。因此,是新的等位基因。與本研究不同的是,是一個(gè)顯性突變,而本研究中的則是一個(gè)隱性突變。的突變位點(diǎn)位于第二個(gè)內(nèi)含子和第三個(gè)外顯子之間的剪接位點(diǎn),導(dǎo)致mRNA剪接發(fā)生改變,產(chǎn)生了2種錯(cuò)誤的轉(zhuǎn)錄本[18]。而的突變位點(diǎn)則是第二外顯子上,僅造成一個(gè)氨基酸的替換。此外,彭永彬[33]利用穗頂端退化突變體,與本研究類似,也是一個(gè)隱性突變。在第三外顯子上發(fā)生了一個(gè)單堿基替換,導(dǎo)致對(duì)應(yīng)編碼氨基酸由精氨酸變?yōu)橘嚢彼?。因此,、和中顯隱性的差異可能是突變方式不同導(dǎo)致的。

    3.2 ROS過度積累引起的PCD是paa21穗頂端退化的原因

    活性氧是PCD的重要觸發(fā)因素,ROS的過度積累使細(xì)胞膜高度氧化,進(jìn)而影響細(xì)胞通透性,最終導(dǎo)致細(xì)胞死亡[23, 34-35]。本研究臺(tái)盼藍(lán)和伊文思藍(lán)染色結(jié)果表明,頂端退化小穗中發(fā)生了PCD(圖3)。DAB染色、H2O2含量測(cè)定發(fā)現(xiàn),頂端退化小穗中積累了大量的ROS(圖4)。此外,穗部的PCD相關(guān)基因和ROS清除相關(guān)基因的表達(dá)水平顯著增加,過氧化氫酶的活性也大幅下降(圖7)。結(jié)果表明,由于ROS在穗頂端的過度積累而引起的PCD可能是穗頂端退化的原因。綜上所述,推測(cè)的突變導(dǎo)致ROS過度積累,最終誘導(dǎo)細(xì)胞死亡。

    4 結(jié)論

    EMS誘變獲得一個(gè)水稻穗頂端退化突變體,突變表型受一對(duì)隱性核基因控制,為突變所致。穗發(fā)育后期積累大量ROS,誘導(dǎo)細(xì)胞死亡,產(chǎn)生穗頂端退化表型。

    [1] XING Y, ZHANG Q. Genetic and molecular bases of rice yield. Annual Review of Plant Biology, 2010, 61: 421-442.

    [2] TEO Z W N, SONG S, WANG Y Q, LIU J, YU H. New insights into the regulation of inflorescence architecture. Trends in Plant Science, 2014, 19(3): 158-165.

    [3] IKEDA-KAWAKATSU K, YASUNO N, OIKAWA T, IIDA S, NAGATO Y, MAEKAWA M, KYOZUKA J. Expression level ofdetermines rice inflorescence form through control of cell proliferation in the meristem. Plant Physiology, 2009, 150(2): 736-747.

    [4] IKEDA K, ITO M, NAGASAWA N, KYOZUKA J, NAGATO Y. Rice, encoding an F-box protein, regulates meristem fate. The Plant Journal, 2007, 51(6): 1030-1040.

    [5] IKEDA K, NAGASAWA N, NAGATO Y.temporally regulates meristem identity in rice. Developmental Biology, 2005, 282(2): 349-360.

    [6] LI M, TANG D, WANG K, WU X, LU L, YU H, GU M, YAN C, CHENG Z. Mutations in the F-box geneimprove the panicle architecture and enhance the grain yield in rice. Plant Biotechnology Journal, 2011, 9(9): 1002-1013.

    [7] HUANG X, QIAN Q, LIU Z, SUN H, HE S, LUO D, XIA G, CHU C, LI J, FU X. Natural variation at thelocus enhances grain yield in rice. Nature Genetics, 2009, 41(4): 494-497.

    [8] KOMATSU M, MAEKAWA M, SHIMAMOTO K, KYOZUKA J. Theandgenes determine the inflorescence architecture of rice by controlling rachis-branch and spikelet development. Developmental Biology, 2001, 231(2): 364-373.

    [9] WANG H, TANG S, ZHI H, XING L, ZHANG H, TANG C, WANG E, ZHAO M, JIA G, FENG B, DIAO X. The boron transporter SiBOR1 functions in cell wall integrity, cellular homeostasis, and panicle development in foxtail millet. the Crop Journal, 2022, 10(2): 342-353.

    [10] PEI Y, DENG Y, ZHANG H, ZHANG Z, LIU J, CHEN Z, CAI D, LI K, DU Y, ZANG J, XIN P, CHU J, CHEN Y, ZHAO L, LIU J, CHEN H. EAR APICAL DEGENERATION1 regulates maize ear development by maintaining malate supply for apical inflorescence. The Plant Cell, 2022, 34(6): 2222-2241.

    [11] BAI J, ZHU X, WANG Q, ZHANG J, CHEN H, DONG G, ZHU L, ZHENG H, XIE Q, NIAN J, CHEN F, FU Y, QIAN Q, ZUO J. Riceencodes a suppressor of cAMP receptor-like protein that is important for actin organization and panicle development. Plant Physiology, 2015, 169(2): 1179-1191.

    [12] 陳惠哲, 朱德峰, 林賢青, 張玉屏. 穗肥施氮量對(duì)兩優(yōu)培九枝梗及穎花分化和退化的影響. 浙江農(nóng)業(yè)學(xué)報(bào), 2008, 20(3): 181-185.

    CHEN H Z, ZHU D F, LIN X Q, ZHANG Y P. Effect of nitrogen levels in spike stage on differentiation and degeneration of branches and spikelet of hybrid rice cultivar Liangyoupeijiu. Acta Agriculturae Zhejiangensis, 2008, 20(3): 181-185. (in Chinese)

    [13] 王亞梁, 張玉屏, 朱德峰, 向鏡, 武輝, 陳惠哲, 張義凱. 水稻穗分化期高溫脅迫對(duì)穎花退化及籽粒充實(shí)的影響. 作物學(xué)報(bào), 2016, 42(9): 1402-1410.

    WANG Y L, ZHANG Y P, ZHU D F, XIANG J, WU H, CHEN H Z, ZHANG Y K. Effect of heat stress on spikelet degeneration and grain filling at panicle initiation period of rice. Acta Agronomica Sinica, 2016, 42(9): 1402-1410. (in Chinese)

    [14] 張興元, 羅勝, 王敏, 叢楠, 趙志超, 程治軍. 與SP1互作的水稻穗頂部退化基因的精細(xì)定位. 中國農(nóng)業(yè)科學(xué), 2015, 48(12): 2287-2295.

    ZHANG X Y, LUO S, WANG M, CONG N, ZHAO Z C, CHENG Z J. Fine mapping of rice panicle apical abortion geneinteracting with. Scientia Agricultura Sinica, 2015, 48(12): 2287-2295. (in Chinese)

    [15] TAN C J, SUN Y J, XU H S, YU S B. Identification of quantitative trait locus and epistatic interaction for degenerated spikelets on the top of panicle in rice. Plant Breeding, 2011, 130(2): 177-184.

    [16] 徐華山, 孫永建, 周紅菊, 余四斌. 構(gòu)建水稻優(yōu)良恢復(fù)系背景的重疊片段代換系及其效應(yīng)分析. 作物學(xué)報(bào), 2007, 33(6): 979-986.

    XU H S, SUN Y J, ZHOU H J, YU S B. Development and characterization of contiguous segment substitution lines with background of an elite restorer line. Acta Agronomica Sinica, 2007, 33(6): 979-986. (in Chinese)

    [17] CHENG Z J, MAO B G, GAO S W, ZHANG L, WANG J L, LEI C L, ZHANG X, WU F Q, GUO X P, WAN J M. Fine mapping of, a gene controlling panicle apical development in rice. Journal of Integrative Plant Biology, 2011, 53(9): 710-718.

    [18] HENG Y, WU C, LONG Y, LUO S, MA J, CHEN J, LIU J, ZHANG H, REN Y, WANG M, TAN J, ZHU S, WANG J, LEI C, ZHANG X, GUO X, WANG H, CHENG Z, WAN J. OsALMT7 maintains panicle size and grain yield in rice by mediating malate transport. The Plant Cell, 2018, 30(4): 889-906.

    [19] ZAFAR S A, PATIL S B, UZAIR M, FANG J, ZHAO J, GUO T, YUAN S, UZAIR M, LUO Q, SHI J, SCHREIBER L, LI X.() encodes a cystathionine β-synthase domain containing protein required for anther cuticle and panicle development in rice. New Phytologist, 2020, 225(1): 356-375.

    [20] PENG Y B, HOU F X, BAI Q, XU P Z, LIAO Y X, ZHANG H Y, GU C J, DENG X S, WU T K, CHEN X Q, ALI A, WU X J. Rice calcineurin b-like protein-interacting protein kinase 31 (OsCIPK31) is involved in the development of panicle apical spikelets. Frontiers in Plant Science, 2018, 9.

    [21] 王中豪. 水稻鈣氫離子交換蛋白基因的圖位克隆和功能分析[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2021.

    WANG Z H. Map-based cloning and functional analysis of the Ca2+/H+exchanger genein rice[D]. Beijing: Chinese Academy of Agricultural Sciences, 2021. (in Chinese)

    [22] LOOR G, KONDAPALLI J, SCHRIEWER J M, CHANDEL N S, VANDEN HOEK T L, SCHUMACKER P T. Menadione triggers cell death through ROS-dependent mechanisms involving PARP activation without requiring apoptosis. Free Radical Biology and Medicine, 2010, 49(12): 1925-1936.

    [23] MITTLER R. ROS are good. Trends in Plant Science, 2017, 22(1): 11-19.

    [24] LI Z, MO W, JIA L, XU Y C, TANG W, YANG W, GUO Y L, LIN R. Rice FLUORESCENT1 is involved in the Regulation of Chlorophyll. Plant Cell Physiology, 2019, 60(19): 2307-2318.

    [25] SU T, WANG P, LI H, ZHAO Y, LU Y, DAI P, REN T, WANG X, LI X, SHAO Q, ZHAO D, ZHAO Y, MA C. Thecatalase triple mutant reveals important roles of catalases and peroxisome- derived signaling in plant development. Journal of Integrative Plant Biology, 2018, 60(7): 591-607.

    [26] DENG M, BIAN H, XIE Y, KIM Y, WANG W, LIN E, ZENG Z, GUO F, PAN J, HAN N, WANG J, QIAN Q, ZHU M. Bcl-2 suppresses hydrogen peroxide-induced programmed cell death via OsVPE2 and OsVPE3, but not via OsVPE1 and OsVPE4, in rice. The FEBS Journal, 2011, 278(24): 4797-4810.

    [27] LU W Y, DENG M J, GUO F, WANG M Q, ZENG Z H, HAN N, YANG Y N, ZHU M Y, BIAN H W. Suppression ofenhances salt tolerance by attenuating vacuole rupture during programmed cell death and affects stomata development in rice. Rice, 2016, 9.

    [28] DAI D, ZHANG H, HE L, CHEN J, DU C, LIANG M, ZHANG M, WANG H, MA L. Panicle apical abortion 7 regulates panicle development in rice (L.). International Journal of Molecular Sciences, 2022, 23(16): 9487.

    [29] HU P, TAN Y Q, WEN Y, FANG Y X, WANG Y Y, WU H, WANG J G, WU K X, CHAI B Z, ZHU L, ZHANG G H, GAO Z Y, REN D Y, ZENG D L, SHEN L, XUE D W, QIAN Q, HU J. LMPA regulates lesion mimic leaf and panicle development through ROS-induced PCD in rice. Frontiers in Plant Science, 2022, 13.

    [30] YANG F, XIONG M, HUANG M, LI Z, WANG Z, ZHU H, CHEN R, LU L, CHENG Q, WANG Y, TANG J, ZHUANG H, LI Y. Panicle apical abortion 3 controls panicle development and seed size in rice. Rice, 2021, 14(1): 68.

    [31] WANG Q L, SUN A Z, CHEN S T, CHEN L S, GUO F Q. SPL6 represses signalling outputs of ER stress in control of panicle cell death in rice. Nature Plants, 2018, 4(5): 280-288.

    [32] ALI A, XU P, RIAZ A, WU X. Current advances in molecular mechanisms and physiological basis of panicle degeneration in rice. International Journal of Molecular Sciences, 2019, 20(7): 1613.

    [33] 彭永彬. 水稻穗頂退化突變體和的基因克隆與功能分析[D]. 成都: 四川農(nóng)業(yè)大學(xué), 2018.

    PENG Y B. Cloning and functional charactization ofand

    [34] WASZCZAK C, CARMODY M, KANGASJ?RVI J. Reactive oxygen species in plant signaling. Annual Review of Plant Biology, 2018, 69(1): 209-236.

    [35] MHAMDI A, VAN BREUSEGEM F. Reactive oxygen species in plant development. Development, 2018, 145(15).

    Phenotypic analysis and gene cloning of rice panicle apical abortion mutant

    HE Lei1, LU Kai1, ZHAO ChunFang1, YAO Shu1, ZHOU LiHui1, ZHAO Ling1, CHEN Tao1, ZHU Zhen1, ZHAO QingYong1, LIANG WenHua1, WANG CaiLin1, ZHU Li2, ZHANG YaDong1

    1Institute of Food Crops, Jiangsu Academy of Agricultural Sciences/East China Branch of National Center of Technology Innovation for Saline-Alkali Tolerant Rice/Jiangsu High Quality Rice R&D Center/Nanjing Branch of China National Center for Rice Improvement/Key laboratory of Jiangsu Province for Agrobiology, Nanjing 210014;2China National Rice Research Institute/State Key Laboratory of Rice Biology, Hangzhou 310006

    Rice panicle apical abortion affects yield. Identification and cloning of genes related to rice panicle apical abortion can enrich the molecular mechanism of rice panicle development regulation, and provide theoretical basis and genetic resources for rice high-yield molecular design breeding.Here, a stably inherited() mutant was screened from EMS mutant library of the japonica rice variety "Wuyunjing 30". Agronomic traits, such as ratio of degraded primary branches, degraded apical spikelets, grains per panicle, plant height, panicle length, and grain yield per plant, were statistically analyzed. Trypan blue and Evans blue staining were used to detect whether programmed cell death occurred in the apical spikelets. H2O2content in young panicles at different development stages and different panicle parts of WT andwas determined. Genetic analysis was carried out by reciprocal cross ofwith indica rice II-32B and 9311 respectively. The F2population constructed by crossingwith indica rice II-32B was used for gene mapping and cloning. The three-dimensional structure of wild-type andproteins were predicted using SWISS-MODEL website. The expression levels of ROS response marker genes, programmed cell death related genes and catalase related genes were analyzed by RT-qPCR.produced panicle apical abortion phenotype and the degenerated spikelets were mainly located on the primary branches at the apical panicle. The plant height, grain number per panicle, panicle length and grain yield per plant ofwere lower than those of WT. After observing the young panicles at different development stages, we found that themutant had a panicle apical abortion phenotype when panicle developed to 12 cm. Trypan blue and Evans blue staining results showed that the apical spikelets of themutant had programmed cell death. Stronger DAB staining was observed in the degenerated apical spikelets ofthan WT. The results of H2O2content determination showed that higher level of ROS was accumulated in panicle ofcompared with WT. Genetic analysis suggested thatmutant phenotype is controlled by a pair of recessive nuclear genes. The results of map-based cloning showed that a C to T mutation occurred in the second exon of, resulting in the mutation of alanine to valine. This gene encodes an aluminum activated malate transporter, ALMT7. The mutation site was located at the fourth transmembrane helix. SWISS-MODEL prediction results showed that the mutation site did not significantly affect the three-dimensional structure of the mutant protein. The expression level of ROS response marker genes,andwas significantly higher than that in WT when the young spike developed to 10 cm. Compared with WT, the expression level of programmed cell death related genesandincreased significantly in. The expression level of,andwhich encode catalase in 10 cm young panicle ofwas significantly higher than that of WT. The activity of CAT in10 cm young spikelet was significantly lower than that of WT.accumulate excess ROS in the apical spikelet at late stage of panicle development, resulting in programmed cell death, which eventually leads to the degeneration of the apical spikelet. These results lay a good foundation for further enriching the genetic regulatory network of panicle development.

    rice; panicle apical abortion; reactive oxygen species; programmed cell death

    2022-09-27;

    2022-10-24

    江蘇省種業(yè)振興揭榜掛帥項(xiàng)目(JBGS[2021]001)、江蘇省重點(diǎn)研發(fā)計(jì)劃(BE2022336)、現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)資金(CARS-01)

    赫磊,E-mail:helei@jaas.ac.cn。通信作者朱麗,E-mail:zhuli05@caas.cn。通信作者張亞東,E-mail:zhangyd@jaas.ac.cn

    (責(zé)任編輯 李莉)

    猜你喜歡
    小穗枝梗幼穗
    多穗型與大穗型超級(jí)稻品種幼穗分化期碳氮代謝與產(chǎn)量構(gòu)成
    播期對(duì)機(jī)插雜交秈稻不同莖蘗部位稻穗枝梗數(shù)和穎花數(shù)的影響
    粳稻鹽粳188 航天誘變后代穗部性狀多樣性及相關(guān)性分析
    北方水稻(2021年2期)2021-04-26 08:25:56
    水稻不同穗型品種幼穗分化相關(guān)基因表達(dá)分析
    最新研究表明,不育小穗有助于提高高粱等禾本植物的產(chǎn)量
    我科學(xué)家成功克隆水稻小穗發(fā)育新基因
    水稻小穗發(fā)育新基因克隆成功
    四分體期低溫脅迫對(duì)2個(gè)小麥品種幼穗結(jié)實(shí)率及生理特性的影響
    不同穗型粳稻品種(系)稻曲病粒在穗內(nèi)粒位上的分布
    高粱幼穗再生體系建立的研究
    亚洲精品亚洲一区二区| 麻豆久久精品国产亚洲av| 欧美另类亚洲清纯唯美| 色综合站精品国产| 亚洲av二区三区四区| 日日摸夜夜添夜夜添av毛片 | 欧美绝顶高潮抽搐喷水| 亚洲精华国产精华精| 免费av毛片视频| 不卡一级毛片| 国产在线男女| 欧美色欧美亚洲另类二区| 制服丝袜大香蕉在线| 少妇人妻一区二区三区视频| 日韩精品中文字幕看吧| 精品久久久久久久久久久久久| 国产精华一区二区三区| 亚洲av成人av| 神马国产精品三级电影在线观看| 国产伦一二天堂av在线观看| 亚洲avbb在线观看| 国产成人影院久久av| 九色成人免费人妻av| 国产av在哪里看| 又黄又爽又免费观看的视频| 久久精品国产亚洲av天美| 桃红色精品国产亚洲av| 亚洲精品粉嫩美女一区| 日韩欧美国产一区二区入口| 日韩国内少妇激情av| www.熟女人妻精品国产| 午夜精品在线福利| 久久久久久大精品| 午夜福利18| 窝窝影院91人妻| 国产亚洲精品久久久com| 久久欧美精品欧美久久欧美| 少妇熟女aⅴ在线视频| 最近最新中文字幕大全电影3| 欧美日韩福利视频一区二区| 制服丝袜大香蕉在线| 国产91精品成人一区二区三区| 我的老师免费观看完整版| 在线观看舔阴道视频| 直男gayav资源| 日韩 亚洲 欧美在线| 欧美成人性av电影在线观看| 免费电影在线观看免费观看| 国内少妇人妻偷人精品xxx网站| 搞女人的毛片| 国产高清三级在线| 男人狂女人下面高潮的视频| 久久欧美精品欧美久久欧美| 国产白丝娇喘喷水9色精品| 男女那种视频在线观看| 亚洲人成伊人成综合网2020| 日本黄色片子视频| 精品午夜福利视频在线观看一区| 欧美成狂野欧美在线观看| 日韩成人在线观看一区二区三区| 99久久九九国产精品国产免费| 亚洲在线自拍视频| 首页视频小说图片口味搜索| 最后的刺客免费高清国语| 一区二区三区四区激情视频 | 免费在线观看影片大全网站| 99国产精品一区二区三区| 观看免费一级毛片| 免费在线观看成人毛片| 久久午夜福利片| 午夜福利免费观看在线| 国产蜜桃级精品一区二区三区| 大型黄色视频在线免费观看| av福利片在线观看| 在线播放无遮挡| 尤物成人国产欧美一区二区三区| 免费在线观看亚洲国产| 中国美女看黄片| 国产精品女同一区二区软件 | 亚洲美女视频黄频| 亚洲av成人不卡在线观看播放网| h日本视频在线播放| 少妇的逼水好多| 在线播放国产精品三级| 日本撒尿小便嘘嘘汇集6| 久久久久九九精品影院| 免费观看人在逋| 波多野结衣高清作品| 欧美bdsm另类| 欧美+亚洲+日韩+国产| av福利片在线观看| 91狼人影院| 亚洲三级黄色毛片| 亚洲精品456在线播放app | 日韩高清综合在线| 免费看美女性在线毛片视频| 免费高清视频大片| 男女下面进入的视频免费午夜| 亚洲专区中文字幕在线| 小说图片视频综合网站| 一本综合久久免费| 免费观看精品视频网站| 亚洲精品乱码久久久v下载方式| 97超级碰碰碰精品色视频在线观看| 青草久久国产| 欧美激情在线99| 亚洲性夜色夜夜综合| 日本精品一区二区三区蜜桃| 在线观看一区二区三区| 午夜免费成人在线视频| 免费在线观看影片大全网站| 亚洲专区国产一区二区| 亚洲,欧美,日韩| 国产伦精品一区二区三区四那| 中文在线观看免费www的网站| 首页视频小说图片口味搜索| 18禁黄网站禁片午夜丰满| 嫁个100分男人电影在线观看| 69人妻影院| 久久久精品欧美日韩精品| 亚洲中文字幕日韩| 成年版毛片免费区| 亚洲,欧美,日韩| 人人妻人人看人人澡| 看免费av毛片| 最近最新中文字幕大全电影3| 亚洲久久久久久中文字幕| 欧美中文日本在线观看视频| bbb黄色大片| 日韩亚洲欧美综合| 高清毛片免费观看视频网站| 久久亚洲真实| 午夜亚洲福利在线播放| 色综合站精品国产| 九色成人免费人妻av| 亚州av有码| 我的女老师完整版在线观看| 国内揄拍国产精品人妻在线| 色在线成人网| 此物有八面人人有两片| 亚洲 欧美 日韩 在线 免费| 免费av不卡在线播放| 一进一出抽搐动态| 亚洲精品粉嫩美女一区| 天堂动漫精品| 日韩欧美三级三区| 欧美日韩亚洲国产一区二区在线观看| aaaaa片日本免费| 国产真实乱freesex| 亚洲精品久久国产高清桃花| 男女视频在线观看网站免费| 天美传媒精品一区二区| 嫩草影院精品99| 亚洲最大成人av| or卡值多少钱| 一个人免费在线观看的高清视频| 色哟哟·www| www.999成人在线观看| 88av欧美| 波多野结衣高清作品| 男女床上黄色一级片免费看| 色吧在线观看| 久久中文看片网| 亚洲激情在线av| 久久精品国产清高在天天线| 日韩欧美一区二区三区在线观看| 一级a爱片免费观看的视频| 免费看a级黄色片| 午夜免费成人在线视频| 日韩精品青青久久久久久| 欧美一区二区精品小视频在线| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美 国产精品| 三级国产精品欧美在线观看| 国产精品一区二区性色av| 国产精品一区二区性色av| 亚洲精品在线美女| 国产一区二区三区视频了| 国产一区二区三区视频了| 国产精品一区二区性色av| 成人av一区二区三区在线看| 99久久久亚洲精品蜜臀av| 久久久久久九九精品二区国产| 12—13女人毛片做爰片一| 一进一出抽搐gif免费好疼| 午夜福利高清视频| 97碰自拍视频| 人人妻人人澡欧美一区二区| 久久国产乱子伦精品免费另类| 亚洲国产精品sss在线观看| eeuss影院久久| 九九久久精品国产亚洲av麻豆| 国产伦精品一区二区三区视频9| 国产蜜桃级精品一区二区三区| 香蕉av资源在线| 精品午夜福利在线看| 变态另类丝袜制服| 午夜福利免费观看在线| 欧美性感艳星| 一个人免费在线观看电影| 嫩草影院新地址| 国产精品久久久久久人妻精品电影| 免费搜索国产男女视频| 精品久久久久久久末码| 免费观看的影片在线观看| 午夜精品久久久久久毛片777| bbb黄色大片| 国产69精品久久久久777片| 久久久国产成人精品二区| 中文字幕免费在线视频6| 一本精品99久久精品77| 国产色婷婷99| 精品久久久久久久人妻蜜臀av| 91狼人影院| 国产白丝娇喘喷水9色精品| 久久香蕉精品热| 999久久久精品免费观看国产| 久久热精品热| 少妇人妻一区二区三区视频| 搡女人真爽免费视频火全软件 | 欧美xxxx黑人xx丫x性爽| 久久久久九九精品影院| 能在线免费观看的黄片| 日本黄色片子视频| 亚洲不卡免费看| 亚洲在线观看片| 日韩欧美在线乱码| 在线十欧美十亚洲十日本专区| 婷婷精品国产亚洲av| 国产精品一区二区性色av| av中文乱码字幕在线| 最新中文字幕久久久久| 88av欧美| av在线天堂中文字幕| 欧美成人a在线观看| 麻豆国产av国片精品| 最新在线观看一区二区三区| 亚洲欧美精品综合久久99| 十八禁网站免费在线| 亚洲国产欧洲综合997久久,| 黄色视频,在线免费观看| 欧美成人a在线观看| 99在线视频只有这里精品首页| 国产视频内射| 亚洲美女搞黄在线观看 | 国产精品,欧美在线| 床上黄色一级片| 在线观看美女被高潮喷水网站 | 内射极品少妇av片p| 国产毛片a区久久久久| 在线天堂最新版资源| 男人舔奶头视频| 亚洲片人在线观看| 九色成人免费人妻av| 我要搜黄色片| 男插女下体视频免费在线播放| 亚洲人成伊人成综合网2020| 老司机福利观看| 中文字幕久久专区| 伊人久久精品亚洲午夜| 亚洲国产高清在线一区二区三| 欧美又色又爽又黄视频| 久久这里只有精品中国| 亚洲精品粉嫩美女一区| 午夜精品在线福利| 久久久久久久亚洲中文字幕 | 成人美女网站在线观看视频| 国产 一区 欧美 日韩| 中文字幕久久专区| 欧美日韩国产亚洲二区| 九九在线视频观看精品| 久久欧美精品欧美久久欧美| 久久这里只有精品中国| 久久久久久国产a免费观看| 日韩av在线大香蕉| 最近在线观看免费完整版| 又爽又黄无遮挡网站| 精品久久久久久久久久免费视频| 婷婷精品国产亚洲av| 免费av毛片视频| 成年免费大片在线观看| 亚洲最大成人中文| 88av欧美| 亚洲欧美精品综合久久99| 国产成人av教育| 亚洲人成网站在线播放欧美日韩| 亚洲 欧美 日韩 在线 免费| 亚洲国产日韩欧美精品在线观看| 极品教师在线视频| 熟女人妻精品中文字幕| 性欧美人与动物交配| 可以在线观看的亚洲视频| 欧美xxxx性猛交bbbb| 老鸭窝网址在线观看| 国产美女午夜福利| 亚洲人成网站在线播| 久久久久久久久中文| 午夜福利在线观看免费完整高清在 | 国产精品一及| 亚洲精品在线观看二区| 高清日韩中文字幕在线| 国产精品女同一区二区软件 | or卡值多少钱| 午夜福利在线观看吧| 亚洲第一欧美日韩一区二区三区| 亚洲三级黄色毛片| 国产精品一区二区三区四区免费观看 | av在线蜜桃| 搡女人真爽免费视频火全软件 | 色精品久久人妻99蜜桃| 非洲黑人性xxxx精品又粗又长| 亚洲欧美精品综合久久99| 日日摸夜夜添夜夜添av毛片 | 国产免费男女视频| 久久久久久九九精品二区国产| 在线观看一区二区三区| 俄罗斯特黄特色一大片| 久9热在线精品视频| 丰满人妻熟妇乱又伦精品不卡| 搞女人的毛片| 神马国产精品三级电影在线观看| 午夜福利在线在线| 国产一区二区在线观看日韩| www.色视频.com| 久久99热这里只有精品18| 亚洲av一区综合| 久久久久性生活片| 啦啦啦韩国在线观看视频| 亚洲电影在线观看av| 99热只有精品国产| 日韩av在线大香蕉| 亚洲av免费高清在线观看| 内地一区二区视频在线| 国产亚洲精品av在线| 午夜福利在线观看免费完整高清在 | 成年女人毛片免费观看观看9| 99精品在免费线老司机午夜| 精品国产亚洲在线| 黄片小视频在线播放| 日本一二三区视频观看| 精品久久久久久久久久久久久| 久久久久国产精品人妻aⅴ院| 美女大奶头视频| 精品免费久久久久久久清纯| 麻豆av噜噜一区二区三区| 国产亚洲精品av在线| 动漫黄色视频在线观看| 亚洲自偷自拍三级| 青草久久国产| 成人亚洲精品av一区二区| 91狼人影院| 少妇人妻一区二区三区视频| а√天堂www在线а√下载| 精品人妻熟女av久视频| 国产在视频线在精品| 欧美不卡视频在线免费观看| 午夜福利18| 超碰av人人做人人爽久久| 无人区码免费观看不卡| 可以在线观看毛片的网站| 精品人妻偷拍中文字幕| 有码 亚洲区| 欧美成狂野欧美在线观看| 99久久成人亚洲精品观看| 看十八女毛片水多多多| 婷婷六月久久综合丁香| 色精品久久人妻99蜜桃| 欧美日本亚洲视频在线播放| 露出奶头的视频| av女优亚洲男人天堂| 欧美一级a爱片免费观看看| 精品国产亚洲在线| 国产伦精品一区二区三区视频9| 啪啪无遮挡十八禁网站| 特大巨黑吊av在线直播| 国内精品一区二区在线观看| 国产主播在线观看一区二区| 一边摸一边抽搐一进一小说| 国产在视频线在精品| 国产精品亚洲一级av第二区| 亚洲av美国av| 国产欧美日韩一区二区精品| 亚洲av中文字字幕乱码综合| 此物有八面人人有两片| 日本免费a在线| 国产精品久久电影中文字幕| 99视频精品全部免费 在线| 一二三四社区在线视频社区8| 欧美xxxx性猛交bbbb| 小蜜桃在线观看免费完整版高清| 最好的美女福利视频网| 成人午夜高清在线视频| 婷婷丁香在线五月| 欧美最新免费一区二区三区 | 舔av片在线| 最新在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 久久精品久久久久久噜噜老黄 | 91狼人影院| 九九久久精品国产亚洲av麻豆| 久久午夜亚洲精品久久| 国产av麻豆久久久久久久| 午夜福利在线在线| 日本a在线网址| 天堂影院成人在线观看| 成人毛片a级毛片在线播放| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 国产精品日韩av在线免费观看| 国产成人av教育| 小说图片视频综合网站| 他把我摸到了高潮在线观看| 国产精品久久久久久人妻精品电影| 观看免费一级毛片| 九九在线视频观看精品| 美女高潮喷水抽搐中文字幕| 久9热在线精品视频| 亚洲国产日韩欧美精品在线观看| 色噜噜av男人的天堂激情| 最后的刺客免费高清国语| 99精品在免费线老司机午夜| 婷婷色综合大香蕉| 最后的刺客免费高清国语| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 亚洲精品影视一区二区三区av| 免费观看精品视频网站| 亚洲国产精品999在线| 久久天躁狠狠躁夜夜2o2o| 99精品在免费线老司机午夜| 少妇人妻一区二区三区视频| 一边摸一边抽搐一进一小说| av在线老鸭窝| 免费看a级黄色片| 一区二区三区免费毛片| av视频在线观看入口| 国产精品一区二区免费欧美| 波野结衣二区三区在线| 性色av乱码一区二区三区2| 国产成人a区在线观看| 国产精品影院久久| av中文乱码字幕在线| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 久久精品国产亚洲av香蕉五月| 国产伦人伦偷精品视频| 午夜精品久久久久久毛片777| 亚洲av不卡在线观看| 看免费av毛片| 国产免费一级a男人的天堂| 成人午夜高清在线视频| 国产在线男女| 精品一区二区免费观看| 无人区码免费观看不卡| 熟女电影av网| 超碰av人人做人人爽久久| 脱女人内裤的视频| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 精品久久久久久久人妻蜜臀av| 免费看a级黄色片| 久久亚洲真实| 精品国内亚洲2022精品成人| 他把我摸到了高潮在线观看| 国内精品久久久久精免费| 成年女人永久免费观看视频| 国产亚洲精品av在线| 亚洲av第一区精品v没综合| 亚洲色图av天堂| 青草久久国产| 国产精品久久久久久亚洲av鲁大| 精品日产1卡2卡| 色综合亚洲欧美另类图片| 国产成人a区在线观看| 精品久久国产蜜桃| 国产探花极品一区二区| 久久久久国内视频| 好男人电影高清在线观看| 欧美丝袜亚洲另类 | av欧美777| av黄色大香蕉| 长腿黑丝高跟| 日韩欧美在线乱码| 51午夜福利影视在线观看| 伦理电影大哥的女人| 亚洲精品久久国产高清桃花| 日本黄大片高清| 黄片小视频在线播放| 国产精品精品国产色婷婷| 久久久久久久精品吃奶| 三级毛片av免费| av黄色大香蕉| 免费人成在线观看视频色| 一本精品99久久精品77| 真人做人爱边吃奶动态| 不卡一级毛片| 成人毛片a级毛片在线播放| 日韩欧美一区二区三区在线观看| 亚洲成人中文字幕在线播放| 久久久久国内视频| 久久精品国产亚洲av涩爱 | 99在线人妻在线中文字幕| 丁香六月欧美| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩东京热| 51午夜福利影视在线观看| 国产视频一区二区在线看| 日韩成人在线观看一区二区三区| 美女大奶头视频| 亚洲最大成人中文| 亚洲中文字幕日韩| 日韩欧美三级三区| 日韩av在线大香蕉| 蜜桃亚洲精品一区二区三区| 在线国产一区二区在线| 三级男女做爰猛烈吃奶摸视频| 日韩欧美精品v在线| 亚洲最大成人av| 国产欧美日韩一区二区精品| av黄色大香蕉| 不卡一级毛片| 久久精品久久久久久噜噜老黄 | 精品久久久久久久人妻蜜臀av| 免费av观看视频| 每晚都被弄得嗷嗷叫到高潮| 国产高潮美女av| 久久国产乱子免费精品| 少妇裸体淫交视频免费看高清| 中亚洲国语对白在线视频| 欧美中文日本在线观看视频| 9191精品国产免费久久| 国产精品久久久久久亚洲av鲁大| 天堂√8在线中文| 色av中文字幕| 亚州av有码| 国产亚洲精品综合一区在线观看| 欧美成人性av电影在线观看| 久久香蕉精品热| 国产成人aa在线观看| 国产三级在线视频| 日韩欧美在线二视频| 国产91精品成人一区二区三区| 高潮久久久久久久久久久不卡| 国产欧美日韩精品亚洲av| 午夜日韩欧美国产| 性插视频无遮挡在线免费观看| 国产一级毛片七仙女欲春2| 亚洲欧美激情综合另类| 久久99热这里只有精品18| 脱女人内裤的视频| 91在线精品国自产拍蜜月| 亚洲成人免费电影在线观看| 国产中年淑女户外野战色| 国产精品亚洲一级av第二区| 国产精品一及| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 色播亚洲综合网| netflix在线观看网站| ponron亚洲| 麻豆国产97在线/欧美| 国产精品久久久久久久久免 | 婷婷精品国产亚洲av| 偷拍熟女少妇极品色| 欧美极品一区二区三区四区| 身体一侧抽搐| 丰满人妻一区二区三区视频av| 九色成人免费人妻av| 美女cb高潮喷水在线观看| 窝窝影院91人妻| 久久久成人免费电影| 国产精品嫩草影院av在线观看 | 国产精品不卡视频一区二区 | 国产精品爽爽va在线观看网站| 久99久视频精品免费| 欧美高清成人免费视频www| 婷婷六月久久综合丁香| 欧美成狂野欧美在线观看| 久久人妻av系列| 亚洲色图av天堂| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 人人妻人人澡欧美一区二区| 偷拍熟女少妇极品色| 十八禁国产超污无遮挡网站| 欧美一区二区精品小视频在线| 国产精品自产拍在线观看55亚洲| 欧美丝袜亚洲另类 | 精品免费久久久久久久清纯| 嫩草影院新地址| 精品午夜福利视频在线观看一区| 国产精品亚洲美女久久久| 亚洲专区国产一区二区| 色吧在线观看| 午夜精品久久久久久毛片777| 欧美潮喷喷水| 中文字幕av成人在线电影| 国产精品av视频在线免费观看| 免费人成视频x8x8入口观看| 高清在线国产一区| 女同久久另类99精品国产91| 国产91精品成人一区二区三区| 男插女下体视频免费在线播放| 国产精品亚洲美女久久久| 18美女黄网站色大片免费观看| 国产精品影院久久| 观看免费一级毛片| 在线观看美女被高潮喷水网站 | 女人十人毛片免费观看3o分钟| 精品一区二区三区人妻视频| 欧美日韩亚洲国产一区二区在线观看| 变态另类丝袜制服| 少妇高潮的动态图| bbb黄色大片| 男女做爰动态图高潮gif福利片|