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

    基于線粒體COⅠ基因的竹筍夜蛾親緣關(guān)系*

    2017-06-05 09:07:55莫有迪張亞波舒金平王浩杰
    林業(yè)科學(xué) 2017年4期
    關(guān)鍵詞:親緣進(jìn)化樹(shù)堿基

    張 琴 莫有迪 張亞波 舒金平 王浩杰 吳 鴻

    (1.中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所 富陽(yáng) 311400; 2.杭州市林業(yè)科學(xué)研究院 杭州 311300;3.廣東省肇慶市懷集縣林業(yè)局 肇慶 526400)

    ?

    基于線粒體COⅠ基因的竹筍夜蛾親緣關(guān)系*

    張 琴1,2莫有迪3張亞波1舒金平1王浩杰1吳 鴻1

    (1.中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所 富陽(yáng) 311400; 2.杭州市林業(yè)科學(xué)研究院 杭州 311300;3.廣東省肇慶市懷集縣林業(yè)局 肇慶 526400)

    【目的】 基于mtDNA COⅠ基因探究竹筍夜蛾的親緣關(guān)系,以確定其分類(lèi)地位,為系統(tǒng)進(jìn)化軌跡提供分子理論基礎(chǔ)?!痉椒ā?以分布于浙江杭州周邊竹子產(chǎn)區(qū)的4種竹筍夜蛾為研究對(duì)象,使用通用引物分別克隆4種竹筍夜蛾線粒體細(xì)胞色素C 氧化酶亞基Ⅰ(mtDNA COⅠ) 基因5′端709 bp 片段序列,通過(guò)構(gòu)建系統(tǒng)進(jìn)化樹(shù),分析4種竹筍夜蛾親緣關(guān)系。【結(jié)果】 1) 4種竹筍夜蛾遺傳距離較??; 2)與比對(duì)昆蟲(chóng)相比,4種竹筍夜蛾同屬一個(gè)進(jìn)化群,分支支持率達(dá)到100%?!窘Y(jié)論】 4種竹筍夜蛾親緣關(guān)系接近,除薩夜蛾外,建議將其他3種竹筍夜蛾歸為一個(gè)屬。4種竹筍夜蛾形態(tài)差異較小,生態(tài)習(xí)性相似,分布區(qū)域基本一致,可作為理想的同域物種進(jìn)化研究素材。本研究可為4種竹筍夜蛾分類(lèi)地位的確定提供初步分子基礎(chǔ),同時(shí)可為竹筍夜蛾DNA條形碼鑒定提供理論基礎(chǔ)。

    竹筍夜蛾; COⅠ基因; 序列比較; 親緣關(guān)系

    竹筍夜蛾屬鱗翅目(Lepidoptera)夜蛾科(Noctuidae),是竹子筍期的主要害蟲(chóng)之一。竹筍夜蛾鉆蛀取食竹筍,活動(dòng)場(chǎng)所隱蔽,且幼蟲(chóng)期較長(zhǎng),竹筍被危害嚴(yán)重時(shí),被害率高達(dá)95%以上。被害竹筍能生長(zhǎng)成竹者,亦為斷頭、折梢、蟲(chóng)孔累累、材質(zhì)硬脆,竹材利用率低(徐天森, 1962; 徐天森等, 2004)。被害竹林大量退筍,使竹林日趨退化,影響筍用林產(chǎn)業(yè)的發(fā)展(胡正堅(jiān)等, 1992; 梁光紅等; 2003; Yoshimatsuetal., 2005; 舒金平等, 2012; 葉碧歡等, 2014)。在我國(guó),研究報(bào)道的竹筍夜蛾主要有基夜蛾(Kumasiakumaso)、竹筍禾夜蛾(Oligiavulgaris)、筍秀夜蛾(Apameaapameoides)和薩夜蛾(Sapporiarepetita)(徐天森等, 2004; 舒金平等, 2012)。這4種竹筍夜蛾均鉆蛀取食竹筍,屬寡食性昆蟲(chóng),生活史極其相似,生態(tài)位嚴(yán)重重疊,在時(shí)間、空間和營(yíng)養(yǎng)等資源利用上存在激烈競(jìng)爭(zhēng)(羅梅浩等, 2002; 黃瓊瑤等, 2009),且在同一地區(qū)均有分布(黃煥華等, 2003)。此外,4種竹筍夜蛾成蟲(chóng)體色差異小,大小接近,實(shí)際生產(chǎn)生活中難以區(qū)分,人們常將其混淆,統(tǒng)稱為竹筍夜蛾,導(dǎo)致研究和生產(chǎn)防治中對(duì)象模糊,目標(biāo)不清。

    傳統(tǒng)昆蟲(chóng)分類(lèi)學(xué)對(duì)這4種竹筍夜蛾的分類(lèi)一直存在爭(zhēng)議,拉丁學(xué)名一再變更,屬的歸屬難以確定(徐天森等, 2004)。例如: 基夜蛾,曾用名淡竹筍夜蛾,其拉丁學(xué)名曾為Apamea(秀夜蛾屬)kumaso(Sugi)(朱弘復(fù), 1964; Inoueetal,1982; 陳一心, 1985; 1992);筍秀夜蛾,曾用名筍秀禾夜蛾,其拉丁學(xué)名曾為Oligia(禾夜蛾屬)apameoidis(Draudt)(陳一心, 1985);竹筍禾夜蛾,曾用名竹筍夜蛾,其拉丁學(xué)名曾為Polydesma(紛夜蛾屬)vulgaris(Butler)、Oligiavulgaris(Hampson)、Bambusiphilavulgaris(Sugi)(朱弘復(fù), 1964; Inoueetal., 1982; 陳一心, 1985; 1992);薩夜蛾,曾用名筍連秀夜蛾,其拉丁學(xué)名曾為Apamea(秀夜蛾屬)repetita(Butler)、Trachea(陌夜蛾屬)repetita(Hampson)、Apamea(秀夜蛾屬)repetita(Chen)、Sapporia(薩夜蛾屬)repetita(Sugi)(Inoueetal., 1982; 陳一心, 1985; 1992)。在研究中更顯混亂,竹筍夜蛾的學(xué)名難以統(tǒng)一。

    隨著分子生物學(xué)的快速發(fā)展,基于DNA所攜帶的遺傳信息進(jìn)行生物分類(lèi)已逐漸成為研究熱點(diǎn)。其中,昆蟲(chóng)線粒體基因組中細(xì)胞色素氧化酶亞基Ⅰ(mtDNA COⅠ)基因序列因既具有高度保守序列,又具有足夠變異位點(diǎn)的特性,被廣泛用于昆蟲(chóng)分類(lèi)、種群和進(jìn)化關(guān)系研究(Luntetal., 1996)。2002年,德國(guó)進(jìn)化生物學(xué)家Tautz等首先提出將DNA作為生物分類(lèi)研究平臺(tái)的觀點(diǎn),即以DNA 序列為基礎(chǔ)建立物種鑒定體系(Tautzetal., 2003)。2003年,加拿大學(xué)者Hebert 等提出建立DNA條形碼鑒定系統(tǒng),即以線粒體細(xì)胞色素氧化酶亞基I(COⅠ) 基因5′端648 bp 的序列多樣性為基礎(chǔ)進(jìn)行物種鑒定。Hebert 等(2003)還發(fā)現(xiàn)這段序列在鱗翅目200個(gè)近緣種中鑒定的成功率為100%。本研究針對(duì)4種竹筍夜蛾近緣種生物學(xué)和生態(tài)學(xué)特性相似、分類(lèi)地位上存在爭(zhēng)議的問(wèn)題,克隆4種竹筍夜蛾線粒體COⅠ基因5′端709 bp片段系列,通過(guò)構(gòu)建系統(tǒng)進(jìn)化樹(shù),分析其親緣關(guān)系,為其進(jìn)一步明確其分類(lèi)地位提供分子理論基礎(chǔ)。

    1 材料與方法

    1.1 供試蟲(chóng)源

    2013年4月—2014年4月,在浙江省安吉市竹博園和富陽(yáng)竹子產(chǎn)區(qū)采集4種供試竹筍夜蛾幼蟲(chóng)各200頭,室內(nèi)飼養(yǎng)至成蟲(chóng)羽化,拍照,搜集活體成蟲(chóng)浸泡于無(wú)水乙醇中,-20 ℃保存?zhèn)溆谩?/p>

    1.2 試驗(yàn)材料

    TAKALA MiniBEST Universal Genomic DNA Extraction Kit VER.5.0 試劑盒,TaKaRa MiniBest DNA Fragment Purification Kit Ver 3.0試劑盒,PMD-19-T載體,大腸桿菌(Escherichiacoli)感受態(tài)細(xì)胞,其中擴(kuò)增4種竹筍夜蛾的線粒體COⅠ基因所用引物來(lái)自通用引物(Folmeretal., 1994):

    L1490-5’-GGTCAACAAATCATAAAGATATT G-3’;

    H2198-5’-TAAACTTCAGGGTGACCAAAA AA T-3’。

    引物由生工生物工程(上海)股份有限公司合成。

    1.3 試驗(yàn)方法

    1.3.1 4種竹筍夜蛾總基因組DNA的提取 分別選取4種竹筍夜蛾成蟲(chóng)樣本,按照TAKALA MiniBEST Universal Genomic DNA Extraction Kit VER.5.0 試劑盒操作步驟提取4種竹筍夜蛾總基因組DNA。產(chǎn)物經(jīng)1%瓊脂糖凝膠進(jìn)行檢測(cè),使用紫外成像系統(tǒng)對(duì)電泳條帶進(jìn)行觀察和拍照。

    1.3.2 4種竹筍夜蛾線粒體COⅠ基因的擴(kuò)增 以提取到的竹筍夜蛾總基因組DNA為模板,引物為上述通用引物,按照如下反應(yīng)體系和反應(yīng)條件對(duì)4種竹筍夜蛾線粒體COⅠ基因進(jìn)行擴(kuò)增。反應(yīng)體系為:25 μL,其中上下游引物各1 μL,10×擴(kuò)增緩沖液 2.5 μL,4種dNTP混合物2 μL,DNA模板5 μL,Taq DNA聚合酶0.3 μL,ddH2O 13.2 μL。反應(yīng)條件為: 預(yù)變性94 ℃ 2 min,變性94 ℃1 min,退火50 ℃1 min,延伸72 ℃1 min,進(jìn)行35個(gè)循環(huán)。取3 μL PCR擴(kuò)增產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測(cè),使用紫外成像系統(tǒng)對(duì)電泳條帶進(jìn)行觀察和拍照。

    1.3.3 4種竹筍夜蛾線粒體COⅠ基因PCR產(chǎn)物純化、克隆和測(cè)序 產(chǎn)物經(jīng)電泳后,按照TaKaRa MiniBest DNA Fragment Purification Kit Ver 3.0試劑盒操作步驟對(duì)4種竹筍夜蛾線粒體COⅠ基因PCR產(chǎn)物進(jìn)行純化,純化的PCR 產(chǎn)物與PMD-19-T載體進(jìn)行連接反應(yīng),16 ℃過(guò)夜。將連接反應(yīng)的產(chǎn)物轉(zhuǎn)化至大腸桿菌JM109 感受態(tài)細(xì)胞,待菌落長(zhǎng)出來(lái)后,挑選單個(gè)陽(yáng)性菌落進(jìn)行菌落PCR驗(yàn)證,鑒定重組質(zhì)粒。測(cè)序公司為生工生物工程(上海)測(cè)序部,要求雙向測(cè)通并進(jìn)行拼接。

    1.3.4 數(shù)據(jù)處理 利用MEGA 5.0軟件分析4種竹筍夜蛾的堿基組成、結(jié)構(gòu)特點(diǎn)、堿基替換模式和遺傳距離等特性。登錄NCBI對(duì)4種竹筍夜蛾線粒體COⅠ基因進(jìn)行同源性搜索,選取其他21條(15屬21種)線粒體COⅠ基因(表1)的氨基酸序列,分別采用鄰接法(NJ)、最大簡(jiǎn)約法(ML)和最大似然法(MP)構(gòu)建系統(tǒng)發(fā)育樹(shù)。分析其系統(tǒng)進(jìn)化關(guān)系,其中NJ法選用的模型是P-Distance。所有方法均運(yùn)用Bootstrap檢驗(yàn)重復(fù),抽樣1 000次,其余參數(shù)為默認(rèn)。

    表1 21種蛾類(lèi)昆蟲(chóng)COⅠ 基因信息①Tab.1 Information of COⅠ genes of 21 noctuid insects

    ① *:無(wú)中文屬名No Chinese name.

    2 結(jié)果與分析

    2.1 4種竹筍夜蛾成蟲(chóng)形態(tài)特征

    由圖1可知,4種竹筍夜蛾成蟲(chóng)形態(tài)基本接近,主要區(qū)別在于其體色和斑紋。4種竹筍夜蛾成蟲(chóng)形態(tài)相關(guān)指標(biāo)見(jiàn)表2。

    2.2 4種竹筍夜蛾COⅠ基因目的片段序列的克隆與鑒定

    對(duì)4種竹筍夜蛾COⅠ基因片段克隆產(chǎn)物進(jìn)行電泳檢測(cè),每個(gè)擴(kuò)增產(chǎn)物均出現(xiàn)1條700 bp左右的特異性條帶(圖2),與試驗(yàn)?zāi)康南喾?。樣品?jīng)測(cè)序公司測(cè)序拼接后,獲得4種竹筍夜蛾mtDNA COⅠ基因片段堿基序列,每個(gè)樣品重復(fù)測(cè)序3次,結(jié)果相同,長(zhǎng)度均為709 bp(登錄號(hào)分別為基夜蛾KR780071、薩夜蛾KR780072、竹筍禾夜蛾KR780073、筍秀夜蛾KR780074)。將該片段序列在NCBI上進(jìn)行同源性搜索,結(jié)果同比對(duì)的夜蛾科中Noctuaatlantic(AY600452.1)等昆蟲(chóng)的線粒體COⅠ基因序列相似度極高(94%以上)。

    圖1 4種竹筍夜蛾成蟲(chóng)形態(tài)Fig.1 Adult morphological characteristic of 4 bamboo shoot-boring noctuids

    表2 4種竹筍夜蛾成蟲(chóng)形態(tài)特征Tab.2 Adult morphological characteristic of 4 bamboo shoot-boring noctuids

    圖2 竹筍夜蛾mtDNA COⅠ基因克隆陽(yáng)性菌落PCR電泳Fig 2 MtDNA COⅠgene’s positive colony PCR electrophoregram of bamboo shoot borers

    2.3 4種竹筍夜蛾及其比對(duì)昆蟲(chóng)mtDNA COⅠ基因片段的序列特征

    4種竹筍夜蛾mtDNA COⅠ基因片段堿基序列同源性達(dá)98.2%。在4種竹筍夜蛾樣本mtDNA COⅠ基因709 bp序列片段中,多態(tài)性位點(diǎn)為30個(gè),僅占整個(gè)序列的4.2%,與其他物種同源基因序列進(jìn)行堿基比對(duì)發(fā)現(xiàn),獲得的4條樣本基因片段均為編碼區(qū)序列,并且沒(méi)有堿基插入和缺失。氨基酸序列比對(duì)發(fā)現(xiàn),4種竹筍夜蛾的mtDNA COⅠ基因片段均從第3個(gè)堿基開(kāi)始翻譯,無(wú)移碼現(xiàn)象,編碼235個(gè)氨基酸,同源性高達(dá)99.89%,其中僅薩夜蛾第162個(gè)氨基酸與其他3種夜蛾不同,其余3種夜蛾氨基酸片段完全相同。

    2.4 4種竹筍夜蛾mtDNA COⅠ基因片段的遺傳變異分析

    使用MEG5.0軟件計(jì)算4種竹筍夜蛾及其比對(duì)昆蟲(chóng)mtDNA COⅠ基因序列的遺傳距離,計(jì)算模型為p-Distance。結(jié)果表明,筍秀夜蛾與基夜蛾間的遺傳距離為0.014,與竹筍禾夜蛾間的遺傳距離為0.006,與薩夜蛾間的遺傳距離為0.069; 基夜蛾與竹筍禾夜蛾間的遺傳距離為0.008,與薩夜蛾間的遺傳距離為0.067; 竹筍禾夜蛾與薩夜蛾間的遺傳距離為0.063(表3)。4種竹筍夜蛾中,除了薩夜蛾與其他3種夜蛾的種間遺傳距離較大,大于所有比對(duì)昆蟲(chóng)的平均種間距離外,其他3種種間遺傳距離均較小,與比對(duì)昆蟲(chóng)的種內(nèi)距離相近。

    4種竹筍夜蛾及其比對(duì)昆蟲(chóng)的堿基替換模式為: Ts >Tv,R值為12.6(表4)。TC間發(fā)生堿基轉(zhuǎn)換次數(shù)最多,達(dá)16次,占替換總數(shù)(19)的84.2%; TA發(fā)生堿基顛換次數(shù)最多,達(dá)13次,占顛換總數(shù)(16)的81.25%, CA和TG間的顛換次數(shù)較少,分別為2次和1次,CG間沒(méi)有替換發(fā)生。

    表3 4種竹筍夜蛾種間遺傳距離Tab 3 Interspecific genetic distance between the our kinds of bamboo shoot borers

    表4 4種竹筍夜蛾及其比對(duì)昆蟲(chóng)mtDNA COⅠ基因的堿基替換頻次①Tab 4 Base substitution frequencies in mitochondrial DNA COⅠ gene from bamboo Shoot borers and other lepidopterous insects

    ① Ts堿基轉(zhuǎn)換數(shù)Transition;Tv:堿基顛換數(shù)Transversion.R=Ts/Tv.

    圖3 基于mtDNA COⅠ基因構(gòu)建的4種竹筍夜蛾系統(tǒng)進(jìn)化樹(shù)(NJ法)Fig.3 Molecular phylogenetic tree of mtDNA COⅠ for bamboo shoot borers sequence with methods of Neighbor-Joining

    在4種竹筍夜蛾樣本mtDNA COⅠ基因709 bp序列片段中,檢測(cè)到保守性位點(diǎn)(conserved sites)555個(gè),多態(tài)性位點(diǎn)(variable sites)154個(gè),占總檢測(cè)位點(diǎn)的21.72%; 簡(jiǎn)約信息位點(diǎn)(parsimony-informative sites)118個(gè),占16.64%; 單突變位點(diǎn)(singleton site)36個(gè),占5.07%。

    2.5 4種竹筍夜蛾基于mtDNA COⅠ基因片段的系統(tǒng)進(jìn)化分析

    用4種竹筍夜蛾及其比對(duì)昆蟲(chóng)的mtDNA COⅠ氨基酸序列在Mega5.0中構(gòu)建進(jìn)化樹(shù),結(jié)果如圖3、4、5所示,3種進(jìn)化樹(shù)顯示結(jié)果大致相同。與比對(duì)序列相比,4種竹筍夜蛾分屬同一支,分支支持率達(dá)100%,親緣關(guān)系最近,其中薩夜蛾與其他3種夜蛾分屬2個(gè)分支,分支支持率為70%左右。MP法進(jìn)化樹(shù)中,基夜蛾又與竹筍禾夜蛾和筍秀夜蛾分為2支,分支支持率為30%。3種進(jìn)化樹(shù)均顯示,4種竹筍夜蛾與禾夜蛾屬、秀夜蛾屬的昆蟲(chóng)遺傳距離較遠(yuǎn)。

    圖4 基于mtDNA COⅠ基因構(gòu)建的4種竹筍夜蛾系統(tǒng)進(jìn)化樹(shù)(ML法)Fig.4 Molecular phylogenetic tree of mtDNA COⅠ for bamboo shoot borers sequence with methods of Maximum Likelihood

    3 討論

    Foster等(1997)認(rèn)為,線粒體基因密碼子位點(diǎn)在物種進(jìn)化過(guò)程中由于受選擇壓力影響,會(huì)表現(xiàn)出不同程度的堿基突變。在不同的密碼子位點(diǎn)上,4種竹筍夜蛾及其比對(duì)昆蟲(chóng)mtDNA COⅠ基因片段堿基含量變化規(guī)律一致。但密碼子第2位點(diǎn)的A+T含量顯著高于第1、3 位點(diǎn),說(shuō)明4種竹筍夜蛾密碼子第2位點(diǎn)上的堿基突變由純化選擇導(dǎo)致(陳星等, 2012)。此外,密碼子第2位點(diǎn)的堿基突變大多屬于同義突變,容易在進(jìn)化過(guò)程中通過(guò)自然選擇固定下來(lái)(周繼亮等, 2001)。本研究結(jié)果與櫟黃掌舟蛾(Phaleraassimilis)mtDNA COⅠ基因密碼子位點(diǎn)變異特點(diǎn)不同,其密碼子第3位點(diǎn)的A+T含量顯著高于第1、2 位點(diǎn)(楊瑞生等, 2013),可能是由于夜蛾種類(lèi)不同導(dǎo)致。在堿基替換模式上,本研究中4種竹筍夜蛾的mtDNA COⅠ基因均為堿基轉(zhuǎn)換高于顛換,其R值為12.6,與鱗翅目蛾類(lèi)昆蟲(chóng)的堿基替換模式一致(楊瑞生等, 2008; 2011; 2013)。生物進(jìn)化可以說(shuō)是生物體在基因水平上堿基突變不斷累積的結(jié)果,一旦同種生物不同個(gè)體間的堿基突變累積到一定程度(量變),就會(huì)產(chǎn)生種的分化(質(zhì)變)。Hebert(2003)等對(duì)動(dòng)物界11門(mén)13 320種的COⅠ 基因序列進(jìn)行比較研究發(fā)現(xiàn),超過(guò)98%的物種間的遺傳距離大于0.002,而種內(nèi)遺傳距離多小于0.001,很少超過(guò)0.002,而同屬種之間的mtDNA COⅠ基因平均距離為0.011 3,除刺細(xì)胞動(dòng)物門(mén)的物種外,COⅠ基因的序列多樣性足以區(qū)分所有動(dòng)物中的近緣種。本研究結(jié)果表明,4種竹筍夜蛾及其比對(duì)昆蟲(chóng)種間遺傳距離為0.002~0.093,符合Hebert等(2003)的研究結(jié)果。

    圖5 基于mtDNA COⅠ基因構(gòu)建的4種竹筍夜蛾系統(tǒng)進(jìn)化樹(shù)(MP法)Fig.5 Molecular phylogenetic tree of mtDNA COⅠ for bamboo shoot borers sequence with methods of maximum parsimony(MP)

    基于COⅠ基因建立的系統(tǒng)發(fā)育樹(shù)顯示,4種夜蛾親緣關(guān)系非常接近,而與秀夜蛾屬和禾夜蛾屬其他昆蟲(chóng)親緣關(guān)系較遠(yuǎn),因未能檢索到基夜蛾屬和薩夜蛾屬其他昆蟲(chóng)COⅠ基因,故無(wú)法比較其親緣關(guān)系。分類(lèi)學(xué)上對(duì)4種竹筍夜蛾的歸屬研究一直存在爭(zhēng)議,本研究或可提供部分參考。本研究中筍秀夜蛾與基夜蛾、竹筍禾夜蛾、薩夜蛾的遺傳距離分別為0.014、0.006、0.069,所建進(jìn)化樹(shù)也表明筍秀夜蛾與竹筍禾夜蛾親緣關(guān)系較近,其次為基夜蛾,與薩夜蛾親緣關(guān)系較遠(yuǎn)。綜上所述,本研究建議可將基夜蛾、竹筍禾夜蛾和筍秀夜蛾歸為同一個(gè)屬,薩夜蛾由于其在浙江分布較少,主要分布于廣東茶稈竹上,與這3種竹筍夜蛾親緣關(guān)系相對(duì)較遠(yuǎn),其分類(lèi)地位可進(jìn)一步進(jìn)行確認(rèn)。此外,4種竹筍夜蛾成蟲(chóng)形態(tài)差異較小,生態(tài)習(xí)性相似,分布區(qū)域基本一致,分子進(jìn)化樹(shù)初步確定其親緣關(guān)系接近,可作為理想的同域物種進(jìn)化研究素材。

    昆蟲(chóng)線粒體DNA 相對(duì)于核基因來(lái)說(shuō),具有進(jìn)化速率快、物種間序列保守等特點(diǎn),因而成為篩選條形碼標(biāo)準(zhǔn)基因的首選(Hebertetal., 2003)。研究發(fā)現(xiàn), 昆蟲(chóng)線粒體DNA含有13個(gè)蛋白質(zhì)編碼基因(Boore, 1999; Taanman, 1999),其中因?yàn)閙tDNA COⅠ基因的序列保守性強(qiáng),在種間鑒定上具有較高的可靠性,且易于擴(kuò)增,被絕大多數(shù)研究者選為昆蟲(chóng)的DNA 條形碼基因。利用DNA條形碼鑒定的研究中,63.7%的研究樣本來(lái)自昆蟲(chóng),其中又以鱗翅目的研究居多(李青青等, 2010),其次為雙翅目(Nelsonetal., 2007)、鞘翅目(Monaghanetal., 2005)等。本研究克隆了4種竹筍夜蛾mtDNA COⅠ基因片段,可為竹筍夜蛾DNA條形碼鑒定提供理論基礎(chǔ)。線粒體基因上保守基因眾多,多基因片段聯(lián)合的DNA 條形碼鑒定系統(tǒng)已經(jīng)成為一個(gè)探索的趨勢(shì)(肖金花等, 2004),進(jìn)一步驗(yàn)證DNA條形碼技術(shù)在物種鑒定和物種親緣關(guān)系研究中的適用性。此外,由于4種竹筍夜蛾親緣關(guān)系較近,今后可通過(guò)采集不同地理種群,克隆其mtDNA COⅠ基因及其他保守序列(COⅡ、16sRNA等),研究其進(jìn)化過(guò)程,進(jìn)一步揭示4種竹筍夜蛾的進(jìn)化關(guān)系。

    4 結(jié)論

    本研究利用4種竹筍夜蛾mtDNA COⅠ基因片段構(gòu)建了進(jìn)化樹(shù),結(jié)果表明4種竹筍夜蛾親緣關(guān)系接近,可為4種竹筍夜蛾分類(lèi)地位的確定提供初步分子基礎(chǔ),同時(shí)可為竹筍夜蛾DNA條形碼鑒定提供理論基礎(chǔ)。鑒于本研究克隆的只是4種竹筍夜蛾mtDNA COⅠ基因的5′端部分序列,今后可通過(guò)設(shè)計(jì)特異性引物,對(duì)竹筍夜蛾mtDNA COⅠ序列全長(zhǎng)進(jìn)行擴(kuò)增,以提高mtDNA COⅠ 基因DNA條形碼的可靠性。4種竹筍夜蛾形態(tài)差異較小,生態(tài)習(xí)性相似,分布區(qū)域基本一致,可作為理想的同域物種進(jìn)化研究素材。

    陳一心. 1985.中國(guó)經(jīng)濟(jì)昆蟲(chóng)志.第32冊(cè).鱗翅目夜蛾科.北京:科學(xué)出版社,1-167.

    (Chen Y X. 1985. Economic insect fauna of China. Chapter 32. Lepidoptera: Noctuidae. Beijing: Science Press. [in Chinese])

    陳一心. 1992. 鱗翅目: 夜蛾.長(zhǎng)沙: 湖南科學(xué)技術(shù)出版社, 947-1040.

    (Chen Y X.1992. Lepidoptera: Noctuidae. Changsha: Hunan Science and Technology Publishing House.[in Chinese])

    陳 星, 沈永義, 張亞平. 2012. 線粒體DNA在分子進(jìn)化研究中的應(yīng)用. 動(dòng)物學(xué)研究, 33(6): 566-573.

    (Chen X, Shen Y Y, Zhang Y P. 2012. Review of mtDNA in molecular evolution studies. Zoological Research,33(6): 566-573.[in Chinese])

    胡正堅(jiān), 陳雙林. 1992. 竹筍夜蛾防治試驗(yàn)初報(bào). 竹子研究匯刊 ,11(3): 37-41.

    (Hu Z J,Chen S L.1992.A control trial on bamboo shoot MothOligiavulgaris(butler). Journal of Bamboo Reseacher,11(3): 37-41.[in Chinese])

    黃煥華, 邵識(shí)煩, 陳 拓, 等. 2003.茶稈竹竹筍夜蛾綜合防治技術(shù)研究. 林業(yè)科學(xué)研究, 16(4): 444-448.

    (Huang H H, Shao S F, Chen T,etal.2003.Study on the integrated control technique of the noctuid on the bamboo shoot ofPseudosasaamabilis. Forest Research, 16(4): 444-448. [in Chinese])

    黃瓊瑤, 舒金平, 張愛(ài)良, 等. 2009. 3 種竹筍夜蛾生態(tài)位及其種間競(jìng)爭(zhēng)的研究. 林業(yè)科學(xué)研究, 22(5): 647-651.

    (Huang Q Y, Shu J P, Zhang A L,etal. 2009.Niche of three dominant shoot-boring noctuids and their interspecific competition. Forest Research, 22(5): 647-651. [in Chinese])

    李青青, 李地艷, 段焰青, 等. 2010. DNA條形碼在鱗翅目昆蟲(chóng)中的應(yīng)用. 生命科學(xué), 22(4): 307-312.

    (Li Q Q, Li D Y, Duan Y Q,etal. 2010. Application of DNA barcoding in lepidopteran insects.Chinese Bulletin of Life Sciences, 22(4): 307-312. [in Chinese])

    梁光紅, 林毓銀. 2003. 黃甜竹基夜蛾生物學(xué)特性及其防治. 福建農(nóng)林大學(xué)學(xué)報(bào):自然科學(xué)版, 32(1):36-40.

    (Liang G H, Lin Y Y. 2003. Study on the biological characteristics and control ofKumasiakumaso(Sugi) onAcidosasaedulis. Journal of Fujian Agriculture and Forestry University: Natural Science Edition, 32(1):36-40. [in Chinese])

    羅梅浩, 郭線茹, 鄭曉軍, 等. 2002. 煙青蟲(chóng)和棉鈴蟲(chóng)在煙草上的生態(tài)位及其種間競(jìng)爭(zhēng). 中國(guó)煙草學(xué)報(bào), 8(4): 34-37.

    (Luo M H, Guo X R, Zheng X J,etal. 2002.The niche and interspecific competition of oriental tobacco budworm and cotton bollworm in tobacco plants. Acta Tabacaria Sinica, 8(4): 34-37. [in Chinese])

    舒金平, 滕 瑩, 張愛(ài)良, 等. 2012. 竹筍基夜蛾的求偶及交配行為. 應(yīng)用生態(tài)學(xué)報(bào), 23(12): 3421-3428.

    (Shu J P,Teng Y,Zhang A L,etal. 2012.Calling and mating behaviors of bamboo shoot borerKumasiakumaso. Chinese Journal of Applied Ecology, 23(12): 3421-3428. [in Chinese])

    徐天森. 1962.竹筍夜蛾防治研究. 林業(yè)科學(xué), 7(4): 283-291.

    (Xu TS.1962.A control trial on bamboo shoot Motholigiavulgaris(butler). Scientia Silvae Sinicae, 7(4): 283-291. [in Chinese])

    徐天森, 王浩杰. 2004. 中國(guó)竹子主要害蟲(chóng). 北京: 中國(guó)林業(yè)出版社, 45-52.

    (Xu T S. Wang H J. 2004.Main pests of bamboo in China. Beijing: China Forestry Publishing House. [in Chinese])

    肖金花,肖 暉,黃大衛(wèi).2004. 生物分類(lèi)學(xué)的新動(dòng)向—DNA條形編碼.動(dòng)物學(xué)報(bào),50(5): 852-855.

    (Xiao J H, Xiao H, Huang D W. 2004. DNA barcoding:new approach of biological taxonomy. Acta Zoologica Sinica,50(5): 852-855. [in Chinese])

    楊瑞生, 王振營(yíng), 何康來(lái), 等. 2008. 稈野螟屬部分種的線粒體COⅡ基因序列分析及其分子系統(tǒng)學(xué) (鱗翅目:草螟科). 昆蟲(chóng)學(xué)報(bào), 51(2): 182-189.

    (Yang R S, Wang Z Y, He K L,etal. 2008. Sequence analysis of mtDNA-COⅡ gene and molecular phylogeny ofOstriniaspp. from China (Lepidoptera: Crambidae).Acta Entomologica Sinica, 51(2): 182-189.[in Chinese])

    楊瑞生, 王振營(yíng), 何康來(lái), 等. 2011. 中國(guó)稈野螟屬昆蟲(chóng)線粒體CO I玉基因遺傳多樣性及其分子系統(tǒng)學(xué)研究 (鱗翅目:草螟科). 南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 34(5):73-80.

    (Yang R S, Wang Z Y, He K L,etal. 2011. Genetic diversity and phylogeny of the genusOstrinia(Lepidoptera:Crambidae) inhabiting China inferred from mitochondrial COⅠ gene. Journal of Nanjing Agricultural University, 34(5):73-80. [in Chinese])

    楊瑞生, 鐘 亮, 姜義仁, 等. 2012. 橡實(shí)象蟲(chóng)等25種昆蟲(chóng)線粒體COⅠ基因的遺傳多樣性及系統(tǒng)進(jìn)化分析. 蠶業(yè)科學(xué), 37(6): 985-992.

    (Yang R S,Zhong L,Jiang Y R,etal. 2012. Genetic diversity and phylogenetic analysis of COⅠ genes fromCurculioarakawaiWeevil and 24 other insect species. Science of Sericulture, 37(6): 985-992.[in Chinese])

    楊瑞生, 姜義仁, 石生林, 等. 2013. 櫟黃掌舟蛾線粒體 DNA COⅠ基因的堿基組成及系統(tǒng)進(jìn)化分析. 蠶業(yè)科學(xué), 39(3):433-441.

    (Yang R S, Jiang Y R, Shi S L,etal. 2013. Base composition and phylogenetic analysis of mitochondrial DNA COⅠ gene fromPhaleraassimilis. Science of Sericulture, 39(3):433-441. [in Chinese])

    葉碧歡, 張亞波, 滕 瑩,等. 2014. 筍秀夜蛾的求偶及交配行為. 生態(tài)學(xué)雜志, 33(8): 2136-2141.

    (Ye B H, Zhang Y B, Teng Y,etal. 2014. Courtship and mating behavior of bamboo shoot borer,Apameaapameoides(Draudt): (Lepidoptera: Noctuidae). Chinese Journal of Ecology,33(8): 2136-2141.[in Chinese])

    周繼亮, 張亞平, 黃美華,等. 2001. 蝮亞科蛇線粒體細(xì)胞色素b基因序列分析及其系統(tǒng)發(fā)育. 動(dòng)物學(xué)報(bào), 47(4): 361-366.

    (Zhou J L, Zhang Y P, Huang M H,etal. 2001. Phylogenetic relationships among crotalinae based on mitochondrial cytochrome b gene sequence variations. Acta Zoologica Sinica, 47(4): 361-366. [in Chinese])

    朱弘復(fù). 1964.中國(guó)經(jīng)濟(jì)昆蟲(chóng)志.第6冊(cè).鱗翅目夜蛾科.北京:科學(xué)出版社,1-206.

    (Zhu H F.1964. Economic insect fauna of China. Chapter 6. Lepidoptera: Noctuidae. Beijing:Science Press.[in Chinese])

    Boore J L. 1999. Animal mitochondrial genomes. Nucleic Acids Research, 27(8): 1767-1780.

    Folmer O,Black M, Hoeh W,etal. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit Ⅰ from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology, 3(5): 294-299.

    Foster S P, Harris M O. 1997. Behavioral manipulation methods for insect pest-management. Annual Review of Entomology, 42(1): 123-146.

    Hebert P D N, Cywinska A, Ball S L.2003. Biological identifications through DNA barcodes. Proceedings of the Royal Society of London B: Biological Sciences, 270(1512): 313-321.

    Inoue H,Sugi S,Kuroko H,etal. 1982. Moths of Japan,Vol.Ⅰ-Ⅱ.Tokyo:Kodansha, 1:746.

    Lunt D H, Zhang D X, Szymura J M,etal. 1996. The insect cytochrome oxidase Ⅰ gene: evolutionary patterns and conserved primers for phylogenetic studies. Insect Molecular Biology, 5(3): 153-165.

    Monaghan M T, Balke M, Gregory T R,etal. 2005. DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers. Philosophical Transactions of the Royal Society B: Biological Sciences, 360(1462): 1925-1933.

    Nelson L A, Wallman J F, Dowton M. 2007. Using COⅠ barcodes to identify forensically and medically important blowflies. Medical and Veterinary Entomology, 21(1): 44-52.

    Taanman J W. 1999. The mitochondrial genome: structure, transcription, translation and replication. Biochimica et Biophysica Acta (BBA)-Bioenergetics, 1410(2): 103-123.

    Tautz D, Arctander P, Minelli A,etal. 2003. A plea for DNA taxonomy.Trends in Ecology and Evolution, 18: 70-74

    Yoshimatsu S, Kusigemati K, Gyoutoku N,etal. 2005. Some lepidopterous pests of bamboo and bamboograss shoots in Japan. Japanese Journal of Entomology (New Series), 8(3): 91-97.

    (責(zé)任編輯 朱乾坤)

    Genetic Relationship of Shoot-Boring Noctuids Based on Mitochondrial Cytochrome Oxidase Ⅰ Gene

    Zhang Qin1,2Mo Youdi3Zhang Yabo1Shu Jinping1Wang Haojie1Wu Hong1

    (1.ResearchInstituteofSubtropicalForestry,CAFFuyang311400; 2.HangzhouAcademyofForestryHangzhou311300; 3.HuaijiForestryBureauofZhaoqing,GuangdongProvinceZhaoging526000)

    【Objective】 The bamboo shoot-boring noctuid is a joint name for several species of moths which are hamful to bamboo shoot stage in southern China. There are four dominant species that have sympatric distribution and the same host. Their biological and ecological characteristics are similar. However,the classification status and evolutionary approach are still controversial. To determine their classification status and provide a molecular theory foundation for the system evolution, their genetic relationship was studied in the paper. 【Method】The bamboo shoot borers, distributing in the bamboo production areas surrounding in Hangzhou, Zhejiang Province, were used as the research object, and four 709 bp segments of cytochrome C oxidase subunit Ⅰ(COⅠ) gene were cloned, using universal primers, from mitochondrial DNA of different bamboo shoot borers′ individuals and compared with homologous sequences from other lepidopterous insects to analyze their base composition characteristics and the genetic relationship for accurate identification of bamboo shoot borers at molecular level by building system evolutionary trees.【Result】1)The results indicated that the genetic distance between the four kinds of bamboo shoot borers was small.2)Three phylogenetic trees showed that they got 100% support on the branch.【Conclusion】 The four bamboo shoot borers have closed genetic relationship,among which three species can be classified into one genus exceptSapporiarepetita.Furthermore, the four kinds of bamboo shoot borers could be used as ideal samples for sympatric speciation research for its similarity in body and habits and distribution.This study provides preliminary molecular basis for the genetic relationship and DNA barcoding technology of the 4 kinds of bamboo shoot moths at the same time.

    shoot-boring noctuids; mitochondrial cytochrome oxidaseⅠgene; sequence comparison; genetic relationship

    10.11707/j.1001-7488.20170411

    2015-05-14;

    2016-12-07。

    林業(yè)公益性行業(yè)科研專(zhuān)項(xiàng)(201304403)。

    S718.7

    A

    1001-7488(2017)04-0096-09

    * 舒金平為通訊作者。

    猜你喜歡
    親緣進(jìn)化樹(shù)堿基
    谷子近緣野生種的親緣關(guān)系及其利用研究
    中國(guó)醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    基于心理旋轉(zhuǎn)的小學(xué)生物進(jìn)化樹(shù)教學(xué)實(shí)驗(yàn)報(bào)告
    常見(jiàn)的進(jìn)化樹(shù)錯(cuò)誤概念及其辨析*
    應(yīng)用思維進(jìn)階構(gòu)建模型 例談培養(yǎng)學(xué)生創(chuàng)造性思維
    中國(guó)科學(xué)家創(chuàng)建出新型糖基化酶堿基編輯器
    生命“字母表”迎來(lái)4名新成員
    生命“字母表”迎來(lái)4名新成員
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    小白菜種質(zhì)遺傳多樣性與親緣關(guān)系的SRAP 和SSR分析
    国产男靠女视频免费网站| 99热6这里只有精品| 一本一本综合久久| 香蕉国产在线看| 狠狠狠狠99中文字幕| 久久久国产精品麻豆| 日本精品一区二区三区蜜桃| 亚洲精品色激情综合| 国产成人福利小说| 九色成人免费人妻av| 国产一区二区在线av高清观看| 观看美女的网站| 18美女黄网站色大片免费观看| 午夜精品在线福利| 久久久久性生活片| 老司机午夜福利在线观看视频| 中出人妻视频一区二区| 亚洲国产精品久久男人天堂| 99国产综合亚洲精品| 51午夜福利影视在线观看| 欧美在线一区亚洲| 欧美黑人巨大hd| 亚洲五月婷婷丁香| 午夜久久久久精精品| 亚洲av五月六月丁香网| 99re在线观看精品视频| а√天堂www在线а√下载| 99久久精品国产亚洲精品| 日日夜夜操网爽| 亚洲电影在线观看av| 黄色女人牲交| 免费搜索国产男女视频| 人妻久久中文字幕网| 欧美日本视频| 嫩草影视91久久| 哪里可以看免费的av片| 国产淫片久久久久久久久 | 国产精品自产拍在线观看55亚洲| 国产麻豆成人av免费视频| 中文字幕久久专区| 男人舔奶头视频| 亚洲国产欧美人成| 狂野欧美白嫩少妇大欣赏| 免费在线观看日本一区| av在线蜜桃| 岛国视频午夜一区免费看| 这个男人来自地球电影免费观看| 夜夜躁狠狠躁天天躁| 亚洲精品在线美女| 波多野结衣高清作品| 日本一二三区视频观看| 欧美一级a爱片免费观看看| 在线a可以看的网站| 热99在线观看视频| 日本 欧美在线| 欧美中文综合在线视频| 全区人妻精品视频| 久久久水蜜桃国产精品网| 男人和女人高潮做爰伦理| 最近最新中文字幕大全电影3| av福利片在线观看| 国产成年人精品一区二区| 日韩高清综合在线| 亚洲精品乱码久久久v下载方式 | 亚洲激情在线av| 久久精品91无色码中文字幕| 亚洲av成人精品一区久久| 夜夜看夜夜爽夜夜摸| 国产91精品成人一区二区三区| 日韩高清综合在线| 午夜精品一区二区三区免费看| 身体一侧抽搐| 国产一区二区三区在线臀色熟女| 校园春色视频在线观看| 国产精品一区二区三区四区免费观看 | 最新中文字幕久久久久 | 欧美乱码精品一区二区三区| 亚洲国产精品999在线| 亚洲欧美一区二区三区黑人| 精品久久蜜臀av无| 变态另类成人亚洲欧美熟女| 日本精品一区二区三区蜜桃| 亚洲成a人片在线一区二区| 韩国av一区二区三区四区| 最好的美女福利视频网| 国产高清有码在线观看视频| 日韩欧美国产在线观看| 我的老师免费观看完整版| 亚洲av免费在线观看| tocl精华| 国产成人精品久久二区二区免费| 免费看a级黄色片| 欧美性猛交╳xxx乱大交人| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品av在线| 国产三级在线视频| 91久久精品国产一区二区成人 | 日日夜夜操网爽| 久久久久久人人人人人| 香蕉丝袜av| 两人在一起打扑克的视频| 国产精品亚洲美女久久久| 少妇人妻一区二区三区视频| 18美女黄网站色大片免费观看| 九色国产91popny在线| 999久久久精品免费观看国产| 十八禁网站免费在线| 男人和女人高潮做爰伦理| 一级毛片精品| 亚洲专区中文字幕在线| 偷拍熟女少妇极品色| 老司机午夜福利在线观看视频| 曰老女人黄片| 99在线人妻在线中文字幕| 后天国语完整版免费观看| 中文字幕av在线有码专区| 欧美日本亚洲视频在线播放| 日本a在线网址| 欧美极品一区二区三区四区| 黄色成人免费大全| 哪里可以看免费的av片| www国产在线视频色| 成人国产综合亚洲| 在线永久观看黄色视频| 岛国在线免费视频观看| 欧美成狂野欧美在线观看| 亚洲人成伊人成综合网2020| 女人被狂操c到高潮| xxx96com| 国产精品永久免费网站| 欧美乱色亚洲激情| 国产三级中文精品| 宅男免费午夜| 91久久精品国产一区二区成人 | 嫩草影院入口| 色综合欧美亚洲国产小说| 麻豆一二三区av精品| 亚洲国产高清在线一区二区三| 美女黄网站色视频| 国产真实乱freesex| 一级毛片女人18水好多| 国产精品综合久久久久久久免费| 制服丝袜大香蕉在线| 国产精品久久电影中文字幕| 久久久久久人人人人人| 我要搜黄色片| 中文在线观看免费www的网站| 麻豆久久精品国产亚洲av| 久久久久久大精品| 一卡2卡三卡四卡精品乱码亚洲| 黑人巨大精品欧美一区二区mp4| 哪里可以看免费的av片| 亚洲欧美精品综合久久99| 变态另类成人亚洲欧美熟女| 亚洲精品国产精品久久久不卡| xxxwww97欧美| 久久香蕉精品热| 国产三级黄色录像| 日本三级黄在线观看| 精品一区二区三区视频在线观看免费| 免费看美女性在线毛片视频| 观看免费一级毛片| 18禁观看日本| 亚洲国产欧美一区二区综合| 免费观看精品视频网站| www日本黄色视频网| 十八禁人妻一区二区| 波多野结衣巨乳人妻| 成人三级黄色视频| 好男人在线观看高清免费视频| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 久久久国产成人免费| 熟妇人妻久久中文字幕3abv| 五月玫瑰六月丁香| 亚洲avbb在线观看| 91字幕亚洲| 亚洲精品中文字幕一二三四区| 午夜久久久久精精品| 男人舔奶头视频| 亚洲国产精品成人综合色| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 亚洲黑人精品在线| 久久久久久久久中文| 午夜福利欧美成人| 午夜激情欧美在线| 欧美不卡视频在线免费观看| 国产欧美日韩精品亚洲av| 亚洲精品456在线播放app | av片东京热男人的天堂| 亚洲欧美精品综合一区二区三区| 色综合站精品国产| xxx96com| 热99在线观看视频| 国产视频内射| 日韩欧美 国产精品| 日韩欧美免费精品| 男女那种视频在线观看| av欧美777| 精品一区二区三区视频在线 | 啦啦啦免费观看视频1| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲美女久久久| 亚洲,欧美精品.| 最新在线观看一区二区三区| 我的老师免费观看完整版| 两个人视频免费观看高清| 久久久久久久久免费视频了| 丁香欧美五月| 在线观看免费午夜福利视频| 老司机在亚洲福利影院| av在线天堂中文字幕| 欧美国产日韩亚洲一区| 在线视频色国产色| 午夜激情福利司机影院| 在线观看66精品国产| 久久久久久大精品| 国产高清videossex| 日韩欧美国产一区二区入口| 久久午夜亚洲精品久久| h日本视频在线播放| 九九热线精品视视频播放| 欧美中文综合在线视频| 少妇人妻一区二区三区视频| 一进一出抽搐动态| 亚洲国产高清在线一区二区三| 免费大片18禁| 丁香六月欧美| 丰满人妻一区二区三区视频av | 中文字幕高清在线视频| 少妇丰满av| 精品熟女少妇八av免费久了| 国产精品98久久久久久宅男小说| 嫩草影院精品99| 午夜久久久久精精品| 午夜福利18| 亚洲精品久久国产高清桃花| 亚洲自偷自拍图片 自拍| 午夜免费观看网址| 老司机福利观看| 国产乱人视频| 国产精品久久视频播放| 亚洲最大成人中文| 久久精品综合一区二区三区| 久久伊人香网站| 国产高潮美女av| 色播亚洲综合网| 中文在线观看免费www的网站| 亚洲国产精品sss在线观看| 男人和女人高潮做爰伦理| 岛国在线观看网站| 观看美女的网站| 国产欧美日韩精品一区二区| 美女黄网站色视频| 精品午夜福利视频在线观看一区| av天堂在线播放| a在线观看视频网站| 91在线观看av| 久久久久久大精品| 丁香六月欧美| 国产成年人精品一区二区| 欧美日韩福利视频一区二区| 好看av亚洲va欧美ⅴa在| 岛国视频午夜一区免费看| 午夜精品久久久久久毛片777| 亚洲av成人一区二区三| 亚洲美女视频黄频| 性欧美人与动物交配| 免费在线观看影片大全网站| 91av网站免费观看| 午夜福利成人在线免费观看| 一级毛片精品| 国产一级毛片七仙女欲春2| 日韩欧美一区二区三区在线观看| 国产精华一区二区三区| 免费在线观看日本一区| 人人妻,人人澡人人爽秒播| 在线观看美女被高潮喷水网站 | 亚洲专区中文字幕在线| 亚洲性夜色夜夜综合| 性欧美人与动物交配| 黄色成人免费大全| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 黄色女人牲交| 亚洲人成电影免费在线| 国产欧美日韩精品一区二区| 国产91精品成人一区二区三区| 麻豆av在线久日| 亚洲一区二区三区不卡视频| 噜噜噜噜噜久久久久久91| 91老司机精品| 深夜精品福利| 亚洲精品乱码久久久v下载方式 | 日日干狠狠操夜夜爽| 亚洲av中文字字幕乱码综合| 亚洲精品美女久久av网站| 丰满人妻熟妇乱又伦精品不卡| 天天躁狠狠躁夜夜躁狠狠躁| 久久这里只有精品19| 日本黄大片高清| 性色avwww在线观看| 欧美不卡视频在线免费观看| 91久久精品国产一区二区成人 | 中文亚洲av片在线观看爽| 亚洲国产精品成人综合色| 亚洲自偷自拍图片 自拍| 亚洲狠狠婷婷综合久久图片| 一区二区三区国产精品乱码| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 一进一出好大好爽视频| 蜜桃久久精品国产亚洲av| 午夜免费成人在线视频| 老司机午夜十八禁免费视频| 嫩草影视91久久| 成人av在线播放网站| 夜夜看夜夜爽夜夜摸| 亚洲精华国产精华精| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 啦啦啦观看免费观看视频高清| 色吧在线观看| 久久人妻av系列| 精品国产三级普通话版| 色尼玛亚洲综合影院| 国产又色又爽无遮挡免费看| 精品99又大又爽又粗少妇毛片 | 国产精品一区二区免费欧美| 美女黄网站色视频| 国产在线精品亚洲第一网站| 波多野结衣高清作品| 国产高清三级在线| 日韩欧美三级三区| 丰满人妻熟妇乱又伦精品不卡| av欧美777| 精品久久蜜臀av无| 在线观看一区二区三区| 欧美另类亚洲清纯唯美| av欧美777| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩无卡精品| 成人三级做爰电影| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 亚洲成人中文字幕在线播放| 国产激情久久老熟女| 午夜福利欧美成人| 国产久久久一区二区三区| 国产精品香港三级国产av潘金莲| 成人三级做爰电影| 婷婷六月久久综合丁香| 丰满人妻熟妇乱又伦精品不卡| 草草在线视频免费看| 欧美日韩亚洲国产一区二区在线观看| 久久这里只有精品中国| 高清在线国产一区| 国产精品久久久久久精品电影| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频| 国产精品影院久久| 国产视频内射| 亚洲成人精品中文字幕电影| 亚洲av电影不卡..在线观看| 久久久成人免费电影| 国产精品久久电影中文字幕| 午夜免费观看网址| 午夜精品一区二区三区免费看| 天堂动漫精品| 美女被艹到高潮喷水动态| 国产午夜精品久久久久久| 欧美一级a爱片免费观看看| 亚洲熟妇中文字幕五十中出| 99热精品在线国产| 精品久久久久久久毛片微露脸| 99久久99久久久精品蜜桃| 国产精品自产拍在线观看55亚洲| 熟妇人妻久久中文字幕3abv| 色播亚洲综合网| 一个人看的www免费观看视频| 精品久久久久久久末码| 亚洲熟女毛片儿| 亚洲精品久久国产高清桃花| 人妻夜夜爽99麻豆av| 精品熟女少妇八av免费久了| 精品人妻1区二区| 欧美黄色淫秽网站| 午夜a级毛片| 欧美国产日韩亚洲一区| 岛国视频午夜一区免费看| 午夜精品一区二区三区免费看| 天堂√8在线中文| 欧美最黄视频在线播放免费| 成人特级av手机在线观看| 欧美最黄视频在线播放免费| 国产成人福利小说| 国产熟女xx| 偷拍熟女少妇极品色| 日本熟妇午夜| 免费看十八禁软件| 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 久久久久久人人人人人| 夜夜夜夜夜久久久久| 亚洲美女黄片视频| 成熟少妇高潮喷水视频| 男女午夜视频在线观看| 精品久久久久久成人av| 香蕉av资源在线| 久久国产精品人妻蜜桃| 哪里可以看免费的av片| 亚洲无线观看免费| 99久久无色码亚洲精品果冻| 国产亚洲av高清不卡| 久久午夜亚洲精品久久| 国产精品98久久久久久宅男小说| 香蕉国产在线看| 18禁观看日本| 成年女人毛片免费观看观看9| 最新中文字幕久久久久 | 色在线成人网| 欧美另类亚洲清纯唯美| 亚洲av免费在线观看| 一级毛片精品| 午夜亚洲福利在线播放| 在线观看66精品国产| 欧美成狂野欧美在线观看| 日韩中文字幕欧美一区二区| bbb黄色大片| www.自偷自拍.com| 老司机午夜十八禁免费视频| 少妇丰满av| 久久午夜综合久久蜜桃| 亚洲中文av在线| 欧美日韩国产亚洲二区| 亚洲五月天丁香| 九九在线视频观看精品| 黄片大片在线免费观看| 国产精品一区二区三区四区免费观看 | 久9热在线精品视频| 大型黄色视频在线免费观看| 成人鲁丝片一二三区免费| 久久精品夜夜夜夜夜久久蜜豆| 久久久久精品国产欧美久久久| 97超视频在线观看视频| 看片在线看免费视频| 嫩草影院精品99| 99国产精品一区二区三区| 狠狠狠狠99中文字幕| 国产熟女xx| 噜噜噜噜噜久久久久久91| 欧美大码av| 国产真人三级小视频在线观看| 国产成人av教育| 国产一区二区三区视频了| 亚洲五月天丁香| 不卡av一区二区三区| 国内毛片毛片毛片毛片毛片| 精品国产乱子伦一区二区三区| 男女做爰动态图高潮gif福利片| 亚洲精品456在线播放app | 午夜福利视频1000在线观看| 欧美精品啪啪一区二区三区| 十八禁人妻一区二区| 日韩欧美 国产精品| 国产97色在线日韩免费| 色综合亚洲欧美另类图片| 啦啦啦韩国在线观看视频| 99国产综合亚洲精品| h日本视频在线播放| 亚洲欧美日韩高清专用| 女人高潮潮喷娇喘18禁视频| 午夜福利成人在线免费观看| 国产高清视频在线观看网站| 国产伦在线观看视频一区| 亚洲av熟女| av天堂中文字幕网| 男人的好看免费观看在线视频| 伊人久久大香线蕉亚洲五| 又爽又黄无遮挡网站| 欧美午夜高清在线| 天堂av国产一区二区熟女人妻| 18禁美女被吸乳视频| 一级黄色大片毛片| 成人性生交大片免费视频hd| 国产成人系列免费观看| 久久久久久久久免费视频了| 久久中文字幕一级| 变态另类成人亚洲欧美熟女| 热99re8久久精品国产| 亚洲成人中文字幕在线播放| 国产成人av激情在线播放| 国产一区二区三区在线臀色熟女| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区精品| 久久人人精品亚洲av| 午夜成年电影在线免费观看| 亚洲精品国产精品久久久不卡| 超碰成人久久| a在线观看视频网站| 日韩欧美国产一区二区入口| 日本在线视频免费播放| 欧美一区二区国产精品久久精品| 成年版毛片免费区| 国产精品美女特级片免费视频播放器 | 精品欧美国产一区二区三| 国产精品久久久久久亚洲av鲁大| 色尼玛亚洲综合影院| 亚洲精品美女久久av网站| 一个人看视频在线观看www免费 | 午夜精品一区二区三区免费看| 国产精品女同一区二区软件 | 精品欧美国产一区二区三| 国产男靠女视频免费网站| 成人18禁在线播放| 淫秽高清视频在线观看| 三级男女做爰猛烈吃奶摸视频| 操出白浆在线播放| 丁香欧美五月| 中文字幕久久专区| 两个人视频免费观看高清| 麻豆久久精品国产亚洲av| 久久久久国产一级毛片高清牌| 美女高潮喷水抽搐中文字幕| 亚洲中文日韩欧美视频| 18美女黄网站色大片免费观看| www.自偷自拍.com| 国产av一区在线观看免费| 国产精品日韩av在线免费观看| 国产伦精品一区二区三区四那| 亚洲国产精品成人综合色| 日本免费一区二区三区高清不卡| 日韩欧美在线二视频| 欧美性猛交黑人性爽| 亚洲色图av天堂| 国产麻豆成人av免费视频| 男女做爰动态图高潮gif福利片| 一区二区三区国产精品乱码| 麻豆av在线久日| 国产精品久久久久久亚洲av鲁大| 日韩欧美国产在线观看| 香蕉丝袜av| 欧美日韩黄片免| 女人被狂操c到高潮| 国产精品免费一区二区三区在线| 欧美乱码精品一区二区三区| 久久久久国内视频| 精品人妻1区二区| 最近最新中文字幕大全电影3| 精品国产乱码久久久久久男人| 国产极品精品免费视频能看的| 老司机午夜福利在线观看视频| 长腿黑丝高跟| 欧美日韩精品网址| 亚洲色图 男人天堂 中文字幕| 国产一区二区在线观看日韩 | 亚洲自偷自拍图片 自拍| tocl精华| 久久精品夜夜夜夜夜久久蜜豆| 少妇的逼水好多| 一个人观看的视频www高清免费观看 | 亚洲精品美女久久av网站| 曰老女人黄片| 国内精品久久久久久久电影| 色综合亚洲欧美另类图片| 久久久久免费精品人妻一区二区| 久久午夜亚洲精品久久| 国产又色又爽无遮挡免费看| 成人精品一区二区免费| 国产精品久久视频播放| 2021天堂中文幕一二区在线观| 久久人妻av系列| 狠狠狠狠99中文字幕| 精品电影一区二区在线| svipshipincom国产片| 精品久久久久久久毛片微露脸| 日韩中文字幕欧美一区二区| 99久久综合精品五月天人人| 久久久久久久午夜电影| 亚洲熟女毛片儿| 亚洲欧美精品综合久久99| 琪琪午夜伦伦电影理论片6080| 在线观看免费视频日本深夜| www国产在线视频色| 亚洲精品国产精品久久久不卡| 又紧又爽又黄一区二区| 国产单亲对白刺激| 久久精品影院6| 亚洲专区国产一区二区| 欧美3d第一页| 久久香蕉国产精品| 黄色日韩在线| 欧美中文日本在线观看视频| 老汉色av国产亚洲站长工具| 欧美不卡视频在线免费观看| xxxwww97欧美| 哪里可以看免费的av片| 一级黄色大片毛片| 成年女人毛片免费观看观看9| 亚洲国产欧美人成| 亚洲成人免费电影在线观看| 色av中文字幕| 国产精品自产拍在线观看55亚洲| 亚洲一区二区三区色噜噜| 国产精品 欧美亚洲| 久久久精品欧美日韩精品| 国内精品美女久久久久久| 日本在线视频免费播放| 日本免费a在线| 国产乱人视频| 99热这里只有是精品50| 亚洲国产欧美一区二区综合|