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

    基于線粒體COⅠ和COⅡ基因的沙果小食心蟲與梨小食心蟲的分子鑒定

    2016-06-15 03:14:02李玉婷段辛樂陳茂華
    關(guān)鍵詞:沙果登錄號種間

    王 康,李玉婷,鄭 燕,段辛樂,張 蒙,彭 雄,陳茂華

    (西北農(nóng)林科技大學(xué) 植物保護學(xué)院,農(nóng)業(yè)部西北黃土高原作物有害生物綜合治理重點實驗室,植保資源與病蟲害治理教育部重點實驗室,陜西 楊凌 712100)

    ?

    基于線粒體COⅠ和COⅡ基因的沙果小食心蟲與梨小食心蟲的分子鑒定

    王康,李玉婷,鄭燕,段辛樂,張蒙,彭雄,陳茂華

    (西北農(nóng)林科技大學(xué) 植物保護學(xué)院,農(nóng)業(yè)部西北黃土高原作物有害生物綜合治理重點實驗室,植保資源與病蟲害治理教育部重點實驗室,陜西 楊凌 712100)

    [摘要]【目的】 建立基于線粒體COⅠ和COⅡ基因序列的沙果小食心蟲、梨小食心蟲分子鑒定方法?!痉椒ā?田間采集試蟲,利用PCR和基因測序技術(shù),擴增2種食心蟲的線粒體COⅠ和COⅡ基因,將獲得的COⅠ和COⅡ基因序列在GenBank中進行BLAST比對,并通過GeneDoc軟件分析基因序列相似性,利用Mega 5.05軟件分別計算基于COⅠ和COⅡ基因序列的遺傳距離,并分別構(gòu)建COⅠ和COⅡ基因的系統(tǒng)發(fā)育樹?!窘Y(jié)果】 在分析的2種食心蟲樣本中,沙果小食心蟲COⅠ基因的種內(nèi)相似度在99.4%以上,梨小食心蟲COⅠ基因的種內(nèi)相似度在99.6%以上,2種食心蟲COⅠ基因的種間相似度為95.0%~95.6%;沙果小食心蟲COⅡ基因的種內(nèi)相似度在 99.0% 以上,梨小食心蟲COⅡ基因的種內(nèi)相似度在99.3%以上, 2種食心蟲COⅡ基因的種間相似度為94.6%~95.5%;2種食心蟲COⅠ基因種間存在30個穩(wěn)定變異位點,COⅡ基因種間有26個穩(wěn)定變異位點。2種食心蟲種內(nèi)遺傳距離為 0.001~0.006(COⅠ)和0.001~0.010(COⅡ),種間遺傳距離為0.046~0.052(COⅠ)和0.047~0.057(COⅡ),基于COⅠ和COⅡ基因的種間遺傳距離均顯著大于種內(nèi)遺傳距離。分別構(gòu)建COⅠ和COⅡ基因單倍型的系統(tǒng)發(fā)育樹,結(jié)果顯示,在2個基因的系統(tǒng)發(fā)育樹上,2種食心蟲的基因序列分別位于不同的進化支,且置信度均達到100%,同種食心蟲的基因序列位于相同的進化枝?!窘Y(jié)論】 可以根據(jù)本研究所用的COⅠ和COⅡ基因序列的差異性,進行沙果小食心蟲和梨小食心蟲的分子鑒定。

    [關(guān)鍵詞]沙果小食心蟲;梨小食心蟲;COⅠ基因;COⅡ基因;單倍型;分子鑒定

    沙果小食心蟲(GrapholitadimorphaKomai)和梨小食心蟲(GrapholitamolestaBusck)是重要的蛀果害蟲,二者均屬于鱗翅目(Lepidoptera)卷蛾科(Tortricidae)。沙果小食心蟲寄主植物為蘋果、沙果、李等,在我國部分地區(qū)報道有發(fā)生[1-2]。梨小食心蟲又名東方果蛀蛾,寄主植物主要為桃、蘋果、梨、杏、李、山楂等,是世界性的主要蛀果害蟲之一,在我國分布廣泛[3]。沙果小食心蟲和梨小食心蟲經(jīng)?;旌习l(fā)生,而且發(fā)生期部分重疊,二者幼蟲和成蟲外部形態(tài)都非常相似,而且性信息素主要成分亦相似[4-5],因此在田間調(diào)查和監(jiān)測種群動態(tài)時,準(zhǔn)確鑒定這2種食心蟲十分困難[1,6-7]。

    基于PCR技術(shù)的分子鑒定可以檢測不同物種間及相同物種內(nèi)不同個體間的細微差異,這為物種鑒定提供了新的手段[8-10]。線粒體DNA具有基因序列相對保守、無重組和單拷貝、進化速率適中、易于檢測等優(yōu)點,因此常被應(yīng)用于昆蟲的分子鑒定與系統(tǒng)進化研究之中,其中線粒體細胞色素氧化酶亞基Ⅰ(COⅠ)和線粒體細胞色素氧化酶亞基Ⅱ(COⅡ)等基因應(yīng)用較為廣泛[11-12]。常虹等[13]基于線粒體COⅠ基因構(gòu)建了齒小蠹屬(Ips)昆蟲的快速鑒定方法,很好地解決了植物檢疫中截獲的齒小蠹幼蟲或殘缺成蟲的鑒定難問題;游中華等[14]報道了基于線粒體COⅠ基因的9種薊馬的分子鑒定方法;通過COⅠ基因序列的分析發(fā)現(xiàn),弄蝶(Astraptesfulgerator)是由多個物種組成的混合種[15];王戎疆等[16]利用線粒體COⅡ基因構(gòu)建蛺蝶屬(Polyura) 5種蝴蝶的系統(tǒng)發(fā)育樹并進行聚類關(guān)系分析,結(jié)果和形態(tài)學(xué)結(jié)果相一致;楊飛龍等[17]通過比較大量中華按蚊(Anophelessinensis)標(biāo)本的COⅠ、COⅡ基因及幾種核糖體基因序列,建立了該蟲基于COⅡ等基因的分子鑒定標(biāo)準(zhǔn)。

    本研究從田間采集沙果小食心蟲和梨小食心蟲試蟲,提取試蟲的基因組DNA,PCR擴增試蟲的線粒體COⅠ和COⅡ基因序列,并對擴增的序列進行測序,以此建立沙果小食心蟲和梨小食心蟲分子鑒定方法,為這2種害蟲的準(zhǔn)確鑒定及預(yù)測預(yù)報提供了重要依據(jù)。

    1材料與方法

    1.1樣本采集

    利用沙果小食心蟲(GrapholitadimorphaKomai)和梨小食心蟲(GrapholitamolestaBusck)信息素誘捕器(誘捕器為三角膠粘式,購自北京中捷四方生物科技公司),在2種害蟲的發(fā)生期,于陜西省的5個蘋果園,分別誘集2種食心蟲成蟲樣本(表1)。誘芯的載體中空,由硅橡膠制成,每個質(zhì)量0.3~0.5 g,每個誘芯的性信息素含量不低于0.001 g,純度90%~97%。選擇較為完整的試蟲,在實驗室利用體視顯微鏡進行準(zhǔn)確的形態(tài)學(xué)鑒定,鑒定后的樣本保存于無水乙醇中,置于-20 ℃冰箱中保存?zhèn)溆谩?/p>

    表 1 沙果小食心蟲和梨小食心蟲樣本的采集信息

    1.2基因組DNA提取

    在采集的樣本中,隨機選取沙果小食心蟲和梨小食心蟲樣本,利用EasyPureTMGenomic DNA Kit(北京全式金生物技術(shù)有限公司)提取單頭試蟲基因組DNA,基因組DNA的提取方法按照試劑盒說明書步驟進行,用1%瓊脂糖凝膠電泳檢測提取的DNA質(zhì)量,提取的基因組DNA置于-20 ℃冰箱保存?zhèn)溆谩?/p>

    1.3PCR擴增及序列測定

    參照已有的報道,選用擴增2種食心蟲COⅠ和COⅡ基因的引物[18],稍作修改。擴增2種食心蟲COⅠ基因的上、下游引物序列分別為:5′-CAA-CATTTATTTTGATTTTTTGG-3′,5′-TCCAATGCACTAATCTGCCATATTA-3′;擴增2種食心蟲COⅡ基因的上、下游引物序列分別為:5′-ATG-GCAGATTATATGTAATGG-3′,5′-GTTTAAG-AGACCAGTACTTG-3′。

    PCR擴增反應(yīng)總體系為25 μL,其中含2×PCR master mix 12.5 μL,上、下游引物各3 μL,DNA模板2 μL(含20~50 ng DNA),加ddH2O至25 μL。PCR擴增程序為:94 ℃預(yù)變性2 min;94 ℃變性30 s,(CO Ⅰ 基因:52 ℃/CO Ⅱ 基因:57 ℃)退火延伸30 s,72 ℃延伸30 s,共35個循環(huán);循環(huán)結(jié)束后72 ℃延伸5 min。利用1%瓊脂糖凝膠電泳檢測PCR擴增產(chǎn)物。

    將擴增的PCR產(chǎn)物直接進行測序,測序由生工生物工程(上海)有限公司完成。

    1.4序列分析

    將所測2種食心蟲COⅠ和COⅡ基因在GenBank中進行BLAST搜索,并將獲得的COⅠ和COⅡ基因登錄到GenBank。通過GeneDoc軟件分析基因序列相似性,利用Mega 5.05軟件中的Kimura2-parameter(K2P)雙參數(shù)模型計算遺傳距離,并利用鄰接法(Neighbor-Joining,NJ)構(gòu)建系統(tǒng)發(fā)育樹[19]。

    2結(jié)果與分析

    2.12種試蟲線粒體COⅠ和COⅡ基因序列分析

    在所有沙果小食心蟲和梨小食心蟲試蟲中都成功擴增到線粒體COⅠ和COⅡ基因序列,其中COⅠ基因序列長度為720 bp,COⅡ基因序列長度為682 bp,所有序列均無堿基插入和缺失。將所獲得的序列在GenBank中進行BLAST分析,結(jié)果顯示,本研究獲得的沙果小食心蟲COⅠ基因和COⅡ基因序列與GenBank中登錄的沙果小食心蟲COⅠ基因(登錄號:AB603522)和COⅡ基因(登錄號:JQ404438)對應(yīng)序列的相似性分別為99.17%~99.31%和 99.19%~100%;測序獲得的梨小食心蟲COⅠ基因和COⅡ基因序列與GenBank中登錄的梨小食心蟲線粒體全基因組(登錄號:HQ116416)上對應(yīng)的COⅠ基因和COⅡ基因序列的相似性分別為 99.72%~100%和99.56%~100%。

    沙果小食心蟲和梨小食心蟲COⅠ和COⅡ基因序列的比對結(jié)果分別見圖1和圖2。供試樣本共獲得3個沙果小食心蟲COⅠ基因單倍型(GenBank登錄號:KC923197~KC923199),這3個單倍型在4個堿基位點上具有多態(tài)性,3個單倍型之間的相似度在99.4%以上。供試樣本共獲得6個梨小食心蟲COⅠ基因單倍型(GenBank登錄號:KC923200~KC923205),6個單倍型在5個堿基位點上具有多態(tài)性,6個單倍型之間的相似度在99.6%以上。沙果小食心蟲和梨小食心蟲COⅠ基因種間的相似度在95.0%~95.6%,COⅠ基因在2種食心蟲之間存在30個變異位點(圖1,表2)。

    供試樣本共獲得4個沙果小食心蟲COⅡ基因單倍型(GenBank登錄號:KC923206~KC923209),這4個單倍型在7個堿基位點上具有多態(tài)性,4個單倍型之間的相度在99.0%以上(表3)。

    圖 1沙果小食心蟲和梨小食心蟲COⅠ基因序列比對

    黑色顯示種間變異位點,下圖同

    Fig.1Alignment of COⅠ sequences ofGrapholitadimorphaandGrapholitamolesta

    The interspecific variation loci are indicated by black color.the same below

    圖 2沙果小食心蟲和梨小食心蟲COⅡ基因序列比對

    Fig.2Alignment of COⅡ sequences ofGrapholitadimorphaandGrapholitamolesta

    供試樣本共獲得8個梨小食心蟲COⅡ基因單倍型(GenBank登錄號:KC923210~KC923217),這8個單倍型在8個堿基位點上具有多態(tài)性,8個單倍型之間的相似度在99.3%以上。沙果小食心蟲和梨小食心蟲COⅡ基因種間的相似度在94.6%~95.5%,種間存在26個變異位點(圖2,表3)。

    2.22種試蟲遺傳距離及系統(tǒng)發(fā)育樹的構(gòu)建

    由表2可以看出,沙果小食心蟲COⅠ基因序列間的遺傳距離在0.001~0.006,梨小食心蟲COⅠ基因序列間的遺傳距離在0.001~0.004,2種食心蟲COⅠ基因種間遺傳距離為0.046~0.052,種間遺傳距離顯著高于種內(nèi)遺傳距離(df=29,t=2.05,P<0.001)。

    由表3可以看出,COⅡ基因中,沙果小食心蟲序列間遺傳距離在0.001~0.010,梨小食心蟲序列間遺傳距離在0.001~0.007,2種食心蟲COⅡ基因種間序列遺傳距離在0.047~0.057,種間遺傳距離顯著高于種內(nèi)遺傳距離(df=58,t=2.00,P<0.001)。

    表 2 沙果小食心蟲和梨小食心蟲COⅠ基因單倍型間的相似度和遺傳距離

    注:對角線上方為單倍型間相似度(%),對角線下方為單倍型間遺傳距離。GD1~GD3為沙果小食心蟲COⅠ基因單倍型,其GenBank登錄號依次為KC923197~KC923199;GM1~GM6是梨小食心蟲COⅠ基因單倍型,其GenBank登錄號依次為KC923200~KC923205。

    Note:The similarities between haplotypes (%) are above the diagonal,and the genetic distances between haplotypes are below the diagonal.GD1-GD3 are haplotypes ofG.dimorphaCOⅠ gene (GenBank accession No:KC923197-KC923199),and GM1-GM6 are the haplotypes ofG.molestaCOⅠ gene (GenBank accession No:KC923200-KC923205).

    表 3 沙果小食心蟲和梨小食心蟲COⅡ基因單倍型間的相似度和遺傳距離

    注:對角線上方為單倍型間相似度(%),對角線下方為單倍型間遺傳距離。GD1~GD4為沙果小食心蟲COⅡ基因單倍型,其GenBank登錄號依次為KC923206~KC923209;GM1~GM8是梨小食心蟲COⅡ基因單倍型,GenBank登錄號依次為KC923210~KC923217。

    Note:The similarities between haplotypes (%) are above the diagonal,and the genetic distances between haplotypes are below the diagonal.GD1-GD4 are haplotypes ofG.dimorphaCOⅡ gene (GenBank accession No:KC923206-KC923209),and GM1-GM8 are the haplotypes ofG.molestaCOⅡ gene (GenBank accession No:KC923210-KC923217).

    以蘋果蠹蛾(Cydiapomonella)COⅠ基因序列(GenBank登錄號:JQ928797)作外類群,將本研究獲得的沙果小食心蟲和梨小食心蟲COⅠ基因單倍型序列和1個已登錄的沙果小食心蟲COⅠ基因序列(GenBank登錄號:AB603522),以及3個已登錄的梨小食心蟲COⅠ基因序列(GenBank登錄號:AB603521、FJ449660、HQ700338)構(gòu)建系統(tǒng)發(fā)育樹,結(jié)果(圖3)顯示,沙果小食心蟲和梨小食心蟲的COⅠ基因單倍型各自聚在一支。

    以蘋果蠹蛾(C.pomonella)線粒體全基因(GenBank登錄號:JX407107)中COⅡ基因序列作外類群,將本研究獲得的2種食心蟲COⅡ基因單倍型序列和已登錄的相應(yīng)沙果小食心蟲COⅡ基因序列(GenBank登錄號:JQ404438)和梨小食心蟲COⅡ基因序列(GenBank登錄號:EU671060)構(gòu)建系統(tǒng)發(fā)育樹,結(jié)果(圖4)表明,所有的沙果小食心蟲COⅡ基因單倍型聚為一支,所有的梨小食心蟲單倍型聚為另外一支。

    圖 3 基于沙果小食心蟲和梨小食心蟲COⅠ基因單倍型的系統(tǒng)發(fā)育樹GD1~GD3和GM1~GM6為本研究獲得的沙果小食心蟲和梨小食心蟲COⅠ基因單倍型;AB603522為GenBank已登錄的相應(yīng)沙果小食心蟲COⅠ基因序列;AB603521、FJ449660、HQ700338為GenBank已登錄的相應(yīng)梨小食心蟲COⅠ序列;JQ928797為GenBank已登錄的相應(yīng)蘋果蠹蛾COⅠ基因序列

    3討論

    線粒體COⅠ和COⅡ基因序列目前已廣泛地應(yīng)用到昆蟲近緣種的鑒定之中[20]。在利用線粒體基因進行分子鑒定時,不同物種間相同的線粒體基因核苷酸序列差異是最重要的標(biāo)準(zhǔn)[21-22]。Hebert等[23]通過分析比較200種鱗翅目昆蟲COⅠ基因序列后發(fā)現(xiàn),COⅠ基因核苷酸序列3%的差異可以作為鑒定不同物種的標(biāo)準(zhǔn),即不同物種COⅠ基因差異大于3%,同一物種COⅠ基因序列差異小于3%,目前大量的相關(guān)研究證實了這一標(biāo)準(zhǔn)的準(zhǔn)確性[24-25]。本研究中,2種食心蟲COⅠ和COⅡ基因種內(nèi)核苷酸差異小于1.0%,而種間核苷酸差異大于4.4%,因此本研究的COⅠ和COⅡ基因序列可以準(zhǔn)確地應(yīng)用于沙果小食心蟲和梨小食心蟲的分子鑒定中。

    Aliabadian等[26]研究指出,當(dāng)種間遺傳距離和種內(nèi)遺傳距離不存在重疊區(qū)域,且分離程度較高時,目標(biāo)片段能準(zhǔn)確鑒定物種。本研究結(jié)果顯示,基于COⅠ和COⅡ基因的沙果小食心蟲和梨小食心蟲的種內(nèi)遺傳距離均小于0.010,而種間遺傳距離大于0.046,基于COⅠ和COⅡ基因的種內(nèi)遺傳距離和種間遺傳距離不存在重疊區(qū)域,并且種內(nèi)遺傳距離與種間遺傳距離差異顯著(P<0.001)。另外,有分析認為,種間平均遺傳距離一般為種內(nèi)平均遺傳距離的10倍以上[27],本研究基于COⅠ基因的種間平均遺傳距離(0.050 0)為種內(nèi)平均遺傳距離(0.002 7)的18.2倍,基于COⅡ基因的種間平均遺傳距離(0.050 6)為種內(nèi)平均遺傳距離(0.003 8)的13.2倍。因此,基于COⅠ和COⅡ基因的沙果小食心蟲和的梨小食心蟲種間和種內(nèi)遺傳距離分析結(jié)果表明,COⅠ和COⅡ基因可以用于這2種食心蟲的分子鑒定。另外,穩(wěn)定的種間差異位點是確定種的重要手段。朱振華等[28]通過比較6種果實蠅mtCytb基因序列,發(fā)現(xiàn)種間存在30個比較穩(wěn)定的變異位點,并以此作為鑒定不同種的標(biāo)準(zhǔn)。本試驗發(fā)現(xiàn),2種食心蟲在COⅠ基因序列上有30個穩(wěn)定變異位點,COⅡ基因序列上有26個穩(wěn)定差異位點,這些位點的堿基及其變化可以作為梨小食心蟲和沙果小食心蟲分子鑒定時的重要參考。

    4結(jié)論

    本研究結(jié)果表明,沙果小食心蟲、梨小食心蟲的COⅠ和COⅡ基因種內(nèi)的相似度在99.0%以上,種間相似度在95.6%以下,且2個基因的種間遺傳距離顯著高于種內(nèi)遺傳距離,2種食心蟲的COⅠ基因單倍型和COⅡ基因單倍型分別在系統(tǒng)發(fā)育樹上明顯各自聚為一支,因此,本研究采用的沙果小食心蟲、梨小食心蟲的COⅠ和COⅡ基因可以用于這2種害蟲的分子鑒定之中。

    [參考文獻]

    [1]嚴(yán)善春,劉友樵,李明文.我國果樹一新害蟲:沙果小食心蟲 [J].森林病害通訊,1999(4):15-16.

    Yan S C,Liu Y Q,Li M W.Grapholita dimorpha:A new record pest damaging fruit trees of China [J].Forest Pest and Disease,1999(4):15-16.(in Chinese)

    [2]Jung C R,Kim Y.Different types of fruit damages of three internal apple feeders diagnosed with mitochondrial molecular markers [J].Journal of Asia-Pacific Entomology,2013,16(2):189-197.

    [3]陳梅香,駱有慶,趙春江,等.梨小食心蟲研究進展 [J].北方園藝,2009(8):144-147.

    Chen M X,Luo Y Q,Zhao C J,et al.Research advance onGrapholitamolestaBusck [J].Northern Horticulture,2009(8):144-147.(in Chinese)

    [4]Choi K H,Lee D H,Byun B K,et al.Occurrence ofGrapholitadimorphaKomai (Lepidoptera:Tortricidae),a new insect pest in apple orchards of Korea [J].Korean Journal of Applied Entomology,2009,48(4):417-421.

    [5]Jung C R,Jung J K,Kim Y.Effects of different sex pheromone compositions and host plants on the mating behavior of twoGrapholitaspecies[J].Journal of Asia-Pacific Entomology,2013,16(4):507-512.

    [6]Komai F.A new species of the genusGrapholitaTreitschke from Japan allied to the oriental fruit moth,Grapholitamolesta(Busck) (Lepidoptera:Tortricidae) [J].Applied Entomology and Zoology,1979,14(2):133-136.

    [7]Hada H,Sekine K T.A diagnostic multiplex polymerase chain reaction method to identify Japanese internal apple-feeding Lepidopteran pests:Grapholitamolesta,Grapholitadimorpha(Lepidoptera:Tortricidae),andCarposinasasakii(Lepidoptera:Carposinidae) [J].Applied Entomology and Zoology,2011,46(2):287-291.

    [8]Tautz D,Arctander P,Minelli A,et al.DNA points the way ahead in taxonomy [J].Nature,2002,418(6897):479-479.

    [9]肖金花,肖暉,黃大衛(wèi).生物分類學(xué)的新動向:DNA條形編碼 [J].動物學(xué)報,2004,50(5):852-855.

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

    [10]莫幫輝,屈莉,韓松,等.DNA條形碼識別:Ⅰ.DNA條形碼研究進展及應(yīng)用前景 [J].四川動物,2008,27(2):303-306.

    Mo B H,Qu L,Han S,et al.DNA barcoding identification:Ⅰ.Research progress and applied perspective of DNA barcoding [J].Sichuan Journal of Zoology,2008,27(2):303-306.(in Chinese)

    [11]Bibb M J,Van Ette R A,Wright C T,et al.Sequence and gene organization of mouse mitochondrial DNA [J].Cell,1981,26:167-180.

    [12]牛屹東,李明,魏輔文,等.線粒體DNA用作分子標(biāo)記的可靠性和研究前景 [J].遺傳,2001,23(6):593-598.

    Niu Y D,Li M,Wei F W,et al.Reliability of mtDNA as molecular marker and its perspective [J].Hereditas,2001,23(6):593-598.(in Chinese)

    [13]常虹,郝德君,肖榮堂,等.基于線粒體COⅠ基因的齒小蠹屬昆蟲DNA條形碼研究 [J].昆蟲學(xué)報,2012,55(9):1075-1081.

    Chang H,Hao D J,Xiao R T,et al.DNA barcoding based on the mitochondrial COⅠ gene sequences forIpsspecies ( Coleoptera: Scolytidae) [J].Acta Entomologica Sinica,2012,55(9):1075-1081.(in Chinese)

    [14]游中華,路虹,張憲省,等.入侵害蟲西花薊馬及其他8種常見薊馬的分子鑒定 [J].昆蟲學(xué)報,2007,50(7):720-726.

    You Z H,Lu H,Zhang X S,et al.Molecular identification of the introduced western flower thrips,Frankliniellaoccidentals(Pergande) and other eight common thrips species (Thysanoptera: Thripidae) [J].Acta Entomologica Sinica,2007,50(7):720-726.(in Chinese)

    [15]Hebert P D N,Penton E H,Burns J M,et al.Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterflyAstraptesfulgerator[J].Proceedings of the National Academy of Sciences of the United States of America,2004a,101(41):14812-14817.

    [16]王戎疆,萬宏,雷光春,等.利用線粒體COⅡ基因序列對中國尾蛺蝶屬系統(tǒng)分化的研究(鱗翅目:蛺蝶科) [J].昆蟲學(xué)報,2004,47(2):243-247.

    Wang R J,Wan H,Lei G C,et al.Phylogenetic analysis ofPolyurain China inferred from mitochondrial COⅡ sequences (Lepidoptera:Nymphalidae) [J].Acta Entomologica Sinica,2004,47(2):243-247.(in Chinese)

    [17]楊飛龍,李旭東,閆振天,等.中華按蚊的分子鑒定標(biāo)準(zhǔn) [J].重慶師范大學(xué)學(xué)報:自然科學(xué)版,2014,31(4):40-44.

    Yang F L,Li X D,Yan Z T,et al.The molecular identification markers ofAnophelessinensis[J].Journal of Chongqing Normal University:Natural Science Edition,2014,31(4):40-44.(in Chinese)

    [18]Simon C,F(xiàn)raiti F,Beckenbach A,et al.Evolution,weighting,and phylogenetic utility of mitochondrial gene sequences and a compilation conserved polymerase chain reaction primer [J].Annals of the Entomological Society of America,1994,87:651-701.

    [19]Kumar S,Nei M,Dudley J,et al.Mega:A biologist-centtric software for evolutionary analysis of DNA and protein sequence [J].Briefings in Bioinformatics,2008,9(4):299-306.

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

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

    [21]Roe A D,Sperling F A H.Patterns of evolution of mitochondrial cytochromecoxidase Ⅰ and Ⅱ DNA and implications for DNA barcoding [J].Molecular Phylogenetics and Evolution,2007,44(1):325-345.

    [22]Tautz D,Arctander P,Minelli A,et al.A plea for DNA taxonomy [J].Trends in Ecology & Evolution,2003,18(2):70-74.

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

    [24]Ball S L,Hebert P D N,Burian S K,et al.Biological identifications of mayflies (Ephemeroptera) using DNA barcodes [J].Journal of the North American Benthological Society,2005,24(3):508-524.

    [25]Hajibabaei M,Janzen D H,Burns J M,et al.DNA barcodes distinguish species of tropical Lepidoptera [J].Proceedings of the National Academy of Sciences of the United States of America,2006,103(4):968-971.

    [26]Aliabadian M,Kaboli M,Nijman V,et al.Molecular identification of birds:Performance of distance-based DNA barcoding in three genes to delimit parapatric species [J].PLoS One,2009,4(1):e4119.

    [27]Hebert P D N,Stoeckle M Y,Zemlak T S,et al.Identification of birds through DNA barcodes [J].PLoS Biology,2004b,2(10):e312.

    [28]朱振華,葉輝,張智英.基于mtDNA Cytb的六種果實蠅的分子鑒定 [J].昆蟲學(xué)報,2005,48(3):386-390.

    Zhu Z H,Ye H,Zhang Z Y.Molecular identification of sixBactroceraspecies (Diptera:Tephritidae) based on mtDNA [J].Acta Entomologica Sinica,2005,48(3):386-390.(in Chinese)

    Molecular identification ofGrapholitadimorphaKomai (Lepidoptera:Tortricidae) andGrapholitamolestaBusck (Lepidoptera:Tortricidae) based on mitochondrial COⅠ and COⅡ genes

    WANG Kang,LI Yu-ting,ZHENG Yan,DUAN Xin-le,ZHANG Meng,PENG Xiong,CHEN Mao-hua

    (CollegeofPlantProtection,KeyLaboratoryofCropPestIntegratedPestManagementontheLoessPlateauofMinistryofAgriculture,KeyLaboratoryofPlantProtectionResourcesandPestManagementofMinistryofEducation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

    Abstract:【Objective】 This study aimed to establish a method for identifing Grapholita dimorpha Komai and Grapholita molesta Busck based on mitochondrial COⅠ and COⅡ genes.【Method】 PCR and gene sequencing methods were used to obtain COⅠ and COⅡ sequences from the samples of the two species collected from apple orchards.The acquired COⅠ and COⅡ sequences were verified by BLAST program in GenBank.GeneDoc was also used to analyze the similarity between the obtained sequences.Mega 5.05 was used to analyze the genetic distances between sequences and to construct the phylogenetic trees of COⅠ and COⅡ gene haplotypes,respectively.【Result】 The intraspecific similarities of COⅠ sequences of G.dimorpha and G.molesta were above 99.4% and 99.6%,respectively,whilst the interspecific similarity of COⅠ gene sequences between them was 95.0%-95.6%.The intraspecific similarities of COⅡ sequences of G.dimorpha and G.molesta were above 99.0% and 99.3%,while its interspecific similarity between them was 94.6%-95.5%.A total of 30 variable sites were detected in COⅠ gene sequences from samples of the two species,and 26 were detected in COⅡ sequences.The intraspecific genetic distances were 0.001-0.006 (COⅠ) and 0.001-0.010 (COⅡ),whereas the interspecific genetic distances were 0.046-0.052 (COⅠ) and 0.047-0.057 (COⅡ).The interspecific genetic distance was significantly higher than intraspecific genetic distance.Phylogenetic trees based on COⅠ gene haplotypes and COⅡ gene haplotypes indicated that gene sequences from the same species were clustered in the same branch,whereas gene sequences from different species were clustered in different branches, with the bootstrap value of 100%.【Conclusion】 The COⅠ and COⅡ gene sequences used in this study can be applied in the molecular identification of G.dimorpha and G.molesta.

    Key words:Grapholita dimorpha;Grapholita molesta;COⅠ gene;COⅡ gene;haplotypes;molecular identification

    DOI:網(wǎng)絡(luò)出版時間:2016-01-0810:2210.13207/j.cnki.jnwafu.2016.02.022

    [收稿日期]2014-06-30

    [基金項目]國家自然科學(xué)基金項目(31071687);教育部高等學(xué)校博士學(xué)科點專項科研基金(博導(dǎo)類)項目(20110204110001);西北農(nóng)林科技大學(xué)人才引進項目;國家科技支撐計劃項目(2012BAK11B03)

    [作者簡介]王康(1989-),男,陜西咸陽人,在讀碩士,主要從事農(nóng)業(yè)昆蟲和害蟲防治研究。E-mail:kangwang@nwsuaf.edu.cn[通信作者]陳茂華(1971-),男,湖北武漢人,教授,博士生導(dǎo)師,主要從事昆蟲種群遺傳學(xué)和分子毒理學(xué)研究。

    [中圖分類號]S763.42;Q78

    [文獻標(biāo)志碼]A

    [文章編號]1671-9387(2016)02-0156-09

    E-mail:maohua.chen@nwsuaf.edu.cn

    猜你喜歡
    沙果登錄號種間
    干旱條件對鬼針草和醉魚草種間相互作用及生長的影響
    植物研究(2023年5期)2023-09-09 08:01:22
    家鄉(xiāng)的沙果
    三峽庫區(qū)支流花溪河浮游植物種間關(guān)聯(lián)及影響因子分析
    基于DNA條形碼進行金花茶組種間鑒別
    種子(2021年2期)2021-03-31 09:23:32
    含氟新農(nóng)藥的研究進展
    浙江化工(2019年5期)2019-06-04 08:33:34
    圖牧吉的沙果紅了
    沙果果醋綠色產(chǎn)品價值及市場分析*
    小反芻獸疫病毒F基因的序列分析
    電子版館藏數(shù)據(jù)回溯中的交叉、重疊、阻滯問題處理
    江蘇省宜興市茶園秋季雜草種間生態(tài)關(guān)系及群落分類
    久9热在线精品视频| 午夜久久久在线观看| 色综合婷婷激情| 久久香蕉精品热| 国产精品香港三级国产av潘金莲| 久久人妻福利社区极品人妻图片| 日本vs欧美在线观看视频| 日日爽夜夜爽网站| 亚洲一区二区三区色噜噜 | 叶爱在线成人免费视频播放| 亚洲免费av在线视频| 国产99久久九九免费精品| 少妇 在线观看| 亚洲国产精品sss在线观看 | 一区二区三区国产精品乱码| 男人舔女人的私密视频| 90打野战视频偷拍视频| 热99re8久久精品国产| 国产高清videossex| 国产精品国产高清国产av| 午夜影院日韩av| 女生性感内裤真人,穿戴方法视频| 久久久久国内视频| 国产片内射在线| 亚洲av五月六月丁香网| 亚洲第一av免费看| 一个人观看的视频www高清免费观看 | 五月开心婷婷网| 两性夫妻黄色片| 在线观看日韩欧美| 国产av一区在线观看免费| 一区二区三区国产精品乱码| 亚洲专区国产一区二区| 免费在线观看完整版高清| 久久天躁狠狠躁夜夜2o2o| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲一级av第二区| 十八禁网站免费在线| 桃红色精品国产亚洲av| 看免费av毛片| 波多野结衣高清无吗| 亚洲精品粉嫩美女一区| 久久久国产成人免费| 午夜免费鲁丝| 好男人电影高清在线观看| 亚洲av成人av| 精品久久久久久电影网| 在线观看舔阴道视频| 波多野结衣高清无吗| 91成年电影在线观看| 99热国产这里只有精品6| 精品国产乱子伦一区二区三区| 亚洲成人免费电影在线观看| 亚洲视频免费观看视频| 男人操女人黄网站| 色综合婷婷激情| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 免费在线观看黄色视频的| xxxhd国产人妻xxx| 99国产极品粉嫩在线观看| 国产免费现黄频在线看| 亚洲成国产人片在线观看| 免费在线观看视频国产中文字幕亚洲| 久热这里只有精品99| 电影成人av| 老司机午夜福利在线观看视频| 久久 成人 亚洲| 国产视频一区二区在线看| 涩涩av久久男人的天堂| 88av欧美| cao死你这个sao货| 极品人妻少妇av视频| 久久久久精品国产欧美久久久| 久久国产亚洲av麻豆专区| 国产精品久久久久成人av| 老汉色∧v一级毛片| 一区福利在线观看| 国产免费现黄频在线看| 青草久久国产| 久久精品国产99精品国产亚洲性色 | 午夜福利在线免费观看网站| 中亚洲国语对白在线视频| 操美女的视频在线观看| 国产在线精品亚洲第一网站| 激情在线观看视频在线高清| 亚洲伊人色综图| 亚洲久久久国产精品| 亚洲国产精品sss在线观看 | 欧美午夜高清在线| 久久天躁狠狠躁夜夜2o2o| 高清av免费在线| 一本综合久久免费| 18禁观看日本| 亚洲av美国av| 天堂中文最新版在线下载| 欧美日韩福利视频一区二区| 正在播放国产对白刺激| 国产精品99久久99久久久不卡| 一区二区日韩欧美中文字幕| 免费搜索国产男女视频| 校园春色视频在线观看| a级毛片黄视频| 亚洲成人免费av在线播放| 91精品三级在线观看| 国产黄色免费在线视频| 亚洲第一欧美日韩一区二区三区| 日韩免费高清中文字幕av| 亚洲人成网站在线播放欧美日韩| 日本黄色视频三级网站网址| 正在播放国产对白刺激| 在线观看舔阴道视频| 亚洲精品国产区一区二| 国产野战对白在线观看| 久热这里只有精品99| 免费av毛片视频| 国产精品爽爽va在线观看网站 | 黄色女人牲交| 亚洲,欧美精品.| 免费高清在线观看日韩| 久久久精品欧美日韩精品| 国产亚洲av高清不卡| 欧美日韩精品网址| 日韩免费av在线播放| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 免费人成视频x8x8入口观看| 老司机午夜福利在线观看视频| 香蕉久久夜色| 日韩欧美国产一区二区入口| 日韩欧美三级三区| 高清黄色对白视频在线免费看| tocl精华| 国产成人欧美在线观看| 亚洲自偷自拍图片 自拍| 色哟哟哟哟哟哟| 在线观看一区二区三区激情| 亚洲久久久国产精品| 99国产精品一区二区蜜桃av| 高清av免费在线| www日本在线高清视频| 乱人伦中国视频| av有码第一页| 日韩高清综合在线| 很黄的视频免费| 久久热在线av| 色老头精品视频在线观看| 大型黄色视频在线免费观看| 亚洲av成人av| 黄片大片在线免费观看| 正在播放国产对白刺激| 亚洲情色 制服丝袜| 久久热在线av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品久久成人aⅴ小说| av网站免费在线观看视频| 国产亚洲欧美精品永久| 欧美激情久久久久久爽电影 | 国产av在哪里看| 琪琪午夜伦伦电影理论片6080| 免费在线观看亚洲国产| av电影中文网址| 最好的美女福利视频网| 琪琪午夜伦伦电影理论片6080| 成人手机av| 一级毛片精品| 一进一出抽搐gif免费好疼 | 50天的宝宝边吃奶边哭怎么回事| 精品少妇一区二区三区视频日本电影| 精品人妻在线不人妻| 69av精品久久久久久| 老司机福利观看| 亚洲一区二区三区不卡视频| 免费在线观看视频国产中文字幕亚洲| 日韩免费av在线播放| 亚洲激情在线av| 国产精品 国内视频| 中出人妻视频一区二区| netflix在线观看网站| 亚洲中文字幕日韩| 欧美黑人精品巨大| 视频区欧美日本亚洲| 十八禁人妻一区二区| 99国产精品一区二区蜜桃av| 丝袜美腿诱惑在线| 国产欧美日韩一区二区精品| 日日爽夜夜爽网站| 国产av在哪里看| 中文字幕人妻丝袜制服| 黄色片一级片一级黄色片| 一级毛片高清免费大全| 黑人巨大精品欧美一区二区蜜桃| 美女高潮喷水抽搐中文字幕| 亚洲五月天丁香| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看 | 欧美激情极品国产一区二区三区| tocl精华| 国产激情欧美一区二区| 欧美亚洲日本最大视频资源| 美女高潮到喷水免费观看| 制服诱惑二区| 国产激情欧美一区二区| 日韩精品青青久久久久久| 日本免费a在线| 大码成人一级视频| 亚洲欧美激情在线| 天堂动漫精品| 精品国产一区二区久久| 又黄又爽又免费观看的视频| 一进一出好大好爽视频| www.www免费av| 两个人看的免费小视频| 天天躁夜夜躁狠狠躁躁| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 麻豆成人av在线观看| 亚洲熟妇中文字幕五十中出 | 日本免费一区二区三区高清不卡 | 久久九九热精品免费| 亚洲中文日韩欧美视频| av在线播放免费不卡| 一区二区三区精品91| 免费在线观看日本一区| 女人高潮潮喷娇喘18禁视频| 亚洲av熟女| 成人18禁高潮啪啪吃奶动态图| 成人国产一区最新在线观看| 国产精品爽爽va在线观看网站 | 少妇 在线观看| 国产一区二区在线av高清观看| 精品一区二区三区四区五区乱码| 757午夜福利合集在线观看| 丁香欧美五月| 亚洲精华国产精华精| 亚洲国产中文字幕在线视频| 在线观看免费视频日本深夜| 19禁男女啪啪无遮挡网站| 国产麻豆69| 午夜激情av网站| 黄片小视频在线播放| 如日韩欧美国产精品一区二区三区| 国产99久久九九免费精品| 国产精品久久久av美女十八| 久久久久久久久久久久大奶| 黄片播放在线免费| av中文乱码字幕在线| 日日干狠狠操夜夜爽| 高清av免费在线| 国产无遮挡羞羞视频在线观看| 久久精品影院6| 国产亚洲欧美98| 亚洲一区二区三区不卡视频| av欧美777| 18禁裸乳无遮挡免费网站照片 | 一区二区日韩欧美中文字幕| 国产成人影院久久av| 免费观看精品视频网站| 女人被躁到高潮嗷嗷叫费观| 国产一区二区激情短视频| 热99re8久久精品国产| 亚洲专区字幕在线| 亚洲国产精品sss在线观看 | 最新美女视频免费是黄的| 亚洲精品国产区一区二| 成人免费观看视频高清| 精品久久久久久电影网| 757午夜福利合集在线观看| 在线观看免费视频网站a站| 狠狠狠狠99中文字幕| 欧美日韩福利视频一区二区| 欧美 亚洲 国产 日韩一| 国产高清视频在线播放一区| 日韩欧美三级三区| 99精国产麻豆久久婷婷| 黄片播放在线免费| 日韩欧美三级三区| 久久久国产一区二区| 99国产极品粉嫩在线观看| tocl精华| 精品国产超薄肉色丝袜足j| 久久99一区二区三区| 亚洲欧美一区二区三区久久| 国产精品爽爽va在线观看网站 | 麻豆国产av国片精品| 一二三四在线观看免费中文在| 久久中文字幕一级| 真人做人爱边吃奶动态| 精品乱码久久久久久99久播| 高清在线国产一区| 午夜福利欧美成人| 操美女的视频在线观看| 国产av在哪里看| 在线观看日韩欧美| 国产99白浆流出| 电影成人av| 热99国产精品久久久久久7| 国产免费现黄频在线看| 中文字幕高清在线视频| 天堂动漫精品| 久久人妻av系列| 精品一品国产午夜福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品在线观看二区| 99久久国产精品久久久| 久久久久久人人人人人| 久久香蕉激情| 身体一侧抽搐| 麻豆国产av国片精品| 亚洲欧美激情在线| 黄片大片在线免费观看| 亚洲三区欧美一区| 日韩高清综合在线| 国产亚洲av高清不卡| 免费搜索国产男女视频| 国产乱人伦免费视频| 欧美日韩亚洲综合一区二区三区_| 国产成人免费无遮挡视频| 十八禁网站免费在线| 国产精品一区二区在线不卡| 亚洲成人国产一区在线观看| √禁漫天堂资源中文www| 在线视频色国产色| 日本撒尿小便嘘嘘汇集6| 12—13女人毛片做爰片一| 成人手机av| 亚洲人成网站在线播放欧美日韩| www.www免费av| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 亚洲欧美精品综合一区二区三区| 亚洲黑人精品在线| 777久久人妻少妇嫩草av网站| 高清黄色对白视频在线免费看| 又黄又粗又硬又大视频| 免费观看人在逋| 男男h啪啪无遮挡| 亚洲精品国产精品久久久不卡| 黑人欧美特级aaaaaa片| 久久天堂一区二区三区四区| www.熟女人妻精品国产| a级片在线免费高清观看视频| 国产精品永久免费网站| 少妇被粗大的猛进出69影院| 久久午夜综合久久蜜桃| 在线永久观看黄色视频| 搡老乐熟女国产| 亚洲视频免费观看视频| 在线播放国产精品三级| 国产精品免费一区二区三区在线| 91大片在线观看| 国产精品国产av在线观看| 成熟少妇高潮喷水视频| 99re在线观看精品视频| 午夜成年电影在线免费观看| 精品人妻在线不人妻| www.精华液| 精品一品国产午夜福利视频| 国产av一区在线观看免费| 日韩欧美一区二区三区在线观看| 国内久久婷婷六月综合欲色啪| 日韩成人在线观看一区二区三区| 91麻豆av在线| 婷婷六月久久综合丁香| 亚洲欧洲精品一区二区精品久久久| 国产av一区二区精品久久| 人妻久久中文字幕网| 国产又爽黄色视频| 校园春色视频在线观看| 久热这里只有精品99| 亚洲精品国产一区二区精华液| 久久国产精品男人的天堂亚洲| 夜夜夜夜夜久久久久| 伊人久久大香线蕉亚洲五| 一进一出抽搐动态| 男人舔女人的私密视频| 夜夜夜夜夜久久久久| 国产av在哪里看| 亚洲免费av在线视频| 99久久人妻综合| 国产成人一区二区三区免费视频网站| 欧美日本亚洲视频在线播放| 久久精品国产99精品国产亚洲性色 | 99国产精品99久久久久| 亚洲狠狠婷婷综合久久图片| 无限看片的www在线观看| 天天添夜夜摸| 亚洲专区字幕在线| 免费人成视频x8x8入口观看| 欧美成人午夜精品| 午夜精品在线福利| 欧美+亚洲+日韩+国产| 国产极品粉嫩免费观看在线| 神马国产精品三级电影在线观看 | 国产一区在线观看成人免费| 免费久久久久久久精品成人欧美视频| 久久精品国产亚洲av香蕉五月| 久久久国产成人精品二区 | 欧美日韩福利视频一区二区| 99国产精品一区二区三区| 深夜精品福利| av片东京热男人的天堂| 国产精品一区二区精品视频观看| 中文字幕另类日韩欧美亚洲嫩草| 制服人妻中文乱码| 国产又爽黄色视频| 日本vs欧美在线观看视频| 久久久久久人人人人人| 真人一进一出gif抽搐免费| 水蜜桃什么品种好| 国产亚洲精品久久久久5区| 一级毛片高清免费大全| 99精品欧美一区二区三区四区| 精品午夜福利视频在线观看一区| 亚洲成人免费av在线播放| 国产精品秋霞免费鲁丝片| 12—13女人毛片做爰片一| 亚洲人成电影观看| 女性被躁到高潮视频| 久久中文字幕人妻熟女| 亚洲精品久久午夜乱码| 国产蜜桃级精品一区二区三区| 欧美成人免费av一区二区三区| 人人妻,人人澡人人爽秒播| av福利片在线| 午夜两性在线视频| 国产精品久久视频播放| 超碰成人久久| 久久中文字幕一级| 午夜日韩欧美国产| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 久久精品国产亚洲av香蕉五月| 中出人妻视频一区二区| www.自偷自拍.com| 国产欧美日韩一区二区精品| 天堂俺去俺来也www色官网| 午夜福利在线观看吧| 日韩精品中文字幕看吧| 日韩欧美免费精品| 无限看片的www在线观看| 天天添夜夜摸| 欧美在线黄色| 久久人人精品亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| 色哟哟哟哟哟哟| 真人一进一出gif抽搐免费| 亚洲精品成人av观看孕妇| 大型av网站在线播放| 热99国产精品久久久久久7| 亚洲av第一区精品v没综合| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| 国内久久婷婷六月综合欲色啪| 久久亚洲精品不卡| 免费在线观看亚洲国产| 99国产极品粉嫩在线观看| 热99re8久久精品国产| 久久人人97超碰香蕉20202| 欧美日韩视频精品一区| 亚洲精品中文字幕在线视频| 欧美成狂野欧美在线观看| 亚洲色图av天堂| 欧美日韩乱码在线| 国产亚洲精品一区二区www| 99国产综合亚洲精品| 亚洲av片天天在线观看| 亚洲精品国产精品久久久不卡| 亚洲 欧美一区二区三区| 国产亚洲欧美在线一区二区| 免费一级毛片在线播放高清视频 | 日韩欧美三级三区| 两个人看的免费小视频| 人人妻人人爽人人添夜夜欢视频| 国产av一区在线观看免费| 1024香蕉在线观看| 精品日产1卡2卡| 身体一侧抽搐| 免费搜索国产男女视频| netflix在线观看网站| 岛国视频午夜一区免费看| 91麻豆av在线| 韩国精品一区二区三区| 亚洲一区二区三区欧美精品| 国产亚洲精品一区二区www| 男女床上黄色一级片免费看| 精品欧美一区二区三区在线| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 欧美日韩精品网址| 视频在线观看一区二区三区| 99国产精品一区二区三区| а√天堂www在线а√下载| 又黄又粗又硬又大视频| 男女之事视频高清在线观看| 成人精品一区二区免费| 90打野战视频偷拍视频| 9191精品国产免费久久| 黑丝袜美女国产一区| 欧美乱色亚洲激情| 90打野战视频偷拍视频| 搡老熟女国产l中国老女人| 叶爱在线成人免费视频播放| 香蕉久久夜色| 国产成人精品久久二区二区91| 不卡av一区二区三区| 久久精品国产亚洲av高清一级| 一边摸一边抽搐一进一出视频| 亚洲aⅴ乱码一区二区在线播放 | 他把我摸到了高潮在线观看| 久久伊人香网站| 精品国产美女av久久久久小说| 熟女少妇亚洲综合色aaa.| 国产不卡一卡二| 法律面前人人平等表现在哪些方面| 天天影视国产精品| 91麻豆av在线| 国产不卡一卡二| svipshipincom国产片| 精品一区二区三卡| 国产精品 欧美亚洲| 亚洲,欧美精品.| 国产精品一区二区免费欧美| 深夜精品福利| 欧美久久黑人一区二区| 亚洲男人天堂网一区| 在线观看免费高清a一片| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 最近最新免费中文字幕在线| 久热爱精品视频在线9| 欧美精品一区二区免费开放| 午夜视频精品福利| 国产欧美日韩综合在线一区二区| 日本五十路高清| 亚洲人成电影免费在线| 又紧又爽又黄一区二区| 天天添夜夜摸| 看免费av毛片| 黄色女人牲交| 国产亚洲欧美98| 女人被躁到高潮嗷嗷叫费观| 午夜福利在线观看吧| 午夜福利影视在线免费观看| 男女床上黄色一级片免费看| 美女扒开内裤让男人捅视频| 大型黄色视频在线免费观看| 亚洲人成电影观看| 老熟妇乱子伦视频在线观看| 亚洲成av片中文字幕在线观看| svipshipincom国产片| 久久国产乱子伦精品免费另类| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 亚洲九九香蕉| 日日摸夜夜添夜夜添小说| 男女做爰动态图高潮gif福利片 | 国产蜜桃级精品一区二区三区| 亚洲成国产人片在线观看| 18禁观看日本| 久久午夜亚洲精品久久| 国产精品98久久久久久宅男小说| 免费观看人在逋| 亚洲av成人一区二区三| 久久香蕉国产精品| 成人精品一区二区免费| 男女之事视频高清在线观看| 亚洲专区国产一区二区| 欧美人与性动交α欧美精品济南到| 日本五十路高清| 校园春色视频在线观看| 国产精品国产av在线观看| 免费在线观看日本一区| 亚洲九九香蕉| 一二三四在线观看免费中文在| 男人操女人黄网站| av中文乱码字幕在线| 免费久久久久久久精品成人欧美视频| 欧美黑人欧美精品刺激| 国产亚洲精品综合一区在线观看 | 女人被狂操c到高潮| 日韩 欧美 亚洲 中文字幕| 制服人妻中文乱码| 国产黄a三级三级三级人| 咕卡用的链子| cao死你这个sao货| 国产视频一区二区在线看| 久9热在线精品视频| 最好的美女福利视频网| 美女扒开内裤让男人捅视频| 午夜福利一区二区在线看| 国产片内射在线| 男女床上黄色一级片免费看| 夫妻午夜视频| 国产av在哪里看| 无人区码免费观看不卡| 国产成人欧美| 国产极品粉嫩免费观看在线| 日韩国内少妇激情av| 午夜福利在线观看吧| 国产欧美日韩一区二区三区在线| 老司机午夜福利在线观看视频| av视频免费观看在线观看| 精品久久久久久久久久免费视频 | 亚洲人成电影免费在线| 日韩成人在线观看一区二区三区| 亚洲国产中文字幕在线视频| 亚洲成人国产一区在线观看| 99久久99久久久精品蜜桃| 日韩av在线大香蕉| 丝袜在线中文字幕| 日本vs欧美在线观看视频| 国产精品久久久人人做人人爽| 黄色毛片三级朝国网站| 亚洲黑人精品在线|