• <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)用探究
    午夜激情久久久久久久| 国产精品99久久久久久久久| 国产高清有码在线观看视频| 国产一区二区三区综合在线观看 | 午夜激情久久久久久久| 大香蕉久久网| 亚洲精品成人av观看孕妇| a级毛片免费高清观看在线播放| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 黑人高潮一二区| 亚洲成人av在线免费| 十八禁高潮呻吟视频 | 久久青草综合色| 国产免费一级a男人的天堂| 精品亚洲乱码少妇综合久久| 国产免费视频播放在线视频| 蜜臀久久99精品久久宅男| 国产免费又黄又爽又色| 99久久综合免费| 少妇被粗大猛烈的视频| 中国美白少妇内射xxxbb| 热re99久久精品国产66热6| 精品午夜福利在线看| 成人亚洲欧美一区二区av| 十八禁高潮呻吟视频 | 免费观看的影片在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久久久久大av| 女性生殖器流出的白浆| 成年av动漫网址| 成人无遮挡网站| 久久久久久久久大av| 天堂中文最新版在线下载| 精品亚洲乱码少妇综合久久| 亚洲高清免费不卡视频| 亚洲综合色惰| 日韩,欧美,国产一区二区三区| 如何舔出高潮| 黄色毛片三级朝国网站 | 中文资源天堂在线| 韩国av在线不卡| 狂野欧美激情性xxxx在线观看| av网站免费在线观看视频| 国产色婷婷99| 久久人人爽人人片av| 国产亚洲午夜精品一区二区久久| 人妻夜夜爽99麻豆av| 男人添女人高潮全过程视频| 久久久久精品性色| 久久久久久久精品精品| 久久久久精品久久久久真实原创| av有码第一页| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产精品成人久久小说| 国产精品女同一区二区软件| freevideosex欧美| 日韩成人伦理影院| 能在线免费看毛片的网站| 亚洲国产精品一区二区三区在线| 午夜视频国产福利| 建设人人有责人人尽责人人享有的| 亚洲av男天堂| 久久99一区二区三区| 久久婷婷青草| 美女福利国产在线| 久久99蜜桃精品久久| 中文资源天堂在线| av在线app专区| 在线 av 中文字幕| 日韩成人伦理影院| 午夜老司机福利剧场| 亚洲国产欧美在线一区| 美女大奶头黄色视频| kizo精华| 秋霞在线观看毛片| 综合色丁香网| 国产黄频视频在线观看| 夫妻性生交免费视频一级片| 在线观看免费视频网站a站| 狂野欧美激情性bbbbbb| 欧美国产精品一级二级三级 | 少妇人妻 视频| 99视频精品全部免费 在线| 国产日韩欧美亚洲二区| 插阴视频在线观看视频| 中国三级夫妇交换| 男女边摸边吃奶| 日日摸夜夜添夜夜添av毛片| 麻豆精品久久久久久蜜桃| 午夜av观看不卡| 狠狠精品人妻久久久久久综合| 午夜激情久久久久久久| 在线观看av片永久免费下载| 在线播放无遮挡| 色视频www国产| 午夜免费鲁丝| av福利片在线| 亚洲欧美一区二区三区黑人 | 亚洲国产精品999| 欧美成人午夜免费资源| 大香蕉97超碰在线| 亚洲国产欧美在线一区| a级毛片在线看网站| 国产又色又爽无遮挡免| 亚洲av不卡在线观看| 亚洲欧美日韩卡通动漫| 亚洲欧美精品专区久久| 久久久欧美国产精品| 国产成人aa在线观看| 最近中文字幕高清免费大全6| 啦啦啦在线观看免费高清www| 亚洲一区二区三区欧美精品| 午夜免费男女啪啪视频观看| 晚上一个人看的免费电影| 午夜福利视频精品| 日韩中文字幕视频在线看片| 亚洲av中文av极速乱| 黑丝袜美女国产一区| 亚洲国产精品999| 97在线人人人人妻| 国产 精品1| 99久久综合免费| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人爽人人夜夜| 日韩 亚洲 欧美在线| 啦啦啦视频在线资源免费观看| 亚洲精品色激情综合| 18+在线观看网站| 极品教师在线视频| 国语对白做爰xxxⅹ性视频网站| 欧美bdsm另类| 亚洲第一区二区三区不卡| 看免费成人av毛片| 免费少妇av软件| 亚洲精品成人av观看孕妇| 天美传媒精品一区二区| 国产伦精品一区二区三区视频9| 欧美最新免费一区二区三区| 赤兔流量卡办理| 狠狠精品人妻久久久久久综合| 能在线免费看毛片的网站| 美女福利国产在线| 男女国产视频网站| av播播在线观看一区| 99精国产麻豆久久婷婷| 久久国产精品男人的天堂亚洲 | 这个男人来自地球电影免费观看 | 国产精品不卡视频一区二区| 欧美成人精品欧美一级黄| 我要看黄色一级片免费的| 亚洲国产精品一区三区| 国产男女超爽视频在线观看| 啦啦啦啦在线视频资源| 一区二区三区四区激情视频| 欧美少妇被猛烈插入视频| 亚洲精品aⅴ在线观看| 欧美日本中文国产一区发布| 啦啦啦啦在线视频资源| 国产亚洲精品久久久com| 精品卡一卡二卡四卡免费| 国产亚洲午夜精品一区二区久久| 久久99蜜桃精品久久| 精品酒店卫生间| 少妇的逼水好多| av在线老鸭窝| 18禁动态无遮挡网站| 亚洲欧美一区二区三区国产| 一区二区三区乱码不卡18| 国产av国产精品国产| 青青草视频在线视频观看| 少妇丰满av| 欧美日韩视频精品一区| 精品国产一区二区三区久久久樱花| 久久av网站| 国产免费又黄又爽又色| 又爽又黄a免费视频| 久久精品国产亚洲网站| 自线自在国产av| 精品久久久噜噜| 中国美白少妇内射xxxbb| 黄色配什么色好看| 天堂中文最新版在线下载| 欧美亚洲 丝袜 人妻 在线| 草草在线视频免费看| 成人国产av品久久久| 久久国产精品男人的天堂亚洲 | 国产精品秋霞免费鲁丝片| 亚洲精品自拍成人| 你懂的网址亚洲精品在线观看| 丁香六月天网| 99热这里只有精品一区| 欧美 亚洲 国产 日韩一| 七月丁香在线播放| 亚洲精品亚洲一区二区| 人人妻人人爽人人添夜夜欢视频 | 亚洲欧美成人精品一区二区| 午夜91福利影院| 亚洲精品,欧美精品| 久久99蜜桃精品久久| 日韩中文字幕视频在线看片| 国产在视频线精品| 亚洲成人一二三区av| 国产91av在线免费观看| 丰满迷人的少妇在线观看| 国产视频内射| 久久6这里有精品| 久久精品国产亚洲网站| 久久婷婷青草| 我的老师免费观看完整版| 制服丝袜香蕉在线| 亚州av有码| 能在线免费看毛片的网站| 国产69精品久久久久777片| 欧美老熟妇乱子伦牲交| 久久久久久久久久久久大奶| 夫妻午夜视频| 91久久精品国产一区二区成人| 少妇猛男粗大的猛烈进出视频| 欧美丝袜亚洲另类| 亚洲欧美成人精品一区二区| 国产色婷婷99| 久久女婷五月综合色啪小说| 久久精品国产自在天天线| 青青草视频在线视频观看| 久久精品国产亚洲网站| 丝袜喷水一区| 一级毛片黄色毛片免费观看视频| 街头女战士在线观看网站| 精品一区在线观看国产| 国产精品一二三区在线看| 国产成人午夜福利电影在线观看| 日韩一本色道免费dvd| 午夜久久久在线观看| 免费观看的影片在线观看| 如何舔出高潮| 波野结衣二区三区在线| 人人妻人人澡人人看| 亚洲,一卡二卡三卡| 国产亚洲精品久久久com| 免费黄频网站在线观看国产| 人人妻人人澡人人爽人人夜夜| av福利片在线| 啦啦啦中文免费视频观看日本| 婷婷色av中文字幕| 欧美精品高潮呻吟av久久| 久久久久久久久久人人人人人人| av免费观看日本| 久久精品国产自在天天线| 日本午夜av视频| 伊人亚洲综合成人网| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| 国产成人免费观看mmmm| 成人漫画全彩无遮挡| 国模一区二区三区四区视频| 欧美丝袜亚洲另类| 午夜免费男女啪啪视频观看| 中文资源天堂在线| 久久久久久久久久久免费av| .国产精品久久| 两个人的视频大全免费| 色婷婷av一区二区三区视频| 校园人妻丝袜中文字幕| 秋霞在线观看毛片| 精品亚洲成国产av| 少妇被粗大猛烈的视频| 日韩一区二区三区影片| 亚洲国产精品一区二区三区在线| 好男人视频免费观看在线| 国产又色又爽无遮挡免| 大片电影免费在线观看免费| 高清在线视频一区二区三区| 日韩中字成人| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 精品卡一卡二卡四卡免费| 亚洲精华国产精华液的使用体验| 国产精品人妻久久久影院| 欧美日韩亚洲高清精品| 亚洲真实伦在线观看| 女性被躁到高潮视频| 黄色配什么色好看| 亚洲情色 制服丝袜| 亚洲av成人精品一区久久| 日韩伦理黄色片| 少妇 在线观看| 国产av码专区亚洲av| av在线播放精品| 国产av一区二区精品久久| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 国产一区二区三区av在线| 国产淫片久久久久久久久| 色5月婷婷丁香| 久久免费观看电影| 日韩av免费高清视频| 精品人妻熟女av久视频| 又爽又黄a免费视频| 日韩欧美 国产精品| 国产精品欧美亚洲77777| 王馨瑶露胸无遮挡在线观看| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 成人亚洲精品一区在线观看| av福利片在线| 一级二级三级毛片免费看| 亚洲欧美清纯卡通| 午夜影院在线不卡| 亚洲av日韩在线播放| 丰满乱子伦码专区| 51国产日韩欧美| 少妇的逼水好多| 国产亚洲5aaaaa淫片| 欧美另类一区| 男人爽女人下面视频在线观看| 国产亚洲av片在线观看秒播厂| 欧美区成人在线视频| 综合色丁香网| 久久av网站| 国产亚洲午夜精品一区二区久久| 亚洲三级黄色毛片| 日本色播在线视频| 男人狂女人下面高潮的视频| 国产一区二区三区av在线| 亚洲欧美日韩另类电影网站| 国产永久视频网站| 欧美日韩国产mv在线观看视频| 亚洲av综合色区一区| 国产精品蜜桃在线观看| 国产毛片在线视频| 免费av不卡在线播放| 国产成人精品无人区| 亚洲伊人久久精品综合| 日韩成人av中文字幕在线观看| 亚洲精品中文字幕在线视频 | 国产深夜福利视频在线观看| tube8黄色片| 国产av精品麻豆| 精品国产一区二区久久| 免费黄频网站在线观看国产| 人体艺术视频欧美日本| 日韩大片免费观看网站| 99久久精品国产国产毛片| 美女中出高潮动态图| av网站免费在线观看视频| 亚洲成人手机| 成人免费观看视频高清| 精品少妇内射三级| 黄色日韩在线| 精品卡一卡二卡四卡免费| 伦理电影大哥的女人| 日本av手机在线免费观看| 日韩在线高清观看一区二区三区| 亚洲精品国产色婷婷电影| 两个人免费观看高清视频 | 99久久精品一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 成人美女网站在线观看视频| 99久久中文字幕三级久久日本| 黑人巨大精品欧美一区二区蜜桃 | 国产女主播在线喷水免费视频网站| 午夜老司机福利剧场| 免费大片黄手机在线观看| 国产综合精华液| 尾随美女入室| 噜噜噜噜噜久久久久久91| 日韩视频在线欧美| 高清av免费在线| 18禁动态无遮挡网站| 免费高清在线观看视频在线观看| 久久99热6这里只有精品| 免费高清在线观看视频在线观看| 在线观看www视频免费| 日日啪夜夜撸| 国产片特级美女逼逼视频| 久久人妻熟女aⅴ| 久久国内精品自在自线图片| 国产成人精品无人区| 五月天丁香电影| 国产成人精品久久久久久| av不卡在线播放| 欧美三级亚洲精品| 久久久久网色| 亚洲av欧美aⅴ国产| 老司机亚洲免费影院| 精品熟女少妇av免费看| 麻豆成人午夜福利视频| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 日韩av不卡免费在线播放| 国产黄色视频一区二区在线观看| 蜜臀久久99精品久久宅男| 建设人人有责人人尽责人人享有的| 国产黄色免费在线视频| 建设人人有责人人尽责人人享有的| av在线观看视频网站免费| a级毛色黄片| 黄片无遮挡物在线观看| 香蕉精品网在线| 哪个播放器可以免费观看大片| 香蕉精品网在线| 久久精品国产自在天天线| 国产91av在线免费观看| 成人国产av品久久久| kizo精华| 日本爱情动作片www.在线观看| www.av在线官网国产| 嫩草影院入口| 永久网站在线| 亚洲欧美日韩另类电影网站| 成人黄色视频免费在线看| 午夜精品国产一区二区电影| 亚洲综合色惰| 男人添女人高潮全过程视频| 丰满乱子伦码专区| 精品少妇黑人巨大在线播放| 亚洲婷婷狠狠爱综合网| 男女免费视频国产| 97精品久久久久久久久久精品| 午夜日本视频在线| 极品教师在线视频| 成人国产麻豆网| 乱人伦中国视频| 亚洲av男天堂| 国产欧美日韩精品一区二区| 美女xxoo啪啪120秒动态图| 在线精品无人区一区二区三| 91精品国产国语对白视频| av福利片在线观看| 亚洲国产精品成人久久小说| 日产精品乱码卡一卡2卡三| 国产乱来视频区| 国产午夜精品久久久久久一区二区三区| 国产一级毛片在线| 亚洲av男天堂| 久热这里只有精品99| 日韩强制内射视频| 成年人免费黄色播放视频 | 伦精品一区二区三区| 久久青草综合色| av福利片在线观看| 国产精品久久久久久精品电影小说| 在线观看av片永久免费下载| 男人狂女人下面高潮的视频| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区 | 国产精品偷伦视频观看了| 欧美日韩亚洲高清精品| 日韩大片免费观看网站| 少妇被粗大的猛进出69影院 | 亚洲精品日韩av片在线观看| 国产日韩欧美亚洲二区| 夫妻午夜视频| 日本黄色日本黄色录像| 成人无遮挡网站| 亚洲欧美一区二区三区黑人 | 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 九色成人免费人妻av| 亚洲精品一二三| 久久久久久久精品精品| 国产高清不卡午夜福利| 国产精品久久久久久精品电影小说| 欧美亚洲 丝袜 人妻 在线| 亚洲精品中文字幕在线视频 | 国产精品嫩草影院av在线观看| 天天躁夜夜躁狠狠久久av| 成年人午夜在线观看视频| 国产极品天堂在线| 国产欧美日韩一区二区三区在线 | 一本大道久久a久久精品| 欧美精品人与动牲交sv欧美| 免费看av在线观看网站| 日韩制服骚丝袜av| 男女无遮挡免费网站观看| 免费大片18禁| 国产av精品麻豆| 水蜜桃什么品种好| 午夜av观看不卡| 亚洲国产精品国产精品| 少妇的逼好多水| 久久av网站| 好男人视频免费观看在线| 亚洲精品aⅴ在线观看| 久久99热6这里只有精品| 精品国产乱码久久久久久小说| 又黄又爽又刺激的免费视频.| 丰满乱子伦码专区| 成人免费观看视频高清| 午夜91福利影院| 午夜福利在线观看免费完整高清在| 一级毛片我不卡| 熟女电影av网| 精品一区在线观看国产| 成人无遮挡网站| 成人亚洲精品一区在线观看| 色视频在线一区二区三区| 91精品一卡2卡3卡4卡| 国产精品欧美亚洲77777| 欧美97在线视频| 日日撸夜夜添| 精品亚洲成国产av| 女人精品久久久久毛片| 另类精品久久| 自拍偷自拍亚洲精品老妇| 久久久久精品性色| 亚洲欧美一区二区三区黑人 | 韩国高清视频一区二区三区| 91午夜精品亚洲一区二区三区| 国产精品久久久久久精品电影小说| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 久久久久久久久久久免费av| 中文字幕亚洲精品专区| 日本午夜av视频| 特大巨黑吊av在线直播| 久久久久精品久久久久真实原创| 精品人妻一区二区三区麻豆| 一级毛片电影观看| 国产精品.久久久| 国产 精品1| 久久精品久久久久久噜噜老黄| 久久狼人影院| 91久久精品国产一区二区成人| 亚洲精品aⅴ在线观看| 国产精品蜜桃在线观看| 国产精品免费大片| 女性生殖器流出的白浆| 欧美日韩av久久| 又大又黄又爽视频免费| h视频一区二区三区| 91久久精品国产一区二区成人| 男的添女的下面高潮视频| 在线亚洲精品国产二区图片欧美 | 狂野欧美激情性bbbbbb| 一区二区av电影网| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 2022亚洲国产成人精品| 熟妇人妻不卡中文字幕| 97在线视频观看| 夫妻性生交免费视频一级片| 亚洲精品成人av观看孕妇| 精品国产一区二区久久| 亚洲精华国产精华液的使用体验| 人妻人人澡人人爽人人| 欧美精品一区二区免费开放| 国产高清三级在线| 一级毛片久久久久久久久女| 午夜福利在线观看免费完整高清在| 国产一级毛片在线| 亚洲内射少妇av| 看十八女毛片水多多多| 99热6这里只有精品| 日本vs欧美在线观看视频 | 欧美日韩国产mv在线观看视频| 日产精品乱码卡一卡2卡三| 在线观看人妻少妇| 久久久久久久精品精品| 91久久精品电影网| 久久精品国产亚洲av天美| 亚洲精品国产成人久久av| 欧美区成人在线视频| 欧美3d第一页| 高清av免费在线| 国产白丝娇喘喷水9色精品| 亚洲精品第二区| 少妇人妻久久综合中文| 国产伦在线观看视频一区| 久久国内精品自在自线图片| 精品久久久久久久久亚洲| 两个人免费观看高清视频 | 亚洲av中文av极速乱| 91成人精品电影| 久久久久久久精品精品| 国内少妇人妻偷人精品xxx网站| 国产在线一区二区三区精| 一级二级三级毛片免费看| 大话2 男鬼变身卡| 九九爱精品视频在线观看| 国产日韩欧美在线精品| 日韩成人av中文字幕在线观看| 成人毛片60女人毛片免费| 欧美成人午夜免费资源| 欧美性感艳星| 亚洲欧美精品自产自拍| 中文字幕制服av| 欧美少妇被猛烈插入视频| 精品国产露脸久久av麻豆| 国产69精品久久久久777片| 2018国产大陆天天弄谢| 久久久国产精品麻豆| 久久99热6这里只有精品| 啦啦啦在线观看免费高清www| 精品酒店卫生间| 亚洲精品成人av观看孕妇| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区二区三区四区免费观看| 秋霞在线观看毛片| 亚洲国产欧美在线一区| 黄色日韩在线| 欧美3d第一页| √禁漫天堂资源中文www| 老司机影院毛片| 妹子高潮喷水视频| √禁漫天堂资源中文www| 你懂的网址亚洲精品在线观看| 特大巨黑吊av在线直播| 国产高清国产精品国产三级| 亚洲电影在线观看av| 欧美 日韩 精品 国产| 国产亚洲最大av| 在线看a的网站| 久久精品国产亚洲av天美|