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

    楊樹PtROP家族基因的表達(dá)分析與功能預(yù)測(cè)

    2017-02-23 07:35:44王麗娟盧孟柱
    林業(yè)科學(xué)研究 2017年1期
    關(guān)鍵詞:擬南芥楊樹結(jié)構(gòu)域

    李 煜,張 進(jìn),王麗娟,盧孟柱*

    (1.林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091;2. 國(guó)家林業(yè)局林木培育重點(diǎn)實(shí)驗(yàn)室,中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091)

    楊樹PtROP家族基因的表達(dá)分析與功能預(yù)測(cè)

    李 煜1,2,張 進(jìn)1,2,王麗娟1,2,盧孟柱1,2*

    (1.林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室,中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091;2. 國(guó)家林業(yè)局林木培育重點(diǎn)實(shí)驗(yàn)室,中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所,北京 100091)

    [目的]為探究ROP基因在林木中的功能。[方法]本研究以楊樹為模式,在全基因組水平上對(duì)ROP基因的家族成員、基因結(jié)構(gòu)、保守結(jié)構(gòu)域、氨基酸序列相似性及表達(dá)模式進(jìn)行了分析。[結(jié)果]顯示,楊樹PtROP基因家族包含13個(gè)成員,不同成員間在進(jìn)化上相對(duì)保守,均存在與GTP結(jié)合與水解相關(guān)的結(jié)構(gòu)域。這些基因在不同組織、不同脅迫條件下的表達(dá)具有明顯差異,說明它們參與不同的生物學(xué)過程。通過構(gòu)建功能基因網(wǎng)絡(luò)(functional gene network for poplar)發(fā)現(xiàn)PtROP基因主要參與了信號(hào)轉(zhuǎn)導(dǎo)過程。[結(jié)論] 楊樹PtROP基因家族包含13個(gè)成員參與不同的生物學(xué)過程,可能主要參與了信號(hào)轉(zhuǎn)導(dǎo)過程。

    楊樹;ROP;基因家族;基因表達(dá)模式

    G蛋白是一類GTP結(jié)合蛋白,按照G蛋白分子量和亞基組成可分為異源三聚體G蛋白、小G蛋白和幾種特殊的GTP結(jié)合蛋白[1]。小G蛋白家族又稱為Ras超家族,根據(jù)其成員之間結(jié)構(gòu)和功能的相似性,將Ras超家族分為五個(gè)不同的家族,分別是Ras、Rab、Arf、Ran和Rho家族[2],其中Rho家族又可被分為CDC42、RAC和RHO三個(gè)亞家族[3-4]。這五個(gè)家族成員在細(xì)胞中參與了不同的生物學(xué)過程,例如Ras家族成員參與調(diào)控細(xì)胞增殖,Rab家族和Arf家族成員在介導(dǎo)不同膜細(xì)胞器之間的囊泡運(yùn)輸中發(fā)揮重要作用,Ran家族成員主要調(diào)控與細(xì)胞核相關(guān)的一些細(xì)胞活動(dòng),而Rho家族成員在細(xì)胞骨架重組及信號(hào)轉(zhuǎn)導(dǎo)中行使重要功能[5]。

    與動(dòng)物或真菌Ras超家族不同,植物沒有Ras家族,而植物中Rho形成了獨(dú)特的一類,稱為ROP(Rho-related proteins from plants)[5]。自1993年首次從豌豆中發(fā)現(xiàn)Rho相似性GTPase以來,已在擬南芥、水稻與玉米等多種植物中進(jìn)行了鑒定[6]。ROP主要在與GTP結(jié)合的活性形式和與GDP結(jié)合的非活性形式之間循環(huán),該循環(huán)過程需要多種蛋白的參與,如鳥苷酸交換因子(guanine nucleotide exchange factor,GEF)、GTPase激活蛋白(GTPase activating protein,GAP)、鳥苷酸解離抑制因子(GDP dissociation inhibitor,GDI)等。一旦受體感知到上游信號(hào),鳥苷酸交換因子就催化ROP從GDP結(jié)合的非活性形式轉(zhuǎn)換為GTP結(jié)合的活性形式,具有活性的ROP蛋白與下游效應(yīng)蛋白結(jié)合發(fā)揮其生理功能[7]。ROP蛋白的結(jié)構(gòu)極為保守,主要包含與鳥苷酸結(jié)合和水解有關(guān)的5個(gè)高度保守結(jié)構(gòu)域、1個(gè)效應(yīng)因子結(jié)合區(qū)、1個(gè)插入序列和C端高度可變區(qū)。GTP結(jié)合結(jié)構(gòu)域和效應(yīng)因子結(jié)合區(qū)決定ROP蛋白的功能,插入序列的功能目前尚不清楚,C端高度可變區(qū)決定ROP蛋白的亞細(xì)胞定位[8]。

    ROP在植物發(fā)育與抗逆等生物學(xué)過程中具有重要的功能。在葉表皮鋪板細(xì)胞中,擬南芥ROP2調(diào)節(jié)PIN1的內(nèi)吞過程[9]。在擬南芥胚胎發(fā)育階段,ROP3在維持PIN1和PIN3在細(xì)胞膜上的極性分布方面具有重要作用[10]。擬南芥中多種ROP基因在花粉管尖端表達(dá),其中ROP1對(duì)花粉管的極性生長(zhǎng)起主要作用[11]。ROP還參與根毛的發(fā)育,超表達(dá)擬南芥ROP2,能增加根毛數(shù)量和密度,而AtROP4控制根毛的發(fā)生和尖端的伸長(zhǎng)[12]。此外,擬南芥保衛(wèi)細(xì)胞中的ROP2能夠被光激活來抑制光誘導(dǎo)的氣孔張開[13]。目前,在擬南芥基因組中共發(fā)現(xiàn)了10個(gè)ROP基因成員,水稻和玉米分別發(fā)現(xiàn)7個(gè)和9個(gè)[6]。雖然一些研究都表明ROPs廣泛參與了植物對(duì)激素,環(huán)境等信號(hào)的響應(yīng)及植物生長(zhǎng)發(fā)育過程,但各個(gè)ROP成員的功能還有待進(jìn)一步解析。楊樹全基因組序列雖已公布,但對(duì)ROP基因家族成員尚未進(jìn)行系統(tǒng)分析[14]。

    本研究利用擬南芥ROP氨基酸序列,在楊樹基因組數(shù)據(jù)庫(kù)中進(jìn)行同源比對(duì),共鑒定出13個(gè)PtROP基因序列。在此基礎(chǔ)上,對(duì)楊樹PtROP基因進(jìn)行了系統(tǒng)進(jìn)化、基因結(jié)構(gòu)及保守結(jié)構(gòu)域分析,并對(duì)其在不同組織和不同脅迫條件下的表達(dá)情況進(jìn)行了分析,構(gòu)建了其基因調(diào)控網(wǎng)絡(luò)。研究結(jié)果為將來系統(tǒng)解析楊樹ROP基因的功能奠定了基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    實(shí)驗(yàn)所用銀腺楊(P.alba×P.glandulosa)無性系84K,由中國(guó)林業(yè)科學(xué)研究院林木遺傳育種國(guó)家重點(diǎn)實(shí)驗(yàn)室保存。楊樹幼苗培養(yǎng)于人工培養(yǎng)箱,培養(yǎng)條件為16 h光照/8 h黑暗,晝夜溫度為23℃/20℃。分別取楊樹幼株(株高約80 cm)的頂端、幼葉、成熟葉、莖頂端、莖下部及根,分別提取總RNA。

    1.2 總RNA的提取

    總RNA提取參見QIAGEN RNeasy Plant Kit RNA描述的方法進(jìn)行。

    1.3 轉(zhuǎn)錄組測(cè)序

    用帶有Oligo(dT)的磁珠富集樣品總RNA中的mRNA,將富集的mRNA片段化后用隨機(jī)引物合成第一鏈cDNA,隨后合成雙鏈cDNA,再對(duì)雙鏈cDNA進(jìn)行純化,末端修復(fù)和添加測(cè)序接頭,最后進(jìn)行PCR擴(kuò)增,建好的測(cè)序文庫(kù)用Illumina HiSeqTM2000測(cè)序平臺(tái)進(jìn)行pair-end 100 bp測(cè)序,每個(gè)樣品分別獲得2G clean data。測(cè)序原始數(shù)據(jù)去除接頭及低質(zhì)量reads,獲得的clean reads利用TopHat2比對(duì)到毛果楊基因組上(http://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Ptrichocarpa,release 3.0)。標(biāo)準(zhǔn)化的FPKM值反映基因的表達(dá)水平。

    1.4 數(shù)據(jù)來源

    擬南芥ROP基因及其編碼的蛋白質(zhì)序列來源于TAIR(http://www.arabidopsis.org/index.jsp)數(shù)據(jù)庫(kù)。以擬南芥ROP蛋白序列為搜索對(duì)象,通過blastp程序在毛果楊基因組數(shù)據(jù)庫(kù)(http://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Ptrichocarpa,release 3.0)中搜索相似的蛋白序列,選取E≤10-10的序列作為候選序列,獲得的毛果楊ROP基因命名為PtROP1-13。

    1.5 系統(tǒng)進(jìn)化樹分析

    將擬南芥和楊樹PtROP氨基酸序列保存為FASTA格式,用Clustal X[15]對(duì)其氨基酸序列進(jìn)行多序列比對(duì),分別利用MEGA5.0[16]中Maximum-Likelihood(ML)法和Neighbor-Joining(NJ)法構(gòu)建系統(tǒng)進(jìn)化樹。

    1.6 基因結(jié)構(gòu)分析

    運(yùn)用在楊樹基因組數(shù)據(jù)庫(kù)中下載的基因組DNA序列及CDS序列使用GSDS[17]程序(http://gsds.cbi.pku.edu.cn/)對(duì)楊樹PtROP基因的外顯子和內(nèi)含子進(jìn)行分析。楊樹PtROP基因組DNA序列及CDS序列下載自Phytozome數(shù)據(jù)庫(kù)(http://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Ptrichocarpa, release 3.0)。

    1.7 基因表達(dá)分析

    分析楊樹不同脅迫條件下的PtROP基因成員的表達(dá)模式所用的基因芯片數(shù)據(jù)來自NCBI Gene Expression Omnibus(GEO)[18]數(shù)據(jù)庫(kù)。

    1.8 楊樹PtROP基因功能網(wǎng)絡(luò)分析

    利用網(wǎng)站(http://10.13.200.25/PoplarGene/index.php)對(duì)楊樹ROP基因進(jìn)行功能網(wǎng)絡(luò)分析。

    2 結(jié)果與分析

    2.1 楊樹PtROP基因鑒定及系統(tǒng)進(jìn)化分析

    表1 楊樹PtROP基因家族成員基本信息Table1 Basic information of Populus PtROP gene

    A.最大似然法(ML)構(gòu)建的系統(tǒng)進(jìn)化樹;B.鄰接法(NJ)構(gòu)建的系統(tǒng)進(jìn)化樹A. Phylogenetic tree was constructed by the Maximum Likelihood method; B. Phylogenetic tree was constructed by the neighbor-joining method圖1 楊樹和擬南芥ROP基因家族系統(tǒng)進(jìn)化分析Fig. 1 Phylogenetic relationship between Populus and Arabidopsis ROP families

    2.2PtROP基因結(jié)構(gòu)及保守結(jié)構(gòu)域分析

    利用MEME結(jié)構(gòu)域搜索工具在楊樹PtROP中檢測(cè)到15個(gè)保守的結(jié)構(gòu)域。結(jié)果表明, PtROP10中不存在motif 2,PtROP5不具有motif 5,而其它ROP均具有motif 1-7。Motif 1、2、3及5代表了ROP蛋白與GTP結(jié)合的結(jié)構(gòu)域(圖4),這四個(gè)結(jié)構(gòu)域中都存在與GTP結(jié)合的保守結(jié)合位點(diǎn)。PtROP5雖然不具有motif 5,但它具有其它與GTP結(jié)合的保守結(jié)合位點(diǎn)。對(duì)楊樹PtROP蛋白motif 2中第20位的蘇氨酸、motif 3中第10位的天冬氨酸分別進(jìn)行點(diǎn)突變,能夠改變ROP蛋白與GTP的親和性,使其成為顯性失活形式。如果對(duì)楊樹PtROP蛋白motif 1中將第6

    位的谷氨酰胺突變?yōu)榱涟彼?,ROP蛋白就成為了持續(xù)活性型形式(圖4)。

    圖2 楊樹PtROP基因結(jié)構(gòu)Fig. 2 Gene structures of PtROP genes

    2.3PtROP基因旁系同源基因?qū)Φ膹?fù)制

    通過比較其系統(tǒng)進(jìn)化關(guān)系、基因結(jié)構(gòu)及氨基酸序列相似性,在楊樹PtROP基因家族中共鑒定得到4對(duì)旁系同源基因?qū)?,分別是PtROP1/PtROP4、PtROP2/PtROP6、PtROP8/PtROP11及PtROP9/PtROP13,這些基因?qū)χg在氨基酸序列、內(nèi)含子的插入大小及插入位點(diǎn)都極為保守。PtROP1/PtOP4、PtROP2/PtROP6、PtROP8/PtROP11、PtROP9/PtROP13所在染色體片段區(qū)域的大部分基因都是同源基因,因此這四個(gè)基因?qū)赡苁怯捎谌旧w復(fù)制產(chǎn)生的(圖5)。計(jì)算這四個(gè)基因?qū)Ψ峭x替換率和同義替換率的比值,發(fā)現(xiàn)它們都小于0.4,因此基因?qū)χ械膬蓚€(gè)旁系同源基因復(fù)制后經(jīng)歷了純化選擇在功能上產(chǎn)生了一定程度的差異(表2)。根據(jù)楊樹擴(kuò)張度λ=9.1×10-9,估算出PtROP旁系同源基因?qū)赡苁窃?.74到17.30百萬年前產(chǎn)生的,其中PtROP1/PtROP4產(chǎn)生的最早,約在17.30百萬年前;PtROP8/PtROP11產(chǎn)生的最晚,約在8.74百萬年前。而PtROP2/PtROP6、PtROP9/PtROP13分別是在13.98百萬年前和11.77百萬年前產(chǎn)生的,這個(gè)時(shí)期與楊樹最近一次在13百萬年前發(fā)生的大規(guī)模片段復(fù)制時(shí)間大致相符[14]。

    圖3 楊樹與擬南芥ROP氨基酸序列相似性Fig. 3 Sequence similarites between Populus and Arabidopsis ROP proteins

    圖4 楊樹PtROP蛋白保守結(jié)構(gòu)域分布Fig. 4 Distribution of conserved motifs in PtROP proteins

    圖5 楊樹PtROP旁系同源基因?qū)ig. 5 PtROP paralogous gene pairs in P trichocarpa ROP gene family

    基因1基因2非同義替換率同義替換率非同義替換率/同義替換率時(shí)期(百萬年前)PtROP1PtROP40.0740.3150.23617.30PtROP2PtROP60.0130.2550.05213.98PtROP8PtROP110.0200.1590.1258.74PtROP9PtROP130.0200.2140.09411.77

    2.4 在不同組織、不同脅迫條件下的表達(dá)

    為了探究楊樹PtROP基因可能參與的生物學(xué)過程,利用課題組的RNA-seq數(shù)據(jù)和已公布的楊樹基因芯片數(shù)據(jù),分析了楊樹PtROP基因在不同組織及不同脅迫條件下的表達(dá)模式。與幼葉、老葉、幼莖、木質(zhì)化莖及根相比,楊樹PtROP1、PtROP2、PtROP8、PtROP11及PtROP12在頂點(diǎn)的表達(dá)水平較高,表明它們可能參與了頂端發(fā)育過程。PtROP7、PtROP9及ROP13特異性地在已木質(zhì)化的莖中表達(dá),說明這些基因可能參與楊樹次生生長(zhǎng)。PtROP1、PtROP8、PtROP9及PtROP12在老葉中的表達(dá)量最低,表明它們可能不參與光合作用及次生物質(zhì)的代謝過程(圖6)。

    圖6 楊樹PtROP基因在不同組織中的表達(dá)模式Fig.6 Expression analyses of PtROPs in various tissues

    在應(yīng)拉木和應(yīng)力木中PtROP7、PtROP9、PtROP10及PtROP13表達(dá)都顯著下調(diào),且PtROP7、PtROP9及PtROP13在次生莖中特異性高表達(dá),表明這三個(gè)基因的表達(dá)受外力的調(diào)控。對(duì)高溫、低溫、干旱和鹽脅迫條件下PtROP基因的表達(dá)分析發(fā)現(xiàn),13個(gè)基因中僅有PtROP7在低溫條件下下調(diào)表達(dá),表明大部分PtROP基因可能不直接對(duì)高、低溫、干旱和鹽等脅迫響應(yīng),或者其在脅迫條件下的穩(wěn)定表達(dá)對(duì)細(xì)胞功能的維持具有重要作用(圖7)。

    2.5 楊樹PtROP基因功能基因網(wǎng)絡(luò)分析

    為了進(jìn)一步解析PtROP基因的生物學(xué)功能,構(gòu)建了PtROP的功能基因網(wǎng)絡(luò)。發(fā)現(xiàn)與PtROP基因緊密聯(lián)系的基因主要分為四類,即Rho GDP-解離抑制活性、信號(hào)轉(zhuǎn)導(dǎo)、tRNA修飾及Actin細(xì)胞骨架(表3)。其中,在整個(gè)網(wǎng)絡(luò)中大約50%的基因與信號(hào)轉(zhuǎn)導(dǎo)相關(guān),并且這些基因與楊樹PtROP基因的相關(guān)性都比較大,這與擬南芥中ROP參與了信號(hào)轉(zhuǎn)導(dǎo)的結(jié)果相符合,說明PtROP基因具有功能保守性。13個(gè)基因中,與PtROP2和PtROP3相關(guān)聯(lián)的基因最多,表明二者在整個(gè)調(diào)控網(wǎng)絡(luò)中發(fā)揮主要作用。在網(wǎng)絡(luò)中,與tRNA修飾相關(guān)的基因占50%以上,但僅有PtROP4與之相關(guān)聯(lián),并且它們之間的相關(guān)性較小,可能PtROP4廣泛地參與到tRNA修飾中,但不在其中起主要作用,另外,PtROP4基因在13個(gè)PtROP基因中僅與PtROP1存在關(guān)聯(lián),顯示楊樹PtROP4可能在進(jìn)化過程中出現(xiàn)了功能分化(圖8)。

    圖7 楊樹PtROP基因不同脅迫下的表達(dá)(PtROPs在不同脅迫下的變化倍數(shù)的log2轉(zhuǎn)換值)Fig.7 Expression analyses of PtROPs under various stresses

    表3 楊樹PtROP功能基因網(wǎng)絡(luò)中基因相關(guān)信息Table 3 Genes information in the PtROP gene functional network

    圖8 楊樹PtROP基因調(diào)控網(wǎng)絡(luò)Fig. 8 PtROP gene functional network

    3 討論

    楊樹是全球廣泛分布的重要經(jīng)濟(jì)樹種,具有生長(zhǎng)迅速,生物量大、輪伐期短等特點(diǎn),隨著楊樹全基因組測(cè)序的完成,楊樹已成為研究多年生林木生長(zhǎng)發(fā)育的模式樹木。在不同物種中對(duì)ROP基因的研究發(fā)現(xiàn)ROP基因廣泛地參與了植物的生長(zhǎng)發(fā)育和抗逆過程,而且ROP基因在初生細(xì)胞壁和次生細(xì)胞壁的沉積過程中具有重要作用[20]。因此對(duì)于楊樹ROP基因的研究,不僅在揭示木本植物發(fā)育和抗逆機(jī)理上具有重要意義,而且具有重要的應(yīng)用價(jià)值。

    本研究利用已知的擬南芥ROP氨基酸序列在楊樹全基因組數(shù)據(jù)庫(kù)中進(jìn)行相似性序列搜索,共檢測(cè)到13個(gè)PtROP基因成員,其中包含四個(gè)旁系同源基因?qū)?;進(jìn)化分析表明,這些基因可能來自全基因組染色體復(fù)制。楊樹PtROP基因結(jié)構(gòu)分析、系統(tǒng)進(jìn)化分析及氨基酸序列相似性分析表明,楊樹PtROP基因家族成員在結(jié)構(gòu)與氨基酸序列上極為保守。楊樹PtROP基因與擬南芥相同,分為四類,它們相對(duì)應(yīng)的類別,氨基酸序列相似度很高。對(duì)楊樹PtROP的motif研究表明,楊樹和擬南芥ROP一樣,都具有相同的與GTP結(jié)合的保守結(jié)構(gòu)域。這些研究結(jié)果預(yù)示著ROPs在楊樹中具有類似的、保守的功能。

    對(duì)楊樹PtROP基因家族成員在不同組織、不同脅迫條件下的表達(dá)情況進(jìn)行了分析,發(fā)現(xiàn)不同的PtROP基因在不同的組織中表達(dá)存在差異,表明它們可能參與了不同的生物學(xué)過程。研究發(fā)現(xiàn)擬南芥ROP11影響次生細(xì)胞壁的沉積,在超表達(dá)持續(xù)活性型ROP11擬南芥中,正在分化的木質(zhì)部細(xì)胞次生壁上沒有紋孔的形成,而ROP效應(yīng)因子ROPGEF4的缺失導(dǎo)致次生壁上紋孔數(shù)量減少[21]。桉樹ROP1也能夠影響次生細(xì)胞壁的沉積,其在形成層和正在分化的木質(zhì)部細(xì)胞中都有很高的表達(dá)水平,在擬南芥中超表達(dá)持續(xù)活性型的桉樹ROP1基因抑制次生木質(zhì)部的木質(zhì)化和木葡聚糖沉積[22]。本研究發(fā)現(xiàn)楊樹PtROP7、PtROP9與PtROP13特異性地在莖的次生長(zhǎng)生中高豐度表達(dá),表明這三個(gè)基因可能參與了楊樹次生生長(zhǎng)過程。PtROP7、PtROP9與PtROP13同屬于楊樹PtROP基因家族的第三類,可能存在功能上的冗余,以保障在莖的次生維管發(fā)育中發(fā)揮重要作用。

    已有研究表明,植物ROP蛋白參與了多種逆境脅迫。例如,油菜ROP5和ROP9在高溫脅迫后表達(dá)發(fā)生明顯變化[23],玉露桃冷害脅迫后,果實(shí)中的ROP基因表達(dá)量顯著降低[24]。本研究通過對(duì)高溫、低溫、干旱與鹽等脅迫條件下楊樹PtROP的表達(dá)分析發(fā)現(xiàn),PtROP基因?qū)@些脅迫沒有產(chǎn)生顯著響應(yīng),可能是由于楊樹PtROP基因功能的特化,趨于基本維持細(xì)胞正常的信息與物質(zhì)的流動(dòng)。通過楊樹PtROP的功能基因網(wǎng)絡(luò)分析,發(fā)現(xiàn)楊樹PtROP基因主要參與了信號(hào)轉(zhuǎn)導(dǎo)過程,這與上述推測(cè)是一致的[25]。

    4 結(jié)論

    本研究以楊樹為模式,在全基因組水平上對(duì)ROP基因的家族成員、基因結(jié)構(gòu)、保守結(jié)構(gòu)域、氨基酸序列相似性及表達(dá)模式進(jìn)行了分析。發(fā)現(xiàn)楊樹中共有13個(gè)PtROP基因,它們的基因結(jié)構(gòu)相似,功能結(jié)構(gòu)域是保守的。根據(jù)氨基酸序列相似性及C端高度可變區(qū)結(jié)構(gòu),楊樹PtROP被分為四類:第一類(PtROP10)、第二類(PtROP8、PtROP11和PtROP12)、第三類(PtROP7、PtROP9和PtROP13)及第四類(PtROP1-6)。對(duì)PtROP基因的復(fù)制情況進(jìn)行分析,發(fā)現(xiàn)PtROP基因家族中存在四對(duì)旁系同源基因?qū)?PtROP1/PtROP4、PtROP2/PtROP6、PtROP8/PtROP11及PtROP9/PtROP13),它們都是通過全基因組復(fù)制產(chǎn)生的。四個(gè)基因?qū)Χ冀?jīng)歷了純化選擇。PtROP旁系同源基因?qū)赡苁窃?.74到17.30百萬年前產(chǎn)生的。對(duì)PtROP基因家族成員在不同組織和不同脅迫響應(yīng)進(jìn)行分析,發(fā)現(xiàn)PtROP基因家族成員表達(dá)模式具有明顯差異。通過構(gòu)建PtROP功能基因網(wǎng)絡(luò),發(fā)現(xiàn)PtROP可能參與了信號(hào)轉(zhuǎn)導(dǎo)、tRNA修飾及細(xì)胞骨架形成等過程。

    本研究系統(tǒng)分析了楊樹PtROP家族的成員及其進(jìn)化關(guān)系,研究PtROP基因在木本植物發(fā)育過程中的表達(dá)特性和可能的功能,為進(jìn)一步研究其在木本植物特別是次生生長(zhǎng)過程中的作用奠定了基礎(chǔ)。

    [1] Assmann S M. Heterotrimeric and unconventional GTP binding proteins in plant cell signaling[J]. The Plant Cell, 2002, 149(suppl 1): S355-S373.

    [2] Kahn R A, Der C J, Bokoch G M. The ras superfamily of GTP-binding proteins: guidelines on nomenclature[J]. The FASEB Journal, 1992, 6(8): 2512-2513.

    [3] Chant J, Stowers L. GTPase cascades choreographing cellular behavior: movement, morphogenesis, and more[J]. Cell, 1995, 81(1): 1-4.[4] Hall A. Rho GTPases and the actin cytoskeleton[J]. Science, 1998, 279(5350): 509-514.

    [5] Vernoud V, Horton A C, Yang Z,etal. Analysis of the small GTPase gene superfamily ofArabidopsis[J]. Plant Physiology, 2003, 131(3): 1191-1208.

    [6] Pandey G K, Sharma M, Pandey A,etal. GTPases: Versatile Regulators of Signal Transduction in Plants[M]. Springer, 2015.

    [7] Ren H, Lin D. ROP GTPase Regulation of Auxin Transport inArabidopsis[J]. Molecular Plant, 2015, 8(2): 193-195.

    [8] Berken A, Wittinghofer A. Structure and function of Rho-type molecular switches in plants[J]. Plant Physiology and Biochemistry, 2008, 46(3): 380-393.

    [9] Nagawa S, Xu T, Lin D,etal. ROP GTPase-dependent actin microfilaments promote PIN1 polarization by localized inhibition of clathrin-dependent endocytosis[J]. PLoS Biol, 2012, 10(4): e1001299.

    [10] Huang J, Liu H, Chen M,etal. ROP3 GTPase contributes to polar auxin transport and auxin responses and is important for embryogenesis and seedling growth inArabidopsis[J]. The Plant Cell, 2014, 26(9): 3501-3518.

    [11] Li H, Wu G, Ware D,etal. Arabidopsis Rho-related GTPases: differential gene expression in pollen and polar localization in fission yeast[J]. Plant Physiology, 1998, 118(2): 407-417.

    [12] Molendijk A J, Bischoff F, Rajendrakumar C S V,etal.Arabidopsisthaliana Rop GTPases are localized to tips of root hairs and control polar growth[J]. The EMBO Journal, 2001, 20(11): 2779-2788.[13] Jeon B W, Hwang J U, Hwang Y,etal. TheArabidopsissmall G protein ROP2 is activated by light in guard cells and inhibits light-induced stomatal opening[J]. The Plant Cell, 2008, 20(1): 75-87.[14] Tuskan G A, Difazio S, Jansson S,etal. The genome of black cottonwood,Populustrichocarpa(Torr. & Gray)[J].Science, 2006, 313(5793): 1596-1604.

    [15] Thompson J D, Gibson T J, Plewniak F,etal. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Research, 1997, 25(24): 4876-4882.

    [16] Tamura K, Peterson D, Peterson N,etal. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Molecular Biology and Evolution, 2011, 28(10): 2731-2739.

    [17] Guo A Y, Zhu Q H, Chen X,etal. [GSDS: a gene structure display server[J]. Yi Chuan= Hereditas/Zhongguo Yi Chuan Xue Hui Bian Ji, 2007, 29(8): 1023-1026.

    [18] Barrett T, Edgar R. [19] Gene Expression Omnibus: Microarray Data Storage, Submission, Retrieval, and Analysis[J]. Methods in enzymology, 2006, 411: 352-369.

    [19]Zheng Z L, Yang Z. The Rop GTPase switch turns on polar growth in pollen[J]. Trends in plant science, 2000, 5(7): 298-303.

    [20] Oda Y, Fukuda H. Emerging roles of small GTPases in secondary cell wall development[J]. Frontiers in Plant Science, 2014, 5.

    [21] Oda Y, Fukuda H. Initiation of cell wall pattern by a Rho-and microtubule-driven symmetry breaking[J]. Science, 2012, 337(6100): 1333-1336.

    [22] Foucart C, Jauneau A, Gion J M,etal. Overexpression ofEgROP1, a Eucalyptus vascular-expressed Rac-like small GTPase, affects secondary xylem formation inArabidopsisthaliana[J]. New Phytologist, 2009, 183(4): 1014-1029.

    [23] Chan J, Pauls K P. Brassica napus Rop GTPases and their expression in microspore cultures[J]. Planta, 2007, 225(2): 469-484.[24] Jin W, Xu C, Li X,etal. Expression ofROP/RACGTPasegenes in postharvest loquat fruit in association with senescence and cold regulated lignification[J]. Postharvest Biology and Technology, 2009, 54(1): 9-14.

    [25] Miyawaki K N, Yang Z. Extracellular signals and receptor-like kinases regulating ROP GTPases in plants[J]. Frontiers in Plant Science, 2014, 5.

    (責(zé)任編輯:張 研)

    Expression and Functional Analysis ofROPGene Family inPopulus

    LIYu1,2,ZHANGJin1,2,WANGLi-juan1,2,LUMeng-zhu1,2

    (1.State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China;2.Key Laboratory of Tree Breeding and Cultivation, State Forestry Administration, Research Institute of Forestry, Chinese Academy of Forestry,Beijing 100091, China)

    [Objective]To explore the putative function ofROPsin forest trees. [Method]A genome-wide analysis ofROPswas performed, including the phylogeny, gene structure, conserved motifs, ROP amino sequences similarity betweenPopulusandArabidopsisROP proteins and the expression patterns usingPopulusas a model. [Result]The results showed that there were 13 members ofPtROPgenes inPopulusand thePtROPgenes were conserved during evolution of species, all of them containing GTP binding and hydrolysis-related domains. Furthermore, the expression profiles revealed thatPopulusPtROPhad distinct expression pattern across different tissues and different stress conditions, suggesting thatPtROPgenes were involved in different biological processes. In addition, the analysis onPtROPgene functional network was conducted. It is predicted thatPopulusPtROPgenes were mainly involved in signal transduction.[Conclude]In this study, systematically bioinformatics analysis onPtROPwere conducted which laid the foundation for the further exploration of poplarPtROPfunctions.

    Populus; ROP; gene family; gene expression pattern

    10.13275/j.cnki.lykxyj.2017.01.001

    2016-02-01

    國(guó)家自然科學(xué)基金面上項(xiàng)目(31270699)。

    李 煜,男,在讀碩士。主要方向:林木極性生長(zhǎng)調(diào)控。E-mail:linxueliyu@163.com
    * 通訊作者:研究員,博士生導(dǎo)師,從事分子生物學(xué)研究。E-mail:lumz@caf.ac.cn

    S792.11

    A

    1001-1498(2017)01-0001-09

    猜你喜歡
    擬南芥楊樹結(jié)構(gòu)域
    擬南芥:活得粗糙,才讓我有了上太空的資格
    蛋白質(zhì)結(jié)構(gòu)域劃分方法及在線服務(wù)綜述
    跟蹤導(dǎo)練(五)
    尿黑酸對(duì)擬南芥酪氨酸降解缺陷突變體sscd1的影響
    口算本失蹤記
    楊樹山漫畫作品欣賞
    喜劇世界(2017年10期)2017-06-01 12:39:29
    兩種LED光源作為擬南芥生長(zhǎng)光源的應(yīng)用探究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    重組綠豆BBI(6-33)結(jié)構(gòu)域的抗腫瘤作用分析
    組蛋白甲基化酶Set2片段調(diào)控SET結(jié)構(gòu)域催化活性的探討
    国产精品香港三级国产av潘金莲| 国产免费男女视频| 久久热在线av| 中文字幕av电影在线播放| 国产亚洲精品av在线| 91九色精品人成在线观看| 亚洲午夜理论影院| 美国免费a级毛片| 激情在线观看视频在线高清| 国产亚洲精品久久久久久毛片| 变态另类成人亚洲欧美熟女| 亚洲电影在线观看av| 国产精品综合久久久久久久免费| av中文乱码字幕在线| 在线永久观看黄色视频| 国产不卡一卡二| 嫩草影视91久久| 久久久久久免费高清国产稀缺| 久久久久久九九精品二区国产 | 波多野结衣巨乳人妻| 久久精品国产综合久久久| 欧美乱色亚洲激情| 两性午夜刺激爽爽歪歪视频在线观看 | 1024手机看黄色片| 好男人在线观看高清免费视频 | 午夜福利在线在线| 黄色视频,在线免费观看| 免费在线观看黄色视频的| 久久久久久久久免费视频了| 午夜视频精品福利| 老汉色av国产亚洲站长工具| 国产av又大| 久久人妻av系列| 久久国产精品影院| 免费在线观看完整版高清| xxxwww97欧美| 无限看片的www在线观看| 婷婷精品国产亚洲av在线| 啦啦啦 在线观看视频| tocl精华| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 在线永久观看黄色视频| 午夜精品在线福利| 久久 成人 亚洲| 日韩欧美 国产精品| 国产精品二区激情视频| 国产极品粉嫩免费观看在线| 亚洲一区二区三区不卡视频| 国产男靠女视频免费网站| 免费av毛片视频| 久久久国产精品麻豆| 91成年电影在线观看| 听说在线观看完整版免费高清| 亚洲国产欧美一区二区综合| 男人的好看免费观看在线视频 | 精品欧美一区二区三区在线| 久久久久精品国产欧美久久久| 日日夜夜操网爽| 又黄又爽又免费观看的视频| 亚洲国产高清在线一区二区三 | 日本免费a在线| 国产色视频综合| 久久草成人影院| 夜夜躁狠狠躁天天躁| 国产一级毛片七仙女欲春2 | 国产av在哪里看| 欧美性猛交黑人性爽| 女警被强在线播放| 99国产精品99久久久久| 97人妻精品一区二区三区麻豆 | 亚洲av片天天在线观看| 亚洲第一av免费看| 亚洲黑人精品在线| 久久人妻av系列| 欧美日韩黄片免| АⅤ资源中文在线天堂| 桃色一区二区三区在线观看| 国产精品一区二区精品视频观看| 天堂动漫精品| 欧美性猛交╳xxx乱大交人| 免费搜索国产男女视频| 夜夜躁狠狠躁天天躁| 一进一出抽搐gif免费好疼| 日韩三级视频一区二区三区| 可以在线观看的亚洲视频| 成人精品一区二区免费| 亚洲自拍偷在线| 51午夜福利影视在线观看| 久久久久国产一级毛片高清牌| 欧美国产精品va在线观看不卡| 啦啦啦免费观看视频1| 69av精品久久久久久| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 亚洲男人天堂网一区| 最近在线观看免费完整版| 国产av一区在线观看免费| 99国产精品99久久久久| 欧美绝顶高潮抽搐喷水| 熟女电影av网| 久久精品人妻少妇| www日本黄色视频网| 天天躁夜夜躁狠狠躁躁| 日本一本二区三区精品| 亚洲av片天天在线观看| 在线观看免费视频日本深夜| 午夜视频精品福利| 久久亚洲精品不卡| 亚洲熟妇中文字幕五十中出| 嫩草影院精品99| 色播亚洲综合网| 亚洲人成伊人成综合网2020| 亚洲欧洲精品一区二区精品久久久| 亚洲精品av麻豆狂野| 欧美激情久久久久久爽电影| 亚洲精品在线美女| 久久久久亚洲av毛片大全| 久久亚洲精品不卡| 黑丝袜美女国产一区| 久久久久九九精品影院| 国产成人一区二区三区免费视频网站| 午夜久久久久精精品| 成人永久免费在线观看视频| 免费在线观看影片大全网站| 欧美黑人精品巨大| 极品教师在线免费播放| 欧美乱妇无乱码| 一卡2卡三卡四卡精品乱码亚洲| 国产成人精品久久二区二区免费| 国产精品98久久久久久宅男小说| 成人av一区二区三区在线看| 国产精品香港三级国产av潘金莲| 国产精品爽爽va在线观看网站 | 国产极品粉嫩免费观看在线| 一本精品99久久精品77| 欧美乱色亚洲激情| 久久久久久九九精品二区国产 | 国产一卡二卡三卡精品| 国产黄色小视频在线观看| 欧美色视频一区免费| 中文在线观看免费www的网站 | 成年人黄色毛片网站| 精品久久蜜臀av无| 色播在线永久视频| 亚洲精品国产区一区二| 午夜a级毛片| 香蕉丝袜av| 一进一出抽搐动态| 日韩大尺度精品在线看网址| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 久久香蕉国产精品| 大型av网站在线播放| 热99re8久久精品国产| 精品熟女少妇八av免费久了| 美女扒开内裤让男人捅视频| 亚洲精品国产一区二区精华液| 成人一区二区视频在线观看| 欧美成人午夜精品| 亚洲国产看品久久| 免费观看人在逋| bbb黄色大片| 可以在线观看毛片的网站| 亚洲中文字幕日韩| 亚洲精品久久成人aⅴ小说| 精品国产乱码久久久久久男人| 18禁黄网站禁片午夜丰满| 亚洲人成网站高清观看| 一级黄色大片毛片| 亚洲国产精品999在线| 亚洲av成人av| 嫁个100分男人电影在线观看| www.精华液| 国产免费男女视频| 女人高潮潮喷娇喘18禁视频| 成人免费观看视频高清| 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 黄色成人免费大全| 一区福利在线观看| 大型av网站在线播放| 少妇的丰满在线观看| 久久久久免费精品人妻一区二区 | 91成年电影在线观看| 亚洲人成电影免费在线| 国产成人av激情在线播放| 在线观看免费视频日本深夜| 亚洲精品久久成人aⅴ小说| 69av精品久久久久久| 老司机福利观看| 免费在线观看视频国产中文字幕亚洲| av天堂在线播放| 视频区欧美日本亚洲| 欧美人与性动交α欧美精品济南到| 精品第一国产精品| 在线观看免费视频日本深夜| 国产高清激情床上av| 69av精品久久久久久| 色av中文字幕| 久久久久国产精品人妻aⅴ院| 日韩成人在线观看一区二区三区| 国产精品自产拍在线观看55亚洲| 一级作爱视频免费观看| 看免费av毛片| 一区二区三区高清视频在线| 亚洲熟妇熟女久久| 欧美色视频一区免费| 久久青草综合色| 日韩欧美 国产精品| 欧美激情极品国产一区二区三区| 久9热在线精品视频| 国产99久久九九免费精品| 身体一侧抽搐| 哪里可以看免费的av片| 久久午夜亚洲精品久久| 高潮久久久久久久久久久不卡| 精品国产亚洲在线| 久久久久久久久久黄片| 午夜福利视频1000在线观看| 欧美大码av| 亚洲电影在线观看av| 亚洲久久久国产精品| 久久久国产成人免费| 看免费av毛片| 好男人电影高清在线观看| 亚洲av日韩精品久久久久久密| 99久久国产精品久久久| 亚洲精品av麻豆狂野| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 久久这里只有精品19| 中文字幕高清在线视频| 亚洲 国产 在线| 精品久久久久久久毛片微露脸| 久久久精品国产亚洲av高清涩受| 色综合婷婷激情| 午夜福利18| 国产伦人伦偷精品视频| 人妻久久中文字幕网| 婷婷精品国产亚洲av| 久久精品国产清高在天天线| 精品欧美国产一区二区三| 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| av在线播放免费不卡| 精品久久久久久久末码| 国产片内射在线| 欧美日韩瑟瑟在线播放| 午夜福利欧美成人| а√天堂www在线а√下载| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www| 波多野结衣高清无吗| 女生性感内裤真人,穿戴方法视频| 久久午夜亚洲精品久久| 99精品在免费线老司机午夜| 国产午夜福利久久久久久| 午夜福利一区二区在线看| cao死你这个sao货| 久久午夜综合久久蜜桃| 日本成人三级电影网站| 婷婷六月久久综合丁香| 国产精品久久视频播放| 一本精品99久久精品77| 一进一出抽搐gif免费好疼| 亚洲av电影在线进入| 在线看三级毛片| 老司机深夜福利视频在线观看| 精华霜和精华液先用哪个| 九色国产91popny在线| 欧美成人午夜精品| 亚洲精品国产精品久久久不卡| 久久九九热精品免费| 国内精品久久久久精免费| 夜夜看夜夜爽夜夜摸| АⅤ资源中文在线天堂| 久久这里只有精品19| 欧美性长视频在线观看| 久久婷婷人人爽人人干人人爱| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品av在线| 亚洲国产毛片av蜜桃av| 亚洲av熟女| 午夜福利成人在线免费观看| 精品福利观看| 黄片小视频在线播放| 国产国语露脸激情在线看| 亚洲一区中文字幕在线| 黄色女人牲交| 别揉我奶头~嗯~啊~动态视频| 久久国产亚洲av麻豆专区| 久久人人精品亚洲av| cao死你这个sao货| 级片在线观看| 午夜福利高清视频| 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲国产一区二区在线观看| 一边摸一边做爽爽视频免费| 国产精品二区激情视频| 亚洲人成伊人成综合网2020| 岛国在线观看网站| 免费在线观看完整版高清| 国产麻豆成人av免费视频| 无限看片的www在线观看| 久久99热这里只有精品18| 最近在线观看免费完整版| 国产亚洲精品av在线| av欧美777| 真人一进一出gif抽搐免费| 久久精品国产综合久久久| avwww免费| tocl精华| 国产伦在线观看视频一区| 久久亚洲精品不卡| 久久久国产成人免费| 无限看片的www在线观看| 国产精品久久久久久精品电影 | 一级黄色大片毛片| av有码第一页| 国产v大片淫在线免费观看| 又紧又爽又黄一区二区| 国产三级黄色录像| 欧美一级毛片孕妇| 国产午夜福利久久久久久| 国产精品自产拍在线观看55亚洲| 国产激情久久老熟女| 亚洲一区二区三区色噜噜| 丝袜美腿诱惑在线| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 免费观看人在逋| 天堂影院成人在线观看| 嫩草影视91久久| 国产精品久久久久久人妻精品电影| 黑人欧美特级aaaaaa片| av中文乱码字幕在线| av超薄肉色丝袜交足视频| 男男h啪啪无遮挡| 久久婷婷人人爽人人干人人爱| 三级毛片av免费| 成年版毛片免费区| 国产欧美日韩精品亚洲av| 亚洲欧美精品综合久久99| 国产精品1区2区在线观看.| 在线天堂中文资源库| 亚洲欧美精品综合久久99| 在线观看www视频免费| 久久中文字幕一级| 91大片在线观看| 最好的美女福利视频网| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 男人舔女人下体高潮全视频| 在线观看免费午夜福利视频| 一本综合久久免费| 91麻豆精品激情在线观看国产| 一区二区三区精品91| 国产成+人综合+亚洲专区| 黄片小视频在线播放| 国产v大片淫在线免费观看| 一区二区三区高清视频在线| 亚洲av熟女| 99精品在免费线老司机午夜| 男女下面进入的视频免费午夜 | 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久精品电影 | 人人妻人人澡欧美一区二区| 午夜福利视频1000在线观看| 国产成+人综合+亚洲专区| 长腿黑丝高跟| 成人三级做爰电影| 自线自在国产av| 国产成人影院久久av| 欧美绝顶高潮抽搐喷水| 日韩视频一区二区在线观看| 国产1区2区3区精品| 十八禁网站免费在线| 久久国产精品人妻蜜桃| 国产又色又爽无遮挡免费看| 中文字幕人妻熟女乱码| 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区| 女人被狂操c到高潮| 一边摸一边做爽爽视频免费| 免费看美女性在线毛片视频| 两个人看的免费小视频| 亚洲国产精品sss在线观看| 1024手机看黄色片| 国产乱人伦免费视频| 亚洲狠狠婷婷综合久久图片| 亚洲欧美精品综合久久99| 老汉色∧v一级毛片| 国产精品久久久久久人妻精品电影| 久久 成人 亚洲| 精品国产美女av久久久久小说| 十分钟在线观看高清视频www| 一级黄色大片毛片| 国产亚洲精品久久久久5区| 脱女人内裤的视频| 少妇被粗大的猛进出69影院| 亚洲欧美激情综合另类| 亚洲精品在线美女| 亚洲黑人精品在线| 91九色精品人成在线观看| 亚洲成人国产一区在线观看| 少妇的丰满在线观看| 国内少妇人妻偷人精品xxx网站 | 美女大奶头视频| 亚洲人成电影免费在线| 麻豆成人av在线观看| 99国产综合亚洲精品| 亚洲欧美日韩高清在线视频| 亚洲片人在线观看| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看 | 久久青草综合色| 两个人免费观看高清视频| 亚洲免费av在线视频| 老汉色∧v一级毛片| 成人亚洲精品av一区二区| 亚洲一区二区三区色噜噜| 夜夜爽天天搞| 亚洲一区二区三区不卡视频| 国产av不卡久久| 伦理电影免费视频| 美女免费视频网站| 无人区码免费观看不卡| 精品不卡国产一区二区三区| a级毛片a级免费在线| 久久国产精品人妻蜜桃| avwww免费| 夜夜夜夜夜久久久久| 国产一区二区三区视频了| 亚洲精品中文字幕一二三四区| 国内精品久久久久精免费| 国产主播在线观看一区二区| 欧美丝袜亚洲另类 | netflix在线观看网站| 宅男免费午夜| 国产激情久久老熟女| 午夜福利18| 18禁黄网站禁片午夜丰满| 妹子高潮喷水视频| 一级片免费观看大全| 欧美国产精品va在线观看不卡| 亚洲性夜色夜夜综合| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| 国产片内射在线| 婷婷六月久久综合丁香| 婷婷丁香在线五月| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区| 自线自在国产av| 给我免费播放毛片高清在线观看| 午夜两性在线视频| 亚洲五月天丁香| 亚洲最大成人中文| 国产精品av久久久久免费| 国产1区2区3区精品| 国产精品九九99| 一级毛片高清免费大全| 美女高潮到喷水免费观看| 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 色综合欧美亚洲国产小说| 丁香欧美五月| 首页视频小说图片口味搜索| 老熟妇仑乱视频hdxx| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| 99在线人妻在线中文字幕| 国产亚洲欧美精品永久| 午夜福利18| 国产日本99.免费观看| 日韩免费av在线播放| 可以在线观看的亚洲视频| 老司机靠b影院| 国产av又大| 国产成人欧美| 可以在线观看的亚洲视频| 亚洲成人久久性| 大香蕉久久成人网| 变态另类成人亚洲欧美熟女| av超薄肉色丝袜交足视频| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| 美国免费a级毛片| 精品久久久久久久毛片微露脸| 亚洲av中文字字幕乱码综合 | 成人一区二区视频在线观看| 美女高潮喷水抽搐中文字幕| 久久伊人香网站| 亚洲自拍偷在线| 亚洲一区高清亚洲精品| 免费看a级黄色片| 日韩一卡2卡3卡4卡2021年| 自线自在国产av| 99精品久久久久人妻精品| 中文亚洲av片在线观看爽| 俺也久久电影网| 午夜两性在线视频| 亚洲五月婷婷丁香| 变态另类成人亚洲欧美熟女| 99久久无色码亚洲精品果冻| 成人国产综合亚洲| 欧美性猛交╳xxx乱大交人| 中文字幕av电影在线播放| 日本精品一区二区三区蜜桃| 黑人欧美特级aaaaaa片| 美女高潮喷水抽搐中文字幕| 成人18禁高潮啪啪吃奶动态图| 亚洲激情在线av| 国产片内射在线| 香蕉av资源在线| 最近最新中文字幕大全电影3 | 亚洲一卡2卡3卡4卡5卡精品中文| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 欧美亚洲日本最大视频资源| tocl精华| 欧美一区二区精品小视频在线| 久99久视频精品免费| 日韩免费av在线播放| 一级作爱视频免费观看| 亚洲成人久久性| 日本一区二区免费在线视频| 97碰自拍视频| 国产日本99.免费观看| 2021天堂中文幕一二区在线观 | 国产在线观看jvid| 国产激情久久老熟女| 国产免费av片在线观看野外av| 级片在线观看| 亚洲精品色激情综合| 怎么达到女性高潮| 国产一区二区在线av高清观看| 国产激情欧美一区二区| 国产一区二区三区在线臀色熟女| 啦啦啦观看免费观看视频高清| 1024香蕉在线观看| 亚洲在线自拍视频| 哪里可以看免费的av片| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 久久 成人 亚洲| 精品久久久久久久毛片微露脸| 91在线观看av| 俄罗斯特黄特色一大片| 日韩国内少妇激情av| 人成视频在线观看免费观看| 亚洲专区国产一区二区| 午夜福利视频1000在线观看| 丝袜在线中文字幕| 久久这里只有精品19| 久久人人精品亚洲av| 成人三级做爰电影| 丰满的人妻完整版| 黄片播放在线免费| 女生性感内裤真人,穿戴方法视频| 看黄色毛片网站| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 最近最新免费中文字幕在线| 国产亚洲精品久久久久久毛片| 国产av又大| 亚洲av成人不卡在线观看播放网| 桃色一区二区三区在线观看| 搡老岳熟女国产| 亚洲第一电影网av| 国产精品二区激情视频| 精品久久久久久久人妻蜜臀av| av在线天堂中文字幕| 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 麻豆久久精品国产亚洲av| 国产成人一区二区三区免费视频网站| 国产亚洲精品第一综合不卡| 午夜福利成人在线免费观看| 久久天堂一区二区三区四区| 人妻丰满熟妇av一区二区三区| 中文字幕人成人乱码亚洲影| 欧美国产日韩亚洲一区| 黄片小视频在线播放| 黄色毛片三级朝国网站| 成在线人永久免费视频| 色老头精品视频在线观看| a在线观看视频网站| 很黄的视频免费| 午夜福利18| 久久伊人香网站| 一二三四社区在线视频社区8| 一区二区三区国产精品乱码| 欧美日韩乱码在线| 超碰成人久久| 日韩大码丰满熟妇| 日韩免费av在线播放| 成人国语在线视频| 日韩欧美一区二区三区在线观看| 人妻丰满熟妇av一区二区三区| 免费高清在线观看日韩| 亚洲最大成人中文| 久久久久九九精品影院| 亚洲全国av大片| 欧美一级毛片孕妇| 精品乱码久久久久久99久播| 一级毛片精品| 亚洲精品中文字幕一二三四区| 一边摸一边做爽爽视频免费| 可以在线观看的亚洲视频| 国产精品久久久av美女十八| 一本精品99久久精品77| 视频区欧美日本亚洲| 岛国在线观看网站|