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

    楊樹(shù)MYB轉(zhuǎn)錄因子家族基因應(yīng)答鹽脅迫表達(dá)特性分析

    2017-11-10 10:43:26李晰妍張雪梅周博如
    植物研究 2017年3期
    關(guān)鍵詞:亞族亞組擬南芥

    李晰妍 趙 凱 張雪梅 周博如

    (東北林業(yè)大學(xué)林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150040)

    楊樹(shù)MYB轉(zhuǎn)錄因子家族基因應(yīng)答鹽脅迫表達(dá)特性分析

    李晰妍 趙 凱 張雪梅 周博如*

    (東北林業(yè)大學(xué)林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室,哈爾濱 150040)

    MYB轉(zhuǎn)錄因子家族是植物中數(shù)量最多的轉(zhuǎn)錄因子家族之一,在植物次生代謝調(diào)節(jié)、信號(hào)轉(zhuǎn)導(dǎo)和抗逆等生物過(guò)程起重要作用。根據(jù)MYB轉(zhuǎn)錄因子結(jié)構(gòu)域組成差異可分4個(gè)亞家族:即1R-MYB(MYB-relaed)、R2R3-MYB、3R-MYB和4R-MYB。其中,R2R3-MYB亞家族數(shù)量最多,可進(jìn)一步分為22個(gè)亞組;利用生物信息學(xué)分析楊樹(shù)MYB轉(zhuǎn)錄因子蛋白序列的保守結(jié)構(gòu)域、系統(tǒng)發(fā)生、基因組定位、氨基酸組成和理化性質(zhì)等;參照擬南芥MYB轉(zhuǎn)錄因子功能,預(yù)測(cè)楊樹(shù)MYB轉(zhuǎn)錄因子功能;基于84K楊轉(zhuǎn)錄組測(cè)序和RT-qPCR分析,從301個(gè)楊樹(shù)MYB轉(zhuǎn)錄因子基因中篩選出69個(gè)應(yīng)答鹽脅迫基因(P≤0.05)。其中,上調(diào)表達(dá)基因32個(gè),下調(diào)表達(dá)基因37個(gè)。該研究可為進(jìn)一步研究楊樹(shù)MYB家族基因功能提供參考依據(jù)。

    楊樹(shù);MYB轉(zhuǎn)錄因子;生物信息學(xué);RT-qPCR

    轉(zhuǎn)錄因子(Transcription factor)也稱(chēng)為反式作用因子,能夠與真核生物基因上游啟動(dòng)子序列中順式作用元件特異性結(jié)合,通過(guò)轉(zhuǎn)錄激活或抑制調(diào)控下游基因的表達(dá),對(duì)其下游基因的表達(dá)充當(dāng)開(kāi)關(guān)的角色。不同轉(zhuǎn)錄因子一般具有不同的功能區(qū)域,植物轉(zhuǎn)錄因子按照其結(jié)構(gòu)域的差異可分為WRKY、bZIP、EIL、MYB、NAC、GATA、AP2/ERF、SBP等轉(zhuǎn)錄因子家族[1]。

    MYB轉(zhuǎn)錄因子家族作為植物中最大的轉(zhuǎn)錄因子家族之一[2]。MYB轉(zhuǎn)錄因子蛋白C端含有轉(zhuǎn)錄激活域可與DNA或其他轉(zhuǎn)錄因子相互作用,N端具有由51-52個(gè)高度保守的氨基酸組成的MYB-DNA結(jié)合域。MYB-DNA結(jié)合域包含3個(gè)α螺旋,第二個(gè)和第三個(gè)螺旋組成helix-turn-helix(HTH)結(jié)構(gòu),形成一個(gè)疏水核心,對(duì)DNA序列的識(shí)別具有重要影響。MYB轉(zhuǎn)錄因子的MYB-DNA結(jié)合域通常由1~4個(gè)重復(fù)組成,依據(jù)其特點(diǎn)可分為4個(gè)亞族,即1R-MYB(MYB-related),R2R3-MYB,3R-MYB和4R-MYB[3~5]。其中1R-MYB亞族類(lèi)僅含有單一的MYB結(jié)構(gòu)域,在不同植物中功能有所差異。如同屬于MYB家族的擬南芥LHY,CCA1基因與擬南芥生物鐘與激素合成調(diào)節(jié)有關(guān)[6],蠟梅EU11603基因與植物抗逆相關(guān),參與植物的抗逆代謝途徑中ABA傳導(dǎo)[7];馬鈴薯過(guò)表達(dá)StMYB1R-1基因可加快氣孔關(guān)閉速度,增強(qiáng)植株抗干旱能力[8]。R2R3-MYB亞族蛋白質(zhì)是MYB轉(zhuǎn)錄因子家族中數(shù)量最多的一類(lèi),在調(diào)控次生代謝[9~10]、響應(yīng)赤霉素(GA)信號(hào)調(diào)節(jié)[11]、器官形成過(guò)程中發(fā)揮作用[12~13]。玉米ZmC1基因調(diào)控查耳酮合成及類(lèi)黃酮生物合成途徑[14],彩艷龍頭RCP1基因在成花過(guò)程中調(diào)控類(lèi)胡蘿卜素生物合成[15],金魚(yú)草MIXTA基因影響金魚(yú)草花瓣形態(tài)建成和花青素合成[16]。3R-MYB基因在動(dòng)物MYB家族數(shù)量較大,在植物中數(shù)量較少,1941年從禽急性成髓細(xì)胞白血病病毒AMV和E26分離出第一個(gè)MYB基因v-MYB即是3R-MYB基因[17]。植物中的也有類(lèi)似基因,如煙草NtMYB-A1,NtMYB-A2,NtMYB-B基因在細(xì)胞分裂G2期與M期發(fā)揮重要作用[18];水稻OsMYB3R-2基因應(yīng)答冷、干旱和鹽脅迫,過(guò)表達(dá)OsMYB3R-2基因?qū)е轮仓暧薪z分裂指數(shù)增加[19~20]。4R-MYB亞家族基因數(shù)量最少,擬南芥中雖然有該類(lèi)基因,但研究較少[21]。

    楊樹(shù)分布廣、適應(yīng)性強(qiáng)、生長(zhǎng)速度快,是我國(guó)重要的造林和用材樹(shù)種,也是研究林木生理和基因工程的模式植物。84K楊(Populusalba×Populusglandulosa)來(lái)自銀白楊和腺毛楊雜交,遺傳轉(zhuǎn)化率較高,是楊樹(shù)基因工程的重要受體樹(shù)種之一。為探明楊樹(shù)MYB轉(zhuǎn)錄因子的結(jié)構(gòu)特征和應(yīng)答鹽脅迫情況,通過(guò)植物轉(zhuǎn)錄因子數(shù)據(jù)庫(kù)(http://planttfdb.cbi.pku.edu.cn/index.php)篩選出楊樹(shù)301個(gè)MYB轉(zhuǎn)錄因子,利用生物信息學(xué)分析其結(jié)構(gòu)特點(diǎn);通過(guò)轉(zhuǎn)錄組測(cè)序分析84K楊中MYB轉(zhuǎn)錄因子基因應(yīng)答鹽脅迫情況,利用RT-qPCR對(duì)候選基因應(yīng)答鹽脅迫情況進(jìn)行驗(yàn)證,為篩選楊樹(shù)耐鹽關(guān)鍵基因和改良楊樹(shù)抗逆性提供參考依據(jù)。

    1 材料與方法

    從植物轉(zhuǎn)錄因子數(shù)據(jù)庫(kù)(http://planttfdb.cbi.pku.edu.cn/index.php)獲得301個(gè)楊樹(shù)MYB轉(zhuǎn)錄因子的cDNA及其編碼氨基酸序列,以及14個(gè)功能已知的擬南芥MYB轉(zhuǎn)錄因子氨基酸序列;利用BioEdit軟件分別進(jìn)行序列比對(duì),利用Weblogo 3.0對(duì)比對(duì)結(jié)果進(jìn)行基序分析(a;b),利用MEGA 5.0構(gòu)建進(jìn)化樹(shù)(NJ),設(shè)置Bootstrap為1 000次重復(fù);利用ProtParam在線分析,分析301個(gè)楊樹(shù)MYB轉(zhuǎn)錄因子家族蛋白氨基酸數(shù)目、分子量、理論等電點(diǎn)、脂肪族氨基酸指數(shù)及蛋白質(zhì)疏水性等理化性質(zhì);利用TargetP1.1 Server在線序列分析及Cello v2.5在線分析預(yù)測(cè)MYB轉(zhuǎn)錄因子亞細(xì)胞定位情況。

    將培養(yǎng)1個(gè)月的84K楊(Populusalba×Populusglandulosa)組培苗洗掉培養(yǎng)基,置于水中,在25℃室溫,14 h·d-1光照條件下培養(yǎng)1個(gè)月后,然后將生長(zhǎng)狀態(tài)相似的植株分為8組,分別用水和0.15 mol·L-1NaCl處理24 h,每個(gè)處理重復(fù)4次。采集葉片樣本,保存于-80℃冰箱,委托金唯智(GENEWIZ)公司進(jìn)行RNA提取和轉(zhuǎn)錄組測(cè)序。

    用瓊脂糖凝膠電泳檢測(cè)轉(zhuǎn)錄組測(cè)序RNA樣品完整性,用Takara的PrimeScriptTMRT reagent kit with gDNA Eraser試劑盒進(jìn)行cDNA合成。根據(jù)69個(gè)MYB候選基因序列,利用primer5.0設(shè)計(jì)定量PCR引物(表1)。其中FN356200(EF1),JM986590(ACT)為內(nèi)參基因[22]。按照TaKaRa公司SYBR?RPremix Ex TaqaTMП(Tli RNaseH Plus)實(shí)時(shí)定量試劑盒進(jìn)行實(shí)時(shí)定量RT-qPCR分析,所用程序?yàn)椋?5.0℃ 30 s,(95.0℃ 5 s、60.0℃ 34 s、95.0℃ 15 s)35個(gè)循環(huán),60.0℃ 60 s,95.0℃ 15 s,利用-△△Ct法進(jìn)行基因表達(dá)分析。

    2 結(jié)果與分析

    2.1 楊樹(shù)MYB轉(zhuǎn)錄因子分類(lèi)

    根據(jù)MYB轉(zhuǎn)錄因子保守結(jié)構(gòu)域差異將301個(gè)楊樹(shù)MYB轉(zhuǎn)錄因子蛋白分為4個(gè)亞族:1R-MYB(MYB-relaed)、 R2R3-MYB、3R-MYB和4R-MYB。

    表1 RT-qPCR引物

    其中,1R-MYB(MYB-relaed)亞族有89個(gè)成員、R2R3-MYB亞族有206個(gè)成員、3R-MYB亞族有5個(gè)成員、4R-MYB亞族僅有1個(gè)成員。將206個(gè)R2R3-MYB亞家族轉(zhuǎn)錄因子的保守結(jié)構(gòu)域(DNA binding domain)與14個(gè)功能已知的擬南芥MYB轉(zhuǎn)錄因子家族保守結(jié)構(gòu)域(DNA binding domain)進(jìn)行比對(duì),將R2R3-MYB亞族分為22個(gè)組,命名為R2R3-1~R2R3-22(圖1)。其中,有7個(gè)轉(zhuǎn)錄因子與擬南芥AtMYB61聚在一起,命名為R2R3-1亞組,該亞組基因多與氣孔調(diào)節(jié)和光合作用有關(guān)[23];有12個(gè)轉(zhuǎn)錄因子與擬南芥AtMYB4聚在一起,命名為R2R3-10亞組,該亞組基因多參與調(diào)節(jié)植物耐紫外輻射過(guò)程[24];有3個(gè)轉(zhuǎn)錄因子與擬南芥AtMYB11聚在一起,命名為R2R3-11亞組,該亞組基因可能參與與細(xì)胞發(fā)育和形態(tài)發(fā)生相關(guān)蛋白質(zhì)的功能分化[25];有9個(gè)轉(zhuǎn)錄因子與擬南芥AtMYB75、AtMYB90和AtMYB113聚在一起,命名為R2R3-14亞組,該亞組基因可能參與調(diào)節(jié)花青素的合成[26~27];有9個(gè)轉(zhuǎn)錄因子與擬南芥AtMYB59聚在一起,命名為2R3R-15亞組,該亞組基因可能參與調(diào)控植物細(xì)胞周期及根的發(fā)育過(guò)程[28]。

    圖1 楊樹(shù)R2R3-MYB亞族保守結(jié)構(gòu)域系統(tǒng)發(fā)生樹(shù) ▼表示14個(gè)功能已知的擬南芥MYB轉(zhuǎn)錄因子保守結(jié)構(gòu)域。Fig.1 Phylogenetic tree of Poplar R2R3-MYB subfamily DNA binding domains based on NJ method ▼represents DNA binding domains of the known functions Arabidopsis thaliana MYB transcription factors.

    2.2 楊樹(shù)MYB轉(zhuǎn)錄因子保守結(jié)構(gòu)域預(yù)測(cè)

    不同楊樹(shù)MYB轉(zhuǎn)錄因子亞家族成員中含有的保守結(jié)構(gòu)域數(shù)目不同,1R-MYB亞族成員有1個(gè)保守結(jié)構(gòu)域,R2R3-MYB亞族成員有2個(gè)保守結(jié)構(gòu)域,3R-MYB亞族成員有3個(gè)保守結(jié)構(gòu)域,4R-MYB亞族成員有4個(gè)保守結(jié)構(gòu)域。將301個(gè)楊樹(shù)MYB轉(zhuǎn)錄因子的第一結(jié)構(gòu)域和第二結(jié)構(gòu)域進(jìn)行比對(duì),可以發(fā)現(xiàn)其保守性均較高,但不同氨基酸的保守性存在差異。圖2MYB轉(zhuǎn)錄因子結(jié)構(gòu)域中,代表氨基酸的字母越大表示該氨基酸的保守程度越高,從圖2a中可以發(fā)現(xiàn)MYB第一結(jié)構(gòu)域的第4位,第32位的色氨酸(W)殘基極保守,但第4位的色氨酸(W)殘基極少數(shù)被苯丙氨酸(F)和甲硫氨酸(M)所取代;第32位的色氨酸(W)殘基極少數(shù)被苯丙氨酸(F)和亮氨酸(L)所取代。從圖2b可以發(fā)現(xiàn)MYB第二結(jié)構(gòu)域的第32位色氨酸(W)殘基始終極保守。總體而言,第一和第二結(jié)構(gòu)域內(nèi)其他位點(diǎn)的色氨酸(W)也較為保守,只有少數(shù)被苯丙氨酸(F),異亮氨酸(I),酪氨酸(Y)和丙氨酸(A)所取代,除此之外,保守結(jié)合域內(nèi)部還有很多保守的氨基酸殘基,比如賴(lài)氨酸(K),精氨酸(R),甘氨酸(G),天冬氨酸(D),谷氨酸(E),天冬酰胺(N),異亮氨酸(I),酪氨酸(Y)以及亮氨酸(L),共同維持著MYB轉(zhuǎn)錄因子家族蛋白結(jié)構(gòu)域的基本穩(wěn)定。

    圖2 楊樹(shù)MYB轉(zhuǎn)錄因子結(jié)構(gòu)域基序分析 a. 1R-MYB,R2R3-MYB,3R-MYB和4R-MYB亞族第一結(jié)構(gòu)域基序;b. R2R3-MYB,3R-MYB和4R-MYB亞族第二結(jié)構(gòu)域基序Fig.2 Sequence alignment of Poplar MYB transcription factors motif in the first domain and the second domaina. The first domain motif of 1R-MYB,R2R3-MYB,3R-MYB and 4R-MYB subfamilies; b. The second domain motif of R2R3-MYB,3R-MYB and 4R-MYB subfamilies

    2.3 楊樹(shù)MYB轉(zhuǎn)錄因子理化性質(zhì)分析

    1R-MYB亞族轉(zhuǎn)錄因子平均長(zhǎng)度為392 aa,平均分子量為43 558.7 kDa,理論等電點(diǎn)多在堿性范圍,細(xì)胞核定位占整個(gè)亞族的88%。2R3R-MYB亞族轉(zhuǎn)錄因子平均長(zhǎng)度為346 aa,平均分子量為38 861 kDa,其中最小的是2R3R-14亞組,平均長(zhǎng)度只有234.89 aa,平均分子量為26 830.79 kDa;最大的為2R3R-21亞組,平均長(zhǎng)度為1 389.25 aa,平均分子量為153 708.10 kDa,除第3、5、21亞組轉(zhuǎn)錄因子的等電點(diǎn)都為酸性外,其他2R3R-MYB亞組轉(zhuǎn)錄因子的等電點(diǎn)均為堿性。第1、3、17、19、20亞組多數(shù)為細(xì)胞外定位,過(guò)氧化酶體和高爾基體定位,其他亞組多數(shù)為細(xì)胞核定位,細(xì)胞核定位占整個(gè)亞族的79.6%。3R-MYB亞族、4R-MYB亞族多為堿性氨基酸,多數(shù)定位在細(xì)胞核外,過(guò)氧化酶體和高爾基體。脂肪族氨基酸指數(shù)是衡量蛋白質(zhì)穩(wěn)定性的指標(biāo),該指數(shù)越高蛋白質(zhì)穩(wěn)定性越強(qiáng),MYB轉(zhuǎn)錄因子蛋白的穩(wěn)定性均較強(qiáng)。蛋白質(zhì)的疏水性負(fù)值越大親水性越強(qiáng),不同MYB轉(zhuǎn)錄因子氨基酸疏水性值雖存在差距,但都屬親水蛋白。

    2.4 楊樹(shù)MYB基因應(yīng)答鹽脅迫分析

    根據(jù)84K楊轉(zhuǎn)錄組測(cè)序數(shù)據(jù),分析楊樹(shù)MYB基因應(yīng)答鹽脅迫情況,結(jié)果有69個(gè)基因表達(dá)發(fā)生了變化,并且P≤0.05,這些基因定義為鹽脅迫應(yīng)答基因(圖3),其中包括37個(gè)下調(diào)表達(dá)基因和32個(gè)上調(diào)表達(dá)基因,分別占家族基因總數(shù)的12.33%和10.67%。在37個(gè)下調(diào)表達(dá)基因中,14個(gè)屬1R-MYB亞族,23個(gè)屬2R3R-MYB亞族,下調(diào)表達(dá)超過(guò)15倍的基因有7個(gè):分別是Potri.002G260000,Potri.004G102600,Potri.005G112000,Potri.007G023800,Potri.009G116600,Potri.013G001000,Potri.018G049600;在32個(gè)上調(diào)表達(dá)基因中,8個(gè)屬1R-MYB亞族,24個(gè)屬2R3R-MYB亞族,上調(diào)表達(dá)超過(guò)15倍的基因有4個(gè):分別是Potri.008G122100,Potri.009G007100,Potri.010G149900,Potri.015G046200。

    圖3 楊樹(shù)MYB應(yīng)答鹽脅迫表達(dá)情況Fig.3 Expression of Poplar MYB genes in response to salt stress

    圖4 69個(gè)楊樹(shù)MYB基因的RT-qPCR分析Fig.4 Expression analysis of 69 Poplar MYB genes by RT-qPCR

    2.5 楊樹(shù)應(yīng)答鹽脅迫MYB基因的RT-qPCR驗(yàn)證

    用RT-qPCR驗(yàn)證69個(gè)應(yīng)答鹽脅迫的MYB基因的表達(dá)情況(圖4),結(jié)果表明RT-qPCR結(jié)果與轉(zhuǎn)錄組測(cè)序分析結(jié)果總體趨勢(shì)基本一致。如在轉(zhuǎn)錄組測(cè)序分析和RT-qPCR結(jié)果中,Potri.007G023800、Potri.018G049600下調(diào)表達(dá)量均超過(guò)15倍;Potri.008G122100、Potri.009G007100、Potri.010G149900、Potri.012G055600、Potri.015G046200上調(diào)表達(dá)量均超過(guò)15倍。但兩種方法檢測(cè)到的基因相對(duì)表達(dá)水平有所差異,如Potri.012G080400在轉(zhuǎn)錄組測(cè)序分析中下調(diào)表達(dá)1.11倍,RT-qPCR分析中下調(diào)表達(dá)11.14倍;而Potri.008G101400在轉(zhuǎn)錄組測(cè)序分析中上調(diào)表達(dá)3.77倍,在RT-qPCR分析中上調(diào)表達(dá)47.78倍。

    我們發(fā)現(xiàn)下調(diào)表達(dá)超過(guò)15倍的Potri.018G049600基因與擬南芥AtMYB7基因聚于2R3R-10亞組。有實(shí)驗(yàn)證明擬南芥AtMYB7基因在鹽脅迫下呈明顯下調(diào)表達(dá)[29];上調(diào)表達(dá)倍數(shù)超過(guò)15倍的Potri.008G122100,Potri.010G149900,Potri.012G055600,Potri.015G046200基因與擬南芥AtMYB2,AtMYB112基因同屬于2R3R-16亞組,擬南芥AtMYB2和AtMYB112基因在高鹽刺激下表達(dá)量上升,過(guò)量表達(dá)的轉(zhuǎn)基因植物耐鹽能力明顯提高[30~31]。Potri.009G007100基因與擬南芥AtMYB68基因?qū)儆?R3R-4亞組,擬南芥AtMYB68基因響應(yīng)環(huán)境脅迫[32];Potri.013G067500、Potri.014G081200、Potri.017G082500、Potri.019G045900基因與擬南芥AtMYB60、AtMYB94基因聚于2R3R-8亞組,擬南芥AtMYB60、AtMYB94基因參與鹽脅迫和干旱脅迫的生理調(diào)節(jié)過(guò)程[33~34];Potri.002G122600、Potri.002G128900、Potri.007G048900與擬南芥AtMYB44基因?qū)儆?R3R-18亞組,過(guò)表達(dá)AtMYB44通過(guò)SOS2負(fù)反饋環(huán)增強(qiáng)植株的耐鹽能力[35]。

    4 結(jié)論

    根據(jù)MYB轉(zhuǎn)錄因子結(jié)構(gòu)域數(shù)量差異將301個(gè)MYB轉(zhuǎn)錄因子基因分為4個(gè)亞族:1R-MYB(MYB-relaed)、R2R3-MYB、3R-MYB和4R-MYB。通過(guò)84K楊轉(zhuǎn)錄組測(cè)序篩選出69個(gè)應(yīng)答鹽脅迫的MYB轉(zhuǎn)錄因子基因,其中下調(diào)表達(dá)基因37個(gè),上調(diào)表達(dá)基因32個(gè)。RT-qPCR分析結(jié)果與轉(zhuǎn)錄組測(cè)序結(jié)果基本一致。本研究對(duì)明確MYB轉(zhuǎn)錄因子基因的應(yīng)答鹽脅迫表達(dá)特點(diǎn)和探明楊樹(shù)耐鹽分子機(jī)制具有參考價(jià)值。

    1.Kasuga M,Liu Q,Miura S,et al.Improving plant drought,salt,and freezing tolerance by gene transfer of a single stress-inducible transcription factor[J].Nature Biotechnology,1999,17(3):287-291.

    2.Frampton J.MYB transcription factors:their role in growth,differentiation and disease:proteins and cell regulation[M].Netherlands:Springer,2004.

    3.Rabinowicz P D,Braun E L,Wolfe A D,et al.Maize R2R3 Myb genes:sequence analysis reveals amplification in the higher plants[J].Genetics,1999,153(1):427-444.

    4.Chen Y H,Yang X Y,He K,et al.The MYB transcription factor superfamily ofArabidopsis:expression analysis and phylogenetic comparison with the Rice MYB Family[J].Plant Molecular Biology,2006,60(1):107-124.

    5.Ogate K,Morikawa S,Nakamura H,et al.Comparison of the free and DNA-complexed forms of the DNA-binding domain from c-MYB[J].Nature Structural Biology,1995,2(4):309-320.

    6.Carré I A,Kim J Y.MYB transcription factors in theArabidopsiscircadian clock[J].Journal of Experimental Botany,2002,53(374):1551-1557.

    7.羅莉莉,劉道鳳,茅允彬,等.蠟梅轉(zhuǎn)錄因子1R-Myb1的cDNA克隆、分子特性與表達(dá)載體的構(gòu)建[C].//中國(guó)園藝學(xué)會(huì)第八屆青年學(xué)術(shù)討論會(huì)論文集.上海:中國(guó)園藝學(xué)會(huì),2008.

    Luo L L,Liu D F,Mao Y B,et al.Chimonanthuspraecox(L.) link gene 1R-Myb1 with plant-specific SANT/MYB domain protein cDNA clones,molecular character and constructing expression vector[C].//China Horticultural Society Youth Seminar.Shanghai:Chinese Society for Horticultural Science,2008.

    8.Shin D,Moon S J,Han S,et al.Expression ofStMYB1R-1,a novel potato single MYB-like domain transcription factor,increases drought to lerance[J].Plant Physiology,2011,155(1):421-432.

    9.Zhang P F,Chopra S,Peterson T.A segmental gene duplication generated differentially expressed MYB-homologous genes inMaize[J].The Plant Cell,2000,12(12):2311-2322.

    10.Nesi N,Jond C,Debeaujon I,et al.TheArabidopsisTT2 gene encodes an R2R3 MYB domain protein that acts as a key determinant for proanthocyanidin accumulation in developing seed[J].The Plant Cell,2001,13(9):2099-2114.

    11.Gubler F,Kalla R,Roberts J K,et al.Gibberellin-regulated expression of a MYB gene in barley aleurone cells:evidence for MYB transactivation of a high-pI α-amylase gene promoter[J].The Plant Cell,1995,7(11):1879-1891.

    12.Steiner-lange S,Unte U S,Eckstein L,et al.Disruption ofArabidopsisthalianaMYB26 results in male sterility due to non-dehiscent anthers[J].The Plant Journal,2003,34(4):519-528.

    13.Penfield S,Meissner R C,Shoue D A,et al.MYB61 is required for mucilage deposition and extrusion in theArabidopsisseed coat[J].The Plant Cell,2001,13(12):2777-2791.

    14.Groteword E,Drummond B J,Bowen B,et al.The MYB-homologousPgene controls phlobaphene pigmentation inMaizefloral organs by directly activating a flavonoid biosynthetic gene subset[J].Cell,1994,76(3):543-553.

    15.Sagawa J M,Stanley L E,Lafountain A M,et al.An R2R3-MYB transcription factor regulates carotenoid pigmentation inMimuluslewisiiflowers[J].New Phytologist,2016,209(3):1049-1057.

    16.張旸,吳佳巖,吳雅妮,等.MYB轉(zhuǎn)錄因子基因MIXTA及其同源基因功能的研究進(jìn)展[J].中國(guó)農(nóng)業(yè)科學(xué),2016,49(7):1230-1241.

    Zhang Y,Wu J Y,Wu Y N,et al.Progresses and perspective of the function of MYB Transcription factorMIXTAand its orthologous gene[J].Scientia Agricultura Sinica,2016,49(7):1230-1241.

    17.Klempnauer K H,Gonda T J,Bishop J M.Nucleotide sequence of the retroviral leukemia genev-MYBand its cellular progenitorc-MYB:the architecture of a transduced oncogene[J].Cell,1982,31(2):453-463.

    18.Ito M,Araki S,Matsunaga S,et al.G2/M-phase-specific transcription during the plant cell cycle is mediated by c-MYB-like transcription factors[J].The Plant Cell,2001,13(8):1891-1905.

    19.Ma Q B,Dai X Y,Xu Y Y,et al.Enhanced tolerance to chilling stress inOsMYB3R-2 transgenic Rice is mediated by alteration in cell cycle and ectopic expression of stress genes[J].Plant Physiology,2009,150(1):244-256.

    20.Dai X Y,Xu Y Y,Ma Q B,et al.Overexpression of an R1R2R3 MYB gene,OsMYB3R-2,increases tolerance to freezing,Drought,and salt stress in transgenicArabidopsis[J].Plant Physiology,2007,143(4):1739-1751.

    21.Cai H S,Tian S,Dong H S.Large scaleinsilicoidentification ofMYBfamily genes from Wheat expressed sequence tags[J].Molecular Biotechnology,2012,52(2):184-192.

    22.Wang S J,Yao W J,Wei H R,et al.Expression patterns of ERF genes underlying abiotic stresses in Di-haploidPopulussimonii×P.nigra[J].The Scientific World Journal,2014,2014:745091.

    23.Liang Y K,Dubos C,Dodd I C,et al.AtMYB61,an R2R3-MYB transcription factor controlling stomatal aperture inArabidopsisthaliana[J].Current Biology,2005,15(13):1201-1206.

    24.Hemm M R,Herrmann K M,Chapple C.AtMYB4:a transcription factor general in the battle against UV[J].Trends in Plant Science,2001,6(4):135-136.

    25.Petroni K,Falasca G,Calvenzani V,et al.The AtMYB11 gene fromArabidopsisis expressed in meristematic cells and modulates growth in planta and organogenesis in vitro[J].Journal of Experimental Botany,2008,59(6):1201-1213.

    26.Bhargava A,Mansfield S D,Hall H C,et al.MYB75 functions in regulation of secondary cell wall formation in theArabidopsisinflorescence stem[J].Plant Physiology,2010,154(3):1428-1438.

    27.Rajagopalan R,Vaucheret H,Trejo J,et al.A diverse and evolutionarily fluid set of microRNAs inArabidopsisthaliana[J].Genes & Development,2006,20(24):3407-3425.

    28.Mu R L,Cao Y R,Liu Y F,et al.An R2R3-type transcription factor geneAtMYB59 regulates root growth and cell cycle progression inArabidopsis[J].Cell Research,2009,19(11):1291-1304.

    29.Kim J H,Hyun W Y,Nguyen H N,et al.AtMyb7,a subgroup 4 R2R3 Myb,negatively regulates ABA-induced inhibition of seed germination by blocking the expression of the bZIP transcription factor ABI5[J].Plant,Cell & Environment,2015,38(3):559-571.

    30.Urao T,Noji M A,Yamaguchi-shinozaki K,et al.A transcriptional activation domain of ATMYB2,a drought-inducibleArabidopsisMYB-related protein[J].The Plant Journal,1996,10(6):1145-1148.

    31.Lotkowska M E,Tohge T,Fernie A R,et al.TheArabidopsistranscription factor MYB112 promotes anthocyanin formation during salinity and under high light stress[J].Plant Physiology,2015,169(3):1862-1880.

    32.Feng C P,Andreasson E,Maslak A,et al.ArabidopsisMYB68 in development and responses to environmental cues[J].Plant Science,2004,167(5):1099-1107.

    33.Cominelli E,Galbiati M,Vavasseur A,et al.A guard-cell-specific MYB transcription factor regulates stomatal movements and plant drought tolerance[J].Current Biology Cb,2005,15(13):1196-1200.

    34.Lee S B,Kim J,Suh M C.Cuticular Wax biosynthesis is up-regulated by the MYB94 transcription factor inArabidopsis[J].Plant & Cell Physiology,2015,56(1):48-60.

    35.Jung C,Seo J S,Han S W,et al.Overexpression ofAtMYB44 enhances stomatal closure to confer abiotic stress tolerance in transgenicArabidopsis[J].Plant Physiology,2008,146(2):623-635.

    National 863 project funding(2013AA102701)

    introduction:LI Xi-Yan(1992—),female,graduate,major in forest resistance mechanism.

    date:2017-02-16

    ExpressionAnalysisofPoplarMYBTranscriptionFactorGeneFamilyinResponsetoSaltStress

    LI Xi-Yan ZHAO Kai ZHANG Xue-Mei ZHOU Bo-Ru*

    (State Key Laboratory of Tree Genetics and Breeding,Northeast Forestry University,Harbin 150040)

    MYB transcription factor family is one of the most important transcription factors in plants, and plays an important role in the process of plant secondary metabolism regulation, signal transduction and stress resistance. According to the difference of domain structure of MYB transcription factor, it can be divided into 4 subfamilies: 1R-MYB(MYB-relaed), R2R3-MYB, 3R-MYB and 4R-MYB. Among them, the R2R3-MYB sub family with the largest number of members can be further divided into 22 subgroups. Using bioinformatics, the characterizations of conserved domain ofPoplarMYB transcription factors, genomic location, amino acid composition and physicochemical properties were analyzed. According to the Arabidopsis MYB transcription factor, the functions ofPoplarMYB transcription factor were predicted. Based on the transcriptome sequencing and RT-qPCR, 69 salt-stress responding genes(P≤0.05) were screened from 301 MYB genes, including 32 up-regulated genes and 37 down regulated genes. The study could provide a reference for further study on the function of MYB family genes inPoplar.

    Poplar;MYB transcription factors;bioinformatics;RT-qPCR

    國(guó)家863課題資助(2013AA102701)

    李晰妍(1992—),女,碩士研究生,主要從事林木抗逆機(jī)理方面研究。

    * 通信作者:E-mail:boruzhou@yahoo.com

    2017-02-16

    * Corresponding author:E-mail:boruzhou@yahoo.com

    S792.11

    A

    10.7525/j.issn.1673-5102.2017.03.013

    猜你喜歡
    亞族亞組擬南芥
    基于Meta分析的黃酮類(lèi)化合物對(duì)奶牛生產(chǎn)性能和血清免疫指標(biāo)影響的研究
    慢性阻塞性肺疾病患者膈肌移動(dòng)度分析
    擬南芥:活得粗糙,才讓我有了上太空的資格
    二穗短柄草CYP72A亞族成員表達(dá)分析及亞細(xì)胞定位
    槭葉鐵線蓮亞組的研究進(jìn)展
    園林科技(2021年3期)2022-01-19 03:17:32
    苦蕎蛋白磷酸酶2C家族的鑒定及表達(dá)分析
    辣椒HD-Zip基因家族鑒定、系統(tǒng)進(jìn)化及表達(dá)分析
    冠心病患者腸道菌群變化的研究 (正文見(jiàn)第45 頁(yè))
    尿黑酸對(duì)擬南芥酪氨酸降解缺陷突變體sscd1的影響
    兩種LED光源作為擬南芥生長(zhǎng)光源的應(yīng)用探究
    黄片小视频在线播放| 日本爱情动作片www.在线观看| 视频在线观看一区二区三区| 香蕉丝袜av| 大陆偷拍与自拍| 午夜福利一区二区在线看| 亚洲三级黄色毛片| 成年动漫av网址| 亚洲精品久久成人aⅴ小说| 精品国产一区二区三区四区第35| av视频免费观看在线观看| 日韩欧美一区视频在线观看| 一区二区日韩欧美中文字幕| 亚洲,欧美精品.| 国产精品二区激情视频| 久久久国产精品麻豆| 中文字幕人妻丝袜一区二区 | 久久热在线av| 日韩三级伦理在线观看| 大码成人一级视频| 成年女人毛片免费观看观看9 | 久久久久久久久久久免费av| 亚洲欧美成人精品一区二区| 日韩一区二区视频免费看| 日韩一区二区视频免费看| 成年女人在线观看亚洲视频| tube8黄色片| 97在线人人人人妻| 午夜福利视频在线观看免费| 纯流量卡能插随身wifi吗| 久久久久久人妻| 97在线人人人人妻| 亚洲精品在线美女| 欧美亚洲日本最大视频资源| 侵犯人妻中文字幕一二三四区| 水蜜桃什么品种好| 97人妻天天添夜夜摸| 久久久久国产精品人妻一区二区| 一二三四在线观看免费中文在| 欧美成人午夜精品| 国产成人aa在线观看| 欧美日韩精品网址| 97精品久久久久久久久久精品| 老司机影院成人| 亚洲国产色片| 黄色视频在线播放观看不卡| 春色校园在线视频观看| av电影中文网址| 欧美日韩国产mv在线观看视频| 精品少妇内射三级| 久久精品国产综合久久久| 成人漫画全彩无遮挡| 青春草亚洲视频在线观看| 黄片播放在线免费| 日本爱情动作片www.在线观看| 高清不卡的av网站| 婷婷色综合www| videossex国产| 中文欧美无线码| 最近最新中文字幕大全免费视频 | 日日撸夜夜添| 欧美日韩视频精品一区| 久久99精品国语久久久| 国产一区二区三区av在线| 精品国产一区二区久久| 最新中文字幕久久久久| 99久久中文字幕三级久久日本| 精品一区二区免费观看| 亚洲综合色网址| 国产精品人妻久久久影院| 欧美另类一区| 一级毛片电影观看| 午夜激情久久久久久久| 日韩制服骚丝袜av| 99久久中文字幕三级久久日本| 国产激情久久老熟女| 中文字幕最新亚洲高清| 校园人妻丝袜中文字幕| 亚洲色图综合在线观看| 啦啦啦视频在线资源免费观看| 午夜日韩欧美国产| 高清欧美精品videossex| 下体分泌物呈黄色| 亚洲国产精品成人久久小说| 国产精品免费视频内射| 久久久久国产精品人妻一区二区| 青春草国产在线视频| 亚洲美女搞黄在线观看| 亚洲国产精品国产精品| 一级爰片在线观看| 国产午夜精品一二区理论片| 91精品三级在线观看| 精品国产露脸久久av麻豆| 制服丝袜香蕉在线| 国产高清不卡午夜福利| 亚洲av日韩在线播放| 街头女战士在线观看网站| 久久久精品国产亚洲av高清涩受| 中文字幕色久视频| 色哟哟·www| 青青草视频在线视频观看| 啦啦啦中文免费视频观看日本| 黄片播放在线免费| av卡一久久| 五月天丁香电影| 国产精品免费视频内射| 亚洲av综合色区一区| www.自偷自拍.com| 免费黄频网站在线观看国产| 秋霞伦理黄片| 十八禁高潮呻吟视频| 晚上一个人看的免费电影| 丝袜喷水一区| 日韩精品免费视频一区二区三区| 在线精品无人区一区二区三| 大陆偷拍与自拍| 国产免费福利视频在线观看| 免费人妻精品一区二区三区视频| 看非洲黑人一级黄片| 在线观看三级黄色| 中国国产av一级| 9191精品国产免费久久| 在线 av 中文字幕| 精品少妇一区二区三区视频日本电影 | 亚洲三区欧美一区| 女人被躁到高潮嗷嗷叫费观| 国产精品麻豆人妻色哟哟久久| 亚洲精华国产精华液的使用体验| 久久婷婷青草| 免费看不卡的av| 五月开心婷婷网| 久久99一区二区三区| 欧美xxⅹ黑人| 成人毛片a级毛片在线播放| 母亲3免费完整高清在线观看 | 成年女人毛片免费观看观看9 | 日韩不卡一区二区三区视频在线| 日韩大片免费观看网站| 亚洲欧洲精品一区二区精品久久久 | 欧美日本中文国产一区发布| 亚洲成人av在线免费| 久久久久久久久久久久大奶| 哪个播放器可以免费观看大片| 精品久久蜜臀av无| 黑人欧美特级aaaaaa片| 日韩欧美精品免费久久| 欧美 亚洲 国产 日韩一| 一级黄片播放器| 丝袜在线中文字幕| 国产色婷婷99| 日韩,欧美,国产一区二区三区| 免费观看在线日韩| 只有这里有精品99| 精品人妻偷拍中文字幕| 中文字幕人妻熟女乱码| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲色图综合在线观看| 在线 av 中文字幕| 乱人伦中国视频| 免费看不卡的av| 免费观看在线日韩| 亚洲色图 男人天堂 中文字幕| 各种免费的搞黄视频| 亚洲情色 制服丝袜| 欧美 亚洲 国产 日韩一| 国产一级毛片在线| 国产不卡av网站在线观看| 老司机影院毛片| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 99九九在线精品视频| 黄片播放在线免费| 国产精品久久久久久av不卡| 亚洲精品av麻豆狂野| 一区二区三区激情视频| 少妇精品久久久久久久| 日本黄色日本黄色录像| 久久久亚洲精品成人影院| 美女国产视频在线观看| 91在线精品国自产拍蜜月| 高清欧美精品videossex| 香蕉精品网在线| 国产午夜精品一二区理论片| 久久精品国产亚洲av高清一级| 免费在线观看黄色视频的| 亚洲一级一片aⅴ在线观看| 日本爱情动作片www.在线观看| 2021少妇久久久久久久久久久| 国产精品国产三级国产专区5o| 男人舔女人的私密视频| 国产又色又爽无遮挡免| 美女脱内裤让男人舔精品视频| 午夜日本视频在线| 日日撸夜夜添| 高清不卡的av网站| 韩国av在线不卡| 日韩视频在线欧美| 性高湖久久久久久久久免费观看| 久久免费观看电影| 99热网站在线观看| 国产精品三级大全| 日日爽夜夜爽网站| 久久精品久久久久久久性| 免费少妇av软件| 国产成人精品福利久久| 大陆偷拍与自拍| 99热国产这里只有精品6| 丰满迷人的少妇在线观看| 欧美人与性动交α欧美软件| 超碰97精品在线观看| 午夜免费观看性视频| 久久国产精品大桥未久av| 免费在线观看黄色视频的| 性少妇av在线| 观看美女的网站| 女性被躁到高潮视频| 丁香六月天网| 国产一区二区激情短视频 | 99国产综合亚洲精品| 久久精品国产自在天天线| 国产人伦9x9x在线观看 | 熟女电影av网| 又大又黄又爽视频免费| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 国产免费又黄又爽又色| videosex国产| 日本av免费视频播放| 午夜久久久在线观看| 婷婷色综合www| videos熟女内射| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 久久国产亚洲av麻豆专区| av女优亚洲男人天堂| av在线老鸭窝| 免费观看av网站的网址| 成人国产av品久久久| 国产av码专区亚洲av| 日本wwww免费看| 午夜91福利影院| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看| 观看av在线不卡| 在线观看免费高清a一片| 晚上一个人看的免费电影| 午夜免费观看性视频| 电影成人av| 国产成人精品婷婷| 热99国产精品久久久久久7| 国产亚洲精品第一综合不卡| 成人亚洲精品一区在线观看| 麻豆av在线久日| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品aⅴ在线观看| 女人精品久久久久毛片| 99久久综合免费| 午夜精品国产一区二区电影| 看免费成人av毛片| 日韩中字成人| xxxhd国产人妻xxx| 妹子高潮喷水视频| 国产精品成人在线| 亚洲精品乱久久久久久| 国产成人免费无遮挡视频| 亚洲精品国产av蜜桃| 亚洲久久久国产精品| 亚洲美女黄色视频免费看| 9色porny在线观看| 啦啦啦啦在线视频资源| 久久精品aⅴ一区二区三区四区 | 亚洲伊人久久精品综合| 9191精品国产免费久久| 久久精品国产a三级三级三级| 国产高清不卡午夜福利| 日本vs欧美在线观看视频| 一区二区三区乱码不卡18| 亚洲av在线观看美女高潮| 国产高清不卡午夜福利| 中国三级夫妇交换| 日本欧美国产在线视频| 精品国产一区二区三区四区第35| 在线天堂最新版资源| 美女国产视频在线观看| 嫩草影院入口| av在线老鸭窝| 人体艺术视频欧美日本| 看非洲黑人一级黄片| 黄色毛片三级朝国网站| 天天躁夜夜躁狠狠躁躁| 国产精品.久久久| 五月伊人婷婷丁香| 91精品伊人久久大香线蕉| 999精品在线视频| 99re6热这里在线精品视频| 美女主播在线视频| 精品少妇内射三级| 国产日韩欧美亚洲二区| 欧美xxⅹ黑人| 午夜久久久在线观看| 人妻少妇偷人精品九色| 久久精品熟女亚洲av麻豆精品| 纵有疾风起免费观看全集完整版| 免费久久久久久久精品成人欧美视频| 免费日韩欧美在线观看| 久久久精品免费免费高清| 只有这里有精品99| 欧美日韩视频高清一区二区三区二| 久久久久久久大尺度免费视频| 亚洲中文av在线| 欧美bdsm另类| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 精品一区二区三区四区五区乱码 | 国产精品香港三级国产av潘金莲 | 蜜桃国产av成人99| 伦精品一区二区三区| 免费高清在线观看视频在线观看| 丰满乱子伦码专区| 黄色视频在线播放观看不卡| 天天影视国产精品| 国产成人精品久久久久久| 亚洲成国产人片在线观看| av卡一久久| 一区二区日韩欧美中文字幕| 国产福利在线免费观看视频| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说| 亚洲视频免费观看视频| 日韩av不卡免费在线播放| 久久影院123| 另类精品久久| 亚洲男人天堂网一区| 欧美亚洲 丝袜 人妻 在线| 久久久久人妻精品一区果冻| 国产黄色视频一区二区在线观看| 成年动漫av网址| 中文欧美无线码| 9色porny在线观看| 如何舔出高潮| 夫妻性生交免费视频一级片| 国产一区亚洲一区在线观看| 春色校园在线视频观看| 一级爰片在线观看| 久久久久久久久久久久大奶| 中国国产av一级| 国产成人欧美| 免费高清在线观看日韩| 国产一区亚洲一区在线观看| 97精品久久久久久久久久精品| 亚洲第一青青草原| 最新中文字幕久久久久| 嫩草影院入口| 新久久久久国产一级毛片| 男女免费视频国产| 国产人伦9x9x在线观看 | 国产黄频视频在线观看| 国产av国产精品国产| 国产在线一区二区三区精| 亚洲欧美一区二区三区黑人 | 精品视频人人做人人爽| 91午夜精品亚洲一区二区三区| 精品少妇久久久久久888优播| 久久精品国产a三级三级三级| 纵有疾风起免费观看全集完整版| 午夜久久久在线观看| 中文字幕制服av| 狠狠婷婷综合久久久久久88av| 亚洲精品国产一区二区精华液| 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| 久久久久视频综合| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱| 日日撸夜夜添| 一级a爱视频在线免费观看| 深夜精品福利| 日本vs欧美在线观看视频| 9热在线视频观看99| 91精品三级在线观看| av片东京热男人的天堂| 免费黄色在线免费观看| 久久精品国产亚洲av涩爱| 七月丁香在线播放| 99久久综合免费| 久久久国产欧美日韩av| 欧美日韩成人在线一区二区| 少妇猛男粗大的猛烈进出视频| 在线观看一区二区三区激情| 国产成人精品无人区| 久久久精品94久久精品| 亚洲少妇的诱惑av| 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 免费黄网站久久成人精品| 熟女电影av网| 青春草亚洲视频在线观看| 亚洲精品一区蜜桃| 成年美女黄网站色视频大全免费| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 一区二区日韩欧美中文字幕| 亚洲熟女精品中文字幕| 精品卡一卡二卡四卡免费| 巨乳人妻的诱惑在线观看| 免费高清在线观看视频在线观看| 中文字幕av电影在线播放| 丝袜喷水一区| 亚洲人成电影观看| 啦啦啦中文免费视频观看日本| 久久久久久久亚洲中文字幕| 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜爱| 精品国产一区二区三区四区第35| 亚洲一区二区三区欧美精品| 一级毛片我不卡| 国产成人午夜福利电影在线观看| 亚洲av男天堂| 国产精品久久久久久精品古装| 天堂中文最新版在线下载| 中文字幕av电影在线播放| 欧美变态另类bdsm刘玥| 日韩大片免费观看网站| 婷婷色av中文字幕| 日韩人妻精品一区2区三区| 久久精品国产综合久久久| 91精品三级在线观看| 欧美中文综合在线视频| 精品亚洲乱码少妇综合久久| 久久免费观看电影| 王馨瑶露胸无遮挡在线观看| 欧美黄色片欧美黄色片| 国产亚洲午夜精品一区二区久久| 男女午夜视频在线观看| 如日韩欧美国产精品一区二区三区| 成年人午夜在线观看视频| 少妇人妻精品综合一区二区| 亚洲第一av免费看| av免费观看日本| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| www.精华液| 午夜福利视频在线观看免费| 国产 一区精品| 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 黄色怎么调成土黄色| av视频免费观看在线观看| 免费不卡的大黄色大毛片视频在线观看| 18禁观看日本| 欧美亚洲 丝袜 人妻 在线| 男人爽女人下面视频在线观看| 成人亚洲欧美一区二区av| av国产精品久久久久影院| 国产伦理片在线播放av一区| 有码 亚洲区| 免费少妇av软件| 日日撸夜夜添| 如何舔出高潮| 国产黄色免费在线视频| 亚洲国产精品成人久久小说| 日韩欧美一区视频在线观看| 一区二区av电影网| 超碰97精品在线观看| 日韩一本色道免费dvd| 亚洲四区av| 国产免费现黄频在线看| 久久久久久久精品精品| 久久亚洲国产成人精品v| 国产精品欧美亚洲77777| 在线观看一区二区三区激情| 精品一区二区免费观看| 中文字幕制服av| 街头女战士在线观看网站| 侵犯人妻中文字幕一二三四区| 国产成人精品久久久久久| 午夜日韩欧美国产| 国产毛片在线视频| 欧美另类一区| 午夜福利一区二区在线看| 亚洲国产欧美在线一区| 18禁观看日本| 欧美人与性动交α欧美精品济南到 | 美女xxoo啪啪120秒动态图| 亚洲欧美色中文字幕在线| 日韩伦理黄色片| 在线观看www视频免费| 男女啪啪激烈高潮av片| 曰老女人黄片| 亚洲一码二码三码区别大吗| 一本大道久久a久久精品| av有码第一页| 最近中文字幕高清免费大全6| 两个人免费观看高清视频| 18禁裸乳无遮挡动漫免费视频| 免费av中文字幕在线| 亚洲国产色片| 国产又色又爽无遮挡免| 国产亚洲欧美精品永久| 飞空精品影院首页| 两性夫妻黄色片| www日本在线高清视频| 亚洲视频免费观看视频| 成人国语在线视频| 国产成人欧美| 亚洲色图 男人天堂 中文字幕| 春色校园在线视频观看| 国产伦理片在线播放av一区| 有码 亚洲区| www.自偷自拍.com| 制服诱惑二区| 婷婷色综合www| 成年人免费黄色播放视频| 丰满迷人的少妇在线观看| 亚洲国产毛片av蜜桃av| 久久精品亚洲av国产电影网| 精品少妇一区二区三区视频日本电影 | 欧美日韩一级在线毛片| 亚洲国产欧美日韩在线播放| 亚洲在久久综合| 在线观看国产h片| 久久精品国产a三级三级三级| 最新中文字幕久久久久| 成年av动漫网址| 国产精品无大码| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| 老司机影院成人| 看免费成人av毛片| 久久青草综合色| 中文字幕亚洲精品专区| 69精品国产乱码久久久| 麻豆乱淫一区二区| 亚洲精品日韩在线中文字幕| 丝袜美足系列| 久久精品人人爽人人爽视色| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 成人影院久久| 精品人妻一区二区三区麻豆| 一级片免费观看大全| 午夜激情av网站| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 亚洲成av片中文字幕在线观看 | 免费高清在线观看日韩| 精品久久蜜臀av无| 午夜激情av网站| 午夜免费观看性视频| 国产有黄有色有爽视频| 亚洲综合精品二区| 国产成人精品福利久久| av天堂久久9| 一区福利在线观看| 国产亚洲午夜精品一区二区久久| 中文字幕亚洲精品专区| 日韩av不卡免费在线播放| 亚洲欧美日韩另类电影网站| 日本-黄色视频高清免费观看| 夜夜骑夜夜射夜夜干| 天天操日日干夜夜撸| 欧美日韩国产mv在线观看视频| 国产精品不卡视频一区二区| 国产精品女同一区二区软件| 成人午夜精彩视频在线观看| 成人国产av品久久久| 黄片无遮挡物在线观看| av女优亚洲男人天堂| 精品人妻偷拍中文字幕| 午夜激情av网站| 国产精品偷伦视频观看了| 一级黄片播放器| 80岁老熟妇乱子伦牲交| 桃花免费在线播放| 三级国产精品片| 久久99一区二区三区| 日韩 亚洲 欧美在线| 精品亚洲成a人片在线观看| 久久久久视频综合| 精品人妻一区二区三区麻豆| av片东京热男人的天堂| 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 最近的中文字幕免费完整| 男的添女的下面高潮视频| 精品少妇一区二区三区视频日本电影 | av网站免费在线观看视频| 啦啦啦视频在线资源免费观看| 成人亚洲精品一区在线观看| 综合色丁香网| 亚洲精品美女久久av网站| 亚洲精品国产色婷婷电影| videos熟女内射| 又粗又硬又长又爽又黄的视频| 99久久人妻综合| 街头女战士在线观看网站| 久久久久国产精品人妻一区二区| 久久99蜜桃精品久久| 9热在线视频观看99| 国产极品天堂在线| 五月天丁香电影| 亚洲精品久久久久久婷婷小说| 国产精品99久久99久久久不卡 | 叶爱在线成人免费视频播放| 观看美女的网站| 久久女婷五月综合色啪小说| 久久 成人 亚洲| 亚洲第一区二区三区不卡| 好男人视频免费观看在线| 国产福利在线免费观看视频| 亚洲国产最新在线播放| 欧美精品一区二区大全| 夜夜骑夜夜射夜夜干|