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

    月季RhMYB61基因的克隆及表達(dá)特性分析

    2017-11-04 07:12:08李紹翠姜新強(qiáng)丁愛琴劉慶超王奎玲劉慶華
    華北農(nóng)學(xué)報(bào) 2017年5期
    關(guān)鍵詞:切花基序級(jí)數(shù)

    李紹翠,姜新強(qiáng),丁愛琴,劉慶超,王奎玲,李 偉,劉慶華

    (青島農(nóng)業(yè)大學(xué) 園林與林學(xué)院,山東 青島 266109)

    月季RhMYB61基因的克隆及表達(dá)特性分析

    李紹翠,姜新強(qiáng),丁愛琴,劉慶超,王奎玲,李 偉,劉慶華

    (青島農(nóng)業(yè)大學(xué) 園林與林學(xué)院,山東 青島 266109)

    為了研究月季MYB類轉(zhuǎn)錄因子在月季花朵開放中的分子特征和表達(dá)特性,利用轉(zhuǎn)錄組測(cè)序獲得的序列信息,結(jié)合RACE技術(shù),從切花月季薩蔓莎中分離獲得一個(gè)MYB類型的轉(zhuǎn)錄因子基因,命名為RhMYB61(登錄號(hào)KY921844)。該基因ORF區(qū)包含1 323 bp,編碼441個(gè)aa,分子量為49.1 kDa,等電點(diǎn)為7.66,公式C2123H3280N620O688S18。系統(tǒng)進(jìn)化樹分析表明,RhMYB61與桃樹中的PpMYB86和枇杷EjMYB8聚為一類,屬于R2R3-MYB類型轉(zhuǎn)錄因子,進(jìn)一步的蛋白序列多重比較發(fā)現(xiàn),RhMYB61在N端具有保守的R2-MYB和R3-MYB區(qū)域,包含有7個(gè)保守的基序,在R3-MYB區(qū)域包含有核定位序列。利用實(shí)時(shí)定量PCR技術(shù)分析了RhMYB61的時(shí)空表達(dá)模式,結(jié)果顯示,隨著開花級(jí)數(shù)的增加,RhMYB61的表達(dá)特性呈逐漸增加的趨勢(shì),在開花級(jí)數(shù)5級(jí)時(shí)表達(dá)倍數(shù)最高;失水12 h處理提高了RhMYB61的表達(dá)倍數(shù);與對(duì)照相比較,乙烯處理6,12,24 h顯著提高了RhMYB61的表達(dá)特性,1-MCP則顯著抑制了其表達(dá)。上述結(jié)果為月季RhMYB61生物學(xué)功能的研究奠定了基礎(chǔ),同時(shí)也為今后月季的分子育種研究提供了理論參考。

    切花月季;RhMYB61;qRT-PCR;表達(dá)特性

    MYB轉(zhuǎn)錄因子是植物體內(nèi)最大的轉(zhuǎn)錄因子基因家族之一,成員眾多、功能多樣,廣泛參與了植物的次生代謝調(diào)控、激素和環(huán)境因子的應(yīng)答,并對(duì)細(xì)胞分化、細(xì)胞周期、器官的形成以及植物葉片的形態(tài)建成具有重要的調(diào)節(jié)作用[1]。眾多植物中均分離獲得了與非生物脅迫應(yīng)答相關(guān)的MYB類型轉(zhuǎn)錄因子,在草本植物擬南芥AtMYB32[2]、水稻[3]、菊花CmMYB1、CmMYB2[4],木本植物蘋果MdSIMYB1[5]、楊樹PtoMYB170[6]中的MYB轉(zhuǎn)錄因子均能夠?qū)Ψ巧锩{迫如高鹽、干旱、冷脅迫等做出積極的響應(yīng),MYB轉(zhuǎn)錄因子在對(duì)非生物脅迫中的應(yīng)答方面作用突出。

    月季(Rosahybrida)是我國(guó)傳統(tǒng)名花也是世界四大切花之一,其切花銷售量和銷售額在全球花卉貿(mào)易中均居切花之首[7]。因其花期長(zhǎng)、花色豐富等特點(diǎn),已逐漸成為花朵開放等方面研究的模式材料[8]。我國(guó)月季切花主產(chǎn)區(qū)為云南,主要消費(fèi)區(qū)域集中在北京、上海和廣州等城市,需要經(jīng)歷長(zhǎng)時(shí)間長(zhǎng)距離的干式運(yùn)輸,極易遭受失水脅迫和乙烯誘導(dǎo)的花朵開放,流通損耗在30%~50%,易造成僵花、僵蕾、花朵不能正常開放和衰老過速,制約了商品價(jià)值的發(fā)揮,對(duì)花卉生產(chǎn)者造成巨大的經(jīng)濟(jì)損失[9]。因此,系統(tǒng)研究月季切花失水脅迫耐性和乙烯作用調(diào)節(jié)的分子機(jī)理,對(duì)切花月季品質(zhì)提高具有重要的理論和實(shí)踐意義。

    迄今為止,在月季中獲得了3個(gè)MYB類型的轉(zhuǎn)錄因子,分別為RhMYB1[10]、RhMYBs4-1和RhMYBs6-1,它們?cè)诓煌ㄉ录净ㄇ嗨睾铣珊驼{(diào)控花瓣顏色方面發(fā)揮著調(diào)控作用[11],但對(duì)于數(shù)量眾多、功能多樣的月季MYBs還未充分挖掘。鑒于此,本試驗(yàn)選取切花月季薩曼莎為材料,結(jié)合前期轉(zhuǎn)錄組測(cè)序結(jié)果,從花瓣cDNA中克隆得到一個(gè)MYB基因,命名為RhMYB61。采用生物信息學(xué)的方法對(duì)該基因的氨基酸序列進(jìn)行序列比對(duì)、理化性質(zhì)和親/疏水性和二級(jí)結(jié)構(gòu)等方面的分析;并利用熒光定量PCR方法檢測(cè)了RhMYB61在月季不同開花級(jí)數(shù)、失水脅迫和乙烯處理下的轉(zhuǎn)錄本水平變化。為進(jìn)一步明確RhMYB61基因特性及其對(duì)非生物脅迫的響應(yīng)機(jī)制提供了理論基礎(chǔ)。

    1 材料和方法

    1.1試驗(yàn)材料

    供試材料為切花月季薩蔓莎(RosahybridaSamantha),開花級(jí)數(shù)2級(jí),按花枝長(zhǎng)30 cm、留3復(fù)葉的基準(zhǔn)剪切并插入蒸餾水中,復(fù)水1 h后備用。

    1.2試驗(yàn)處理

    開花級(jí)數(shù)參考馬男等[12]對(duì)月季開花級(jí)數(shù)的劃分并進(jìn)行取樣。

    失水脅迫處理參照按照Dai等[13]進(jìn)行,在溫度為22~25 ℃、相對(duì)濕度為30%~50%、140 μmol/(m2·s)熒光燈(照光12 h/d)照光的條件下水平干置0~24 h。在0,1,3,6,12,24 h分別取樣。

    乙烯處理參照馬男等[12]進(jìn)行,將修剪后的花枝置于去離子水中進(jìn)行處理。乙烯處理時(shí),在64 L密閉容器內(nèi)注入乙烯使其終濃度為10 μL/L,處理24 h,以密閉于不含乙烯氣體的64 L密封箱中的花枝為對(duì)照。

    各處理結(jié)束后將花瓣取下立即進(jìn)行液氮冷凍,存于-80 ℃冰箱備用。

    1.3試驗(yàn)方法

    1.3.1 總RNA提取和cDNA的合成 以1.2中處理后的切花月季花瓣為材料,Hot borate法提取不同處理下花瓣總RNA[14]。利用Nanodrop(美國(guó)Effendorf公司)檢測(cè)提取的RNA樣品濃度,1.2%瓊脂糖凝膠電泳檢測(cè)提取質(zhì)量。cDNA反轉(zhuǎn)錄采用PrimeScript RT reagent Kit With gDNA Erase(大連TaKaRa公司)試劑盒,將反轉(zhuǎn)錄產(chǎn)物稀釋4倍后-20 ℃保存?zhèn)溆谩?/p>

    1.3.2RhMYB61基因的克隆 利用月季EST數(shù)據(jù)庫(kù)[13]的Unigene序列為模板設(shè)計(jì)5′Race擴(kuò)增和3′Race擴(kuò)增特異性引物(表1)并進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系為上下游引物 (10 μmol/L) 各 1 μL,ExTaq聚合酶 (5 U/μL) 0.1 μL,10×PCR Buffer 2 μL,dNTP(2.5 mmol/L) 1.6 μL,總cDNA (1 μg/μL) 1 μL,ddH2O補(bǔ)足20 μL。PCR反應(yīng)程序?yàn)椋?4 ℃預(yù)變性5 min;94 ℃變性30 s,60 ℃退火30 s,72 ℃延伸1 min,30個(gè)循環(huán);72 ℃延伸5 min,4 ℃保存。將PCR產(chǎn)物進(jìn)行瓊脂糖凝膠電泳檢測(cè),利用TaKaRa Mini BESTAgarose Gel DNA Extraction Kit試劑盒(大連TaKaRa公司)回收目的條帶,并與pMD18-T載體連接。連接產(chǎn)物轉(zhuǎn)化大腸桿菌,藍(lán)白斑篩選,挑取白色單菌落鑒定,根據(jù)菌落PCR結(jié)果,將鑒定出的陽性克隆交由青島擎科梓熙生物工程有限公司測(cè)序。

    1.3.3 RhMYB61的生物信息學(xué)分析 采用ProtParam (http://web.Expasy.org/compute_pi/) 預(yù)測(cè)RhMYB61分子量和理論等電點(diǎn)。通過SignalP軟件進(jìn)行信號(hào)肽預(yù)測(cè)[15];利用ProtScale (http://web.Expasy.org/protscale/) 分析RhMYB61氨基酸序列的疏水性/親水性;利用TMHMM Server v. 2.0軟件(http://www.cbs.dtu.dk/services/TMHMM/)對(duì)RhMYB61進(jìn)行蛋白跨膜性分析[16]。利用ProtFun(http://www.cbs.dtu.dk/services/ProtFun/) 和PRABI(https://npsa-prabi.ibcp.fr/cgi-bin/secpred_sopma.pl)對(duì)蛋白質(zhì)功能和二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)[17]。利用Motif分析工具[18]進(jìn)行RhMYB61基序保守性分析。用DNAMAN軟件進(jìn)行序列比對(duì)及翻譯,利用MEGA 6.0軟件構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹[11]。

    表1 切花月季RhMYB61基因表達(dá)檢測(cè)引物Tab.1 Primers used to determine the expression of RhMYB61 in cut rose Samantha

    1.3.4RhMYB61基因表達(dá)特性分析 采用實(shí)時(shí)熒光定量PCR(Quantitative real-time PCR,qRT-PCR)分析月季切花RhMYB61基因在不同開花級(jí)數(shù)及失水、乙烯和1-MCP處理下的表達(dá)模式。不同處理?xiàng)l件總RNA提取和cDNA的合成如1.3.1,按照SYBR PremixExTaq試劑盒(大連TaKaRa公司)的操作說明進(jìn)行qRT-PCR,所有反應(yīng)均在Setp One Plus熒光定量PCR平臺(tái)(美國(guó)ABI公司)進(jìn)行,RhUBI1為內(nèi)參基因[19](表1)。反應(yīng)條件為:95 ℃預(yù)變性60 s;95 ℃變性10 s,58 ℃退火20 s,72 ℃延伸15 s,共40個(gè)循環(huán)。Melt曲線從72 ℃至99 ℃,第1步維持45 s,以后每升高1 ℃維持5 s。采用2-ΔΔCT法分析數(shù)據(jù)[20]。

    2 結(jié)果與分析

    2.1RhMYB61克隆

    以開花級(jí)數(shù)2級(jí)的切花月季薩蔓莎進(jìn)行失水12 h處理,提取總RNA并反轉(zhuǎn)錄成cDNA模板,根據(jù)轉(zhuǎn)錄組EST序列信息設(shè)計(jì)5′和3′擴(kuò)增引物(表1),分別進(jìn)行3′RACE和5′RACE,將上述片段與pMD-18T載體連接、轉(zhuǎn)化,并進(jìn)行藍(lán)白斑篩選,將鑒定正確的樣品測(cè)序。所得結(jié)果如圖1所示,月季RhMYB61登錄號(hào)為KY921844,全長(zhǎng)為1 727 bp,5′端非編碼區(qū)(5′-UTR)為178 bp,開放閱讀框?yàn)? 323 bp,編碼441個(gè)氨基酸,3′端非編碼區(qū)(3′-UTR)為226 bp。

    M. DL2000 Marker;A. 3′RACE 產(chǎn)物;B.5′RACE 產(chǎn)物;C. RhMYB61開放閱讀框片段克隆;D. RhMYB61全長(zhǎng)片段克隆。M. DL2000 Marker;A. 3′RACE fragment;B.5′RACE fragment;C. Open reading frame of RhMYB61;D. Full length of RhMYB61.

    2.2RhMYB61生物信息學(xué)分析

    2.2.1 RhMYB61基本特性分析 利用生物信息學(xué)軟件ProtParam等對(duì)RhMYB61的蛋白組成成分和理化性質(zhì)進(jìn)行了分析[21]。RhMYB61蛋白相對(duì)分子量49.1 kDa,理論等電點(diǎn)(pI)為7.66,公式C2123H3280N620O688S18,原子量6 729。RhMYB61蛋白為親水蛋白,其N-端頭部較親水,中間氨基酸也表現(xiàn)出親水性,蛋白在第81位疏水性最強(qiáng),最高值為2.189,第177位親水性最強(qiáng),最低值為-3.111。蛋白峰值的分布在0.5以上的比在-0.5以下的要少,表明蛋白為親水性蛋白。信號(hào)肽結(jié)果分析表明,RhMYB61不屬于分泌蛋白。蛋白跨膜性分析顯示RhMYB61氨基酸序列中未發(fā)現(xiàn)含有跨膜區(qū)域,不屬于膜蛋白;對(duì)RhMYB61蛋白進(jìn)行二級(jí)結(jié)構(gòu)預(yù)測(cè)顯示,其氨基酸組成α-螺旋(Hh)為137個(gè)(31.07%),延伸鏈(Ee)為53個(gè)(12.02%),無規(guī)則卷曲(Cc)為219個(gè)(49.66%),屬于不規(guī)則結(jié)構(gòu)。生物學(xué)功能預(yù)測(cè)顯示RhMYB61具有合成輔酶因子、嘌呤、嘧啶和脂肪酸代謝等方面的功能(表2)。

    表2 RhMYB61蛋白質(zhì)功能預(yù)測(cè)Tab.2 Biological function prediction of RhMYB61

    2.2.2 RhMYB61保守基序和氨基酸序列比對(duì) 選取擬南芥(Arabidopsisthaliana)、大豆(Glycinemax)、蘋果(Malusdomestica)、枇杷(Eriobotryajaponica)、白樺(Betulaplatyphylla)、菊花(Chrysanthemummorifolium)、厚葉蒴苣苔(Boeacrassifolia)、黑楊(Populusnigra)、桃(Prunuspercica)、水稻(Oryzasativa)和小麥(Triticumaestivum)中的MYB蛋白序列與月季RhMYB61進(jìn)行基序保守分析,結(jié)果如圖2所示。由圖可知,不同植物MYB蛋白保守基序數(shù)量和位置不盡相同,所有植物均含有基序1、基序2和基序3。RhMYB61與AtMYB61、PpMYB86、MdMYB、EjMYB8、PnMYB和GmMYB86氨基酸序列類似,包含有7個(gè)保守基序,在N-端含有5個(gè)保守基序,并呈串聯(lián)排列,C-端含有2個(gè)保守基序,分別為基序4和基序8。基序7僅在TaMYB3R1和OsMYB3R2中有,其他植物中未發(fā)現(xiàn)。

    方框表示MYB蛋白的保守基序;灰線條顯示非保守序列;標(biāo)尺代表100個(gè)氨基酸。Each motif is represented by a number in the box and the grey lines represent the nonconserved sequences;The scale bar represents 100 amino acids.

    進(jìn)一步選取擬南芥、大豆、蘋果、枇杷和白樺中的MYB蛋白與月季RhMYB61蛋白進(jìn)行序列多重比對(duì)分析(圖3)。結(jié)果顯示,5個(gè)不同物種MYB蛋白氨基酸序列與月季RhMYB61的一致性分別為42.63%,47.78%,63.19%,62.61%和44.62%,在12-60 aa和61-111 aa處含有2個(gè)保守的MYB結(jié)構(gòu)域,上述不同植物MYB蛋白序列都屬于2R類型MYB轉(zhuǎn)錄因子,其N-端都含有非常保守的R2-MYB和R3-MYB結(jié)構(gòu)域,具有 MYB家族典型的結(jié)構(gòu)特征,SeqNLS[22]分析顯示RhMYB61的R3-MYB結(jié)構(gòu)域中還包含有3處核定位序列(NLS,圖3方框部分),所有植物MYB蛋白C-端氨基酸序列高度分化,推測(cè)不同物種中MYB蛋白功能的多樣性。

    At.擬南芥;Gm.大豆;Md.蘋果;Ej.枇杷;Bp.白樺;Rh.月季。箭頭表示保守的R2和R3結(jié)構(gòu)域;NLS表示核定位序列。At. Arabidopsis thaliana;Gm.Glycine max;Md.Malus domestica;Ej.Eriobotrya japonica;Bp.Betula platyphylla;Rh.Rosa hybrid.The two conserved MYB domain repeat (R2 and R3) are indicated with arrows;Nuclear localization sites was shown with boxed.

    2.2.3 RhMYB61蛋白進(jìn)化樹分析 為進(jìn)一步探究月季RhMYB61與其他物種的同源性進(jìn)化關(guān)系,選用擬南芥、水稻、厚葉旋蒴苣苔、蘋果、菊花、大豆、枇杷、桃、黑楊、白樺和小麥中的22個(gè)MYB蛋白質(zhì)序列與月季RhMYB61構(gòu)建系統(tǒng)進(jìn)化關(guān)系圖(圖4)。結(jié)果顯示,月季RhMYB61與桃PpMYB86、枇杷EjMYB8、蘋果MdMYB、大豆GmMYB86、擬南芥AtMYB61親緣關(guān)系較近,與水稻OsMYB4、蘋果MdoMYB121、厚葉旋蒴苣苔BcMYB1親緣關(guān)系較遠(yuǎn),上述MYB蛋白均屬于R2R3-MYB類蛋白。單子葉植物小麥TaMYB3R1和水稻OsMYB3R1、OsMYB3R2屬于R1R2R3-MYB類蛋白,與月季RhMYB61親緣關(guān)系較遠(yuǎn)。上述結(jié)果暗示了MYB類蛋白在雙子葉植物和單子葉植物進(jìn)化上的保守性和多樣性。

    2.3RhMYB61表達(dá)特異性分析

    2.3.1RhMYB61在不同開花級(jí)數(shù)中的表達(dá) 為了明確RhMYB61在月季花朵發(fā)育中的作用,采用qRT-PCR的方法檢測(cè)RhMYB61在月季不同開花級(jí)數(shù)中的表達(dá)情況(圖5)。由圖5可知,隨著月季開花過程的進(jìn)行,花瓣中RhMYB61的表達(dá)量呈現(xiàn)逐漸增高的趨勢(shì),在開花級(jí)數(shù)0級(jí)和1級(jí)時(shí)相對(duì)表達(dá)量差別不大,分別為1.00和0.81,開花級(jí)數(shù)2級(jí)時(shí)表達(dá)量升高至3.00,顯著高于開花級(jí)數(shù)0級(jí)和1級(jí),開花級(jí)數(shù)5級(jí)時(shí)相對(duì)表達(dá)量最高,為13.60,顯著高于其他開花級(jí)數(shù)。RhMYB61在月季不同開花級(jí)數(shù)中的表達(dá)量暗示著其參與了花朵開放的過程。

    2.3.2 不同處理對(duì)RhMYB61表達(dá)量的影響 用定量PCR檢測(cè)了RhMYB61在轉(zhuǎn)錄本水平上對(duì)乙烯處理和失水脅迫處理的響應(yīng)(圖6)。結(jié)果表明,在正常瓶插條件下,花瓣RhMYB61瓶插1,6,12 h只有微弱的表達(dá),表達(dá)量分別為0.73,1.50,1.80,瓶插24 h后表達(dá)量升高到5.69,乙烯處理明顯誘導(dǎo)了RhMYB61的表達(dá),乙烯處理6,12 h后,相對(duì)表達(dá)量分別為7.92,7.72,高于同一時(shí)間正常瓶插下的表達(dá)量,乙烯處理24 hRhMYB61相對(duì)表達(dá)量最高為8.37。乙烯抑制劑1-MCP處理降低了RhMYB61的表達(dá)水平,1-MCP處理1,6,12,24 h后RhMYB61表達(dá)量分別為0.21,0.46,0.57,1.07。

    At.擬南芥;Gm.大豆;Md.蘋果;Ej.枇杷;Bp.白樺;Cm.菊花;Pp.桃;Pn.黑楊;Bc.厚葉旋蒴苣苔;Ta.小麥;Os.水稻。At.Arabidopsis thaliana;Gm.Glycine max;Md.Malus domestica;Ej.Eriobotrya japonica;Bp.Betula platyphylla;Cm.Chrysanthemum morifolium;Pp.Prunus percica;Pn.Populus nigra;Bc.Boea crassifolia;Ta.Triticum aestivum;Os.Oryza sativa.

    圖中不同小寫字母表示處理間差異顯著(P<0.05)。圖6-7同。Bars with different lowercase letters indicate significant differences by Turkey HSD test(P<0.05).The same as Fig.6-7.

    對(duì)RhMYB61在不同失水情況下的表達(dá)特性進(jìn)行了分析,結(jié)果如圖7所示,隨著失水脅迫時(shí)間的延長(zhǎng),RhMYB61的表達(dá)呈現(xiàn)出先上升后下降的趨勢(shì),與正常瓶插植株相比,失水脅迫1 h的RhMYB61轉(zhuǎn)錄本水平差異不大,失水脅迫3,6 h誘導(dǎo)了RhMYB61表達(dá)水平的提高,相對(duì)表達(dá)量分別為1.56和1.87,隨后RhMYB61表達(dá)量逐漸降低,失水脅迫24 h表達(dá)量最低,為0.48。

    圖6 不同時(shí)間乙烯處理RhMYB61的表達(dá)Fig.6 Expression analysis of RhMYB61 in response to ethylene in rose petals

    圖7 RhMYB61在切花月季不同失水處理的表達(dá)Fig.7 Expression analysis of RhMYB61 expression under different dehydration time

    上述結(jié)果表明,乙烯處理顯著誘導(dǎo)了RhMYB61的表達(dá),失水脅迫處理對(duì)RhMYB61的表達(dá)影響不大。

    3 討論

    已有研究表明,MYB是植物體內(nèi)最大的一類轉(zhuǎn)錄因子家族之一,在植物次生代謝調(diào)控、激素和環(huán)境因子的應(yīng)答、細(xì)胞分化、細(xì)胞周期、器官的形成以及植物葉片的形態(tài)建成中發(fā)揮了重要作用。對(duì)多種植物中的MYB轉(zhuǎn)錄因子的功能研究發(fā)現(xiàn),R2R3類型MYB轉(zhuǎn)錄因子在植物響應(yīng)非生物脅迫逆境方面發(fā)揮了非常關(guān)鍵的作用,擬南芥AtMYB61[23]、水稻OsMYB91[24]、小麥TaMYB3R1[25]和大豆GmMYBJ7[26]中報(bào)道的R2R3類型MYB轉(zhuǎn)錄因子均能夠提高植株對(duì)干旱、高鹽、低溫等非生物逆境脅迫的抵抗能力。而在切花月季中,未見有對(duì)花瓣MYB轉(zhuǎn)錄因子參與非生物脅迫逆境的研究報(bào)道。本研究中,通過高通量測(cè)序和RACE技術(shù)克隆了月季切花中R2R3類型MYB轉(zhuǎn)錄因子RhMYB61,并利用定量PCR檢測(cè)了其在不同開花級(jí)數(shù)、乙烯處理和失水脅迫過程中的表達(dá)變化。結(jié)果表明,RhMYB61廣泛參與了切花月季花朵開放的過程,隨著開花級(jí)數(shù)的逐漸增加,表達(dá)量逐漸升高,開花級(jí)數(shù)5級(jí)時(shí)表達(dá)量最高,暗示著RhMYB61在月季花朵開放和花瓣擴(kuò)展方面發(fā)揮著重要的作用,特別是花朵開放后期,參與了月季花朵的衰老進(jìn)程。對(duì)RhMYB61的生物信息學(xué)分析表明,其生物學(xué)功能多集中在參與合成輔酶因子、嘌呤、嘧啶和脂肪酸代謝等方面,其參與月季花瓣衰老進(jìn)程是否通過上述過程進(jìn)行還有待進(jìn)一步研究。

    切花月季在采收后,從花農(nóng)生產(chǎn)者到消費(fèi)者手中要經(jīng)歷長(zhǎng)時(shí)間的干式遠(yuǎn)距離運(yùn)輸,運(yùn)輸過程中的外源乙烯積累和失水脅迫極易導(dǎo)致切花開放異常,對(duì)切花品質(zhì)和商品價(jià)值產(chǎn)生影響,是月季切花采后損耗的重要因素,因此了解不同切花月季品種對(duì)外源乙烯和失水脅迫的響應(yīng)機(jī)制,對(duì)于切花月季采后保鮮、提高流通質(zhì)量、改良現(xiàn)代月季育種技術(shù)等均具有重要的現(xiàn)實(shí)意義。切花月季薩曼莎對(duì)外源乙烯非常敏感,外源乙烯處理能夠促進(jìn)月季花朵的開放,乙烯抑制劑1-MCP能夠通過與乙烯受體的結(jié)合,阻斷乙烯的感受,進(jìn)而抑制花瓣內(nèi)源乙烯的自我催化,延緩花瓣的衰老啟動(dòng),有效抑制花朵乙烯的生成。本試驗(yàn)中,乙烯處理顯著增加了RhMYB61的轉(zhuǎn)錄本水平表達(dá),乙烯抑制劑1-MCP處理則抑制了RhMYB61的表達(dá)水平,暗示著RhMYB61參與了乙烯誘導(dǎo)的月季花朵開放過程。研究表明,月季乙烯生物合成關(guān)鍵基因(RhACSs和RhACO1)和乙烯信號(hào)轉(zhuǎn)導(dǎo)關(guān)鍵基因(RhETRs和RhCTRs)在月季花朵開放過程中發(fā)揮著關(guān)鍵作用[27-28],RhMYB61是否通過調(diào)控上述基因或者通過調(diào)控其他功能蛋白進(jìn)而參與對(duì)花瓣擴(kuò)展的調(diào)控,明確RhMYB61在月季花朵開放中的調(diào)節(jié)機(jī)理將是后續(xù)研究中需要解決的重要問題。

    前人研究表明,失水脅迫對(duì)乙烯具有誘導(dǎo)作用[29]。切花月季薩曼莎經(jīng)過乙烯處理和失水脅迫處理后花朵開放狀態(tài)非常類似,切花月季中可能存在水分脅迫和乙烯信號(hào)的交叉[30]?;诒驹囼?yàn)結(jié)果,失水脅迫處理3,6 h提高了RhMYB61的轉(zhuǎn)錄本水平,而失水脅迫處理12,24 h則降低了其表達(dá),說明RhMYB61至少在轉(zhuǎn)錄水平上未參與月季切花失水脅迫與乙烯信號(hào)的交叉。本研究克隆了月季RhMYB61基因,對(duì)其進(jìn)行了初步研究,后期可進(jìn)一步對(duì)RhMYB61的功能和在切花月季響應(yīng)失水脅迫影響花瓣擴(kuò)展中的作用機(jī)制進(jìn)行研究,對(duì)于明確花瓣開放機(jī)理和參與非生物脅迫耐性,利用現(xiàn)代生物育種技術(shù)改良現(xiàn)代切花月季品質(zhì)將具有一定的意義。

    [1] 王 棟,李利紅,陳志玲,等. 擬南芥根特異表達(dá)轉(zhuǎn)錄因子AtMYB305的鑒定及功能研究[J]. 科學(xué)通報(bào),2001,46(21):1804-1809.

    [2] 馮盼盼.AtMYB32對(duì)AtGA20ox1基因表達(dá)的調(diào)控及在植物抗脅迫中的作用研究[D]. 長(zhǎng)沙:湖南大學(xué),2016.

    [3] Smita S,Katiyar A,Chinnusamy V,et al. Transcriptional regulatory network analysis of MYB transcription factor family genes in rice[J]. Frontiers in Plant Science,2015,6(682):1157.

    [4] 單 紅. 菊花CmMYB1和CmMYB2轉(zhuǎn)錄因子基因的克隆及功能鑒定[D]. 南京:南京農(nóng)業(yè)大學(xué),2011.

    [5] Wang R K,Cao Z H,Hao Y J. Overexpression of a R2R3 MYB geneMdSIMYB1 increases tolerance to multiple stresses in transgenic tobacco and apples[J]. Physiologia Plantarum,2014,150(1):76-87.

    [6] 劉 瑞. 楊樹PtoMYB170轉(zhuǎn)錄因子在次生壁合成調(diào)控及抗旱中的功能研究[D]. 重慶:西南大學(xué),2016.

    [7] Zuker A,Tzfira T,Vainstein A. Genetic engineering for cut-flower improvement[J]. Biotechnology Advances,1998,16(1):33-79.

    [8] Debener T,Linde M. Exploring complex ornamental genomes:the rose as a model plant[J]. Critical Reviews in Plant Sciences,2009,28(4):267-280.

    [9] 呂培濤. HD-Zip I家族成員RhHB1在月季花瓣衰老過程中的功能分析[D]. 北京:中國(guó)農(nóng)業(yè)大學(xué),2014.

    [10] Yan H,Zhang H,Wang Q,et al. Isolation and identification of a putative scent-related geneRhMYB1 from rose[J]. Molecular Biology Reports,2011,38(7):4475-4482.

    [11] 趙 佳,劉 榮,楊 帆,等. 月季花青素苷相關(guān)R2R3-MYB 蛋白基因的克隆和表達(dá)分析[J]. 中國(guó)農(nóng)業(yè)科學(xué),2015,48(7):1392-1404.

    [12] 馬 男,蔡 蕾,陸旺金,等. 外源乙烯對(duì)月季(Rosahybrida)切花花朵開放的影響與乙烯生物合成相關(guān)基因表達(dá)的關(guān)聯(lián)[J]. 中國(guó)科學(xué)C輯,2005,35(2):104-114.

    [13] Dai F W,Zhang C Q,Jiang X Q,et al.RhNAC2 andRhEXPA4 are involved in the regulation of dehydration tolerance during the expansion of rose petals[J]. Plant Physiology,2012,160(4):2064-2082.

    [14] Wan C Y,Wilkins T A. A modified hot borate method significantly enhances the yield of high-quality RNA from cotton (GossypiumhirsutumL.) [J]. Analytical Biochemistry,1994,223(1):7-12.

    [15] Petersen T N,Brunak S,Von Heijne G,et al. SignalP 4.0:discriminating signal peptides from transmembrane regions[J]. Nature Methods,2011,8(10):785-786.

    [16] 徐 麗,陳 新,宗曉娟,等. 甜櫻桃砧木PcDHN1的克隆及其對(duì)非生物脅迫的響應(yīng)[J]. 核農(nóng)學(xué)報(bào),2017,31(1):14-20.

    [17] 董慧霞,理永霞,賈秀貞,等. 楊樹SABP2基因的克隆與分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào):自然科學(xué)版,2017,45(3):82-88,95.

    [18] Ma W,Noble W S,Bailey T L. Motif-based analysis of large nucleotide data sets using MEME-ChIP[J]. Nature Protocols,2014,9(6):1428-1450.

    [19] Jiang X Q,Zhang C Q,Lu P T,et al.RhNAC3,a stress-associated NAC transcription factor,has a role in dehydration tolerance through regulating osmotic stress-related genes in rose petals[J]. Plant Biotechnology Journal,2014,12(1):38-48.

    [20] Livak K J,Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method[J]. Methods,2001,25(4):402-408.

    [21] 鄭伶杰,馬 宏,張 媛,等. 蘋果矮化砧木GA20氧化酶基因的克隆及其表達(dá)分析[J]. 華北農(nóng)學(xué)報(bào),2017,32(1):53-59.

    [22] Lin J R,Hu J. SeqNLS:nuclear localization signal prediction based on frequent pattern mining and linear motif scoring[J]. PLoS One,2013,8(10):e76864.

    [23] 王 云,任永兵,楊 碩,等. 擬南芥AtMYB61基因的克隆及表達(dá)載體構(gòu)建[J]. 合肥工業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版,2011,34(5):757-761.

    [24] 朱 寧. MYB類轉(zhuǎn)錄因子OsMYB91在水稻發(fā)育和抗逆中的功能研究[D]. 武漢:華中農(nóng)業(yè)大學(xué),2014.

    [25] Cai H,Tian S,Dong H,et al. Pleiotropic effects ofTaMYB3R1 on plant development and response to osmotic stress in transgenicArabidopsis[J]. Gene,2015,558(2):227-234.

    [26] 楊文杰,吳燕民,唐益雄. 大豆轉(zhuǎn)錄因子基因GmMYBJ7的表達(dá)及功能分析[J]. 華北農(nóng)學(xué)報(bào),2012,27(6):24-29.

    [27] Ma N,Cai L,Lu W,et al. Exogenous ethylene influences flower opening of cut roses (Rosahybrida) by regulating the genes encoding ethylene biosynthesis enzymes[J]. Science in China Series C-life Sciences,2005,48(5):434-444.

    [28] 譚 輝. 外源乙烯通過調(diào)節(jié)乙烯受體和信號(hào)轉(zhuǎn)導(dǎo)元件基因表達(dá)影響月季切花開放[D]. 北京:中國(guó)農(nóng)業(yè)大學(xué),2005.

    [29] Nakano R,Ogura E,Kubo Y,et al. Ethylene biosynthesis in detached young persimmon fruit is initiated in calyx and modulated by water loss from the fruit[J]. Plant Physiology,2003,131(1):276-286.

    [30] 王 婷,仝 征,馬 男,等. 切花月季Rh-DREB1基因的克隆及其在乙烯和失水脅迫下的表達(dá)[J]. 園藝學(xué)報(bào),2009,36(1):65-72.

    CloningandExpressionCharacteristicsofRhMYB61GenefromRosahybrida

    LI Shaocui,JIANG Xinqiang,DING Aiqin,LIU Qingchao,WANG Kuiling,LI Wei,LIU Qinghua

    (College of Landscape Architecture and Forestry,Qingdao Agricultural University,Qingdao 266109,China)

    In order to determine the molecular characterization and expression patterns of MYB transcription factors involved in flower opening in Rose. According to the expressed sequence tags from transcriptome assembly and RACE PCR methods,we obtained a new MYB transcription factor gene,and namedRhMYB61(GenBank accession number is KY921844 ),inRosahybridaSamantha petals. The ORF ofRhMYB61 was 1 323 bp,which encoded 441 amino acids. Bioinformatics analysis showed that the molecular weight and theoretical isoelectric point of RhMYB61 was 49.1 kDa and 7.66,respectively,and the formula of RhMYB61 was C2123H3280N620O688S18. Phylogenic tree analysis showed that RhMYB61 was clustered with PpMYB86 and EjMYB8,and belonged to the R2R3-MYB type transcription factor. Multi-alignment of RhMYB61 of several plants found that RhMYB61 had the seven characteristic conservative motifs.Two conserved domains,including R2-MYB domain and R3-MYB domain,were found in its N-terminal. The nuclear localization sites was found in R3-MYB domain. Qualitative PCR analysis showed thatRhMYB61 expression was induced as the rose flower opening.RhMYB61 expression was accumulated at 12 h dehydration treatment,compared with control. In addition,RhMYB61 transcript was rapidly and dramatically induced by 6,12,24 h ethylene treatment compared with control,and was inhibited following 1-methylcyclopropene (1-MCP) treatment.These results above gave a clear data for understanding the biological function ofRhMYB61, and provided a foundation for further research of molecular breeding of rose plant.

    Cut roses;RhMYB61;qRT-PCR;Expression characteristics

    2017-07-21

    國(guó)家自然科學(xué)基金項(xiàng)目(31501798);山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金項(xiàng)目(2014BSB01563)

    李紹翠(1992-),女,山東招遠(yuǎn)人,在讀碩士,主要從事園林植物種質(zhì)資源創(chuàng)新和利用研究。

    劉慶華(1962-),男,山東乳山人,教授,博士,博士生導(dǎo)師,主要從事園林植物種質(zhì)資源創(chuàng)新及園林植物種植設(shè)計(jì)研究。

    Q78;S685.03

    A

    1000-7091(2017)05-0061-08

    10.7668/hbnxb.2017.05.010

    猜你喜歡
    切花基序級(jí)數(shù)
    納米銀和1-MCP處理對(duì)月季切花的保鮮作用
    EPIYA 基序與幽門螺桿菌感染相關(guān)胃病關(guān)系的研究進(jìn)展
    帶TRS基序突變的新型冠狀病毒威脅更大
    芥藍(lán)Aux/IAA家族基因生物信息學(xué)與表達(dá)分析
    Dirichlet級(jí)數(shù)及其Dirichlet-Hadamard乘積的增長(zhǎng)性
    切花生產(chǎn)技術(shù)
    幾個(gè)常數(shù)項(xiàng)級(jí)數(shù)的和
    遼西劍蘭切花生產(chǎn)技術(shù)要點(diǎn)
    p級(jí)數(shù)求和的兩種方法
    Dirichlet級(jí)數(shù)的Dirichlet-Hadamard乘積
    日韩av在线大香蕉| 美女免费视频网站| 19禁男女啪啪无遮挡网站| 婷婷丁香在线五月| 国产精品av视频在线免费观看| 免费大片18禁| 蜜桃久久精品国产亚洲av| 国产真实乱freesex| 免费观看人在逋| 国产精品九九99| 亚洲精品色激情综合| 一区二区三区高清视频在线| 亚洲av五月六月丁香网| 制服人妻中文乱码| 午夜激情福利司机影院| 日韩有码中文字幕| 亚洲av第一区精品v没综合| 国产三级黄色录像| 国产伦精品一区二区三区视频9 | 国产精华一区二区三区| 嫩草影院精品99| 日韩欧美三级三区| 日韩 欧美 亚洲 中文字幕| 国产午夜福利久久久久久| 国产乱人伦免费视频| 一本精品99久久精品77| 国产日本99.免费观看| 国产视频内射| 欧美激情在线99| 久久久成人免费电影| 久久精品aⅴ一区二区三区四区| 亚洲一区二区三区色噜噜| 欧美最黄视频在线播放免费| 老司机午夜十八禁免费视频| 中亚洲国语对白在线视频| АⅤ资源中文在线天堂| 最新中文字幕久久久久 | 国内揄拍国产精品人妻在线| 精品久久久久久,| 国产不卡一卡二| 日本与韩国留学比较| 老熟妇仑乱视频hdxx| 国产乱人视频| 亚洲欧美日韩无卡精品| 午夜免费观看网址| 国产精品一及| 女人高潮潮喷娇喘18禁视频| 欧美+亚洲+日韩+国产| 国产伦精品一区二区三区视频9 | 中文资源天堂在线| 日日夜夜操网爽| 久久久精品大字幕| 免费一级毛片在线播放高清视频| 精品一区二区三区av网在线观看| 两个人看的免费小视频| 精品久久久久久久久久免费视频| 窝窝影院91人妻| 美女被艹到高潮喷水动态| 成人欧美大片| 精品国内亚洲2022精品成人| 亚洲欧洲精品一区二区精品久久久| www.精华液| 亚洲男人的天堂狠狠| 久久久久性生活片| 黄色片一级片一级黄色片| 午夜精品久久久久久毛片777| 丝袜人妻中文字幕| 天天躁日日操中文字幕| 精品国产乱码久久久久久男人| 18禁美女被吸乳视频| 欧美一区二区国产精品久久精品| 国产高潮美女av| 99视频精品全部免费 在线 | 成年女人永久免费观看视频| 亚洲精品色激情综合| 国产精品久久久人人做人人爽| 热99在线观看视频| 国产精品av视频在线免费观看| 男人舔奶头视频| 1000部很黄的大片| 婷婷亚洲欧美| 香蕉久久夜色| 久久久久久久久中文| 国产成人av激情在线播放| 两个人看的免费小视频| 亚洲人成网站高清观看| 国产高清videossex| 中出人妻视频一区二区| e午夜精品久久久久久久| 久久这里只有精品19| 久久久久国内视频| 亚洲黑人精品在线| 国产一区二区三区在线臀色熟女| www.999成人在线观看| 国产精品一及| www日本黄色视频网| 美女大奶头视频| 色在线成人网| 亚洲精品在线观看二区| 国产又色又爽无遮挡免费看| 日韩欧美精品v在线| 五月伊人婷婷丁香| 国产野战对白在线观看| 波多野结衣高清无吗| 亚洲 国产 在线| 日本免费一区二区三区高清不卡| 网址你懂的国产日韩在线| 亚洲av五月六月丁香网| 欧美绝顶高潮抽搐喷水| 亚洲国产日韩欧美精品在线观看 | 国产v大片淫在线免费观看| 午夜久久久久精精品| 在线播放国产精品三级| 看黄色毛片网站| 欧美绝顶高潮抽搐喷水| 九九在线视频观看精品| 男人舔女人的私密视频| 少妇人妻一区二区三区视频| www.999成人在线观看| 99国产精品一区二区蜜桃av| 男人的好看免费观看在线视频| 麻豆成人午夜福利视频| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 不卡av一区二区三区| 亚洲国产高清在线一区二区三| 99国产精品一区二区三区| 成人一区二区视频在线观看| 人妻久久中文字幕网| 成年免费大片在线观看| 久久久色成人| 麻豆av在线久日| 午夜福利视频1000在线观看| 丝袜人妻中文字幕| 日本 欧美在线| 91久久精品国产一区二区成人 | 国产野战对白在线观看| 少妇人妻一区二区三区视频| 国产精品久久久久久精品电影| 亚洲九九香蕉| 国产成人影院久久av| 韩国av一区二区三区四区| 中文字幕av在线有码专区| 看黄色毛片网站| 日本熟妇午夜| 成人av一区二区三区在线看| 精品久久久久久,| 久久精品国产清高在天天线| 人人妻,人人澡人人爽秒播| 一进一出抽搐gif免费好疼| 999久久久国产精品视频| 在线视频色国产色| 久久精品aⅴ一区二区三区四区| av欧美777| 91av网一区二区| 在线视频色国产色| av欧美777| www日本黄色视频网| 18禁黄网站禁片免费观看直播| 他把我摸到了高潮在线观看| 国产av不卡久久| 国产精品一区二区三区四区久久| 久久久久久久午夜电影| www日本黄色视频网| 日韩成人在线观看一区二区三区| 91av网一区二区| 亚洲人成网站在线播放欧美日韩| 欧美日韩国产亚洲二区| 欧美另类亚洲清纯唯美| www.熟女人妻精品国产| 脱女人内裤的视频| 国内精品一区二区在线观看| 不卡av一区二区三区| 又紧又爽又黄一区二区| 黄色 视频免费看| 神马国产精品三级电影在线观看| 巨乳人妻的诱惑在线观看| 久久精品国产亚洲av香蕉五月| 日本精品一区二区三区蜜桃| 俄罗斯特黄特色一大片| 色播亚洲综合网| 男人和女人高潮做爰伦理| 亚洲成a人片在线一区二区| 啦啦啦韩国在线观看视频| 在线十欧美十亚洲十日本专区| 在线观看66精品国产| 午夜视频精品福利| 国产精品久久久久久精品电影| 久久久国产精品麻豆| 亚洲精华国产精华精| 午夜激情福利司机影院| 婷婷丁香在线五月| 99热精品在线国产| 国产黄片美女视频| 欧美日本亚洲视频在线播放| 国产亚洲欧美在线一区二区| 最近最新中文字幕大全免费视频| 久久这里只有精品19| 亚洲天堂国产精品一区在线| 超碰成人久久| 欧美一区二区精品小视频在线| 国产一区在线观看成人免费| 九九热线精品视视频播放| av视频在线观看入口| 国产成人av教育| 亚洲av成人精品一区久久| 亚洲欧美日韩东京热| 婷婷精品国产亚洲av| 一区二区三区高清视频在线| av欧美777| 国产三级黄色录像| 成人午夜高清在线视频| 亚洲成人久久爱视频| 99久久成人亚洲精品观看| 国产av麻豆久久久久久久| av天堂中文字幕网| 午夜影院日韩av| 一二三四社区在线视频社区8| 久久中文字幕人妻熟女| 色播亚洲综合网| 老司机午夜福利在线观看视频| 丁香六月欧美| 变态另类成人亚洲欧美熟女| 好男人电影高清在线观看| 在线永久观看黄色视频| 久久香蕉国产精品| 狠狠狠狠99中文字幕| 丝袜人妻中文字幕| 夜夜夜夜夜久久久久| 国产精品久久电影中文字幕| 精品久久久久久久人妻蜜臀av| 国产精品99久久99久久久不卡| 无遮挡黄片免费观看| 五月玫瑰六月丁香| 99riav亚洲国产免费| svipshipincom国产片| 亚洲无线观看免费| 国产精品一区二区三区四区免费观看 | 香蕉国产在线看| 亚洲一区高清亚洲精品| 两个人看的免费小视频| 757午夜福利合集在线观看| 欧美日韩黄片免| 在线视频色国产色| 最近视频中文字幕2019在线8| 久久精品国产综合久久久| 国产又黄又爽又无遮挡在线| 免费观看精品视频网站| 欧美激情久久久久久爽电影| 欧美一区二区精品小视频在线| 免费电影在线观看免费观看| 91av网站免费观看| 18美女黄网站色大片免费观看| 精品熟女少妇八av免费久了| 国内久久婷婷六月综合欲色啪| 国产激情偷乱视频一区二区| 香蕉久久夜色| 三级男女做爰猛烈吃奶摸视频| 亚洲成人久久性| 久久久久久人人人人人| 国产精品久久久av美女十八| 男女做爰动态图高潮gif福利片| 免费高清视频大片| а√天堂www在线а√下载| 国产欧美日韩一区二区精品| 欧美一级毛片孕妇| 亚洲 欧美一区二区三区| 欧美乱色亚洲激情| 国产极品精品免费视频能看的| 亚洲美女黄片视频| 搡老岳熟女国产| 国产精品免费一区二区三区在线| 国产亚洲欧美在线一区二区| 天天躁日日操中文字幕| 在线免费观看不下载黄p国产 | 一级a爱片免费观看的视频| 国产一区二区激情短视频| 欧美zozozo另类| 别揉我奶头~嗯~啊~动态视频| 免费看日本二区| 级片在线观看| 岛国在线观看网站| 亚洲精品在线观看二区| 精品乱码久久久久久99久播| 成人午夜高清在线视频| 亚洲成人久久爱视频| 高清在线国产一区| 香蕉久久夜色| 国产三级中文精品| 99热只有精品国产| 免费人成视频x8x8入口观看| 国产真实乱freesex| 一卡2卡三卡四卡精品乱码亚洲| 小蜜桃在线观看免费完整版高清| 麻豆久久精品国产亚洲av| 亚洲成人免费电影在线观看| 最新中文字幕久久久久 | 亚洲七黄色美女视频| 国产日本99.免费观看| 日日摸夜夜添夜夜添小说| 每晚都被弄得嗷嗷叫到高潮| 国产精品一及| 免费人成视频x8x8入口观看| 免费看日本二区| 亚洲专区字幕在线| 久久人人精品亚洲av| 欧美性猛交╳xxx乱大交人| x7x7x7水蜜桃| 久久精品国产综合久久久| 91av网一区二区| 91久久精品国产一区二区成人 | 十八禁人妻一区二区| 精品一区二区三区av网在线观看| 亚洲国产高清在线一区二区三| 欧美黄色片欧美黄色片| 国模一区二区三区四区视频 | 美女大奶头视频| 亚洲九九香蕉| 国产视频内射| 国产伦一二天堂av在线观看| 美女免费视频网站| 日韩av在线大香蕉| 中文在线观看免费www的网站| 亚洲中文字幕一区二区三区有码在线看 | 精品国产乱码久久久久久男人| 99热6这里只有精品| 9191精品国产免费久久| 男人舔奶头视频| 午夜福利免费观看在线| 国产成人aa在线观看| 欧美日韩黄片免| 国产又色又爽无遮挡免费看| 亚洲黑人精品在线| 国产成人av激情在线播放| 亚洲精品在线观看二区| 免费高清视频大片| 国产97色在线日韩免费| 看黄色毛片网站| 亚洲自拍偷在线| 免费人成视频x8x8入口观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区三区四区免费观看 | 老汉色av国产亚洲站长工具| tocl精华| 啦啦啦韩国在线观看视频| 日韩成人在线观看一区二区三区| av女优亚洲男人天堂 | 18美女黄网站色大片免费观看| netflix在线观看网站| 91字幕亚洲| 亚洲国产精品成人综合色| 中文字幕久久专区| 亚洲专区中文字幕在线| 亚洲国产精品成人综合色| 91九色精品人成在线观看| 长腿黑丝高跟| 丰满人妻熟妇乱又伦精品不卡| 脱女人内裤的视频| 91字幕亚洲| 18美女黄网站色大片免费观看| 国产精品亚洲美女久久久| 精品久久久久久,| 99在线视频只有这里精品首页| 黄色日韩在线| 国产精品一区二区三区四区久久| 欧美日韩一级在线毛片| 国产精品国产高清国产av| 亚洲熟妇熟女久久| aaaaa片日本免费| 国产黄片美女视频| 极品教师在线免费播放| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 国产美女午夜福利| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 蜜桃久久精品国产亚洲av| 日韩人妻高清精品专区| 欧美性猛交黑人性爽| 欧美一级毛片孕妇| 久久国产精品影院| 国产精品女同一区二区软件 | 久久热在线av| 男人舔奶头视频| 日本免费a在线| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产| 制服丝袜大香蕉在线| 久久人妻av系列| 久久久久久九九精品二区国产| 免费av毛片视频| 久久天躁狠狠躁夜夜2o2o| 高潮久久久久久久久久久不卡| 99国产极品粉嫩在线观看| 亚洲性夜色夜夜综合| 男人的好看免费观看在线视频| 真人做人爱边吃奶动态| 精品国产超薄肉色丝袜足j| 99国产极品粉嫩在线观看| 精品国产乱码久久久久久男人| 99热这里只有精品一区 | 一a级毛片在线观看| 国产成人精品久久二区二区91| 法律面前人人平等表现在哪些方面| 成人特级黄色片久久久久久久| 天堂动漫精品| 色av中文字幕| 一区二区三区高清视频在线| 免费看美女性在线毛片视频| 亚洲在线自拍视频| 国产精品影院久久| av视频在线观看入口| 成人18禁在线播放| 精品国产超薄肉色丝袜足j| 国产精品日韩av在线免费观看| 精品日产1卡2卡| 在线观看66精品国产| 女生性感内裤真人,穿戴方法视频| 99久久国产精品久久久| 亚洲国产欧美网| 欧美av亚洲av综合av国产av| 欧美xxxx黑人xx丫x性爽| 高清在线国产一区| 麻豆久久精品国产亚洲av| a级毛片a级免费在线| 亚洲精品456在线播放app | 真人做人爱边吃奶动态| 最好的美女福利视频网| 国产高清有码在线观看视频| 在线播放国产精品三级| 亚洲 欧美 日韩 在线 免费| 视频区欧美日本亚洲| 欧美中文日本在线观看视频| 国产黄片美女视频| 男女床上黄色一级片免费看| 日韩有码中文字幕| 日韩国内少妇激情av| 久久天堂一区二区三区四区| 精品一区二区三区视频在线 | 久久久久久久久免费视频了| 国产精品一区二区三区四区免费观看 | 成在线人永久免费视频| 欧美日韩亚洲国产一区二区在线观看| 午夜福利18| 热99在线观看视频| 后天国语完整版免费观看| 欧美乱码精品一区二区三区| 午夜a级毛片| 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 国产三级黄色录像| 亚洲人成网站高清观看| 99国产极品粉嫩在线观看| 香蕉国产在线看| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 又大又爽又粗| 国产成+人综合+亚洲专区| 亚洲美女视频黄频| 非洲黑人性xxxx精品又粗又长| 亚洲精品一卡2卡三卡4卡5卡| АⅤ资源中文在线天堂| 老司机午夜福利在线观看视频| 亚洲电影在线观看av| а√天堂www在线а√下载| 久久久精品大字幕| 色在线成人网| 国产精华一区二区三区| 在线观看舔阴道视频| 国产成年人精品一区二区| 狂野欧美激情性xxxx| 亚洲专区中文字幕在线| 18禁裸乳无遮挡免费网站照片| 特级一级黄色大片| 十八禁网站免费在线| 色综合欧美亚洲国产小说| 12—13女人毛片做爰片一| 一边摸一边抽搐一进一小说| 小蜜桃在线观看免费完整版高清| 久久99热这里只有精品18| 免费在线观看成人毛片| 国产爱豆传媒在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲av电影在线进入| 欧美大码av| 日韩欧美精品v在线| 欧美色视频一区免费| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 亚洲精品乱码久久久v下载方式 | 国产精品亚洲美女久久久| 美女高潮喷水抽搐中文字幕| 欧美色欧美亚洲另类二区| 欧美三级亚洲精品| 一个人免费在线观看电影 | 他把我摸到了高潮在线观看| 亚洲人成伊人成综合网2020| 亚洲专区国产一区二区| 亚洲性夜色夜夜综合| 制服人妻中文乱码| 91av网一区二区| 一级作爱视频免费观看| 国产免费男女视频| 欧美日本亚洲视频在线播放| 国产精华一区二区三区| 99久久99久久久精品蜜桃| 午夜福利在线在线| 老鸭窝网址在线观看| 黑人操中国人逼视频| 18禁国产床啪视频网站| 校园春色视频在线观看| 午夜福利高清视频| 十八禁人妻一区二区| 不卡av一区二区三区| 成年女人永久免费观看视频| 一本精品99久久精品77| 搡老岳熟女国产| 99久久国产精品久久久| 成人特级黄色片久久久久久久| 91av网站免费观看| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 欧美中文日本在线观看视频| 久久久久久人人人人人| 在线观看午夜福利视频| 不卡av一区二区三区| 综合色av麻豆| 日韩大尺度精品在线看网址| 久久精品aⅴ一区二区三区四区| 亚洲成a人片在线一区二区| 免费在线观看影片大全网站| 老鸭窝网址在线观看| 香蕉久久夜色| 琪琪午夜伦伦电影理论片6080| 国产高清三级在线| 国产一区在线观看成人免费| 97超级碰碰碰精品色视频在线观看| 麻豆久久精品国产亚洲av| 国产精品久久视频播放| 免费看光身美女| 亚洲av电影在线进入| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 在线国产一区二区在线| 亚洲18禁久久av| 国产高清videossex| 五月伊人婷婷丁香| 亚洲成a人片在线一区二区| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站| 国产精品国产高清国产av| 叶爱在线成人免费视频播放| 波多野结衣巨乳人妻| 一二三四社区在线视频社区8| 亚洲第一电影网av| 色噜噜av男人的天堂激情| 色视频www国产| 91字幕亚洲| 男人和女人高潮做爰伦理| 精品不卡国产一区二区三区| 亚洲欧美日韩卡通动漫| 又黄又粗又硬又大视频| 天天添夜夜摸| 看黄色毛片网站| 九九热线精品视视频播放| 欧美最黄视频在线播放免费| 免费av毛片视频| 在线播放国产精品三级| 国产主播在线观看一区二区| 久久久久国产一级毛片高清牌| 久久久国产欧美日韩av| 桃色一区二区三区在线观看| 亚洲精品在线观看二区| 国产麻豆成人av免费视频| 国产高清视频在线播放一区| 欧美乱妇无乱码| 精品国产亚洲在线| 久久久久九九精品影院| 精品不卡国产一区二区三区| 国产精品 国内视频| 九九在线视频观看精品| 久久天躁狠狠躁夜夜2o2o| 日本成人三级电影网站| a在线观看视频网站| av福利片在线观看| 亚洲国产欧美一区二区综合| 熟女电影av网| 欧美乱妇无乱码| 岛国在线观看网站| 国产精品 国内视频| 国产亚洲精品久久久久久毛片| 午夜福利视频1000在线观看| 久久香蕉精品热| 1024手机看黄色片| 12—13女人毛片做爰片一| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| 麻豆国产av国片精品| 日韩免费av在线播放| 亚洲人与动物交配视频| 99久久久亚洲精品蜜臀av| 后天国语完整版免费观看| 老熟妇仑乱视频hdxx| bbb黄色大片| 嫩草影视91久久| 两人在一起打扑克的视频| 啦啦啦免费观看视频1| 精品久久久久久久久久久久久| 一二三四社区在线视频社区8| 国产欧美日韩精品一区二区| 亚洲av日韩精品久久久久久密| 99国产精品99久久久久| 亚洲狠狠婷婷综合久久图片| 久久久久免费精品人妻一区二区| 国产黄色小视频在线观看|