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

    凡納濱對(duì)蝦Hox基因及其在早期發(fā)育中表達(dá)模式的研究?

    2015-03-18 08:33:39孫曉晴隗健凱袁劍波張曉軍李富花相建海中國科學(xué)院海洋研究所實(shí)驗(yàn)海洋生物學(xué)重點(diǎn)實(shí)驗(yàn)室山東青島26607中國科學(xué)院大學(xué)北京00049
    關(guān)鍵詞:凡納濱幼體果蠅

    孫曉晴,隗健凱,袁劍波,張曉軍??,李富花,相建海(. 中國科學(xué)院海洋研究所實(shí)驗(yàn)海洋生物學(xué)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 26607;2. 中國科學(xué)院大學(xué), 北京 00049)

    ?

    凡納濱對(duì)蝦Hox基因及其在早期發(fā)育中表達(dá)模式的研究?

    孫曉晴1,2,隗健凱1,2,袁劍波1,張曉軍1??,李富花1,相建海1
    (1. 中國科學(xué)院海洋研究所實(shí)驗(yàn)海洋生物學(xué)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071;2. 中國科學(xué)院大學(xué), 北京 100049)

    Hox基因是決定動(dòng)物形態(tài)多樣性的關(guān)鍵基因,同一類型Hox基因的同源異型結(jié)構(gòu)域(Homeodomain)序列及其功能在進(jìn)化上高度保守。甲殼動(dòng)物是僅次于昆蟲的第二大類節(jié)肢動(dòng)物,但有關(guān)其Hox基因研究在深度和廣度上都遠(yuǎn)遠(yuǎn)不如昆蟲。本研究首先通過RNA-Seq測序技術(shù)和生物信息學(xué)方法分析發(fā)現(xiàn)在凡納濱對(duì)蝦(Litopenaeusvannamei)中存在13種Hox基因,分別為Lvlab、Lvpb、Lvhox3、Lvdfd、Lvscr、Lvftz、Lvantp、Lvubx、LvabdA、LvabdB、Lvmsx、Lvmnx和Lvdll;然后根據(jù)RNA-Seq數(shù)據(jù),克隆獲得了這13條基因的開放閱讀框(ORF)的cDNA序列,并對(duì)其序列結(jié)構(gòu)、早期發(fā)育表達(dá)模式及同源關(guān)系進(jìn)行了初步分析。序列分析顯示這些Hox基因的結(jié)構(gòu)非常保守,同源比對(duì)顯示同一類型的Hox基因的氨基酸序列在不同物種之間具有較高的保守性;表達(dá)模式分析結(jié)果發(fā)現(xiàn),不同Hox基因在早期發(fā)育過程中的表達(dá)模式不同,從囊胚期開始,按時(shí)序表達(dá)的先后順序基本為Lvmsx、Lvlab、Lvpb、Lvhox3、Lvdfd、Lvscr、Lvftz、Lvmnx、Lvdll、Lvantp、Lvubx、LvabdA、LvabdB。根據(jù)各Hox基因表達(dá)時(shí)期對(duì)應(yīng)的對(duì)蝦胚胎和幼體發(fā)育情況,初步推測了這些基因在凡納濱對(duì)蝦的發(fā)育過程中可能的功能。

    Hox基因;凡納濱對(duì)蝦;甲殼動(dòng)物;早期發(fā)育;表達(dá)模式

    Hox基因是后生動(dòng)物中普遍存在的一類含有同源框(Homeobox)的基因,它們都具有一段約180 bp的可轉(zhuǎn)錄出約60個(gè)氨基酸的同源異型結(jié)構(gòu)域(Homeodomain)保守蛋白的序列。這類基因在基因組上成簇存在,同一類型Hox基因在進(jìn)化上高度保守。Hox基因調(diào)控生物體軸形成,參與動(dòng)物早期胚胎發(fā)育,它編碼的Homeodomain蛋白是一種轉(zhuǎn)錄調(diào)控因子,通過與靶基因序列特異結(jié)合而實(shí)現(xiàn)對(duì)生物發(fā)育的調(diào)控[1-2]。1984年,Hox基因在果蠅(Drosophilamelanogaster)中首次被發(fā)現(xiàn)[3-4],逐漸成為人們研究生物體形態(tài)進(jìn)化的焦點(diǎn)。Hox基因若發(fā)生突變,則會(huì)產(chǎn)生同源異型表型,即同源異型突變。例如果蠅的BX-C基因簇中的Ubx發(fā)生突變后,第三胸節(jié)可變?yōu)榕c第2胸節(jié)類似的結(jié)構(gòu),原有的平衡棒被一對(duì)小翅所代替,這就是著名的果蠅雙胸突變體[5]。因此,通過調(diào)控Hox基因的轉(zhuǎn)錄表達(dá),可以達(dá)到控制生物形態(tài)發(fā)育的目的,這對(duì)研究生物體的形態(tài)發(fā)育和進(jìn)化有重要意義。

    根據(jù)大量節(jié)肢動(dòng)物的Hox基因序列以及表達(dá)模式的研究報(bào)道,現(xiàn)在普遍認(rèn)為節(jié)肢動(dòng)物一套完整的Hox基因簇(Hoxgene clusters)可能由laial(lab)、proboscipedia(pb)、Hox3/zen、Deformed(Dfd)、Sexcombsreduced(Scr)、fushitarazu(ftz)、Antennapedia(Antp)、Ultrabithorax(Ubx)、abdominal-A(abdA) 和Abdominal-B(abdB) 10個(gè)基因所組成[6-8]。幾乎所有的節(jié)肢動(dòng)物都具有這一套Hox基因簇,這些基因之間通過復(fù)雜的相互作用造就了節(jié)肢動(dòng)物的形態(tài)多樣性。目前,在節(jié)肢動(dòng)物中,有關(guān)Hox基因的研究大量集中于昆蟲,特別是果蠅的Hox基因結(jié)構(gòu)和功能研究已經(jīng)非常深入,各Hox基因的功能和表達(dá)調(diào)控已經(jīng)基本清楚。然而在甲殼動(dòng)物中,Hox基因的研究還很缺乏。雖然甲殼動(dòng)物與昆蟲同屬節(jié)肢動(dòng)物門,但是它們的形態(tài)特征、遺傳發(fā)育及生活習(xí)性有著非常大的差異,因此盡管它們共用一套同源的Hox基因來調(diào)控形態(tài)發(fā)育,但是Hox基因間的調(diào)控表達(dá)方式卻很可能大不相同,從而在不同的體節(jié)產(chǎn)生不同功能的附肢來適應(yīng)不同的生活環(huán)境。在甲殼動(dòng)物的代表種類-對(duì)蝦中,目前尚無有關(guān)Hox基因的報(bào)道。

    對(duì)蝦的發(fā)育過程與果蠅等昆蟲及其他節(jié)肢動(dòng)物有很大的不同,屬于增節(jié)變態(tài)發(fā)育。在無節(jié)幼體階段身體是不分節(jié)的,但已開始頭胸發(fā)育并產(chǎn)生3對(duì)附肢;從無節(jié)幼體后期開始腹部出現(xiàn)分節(jié);至溞狀幼體和糠蝦幼體,初期的3對(duì)附肢逐漸演變成口器、觸角等器官,體節(jié)增多,身體變長并伴隨著附肢的增長和功能分化(例如胸部附肢稱為步足,行使行走功能;腹部附肢稱為游泳足,行使游動(dòng)功能),同時(shí)在發(fā)育過程中一直伴隨著蛻皮現(xiàn)象[9]。對(duì)蝦類的體節(jié)發(fā)育方式復(fù)雜,相關(guān)研究很少。本研究選擇凡納濱對(duì)蝦(L.vannamei)為實(shí)驗(yàn)材料,通過用生物信息學(xué)方法對(duì)其轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行分析篩選,對(duì)相關(guān)Hox基因的開放閱讀框的cDNA序列進(jìn)行分析和克隆驗(yàn)證,同時(shí)與果蠅的Hox基因進(jìn)行比對(duì)分析,并初步分析這些基因在對(duì)蝦早期發(fā)育中的表達(dá)模式,以便嘗試了解Hox基因在對(duì)蝦類發(fā)育過程中的作用。這些研究將豐富甲殼動(dòng)物Hox基因研究,為對(duì)蝦形態(tài)發(fā)育和進(jìn)化研究奠定基礎(chǔ)。

    1 材料和方法

    1.1Hox基因序列的獲得和生物信息學(xué)分析

    在前期研究中,本課題組開展了凡納濱對(duì)蝦不同發(fā)育時(shí)期(胚胎、無節(jié)幼體、溞狀幼體、糠蝦幼體和仔蝦)的RNA-Seq測序,共得到約26 Gb轉(zhuǎn)錄組數(shù)據(jù),經(jīng)組裝得到66 815條單基因簇(Unigenes)[10]。從Fruitfly Base(http://flybase.org/)中下載了果蠅的2個(gè)Hox基因簇ANT-C和BX-C共8個(gè)Hox基因(lab、pb、Dfd、Scr、Antp、Ubx、abdA和abdB)的cDNA全長序列,將這些序列與凡納濱對(duì)蝦的轉(zhuǎn)錄組數(shù)據(jù)(SRR1039534)進(jìn)行比對(duì),得到同源序列數(shù)據(jù)。同時(shí)對(duì)轉(zhuǎn)錄組組裝得到的Unigenes或Transcripts(轉(zhuǎn)錄本)進(jìn)行篩選,收集注釋為Hox基因的序列。將2方面得到的數(shù)據(jù)合并后進(jìn)行分類和篩選,然后將得到的序列在NCBI去冗余蛋白質(zhì)數(shù)據(jù)庫(Nr)中進(jìn)行BlastX比對(duì)檢索,對(duì)候選Hox基因序列進(jìn)行精確分類,找出各自編碼區(qū)和編碼同源異型結(jié)構(gòu)域區(qū),并進(jìn)行序列間的比對(duì)以及與果蠅各Hox基因的比對(duì)分析。

    1.2 PCR克隆驗(yàn)證

    由于不同的Hox基因表達(dá)的時(shí)期不一樣,根據(jù)RNA-seq結(jié)果,選擇凡納濱對(duì)蝦多個(gè)Hox基因集中表達(dá)的無節(jié)幼體、溞狀幼體和糠蝦幼體3個(gè)發(fā)育時(shí)期的樣品,按照北京全式金生物技術(shù)公司的TransZol Up方法提取3個(gè)樣品的總RNA。以3種總RNA為模板,進(jìn)行反轉(zhuǎn)錄,合成cDNA。

    對(duì)篩選得到的各Hox基因序列,用primer 5.0軟件設(shè)計(jì)系列引物,引物由上海生工生物技術(shù)有限公司合成。對(duì)于一個(gè)基因有多個(gè)Transcripts的情況,選取序列最長或包含同源異型結(jié)構(gòu)域區(qū)最完整的一條轉(zhuǎn)錄本設(shè)計(jì)引物。

    以制備的凡納濱對(duì)蝦的3種cDNA為模板,進(jìn)行PCR擴(kuò)增,獲取不同的Hox基因cDNA片段。PCR反應(yīng)體系為25μL:2×PCR Mix 12.5μL,模板cDNA 1μL,4μmol/L正向引物1μL,4μmol/L反向引物1μL,DNA Taq polymerase 0.25μL,ddH2O 9.25μL。PCR反應(yīng)程序?yàn)椋?5℃變性5min;95℃變性30s,退火(溫度因引物而異)30s,72℃延伸45s,40個(gè)循環(huán);72℃延伸10min。

    PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測,然后將產(chǎn)物片段純化克隆并送到上海桑尼生物科技公司測序,得到的序列與轉(zhuǎn)錄組序列進(jìn)行比對(duì)驗(yàn)證。

    1.3Hox基因在早期發(fā)育中表達(dá)模式

    根據(jù)本課題組的凡納濱對(duì)蝦20個(gè)不同早期發(fā)育時(shí)期的表達(dá)譜RNA-Seq測序數(shù)據(jù)[10],對(duì)經(jīng)篩選鑒定的Hox基因序列在不同的早期發(fā)育時(shí)期的數(shù)據(jù)中進(jìn)行檢索,獲得不同時(shí)期各Hox基因的表達(dá)量信息,并繪制熱圖進(jìn)行聚類分析。根據(jù)各基因的表達(dá)量變化,結(jié)合對(duì)應(yīng)的對(duì)蝦胚胎和幼體發(fā)育時(shí)期的特征,分析各Hox基因在凡納濱對(duì)蝦早期發(fā)育中的表達(dá)模式及可能具有的功能。

    2 結(jié)果與討論

    2.1Hox基因序列

    在凡納濱對(duì)蝦轉(zhuǎn)錄組中,共篩選得到27條Unigenes或Transcripts與Hox基因相關(guān)(見表1)。經(jīng)過BlastX比對(duì)進(jìn)一步分類和分析,這27條Unigenes或Transcripts序列歸結(jié)到13條Hox基因,分別是Lvlab、Lvpb、Lvhox3、Lvdfd、Lvscr、Lvftz、Lvantp、Lvubx、LvabdALvabdB、Lvmsx、Lvmnx和Lvdll。其中Lvlab、Lvpb、Lvdfd、Lvscr、Lvftz、Lvantp、Lvubx、LvabdA、LvabdB、Lvmsx和Lvdll擁有完整的ORF,BlastX比對(duì)覆蓋了整個(gè)Homoedomain區(qū)域;其余的Hox基因序列都只包含部分ORF,經(jīng)BlastX比對(duì)后發(fā)現(xiàn),這些部分編碼區(qū)覆蓋不同Hox基因的全部或部分Homoedomain區(qū)域。

    對(duì)13條Hox基因序列設(shè)計(jì)引物,PCR擴(kuò)增得到目的片段,產(chǎn)物經(jīng)克隆測序驗(yàn)證,各序列與轉(zhuǎn)錄組數(shù)據(jù)中的原始序列保持一致(所有13條序列均已上傳至NCBI,可能存在多個(gè)轉(zhuǎn)錄本的基因,選擇包含生物信息最多的一條上傳)(見表1)。

    2.2Hox基因序列分析

    將凡納濱對(duì)蝦的13條Hox基因的蛋白序列進(jìn)行比對(duì)和結(jié)構(gòu)分析,發(fā)現(xiàn)Hox基因編碼的蛋白結(jié)構(gòu)主要由一個(gè)含YPWM的Motif(基序)(或者類似的4個(gè)氨基酸Motif)和一個(gè)高度保守的Homeodomain組成,并且不同類型的Hox基因序列的Homeodomain區(qū)域具有一定的相似性(見圖1)。這說明不同類型的Hox基因可能來源于同一祖先,經(jīng)過漫長的進(jìn)化,逐漸在某些關(guān)鍵的氨基酸部位發(fā)生突變從而獲得新的功能變成新的基因。很多Hox基因的作用因子通過與YPWM motif的結(jié)合來發(fā)揮調(diào)控作用。有研究認(rèn)為,Hox3和ftz的Hox功能逐漸缺失與YPWM基序的丟失和突變有關(guān)[6],凡納濱對(duì)蝦的Hox基因大多含有這個(gè)基序或者類似的4個(gè)氨基酸序列,這可能與不同Hox基因發(fā)揮不同的功能作用密切相關(guān),其中Lvhox3和Lvftz這2個(gè)基因缺失這種基序或者序列已經(jīng)發(fā)生了很大的變化,其Hox功能是否同樣缺失或改變將在下文討論。

    在NCBI中下載與這13條凡納濱對(duì)蝦Hox基因序列同源的果蠅Hox基因序列。經(jīng)過比對(duì),發(fā)現(xiàn)雖然凡納濱對(duì)蝦與果蠅在進(jìn)化上的親緣關(guān)系甚遠(yuǎn),但是Hox基因的序列及結(jié)構(gòu)保守程度非常高。其中2個(gè)物種的Hox3和ftz的序列相似度較低,尤其是凡納濱對(duì)蝦的Hox3與果蠅的由Hox3衍生的3個(gè)基因zen、zen2和bicoid相比,氨基酸序列的相似度逐步降低。在果蠅中,Hox3衍生的3個(gè)基因序列的YPWM Motif已經(jīng)丟失,并且Homeodomain序列的改變很大,而ftz雖然Homeodomain序列仍舊很保守但已經(jīng)丟失了YPWM序列,可見導(dǎo)致這2個(gè)基因完全丟失Hox基因調(diào)控體節(jié)發(fā)育功能的原因可能是由上文提到的YPWM Motif丟失引起[6]。

    表1 凡納濱對(duì)蝦Hox基因信息

    圖1 凡納濱對(duì)蝦13條Hox基因的蛋白序列的同源性比對(duì)Fig.1 The similarity analysis of 13 Hox protein sequence of L. vannamei

    (D.m:黑腹果蠅;L.v:凡納濱對(duì)蝦。D.m: D.melanogaster; L.v: L.vannamei)圖2 凡納濱對(duì)蝦13條Hox基因的蛋白序列與果蠅的同源性比對(duì)Fig.2 Compare the homology of Hox protein sequences in D.melanogaster and L.vannamei

    2.3Hox基因在凡納濱對(duì)蝦幼體不同時(shí)期中的表達(dá)

    根據(jù)凡納濱對(duì)蝦幼體不同發(fā)育時(shí)期的表達(dá)譜RNA-Seq數(shù)據(jù),將注釋為Hox基因的Unigenes或Transcripts進(jìn)行聚類分析,相似表達(dá)形式的Unigenes或Transcripts聚類到一個(gè)模塊中。見圖3,從這個(gè)熱圖可以看出不同的Unigenes或Transcripts在不同時(shí)期的表達(dá)水平各不相同,但總體看來,Lvlab、Lvpb、Lvhox3、Lvdfd、Lvftz、Lvmsx、Lvdll和Lvmnx主要在無節(jié)幼體及之前的時(shí)期表達(dá)活躍,Lvscr、Lvantp、Lvubx、LvabdA和LvadB則在溞狀幼體及之后的時(shí)期活躍表達(dá)。

    2.4 對(duì)凡納濱對(duì)蝦Hox基因功能的推測

    由于Hox基因是一類高度保守的基因,因此可以通過同一類型的Hox基因在其他節(jié)肢動(dòng)物種類中的作用和該基因在凡納濱對(duì)蝦發(fā)育過程中高表達(dá)時(shí)期對(duì)應(yīng)的發(fā)育特征推測Hox基因在凡納濱對(duì)蝦的發(fā)育過程中的功能。

    從圖3中可以看出,凡納濱對(duì)蝦在受精卵處于囊胚以前的時(shí)期,所有的Hox基因基本都沒有表達(dá)或者非常低量表達(dá);從囊胚期開始隨著細(xì)胞分化而開始陸續(xù)表達(dá)Hox基因。

    (左面:聚類結(jié)果;右面:Hox基因的表達(dá)譜數(shù)據(jù);上面:不同的發(fā)育時(shí)期,由左到右依次是:zygo(受精卵),C2(2細(xì)胞期),C4(4細(xì)胞期),C32(32細(xì)胞期),blas(囊胚期),gast(原腸胚期),Lbe1(肢芽早期),Lbe2(肢芽晚期),Lim1(膜內(nèi)無節(jié)幼體早期),Lim2(膜內(nèi)無節(jié)幼體晚期),N1(無節(jié)幼體1期),N3(無節(jié)幼體3期),N6(無節(jié)幼體6期),Z1(溞狀幼體1期),Z2(溞狀幼體2期),Z3(溞狀幼體3期),M1(糠蝦幼體1期),M2(糠蝦幼體2期),M3(糠蝦幼體3期),P1(仔蝦1期)。Left: The result of cluster analysis; Right: Expression level ofHoxgenes; Up: Different developmental stage ofL.vannamei, from the left to right: zygote, 2 cells, 4 cells, 32cells, blastula, gastrulae, limb bud 1, limb bud 2, nauplius in membrane 1, nauplius in membrane 2, Nauplius 1, Nauplius 3, Nauplius 6, Zoeae 1, Zoeae 2, Zoeae 3, Mysis 1, Mysis 2, Mysis 3,Post-larva 1.)

    圖3 凡納濱對(duì)蝦Hox基因在不同發(fā)育時(shí)期的表達(dá)譜
    Fig.3 TheHoxgenes expression profiling ofL.vannameiin different developmental stages

    從肢芽早期到膜內(nèi)無節(jié)幼體晚期,凡納濱對(duì)蝦先后出現(xiàn)3對(duì)附肢:第一觸角、第二觸角和大顎,并且在胚體前端的腹面中央出現(xiàn)紅褐色的單眼[9]。在這期間,Lvlab、Lvhox3、Lvmsx、Lvdfd、Lvftz、Lvmnx和Lvdll7個(gè)基因先后開始表達(dá),表達(dá)量先升高再降低。Lvlab、Lvdfd和Lvmsx3個(gè)基因的表達(dá)趨勢(shì)一致,而Lvhox3的表達(dá)則表現(xiàn)為與其相反的趨勢(shì)。Lvftz與Lvdfd同時(shí)在膜內(nèi)無節(jié)幼體開始表達(dá),但表達(dá)量上調(diào)和下調(diào)的變化很大。Lvmnx和Lvdll開始表達(dá)的時(shí)期較晚,且表達(dá)量一直不高,但表達(dá)時(shí)期延伸到了后面的溞狀幼體1期(見圖4a~f)。因此,推測Lvlab、Lvmsx、Lvhox3、Lvdfd和Lvftz可能參與調(diào)控凡納濱對(duì)蝦觸角、大顎和單眼的發(fā)育,Lvdll和Lvmnx可能與3對(duì)附肢的后期發(fā)育有關(guān),并且參與無節(jié)幼體期附肢的細(xì)節(jié)特征和功能分化。對(duì)蝦無節(jié)幼體期分為6個(gè)時(shí)期,主要特征是軀體增大,3對(duì)附肢相應(yīng)增長增大,細(xì)節(jié)特征增多,到無節(jié)幼體期6期,出現(xiàn)頭胸甲、2對(duì)小顎和2對(duì)顎足,腹部出現(xiàn)多對(duì)附肢原基,口器和消化器官逐漸形成[9]。Hox基因的表達(dá)表現(xiàn)為孵化前期高表達(dá)的基因如Lvlab、Lvmsx、Lvhox3和Lvftz等開始低量表達(dá)或者不表達(dá),而另外一些Hox基因,如Lvpb、Lvscr、Lvantp開始上調(diào)表達(dá),推測Lvpb、Lvscr、Lvantp可能主要參與口器和消化器官的發(fā)育形成,并且Lvdll和Lvmnx可能也起到重要的調(diào)控作用。

    至溞狀幼體時(shí)期具完整的口器和消化器官,身體開始分節(jié)并且在胸節(jié)長出附肢[9]。從Hox基因的表達(dá)上可以看出,Ubx、LvabdA、LvabdB開始表達(dá)。

    凡納濱對(duì)蝦在糠蝦幼體至仔蝦時(shí)期已具有全部的體節(jié)和附肢,且各附肢的細(xì)節(jié)特征逐步得到完善。這段時(shí)期主要表達(dá)LabdA、LvabdB3個(gè)基因(見圖4g~p)。這3個(gè)基因之間LvabdA和Lvubx的表達(dá)趨勢(shì)一致并且同LvabdB的表達(dá)趨勢(shì)相反,呈此消彼長的變化規(guī)律,推測這3個(gè)基因之間相互作用共同調(diào)節(jié)和完善凡納濱對(duì)蝦胸腹節(jié)及其附肢的發(fā)育。

    2.4.1Lvlab在節(jié)肢動(dòng)物Hox基因的研究中,lab總是最先表達(dá)的一個(gè)基因[6-7]。在果蠅幼體期,lab與觸角、視葉原基、中腸和部分由腦發(fā)出的神經(jīng)發(fā)育密切相關(guān);成體時(shí)表達(dá)則僅限于眼部、觸角及其他頭部區(qū)域的發(fā)育[11]。在甲殼動(dòng)物球鼠婦(Porcelliolaevis)中,lab僅在第二對(duì)觸角表達(dá),調(diào)控觸角發(fā)育[12]。Lvlab主要在凡納濱對(duì)蝦肢芽早期開始大量表達(dá),在膜內(nèi)無節(jié)幼體晚期迅速下降,以后表達(dá)量一直很低。推測Lvlab在凡納濱對(duì)蝦發(fā)育過程中的作用至少有一點(diǎn)與其他節(jié)肢動(dòng)物相同,就是與觸角的發(fā)育密切相關(guān)。

    2.4.2Lvhox3和lvftzHox3和ftz都是非常古老的Hox基因之一,在生物漫長的進(jìn)化史中,這2個(gè)基因的序列、表達(dá)區(qū)域以及功能都發(fā)生了改變。在蜈蚣和蜘蛛等多足類和螯肢類動(dòng)物中,這2個(gè)基因仍具有Hox基因的功能,參與顎的發(fā)育[6];在果蠅等昆蟲中,這2個(gè)基因經(jīng)過復(fù)雜的進(jìn)化歷史,已經(jīng)丟失了Hox-like功能,而分別在胚外組織和神經(jīng)元中獲得了不同的新的功能[6]。在甲殼動(dòng)物中,這2個(gè)基因在水溞(Daphniapulex)胚胎發(fā)育的不同時(shí)期表達(dá)位置不同,最終在大顎節(jié)或大顎節(jié)后部表達(dá),因此研究者認(rèn)為ftz和Hox3在水溞中仍具有Hox基因功能[13];球鼠婦中Hox3僅在大顎肢芽中心表達(dá)[8],藤壺(Sacculinacarcini)中,ftz僅在神經(jīng)系統(tǒng)表達(dá)[14]。在凡納濱對(duì)蝦中,Lvhox3在囊胚期迅速大量增加,原腸胚表達(dá)量達(dá)到最大后迅速下降,在其他時(shí)期的表達(dá)水平都較低。Lvftz從原腸胚期開始表達(dá),在肢芽期達(dá)到最大量表達(dá),無節(jié)幼體時(shí)期及之后微量表達(dá)或不表達(dá)。推測,既然這2個(gè)基因在凡納濱對(duì)蝦胚胎早期、變態(tài)發(fā)育開始發(fā)生前大量表達(dá),就可能仍具有Hox基因調(diào)控動(dòng)物形態(tài)發(fā)育的功能,至于具體調(diào)控胚胎和幼體的哪一部分發(fā)育,以及是否具有調(diào)其他新功能,還有待深入研究。

    (每圖縱坐標(biāo)代表Hox基因的表達(dá)量,橫坐標(biāo)是不同的Hox基因,由左到右依次是:Lvlab、Lvpb、Lvhox3、Lvdfd、Lvscr、Lvftz、Lvmsx、Lvdll、Lvmnx、Lvantp、Lvubx、LvabdA、LvabdB。A~p分別是凡納濱對(duì)蝦不同的發(fā)育時(shí)期各Hox基因的表達(dá)情況:a.囊胚期,b.原腸胚期,c.肢芽早期,d. 肢芽晚期,e. 膜內(nèi)無節(jié)幼體早期,f.膜內(nèi)無節(jié)幼體晚期,g. 無節(jié)幼體1期,h. 無節(jié)幼體3期,i.無節(jié)幼體6期,j.溞狀幼體1期,k. 溞狀幼體2期,l.溞狀幼體3期,m.糠蝦幼體1期,n.糠蝦幼體2期,o.糠蝦幼體3期,p.仔蝦1期。Ordinate: Expression level ofHoxgenes; Abscissa: 13Hoxgene, from the left to right:Lvlab,Lvpb,Lvhox3,Lvdfd,Lvscr,Lvftz,Lvmsx,lvdll,Lvmnx,Lvantp,Lvubx,LvabdA,LvabdB, a~p showed different development stage ofL.vannamei: a. blastula, b. gastrulae, c. limb bud 1, d. limb bud 2, e. nauplius in membrane 1, f. nauplius in membrane 2, g. Nauplius 1, h. Nauplius 3, i. Nauplius 6, j. Zoea 1, k. Zoea 2, l. Zoea 3, m. Mysis 1, n. Mysis 2, o. Mysis 3, p. Post-larva 1.)

    圖4 凡納濱對(duì)蝦胚胎及幼體發(fā)育的時(shí)期各Hox基因的表達(dá)變化
    Fig.4 The dynamic expressions of 13Hoxgenes in difference early development stage ofL.vannamei

    2.4.3LvMsxMsx從腔腸動(dòng)物到哺乳動(dòng)物均有表達(dá),是一類高度保守的Hox基因。無脊椎動(dòng)物只有單一的MSH/MSX,對(duì)背腹中樞神經(jīng)系統(tǒng)的形成起著至關(guān)重要的作用,該功能在進(jìn)化上高度保守[15-17]。凡納濱對(duì)蝦的Lvmsx從受精卵即出現(xiàn),至原腸胚時(shí)期開始大量表達(dá),幾乎覆蓋隨后的所有發(fā)育時(shí)期,而高度表達(dá)時(shí)期主要集中在肢芽期到膜內(nèi)無節(jié)幼體時(shí)期,根據(jù)前期文獻(xiàn),推測在Lvmsx的功能很多,可能對(duì)凡納濱對(duì)蝦的神經(jīng)系統(tǒng)的發(fā)育起著至關(guān)重要的作用。在凡納濱對(duì)蝦中,疑似存在多個(gè)Lvmsx轉(zhuǎn)錄本,推測不同的轉(zhuǎn)錄本作為不同的轉(zhuǎn)錄因子或信號(hào)分子調(diào)控不同部分的神經(jīng)發(fā)育[18]。

    2.4.4Lvdfd有研究發(fā)現(xiàn)在球鼠婦中,Dfd僅在大顎節(jié)和小顎節(jié)表達(dá),但是在大顎節(jié)中心受到抑制,而這個(gè)區(qū)域正是Hox3表達(dá)位置,因此推測Lvdfd受Hox3的表達(dá)抑制[12]。在凡納濱對(duì)蝦中,Lvdfd從原腸胚期開始連續(xù)表達(dá),在膜內(nèi)無節(jié)幼體早期達(dá)到最大后開始迅速下降,之后低量表達(dá)。在此期間形成的無節(jié)幼體的身體將來主要發(fā)育為對(duì)蝦的頭部,因此推測,Lvdfd的功能可能是負(fù)責(zé)頭部大小顎的形成發(fā)育。結(jié)合Lvhox3的表達(dá)趨勢(shì)及其在鼠婦中的作用,凡納濱對(duì)蝦中Lvdfd和lvhox3之間可能也存在著相互抑制的作用,共同調(diào)控大顎小顎的發(fā)育。

    2.4.5LvdllDll在節(jié)肢動(dòng)物中的最主要最普遍的功能是調(diào)控附肢的發(fā)育[6]。在昆蟲中,Dll在腹部的附肢發(fā)育時(shí)表達(dá)受到Ubx/abdA的抑制而使腹部無附肢[19-20];而在多足綱和甲殼綱動(dòng)物中,Dll在腹部的表達(dá)和及其對(duì)附肢的調(diào)控并不受Ubx或abdA抑制[21]。在果蠅中,Dll第二個(gè)功能還與觸角的發(fā)育密切相關(guān),由于果蠅的觸角具有聽覺和嗅覺的功能,研究證實(shí),Dll還與嗅覺神經(jīng)系統(tǒng)發(fā)育相關(guān)[22-23]。Lvdll在凡納濱對(duì)蝦肢芽早期表達(dá)迅速上調(diào),膜內(nèi)無節(jié)幼體早期達(dá)到最大量表達(dá),之后表達(dá)量迅速下降,到溞狀幼體時(shí)期以后一直微量表達(dá)。推測Lvdll主要參與調(diào)控凡納濱對(duì)蝦早期3對(duì)附肢的發(fā)育,可能遠(yuǎn)程調(diào)控胸腹部附肢的形成,由于胸部和腹部分別是功能分明的步足和游泳足,因此Lvdll的表達(dá)并不受完全Ubx或LvabdA的抑制,更有可能是一起共表達(dá)調(diào)控胸腹部附肢的發(fā)育,決定附肢的種類和功能。

    2.4.6LvmnxMnx在人類基因中被稱為HB9或HLXB9,對(duì)哺乳動(dòng)物的胰腺發(fā)育和運(yùn)動(dòng)神經(jīng)元的分化合并至關(guān)重要,突變可導(dǎo)致Currarino綜合征[24]。無脊椎動(dòng)物雖然缺乏胰腺,但是在原始的胚胎內(nèi)胚層仍有該基因的表達(dá)[25]。在果蠅中,Mnx還是運(yùn)動(dòng)神經(jīng)細(xì)胞遷移、軸突生長和決定運(yùn)動(dòng)神經(jīng)元亞型種類所必需的調(diào)控因子。在凡納濱對(duì)蝦中,Lvmnx從肢芽早期開始有微量表達(dá),在整個(gè)胚胎及胚后發(fā)育時(shí)期的表達(dá)量一直很低。推測Lvmnx在凡納濱對(duì)蝦中的功能有2種可能,一種是促進(jìn)組織肝胰臟的發(fā)育,另一種可能與運(yùn)動(dòng)神經(jīng)的發(fā)育分化密切相關(guān)。

    2.4.7Lvpbpb在是昆蟲特化結(jié)構(gòu)唇顎發(fā)育的重要調(diào)控基因。在果蠅中主要在胚胎和成體的唇節(jié)顎節(jié)表達(dá),并且分別與Dfd和Scr相互結(jié)合作用決定成體果蠅顎須和唇的特化[26];在螯肢類和多足類動(dòng)物中,Lvpb在從大小顎到步足都有表達(dá)[6];而在球鼠婦中僅在第二觸角表達(dá)[12]。Lvpb在凡納濱對(duì)蝦膜內(nèi)無節(jié)幼體早期開始表達(dá),這正是口形成的時(shí)期,但Lvpb一直到仔蝦所有表達(dá)時(shí)期的有表達(dá),只是表達(dá)量都很低。推測Lvpb可能與對(duì)蝦口器的發(fā)育密切相關(guān),但很可能還有其他的功能。

    2.4.8LvscrScr在不同的節(jié)肢動(dòng)物中的表達(dá)區(qū)域不一樣,且功能呈現(xiàn)多樣化。在果蠅中,可誘導(dǎo)唾液腺的發(fā)育,在家蠶中誘導(dǎo)絲腺和唾液腺的發(fā)育,在其他昆蟲的一些特化附肢中也具有重要作用[6]。在馬利筋長蝽(Oncopeltusfasciatus)中,Scr與pb共同作用決定唇的發(fā)育,與Dfd共同作用決定小顎的發(fā)育[27]。在甲殼動(dòng)物中,Scr與Ubx相互作用調(diào)控顎肢的發(fā)育,兩者動(dòng)態(tài)表達(dá)產(chǎn)生不同的顎肢[7],如在球鼠婦和克原氏螯蝦(Procambarusclarkii)中,Scr在早期胚胎發(fā)育的第一對(duì)小顎中表達(dá)[28-29]。在凡納濱對(duì)蝦中,Lvscr主要在無節(jié)幼體時(shí)期和溞狀幼體早期表達(dá),表達(dá)趨勢(shì)與Lvpb的表達(dá)一致。由于在果蠅等昆蟲中Scr對(duì)pb具有活化的正調(diào)控作用,推測Lvscr可能與Lvpb共同作用協(xié)調(diào)表達(dá)調(diào)控對(duì)蝦口器顎的發(fā)育。

    2.4.9LvantpLvantp已在3種甲殼動(dòng)物中研究報(bào)道:在鹵蟲(Artemiafranciscan)中,Lvantp表達(dá)區(qū)域從第一對(duì)小顎節(jié)后部開始延伸到全部胸節(jié)[29];在球鼠婦中則從第二對(duì)小顎節(jié)開始延伸到整個(gè)胸部[30];在克原氏螯蝦中,Antp的表達(dá)在腹部的表達(dá)較弱,在第二對(duì)小顎節(jié)附肢和3對(duì)螯肢中最強(qiáng)[28]。因此,Antp主要參與決定甲殼動(dòng)物小顎和胸部附肢的發(fā)育。在凡納濱對(duì)蝦中,Antp在無節(jié)幼體時(shí)期表達(dá)量開始快速增加,至溞狀幼體時(shí)期開始下降并趨于穩(wěn)定低量表達(dá)。凡納濱對(duì)蝦自無節(jié)幼體時(shí)期末期開始出現(xiàn)多對(duì)附肢和頭胸甲,在溞狀幼體時(shí)期開始分節(jié)并且胸節(jié)長出附肢。由此推測lvantp主要參與調(diào)控凡納濱對(duì)蝦小顎、胸節(jié)及其附肢的發(fā)育。

    2.4.10Lvubx在果蠅中,Ubx決定后部胸節(jié)和腹節(jié)發(fā)育,Ubx還通過抑制Dll的表達(dá)抑制腹部附肢的發(fā)育[19-20]。在鹵蟲中,Ubx/abdA的表達(dá)與胸腹部附肢的形態(tài)發(fā)育關(guān)系密切,Ubx/abdA在不同物種間的動(dòng)態(tài)表達(dá),誘導(dǎo)決定了各對(duì)附肢的功能,例如捕食、游泳、爬行或戰(zhàn)斗[31]。在凡納濱對(duì)蝦中,Lvubx從無節(jié)幼體時(shí)期開始一直到糠蝦幼體和仔蝦時(shí)期都有表達(dá),這段時(shí)間主要發(fā)生增節(jié)變態(tài)發(fā)育,出現(xiàn)胸節(jié)和腹節(jié)及附肢并且各附肢的分工明確。因此,推測Lvubx在凡納濱對(duì)蝦中的功能作用主要是與LvabdA和LvabdB共同作用在胸腹部決定體節(jié)和附肢的形態(tài)發(fā)育。

    2.4.11LvabdAabdA在節(jié)肢動(dòng)物中的主要功能作用是調(diào)節(jié)腹部附肢的發(fā)育。abdA在已研究的3種甲殼動(dòng)物表達(dá)各異:在鹵蟲中,自第二胸節(jié)的神經(jīng)節(jié)開始表達(dá),覆蓋整個(gè)胸部,與Ubx的表達(dá)區(qū)域重疊[29];在球鼠婦中,abdA在腹部表達(dá)跟Ubx表達(dá)區(qū)域無重疊[30];在克原氏螯蝦胚胎發(fā)育后期,abdA與Ubx表達(dá)位置大多無重疊[28],因此導(dǎo)致這3種甲殼動(dòng)物的胸腹部附肢形態(tài)各異。凡納濱對(duì)蝦的lvabdA在無節(jié)幼體時(shí)期開始連續(xù)表達(dá),呈上調(diào)趨勢(shì),在糠蝦幼體時(shí)期開始穩(wěn)定表達(dá),但表達(dá)量不高。表達(dá)趨勢(shì)同Lvubx的表達(dá)呈現(xiàn)一致,推測LvabdA與Lvubx基因相互作用共同調(diào)控對(duì)蝦胸部和腹部體節(jié)及附肢的發(fā)育。

    2.4.12LvabdB在昆蟲發(fā)育的過程中,AbdB不論在基因簇上的位置還是表達(dá)時(shí)間先后順序都是最后一個(gè),決定昆蟲腹部后部體節(jié)和生殖腺的發(fā)育[7]。而在鹵蟲中AbdB基因只在位于胸節(jié)后部腹節(jié)前部之間的生殖節(jié)中表達(dá),AbdB對(duì)生殖器官的發(fā)育具有重要的調(diào)控作用且該功能非常保守[29]。凡納濱對(duì)蝦LvabdB基因從無節(jié)幼體時(shí)期開始連續(xù)上調(diào)表達(dá),但是表達(dá)量都很低,與Lvubx/LvabdA呈現(xiàn)一種此消彼長的表達(dá)趨勢(shì)。推測,AbdB在凡納濱對(duì)蝦發(fā)育早期主要是與Lvubx和LvabdA共同參與調(diào)控體節(jié)和附肢形成,調(diào)控生殖器官發(fā)育的功能可能在仔蝦后期發(fā)揮作用。

    本研究中發(fā)現(xiàn)Lvlab、lvmsx和lvubx擁有多個(gè)轉(zhuǎn)錄本,這些轉(zhuǎn)錄本開始表達(dá)的時(shí)期、表達(dá)量及表達(dá)上調(diào)下調(diào)的趨勢(shì)各不相同(見圖5)。不同的轉(zhuǎn)錄本作為轉(zhuǎn)錄調(diào)控因子可能調(diào)控不同的靶基因或?qū)ν话谢虻恼{(diào)控方式不同,從而引起不同的效應(yīng),例如造成體節(jié)特異的特征等[18]。由于目前缺乏對(duì)蝦全基因組參考序列,目前還不能確定這些轉(zhuǎn)錄本是來自于可變剪切還是基因家族,隨著研究的深入和對(duì)蝦基因組序列數(shù)據(jù)的完善,這些轉(zhuǎn)錄本的來源和功能可以獲得更深入的闡釋。

    圖5 凡納濱對(duì)蝦多轉(zhuǎn)錄本的Hox基因在不同發(fā)育時(shí)期的表達(dá)Fig.5 The expressions of Hox genes with more than one transcripts in difference early development stage of L. vannamei

    2.5 凡納濱對(duì)蝦Hox基因的同源關(guān)系分析

    現(xiàn)在普遍認(rèn)為所有的Hox基因擁有共同的祖先,經(jīng)歷多次基因復(fù)制或基因組加倍事件后逐漸產(chǎn)生具有不同功能的Hox基因。將13條Hox基因的Homeodomain序列進(jìn)行比對(duì)和聚類分析,發(fā)現(xiàn)在凡納濱對(duì)蝦中,按照節(jié)肢動(dòng)物Hox基因排列順序相鄰的2個(gè)Hox基因首先聚為一類,其次再和相隔基因聚在一起,如Lvpb和Lvhox3聚類后再與Lvlab聚類,Lvubx和LvabdA聚類后與Lvftz和Lvantp聚類然后再與Lvscr聚類(見圖6)。這也在一定程度上支持了Hox基因擁有同一祖先的假設(shè)。聚類和同源性分析研究也發(fā)現(xiàn)存在于凡納濱對(duì)蝦中的3個(gè)Hox基因Lvmsx、Lvmnx和Lvdll與其他Hox基因的關(guān)系較遠(yuǎn),自成一支,這3個(gè)基因或許都屬于比較古老的Hox基因。

    圖6 凡納濱對(duì)蝦Hox基因的同源關(guān)系Fig.6 The homology of Hox homoedomain sequence of L. vannamei

    3 結(jié)語

    作為生物發(fā)育過程中至關(guān)重要的一類基因,Hox基因在漫長的進(jìn)化中,各自的序列和功能既保守又動(dòng)態(tài)變化,它們彼此之間相互作用,在時(shí)空順序上組成了一個(gè)復(fù)雜的調(diào)控網(wǎng)絡(luò),共同決定生物體的組織發(fā)育和形態(tài)特征等。在凡納濱對(duì)蝦中首次發(fā)現(xiàn)了13種Hox基因,序列結(jié)構(gòu)分析和同源比對(duì)結(jié)果顯示這些Hox基因的Homeobox序列和Homeodomain都很保守,不僅是同一類型Hox基因與遠(yuǎn)緣物種-果蠅間保守,并且凡納濱對(duì)蝦的不同類型的Hox基因之間也具有一定的同源性,支持Hox來源域同一祖先的說法。通過對(duì)這些Hox基因在對(duì)蝦早期發(fā)育不同時(shí)期的表達(dá)模式的分析,認(rèn)為這些Hox基因全部參與了凡納濱對(duì)蝦的早期發(fā)育,且按一定的時(shí)空順序先后開始表達(dá)。推測Lvlab、Lvmsx、Lvhox3、Lvdfd和Lvftz可能參與調(diào)控凡納濱對(duì)蝦觸角、大顎的發(fā)育,Lvhox3和Lvftz2個(gè)基因可能仍具有Hox基因的功能。Lvdll和Lvmnx可能與3對(duì)附肢的后期發(fā)育有關(guān),并且與Lvpb、Lvscr、Lvantp共同參與口器和消化器官的發(fā)育形成。同時(shí)發(fā)現(xiàn)這些Hox基因之間可能存在相互抑制共同作用的情況,如Lvubx、LvabdA和LvabdB3個(gè)基因之間相互作用共同調(diào)節(jié)和完善凡納濱對(duì)蝦胸腹節(jié)及其附肢的發(fā)育。根據(jù)不同時(shí)期凡納濱對(duì)蝦的生長發(fā)育情況,推測了Hox基因在這些時(shí)期可能行使的功能,這些結(jié)果為今后的對(duì)蝦體節(jié)發(fā)育和Hox基因的功能驗(yàn)證研究提供了參考,為深入研究甲殼動(dòng)物的進(jìn)化及解釋其形態(tài)多樣性提供了更多的信息。

    [1] Banerjee B S, Ferlanti E S, Ryan J F, et al. The Homeodomain Resource: sequences, structures and genomic information[J]. Nucleic Acids Res, 1999, 27(1): 336-337.

    [2] Banerjee B S, Moreland T, Hsu B J, et al. The Homeodomain Resource: 2003 update[J]. Nucleic Acids Res, 2003, 31(1): 304-306.

    [3] McGinnis W, Levine MS, Hafen E, et al. A conserved DNA sequence in homoeotic genes of theDrosophilaAntennapediaand bithorax complexes[J]. Nature, 1984, 308(5958): 428-433.

    [4] Scott M P, Weiner A J. Structural relationships among genes that control development: sequence homology between theAntennapedia,Ultrabithorax, andfushitarazuloci ofDrosophila[J]. Proc Natl Acad Sci USA, b, 1984, 81(13): 4115-4119.

    [5] Lewis E B. A gene complex controlling segmentation inDrosophila[J]. Nature, 1978, 276(5688): 565-570.

    [6] Hughes C L, Kaufman T C.Hoxgenes and the evolution of the arthropod body plan[J]. Evolution & Development, 2002, 4(6): 459-499.

    [7] Hou L.Hoxgenes and study ofHoxgenes in crustacean[J]. Chinese Journal of Oceanology and Limnology,2004,4(22): 392-398.

    [8] Michalis A, Anastasios P Zacharias K, et al. Evolution of new appendage types by gradual changes inHoxgene expression-the case of crustacean maxillipeds[J]. Palaeodiversity, 2010, 3: 141-155.

    [9] 趙法箴. 海洋水產(chǎn)研究資料, 對(duì)蝦幼體發(fā)育形態(tài)[M]. 北京: 農(nóng)業(yè)出版社, 1965.

    [10] Wei J, Zhang X, Yu Y, et al. Comparative transcriptomic characterization of the early development in Pacific white shrimpLitopenaeusvannamei[J]. PLoS ONE, 2014, 9(9): e106201.

    [11] Merrill V K, Diederich R J, Turner F R, et al. A genetic and developmental analysis of mutations inlabial, a gene necessary for proper head formation inDrosophilamelanogaster[J]. Dev Biol, 1989, 135: 376-391.

    [12] Abzhanov A, Kaufman T C. Homeotic genes and the arthropod head: expression patterns of thelabial,proboscipedia, andDeformedgenes in crustaceans and insects[J]. Proc Natl Acad Sci USA, a, 1999, 96: 10224-10229.

    [13] Papillon D, Telford M J. Evolution ofHox3 andftzin arthropods: insights from the crustaceanDaphniapulex[J]. Dev Genes Evol, 2007, 217(4): 315-322.

    [14] Mouchel V E, Blin M Rigolot C, et al. Expression of a homologue of thefushitarazu(ftz) gene in a cirripede crustacean[J]. Evol Dev, 2002, 4: 76-85.

    [15] Davidson D. The function and evolution ofMsxgenes pointers and paradoxes[J]. Trends Genet, 1995, 11(10): 405-411.

    [16] Cornell R A, Ohlen T V.Vnd/nkx,ind/gsh, andmsh/msx: conserved regulators of dorsoventral neural patterning?[J]. Curr Opin Neurobiol, 2000, 10(1): 63-71.

    [17] Ramos C, Robert B.Msh/Msxgene family in neural development[J]. Trends Genet, 2005, 21(11): 624-632.

    [18] 田曉軒, 謝強(qiáng), 卜文俊. 乳草長蝽Ubx基因克隆及多轉(zhuǎn)錄本分析[J]. 昆蟲學(xué)報(bào), 2011, 54(4): 390-396.

    [19] Vachon G, Cohen B, Pfeiflec, et al. Homeotic genes of the bithorax complex repress limb development in the abdomen of theDrosophilaembryo through the target geneDistal-Less[J]. Cell, 1992, 71(3): 437-450.

    [20] Warren R W, Nagy L, Seleque J. et al. Evolution of homeotic gene regulation and function in flies and butterflies[J]. Nature, 1994, 372: 458-461.

    [21] Panganiban G, Sebring A, Nagy L, et al. The development of crustacean limbs and the evolution of arthropods[J]. Science, 1995, 270: 1363-1366.

    [22] Panganiban G, Rubenstein J L. Developmental functions of theDistal-less/Dlxhomeobox genes[J]. Development, 2002, 129(19): 4371-4386.

    [23] Plavicki J, Mader S, Pueschel E, et al. Homeobox genedistal-lessis required for neuronal differentiation and neurite outgrowth in theDrosophilaolfactory system[J]. Proc Natl Acad Sci USA, 2012, 109(5): 1578-1583.

    [24] Hagan D M, Ross A J, Strachan T, et a1. Mutation analysis and embryonic expression of theHLXB9 Currarino syndrome gene[J]. Am J Hum Genet, 2000, 66(5): 1504-1515.

    [25] Wendik B, Maier E, Meyer D. Zebrafishmnxgenes in endocrine and exocrine pancreas formation[J]. Dev Biol, 2004, 268(2): 372-383.

    [26] Hughes C L, Kaufman T C. RNAi analysis ofDeformed,proboscipediaandSexcombsreducedin the milkweed bugOncopeltusfasciatus: novel roles forHoxgenes in the hemipteran head[J]. Development, 2000, 127: 3683-3694.

    [27] Abzhanov A, Kaufman T C. Novel regulation of the homeotic geneScrassociated with a crustacean leg-to-maxilliped appendage transformation[J]. Development, b, 1999, 126: 1121-1128.

    [28] Abzhanov A, Kaufman T C. Embryonic expression patterns of theHoxgenes of the crayfishProcambarusclarkii(Crustacea, Decapoda) [J]. Evol Dev, b, 2000, 2: 271-283.

    [29] Averof M, Akam M.Hoxgenes and the diversification of insect and crustacean body plans[J]. Nature, 1995, 376: 420-423.

    [30] Abzhanov A, Kaufman T C. Crustacean (malacostracan)Hoxgenes and the evolution of the arthropod trunk[J]. Development, a, 2000, 127: 2239-2249.

    [31] Averof M, Patel N H. Crustacean appendage evolution associated with changes inHoxgene expression[J]. Nature, 1997, 388: 682-686.

    責(zé)任編輯 高 蓓

    HoxGenes and Their Expression Pattern in Early Development ofLitopenaeusvannamei

    SUN Xiao-Qing1,2, WEI Jian-Kai1,2, YUAN Jian-Bo1, ZHANG Xiao-Jun1, LI Fu-Hua1, XIANG Jian-Hai1

    (1. The Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China)

    Hoxgenes involve in early embryonic development of animals and regulate the formation of body segments, tissues and organs. They have been widely studied to explain the morphological development and evolution in insects, especially inDrosophila, but rarely in crustaceans, which is the second largest group of arthropods. In this study, by RNA-seq sequencing and bioinformatics analysis, we found 13Hoxgenes in Pacific white shrimp (Litopenaeusvannamei):Lvlab,lvpb,lvhox3,lvdfd,lvscr,lvftz,lvmsx,lvdll,lvmnx,lvantp,lvubx,lvabdAandlvabdB. Molecular clone and sequence analysis showed the 13Hoxgenes are conserved in sequence structure. By analysis the expression pattern of these genes at the early development stages ofL.vannamei, we found that their expression indicated a dynamic change at different development stages and acoording a temporal order based on different function. Our findings suggested that there is a complicated interaction with theHoxgenes together to determine the morphology of the shrimp. This study provided important information for the morphological development and evolution analysis in crustaceans.

    Hoxgene;Litopenaeusvannamei; crustacean; early development; expression pattern

    國家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2012AA092205;2012AA10A404);國家自然科學(xué)基金項(xiàng)目(31172396;31272683)資助

    2014-09-30;

    2014-12-12

    孫曉晴(1985-),女,碩士生。E-mail:xiaoqingsun@126.com

    ?? 通訊作者: E-mail: xjzhang@qdio.ac.cn

    Q959.223

    A

    1672-5174(2015)08-052-11

    10.16441/j.cnki.hdxb.20140323

    猜你喜歡
    凡納濱幼體果蠅
    果蠅也會(huì)“觸景傷身”
    小果蠅大貢獻(xiàn)
    黃海近岸潮汐鋒海域蟹類幼體日間垂直遷移特征?
    果蠅遇到危險(xiǎn)時(shí)會(huì)心跳加速
    探究“胚是新植物的幼體”實(shí)驗(yàn)創(chuàng)新
    奇異海蟑螂Ligia exotica幼體的生長特征研究
    小果蠅助力治療孤獨(dú)癥
    凡納濱對(duì)蝦白斑綜合征病毒防治研究進(jìn)展(一)
    凡納濱對(duì)蝦與點(diǎn)帶石斑魚的混養(yǎng)模式
    凡納濱對(duì)蝦CTSL基因與生長相關(guān)的SNP位點(diǎn)的特征
    99久久无色码亚洲精品果冻| 777久久人妻少妇嫩草av网站| 久久久久久人人人人人| 18禁国产床啪视频网站| 亚洲精华国产精华精| 国产区一区二久久| 久久精品亚洲精品国产色婷小说| 大型黄色视频在线免费观看| 日韩三级视频一区二区三区| 一级作爱视频免费观看| 一级毛片高清免费大全| 老司机靠b影院| 美女黄网站色视频| 色精品久久人妻99蜜桃| 亚洲国产精品久久男人天堂| 小说图片视频综合网站| 岛国视频午夜一区免费看| 一个人观看的视频www高清免费观看 | 免费高清视频大片| 麻豆成人午夜福利视频| 精品久久久久久久人妻蜜臀av| 国产一级毛片七仙女欲春2| 国产精品久久久人人做人人爽| 叶爱在线成人免费视频播放| 91麻豆av在线| 欧美高清成人免费视频www| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 国产熟女午夜一区二区三区| 可以在线观看毛片的网站| 亚洲中文字幕一区二区三区有码在线看 | 18禁国产床啪视频网站| 精品欧美国产一区二区三| 中文在线观看免费www的网站 | 麻豆成人午夜福利视频| 看黄色毛片网站| 91字幕亚洲| 99在线人妻在线中文字幕| 一本一本综合久久| 亚洲在线自拍视频| 一进一出抽搐gif免费好疼| 久久久久久久久久黄片| 亚洲欧美日韩高清专用| 看黄色毛片网站| 97碰自拍视频| 精品国产乱码久久久久久男人| 欧美色欧美亚洲另类二区| 看黄色毛片网站| 此物有八面人人有两片| 叶爱在线成人免费视频播放| 欧美另类亚洲清纯唯美| 欧美成人午夜精品| 老汉色av国产亚洲站长工具| 国产成人精品无人区| 特大巨黑吊av在线直播| 国语自产精品视频在线第100页| 99re在线观看精品视频| 在线观看午夜福利视频| 久久久精品大字幕| 神马国产精品三级电影在线观看 | 国产精品国产高清国产av| 18禁黄网站禁片免费观看直播| 男女之事视频高清在线观看| 亚洲欧美激情综合另类| 超碰成人久久| 日韩欧美一区二区三区在线观看| 麻豆国产av国片精品| 国内精品久久久久精免费| 999久久久精品免费观看国产| 国产私拍福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 久久久久久久久免费视频了| 97碰自拍视频| 大型av网站在线播放| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 丰满人妻一区二区三区视频av | 99精品欧美一区二区三区四区| 给我免费播放毛片高清在线观看| 一级毛片精品| 男人舔女人的私密视频| 国产高清有码在线观看视频 | 欧美日本亚洲视频在线播放| 欧美日韩乱码在线| 成人18禁在线播放| 在线a可以看的网站| 欧美午夜高清在线| 可以在线观看的亚洲视频| 成人精品一区二区免费| 美女免费视频网站| 日韩av在线大香蕉| 淫妇啪啪啪对白视频| 国产视频一区二区在线看| 国产男靠女视频免费网站| 色综合亚洲欧美另类图片| 午夜福利免费观看在线| 国产精品久久久久久人妻精品电影| 国产精华一区二区三区| 亚洲成a人片在线一区二区| 日韩欧美国产在线观看| 天天躁夜夜躁狠狠躁躁| 一夜夜www| 成人三级做爰电影| 麻豆成人av在线观看| 丰满人妻一区二区三区视频av | 亚洲av成人av| 欧美3d第一页| 小说图片视频综合网站| 91成年电影在线观看| 国产区一区二久久| 国产精品久久久久久久电影 | 一进一出好大好爽视频| 久久九九热精品免费| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 嫩草影院精品99| 少妇的丰满在线观看| 欧美不卡视频在线免费观看 | 波多野结衣巨乳人妻| 国产亚洲精品第一综合不卡| 久久精品91无色码中文字幕| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 最近在线观看免费完整版| 99国产精品99久久久久| 男女下面进入的视频免费午夜| 久久久久国内视频| 亚洲成av人片在线播放无| 巨乳人妻的诱惑在线观看| 免费搜索国产男女视频| 丁香六月欧美| 午夜a级毛片| 国产成人影院久久av| 夜夜躁狠狠躁天天躁| 此物有八面人人有两片| 中国美女看黄片| 国内精品久久久久精免费| 亚洲成人久久性| 国产91精品成人一区二区三区| 色播亚洲综合网| 可以免费在线观看a视频的电影网站| 手机成人av网站| 日韩 欧美 亚洲 中文字幕| 欧美一区二区国产精品久久精品 | 国产69精品久久久久777片 | 香蕉av资源在线| 精品欧美国产一区二区三| 中文字幕av在线有码专区| 精华霜和精华液先用哪个| 在线观看午夜福利视频| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 欧美人与性动交α欧美精品济南到| 99久久99久久久精品蜜桃| 色噜噜av男人的天堂激情| 久久久久国内视频| 好看av亚洲va欧美ⅴa在| 男女下面进入的视频免费午夜| 久久九九热精品免费| 午夜日韩欧美国产| 一区福利在线观看| 给我免费播放毛片高清在线观看| 欧美性长视频在线观看| av有码第一页| 国产高清videossex| 我的老师免费观看完整版| av在线天堂中文字幕| 欧美另类亚洲清纯唯美| 老汉色∧v一级毛片| 国内毛片毛片毛片毛片毛片| 50天的宝宝边吃奶边哭怎么回事| 午夜福利欧美成人| 亚洲成人久久爱视频| 少妇熟女aⅴ在线视频| 丰满人妻一区二区三区视频av | 久久中文字幕一级| √禁漫天堂资源中文www| 亚洲成人久久性| 欧美日韩乱码在线| 成人欧美大片| 国内揄拍国产精品人妻在线| x7x7x7水蜜桃| 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 欧美另类亚洲清纯唯美| 午夜福利在线在线| 黄色a级毛片大全视频| 午夜福利在线在线| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲美女久久久| 男女做爰动态图高潮gif福利片| 亚洲一卡2卡3卡4卡5卡精品中文| 三级男女做爰猛烈吃奶摸视频| 久久久久久久午夜电影| 中文在线观看免费www的网站 | 久久伊人香网站| 成人av一区二区三区在线看| 亚洲国产欧美一区二区综合| 中文字幕高清在线视频| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 国产亚洲欧美在线一区二区| 高清在线国产一区| 欧美又色又爽又黄视频| 免费搜索国产男女视频| 看黄色毛片网站| 国产又色又爽无遮挡免费看| 毛片女人毛片| 欧美精品啪啪一区二区三区| 日韩三级视频一区二区三区| 午夜成年电影在线免费观看| 午夜免费成人在线视频| 日韩高清综合在线| 欧美日韩国产亚洲二区| 亚洲熟妇中文字幕五十中出| 美女扒开内裤让男人捅视频| 婷婷精品国产亚洲av在线| 久久亚洲真实| 99久久99久久久精品蜜桃| 久久婷婷人人爽人人干人人爱| 一级片免费观看大全| 欧美日韩一级在线毛片| 国产伦人伦偷精品视频| 最新在线观看一区二区三区| 欧美乱码精品一区二区三区| 人人妻,人人澡人人爽秒播| 国产精品一区二区三区四区久久| 一本久久中文字幕| 久久人妻av系列| 久久久久久免费高清国产稀缺| 欧美黑人欧美精品刺激| 一级作爱视频免费观看| 熟妇人妻久久中文字幕3abv| 欧美在线黄色| 我要搜黄色片| aaaaa片日本免费| 在线十欧美十亚洲十日本专区| 欧美黑人精品巨大| 午夜福利在线观看吧| 国产高清videossex| 日韩中文字幕欧美一区二区| 亚洲一区高清亚洲精品| 欧美乱色亚洲激情| 啦啦啦免费观看视频1| 久久精品国产亚洲av高清一级| 免费在线观看黄色视频的| 淫妇啪啪啪对白视频| 999久久久精品免费观看国产| 久久久久久久久久黄片| 99riav亚洲国产免费| 观看免费一级毛片| 亚洲熟女毛片儿| 亚洲成人国产一区在线观看| 99riav亚洲国产免费| 亚洲av成人精品一区久久| 人妻夜夜爽99麻豆av| 国产黄a三级三级三级人| 亚洲人成电影免费在线| 日日爽夜夜爽网站| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 色综合欧美亚洲国产小说| 黄色成人免费大全| 19禁男女啪啪无遮挡网站| 亚洲人成伊人成综合网2020| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美 | 夜夜看夜夜爽夜夜摸| 国产高清有码在线观看视频 | 中文字幕人妻丝袜一区二区| 国产精品永久免费网站| 日韩 欧美 亚洲 中文字幕| 免费在线观看成人毛片| 亚洲 国产 在线| 男人的好看免费观看在线视频 | tocl精华| 亚洲国产高清在线一区二区三| 99国产精品一区二区蜜桃av| 日本熟妇午夜| 国内毛片毛片毛片毛片毛片| 欧美日本视频| 欧美性猛交黑人性爽| 精品久久久久久久久久免费视频| 国产私拍福利视频在线观看| 国产片内射在线| 国产在线观看jvid| 国产成人欧美在线观看| 天天躁夜夜躁狠狠躁躁| 精品熟女少妇八av免费久了| 每晚都被弄得嗷嗷叫到高潮| 两性夫妻黄色片| av免费在线观看网站| 制服丝袜大香蕉在线| 久久久久久久久中文| 亚洲国产欧洲综合997久久,| 久久热在线av| 午夜免费观看网址| 在线播放国产精品三级| 国产免费av片在线观看野外av| 女人被狂操c到高潮| 亚洲av成人不卡在线观看播放网| 午夜老司机福利片| 好男人电影高清在线观看| 午夜精品在线福利| 一二三四在线观看免费中文在| 精品免费久久久久久久清纯| 日韩国内少妇激情av| 亚洲性夜色夜夜综合| 丁香六月欧美| 久久久久久久午夜电影| 成人精品一区二区免费| 91成年电影在线观看| 啦啦啦韩国在线观看视频| 久久精品国产清高在天天线| 久久午夜综合久久蜜桃| 他把我摸到了高潮在线观看| 在线播放国产精品三级| 香蕉av资源在线| 法律面前人人平等表现在哪些方面| 国产成人系列免费观看| 婷婷丁香在线五月| av在线播放免费不卡| 可以免费在线观看a视频的电影网站| 婷婷精品国产亚洲av在线| av福利片在线| 欧美+亚洲+日韩+国产| av福利片在线| 好男人电影高清在线观看| 久9热在线精品视频| 欧美成人午夜精品| 禁无遮挡网站| 我的老师免费观看完整版| 青草久久国产| 身体一侧抽搐| 成人av在线播放网站| 动漫黄色视频在线观看| 成人手机av| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| 三级男女做爰猛烈吃奶摸视频| 亚洲黑人精品在线| 欧美精品亚洲一区二区| 国产精品久久视频播放| 国产亚洲av高清不卡| 久久亚洲真实| 一边摸一边做爽爽视频免费| 亚洲欧美日韩高清专用| 久久久久久免费高清国产稀缺| 久久久久九九精品影院| 精品国产亚洲在线| 国产熟女午夜一区二区三区| 999久久久国产精品视频| 99精品久久久久人妻精品| 999久久久精品免费观看国产| 男插女下体视频免费在线播放| 一边摸一边做爽爽视频免费| 欧美乱妇无乱码| 亚洲熟女毛片儿| 日韩大码丰满熟妇| 一区二区三区激情视频| 99在线视频只有这里精品首页| 不卡一级毛片| 午夜精品在线福利| 99久久久亚洲精品蜜臀av| 欧美成狂野欧美在线观看| 午夜两性在线视频| 1024手机看黄色片| 亚洲全国av大片| 日本一区二区免费在线视频| 午夜福利在线在线| 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| 欧美乱妇无乱码| 人妻夜夜爽99麻豆av| 亚洲一卡2卡3卡4卡5卡精品中文| 97碰自拍视频| 成年免费大片在线观看| 最好的美女福利视频网| e午夜精品久久久久久久| 人人妻,人人澡人人爽秒播| 在线看三级毛片| 日韩精品青青久久久久久| 99久久国产精品久久久| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 成人三级黄色视频| 国产激情久久老熟女| 成人国产一区最新在线观看| 俺也久久电影网| 成年人黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看 | 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 叶爱在线成人免费视频播放| 人妻丰满熟妇av一区二区三区| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| 亚洲中文字幕日韩| 精品免费久久久久久久清纯| 一级毛片精品| 国产91精品成人一区二区三区| 国产单亲对白刺激| 一区二区三区高清视频在线| 中亚洲国语对白在线视频| 亚洲专区字幕在线| 久久精品综合一区二区三区| 黑人欧美特级aaaaaa片| avwww免费| 亚洲人成伊人成综合网2020| 欧美日韩黄片免| 青草久久国产| 久久精品人妻少妇| 成年女人毛片免费观看观看9| 男女下面进入的视频免费午夜| 国内久久婷婷六月综合欲色啪| 免费电影在线观看免费观看| 脱女人内裤的视频| 99国产综合亚洲精品| 一本大道久久a久久精品| 一夜夜www| 757午夜福利合集在线观看| 免费看美女性在线毛片视频| 一级毛片高清免费大全| 亚洲av成人不卡在线观看播放网| 在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 我的老师免费观看完整版| 性色av乱码一区二区三区2| 亚洲国产欧美网| 久久香蕉国产精品| 精品国产乱子伦一区二区三区| 亚洲精品久久成人aⅴ小说| 欧美黄色淫秽网站| 久久久久亚洲av毛片大全| 99精品久久久久人妻精品| 亚洲电影在线观看av| 成在线人永久免费视频| 日日夜夜操网爽| 88av欧美| 午夜成年电影在线免费观看| 国产区一区二久久| 亚洲 欧美 日韩 在线 免费| 波多野结衣高清作品| 老司机午夜福利在线观看视频| 国产精品98久久久久久宅男小说| 一区二区三区激情视频| 丝袜美腿诱惑在线| 99国产精品99久久久久| 亚洲精品美女久久久久99蜜臀| 色尼玛亚洲综合影院| 啦啦啦免费观看视频1| 免费在线观看亚洲国产| 亚洲全国av大片| 久久久国产成人免费| 国产午夜福利久久久久久| 欧美在线一区亚洲| 成人永久免费在线观看视频| 午夜福利免费观看在线| 丁香六月欧美| 十八禁网站免费在线| 精品第一国产精品| 国产av一区二区精品久久| 可以免费在线观看a视频的电影网站| 亚洲精品在线观看二区| 久99久视频精品免费| 日韩中文字幕欧美一区二区| 中亚洲国语对白在线视频| 国产免费男女视频| 我要搜黄色片| 国产亚洲精品av在线| 国产一区二区在线av高清观看| 中文字幕最新亚洲高清| 中文字幕人成人乱码亚洲影| 热99re8久久精品国产| 成人av一区二区三区在线看| 午夜精品在线福利| 亚洲人成网站在线播放欧美日韩| 99久久无色码亚洲精品果冻| 亚洲av片天天在线观看| 亚洲av日韩精品久久久久久密| 精品第一国产精品| 亚洲自拍偷在线| 一级作爱视频免费观看| 国产精品亚洲一级av第二区| 国产精品久久视频播放| 国产精品一区二区三区四区免费观看 | 午夜福利在线在线| 狂野欧美白嫩少妇大欣赏| ponron亚洲| av福利片在线| 老熟妇仑乱视频hdxx| 在线观看免费午夜福利视频| 久久久久久免费高清国产稀缺| 久久中文字幕人妻熟女| 久久久久久久久免费视频了| 草草在线视频免费看| 成人国语在线视频| 波多野结衣高清无吗| 国产熟女午夜一区二区三区| 一区二区三区激情视频| 久久精品91无色码中文字幕| 日本熟妇午夜| 舔av片在线| 亚洲精品av麻豆狂野| 亚洲精品色激情综合| 久久伊人香网站| 久久久久久久午夜电影| 高潮久久久久久久久久久不卡| 午夜福利欧美成人| 91麻豆av在线| 久久国产乱子伦精品免费另类| 一进一出抽搐动态| 999精品在线视频| 两性夫妻黄色片| 一边摸一边做爽爽视频免费| 免费搜索国产男女视频| www日本在线高清视频| 国产精品av久久久久免费| 国产成人av激情在线播放| 他把我摸到了高潮在线观看| 老司机福利观看| 日韩 欧美 亚洲 中文字幕| 国产97色在线日韩免费| 97人妻精品一区二区三区麻豆| 俺也久久电影网| 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| 国产99白浆流出| 欧美日韩一级在线毛片| 久久婷婷成人综合色麻豆| 在线观看免费午夜福利视频| 波多野结衣高清无吗| 国产日本99.免费观看| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 日本在线视频免费播放| 婷婷六月久久综合丁香| 久久午夜综合久久蜜桃| 久久久久久久久久黄片| 国产探花在线观看一区二区| 国产一区二区激情短视频| 99久久99久久久精品蜜桃| 制服丝袜大香蕉在线| 欧美日韩精品网址| 18禁裸乳无遮挡免费网站照片| 一级黄色大片毛片| 成年人黄色毛片网站| 12—13女人毛片做爰片一| av福利片在线| 大型黄色视频在线免费观看| 久99久视频精品免费| av天堂在线播放| 欧美中文日本在线观看视频| avwww免费| 久久人妻av系列| 国产午夜精品久久久久久| 国产亚洲精品av在线| 国产伦人伦偷精品视频| 成人高潮视频无遮挡免费网站| 亚洲五月天丁香| 全区人妻精品视频| 久9热在线精品视频| 伊人久久大香线蕉亚洲五| 亚洲一区高清亚洲精品| 在线观看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 美女黄网站色视频| 窝窝影院91人妻| 久久久久久大精品| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看免费日韩欧美大片| 亚洲中文字幕一区二区三区有码在线看 | 国产精品自产拍在线观看55亚洲| 亚洲va日本ⅴa欧美va伊人久久| 亚洲五月天丁香| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 日韩欧美 国产精品| 久久久精品国产亚洲av高清涩受| 在线视频色国产色| 亚洲成a人片在线一区二区| 日韩 欧美 亚洲 中文字幕| 中文在线观看免费www的网站 | 人妻丰满熟妇av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 中文资源天堂在线| 99国产精品一区二区蜜桃av| 日韩中文字幕欧美一区二区| 狂野欧美白嫩少妇大欣赏| 欧美av亚洲av综合av国产av| 国内久久婷婷六月综合欲色啪| 天堂动漫精品| 亚洲av片天天在线观看| 岛国在线免费视频观看| 国产精品久久久久久久电影 | 99国产极品粉嫩在线观看| 欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| 久久国产精品人妻蜜桃| 国产精品日韩av在线免费观看| 欧美绝顶高潮抽搐喷水| 最近最新免费中文字幕在线| 亚洲一区二区三区色噜噜| 国产成人一区二区三区免费视频网站| 少妇人妻一区二区三区视频| 亚洲第一电影网av| 国产精品一区二区三区四区久久| 亚洲国产日韩欧美精品在线观看 | 日本黄色视频三级网站网址| 狂野欧美白嫩少妇大欣赏| 美女午夜性视频免费|