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

    離腹寡毛果實(shí)蠅性別決定基因transformer 2及其在昆蟲中的進(jìn)化

    2015-12-09 09:15:10劉桂清萬方浩
    環(huán)境昆蟲學(xué)報(bào) 2015年4期
    關(guān)鍵詞:毛果實(shí)蠅同義

    劉桂清,萬方浩

    (1.中國農(nóng)業(yè)科學(xué)院植物保護(hù)研究所,植物病蟲害生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100193;2.廣東省昆蟲研究所,廣東省野生動(dòng)物保護(hù)和利用公共實(shí)驗(yàn)室,廣東省農(nóng)業(yè)害蟲綜合治理重點(diǎn)實(shí)驗(yàn)室,廣州 510260;3.青島農(nóng)業(yè)大學(xué)農(nóng)學(xué)與植物保護(hù)學(xué)院,山東青島 266109)

    離腹寡毛實(shí)蠅屬Bactrocera Macquart、小條實(shí)蠅屬Ceratitis Macleay 和按實(shí)蠅屬Anastrepha Schiner是實(shí)蠅科Tephritidae 中寄主范圍最廣,最具經(jīng)濟(jì)重要性的三大類群。世界已知的離腹寡毛果實(shí)蠅有561種,有記載的為害經(jīng)濟(jì)植物的種類有53種,其中具有經(jīng)濟(jì)重要性的種類36種,除橄欖實(shí)蠅B.oleae (Rossi)原產(chǎn)于南歐,現(xiàn)擴(kuò)至非洲、亞洲印度北部和巴基斯坦西北部外,多數(shù)種主產(chǎn)舊大陸的亞洲熱帶、亞熱帶和暖溫帶,以及澳大利亞和南太平洋地區(qū),少數(shù)種類分布于非洲南部和夏威夷群島。該屬包含許多世界性果蔬植物的檢疫害蟲,如桔小實(shí)蠅B.dorsalis (Hendel)、橄欖實(shí)蠅、瓜實(shí)蠅B.cucurbitae (Coquillett)、番石榴實(shí)蠅B.correcta (Bezzi)、桔大實(shí)蠅B.minax (Enderlein)、辣椒實(shí)蠅B.latifrons (Hendel)、昆士蘭果實(shí)蠅B.tryoni (Froggatt)、扎氏果實(shí)蠅Bactrocera jarvisi(Tryon)、桃實(shí)蠅B.zonata (Saunders)、南瓜實(shí)蠅B.tau (Walker)等,可危害水果類寄主植物有2l 科37 屬,蔬菜類寄主植物有4 科16 屬(梁廣勤等,1996)。

    鑒于離腹寡毛實(shí)蠅屬危害范圍廣及對(duì)國內(nèi)外貿(mào)易造成的嚴(yán)重?fù)p失,世界上許多國家和地區(qū)都將離腹寡毛果實(shí)蠅列為重要的危險(xiǎn)性檢疫害蟲,我國農(nóng)業(yè)部2007年5月28 日頒布實(shí)施的《中華人民共和國進(jìn)境植物檢疫性有害生物名錄》中,離腹寡毛實(shí)蠅屬所有物種都被列為檢疫性對(duì)象。除加強(qiáng)檢疫外,昆蟲不育技術(shù) (Sterile insect technique,SIT)是根除重大實(shí)蠅害蟲或預(yù)防其在新地區(qū)定殖的有效方法 (Klassen and Curtis,2005),該技術(shù)的核心是通過大規(guī)模釋放不育雄蟲與野生雌蟲發(fā)生不育交配從而達(dá)到控制該地區(qū)害蟲種群的目的(Knipling,1955)。在輻照處理蛹之前將靶標(biāo)雌、雄蟲分離是不育雄蟲釋放的前提,不育雄蟲釋放不僅可以提高SIT 的防治效果(Rendón et al.,2004),而且是病媒昆蟲如蚊蟲昆蟲不育技術(shù)項(xiàng)目實(shí)施的先決條件 (Catteruccia et al.,2005)。

    通過分析性別決定通路關(guān)鍵基因的進(jìn)化,有助于闡明昆蟲性別分化的分子機(jī)制,而且為發(fā)展雌、雄蟲性別分離策略,提高昆蟲不育昆蟲技術(shù)的應(yīng)用效果和開發(fā)新的害蟲種群遺傳控制方法奠定基礎(chǔ)。目前,黑腹果蠅Drosophila melanogaster Meigen 性別決定分子機(jī)制研究相對(duì)成熟,為離腹寡毛果實(shí)蠅性別決定機(jī)制的比較研究提供了依據(jù)。在果蠅中,遺傳級(jí)聯(lián) (sex lethal >transformer +transformer 2 >double sex)上游基因通過性別特異剪接方式產(chǎn)生不同的產(chǎn)物控制級(jí)聯(lián)下游基因的剪接,從而控制雌雄分化 (Bell et al.,1991;Sánchez,2008)。Tra-2 基因是性別決定通路關(guān)鍵基因,tranformer (tra)需要與trasformer 2 (tra-2)相互結(jié)合形成Tra/Tra-2 復(fù)合體,調(diào)控下游基因double sex (dsx)pre-mRNA 以雌性特異方式剪接產(chǎn)生雌性特異Dsx 蛋白。胚胎顯微注射tra-2 基因雙鏈RNA 會(huì)導(dǎo)致XX 染色體胚胎發(fā)育為雄蟲表型(Salvemini et al.,2009;Schetelig et al.,2012),這種性別反轉(zhuǎn)有利于只釋放不育雄蟲的害蟲防治項(xiàng)目。目前,有關(guān)離腹寡毛果實(shí)蠅tra-2 基因生物信息學(xué)方面的系統(tǒng)研究尚未見報(bào)道,故本研究首先克隆得到瓜實(shí)蠅tra-2 基因,對(duì)7種離腹寡毛果實(shí)蠅tra-2 基因結(jié)構(gòu)進(jìn)行比較研究,并對(duì)18種昆蟲tra-2 基因進(jìn)行系統(tǒng)發(fā)育和分子進(jìn)化分析。

    1 材料與方法

    1.1 瓜實(shí)蠅基因組DNA 和總RNA 的提取和文庫構(gòu)建

    瓜實(shí)蠅為廣東出入境檢驗(yàn)檢疫局檢驗(yàn)檢疫技術(shù)中心室內(nèi)飼養(yǎng)蟲源。成蟲以人工飼料飼養(yǎng),幼蟲以南瓜飼養(yǎng)。飼養(yǎng)條件為27℃±1℃、RH 75%±5%,光周期為14 h∶10 h (光照∶黑暗)。選取初羽化瓜實(shí)蠅雌、雄蟲提取基因組DNA 和總RNA。

    參照Bender et al.(1983)的方法提取瓜實(shí)蠅成蟲基因組DNA。采用TriZol Reagent (Invitrogen)分別提取單頭瓜實(shí)蠅雌、雄成蟲總RNA,1.2%瓊脂糖凝膠電泳檢測(cè) RNA 樣本的完整性,NanophotometerTMP-Class (IMPLEN,德國)檢測(cè)RNA 的純度和濃度,每個(gè)樣本取2 μg 總量的合格RNA,按 照 TransScript TM One-Step gDNA Removal and cDNA Synthesis SuperMix (北京全式金生物技術(shù)有限公司)第一鏈cDNA 合成試劑盒說明書構(gòu)建cDNA 文庫;按照BD SMARTer RACE cDNA Amplification Kit (BD Biosciences)試劑盒說明書分別構(gòu)建雌、雄蟲5'-RACE cDNA 和3'-RACE cDNA 文庫,具體步驟及注意事項(xiàng)參照試劑盒說明書。

    1.2 引物設(shè)計(jì)

    根據(jù)地中海實(shí)蠅Ceratitis capitata、橄欖實(shí)蠅、加勒比按實(shí)蠅Anastrepha suspensa 同源tra-2 基因的氨基酸序列比對(duì)結(jié)果,設(shè)計(jì)簡并引物F30/R530 擴(kuò)增瓜實(shí)蠅tra-2 的中間片段,在獲得的中間片段序列基礎(chǔ)上,分別設(shè)計(jì)5' RACE 和3' RACE 特異引物,擴(kuò)增tra-2 的cDNA 全長。參照番石榴實(shí)蠅同源tra-2 基因的分子結(jié)構(gòu),設(shè)計(jì)3 對(duì)外顯子特異引物F13/R337,F(xiàn)327/R711,F(xiàn)686/R1140 用于擴(kuò)增瓜實(shí)蠅tra-2 的內(nèi)含子序列。本研究所用引物均使用Primer Premier 5.0 設(shè)計(jì),引物序列詳細(xì)信息見表1,所有引物均由生工生物工程(上海)股份有限公司合成。

    表1 瓜實(shí)蠅tra-2 基因克隆引物Table 1 Primers for cloning Bactrocera cucurbitae homolog tra-2

    1.3 瓜實(shí)蠅同源tra-2 基因克隆

    以單頭雌、雄蟲第一鏈cDNA 為模板,采用簡并引物F30/R530 擴(kuò)增tra-2 的中間片段,PCR 反應(yīng)體系為:滅菌ddH2O 水19.5 μL,10×PCR reaction buffer 2.5 μL,dNTP Mix (10mM)0.5 μL,上、下游引物(10 μmol/L)各0.5 μL,TranStart Taq DNA 聚合酶 (2.5U)0.5 μL,cDNA 模板1 μL。反應(yīng)程序?yàn)?5℃變性3 min;94℃,30 s,55℃,30 s,72℃,1 min,共35 個(gè)循環(huán);72℃延伸7 min。PCR 產(chǎn)物經(jīng)1 %瓊脂糖凝膠電泳檢測(cè),割取目的大小條帶并用 EasyPure Quick Gel Extraction Kit 瓊脂糖凝膠DNA 回收試劑盒(北京全式金生物技術(shù)有限公司)說明書回收PCR 產(chǎn)物,參照pEASY-T1 Simple Cloning Kit 克隆試劑盒(北京全式金生物技術(shù)有限公司)說明書將回收產(chǎn)物連接到pEASY-T1 克隆載體,重組質(zhì)粒轉(zhuǎn)化Trans1-T1 感受態(tài)細(xì)胞,然后涂布在含有Ampicillin/X-gal/IPTG 的LB 平板上,37℃倒置培養(yǎng)12 h,藍(lán)白斑篩選,挑取5 個(gè)白色菌落進(jìn)行菌液PCR,并將菌落PCR 陽性克隆株挑入LB 液體培養(yǎng)基(含 Amp) 中過夜培養(yǎng),按照 AxyPrepTMPlasmid Miniprep Kit 質(zhì)粒提取試劑盒(Axygen)說明書提取質(zhì)粒,送至生工生物工程(上海)股份有限公司進(jìn)行測(cè)序。

    為擴(kuò)增tra-2 基因的cDNA 全長,以5' RACE Ready cDNA 為模板,5' RACE out/in 為引物擴(kuò)增tra-2 5' UTR;以3' RACE Ready cDNA 為模板,3'RACE out/in 擴(kuò)增tra-2 3' UTR,RACE-PCR 反應(yīng)程序如下:95℃預(yù)變性3 min;94℃,30 s,72℃,2 min 擴(kuò)增5 個(gè)循環(huán);94℃,30 s;70℃,30 s;72℃,1 min 擴(kuò)增5 個(gè)循環(huán);94℃,30 s;68℃,30 s;72℃,1 min 擴(kuò)增35 個(gè)循環(huán);72℃延伸10 min。最后,以基因組DNA 為模板,分別以引物F13/R337、F327/R711、F686/R1140 擴(kuò)增該基因的內(nèi)含子,PCR 反應(yīng)程序?yàn)?95℃退火3 min,94℃,30 s,60℃,45 s,72℃,1 min 擴(kuò)增35 個(gè)循環(huán),72℃延伸7 min。所有PCR 產(chǎn)物經(jīng)1 %瓊脂糖凝膠電泳檢測(cè),割膠回收,并將回收產(chǎn)物連接至pEASY-T1 克隆載體,重組質(zhì)粒轉(zhuǎn)化Trans1-T1 感受態(tài)細(xì)胞,篩選陽性克隆株,送測(cè)序,具體步驟同上。

    1.4 序列分析

    利用Jemboss 軟件進(jìn)行核苷酸序列翻譯,并預(yù)測(cè)蛋白質(zhì)的分子量、等電點(diǎn)和跨膜結(jié)構(gòu)域。NCBI網(wǎng)站Blast 工具對(duì)序列進(jìn)行相似性比對(duì)分析,利用Spidey 分析基因外顯子和內(nèi)含子的邊界(http://www.ncbi.nlm.nih.gov/spidey/)。利 用 Bioedit(Hall,1999)和Clustal W (Thompson et al.,1994)軟件進(jìn)行氨基酸序列比對(duì)分析,參照地中海實(shí)蠅Tra-2 蛋白結(jié)構(gòu)特征確定Tra-2 蛋白的功能結(jié)構(gòu)域,如N 端和C 端精氨酸/絲氨酸富集區(qū)(Arginine/Serine rich region,RS)、RNA 識(shí)別基序 (RNA recognition motif,RRM) 等 (Salvemini et al.,2009)。

    1.5 分子進(jìn)化分析

    本文中引用了18種昆蟲的Tra-2 序列進(jìn)行分析,物種及登錄號(hào)分別為:雙翅目:瓜實(shí)蠅(GenBank:KR056217)、桔小實(shí)蠅 (GenBank:KP342061)、番石榴實(shí)蠅(GenBank:KM658207)、橄欖實(shí)蠅(GenBank:AJ547623)、昆士蘭果實(shí)蠅(GenBank:KJ443723)、扎式果實(shí)蠅 B.jarvisi(Tryon)(GenBank:KJ443722)、西印度按實(shí)蠅Anastrepha obliqua (Macquart) (GenBank:FN658607)、加勒比按實(shí)蠅A.suspensa (Loew)(GenBank:JN597290)、地中海實(shí)蠅 Ceratitis capitata (Wiedemann)(GenBank:EU999754)、家蠅Musca domestica L.(GenBank:AY847518)、銅綠蠅 Lucilia cuprinia (Wiedemann) (GenBank:FJ461620)、黑腹果蠅 Drosophila melanogaster Meigen (GenBank:M23633)、黑果蠅D.virilis Sturtevant (GenBank:XM_ 002049663)、擬暗果蠅D.pseudoobscura Frolova & Astaurov (GenBank:XM_ 001360568);鱗翅目:家蠶Bombyx mori L.(GenBank:NM_ 001126233);膜翅目:意大利蜜蜂Apis mellifera L.(GenBank:XM_ 001121070)、麗蠅蛹集金小蜂 Nasonia vitripennis Walker(GenBank:XP_ 001601106)。

    利用MEGA 5 軟件進(jìn)行分子進(jìn)化分析,采用Kimura 2 參數(shù)模型(Kimura,1980)估算核苷酸序列分異 (pNT),并計(jì)算轉(zhuǎn)換與顛換的比值(transition/transversion ratio,R)。采用差異位點(diǎn)占所有比較位點(diǎn)的百分?jǐn)?shù)(p-distance)估算Tra-2蛋白及其RS 和RRM 保守結(jié)構(gòu)域的氨基酸序列分異(pAA)(Nei and Kumar,2000)。采用Nei-Gojobori 修正算法 (Zhang et al.,1998),假設(shè)transition/transversion ratio 為1,計(jì)算核苷酸同義突變率(pS)和非同義突變率(pN)。采用競(jìng)爭刪除選項(xiàng)(Compete deletion option)計(jì)算進(jìn)化距離,鄰接法(Neighbor-joining method)(Tamura et al.,2011)構(gòu)建系統(tǒng)發(fā)育樹,重復(fù)構(gòu)建1000 次。同時(shí)利用DnaSP v.5 軟件(Librado and Rozas,2009)分析核苷酸序列分異,所有突變位點(diǎn)經(jīng)Tajima D檢驗(yàn)。并通過計(jì)算Effective Number of Codons(ENC)、Codon BiasIndex (CBI)和Scaled Chisquare (Schi2)分析該基因的密碼子偏好性。

    2 結(jié)果與分析

    2.1 瓜實(shí)蠅tra-2 基因結(jié)構(gòu)

    以第一鏈cDNA 為模板,以F30/R530 為引物,RT-PCR 擴(kuò)增得到片段長度為501 bp 的產(chǎn)物,BlastX 分析發(fā)現(xiàn)推斷的氨基酸序列與橄欖實(shí)蠅、暗色實(shí)蠅A.serpentina (Wiedemann)(GenBank:FN658613)Tra 2 蛋白序列的相似度為97%,推斷得到的501 bp 的片段為瓜實(shí)蠅tra-2 基因序列。通過RACE 技術(shù)擴(kuò)增得到tra-2 cDNA 全長1467 bp,開放讀碼框753 bp,共編碼250 個(gè)氨基酸殘基,分子量為295.59 kD,等電點(diǎn)PI 為6.87,有2 個(gè)跨膜結(jié)構(gòu)域。瓜實(shí)蠅tra-2 與桔大實(shí)蠅同源tra-2 基因的序列一致性最高,達(dá)到90%,其次是與其他離腹寡毛果實(shí)蠅,與地中海實(shí)蠅tra-2 基因的序列一致性為83.6%,與黑腹果蠅的最低,僅為51.6%(表2)。7種離腹寡毛果實(shí)蠅tra-2 基因結(jié)構(gòu)類似,均包含8 個(gè)外顯子和7 個(gè)內(nèi)含子,外顯子和內(nèi)含子大小及其翻譯起始位點(diǎn)見圖1,翻譯起始密碼子均位于第1 個(gè)外顯子末端,與地中海實(shí)蠅tra-2 基因結(jié)構(gòu)相似(Salvemini et al.,2009),與黑腹果蠅tra-2 (翻譯起始密碼子位于第2 個(gè)外顯子)不同。

    圖1 瓜實(shí)蠅tra-2 基因結(jié)構(gòu)及與同屬其他種tra-2 基因結(jié)構(gòu)比較Fig.1 Comparison of genomic organization of Bactrocera cucurbitae tra-2 with that from other Bactrocera species

    2.2 瓜實(shí)蠅和其他離腹寡毛果實(shí)蠅Tra-2 蛋白

    以該基因推導(dǎo)的Tra-2 氨基酸序列具有SR 蛋白家族的典型特征,如RNA 結(jié)合基序 (RNA Recognition Motif,RRM)及其側(cè)翼的2 個(gè)精氨酸/絲氨酸富集結(jié)構(gòu)域(Arginine/Serine rich domain,RS),暗示Tra-2 蛋白可以與其他蛋白結(jié)合發(fā)揮其功能。7種離腹寡毛果實(shí)蠅之間的序列一致性最高,其中扎式果實(shí)蠅與昆士蘭果實(shí)蠅的氨基酸序列一致性最高,達(dá)到99.6%,橄欖實(shí)蠅和瓜實(shí)蠅的最低,為94%,番石榴實(shí)蠅與桔小實(shí)蠅的核苷酸序列一致性最高,達(dá)到97.6%,番石榴實(shí)蠅與瓜實(shí)蠅的最低,為87.4%。7種離腹寡毛果實(shí)蠅與地中海實(shí)蠅的序列一致性較高,與黑腹果蠅的最低,只有45.7% (表2)。在氨基酸殘基數(shù)量上,瓜實(shí)蠅和桔大實(shí)蠅的Tra-2 蛋白均有250 個(gè)氨基酸殘基數(shù),比桔小實(shí)蠅、番石榴實(shí)蠅和扎實(shí)果實(shí)蠅的Tra-2 蛋白少1 個(gè),比昆士蘭果實(shí)蠅的Tra-2 蛋白少5 個(gè),但氨基酸差異主要發(fā)生在RS1 和RS2結(jié)構(gòu)域,所有物種的RNA 識(shí)別基序RRM 結(jié)構(gòu)域均有相同的氨基酸殘基數(shù)(圖2)。

    2.3 Tra-2 系統(tǒng)發(fā)育和分子進(jìn)化分析

    以家蠶、意大利蜜蜂、麗蠅蛹集金小蜂的Tra-2 蛋白作為系統(tǒng)進(jìn)化樹的外圍支,重構(gòu)雙翅目Tra-2 蛋白的進(jìn)化關(guān)系(圖3),結(jié)果表明,家蠶、意大利蜜蜂、麗蠅蛹集金小蜂Tra-2 蛋白代表1 個(gè)基礎(chǔ)分支,雙翅目昆蟲Tra-2 代表另1 個(gè)分支。在雙翅目昆蟲中,果蠅單獨(dú)聚為一支,蠅科的家蠅和麗蠅科的銅綠蠅以及實(shí)蠅科果實(shí)蠅聚為另一支,在實(shí)蠅科果實(shí)蠅中,按實(shí)蠅屬與小條實(shí)蠅屬先聚在一起,然后與離腹寡毛實(shí)蠅屬聚在一起。離腹寡毛實(shí)蠅屬Tra-2 蛋白與小條實(shí)蠅屬和按實(shí)蠅屬的最為接近,其次是麗蠅科和蠅科的Tra-2 蛋白,這些物種Tra-2 蛋白表現(xiàn)為單一起源,而果蠅Tra-2 蛋白表現(xiàn)出實(shí)質(zhì)性的分化。

    表2 Tra-2 氨基酸和cDNA 序列一致性Table 2 Tra-2 amino acids and cDNA sequence percentage identity

    圖2 雙翅目昆蟲Tra-2 氨基酸序列比較Fig.2 Comparison of the putative Tra-2 polypeptides from Diptera species

    圖3 鄰接法構(gòu)建18種昆蟲Tra-2 蛋白系統(tǒng)發(fā)育樹Fig.3 Phylogenetic tree reconstructed from 18 Tra-2 proteins using the neighbor-joining method

    表3 Tra-2 基因的同義和非同義分異Table 3 Synonymous and nonsynonymous divergence of tra-2 genes

    DnaSP5 分析18種昆蟲tra-2 基因核苷酸序列比對(duì)的1152 個(gè)核苷酸位點(diǎn),發(fā)現(xiàn)425 個(gè)核苷酸位點(diǎn)存在變異,其中372 個(gè)為簡約信息位點(diǎn)(Parsimony information sites),53 個(gè)為單突變位點(diǎn)(Singleton sites)。Tra-2 基因序列分異結(jié)果如表3所示,同義分異Ks 均大于非同義分異Ka,其中同義替換位點(diǎn)數(shù)(SS)為133.33-144.50,同義替換核苷酸多樣性均值Pi (s)為0.63599;非同義替換位點(diǎn)數(shù)(NSS)為395.50-406.67,非同義替換核苷酸多樣性均值Pi (a)為0.25439,表明同義變異引起的tra-2 基因核苷酸變異高于非同義變異,核苷酸多樣性在不同的核苷酸區(qū)段也存在差異,如圖4A 所示,tra-2 基因RS1 和RS2 結(jié)構(gòu)域的核苷酸多樣性高于RRM 結(jié)構(gòu)域。所有突變位點(diǎn)經(jīng)Tajima D 檢驗(yàn),Tajima's D 為-0.94782,差異不顯著(P >0.10)(圖4B),表明tra-2 基因核苷酸變異符合中性進(jìn)化模型。

    圖4 18種昆蟲tra-2 基因編碼核苷酸多態(tài)性和Tajima 中性檢測(cè)Fig.4 DNA polymorphism of tra-2 coding sequences from 18 insects and Tajima's test of neutrality

    在氨基酸和核苷酸水平分析18種昆蟲tra-2 基因的變異,同時(shí)分析Tra-2 蛋白R(shí)RM、RS1、RS2結(jié)構(gòu)域的氨基酸和核苷酸變異,結(jié)果如表4 所示,tra-2 基因的總體氨基酸和核苷酸變異分別為0.333±0.019 和0.533±0.022,均高于實(shí)蠅科昆蟲tra-2 基因 (氨基酸和核苷酸變異分別為0.073±0.010 和0.148± 0.010),實(shí)蠅科 蟲tra-2蛋白高度保守。Tra-2 蛋白及其RRM、RS1 和RS2 結(jié)構(gòu)域的同義突變水平均顯著高于非同義突變水平(Z-test,P<0.001),RS1 結(jié)構(gòu)域變異高于RS2,兩者均高于RRM 結(jié)構(gòu)域,說明tra-2 基因變異主要來源于同義突變,受凈化選擇壓力,RRM與RS 結(jié)構(gòu)域相比更為保守。密碼子偏好性分析表明Tra-2 RRM 結(jié)構(gòu)域密碼子偏好性低于RS1 和RS2結(jié)構(gòu)域。

    表4 18種昆蟲Tra-2 蛋白及其保守結(jié)構(gòu)域進(jìn)化參數(shù)Table 4 Evolutionary parameters of the 18 Tra-2 proteins and their conserved domains

    3 結(jié)論與討論

    通過同源基因克隆和RACE 技術(shù)克隆得到瓜實(shí)蠅tra-2 基因組DNA 和cDNA 全長,瓜實(shí)蠅tra-2基因與其他6種寡毛果實(shí)蠅和地中海實(shí)蠅的基因結(jié)構(gòu)一致,均包含8 個(gè)外顯子和7 個(gè)內(nèi)含子,而且存在內(nèi)含子剪接抑制ISS 位點(diǎn)和RNA 結(jié)合蛋白R(shí)BP1 結(jié)合位點(diǎn)。瓜實(shí)蠅tra-2 基因不存在性別特異剪接,只有1 個(gè)轉(zhuǎn)錄本,編碼250 個(gè)氨基酸殘基,編碼氨基酸具有SR 蛋白家族特征,如RS 和RRM 結(jié)構(gòu)域。與黑腹果蠅tra-2 存在4 個(gè)轉(zhuǎn)錄本(Amrein et al.,1990;McGuffin et al.,1998)不同,瓜實(shí)蠅tra-2 與地中海實(shí)蠅(Salvemini et al.,2009)、家蠅 (Burghardt et al.,2005)、銅綠蠅(Concha and Scott,2009)、加勒比實(shí)蠅(Schetelig et al.,2012)等同源tra-2 基因均只有1個(gè)轉(zhuǎn)錄本,并未檢測(cè)到生殖系特異的轉(zhuǎn)錄本。序列分析結(jié)果表明瓜實(shí)蠅Tra-2 與其他離腹寡毛果實(shí)蠅和地中海實(shí)蠅的氨基酸和核苷酸的序列一致性均比黑腹果蠅的Tra-2 高,Tra-2 系統(tǒng)發(fā)育結(jié)果也表明7種離腹寡毛果實(shí)蠅與實(shí)蠅科的地中海實(shí)蠅和按實(shí)蠅、蠅科的家蠅和麗蠅科的銅綠蠅的遺傳距離更近,為單一起源,而果蠅Tra-2 表現(xiàn)出明顯的分異。分子進(jìn)化分析結(jié)果表明tra-2 基因核苷酸變異主要來源于同義變異,Tra-2 蛋白受凈化選擇壓力,以維持RNA 結(jié)合蛋白的典型特征,特別是RNA 識(shí)別基序RRM 的結(jié)構(gòu)。

    在地中海實(shí)蠅、橄欖實(shí)蠅、銅綠蠅和加勒比按實(shí)蠅等昆蟲中,Tra-2 通過與Tra 結(jié)合形成Tra/Tra-2復(fù)合體,在性別決定中發(fā)揮雙重功能:一方面通過與下游基因dsx 雌性特異外顯子結(jié)合,調(diào)控dsx 雌性特異剪接;另一方面通過與tra 雄性特異外顯子結(jié)合,抑制這些外顯子被剪接翻譯為成熟的mRNA (Concha and Scott,2009;Salvemini et al.,2009;Schetelig et al.,2012),Tra-2 內(nèi)含子剪接抑制ISS 位點(diǎn)在其中發(fā)揮關(guān)鍵作用 (Ruiz et al.,2007)。由于Tra 缺乏RNA 結(jié)合蛋白結(jié)構(gòu)域,因此,Tra 發(fā)揮作用必須與RNA 結(jié)合蛋白如Tra-2 相互作用(Black,2003)。tra-2 基因核苷酸序列同義和非同義替換率結(jié)果表明該基因受凈化選擇壓力,通過維持Tra-2 蛋白結(jié)構(gòu)的高度保守發(fā)揮其在昆蟲性別決定中的作用。在果蠅中,Tra 已發(fā)生高水平的中性進(jìn)化(Kulathinal et al.,2003),而Tra-2 精氨酸/絲氨酸富集區(qū)RS 是Tra 和Tra-2相互結(jié)合的位置,因此,tra-2 基因RS 結(jié)構(gòu)域核苷酸變異很可能是為了適應(yīng)tra 核苷酸變異,以保證Tra/Tra-2 復(fù)合體的形成,Tra-2 通過與Tra 協(xié)同進(jìn)化發(fā)揮其在昆蟲性別決定中的作用。

    References)

    Amrein H,Maniatis T,N?thiger R.Alternatively spliced transcripts of the sex determining gene tra-2 of Drosophila encode functional proteins of different size [J].EMBO Journal,1990,9:3619-3629.

    Bell LR,Horabin JI,Schedl P,et al.Positive auto-regulation of sex lethal by alternative splicing maintains the female determined state in Drosophila [J].Cell,1991,65:229-239.

    Bender W,Spierer P,Hogness DS.Chromosomal walking and jumping to isolate DNA from the Ace and rosy loci and the bithorax complex in Drosophila melanogaster [J].Journal of Molecular Biology,1983,168:17-33.

    Black DL.Mechanisms of alternative pre-messenger RNA splicing[J].Annual Review of Biochemistry,2003,72:291-336.

    Burghardt G,Hediger M,Siegenthaler C,et al.The transformer 2 gene in Musca domestica is required for selecting and maintaining the female pathway of development [J].Development Genes and Evolution,2005,215:165-176.

    Catteruccia F,Benton JP,Crisanti A.An anopheles transgenic sexing strain for vector control [J].Nature Biotechnology,2005,23:1414-1417.

    Concha C,Scott MJ.Sexual development in Lucilia cuprina (Diptera,Calliphoridae)is controlled by the transformer gene[J].Genetics,2009,182:785-789.

    Hall TA.BioEdit:a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT [J].Nucleic Acids Symposium Series,1999,41:95-98.

    Kimura M.A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences[J].Journal of Molecular Evolution,1980,16:11-120.

    Klassen W,Curtis CF.History of the sterile insect technique.In:Dyck V,Hendrichs J,Robinson AS,eds.Sterile Insect Technique[M].Netherlands:Springer,2005.

    Knipling EF.Possibilities of insect control or eradication through the use of sexually sterile males [J].Journal of Economic Entomology,1955,48:459e462.

    Kulathinal RJ,Skwarek L,Morton RA,et al.Rapid evolution of the sex determining gene,transformer:Structural diversity and rate heterogeneity among sibling species of Drosophila[J].Molecular Biology Evolution,2003,20:441-452.

    Liang GQ,Yang GH,Liang F,et al.Bactrocera Species in Asia-Pacific Area[M].Guangzhou:Guangdong Science and Technology Press,1996.[梁廣勤,楊國海,梁帆,等.亞太地區(qū)寡毛實(shí)蠅.廣州:廣東科技出版社,1996]

    Librado P,Rozas J.DnaSP v5:A software for comprehensive analysis of DNA polymorphism data [J].Bioinformatics,2009,25 (11):1451-1452.

    McGuffin ME,Chandler D,Somaiya D,et al.Auto-regulation of transformer-2 alternative splicing is necessary for normal male fertility in Drosophila[J].Genetics,1998,149:1 477-1 486.

    Nei M,Kumar S.Molecular Evolution and Phylogenetics [M].New York:Oxford University Press,2000.

    Rendón P,McInnis D,Lance D,et al.Medfly (Diptera:Tephritidae)genetic sexing:large-scale field comparison of males-only and bisexual sterile fly releases in Guatemala[J].Journal of Economic Entomology,2004,97:1547-1553.

    Ruiz MF,Milano A,Salvemini M,et al.The gene transformer of Anastrepha fruit flies (Diptera,tephritidae)and its evolution in insects[J].PLoS ONE,2007,2:e1239.

    Salvemini M,Robertson M,Aronson B,et al.Ceratitis capitata transformer-2 gene is required to establish and maintain the autoregulation of ctra,the master gene for female sex determination[J].International Journal of Developmental Biology,2009,53:109-120.

    Sánchez L.Sex-determining mechanisms in insects[J].International Journal of Developmental Biology,2008,52:837-856.

    Schetelig MF,Milano A,Saccone G,et al.Male only progeny in Anastrepha suspensa by RNAi-induced sex reversion of chromosomal females [J].Insect Biochemistry and Molecular Biology,2012,42 (1):51-57.

    Tamura K,Peterson D,Peterson N,et al.MEGA 5:Molecular evolutionary genetics analysis using maximum likelihood,evolutionary distance,and maximum parsimony methods [J].Molecular Biology and Evolution,2011,28:2731-2739.

    Thompson JD,Higgins DG,Gibson TJ,et al.CLUSTAL W.Improving the sensitivity of progressive multiple sequence alignments through sequence weighting,position specific gap penalties and weight matrix choice [J].Nucleic Acids Research,1994,22:4673-4680.

    Zhang J,Rosenberg HF,Nei M.Positive Darwinian selection after gene duplication in primate ribonuclease genes[J].Proceedings of the National Academy of Sciences of the United States of America,1998,95:3708-3713.

    猜你喜歡
    毛果實(shí)蠅同義
    Dale Carnegie
    住在樓上的老馬(三)
    西夏文《同義》重復(fù)字研究
    西夏學(xué)(2019年1期)2019-02-10 06:22:08
    西夏文《同義》考釋三則
    西夏學(xué)(2018年2期)2018-05-15 11:25:30
    瓜實(shí)蠅的發(fā)生與綜合防治
    蔬菜(2016年8期)2016-10-10 06:49:12
    實(shí)蠅蟲果悶殺袋
    柑橘大實(shí)蠅綜合治理
    微波輻射-溶劑回流法提取/HPLC法測(cè)定倒提壺中天芥菜堿與毛果天芥菜堿的含量
    柑橘實(shí)蠅防治一法
    毛果魚藤黃酮類化合物提取工藝及抑菌作用的研究
    中成藥(2014年11期)2014-02-28 22:30:08
    麻豆成人av在线观看| 成人特级av手机在线观看| 日本爱情动作片www.在线观看 | 给我免费播放毛片高清在线观看| 色尼玛亚洲综合影院| 少妇的逼水好多| av天堂中文字幕网| 久久这里只有精品中国| 欧美成人a在线观看| 成人鲁丝片一二三区免费| 国产av麻豆久久久久久久| 又紧又爽又黄一区二区| 丰满的人妻完整版| 美女cb高潮喷水在线观看| 成人亚洲精品av一区二区| 日韩欧美在线乱码| 欧美一区二区精品小视频在线| 亚州av有码| av视频在线观看入口| 中文字幕精品亚洲无线码一区| 啦啦啦观看免费观看视频高清| 亚洲人成网站在线播放欧美日韩| 99久久精品热视频| 欧美xxxx黑人xx丫x性爽| 最新在线观看一区二区三区| 欧美日韩瑟瑟在线播放| 国产高清激情床上av| 国产黄色小视频在线观看| 国产成人影院久久av| 欧美国产日韩亚洲一区| or卡值多少钱| 少妇的逼好多水| 色吧在线观看| 国产精品人妻久久久影院| 天堂√8在线中文| 美女高潮喷水抽搐中文字幕| 人人妻人人澡欧美一区二区| 亚洲第一区二区三区不卡| 看免费成人av毛片| 亚洲人成网站高清观看| 国产成人a区在线观看| 亚洲在线自拍视频| 少妇的逼水好多| 成人国产综合亚洲| 欧美最新免费一区二区三区| 亚洲一区二区三区色噜噜| 久久久午夜欧美精品| 91av网一区二区| 99久久中文字幕三级久久日本| 国产欧美日韩精品亚洲av| 国产精品自产拍在线观看55亚洲| 欧美一区二区亚洲| 欧美日韩乱码在线| 国产一区二区三区在线臀色熟女| 嫩草影院新地址| 日本一二三区视频观看| videossex国产| 最近在线观看免费完整版| 欧美色视频一区免费| 嫁个100分男人电影在线观看| 国模一区二区三区四区视频| 亚洲欧美清纯卡通| 日本 欧美在线| 好男人在线观看高清免费视频| 久久午夜亚洲精品久久| 欧美xxxx黑人xx丫x性爽| 俺也久久电影网| 欧美成人一区二区免费高清观看| 免费人成视频x8x8入口观看| 成人午夜高清在线视频| 免费观看的影片在线观看| 男插女下体视频免费在线播放| 在线观看午夜福利视频| 一本一本综合久久| 别揉我奶头~嗯~啊~动态视频| 麻豆国产av国片精品| 成年免费大片在线观看| 美女高潮的动态| 一卡2卡三卡四卡精品乱码亚洲| 国产真实乱freesex| 99riav亚洲国产免费| 偷拍熟女少妇极品色| 国产精品爽爽va在线观看网站| 嫩草影院精品99| 深爱激情五月婷婷| 最近视频中文字幕2019在线8| 五月玫瑰六月丁香| 美女大奶头视频| 18禁黄网站禁片免费观看直播| 午夜影院日韩av| 欧美又色又爽又黄视频| 22中文网久久字幕| 欧美三级亚洲精品| 伦理电影大哥的女人| 看免费成人av毛片| 国产亚洲精品av在线| 亚洲成av人片在线播放无| 午夜视频国产福利| 观看免费一级毛片| 精品久久久久久成人av| 一级av片app| 亚洲人成伊人成综合网2020| 日韩亚洲欧美综合| 久久天躁狠狠躁夜夜2o2o| 日本a在线网址| 一夜夜www| 亚洲国产精品成人综合色| 亚洲精品影视一区二区三区av| 欧美成人一区二区免费高清观看| 99riav亚洲国产免费| 精品欧美国产一区二区三| 国产精品一区www在线观看 | 男女那种视频在线观看| 免费一级毛片在线播放高清视频| 国产精品99久久久久久久久| 伦理电影大哥的女人| 又紧又爽又黄一区二区| 久久香蕉精品热| 麻豆久久精品国产亚洲av| 国产av在哪里看| 久久精品夜夜夜夜夜久久蜜豆| 免费av毛片视频| 精品乱码久久久久久99久播| 麻豆国产97在线/欧美| 精品午夜福利在线看| АⅤ资源中文在线天堂| 噜噜噜噜噜久久久久久91| 成人午夜高清在线视频| 中国美女看黄片| 久久久成人免费电影| 三级男女做爰猛烈吃奶摸视频| 欧美不卡视频在线免费观看| 色综合色国产| 色哟哟哟哟哟哟| 少妇被粗大猛烈的视频| 久久国产乱子免费精品| 18禁黄网站禁片免费观看直播| 亚洲av.av天堂| 亚洲精品色激情综合| 高清日韩中文字幕在线| 久久久久久久久久黄片| 免费高清视频大片| 丰满的人妻完整版| 国产一区二区激情短视频| 乱人视频在线观看| 天美传媒精品一区二区| 午夜老司机福利剧场| 国产真实乱freesex| 成人二区视频| 免费观看精品视频网站| 中文字幕av在线有码专区| 欧美潮喷喷水| 波多野结衣巨乳人妻| 国产精品国产高清国产av| 日日摸夜夜添夜夜添av毛片 | 国产精品乱码一区二三区的特点| 欧美日韩国产亚洲二区| av.在线天堂| av在线蜜桃| 国产伦精品一区二区三区视频9| 亚洲精品粉嫩美女一区| 中国美女看黄片| 亚洲av电影不卡..在线观看| 中文资源天堂在线| 亚洲中文字幕一区二区三区有码在线看| 自拍偷自拍亚洲精品老妇| 久久99热这里只有精品18| 欧美日韩中文字幕国产精品一区二区三区| 两个人视频免费观看高清| 久9热在线精品视频| 欧美成人免费av一区二区三区| 欧美一级a爱片免费观看看| 十八禁国产超污无遮挡网站| 亚洲精品456在线播放app | 黄色丝袜av网址大全| 观看免费一级毛片| 在线国产一区二区在线| 赤兔流量卡办理| 熟妇人妻久久中文字幕3abv| 色5月婷婷丁香| 亚洲av五月六月丁香网| 亚洲中文字幕一区二区三区有码在线看| 永久网站在线| 欧美又色又爽又黄视频| av在线观看视频网站免费| 午夜福利在线在线| 免费看a级黄色片| 精品人妻1区二区| av天堂在线播放| 日日啪夜夜撸| 午夜影院日韩av| 色在线成人网| 日日啪夜夜撸| 久久精品国产自在天天线| 一级a爱片免费观看的视频| 村上凉子中文字幕在线| 久久精品久久久久久噜噜老黄 | 国产 一区精品| 久久久久久久久大av| 亚洲自偷自拍三级| 女生性感内裤真人,穿戴方法视频| 欧美性猛交╳xxx乱大交人| 91在线精品国自产拍蜜月| 俺也久久电影网| 午夜精品久久久久久毛片777| 在线看三级毛片| 国产精品伦人一区二区| 欧美日韩亚洲国产一区二区在线观看| 别揉我奶头 嗯啊视频| 亚洲在线观看片| 亚洲国产日韩欧美精品在线观看| 国产精品亚洲美女久久久| 看十八女毛片水多多多| 色哟哟哟哟哟哟| 午夜福利在线在线| 国产不卡一卡二| or卡值多少钱| bbb黄色大片| 99久久成人亚洲精品观看| 亚洲午夜理论影院| 免费在线观看影片大全网站| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| 国产一区二区三区在线臀色熟女| 在线观看午夜福利视频| 成人鲁丝片一二三区免费| 精品一区二区三区av网在线观看| 欧美精品国产亚洲| 亚洲国产欧洲综合997久久,| 国产高清视频在线播放一区| 成人国产一区最新在线观看| 久久精品国产亚洲av香蕉五月| 少妇的逼好多水| 欧美zozozo另类| 亚洲av免费高清在线观看| 91精品国产九色| 国产高清三级在线| 欧美激情久久久久久爽电影| 日韩一本色道免费dvd| 精品乱码久久久久久99久播| 欧美性感艳星| 国产精品一区www在线观看 | 亚洲中文字幕日韩| av在线观看视频网站免费| 在线免费观看的www视频| 国产高清有码在线观看视频| 简卡轻食公司| 欧美成人免费av一区二区三区| 国产精品电影一区二区三区| 国产乱人伦免费视频| 国产三级在线视频| 日韩欧美三级三区| 一个人免费在线观看电影| 18禁裸乳无遮挡免费网站照片| 国产单亲对白刺激| 国产探花在线观看一区二区| 免费无遮挡裸体视频| 91久久精品电影网| 一区福利在线观看| 亚洲欧美日韩无卡精品| 在线免费十八禁| 久久人妻av系列| 日本熟妇午夜| 人人妻,人人澡人人爽秒播| 亚洲自偷自拍三级| 中国美女看黄片| 亚洲美女搞黄在线观看 | 日韩 亚洲 欧美在线| 欧美国产日韩亚洲一区| 色综合婷婷激情| 精品一区二区三区人妻视频| 美女高潮喷水抽搐中文字幕| a在线观看视频网站| 99热这里只有精品一区| 在线国产一区二区在线| 亚洲一区二区三区色噜噜| 18禁黄网站禁片免费观看直播| 免费看美女性在线毛片视频| 国产成人av教育| 国产探花在线观看一区二区| 成人特级黄色片久久久久久久| 午夜a级毛片| 一夜夜www| 午夜激情福利司机影院| 中文字幕久久专区| 国产午夜福利久久久久久| 久久中文看片网| 国模一区二区三区四区视频| 一级av片app| 黄色一级大片看看| 少妇丰满av| 中亚洲国语对白在线视频| 日本欧美国产在线视频| 亚洲人成网站在线播| 亚洲av二区三区四区| 在线观看一区二区三区| 69av精品久久久久久| 亚洲av电影不卡..在线观看| 国产精品乱码一区二三区的特点| 2021天堂中文幕一二区在线观| 久久久久久国产a免费观看| 久久午夜亚洲精品久久| 国产精品一区二区免费欧美| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影| 嫁个100分男人电影在线观看| 别揉我奶头 嗯啊视频| 91在线观看av| 国产高清视频在线播放一区| 色吧在线观看| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 成人永久免费在线观看视频| 97碰自拍视频| 欧美精品国产亚洲| 国产精品一区www在线观看 | 22中文网久久字幕| 免费高清视频大片| 99精品久久久久人妻精品| 窝窝影院91人妻| 国产乱人视频| 九九热线精品视视频播放| 18禁黄网站禁片午夜丰满| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 日韩av在线大香蕉| 国产日本99.免费观看| 久久久久久久久中文| 日韩人妻高清精品专区| 精品久久久久久久久av| 午夜激情福利司机影院| 精品国内亚洲2022精品成人| 小蜜桃在线观看免费完整版高清| 成年女人看的毛片在线观看| 免费不卡的大黄色大毛片视频在线观看 | 深夜精品福利| videossex国产| 国产aⅴ精品一区二区三区波| 三级男女做爰猛烈吃奶摸视频| 精品国产三级普通话版| 91久久精品国产一区二区成人| 一级黄片播放器| 精品久久久久久久人妻蜜臀av| 国产免费一级a男人的天堂| 美女高潮喷水抽搐中文字幕| 成人无遮挡网站| 免费黄网站久久成人精品| 成人高潮视频无遮挡免费网站| 亚洲精品日韩av片在线观看| 精品免费久久久久久久清纯| 热99re8久久精品国产| 国产aⅴ精品一区二区三区波| 在线a可以看的网站| 久久久成人免费电影| 日韩高清综合在线| 国产在视频线在精品| 亚洲人成网站在线播放欧美日韩| 国产免费av片在线观看野外av| 免费人成视频x8x8入口观看| 精品日产1卡2卡| 欧美日韩瑟瑟在线播放| 亚洲成人中文字幕在线播放| 亚洲av一区综合| 久久久国产成人精品二区| 亚洲国产日韩欧美精品在线观看| 在线观看av片永久免费下载| 男女边吃奶边做爰视频| 欧美xxxx性猛交bbbb| 午夜视频国产福利| 中文字幕av成人在线电影| 欧美一区二区国产精品久久精品| 美女 人体艺术 gogo| 国产精品亚洲一级av第二区| 国产一区二区三区在线臀色熟女| 成人二区视频| 欧美人与善性xxx| 亚洲成人精品中文字幕电影| 免费大片18禁| 老司机午夜福利在线观看视频| 特级一级黄色大片| 少妇被粗大猛烈的视频| 亚洲黑人精品在线| 黄色丝袜av网址大全| 中文亚洲av片在线观看爽| 欧美黑人欧美精品刺激| 精品一区二区三区视频在线| 国产精品野战在线观看| 日本五十路高清| 国内毛片毛片毛片毛片毛片| 国产单亲对白刺激| 伦理电影大哥的女人| 免费不卡的大黄色大毛片视频在线观看 | 97超级碰碰碰精品色视频在线观看| 亚洲成av人片在线播放无| 午夜影院日韩av| 免费看美女性在线毛片视频| 色视频www国产| 黄色配什么色好看| 久久久成人免费电影| 久久国内精品自在自线图片| 啦啦啦韩国在线观看视频| 国产三级中文精品| 制服丝袜大香蕉在线| 亚洲男人的天堂狠狠| 国产精品亚洲一级av第二区| 国产男人的电影天堂91| 一进一出抽搐动态| 日日摸夜夜添夜夜添av毛片 | 国产淫片久久久久久久久| 成人精品一区二区免费| 亚洲成人精品中文字幕电影| 丝袜美腿在线中文| 黄色丝袜av网址大全| 国产免费一级a男人的天堂| 九九热线精品视视频播放| 久久久国产成人精品二区| 色综合色国产| 88av欧美| 无人区码免费观看不卡| 一个人看视频在线观看www免费| 日本 av在线| 禁无遮挡网站| 看片在线看免费视频| 亚洲成人精品中文字幕电影| 丝袜美腿在线中文| 十八禁国产超污无遮挡网站| 亚洲专区中文字幕在线| 观看美女的网站| 精品久久久久久久久久久久久| 国产极品精品免费视频能看的| 男女边吃奶边做爰视频| 男人狂女人下面高潮的视频| 在线观看免费视频日本深夜| 国产精品三级大全| 精品日产1卡2卡| 不卡一级毛片| 日韩强制内射视频| 女的被弄到高潮叫床怎么办 | 观看美女的网站| 国产在视频线在精品| 亚洲 国产 在线| 亚洲aⅴ乱码一区二区在线播放| 观看免费一级毛片| 亚洲乱码一区二区免费版| 午夜老司机福利剧场| 看片在线看免费视频| 欧美日韩精品成人综合77777| 久久久久九九精品影院| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 中文字幕精品亚洲无线码一区| 国产探花极品一区二区| 国产欧美日韩精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 一进一出抽搐gif免费好疼| 欧美成人a在线观看| 亚洲va日本ⅴa欧美va伊人久久| x7x7x7水蜜桃| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看| 天堂av国产一区二区熟女人妻| eeuss影院久久| 成人性生交大片免费视频hd| 国产黄片美女视频| 国产一区二区激情短视频| 精品久久久久久久久久久久久| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 国产午夜精品论理片| 长腿黑丝高跟| 亚洲精品乱码久久久v下载方式| 婷婷精品国产亚洲av在线| 午夜老司机福利剧场| 亚洲精华国产精华液的使用体验 | 色视频www国产| 又紧又爽又黄一区二区| 狂野欧美激情性xxxx在线观看| 国产精品永久免费网站| 免费看a级黄色片| 亚洲欧美日韩高清专用| 欧美绝顶高潮抽搐喷水| 国产一区二区在线av高清观看| 91精品国产九色| 91久久精品电影网| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人福利小说| 别揉我奶头 嗯啊视频| 日韩欧美国产一区二区入口| 午夜免费激情av| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 国产精品一区二区免费欧美| 亚洲精品成人久久久久久| 国产日本99.免费观看| 精品不卡国产一区二区三区| 欧美性感艳星| 精品久久久噜噜| 国产高清激情床上av| 久久久国产成人精品二区| 欧美+亚洲+日韩+国产| 久久午夜亚洲精品久久| 波野结衣二区三区在线| 搡老妇女老女人老熟妇| 赤兔流量卡办理| 久久久久久国产a免费观看| av在线老鸭窝| 免费av毛片视频| 国产v大片淫在线免费观看| 97碰自拍视频| 男人和女人高潮做爰伦理| 看十八女毛片水多多多| 欧美性猛交╳xxx乱大交人| 人妻丰满熟妇av一区二区三区| 国产不卡一卡二| 久久99热6这里只有精品| 身体一侧抽搐| 日韩欧美精品免费久久| 免费无遮挡裸体视频| 听说在线观看完整版免费高清| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 99久国产av精品| 国产精品福利在线免费观看| 国产真实乱freesex| 麻豆成人av在线观看| 一区二区三区四区激情视频 | 久久婷婷人人爽人人干人人爱| 男人狂女人下面高潮的视频| 国产精品久久久久久亚洲av鲁大| 1024手机看黄色片| 三级男女做爰猛烈吃奶摸视频| 91在线精品国自产拍蜜月| 亚洲无线观看免费| 国产一区二区亚洲精品在线观看| 国产色爽女视频免费观看| 日本色播在线视频| 国产精品美女特级片免费视频播放器| 日本成人三级电影网站| 免费搜索国产男女视频| 亚洲国产日韩欧美精品在线观看| 亚洲自偷自拍三级| 亚洲天堂国产精品一区在线| 国产精品av视频在线免费观看| 精品午夜福利视频在线观看一区| 精品人妻1区二区| 麻豆成人午夜福利视频| 欧美区成人在线视频| 亚洲黑人精品在线| av在线亚洲专区| 国产精品自产拍在线观看55亚洲| 免费人成在线观看视频色| 99热只有精品国产| 老女人水多毛片| 久久久成人免费电影| 欧美另类亚洲清纯唯美| 色在线成人网| 性欧美人与动物交配| 成人国产一区最新在线观看| 麻豆国产av国片精品| av在线蜜桃| 真人一进一出gif抽搐免费| 成年女人看的毛片在线观看| 毛片女人毛片| 精品一区二区三区视频在线观看免费| 一区二区三区免费毛片| av福利片在线观看| 在线观看午夜福利视频| 国产熟女欧美一区二区| 国产亚洲精品久久久com| 成人特级黄色片久久久久久久| 美女cb高潮喷水在线观看| 国产老妇女一区| 精品久久久久久久末码| 国产乱人伦免费视频| 亚洲精华国产精华精| 黄色一级大片看看| 啪啪无遮挡十八禁网站| 白带黄色成豆腐渣| 一区二区三区免费毛片| 夜夜看夜夜爽夜夜摸| 久久久色成人| 国内精品美女久久久久久| 国产真实乱freesex| 99国产极品粉嫩在线观看| 国产视频内射| 午夜福利18| 成人亚洲精品av一区二区| 少妇丰满av| 99久久成人亚洲精品观看| 亚州av有码| 午夜视频国产福利| 日本免费a在线| 婷婷亚洲欧美| 午夜福利在线观看吧| 亚洲中文字幕日韩| 1024手机看黄色片| 在线观看av片永久免费下载| 无人区码免费观看不卡| 日本五十路高清| 国产v大片淫在线免费观看| 成人国产麻豆网| 国产综合懂色| 国产欧美日韩精品亚洲av| 真人做人爱边吃奶动态| 欧美一区二区国产精品久久精品| 国产精品人妻久久久久久| 我要搜黄色片| 在线观看66精品国产| 国产精品一区www在线观看 | 亚洲黑人精品在线| 动漫黄色视频在线观看| 日韩高清综合在线| 精品福利观看|