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

    基于線粒體COⅠ基因的竹筍夜蛾親緣關系*

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

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

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

    ?

    基于線粒體COⅠ基因的竹筍夜蛾親緣關系*

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

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

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

    竹筍夜蛾; COⅠ基因; 序列比較; 親緣關系

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

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

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

    1 材料與方法

    1.1 供試蟲源

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

    1.2 試驗材料

    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)細胞,其中擴增4種竹筍夜蛾的線粒體COⅠ基因所用引物來自通用引物(Folmeretal., 1994):

    L1490-5’-GGTCAACAAATCATAAAGATATT G-3’;

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

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

    1.3 試驗方法

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

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

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

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

    表1 21種蛾類昆蟲COⅠ 基因信息①Tab.1 Information of COⅠ genes of 21 noctuid insects

    ① *:無中文屬名No Chinese name.

    2 結果與分析

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

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

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

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

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

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

    2.3 4種竹筍夜蛾及其比對昆蟲mtDNA COⅠ基因片段的序列特征

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

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

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

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

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

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

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

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

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

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

    用4種竹筍夜蛾及其比對昆蟲的mtDNA COⅠ氨基酸序列在Mega5.0中構建進化樹,結果如圖3、4、5所示,3種進化樹顯示結果大致相同。與比對序列相比,4種竹筍夜蛾分屬同一支,分支支持率達100%,親緣關系最近,其中薩夜蛾與其他3種夜蛾分屬2個分支,分支支持率為70%左右。MP法進化樹中,基夜蛾又與竹筍禾夜蛾和筍秀夜蛾分為2支,分支支持率為30%。3種進化樹均顯示,4種竹筍夜蛾與禾夜蛾屬、秀夜蛾屬的昆蟲遺傳距離較遠。

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

    3 討論

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

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

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

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

    4 結論

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

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

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

    陳一心. 1992. 鱗翅目: 夜蛾.長沙: 湖南科學技術出版社, 947-1040.

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

    陳 星, 沈永義, 張亞平. 2012. 線粒體DNA在分子進化研究中的應用. 動物學研究, 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])

    胡正堅, 陳雙林. 1992. 竹筍夜蛾防治試驗初報. 竹子研究匯刊 ,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])

    黃煥華, 邵識煩, 陳 拓, 等. 2003.茶稈竹竹筍夜蛾綜合防治技術研究. 林業(yè)科學研究, 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])

    黃瓊瑤, 舒金平, 張愛良, 等. 2009. 3 種竹筍夜蛾生態(tài)位及其種間競爭的研究. 林業(yè)科學研究, 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條形碼在鱗翅目昆蟲中的應用. 生命科學, 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. 黃甜竹基夜蛾生物學特性及其防治. 福建農(nóng)林大學學報:自然科學版, 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. 煙青蟲和棉鈴蟲在煙草上的生態(tài)位及其種間競爭. 中國煙草學報, 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])

    舒金平, 滕 瑩, 張愛良, 等. 2012. 竹筍基夜蛾的求偶及交配行為. 應用生態(tài)學報, 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è)科學, 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. 中國竹子主要害蟲. 北京: 中國林業(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. 生物分類學的新動向—DNA條形編碼.動物學報,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])

    楊瑞生, 王振營, 何康來, 等. 2008. 稈野螟屬部分種的線粒體COⅡ基因序列分析及其分子系統(tǒng)學 (鱗翅目:草螟科). 昆蟲學報, 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])

    楊瑞生, 王振營, 何康來, 等. 2011. 中國稈野螟屬昆蟲線粒體CO I玉基因遺傳多樣性及其分子系統(tǒng)學研究 (鱗翅目:草螟科). 南京農(nóng)業(yè)大學學報, 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. 橡實象蟲等25種昆蟲線粒體COⅠ基因的遺傳多樣性及系統(tǒng)進化分析. 蠶業(yè)科學, 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)進化分析. 蠶業(yè)科學, 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)學雜志, 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. 蝮亞科蛇線粒體細胞色素b基因序列分析及其系統(tǒng)發(fā)育. 動物學報, 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])

    朱弘復. 1964.中國經(jīng)濟昆蟲志.第6冊.鱗翅目夜蛾科.北京:科學出版社,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.

    (責任編輯 朱乾坤)

    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è)科研專項(201304403)。

    S718.7

    A

    1001-7488(2017)04-0096-09

    * 舒金平為通訊作者。

    猜你喜歡
    親緣進化樹堿基
    谷子近緣野生種的親緣關系及其利用研究
    作物學報(2022年2期)2022-11-06 12:11:02
    中國醫(yī)學科學院藥用植物研究所藥用植物親緣學研究中心
    基于心理旋轉的小學生物進化樹教學實驗報告
    常見的進化樹錯誤概念及其辨析*
    生物學通報(2021年9期)2021-07-01 03:24:44
    應用思維進階構建模型 例談培養(yǎng)學生創(chuàng)造性思維
    中國科學家創(chuàng)建出新型糖基化酶堿基編輯器
    生命“字母表”迎來4名新成員
    科學24小時(2019年5期)2019-06-11 08:39:38
    生命“字母表”迎來4名新成員
    菊科藥用植物遺傳多樣性及親緣關系的ISSR分析
    小白菜種質(zhì)遺傳多樣性與親緣關系的SRAP 和SSR分析
    午夜激情福利司机影院| av在线播放精品| 亚洲熟女精品中文字幕| 欧美xxⅹ黑人| 女人久久www免费人成看片| 成人亚洲精品一区在线观看| 亚洲第一av免费看| av在线app专区| 亚洲精品国产av蜜桃| 在线亚洲精品国产二区图片欧美 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产黄频视频在线观看| 边亲边吃奶的免费视频| 男女无遮挡免费网站观看| 热99久久久久精品小说推荐| 午夜视频国产福利| 热re99久久精品国产66热6| 三上悠亚av全集在线观看| 久久久久久久久久久丰满| 亚洲丝袜综合中文字幕| 亚洲国产日韩一区二区| 国产不卡av网站在线观看| 日韩精品免费视频一区二区三区 | 91久久精品国产一区二区成人| 美女cb高潮喷水在线观看| 久久午夜福利片| 色婷婷久久久亚洲欧美| 最近中文字幕高清免费大全6| 最近中文字幕高清免费大全6| 天堂8中文在线网| 免费久久久久久久精品成人欧美视频 | 热re99久久精品国产66热6| 欧美 亚洲 国产 日韩一| 久久久久精品性色| av国产精品久久久久影院| 成人无遮挡网站| 国产精品久久久久成人av| 免费观看的影片在线观看| av在线播放精品| 国内精品宾馆在线| 国产日韩欧美视频二区| 春色校园在线视频观看| 国内精品宾馆在线| 一级黄片播放器| 大香蕉97超碰在线| 啦啦啦啦在线视频资源| 人妻少妇偷人精品九色| 国模一区二区三区四区视频| 日韩 亚洲 欧美在线| 久久久久国产精品人妻一区二区| 午夜精品国产一区二区电影| 老司机影院成人| 99热6这里只有精品| 国产精品久久久久成人av| 黄色怎么调成土黄色| 天天操日日干夜夜撸| 久久狼人影院| 极品少妇高潮喷水抽搐| 高清在线视频一区二区三区| h视频一区二区三区| 18禁观看日本| 久久久久久久久久久丰满| 黑人巨大精品欧美一区二区蜜桃 | 亚洲四区av| 国产精品国产三级国产专区5o| 久久国产亚洲av麻豆专区| 国产成人精品婷婷| 日本猛色少妇xxxxx猛交久久| 少妇丰满av| 亚洲色图 男人天堂 中文字幕 | 国产不卡av网站在线观看| 国产一区有黄有色的免费视频| 国产在线一区二区三区精| 美女福利国产在线| 高清不卡的av网站| 亚洲精品av麻豆狂野| 欧美 亚洲 国产 日韩一| 少妇猛男粗大的猛烈进出视频| 亚洲av福利一区| 亚洲国产精品一区三区| 国产片内射在线| 日韩三级伦理在线观看| 欧美人与性动交α欧美精品济南到 | 美女大奶头黄色视频| 18在线观看网站| 九九在线视频观看精品| 在线观看www视频免费| 性色avwww在线观看| 99九九线精品视频在线观看视频| 一二三四中文在线观看免费高清| 毛片一级片免费看久久久久| 亚洲高清免费不卡视频| 热re99久久国产66热| 丰满少妇做爰视频| 三级国产精品片| 热re99久久精品国产66热6| 精品一区二区三区视频在线| 高清毛片免费看| 99热这里只有是精品在线观看| 国产在视频线精品| 在线观看www视频免费| 亚洲高清免费不卡视频| a级毛色黄片| 亚洲成人一二三区av| 国产黄色视频一区二区在线观看| 丁香六月天网| av视频免费观看在线观看| 大又大粗又爽又黄少妇毛片口| 久久鲁丝午夜福利片| 热99久久久久精品小说推荐| 大陆偷拍与自拍| 久久热精品热| 国产一区二区在线观看日韩| 极品人妻少妇av视频| 亚洲国产av影院在线观看| av天堂久久9| 久久影院123| 久久久久精品性色| 中文字幕av电影在线播放| 高清av免费在线| 黄色怎么调成土黄色| 国产视频首页在线观看| 久久韩国三级中文字幕| 少妇的逼好多水| 好男人视频免费观看在线| 国产精品嫩草影院av在线观看| 99热这里只有是精品在线观看| 欧美精品高潮呻吟av久久| 久久久久国产网址| 成人影院久久| 高清欧美精品videossex| 国产亚洲最大av| 如日韩欧美国产精品一区二区三区 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品久久久久久电影网| 99久久中文字幕三级久久日本| 好男人视频免费观看在线| 99热6这里只有精品| 亚洲少妇的诱惑av| 伊人久久精品亚洲午夜| 午夜日本视频在线| 高清av免费在线| 一级毛片我不卡| 欧美亚洲日本最大视频资源| 欧美精品高潮呻吟av久久| 看十八女毛片水多多多| 国产黄频视频在线观看| 大码成人一级视频| videossex国产| 久久精品国产亚洲网站| 欧美一级a爱片免费观看看| 看免费成人av毛片| 国产极品天堂在线| 777米奇影视久久| 久久午夜福利片| 18禁动态无遮挡网站| 日韩精品免费视频一区二区三区 | 尾随美女入室| 男女啪啪激烈高潮av片| 欧美精品国产亚洲| 搡老乐熟女国产| 国产精品一区www在线观看| 伊人久久国产一区二区| 精品亚洲成国产av| 女性生殖器流出的白浆| 久久热精品热| 在线观看三级黄色| 亚洲国产最新在线播放| 狂野欧美激情性bbbbbb| 日本欧美视频一区| 亚洲人成77777在线视频| 少妇高潮的动态图| 成人黄色视频免费在线看| 国产精品秋霞免费鲁丝片| 热re99久久国产66热| 欧美成人午夜免费资源| 亚洲第一av免费看| 天堂中文最新版在线下载| 国产极品粉嫩免费观看在线 | av黄色大香蕉| 亚洲精品美女久久av网站| 国产女主播在线喷水免费视频网站| 99热国产这里只有精品6| 久久99热这里只频精品6学生| 国产成人免费无遮挡视频| 亚洲不卡免费看| 国产一区二区三区综合在线观看 | 伊人亚洲综合成人网| 成人国产av品久久久| 看十八女毛片水多多多| 夜夜看夜夜爽夜夜摸| 日韩不卡一区二区三区视频在线| 丰满乱子伦码专区| 综合色丁香网| 欧美日韩国产mv在线观看视频| 高清在线视频一区二区三区| 色5月婷婷丁香| 国产成人精品在线电影| 少妇被粗大猛烈的视频| 亚洲国产精品专区欧美| 亚洲精品av麻豆狂野| 超碰97精品在线观看| 蜜桃在线观看..| 亚洲欧美成人精品一区二区| 18在线观看网站| 一级片'在线观看视频| 午夜久久久在线观看| 欧美xxxx性猛交bbbb| 欧美日韩av久久| 内地一区二区视频在线| 秋霞伦理黄片| 欧美老熟妇乱子伦牲交| 男女边摸边吃奶| 色吧在线观看| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久成人av| 亚洲国产色片| 免费少妇av软件| 久久av网站| 制服诱惑二区| 国产精品 国内视频| 黑人猛操日本美女一级片| 高清不卡的av网站| 婷婷成人精品国产| 黄色一级大片看看| 99精国产麻豆久久婷婷| 日韩中字成人| 日本猛色少妇xxxxx猛交久久| 美女cb高潮喷水在线观看| 午夜精品国产一区二区电影| 精品人妻在线不人妻| 国产黄色免费在线视频| 看十八女毛片水多多多| 夫妻性生交免费视频一级片| 人人妻人人澡人人看| 91精品伊人久久大香线蕉| 婷婷色av中文字幕| 亚洲国产精品999| 亚洲熟女精品中文字幕| 各种免费的搞黄视频| 一区二区三区乱码不卡18| tube8黄色片| 午夜福利网站1000一区二区三区| 日本黄大片高清| 99久久中文字幕三级久久日本| 亚洲精品av麻豆狂野| 亚洲精品一区蜜桃| 日韩精品免费视频一区二区三区 | 在线播放无遮挡| 国产白丝娇喘喷水9色精品| 欧美精品一区二区免费开放| 欧美xxⅹ黑人| 又粗又硬又长又爽又黄的视频| 美女xxoo啪啪120秒动态图| 晚上一个人看的免费电影| 亚洲无线观看免费| 日本与韩国留学比较| 人人澡人人妻人| 性色av一级| 十八禁网站网址无遮挡| 女性生殖器流出的白浆| 国产一区二区三区av在线| av卡一久久| 日韩欧美一区视频在线观看| 国产精品无大码| 亚洲成人手机| 精品亚洲乱码少妇综合久久| 色网站视频免费| 一级毛片 在线播放| 一本久久精品| videos熟女内射| 国产成人精品无人区| a级毛片免费高清观看在线播放| 欧美丝袜亚洲另类| 国产精品.久久久| 亚洲综合色惰| 大香蕉97超碰在线| 久久这里有精品视频免费| 婷婷色麻豆天堂久久| 久久综合国产亚洲精品| 18禁观看日本| 狠狠婷婷综合久久久久久88av| 日本猛色少妇xxxxx猛交久久| 国产亚洲一区二区精品| 国产免费现黄频在线看| 午夜激情久久久久久久| 国产视频内射| 亚洲中文av在线| 有码 亚洲区| 97在线人人人人妻| 国产欧美日韩综合在线一区二区| 国产69精品久久久久777片| 午夜福利在线观看免费完整高清在| 亚洲,一卡二卡三卡| 国产又色又爽无遮挡免| 久久久精品94久久精品| av专区在线播放| 爱豆传媒免费全集在线观看| 99热这里只有精品一区| 99热国产这里只有精品6| 久久99热6这里只有精品| 久久99一区二区三区| 一区二区三区免费毛片| 老司机影院毛片| 亚洲,一卡二卡三卡| 国产伦理片在线播放av一区| 亚洲av电影在线观看一区二区三区| 成人综合一区亚洲| 国产一区亚洲一区在线观看| 丁香六月天网| 中国三级夫妇交换| 十八禁高潮呻吟视频| av免费观看日本| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 秋霞在线观看毛片| 国产极品天堂在线| 日韩 亚洲 欧美在线| 亚洲国产精品一区二区三区在线| 日韩欧美精品免费久久| 大香蕉久久网| 91精品伊人久久大香线蕉| 国产在线视频一区二区| 搡老乐熟女国产| 久久这里有精品视频免费| 国产精品人妻久久久影院| 亚洲av在线观看美女高潮| 看非洲黑人一级黄片| 成人二区视频| 免费观看a级毛片全部| 亚洲五月色婷婷综合| av网站免费在线观看视频| 一本大道久久a久久精品| 国产男女超爽视频在线观看| 午夜老司机福利剧场| 免费大片黄手机在线观看| 观看美女的网站| 人妻一区二区av| av一本久久久久| 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 日本91视频免费播放| 精品一区二区免费观看| 新久久久久国产一级毛片| 日韩精品有码人妻一区| 久久久久久人妻| 亚洲精品日韩av片在线观看| 纯流量卡能插随身wifi吗| 91aial.com中文字幕在线观看| 亚洲欧美日韩另类电影网站| 久久99精品国语久久久| 韩国高清视频一区二区三区| 一区二区三区精品91| 狂野欧美白嫩少妇大欣赏| 22中文网久久字幕| 最近的中文字幕免费完整| 22中文网久久字幕| 99热全是精品| 国产一区二区三区av在线| 国产无遮挡羞羞视频在线观看| 日本与韩国留学比较| 建设人人有责人人尽责人人享有的| 久久久亚洲精品成人影院| 九色亚洲精品在线播放| 亚洲精品亚洲一区二区| 青春草国产在线视频| 亚洲av国产av综合av卡| 青春草国产在线视频| 91国产中文字幕| 久久精品久久久久久久性| 国产午夜精品久久久久久一区二区三区| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 中文字幕制服av| 久久精品久久久久久噜噜老黄| 欧美成人精品欧美一级黄| 欧美日韩在线观看h| 亚洲成人一二三区av| 国产 一区精品| 999精品在线视频| 亚洲精品日韩av片在线观看| 亚洲美女搞黄在线观看| 国产精品国产三级国产av玫瑰| 国产色爽女视频免费观看| 成人国产av品久久久| 欧美日韩综合久久久久久| 国产高清国产精品国产三级| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 日日摸夜夜添夜夜添av毛片| 免费观看性生交大片5| 免费大片黄手机在线观看| 亚洲成人一二三区av| 91国产中文字幕| 成人二区视频| 一区二区av电影网| 免费黄频网站在线观看国产| 2018国产大陆天天弄谢| 在线播放无遮挡| 人人澡人人妻人| 国产成人freesex在线| 亚洲精品av麻豆狂野| 亚洲成色77777| 人妻人人澡人人爽人人| 纵有疾风起免费观看全集完整版| 亚洲欧美一区二区三区国产| 91aial.com中文字幕在线观看| 久热久热在线精品观看| 亚洲国产av影院在线观看| av黄色大香蕉| 少妇的逼好多水| 欧美变态另类bdsm刘玥| 十分钟在线观看高清视频www| 另类精品久久| 日本免费在线观看一区| 黑人猛操日本美女一级片| 欧美3d第一页| 2022亚洲国产成人精品| 丝袜脚勾引网站| 狠狠婷婷综合久久久久久88av| av国产久精品久网站免费入址| 免费大片黄手机在线观看| 欧美性感艳星| 欧美变态另类bdsm刘玥| 国产亚洲欧美精品永久| 国产亚洲av片在线观看秒播厂| 3wmmmm亚洲av在线观看| 蜜桃国产av成人99| 街头女战士在线观看网站| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| 亚洲av免费高清在线观看| 99国产综合亚洲精品| 国产高清有码在线观看视频| 国产伦理片在线播放av一区| 国产精品欧美亚洲77777| 精品一区二区三卡| 国产精品人妻久久久影院| 久久久久久人妻| 日本色播在线视频| 天堂俺去俺来也www色官网| 免费高清在线观看视频在线观看| 九色亚洲精品在线播放| 亚洲精品aⅴ在线观看| 中文字幕av电影在线播放| 国产在线视频一区二区| 中国国产av一级| 亚洲精品一二三| 中国国产av一级| 尾随美女入室| 欧美性感艳星| 高清在线视频一区二区三区| 亚洲美女黄色视频免费看| 国产69精品久久久久777片| 欧美xxⅹ黑人| 亚洲av免费高清在线观看| 插逼视频在线观看| 日日爽夜夜爽网站| 日本黄色日本黄色录像| 亚洲综合色网址| 男女无遮挡免费网站观看| 爱豆传媒免费全集在线观看| 丝瓜视频免费看黄片| 欧美精品一区二区免费开放| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 亚洲五月色婷婷综合| 日韩av不卡免费在线播放| 国产成人aa在线观看| 国产女主播在线喷水免费视频网站| 一区二区三区乱码不卡18| 大香蕉久久网| 亚洲精品av麻豆狂野| 国产精品人妻久久久久久| 亚洲美女视频黄频| 国产精品久久久久久精品电影小说| 波野结衣二区三区在线| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 亚洲国产av新网站| 亚洲av不卡在线观看| 日韩av在线免费看完整版不卡| 国产片特级美女逼逼视频| 夜夜看夜夜爽夜夜摸| 91精品国产国语对白视频| 久久亚洲国产成人精品v| www.色视频.com| 久久精品久久久久久久性| 成年美女黄网站色视频大全免费 | 日本欧美视频一区| 亚洲怡红院男人天堂| 在线亚洲精品国产二区图片欧美 | 国产日韩欧美视频二区| 久久久久久久久久人人人人人人| 国产亚洲最大av| 久热这里只有精品99| 亚洲欧洲日产国产| 国产精品国产av在线观看| 各种免费的搞黄视频| 一个人免费看片子| 99热网站在线观看| av在线老鸭窝| av又黄又爽大尺度在线免费看| 男人添女人高潮全过程视频| 男女边吃奶边做爰视频| 欧美变态另类bdsm刘玥| 丰满乱子伦码专区| 日日摸夜夜添夜夜添av毛片| 黄色怎么调成土黄色| 久久精品国产亚洲av涩爱| 中文字幕av电影在线播放| 久久久久久久久久久免费av| 午夜福利影视在线免费观看| 亚洲在久久综合| 夫妻午夜视频| 成人综合一区亚洲| 国国产精品蜜臀av免费| 欧美老熟妇乱子伦牲交| 免费大片18禁| 不卡视频在线观看欧美| 视频中文字幕在线观看| 黄色一级大片看看| 国模一区二区三区四区视频| 国产男女超爽视频在线观看| 亚洲婷婷狠狠爱综合网| 午夜福利在线观看免费完整高清在| 亚洲国产精品成人久久小说| 午夜福利在线观看免费完整高清在| 我要看黄色一级片免费的| 我的老师免费观看完整版| 最近2019中文字幕mv第一页| 母亲3免费完整高清在线观看 | a级毛片在线看网站| 精品一区二区三区视频在线| 日韩一本色道免费dvd| 日本黄色片子视频| 精品久久久久久久久亚洲| 在线亚洲精品国产二区图片欧美 | 老熟女久久久| 中国三级夫妇交换| 国产淫语在线视频| 97在线视频观看| 久久久久久久亚洲中文字幕| 国产亚洲欧美精品永久| 午夜视频国产福利| 久久影院123| 久久av网站| 亚洲精品av麻豆狂野| 欧美最新免费一区二区三区| 人妻人人澡人人爽人人| 成年美女黄网站色视频大全免费 | 亚洲欧美精品自产自拍| 国产精品久久久久久精品电影小说| 欧美精品一区二区大全| 一边摸一边做爽爽视频免费| av福利片在线| 十八禁高潮呻吟视频| 国产一区亚洲一区在线观看| 免费黄频网站在线观看国产| 国国产精品蜜臀av免费| 天美传媒精品一区二区| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 亚洲av二区三区四区| 中文字幕精品免费在线观看视频 | 高清午夜精品一区二区三区| 亚洲四区av| 男女边摸边吃奶| 丝袜在线中文字幕| 精品久久久久久久久亚洲| 新久久久久国产一级毛片| 国产男人的电影天堂91| 精品人妻在线不人妻| a级毛片免费高清观看在线播放| 一区二区日韩欧美中文字幕 | 哪个播放器可以免费观看大片| 亚洲美女搞黄在线观看| 免费不卡的大黄色大毛片视频在线观看| 久久狼人影院| 在线观看国产h片| 国产在线免费精品| 国语对白做爰xxxⅹ性视频网站| 我的老师免费观看完整版| 久久久久国产精品人妻一区二区| 成人手机av| 国产深夜福利视频在线观看| 各种免费的搞黄视频| 欧美变态另类bdsm刘玥| 校园人妻丝袜中文字幕| 久久国产精品大桥未久av| a级毛片黄视频| 国产午夜精品久久久久久一区二区三区| 纯流量卡能插随身wifi吗| 日韩大片免费观看网站| 午夜视频国产福利| 亚洲精品久久午夜乱码| 99精国产麻豆久久婷婷| 欧美日韩av久久| 国产一区二区三区av在线| 亚洲av在线观看美女高潮| 国产国拍精品亚洲av在线观看| 日韩在线高清观看一区二区三区| 精品少妇内射三级| 免费高清在线观看视频在线观看| 欧美一级a爱片免费观看看| 91精品伊人久久大香线蕉| 免费人妻精品一区二区三区视频| 国产一区有黄有色的免费视频| 久久 成人 亚洲| 只有这里有精品99| 久久狼人影院| 一级毛片我不卡| 看非洲黑人一级黄片| 91精品一卡2卡3卡4卡| 国产老妇伦熟女老妇高清|