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

    真核生物中支架蛋白家族(IQGAPfamily)的分子進(jìn)化分析

    2017-10-13 03:36:19熱孜亞吾買爾曹立環(huán)
    關(guān)鍵詞:真核同源結(jié)構(gòu)域

    熱孜亞·吾買爾,曹立環(huán),余 龍

    (1. 復(fù)旦大學(xué) 生命科學(xué)學(xué)院,上海 200438; 2. 復(fù)旦大學(xué) 遺傳工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200438)

    真核生物中支架蛋白家族(IQGAPfamily)的分子進(jìn)化分析

    熱孜亞·吾買爾1,2,曹立環(huán)1,2,余 龍1,2

    (1. 復(fù)旦大學(xué) 生命科學(xué)學(xué)院,上海 200438; 2. 復(fù)旦大學(xué) 遺傳工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200438)

    支架蛋白家族(IQGAP family)廣泛存在于真核生物中,在細(xì)胞的信號(hào)轉(zhuǎn)導(dǎo)、骨架運(yùn)動(dòng)、細(xì)胞分裂的過(guò)程中都有著重要功能.研究分析IQGAP蛋白家族的分子進(jìn)化,有助于深入全面地了解整個(gè)IQGAP蛋白家族在不同物種中功能.通過(guò)使用相關(guān)數(shù)據(jù)庫(kù),在35個(gè)真核生物物種中檢索到了70個(gè)支架蛋白家族成員的序列信息,并對(duì)其進(jìn)行了系統(tǒng)的分子進(jìn)化研究.我們發(fā)現(xiàn): 支架蛋白家族基因廣泛存在于后生動(dòng)物(Metazoan)、真菌(Fungi)、變形蟲界(Amoebozoa)以及其他一些真核生物中,提示支架蛋白存在于最后的真核生物的共同祖先(Last Eukaryotic Common Ancestor, LECA).系統(tǒng)發(fā)育分析表明脊椎動(dòng)物的IQGAP蛋白和非脊椎動(dòng)物中代表性IQGAP基因共享一個(gè)共同的祖先基因,單個(gè)的IQGAP祖先基因在早期脊椎動(dòng)物中(在四足總綱和硬骨魚綱的分離之前),發(fā)生了兩次基因復(fù)制事件,導(dǎo)致現(xiàn)存的脊椎動(dòng)物中3組IQGAP基因(IQGAP1,IQGAP2,IQGAP3)產(chǎn)生.

    真核生物; 支架蛋白家族; 分子進(jìn)化分析

    基因的分子進(jìn)化研究能通過(guò)序列同源性的比較,系統(tǒng)發(fā)育樹的構(gòu)建等方法研究基因在不同物種中的起源,復(fù)制,突變,丟失的信息,這有助于人們了解在不同生物物種中基因結(jié)構(gòu)變化及其功能衍化的規(guī)律,使人們更深入全面地了解基因的功能.

    支架蛋白家族(IQGAP family)是一種進(jìn)化上保守的多結(jié)構(gòu)域蛋白質(zhì),IQGAP蛋白在酵母,魚類,非洲爪蟾和哺乳動(dòng)物中,都有廣泛表達(dá)[1].在真菌中,釀酒酵母(Saccharomycescerevisiae)IQGAP同源蛋白iqg1和裂殖酵母(Schizosaccharomycespombe)IQGAP同源蛋白R(shí)NG2都在胞質(zhì)分裂過(guò)程中發(fā)揮重要作用,兩者都定位于細(xì)胞收縮環(huán)和負(fù)責(zé)招募其他參與細(xì)胞分裂的蛋白[2-3].變形蟲界的盤基網(wǎng)柄菌(Dictyosteliumdiscoideum)的兩種IQGAP同源蛋白R(shí)GAA和GAPA也都被發(fā)現(xiàn)在胞質(zhì)分裂時(shí)發(fā)揮關(guān)鍵作用[4-5].人擁有3個(gè)IQGAP家族的成員(IQGAP1,IQGAP2和IQGAP3),研究發(fā)現(xiàn),IQGAP家族的成員在各組織中有不同表達(dá)模式.人IQGAP1蛋白在人體的組織中廣泛表達(dá),而人IQGAP2僅表達(dá)于肝臟和睪丸,人IQGAP3主要在腦和肺中表達(dá)[6].當(dāng)前,人們對(duì)IQGAP蛋白的研究主要集中在IQGAP1蛋白.IQGAPl分子大小190kD,位于染色體15q26.IQGAPl有鈣調(diào)理蛋白同源結(jié)構(gòu)域(Calponin Homology Domain, CHD),IQ結(jié)構(gòu)域(IQ domain),聚脯氨酸結(jié)合區(qū)(polyproline binding region, Ww),RasGAP結(jié)構(gòu)域(RasGap domain),以及在蛋白C末端的PfamRasGAP_C結(jié)構(gòu)域[7].其中,IQGAPl的鈣調(diào)理蛋白同源結(jié)構(gòu)域能結(jié)合肌動(dòng)蛋白(Actin),聚脯氨酸結(jié)合區(qū)能結(jié)合ERK蛋白,IQGAP1上的RasGap結(jié)構(gòu)域能結(jié)合Cdc42蛋白和Rac1蛋白,PfamRasGAP_C 結(jié)構(gòu)域能結(jié)合鈣黏附蛋白E(E-cadherin)和β-鏈蛋白(β-catenin),IQGAP1被認(rèn)為是一個(gè)介導(dǎo)多個(gè)信號(hào)通路的樞紐,參與細(xì)胞生長(zhǎng)、分化和分裂[8].人IQGAP1和IQGAP3在HeLa細(xì)胞的胞質(zhì)分裂中發(fā)揮重要作用,RNA干擾IQGAP1或IQGAP3地表達(dá)能導(dǎo)致胞質(zhì)分裂受阻,多核細(xì)胞產(chǎn)生[9].

    最近,關(guān)于人IQGAP家族在腫瘤發(fā)生發(fā)展中的作用也引起了人們的注意.在肝癌細(xì)胞中,人IQGAP1蛋白表達(dá)水平顯著上升,而人IQGAP2蛋白表達(dá)水平卻下降[10].同時(shí),肝癌病人中的人IQGAP1,IQGAP2的表達(dá)水平與病人的愈后相關(guān)[11].

    盡管IQGAP蛋白家族分子被發(fā)現(xiàn)和研究已經(jīng)20多年了,但人們對(duì)整個(gè)IQGAP家族在多種生物中的功能研究依然十分有限,對(duì)IQGAP蛋白家族的分子進(jìn)化歷史的了解也不明確.當(dāng)前,隨著越來(lái)越多生物的基因組被完全測(cè)序,使我們有機(jī)會(huì)在相對(duì)廣泛的真核生物物種中探討IQGAP基因的分子進(jìn)化歷史,了解不同物種間的IQGAP分子以及同一物種間不同的IQGAP分子之間的進(jìn)化關(guān)系,幫助我們進(jìn)一步深入地理解IQGAP蛋白的功能.

    1 材料與方法

    1.1真核生物中IQGAP蛋白的檢索

    我們下載了SMART蛋白結(jié)構(gòu)域數(shù)據(jù)庫(kù)[12](http:∥smart.embl-heidelberg.de/)中RasGAP結(jié)構(gòu)域(SM00323)的種子序列(alignment Consensus sequences)的序列,利用HMMER程序在非冗余蛋白序列數(shù)據(jù)庫(kù)進(jìn)行HMM搜索(http:∥hmmer.janelia.org/search/hmmsearch)[13],所得到的序列在計(jì)算機(jī)上用BLAST程序?qū)θ说鞍仔蛄袛?shù)據(jù)庫(kù)進(jìn)行BLAST搜索[14],僅保留其中BEST Hit為人IQGAP蛋白的序列,而其他的BEST Hit為人RASA1, RASL3等蛋白則被丟棄.對(duì)于每一個(gè)潛在的IQGAP家族基因,我們最后只收錄蛋白序列最長(zhǎng)的剪接本,其他剪接本或重復(fù)提交序列都沒(méi)有被最后收錄.最后,我們選擇了35個(gè)物種中的IQGAP蛋白進(jìn)行研究,35個(gè)物種中包含16個(gè)后生生物物種:Homosapiens(人);Musmusculus(小鼠);Rattusnorvegicus(大鼠);Bostaurus(牛);Gallusgallus(原雞);Xenopuslaevis(非洲爪蟾);Xenopustropicalis(熱帶爪蟾);Daniorerio(斑馬魚);Tetraodonnigroviridis(綠河豚);Cionaintestinalis(玻璃海鞘);Cionasavignyi;Branchiostomafloridae(文昌魚);Caenorhabditiselegans(線蟲);Nematostellavectensis(???;Trichoplaxadhaerens(絲盤蟲);Amphimedonqueenslandica(海綿).3個(gè)和后生動(dòng)物(Metazoa)親緣關(guān)系很近的原生生物物種:Monosigabrevicollis;Salpingoecarosetta;Capsasporaowczarzaki.8個(gè)真菌(Fungi)物種:Allomycesmacrogynus;Candidaalbicans(白念珠菌);Coccidioidesimmitis(粗球孢子菌);Coprinopsiscinerea(灰蓋鬼傘);Neurosporacrassa(粗糙鏈孢霉);Saccharomycescerevisiae(釀酒酵母);Schizosaccharomycespombe(裂殖酵母);Spizellomycespunctatus.5個(gè)歸屬于變形蟲界(Amoebozoa)的物種:Acanthamoebacastellanii(卡氏棘阿米巴);Dictyosteliumdiscoideum(盤基網(wǎng)柄菌);Dictyosteliumpurpureum;Entamoebahistolytica(溶組織內(nèi)阿米巴);Polysphondyliumpallidum; 以及其他的3個(gè)物種:Naegleriagruberi(尾刺耐格里原蟲);Thecamonastrahens;Spongosporasubterranean.

    1.2蛋白結(jié)構(gòu)域組成分析

    利用SMART蛋白結(jié)構(gòu)域數(shù)據(jù)庫(kù)[12](http:∥smart.embl-heidelberg.de/)和Pfam蛋白結(jié)構(gòu)域數(shù)據(jù)庫(kù)[15](http:∥pfam.xfam.org/)對(duì)我們所收集的70個(gè)蛋白進(jìn)行結(jié)構(gòu)域組成分析.

    1.3蛋白多序列比對(duì)和系統(tǒng)發(fā)育分析

    我們利用Mafft程序[16],參數(shù)設(shè)定為(L-INS-I),進(jìn)行蛋白序列的多序列比對(duì).我們選擇多序列比對(duì)中的最保守區(qū)段(IQGAP蛋白中的RasGap結(jié)構(gòu)域部分)進(jìn)行進(jìn)一步的系統(tǒng)發(fā)育分析,我們先用ProtTest2.4程序[17],利用Akaike Information Criterion(AIC)標(biāo)準(zhǔn)計(jì)算了相應(yīng)的最優(yōu)蛋白替代模型,然后用最大似然法(Maximum Likelihood, ML)和貝葉斯推斷法(Bayesian Inference, BI)構(gòu)建進(jìn)化樹.我們基于CIPRES Science Gateway V. 3.1[18]利用RAxML 7.2.8程序[19]構(gòu)建ML樹,bootstrap值設(shè)定為1000;用PHYLOBAYES3.3程序[20]構(gòu)建貝葉斯進(jìn)化樹,停止蒙特卡羅馬爾可夫鏈計(jì)算的參數(shù)設(shè)定為Maxdiff<0.1,并使用Dendroscope程序[21]顯示所有進(jìn)化樹.

    2 結(jié) 果

    2.1IQGAP蛋白的檢索

    我們最后檢索到35個(gè)真核物種中的70個(gè)IQGAP蛋白的序列信息,所有IQGAP蛋白合乎下列3個(gè)條件: (1)含有RasGap結(jié)構(gòu)域(SM00323)或pfamRasGap結(jié)構(gòu)域(PF00616);(2)這些蛋白對(duì)人蛋白序列數(shù)據(jù)庫(kù)進(jìn)行BLAST搜索,Best hit為人的IQGAP1,IQGAP2或IQGAP3蛋白;(3)所有蛋白含有PfamRasGap_C(PF03836)結(jié)構(gòu).我們發(fā)現(xiàn)這些蛋白分布在后生動(dòng)物、真菌、變形蟲界和其他的原生生物中,具體結(jié)果見表1.本文中蛋白名稱均依據(jù)Uniprot蛋白數(shù)據(jù)庫(kù).

    表1 代表性真核生物物種中的IQGAP蛋白

    注: 此表的數(shù)據(jù),在匯集本文鑒定的IQGAP家族蛋白的同時(shí),部分整合了本文圖2中系統(tǒng)發(fā)育分析的數(shù)據(jù),其中脊椎動(dòng)物IQGAP被成功地歸類于3個(gè)小組(IQGAP1,IQGAP2,IQGAP3),而其他物種的IQGAP蛋白未能進(jìn)行這三個(gè)小組的分類.用黑體標(biāo)注的蛋白被認(rèn)為是脊椎動(dòng)物IQGAP的直系同源基因(參見正文).

    2.2各物種中IQGAP蛋白的蛋白結(jié)構(gòu)域組成

    首先,我們研究了這些蛋白的蛋白結(jié)構(gòu)域的組成,我們發(fā)現(xiàn),大部分IQGAP蛋白(圖1,蛋白4~13)都含有鈣調(diào)理蛋白同源結(jié)構(gòu)域,IQ結(jié)構(gòu)域,RasGAP結(jié)構(gòu)域,以及在C末端的PfamRasGAP_C結(jié)構(gòu)域.但在相當(dāng)多的真菌和變形蟲來(lái)源的IQGAP蛋白,例如來(lái)源于變形蟲界的L8GLC0_ACACA、L8HI44_ACACA、Q54I64_DICDI、RGAA_DICDI、GAPA_DICDI等蛋白,以及來(lái)源于真菌的0A0L0T9F0_ALLMA、J3KKA8_COCIM、A8NUP6_COPC7、Q7S5G6_NEUCR、AP1_SCHPO等蛋白,缺失鈣調(diào)理蛋白同源結(jié)構(gòu)域和IQ結(jié)構(gòu)域,但這些蛋白中的RasGAP結(jié)構(gòu)域和C末端的PfamRasGAP_C結(jié)構(gòu)域依然存在,限于篇幅,圖1中僅顯示了3個(gè)此類代表性蛋白(蛋白1~3),這顯示RasGAPs和PfamRasGAP_C這兩個(gè)結(jié)構(gòu)域是IQGAP家族最為核心的區(qū)域.另外,我們只在部分IQGAP蛋白(圖1,蛋白7,8,11)中檢測(cè)到聚脯氨酸結(jié)合區(qū)(Ww)蛋白結(jié)構(gòu)域.

    圖1 真核生物中代表性IQGAP蛋白的結(jié)構(gòu)域組成Fig.1 Domain organization of representative IQGAP proteins in Eukaryotes真核生物中的代表性IQGAP蛋白,依次標(biāo)示了蛋白名稱、序列長(zhǎng)度、物種名稱和結(jié)構(gòu)域組成.蛋白名稱依據(jù)Uniprot蛋白數(shù)據(jù)庫(kù),結(jié)構(gòu)域組成依據(jù)SMART蛋白結(jié)構(gòu)域數(shù)據(jù)庫(kù)和Pfam蛋白結(jié)構(gòu)域數(shù)據(jù)庫(kù)的結(jié)果.

    2.3真核生物中IQGAP蛋白的系統(tǒng)發(fā)育分析結(jié)果

    因?yàn)閬?lái)源于線蟲的O17772_CAEEL蛋白,釀酒酵母的IQG1_YEAST蛋白,裂殖酵母的RNG2_SCHPO蛋白,白念珠菌的IQG1_CANAL蛋白盡管是IQGAP的同源蛋白,但這4個(gè)蛋白與其他的IQGAP蛋白序列差異較大,會(huì)干擾我們整個(gè)系統(tǒng)發(fā)育分析的可信度.為了獲得有效的系統(tǒng)發(fā)育分析結(jié)果,我們的分析中沒(méi)有涵蓋上述的4個(gè)蛋白.這樣,我們對(duì)70個(gè)蛋白中余下的66個(gè)蛋白進(jìn)行了系統(tǒng)發(fā)育分析.由于來(lái)源于動(dòng)物、真菌、變形蟲界等的IQGAP蛋白相互之間的進(jìn)化距離較遠(yuǎn),故這些IQGAP蛋白序列比對(duì)中某些區(qū)域的序列相似性也較差.為獲得可信的系統(tǒng)發(fā)育分析結(jié)果,我們選擇IQGAP蛋白中的最保守蛋白序列區(qū)段(IQGAP蛋白中的RasGap domain部分)進(jìn)行多序列比對(duì)和進(jìn)化分析.

    我們對(duì)上述66個(gè)IQGAP家族的基因的蛋白序列用最大似然法(ML)和貝葉斯推斷法(BI)進(jìn)行了系統(tǒng)發(fā)育分析,來(lái)自人的RASA1_HUMAN, RASL3_HUMAN和SYGP1_HUMAN蛋白設(shè)定為我們進(jìn)化分析的外類群(outgroup),結(jié)果顯示最大似然法(ML)和貝葉斯推斷法(BI)產(chǎn)生了幾乎相同的拓?fù)浣Y(jié)構(gòu)(圖2).在圖2中所有的IQGAP蛋白家族也總結(jié)于表1.

    圖2顯示,所有的來(lái)自脊椎動(dòng)物(人,小鼠,大鼠,牛,原雞,非洲爪蟾,熱帶爪蟾,斑馬魚,綠河豚)的IQGAP蛋白都可以很好地歸類到3個(gè)IQGAP組,分別為IQGAP1Group,IQGAP2Group,IQGAP3Group.同時(shí),我們還發(fā)現(xiàn),來(lái)自??腁7SXE3_NEMVE,絲盤蟲的B3RR54_TRIAD,海綿的I1F9L9_AMPQE蛋白和脊椎動(dòng)物中的IQGAP蛋白很好的歸類到一起(圖中標(biāo)示為IQGAPsinAnimalia),這一結(jié)果提示這些基因應(yīng)共享一個(gè)共同的祖先基因.事實(shí)上,從圖中可以看到,一些來(lái)源于后生動(dòng)物進(jìn)化關(guān)系最近的單細(xì)胞生物的IQGAP蛋白(如A9USM2_MONBE,F(xiàn)2TWA6_SALR5),一些來(lái)源于真菌的IQGAP蛋白(圖中標(biāo)示為IQGAPsinFungi)和一些來(lái)源于變形蟲界的IQGAP蛋白(圖中標(biāo)示為IQGAPsinAmoebozoa),都可被認(rèn)為是脊椎動(dòng)物中的IQGAP蛋白直系同源基因.這提示最后的真核生物的共同祖先(Last Eukaryotic Common Ancestor, LECA)已擁有支架蛋白家族基因這一古老的基因.

    盡管脊椎動(dòng)物IQGAP蛋白與真菌的IQGAP蛋白,變形蟲界的IQGAP蛋白共享一個(gè)共同的祖先基因,但脊椎動(dòng)物中3組IQGAP基因(IQGAP1,IQGAP2,IQGAP3)的形成卻發(fā)生于相對(duì)晚近的時(shí)期.我們的分析結(jié)果顯示在早期脊椎動(dòng)物(在四足總綱和硬骨魚綱的分離之前)一個(gè)IQGAP的祖先基因,發(fā)生了兩次基因復(fù)制事件,其中第一次復(fù)制,形成IQGAP1/2的祖先基因和IQGAP3基因,然后,所產(chǎn)生的IQGAP1/2的祖先基因再次發(fā)生基因復(fù)制,產(chǎn)生相應(yīng)的IQGAP1和IQGAP2基因,這兩次基因復(fù)制事件導(dǎo)致脊椎動(dòng)物中3組IQGAP基因(IQGAP1,IQGAP2,IQGAP3)形成(圖3,見第380頁(yè)).

    在進(jìn)一步對(duì)多種脊椎動(dòng)物和非脊椎動(dòng)物中的IQGAP蛋白進(jìn)行多序列比對(duì)時(shí),我們發(fā)現(xiàn)人,小鼠,大鼠,牛的IQGAP2基因的蛋白序列相對(duì)于它們對(duì)應(yīng)的IQGAP1和IQGAP3發(fā)生了長(zhǎng)度為84個(gè)氨基酸的片段缺失(圖4,見第380頁(yè)),其缺失的片斷位于IQGAP蛋白的鈣調(diào)理蛋白同源結(jié)構(gòu)域和IQ結(jié)構(gòu)域之間.因?yàn)榧棺祫?dòng)物中熱帶爪蟾和斑馬魚中IQGAP2基因,以及非脊椎動(dòng)物中相應(yīng)直系同源基因(如海鞘的H2Z719_CIOSA,海綿的I1F9L9_AMPQE基因)未見此類型的蛋白片段缺失.提示部分脊椎動(dòng)物中發(fā)生此類片斷缺失發(fā)生在相對(duì)晚近的時(shí)期.但此類發(fā)生在部分脊椎動(dòng)物中的IQGAP2基因片斷缺失對(duì)IQGAP蛋白的功能影響還未見報(bào)道.

    圖2 真核生物中IQGAP蛋白的系統(tǒng)發(fā)育分析Fig.2 Phylogenetic analysis of IQGAP proteins in Eukaryotic lineages對(duì)真核生物中IQGAP家族的的蛋白序列用最大似然法(ML)和貝葉斯推斷法(BI)進(jìn)行了系統(tǒng)發(fā)育分析,結(jié)果顯示最大似然法(ML)和貝葉斯推斷法(BI)產(chǎn)生了幾乎相同的拓?fù)浣Y(jié)構(gòu).進(jìn)化樹分支中的第一個(gè)數(shù)字表明貝葉斯后驗(yàn)概率(只有這些關(guān)鍵分支被標(biāo)記),進(jìn)化樹分支中的第二個(gè)數(shù)字顯示最大似然法中最大似然比.來(lái)自人的RASA1_HUMAN, RASL3_HUMAN,和SYGP1_HUMAN蛋白設(shè)定為我們進(jìn)化分析的外類群(outgroup).蛋白名稱依據(jù)Uniprot蛋白數(shù)據(jù)庫(kù).

    圖3 IQGAP基因在早期脊椎動(dòng)物中復(fù)制的示意圖Fig.3 Summaries of the gene duplication events for IQGAP genes in early vertebrataIQGAP的祖先基因在早期脊椎動(dòng)物(在四足總綱和硬骨魚綱的分離之前)發(fā)生了兩次基因復(fù)制事件,導(dǎo)致脊椎動(dòng)物中3組IQGAP基因(IQGAP1,IQGAP2,IQGAP3)的形成.

    圖4 脊椎動(dòng)物和非脊椎動(dòng)物IQGAP蛋白的多序列比對(duì)Fig.4 Multiple sequence alignment of IQGAP proteins from vertebrate and invertebrate lineages來(lái)自人、小鼠、牛、斑馬魚、熱帶爪蟾、海鞘、海綿等的IQGAP蛋白進(jìn)行多序列比對(duì)(僅顯示與IQGAP2蛋白序列片斷缺失相近區(qū)域部份).蛋白名稱依據(jù)Uniprot蛋白數(shù)據(jù)庫(kù).

    3 討 論

    我們檢索了整個(gè)真核生物蛋白數(shù)據(jù)庫(kù)中的IQGAP家族蛋白,發(fā)現(xiàn)IQGAP蛋白在植物中缺失;然后,我們選擇了35個(gè)真核生物物種作為我們的研究范圍,對(duì)這些物種的選擇,首先考慮了物種廣泛性,35個(gè)真核生物物種包含了16種動(dòng)物物種,3種和后生動(dòng)物親緣關(guān)系很近的單細(xì)胞物種,8種真菌物種,5種變形蟲界和3種其他的物種.這些物種中包含了在動(dòng)物進(jìn)化歷史上有重要意義的代表性的生物.例如: 扁盤動(dòng)物門(Placozoa)中絲盤蟲(T.adhaerens)是一種原始的后生動(dòng)物,被認(rèn)為是目前發(fā)現(xiàn)的最簡(jiǎn)單的真后生動(dòng)物(Eumetazoan),蟲體由上下兩層細(xì)胞構(gòu)成,無(wú)體腔,無(wú)消化腔,無(wú)神經(jīng)系統(tǒng)[22].刺胞動(dòng)物門(Cnidaria)的???N.vectensis)是一個(gè)非兩側(cè)對(duì)稱動(dòng)物[23].領(lǐng)鞭毛蟲類(Choanoflagellate)的M.brevicollis被認(rèn)為是與后生動(dòng)物進(jìn)化關(guān)系最近的單細(xì)胞生物,在研究后生動(dòng)物的起源和多樣化有重要的價(jià)值[24].盤基網(wǎng)柄菌(D.discoideum)作為一種土壤變形蟲,具有特殊的能力,能以單細(xì)胞形式和多細(xì)胞形式交替的方式存在.在生物的進(jìn)化譜系上,盤基網(wǎng)柄菌和后生生物的分離,早于酵母和后生生物的分離,但晚于植物和后生生物的分離[25].

    在研究方法上,在同源蛋白搜索中,我們結(jié)合了HMMer方法[13]和BLAST方法[14],以保證方法的靈敏性和可靠性.多序列比對(duì)所用的Mafft程序[16],系統(tǒng)發(fā)育樹構(gòu)建所用的最大似然法的RAxML 7.2.8程序[17]和貝葉斯法的PHYLOBAYES3.3程序[18]都是當(dāng)前國(guó)際上進(jìn)行系統(tǒng)發(fā)育研究的主流方法和程序,這也保證了我們結(jié)果的可靠性.

    我們的研究結(jié)果表明,IQGAP基因存在于后生動(dòng)物、真菌、變形蟲界等真核生物中,在Excavata類的真核生物尾刺耐格里原蟲(Naegleriagruberi),在有孔蟲界(Rhizaria)的Spongosporasubterranean中也發(fā)現(xiàn)了IQGAP的同源基因(表1),這充分顯示了IQGAP基因廣泛分布于真核生物中,跟據(jù)上述物種在真核生物中的進(jìn)化關(guān)系[26].我們認(rèn)為支架蛋白存在于LECA,即在后生動(dòng)物,真菌,變形蟲界和有孔蟲界等真核生物的物種分離之前,IQGAP基因就已經(jīng)存在.但是,變形蟲界IQGAP基因一般僅擁有RasGAP結(jié)構(gòu)域,C末端的PfamRasGAP_C結(jié)構(gòu)域.而后生動(dòng)物,真菌中的IQGAP基因除了含有RasGAP結(jié)構(gòu)域,C末端的PfamRasGAP_C結(jié)構(gòu)域外,往往還含有鈣調(diào)理蛋白同源結(jié)構(gòu)域,IQ結(jié)構(gòu)域(圖1).關(guān)于各物種IQGAP基因的結(jié)構(gòu)域組成變化的進(jìn)化歷史,以及結(jié)構(gòu)域組成對(duì)IQGAP基因功能的影響將會(huì)是進(jìn)一步研究的對(duì)象.

    但人的3個(gè)IQGAP基因是在早期脊椎動(dòng)物中的一個(gè)IQGAP的祖先基因通過(guò)發(fā)生了兩次基因復(fù)制事件后形成的.其中第一次基因復(fù)制,形成IQGAP1/2的祖先基因和IQGAP3基因,然后,所產(chǎn)生的IQGAP1/2的祖先基因再次發(fā)生基因復(fù)制,產(chǎn)生相應(yīng)的IQGAP1和IQGAP2基因(圖3).這一結(jié)果與人IQGAP1,IQGAP2,IQGAP3基因相互之間較高的序列相似性相吻合.在人的3個(gè)IQGAP基因中,IQGAP1基因表達(dá)于人體絕大多數(shù)組織中,并受到廣泛的研究,而IQGAP2基因和IQGAP3研究則相對(duì)較少[6,27].人的3個(gè)IQGAP基因,除了在人體各組織的表達(dá)譜不同以外,功能也有所不同[27].人IQGAP1基因被認(rèn)為是一個(gè)癌基因(oncogene),在肝細(xì)胞癌、肺癌、胃癌、乳腺癌表達(dá)上升[28],而人IQGAP2基因在肝細(xì)胞癌,前列腺中表達(dá)下降[29-30].實(shí)驗(yàn)研究還表明敲除IQGAP2基因的小鼠易發(fā)生肝癌,但敲除IQGAP1和IQGAP2雙基因的小鼠則能正常生存,提示IQGAP1基因和IQGAP2基因發(fā)生某種程度的功能拮抗[31].而IQGAP3功能與人IQGAP1有相似性,兩者都能促進(jìn)肝細(xì)胞和乳腺上皮細(xì)胞的增殖[32-33].這樣,IQGAP2相對(duì)IQGAP1和IQGAP3基因,有其功能的特殊性.有趣的是我們發(fā)現(xiàn)人,小鼠,大鼠,牛的IQGAP2基因的蛋白序列相對(duì)于它們對(duì)應(yīng)的IQGAP1和IQGAP3發(fā)生了長(zhǎng)度為84個(gè)氨基酸的片段缺失(見圖4).現(xiàn)在還不清楚是否是IQGAP2上84個(gè)氨基酸的片段缺失造成了IQGAP2基因和IQGAP1基因的功能拮抗,進(jìn)一步研究IQGAP2基因上的84個(gè)氨基酸的片段缺失對(duì)IQGAP功能的影響將給我們以新的提示.

    我們的研究明確揭示了IQGAP家族在早期脊椎動(dòng)物中的進(jìn)化歷史,即在早期脊椎動(dòng)物中發(fā)生了兩次基因復(fù)制事件,導(dǎo)致現(xiàn)存的脊椎動(dòng)物中3組IQGAP基因(IQGAP1,IQGAP2,IQGAP3)產(chǎn)生.人,小鼠,大鼠,牛的IQGAP2基因的蛋白序列相對(duì)于它們對(duì)應(yīng)的IQGAP1和IQGAP3發(fā)生了長(zhǎng)度為84個(gè)氨基酸的片段缺失.這些結(jié)果將有助于我們進(jìn)一步研究人和多種其他生物中IQGAP蛋白的功能.

    [1] SHANNON K. IQGAP family members in yeast, dictyostelium, and mammalian cells [J].IntJCellBiol, 2012, 894817.

    [2] EPP A, CHANT J. An IQGAP-related protein controls actin-ring formation and cytokinesis in yeast [J].CurrBiol, 1997,7: 921-929.

    [3] ENG K, NAQVI N, WONG K,etal. Rng2p, a protein required for cytokinesis in fission yeast, is a component of the actomyosin ring and the spindle pole body [J].CurrBiol, 1998,8(11): 611-621.

    [4] LEE S, ESCALANTE R, FIRTEL R A. A Ras GAP is essential for cytokinesis and spatial patterning in Dictyostelium [J].Development, 1997,124: 983-996.

    [5] ADACHI H, TAKAHASHI Y, HASEBE T,etal. Dictyostelium IQGAP-related protein specifically involved in the completion of cytokinesis [J].JCellBiol, 1997,137: 891-898.

    [6] WANG S, WATANABE T, NORITAKE J. IQGAP3, a novel effector of Rac1 and Cdc42, regulates neurite outgrowth [J].JCellSci, 2007,120(4): 567-577.

    [7] HEDMAN A C, SMITH J M, SACKS D B. The biology of IQGAP proteins: Beyond the cytoskeleton [J].EMBORep, 2015,16, 427-446.

    [8] WHITE C D, ERDEMIR H H, SACKS D B. IQGAP1 and its binding proteins control diverse biological functions [J].CellSignal, 2012,24,826-834.

    [9] ADACHI M, KAWASAKI A, NOJIMA H,etal. Involvement of IQGAP family proteins in the regulation of mammalian cell cytokinesis [J].GenesCells, 2014,19(11): 803-820.

    [10] WHITE C D, KHURANA H, GNATENKO D V,etal. IQGAP1 and IQGAP2 are reciprocally altered in hepatocellular carcinoma [J].BMCGastroenterol, 2010,10: 125.

    [11] XIA F D, WANG Z L, CHEN H X,etal. Differential expression of IQGAP1/2 inHepatocellularcarcinomaand its relationship with clinical outcomes [J].AsianPacJCancerPrev, 2014,15(12): 4951-4956.

    [12] LETUNIC I, DOERKS T, BORK P. SMART: Recent updates, new developments and status in 2015 [J]NucleicAcidsRes, 2015,43(Database issue): 257-260.

    [13] FINN R D, CLEMENTS J, EDDY S R. HMMER Web Server: Interactive sequence similarity searching [J].NucleicAcidsRes, 2011,39: 29-37.

    [14] ALTSCHUL S F, GISH W, MILLER W,etal. Basic local alignment search tool [J].JMolBiol, 1990,215: 403-410.

    [15] FINN R D, COGGILL P, EBERHARDT R Y,etal. The Pfam protein families database: Towards a more sustainable future [J].NucleicAcidsRes, 2016,44(D1): 279-285.

    [16] KATOH K, KUMA K, TOH H, MIYATA T. MAFFT version 5: Improvement in accuracy of multiple sequence alignment[J].NucleicAcidsRes, 2005,33(2): 511-518.

    [17] DARRIBA D, TABOADA G L, DOALLO R,etal. ProtTest 3: Fast selection of best-fit models of protein evolution [J].Bioinformatics, 2011,27(8): 1164-1165.

    [18] STAMATAKIS A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies [J].Bioinformatics, 2014,30(9): 1312-1313.

    [19] MILLER M A, PFEIFFER W, SCHWARTZ T. Creating the CIPRES science gateway for inference of large phylogenetic trees [C]∥Proceedings of the Gateway Computing Environments Workshop(GCE). New Orleans, LA: ACM. 2010: 1-8.

    [20] LARTILLOT N, LEPAGE T, BLANQUART S. PhyloBayes 3: A Bayesian software package for phylogenetic reconstruction and molecular dating [J].Bioinformatics, 2009,25(17): 2286-2288.

    [21] HUSON D H, RICHTER D C, RAUSCH C,etal. Dendroscope: An interactive viewer for large phylogenetic trees [J].BMCBioinforma, 2007,8: 460.

    [22] SRIVASTAVA M, BEGOVIC E, CHAPMAN J,etal. The Trichoplax genome and the nature of placozoans [J].Nature, 2008,454(7207): 955-960.

    [23] PUTNAM N H, SRIVASTAVA M, HELLSTEN U,etal. Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization [J].Science, 2007,317(5834): 86-94.

    [24] KING N, WESTBROOK M J, YOUNG S L,etal. The genome of the choanoflagellateMonosigabrevicollisand the origin of metazoans [J].Nature, 2008,451(7180): 783-788.

    [25] EICHINGER L, PACHEBAT J A, GLOCKNER G,etal. The genome of the social amoebaDictyosteliumdiscoideum[J].Nature, 2005,435(7038): 43-57.

    [26] BURKI F. The eukaryotic tree of life from a global phylogenomic perspective [J].ColdSpringHarbPerspectBiol, 2014,6(5): a016147.

    [27] SMITH J M, HEDMAN A C, SACKS D B. IQGAPs choreograph cellular signaling from the membrane to the nucleus [J].TrendsCellBiol, 2015,25(3): 171-184.

    [28] JOHNSON M, SHARMA M, HENDERSON B R. IQGAP1 regulation and roles in cancer [J].CellSignal, 2009,21(10): 1471-1478.

    [29] GNATENKO D V, XU X, ZHU W,etal. Transcript profiling identifies IQGAP2(-/-) mouse as a model for advanced human hepatocellular carcinoma [J].PLoSOne, 2013,8(8): e71826.

    [30] XIE Y, YAN J, CUTZ J C,etal. IQGAP2, A candidate tumour suppressor of prostate tumorigenesis [J].BiochimBiophysActa, 2012,1822(6): 875-884.

    [31] SCHMIDT V A, CHIARIELLO C S, CAPILLA E,etal. Development of hepatocellular carcinoma in Iqgap2-deficient mice is IQGAP1 dependent [J].MolCellBiol, 2008,28(5): 1489-1502.

    [32] KUNIMOTO K, NOJIMA H, YAMAZAKI Y,etal. Involvement of IQGAP3, a regulator of Ras/ERK-related cascade, in hepatocyte proliferation in mouse liver regeneration and development [J].JCellPhysiol, 2009,220(3): 621-631.

    [33] NOJIMA H, ADACHI M, MATSUI T,etal. IQGAP3 regulates cell proliferation through the Ras/ERK signalling cascade [J].NatCellBiol, 2008,10(8): 971-978.

    Abstract: IQGAP family proteins present widely in eukaryotic lineages, and play important roles in the process of cell signal transduction, cell skeletal movement and cell division. Analysis of the molecular evolution of the IQGAP protein family is helpful in understanding the function of the whole IQGAP protein family in different species. Here we retrieved 70 members of the IQGAP family in 35 eukaryotic species, and carried a comprehensive phylogenetic analysis. We found that IQGAP family proteins exist in Metazoan, Fungi, Amoebozoa, as well as some other eukaryotes, suggesting that the common ancestor of IQGAP family present in the Last Eukaryotic Common the Ancestor (LECA). Phylogenetic analysis also showed that vertebrate IQGAP genes shared a common ancestral gene with several representative non vertebrate IQGAP genes, and a single ancestral IQGAP gene duplicated two times in early vertebrates (before tetrapoda and teleostei divergence), resulting in three groups of IQGAP genes (IQGAP1, IQGAP2, IQGAP3) formed in vertebrates.

    Keywords: Eukaryota;IQGAP family; Phylogenetic analysis

    PhylogeneticAnalysisofIQGAPProteinsinEukaryoticLineages

    REZIYA Wumaier1,2, CAO Lihuan1,2, YU Long1,2

    (1.StateKeyLaboratoryofGeneticEngineering,InstituteofGenetics,Shanghai200438,China; >2.SchoolofLifeSciences,FudanUniversity,Shanghai200438,China)

    Q51

    A

    0427-7104(2017)03-0374-09

    2016-09-09

    國(guó)家自然科學(xué)基金(31071186)

    熱孜亞·吾買爾(1990—),女,碩士研究生;余 龍,男,教授,通信聯(lián)系人,E-mail: longyu@fudan.edu.cn;曹立環(huán),男,副教授,通信聯(lián)系人,E-mail: lihuancao@fudan.edu.cn.

    猜你喜歡
    真核同源結(jié)構(gòu)域
    藥食同源
    ——紫 蘇
    兩岸年味連根同源
    以同源詞看《詩(shī)經(jīng)》的訓(xùn)釋三則
    真核翻譯起始因子-5A2在肝內(nèi)膽管癌中的表達(dá)及意義
    蛋白質(zhì)結(jié)構(gòu)域劃分方法及在線服務(wù)綜述
    虔誠(chéng)書畫乃同源
    重組綠豆BBI(6-33)結(jié)構(gòu)域的抗腫瘤作用分析
    組蛋白甲基化酶Set2片段調(diào)控SET結(jié)構(gòu)域催化活性的探討
    泛素結(jié)合結(jié)構(gòu)域與泛素化信號(hào)的識(shí)別
    人醛縮酶A干擾RNA真核表達(dá)載體的構(gòu)建
    国产精华一区二区三区| 真实男女啪啪啪动态图| av国产免费在线观看| 亚洲乱码一区二区免费版| 深夜a级毛片| 国产一区二区三区av在线| av免费观看日本| 国产真实伦视频高清在线观看| 欧美一区二区精品小视频在线| 久久精品91蜜桃| 国产免费福利视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲熟妇中文字幕五十中出| 亚洲伊人久久精品综合 | 性插视频无遮挡在线免费观看| 免费搜索国产男女视频| 99久久精品国产国产毛片| 久久精品国产亚洲网站| 国产淫语在线视频| 又爽又黄a免费视频| 国产精品蜜桃在线观看| 亚洲国产精品久久男人天堂| 蜜桃亚洲精品一区二区三区| 午夜亚洲福利在线播放| 舔av片在线| 久久99热这里只频精品6学生 | 午夜激情欧美在线| 女人被狂操c到高潮| 国产精品伦人一区二区| 亚洲欧美一区二区三区国产| 人妻少妇偷人精品九色| 日本欧美国产在线视频| 国产人妻一区二区三区在| 久久精品夜色国产| 国产老妇伦熟女老妇高清| 两个人视频免费观看高清| 淫秽高清视频在线观看| 青春草亚洲视频在线观看| 天美传媒精品一区二区| 亚洲精品影视一区二区三区av| 一本一本综合久久| 一级av片app| 91精品伊人久久大香线蕉| 九九热线精品视视频播放| 久久久精品欧美日韩精品| 舔av片在线| 一级毛片我不卡| 91午夜精品亚洲一区二区三区| 久久婷婷人人爽人人干人人爱| 国产精品一及| 少妇的逼好多水| 91在线精品国自产拍蜜月| 国模一区二区三区四区视频| 成人无遮挡网站| 最后的刺客免费高清国语| 久久精品夜色国产| 欧美成人精品欧美一级黄| 国产精品永久免费网站| 老师上课跳d突然被开到最大视频| 色尼玛亚洲综合影院| 亚洲在线自拍视频| 波多野结衣巨乳人妻| 不卡视频在线观看欧美| 国内精品美女久久久久久| 亚洲成人中文字幕在线播放| 99久久无色码亚洲精品果冻| 国产色爽女视频免费观看| 国产老妇女一区| 国产精品久久久久久精品电影| 天堂网av新在线| 色网站视频免费| eeuss影院久久| 国产亚洲91精品色在线| 综合色丁香网| 国模一区二区三区四区视频| 桃色一区二区三区在线观看| 欧美3d第一页| 精品久久久久久成人av| 亚洲在线自拍视频| 97超碰精品成人国产| 日本熟妇午夜| 在线a可以看的网站| 亚洲成人av在线免费| 少妇被粗大猛烈的视频| 亚洲高清免费不卡视频| 最近最新中文字幕大全电影3| 亚洲色图av天堂| 最新中文字幕久久久久| 男插女下体视频免费在线播放| 国产又黄又爽又无遮挡在线| 午夜福利高清视频| 男女边吃奶边做爰视频| 精品人妻熟女av久视频| 亚洲精品乱码久久久v下载方式| 免费电影在线观看免费观看| 91aial.com中文字幕在线观看| 亚洲欧美日韩无卡精品| 国产私拍福利视频在线观看| 国产淫片久久久久久久久| 免费人成在线观看视频色| 欧美又色又爽又黄视频| 美女黄网站色视频| 三级毛片av免费| 亚洲成av人片在线播放无| 亚洲精品乱码久久久久久按摩| 97在线视频观看| 午夜激情欧美在线| 精品酒店卫生间| 九九久久精品国产亚洲av麻豆| 欧美成人一区二区免费高清观看| 亚洲丝袜综合中文字幕| ponron亚洲| 少妇的逼好多水| 久久午夜福利片| 久久精品人妻少妇| 日韩精品有码人妻一区| 九九热线精品视视频播放| 狠狠狠狠99中文字幕| 国产老妇伦熟女老妇高清| 国产91av在线免费观看| 久久这里有精品视频免费| 看免费成人av毛片| 身体一侧抽搐| 91av网一区二区| 美女黄网站色视频| 一个人观看的视频www高清免费观看| 桃色一区二区三区在线观看| 亚洲av成人精品一区久久| 欧美激情国产日韩精品一区| 国产一级毛片七仙女欲春2| 欧美日韩在线观看h| 有码 亚洲区| 能在线免费看毛片的网站| 麻豆乱淫一区二区| 免费观看性生交大片5| 好男人在线观看高清免费视频| 精品久久久久久电影网 | 97热精品久久久久久| 在线观看一区二区三区| 日韩制服骚丝袜av| 国产国拍精品亚洲av在线观看| 久久99热这里只频精品6学生 | av卡一久久| 毛片一级片免费看久久久久| 精品欧美国产一区二区三| 亚洲国产成人一精品久久久| 啦啦啦韩国在线观看视频| 欧美日本视频| 91狼人影院| 黄色配什么色好看| av在线蜜桃| 久久久久久伊人网av| 在线观看美女被高潮喷水网站| 精品一区二区三区视频在线| 亚洲国产成人一精品久久久| 美女被艹到高潮喷水动态| 人人妻人人澡人人爽人人夜夜 | 日韩成人伦理影院| 免费观看在线日韩| 亚洲成人精品中文字幕电影| av免费在线看不卡| 最近视频中文字幕2019在线8| 一区二区三区乱码不卡18| 亚洲欧美精品综合久久99| 麻豆成人av视频| ponron亚洲| 亚洲18禁久久av| 午夜福利视频1000在线观看| 九九爱精品视频在线观看| 亚洲无线观看免费| 熟妇人妻久久中文字幕3abv| 色网站视频免费| 搡女人真爽免费视频火全软件| 男女视频在线观看网站免费| 亚洲美女搞黄在线观看| 夫妻性生交免费视频一级片| 最新中文字幕久久久久| 最近中文字幕高清免费大全6| 午夜激情福利司机影院| 国产成人精品一,二区| 黄色日韩在线| 只有这里有精品99| 97超视频在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 最近手机中文字幕大全| av国产免费在线观看| a级毛片免费高清观看在线播放| av在线播放精品| 久久草成人影院| 真实男女啪啪啪动态图| 黄色日韩在线| 国产av在哪里看| 成人国产麻豆网| 一个人免费在线观看电影| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 伦理电影大哥的女人| 成人鲁丝片一二三区免费| 国产久久久一区二区三区| 久久精品国产鲁丝片午夜精品| 亚洲在久久综合| 日韩中字成人| 美女cb高潮喷水在线观看| 精品久久久久久久人妻蜜臀av| 久久精品国产自在天天线| 大香蕉97超碰在线| 亚洲欧美日韩东京热| 日本免费a在线| 欧美3d第一页| 91久久精品电影网| 少妇裸体淫交视频免费看高清| av线在线观看网站| av在线观看视频网站免费| 欧美成人a在线观看| 国产大屁股一区二区在线视频| 神马国产精品三级电影在线观看| 男插女下体视频免费在线播放| 午夜免费男女啪啪视频观看| 国产一区二区在线观看日韩| 久久久国产成人免费| 亚洲在久久综合| 美女黄网站色视频| 久久久久久久国产电影| 男女啪啪激烈高潮av片| 春色校园在线视频观看| 久久久久久久久久成人| 精品久久国产蜜桃| 国产 一区 欧美 日韩| 久久久午夜欧美精品| 国产精品久久视频播放| 色5月婷婷丁香| av免费在线看不卡| 亚洲一区高清亚洲精品| 成人毛片a级毛片在线播放| 一夜夜www| 国产成人精品一,二区| 深夜a级毛片| 一本一本综合久久| 夫妻性生交免费视频一级片| 舔av片在线| 天堂影院成人在线观看| 精华霜和精华液先用哪个| 中文乱码字字幕精品一区二区三区 | 久久精品久久久久久久性| 亚洲中文字幕日韩| 欧美zozozo另类| 亚洲丝袜综合中文字幕| 美女黄网站色视频| 久久久国产成人免费| 老司机福利观看| 99久国产av精品国产电影| 欧美3d第一页| 麻豆一二三区av精品| 中文天堂在线官网| 国产亚洲午夜精品一区二区久久 | 中文字幕熟女人妻在线| 久久久久性生活片| 日日摸夜夜添夜夜添av毛片| 91午夜精品亚洲一区二区三区| 欧美一区二区国产精品久久精品| 成人无遮挡网站| 在线播放国产精品三级| av福利片在线观看| 国产成人91sexporn| 好男人视频免费观看在线| 内地一区二区视频在线| 成人亚洲欧美一区二区av| 综合色av麻豆| 国产成人精品一,二区| 欧美日韩精品成人综合77777| 老司机影院成人| 黄色配什么色好看| 免费观看人在逋| 99热精品在线国产| 久久久久久久久久黄片| 国产91av在线免费观看| 精品人妻视频免费看| 一本一本综合久久| 精品少妇黑人巨大在线播放 | 大香蕉97超碰在线| 国产91av在线免费观看| av福利片在线观看| 爱豆传媒免费全集在线观看| 久久久色成人| 免费观看精品视频网站| 又粗又爽又猛毛片免费看| 一级毛片久久久久久久久女| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 午夜精品在线福利| 国产av码专区亚洲av| 视频中文字幕在线观看| kizo精华| 精品久久国产蜜桃| 成人综合一区亚洲| 亚洲精品乱久久久久久| 尤物成人国产欧美一区二区三区| 亚洲欧美成人综合另类久久久 | 久久久成人免费电影| 国产精品一二三区在线看| 两个人的视频大全免费| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 偷拍熟女少妇极品色| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| 欧美一区二区精品小视频在线| av国产久精品久网站免费入址| 亚洲无线观看免费| 青春草亚洲视频在线观看| 99热精品在线国产| 国产成人精品婷婷| 真实男女啪啪啪动态图| 色吧在线观看| 亚洲精品乱码久久久v下载方式| 黄片无遮挡物在线观看| 国产亚洲一区二区精品| 汤姆久久久久久久影院中文字幕 | 久久精品国产亚洲网站| 中文天堂在线官网| 亚洲欧美中文字幕日韩二区| av在线蜜桃| 性色avwww在线观看| 国产精品一区二区三区四区免费观看| 国产黄色小视频在线观看| 国产淫语在线视频| 少妇的逼好多水| 国产 一区 欧美 日韩| 成人国产麻豆网| 国产一级毛片在线| 久久国产乱子免费精品| 美女高潮的动态| 日韩av在线免费看完整版不卡| 男插女下体视频免费在线播放| 三级经典国产精品| 91精品伊人久久大香线蕉| 内射极品少妇av片p| 日韩强制内射视频| 国产成人午夜福利电影在线观看| 成年女人永久免费观看视频| 亚洲av中文字字幕乱码综合| 色哟哟·www| 精品久久久噜噜| 国产精品,欧美在线| 亚洲精品乱码久久久v下载方式| 91av网一区二区| 色视频www国产| 亚洲精品日韩在线中文字幕| 少妇人妻精品综合一区二区| 亚洲av中文av极速乱| 人妻少妇偷人精品九色| 日韩人妻高清精品专区| 免费黄色在线免费观看| 男女啪啪激烈高潮av片| 99热6这里只有精品| 我的老师免费观看完整版| 久久久久九九精品影院| 插逼视频在线观看| 午夜福利视频1000在线观看| 亚洲中文字幕一区二区三区有码在线看| 国模一区二区三区四区视频| av天堂中文字幕网| 亚洲av电影在线观看一区二区三区 | 国产探花极品一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品一区蜜桃| 国产精品久久久久久久电影| 少妇的逼好多水| 亚洲成人久久爱视频| 尾随美女入室| 久久精品国产自在天天线| 男人和女人高潮做爰伦理| kizo精华| 99久久精品热视频| 免费黄色在线免费观看| 欧美日韩综合久久久久久| 天美传媒精品一区二区| 黄色一级大片看看| 天天躁日日操中文字幕| av天堂中文字幕网| 春色校园在线视频观看| 精品一区二区三区人妻视频| 亚洲国产色片| 国产女主播在线喷水免费视频网站 | 亚洲精品成人久久久久久| 国产人妻一区二区三区在| 日韩,欧美,国产一区二区三区 | 国产成人精品一,二区| 村上凉子中文字幕在线| 一区二区三区乱码不卡18| 伦精品一区二区三区| 亚洲av成人精品一二三区| 日本欧美国产在线视频| 精品少妇黑人巨大在线播放 | 免费看日本二区| 国产午夜福利久久久久久| 中国国产av一级| 丝袜喷水一区| 精品人妻视频免费看| 岛国在线免费视频观看| 啦啦啦观看免费观看视频高清| www日本黄色视频网| videossex国产| 中文字幕久久专区| 美女大奶头视频| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 亚洲美女视频黄频| 亚洲人与动物交配视频| 色综合色国产| 一级毛片电影观看 | 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久 | 午夜激情福利司机影院| 国产三级在线视频| 亚洲av电影不卡..在线观看| 久久久成人免费电影| 中文在线观看免费www的网站| 99久久中文字幕三级久久日本| 男女那种视频在线观看| 国产亚洲午夜精品一区二区久久 | 国产综合懂色| 99在线人妻在线中文字幕| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 伦精品一区二区三区| 免费观看性生交大片5| 国产精品.久久久| 亚洲婷婷狠狠爱综合网| av天堂中文字幕网| 一本久久精品| 成人毛片a级毛片在线播放| 三级毛片av免费| 欧美日本亚洲视频在线播放| 欧美成人免费av一区二区三区| 欧美bdsm另类| 老师上课跳d突然被开到最大视频| 爱豆传媒免费全集在线观看| 在线播放国产精品三级| 免费观看a级毛片全部| 狂野欧美白嫩少妇大欣赏| 国产精品.久久久| 色综合站精品国产| 国产在线一区二区三区精 | 色综合站精品国产| 日韩欧美国产在线观看| 亚洲成人精品中文字幕电影| 成人av在线播放网站| 波多野结衣高清无吗| 麻豆av噜噜一区二区三区| 国产一区二区三区av在线| 亚洲国产成人一精品久久久| 日韩大片免费观看网站 | 亚洲美女视频黄频| 亚洲在线自拍视频| 国产精品爽爽va在线观看网站| 亚洲性久久影院| 久久6这里有精品| 小说图片视频综合网站| 精品人妻视频免费看| 男插女下体视频免费在线播放| 三级国产精品欧美在线观看| 久久人人爽人人片av| 精品人妻视频免费看| 亚洲美女搞黄在线观看| 亚洲国产欧洲综合997久久,| 亚洲成av人片在线播放无| 两个人视频免费观看高清| 在线a可以看的网站| 亚洲欧美日韩高清专用| 只有这里有精品99| 99热这里只有精品一区| 亚洲一级一片aⅴ在线观看| 可以在线观看毛片的网站| 中国国产av一级| 国产毛片a区久久久久| 性色avwww在线观看| 一区二区三区高清视频在线| 成人鲁丝片一二三区免费| 热99在线观看视频| 大又大粗又爽又黄少妇毛片口| 国产亚洲最大av| 不卡视频在线观看欧美| 午夜免费激情av| 成人毛片a级毛片在线播放| 国产精品一二三区在线看| 免费播放大片免费观看视频在线观看 | 熟妇人妻久久中文字幕3abv| 丰满少妇做爰视频| 91午夜精品亚洲一区二区三区| 亚洲精品456在线播放app| 久久久久久久久大av| 大又大粗又爽又黄少妇毛片口| av在线天堂中文字幕| 亚洲经典国产精华液单| 亚洲综合精品二区| 小说图片视频综合网站| 亚洲精品影视一区二区三区av| 国产成人一区二区在线| 淫秽高清视频在线观看| 国产探花在线观看一区二区| 国产精品,欧美在线| 国产久久久一区二区三区| 国产一区二区在线观看日韩| 又爽又黄无遮挡网站| 久久99蜜桃精品久久| 一级黄片播放器| 男女啪啪激烈高潮av片| 午夜爱爱视频在线播放| 2022亚洲国产成人精品| 欧美成人午夜免费资源| 男女边吃奶边做爰视频| 精品酒店卫生间| 国产老妇伦熟女老妇高清| 中文欧美无线码| 午夜久久久久精精品| 亚洲精品国产成人久久av| 白带黄色成豆腐渣| 麻豆成人午夜福利视频| 五月伊人婷婷丁香| 啦啦啦观看免费观看视频高清| 我要搜黄色片| 亚洲国产高清在线一区二区三| 韩国av在线不卡| av国产免费在线观看| 一级爰片在线观看| 99九九线精品视频在线观看视频| 精品不卡国产一区二区三区| 久久人人爽人人片av| 久久婷婷人人爽人人干人人爱| 成人欧美大片| 69av精品久久久久久| 国产免费一级a男人的天堂| 久久精品国产99精品国产亚洲性色| 日韩大片免费观看网站 | 国产成年人精品一区二区| 欧美日本视频| 亚洲精品国产av成人精品| av在线观看视频网站免费| 午夜福利成人在线免费观看| 男人舔女人下体高潮全视频| 国产91av在线免费观看| 亚洲欧洲国产日韩| 特大巨黑吊av在线直播| 亚洲精品日韩av片在线观看| 又爽又黄无遮挡网站| 亚洲精品日韩av片在线观看| 亚洲欧美精品综合久久99| 黄片wwwwww| 免费黄网站久久成人精品| 久久99热这里只有精品18| 国产激情偷乱视频一区二区| 69av精品久久久久久| 男的添女的下面高潮视频| 欧美极品一区二区三区四区| 男女边吃奶边做爰视频| 人人妻人人澡人人爽人人夜夜 | 91精品国产九色| 国产精品一区二区三区四区久久| 一个人看视频在线观看www免费| 国内精品宾馆在线| 日本午夜av视频| 联通29元200g的流量卡| 亚洲国产精品成人综合色| 久久久亚洲精品成人影院| 亚洲国产精品久久男人天堂| 国产视频首页在线观看| 九九在线视频观看精品| 国产69精品久久久久777片| 国产成年人精品一区二区| 免费看光身美女| 国产精品久久久久久精品电影小说 | 一本久久精品| 国产一区二区三区av在线| 午夜精品一区二区三区免费看| 两个人视频免费观看高清| 两个人的视频大全免费| 国产麻豆成人av免费视频| 国产三级中文精品| 成人鲁丝片一二三区免费| 婷婷色综合大香蕉| 嘟嘟电影网在线观看| 久久精品人妻少妇| 超碰97精品在线观看| 婷婷色综合大香蕉| 日本三级黄在线观看| 精品久久久久久电影网 | 亚洲综合色惰| 91精品一卡2卡3卡4卡| 国产亚洲av嫩草精品影院| 精华霜和精华液先用哪个| 国产欧美日韩精品一区二区| 美女内射精品一级片tv| 午夜福利在线在线| 真实男女啪啪啪动态图| 国产精品爽爽va在线观看网站| 日韩,欧美,国产一区二区三区 | 久久热精品热| 级片在线观看| 少妇裸体淫交视频免费看高清| 国产午夜福利久久久久久| 中文字幕亚洲精品专区| 久久人人爽人人爽人人片va| 日韩成人伦理影院| 丝袜美腿在线中文| 色综合色国产| 国语对白做爰xxxⅹ性视频网站| 欧美激情国产日韩精品一区| 久久久久久久午夜电影| 亚洲,欧美,日韩| 18禁在线无遮挡免费观看视频| 国产精品一区二区性色av| 三级国产精品片| 男女下面进入的视频免费午夜| 老师上课跳d突然被开到最大视频|