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

    基于線粒體COⅠ和COⅡ基因的5種不同寄主植物西花薊馬種群遺傳多樣性研究

    2018-05-14 12:17田虎張蓉王玉生萬方浩張桂芬
    植物保護 2018年1期
    關(guān)鍵詞:薊馬靶標(biāo)品系

    田虎 張蓉 王玉生 萬方浩 張桂芬

    摘要

    西花薊馬Frankliniella occidentalis (Pergande)是一種分布廣、危害大的世界性檢疫害蟲。本研究旨在探討田間常見5種寄主植物上的西花薊馬種群的寄主?;院瓦z傳多樣性。以采自甘肅、寧夏的辣椒、茄子、蜀葵、月季、美人蕉等5種寄主植物上的西花薊馬為對象,以線粒體CO Ⅰ和CO Ⅱ基因為靶標(biāo),應(yīng)用Arlequin 3.5 軟件進行種群遺傳多樣性、遺傳分化、基因流水平及分子變異分析,以MEGA 7.0 和Network 5.0軟件分別構(gòu)建單倍型系統(tǒng)發(fā)育樹和中接網(wǎng)狀樹。結(jié)果顯示,當(dāng)以mtDNA CO Ⅰ和CO Ⅱ基因為靶標(biāo)時,分別檢測到13種和12種單倍型,其中單倍型H_1和H_2為優(yōu)勢單倍型;辣椒和茄子上的單倍型數(shù)量只有3種(CO Ⅰ基因)或2種(CO Ⅱ基因),而蜀葵、月季和美人蕉上的單倍型數(shù)量比較多,為7~9種(CO Ⅰ基因)或8~9種(CO Ⅱ基因);辣椒和茄子上的種群單倍型種類、單倍型多樣性、核苷酸多樣性和核苷酸平均差異數(shù)等均較蜀葵、月季和美人蕉上的種群低。CO Ⅰ和CO Ⅱ基因總寄主種群Fus Fs檢驗結(jié)果分別為2.36和4.06,種群大小整體保持穩(wěn)定;總寄主種群固定指數(shù)FST和基因流Nm分析表明,西花薊馬各寄主植物種群之間基因交流充分,尚未發(fā)生明顯的遺傳分化;AMOVA分析顯示遺傳變異主要來自種群內(nèi)部;基于CO Ⅰ和CO Ⅱ基因序列的單倍型聚類和網(wǎng)狀樹都明顯分為兩支,分別對應(yīng)西花薊馬的溫室品系Glasshouse strain和羽扇豆品系Lupin strain;其中溫室品系在所有寄主植物種群中均有發(fā)生,而羽扇豆品系的單倍型主要來自多年生的蜀葵、月季和美人蕉等植物,在辣椒和茄子上并無羽扇豆品系發(fā)生。研究結(jié)果對西花薊馬種群擴張機制探討、遺傳動態(tài)分析以及有效防控措施的制定具有一定指導(dǎo)意義。

    關(guān)鍵詞

    西花薊馬; 線粒體DNA; CO Ⅰ和CO Ⅱ基因; 寄主植物; 遺傳多樣性; 系統(tǒng)發(fā)育分析; 網(wǎng)狀樹

    中圖分類號:

    S436.3

    文獻標(biāo)識碼: A

    DOI: 10.16688/j.zwbh.2017161

    Genetic diversity among five different host plant populations of Frankliniella

    occidentalis (Pergande) based on COⅠ and COⅡ gene sequences

    TIAN Hu1, ZHANG Rong1, WANG Yusheng1, WAN Fanghao1,2, ZHANG Guifen1,2

    (1. State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant

    Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China; 2. Center for

    Management of Invasive Alien Species, Ministry of Agriculture, Beijing 100193, China)

    Abstract

    Frankliniella occidentalis (Pergande) is a worldwide quarantine pest with wide distribution and can cause serious damage. This study aimed to evaluate the host specialization and genetic diversity among different host plant populations of F.occidentalis. In the present study, the mtDNA COⅠ and COⅡ were selected as molecular markers. Samples were collected from five host plant species, including Capsicum annuum, Solanum melongena, Althaea rosea, Rosa chinensis, and Canna indica in Gansu and Ningxia provinces. The genetic diversity, genetic differentiation, gene flow and molecular variance (AMOVA) were analyzed by using the software Arlequin 3.5. Phylogenetic relationships and medianjoining networks were inferred from MEGA 7.0 and Network 5.0, respectively. The results showed that 13 and 12 types of haplotype were detected when the COⅠ and COⅡ gene sequences analyzed, respectively, and H_1 and H_2 were the dominant haplotypes. Moreover, only 3 (COⅠ) and 2 (COⅡ) types of haplotype were detected in C. annuum (LJ) and S.melongena populations, but 79 (COⅠ) or 89 (COⅡ) types were detected in A.rosea, R.chinensis or C. indica populations. Furthermore, the haplotype diversity (Hd), nucleotide diversity (π) and average number of nucleotide difference (K) of C. annuum and S.melongena populations were all lower than those of A.rosea, R.chinensis and C. indica populations. Analyses of fixation index (FST) and gene flow (Nm) based on COⅠ or COⅡ gene suggested high level of gene flow and no obvious genetic differentiation among host populations of F.occidentalis. Analysis of molecular variance (AMOVA) indicated that the genetic variance was mainly attributes to variation within populations. The phylogenetic tree and haplotype network analyses based on haplotypes of COⅠ and COⅡ genes showed two clear genetic branches, corresponding to the glasshouse strain and lupin strain of F.occidentalis, respectively. The haplotypes standing for glasshouse strain happened on all host plants, while the haplotypes standing for lupin strain were from perennial plant species A.rosea, R.chinensis and C. indica, but not from C. annuum and S.melongena. The present results may help us understand the population expansion mechanism, analyze genetic dynamics and set up effective management measures of F.occidentalis.

    Key words

    Frankliniella occidentalis; mtDNA; COⅠ and COⅡ genes; host plant; genetic diversity; phylogenic analysis; haplotype network

    西花薊馬Frankliniella occidentalis (Pergande)又稱苜蓿薊馬,屬纓翅目Thysanoptera薊馬科Thripidae,是一種世界性檢疫害蟲[12],原產(chǎn)于美國和加拿大的西部山區(qū),20世紀(jì)80年代遍及北美洲,之后隨國際貿(mào)易的快速發(fā)展迅速在世界范圍內(nèi)傳播擴散[3]。我國于2003年在北京郊區(qū)溫室大棚的辣椒上首次發(fā)現(xiàn)其為害[4],隨后迅速擴張到云南、浙江(鮮切花店)、山東、貴州、湖南、新疆等地區(qū)[59],目前已成為我國蔬菜花卉生產(chǎn)中的主要害蟲[10]。

    近年來,我國學(xué)者已對西花薊馬的入侵來源、傳播擴散途徑以及局部暴發(fā)成災(zāi)機制等開展了大量研究[1114],為其綜合防控措施的制定提供了理論支撐。然而,隨著全球氣候變化的逐年加劇、果蔬貿(mào)易的快速發(fā)展以及設(shè)施農(nóng)業(yè)的普遍興起,西花薊馬的發(fā)生為害呈逐年加重的趨勢[10]。

    西花薊馬為雜食性害蟲,且植食時寄主植物范圍十分廣泛,目前已知的寄主植物種類多達62科,500余種[1516]。以往的研究指出西花薊馬存在兩個品系,即,溫室品系(glasshouse strain)和羽扇豆品系(lupin strain)[3]。其中,溫室品系主要為害溫室作物,寄主植物范圍廣泛,并且對常用的殺蟲劑已產(chǎn)生了較強的抗性[17];羽扇豆品系主要為害豆科植物羽扇豆;但兩個品系在形態(tài)上無法準(zhǔn)確區(qū)分[3, 18]。近年來的研究顯示,兩個品系的西花薊馬已同時入侵我國,并發(fā)現(xiàn)羽扇豆品系的寄主植物雖并非只有羽扇豆,但其寄主植物種類相對較少,初步統(tǒng)計約有10種[1114]。通常,植食性昆蟲與寄主植物之間存在明顯的協(xié)同進化關(guān)系[1920],并進而可能導(dǎo)致寄主?;?宗的出現(xiàn)[2122]。此外,地理分布,寄主植物的種類、次生代謝物以及氣候條件等[2325]也會影響昆蟲的種群遺傳多樣性和遺傳結(jié)構(gòu)。然而,不同寄主植物上西花薊馬種群的遺傳多樣性及其與溫室品系和羽扇豆品系的關(guān)系尚不得而知。

    本研究圍繞西花薊馬品系的寄主偏好性及其遺傳特性,以線粒體細(xì)胞色素C氧化酶亞基Ⅰ(mitochondrial cytochrome c oxidase subunitⅠ, mtDNA COⅠ)和亞基Ⅱ (mitochondrial cytochrome c oxidase subunit Ⅱ, mtDNA CO Ⅱ)基因為分子標(biāo)記,以采自我國寧夏和甘肅5種常見寄主植物上的西花薊馬為材料,通過比較不同寄主植物上西花薊馬種群的分子遺傳學(xué)參數(shù),分析種群間遺傳分化和基因交流特點,構(gòu)建系統(tǒng)發(fā)育樹等,明確各寄主植物上的西花薊馬種群的遺傳結(jié)構(gòu)和遺傳多樣性,研究結(jié)果對西花薊馬種群擴張機制探討、遺傳動態(tài)分析以及有效防控措施制定具有一定指導(dǎo)意義。

    1 材料和方法

    1.1 樣本采集

    于田間薊馬高發(fā)期,在西北地區(qū)的內(nèi)蒙古(包頭)、新疆(烏魯木齊、吐魯番)、甘肅(蘭州、張掖、嘉峪關(guān))、寧夏(銀川、吳忠、青銅峽)、青海(西寧、海東)等5省份11個區(qū)域,選擇田間常見的植物(包括茄子、辣椒等蔬菜植物,蜀葵、美人蕉等花卉植物,以及國槐、小葉黃楊等園林植物和棉花等共計12種),分別以盤拍法隨機5點進行樣本采集。采集的所有薊馬保存于含有99.7%乙醇的離心管中;之后,于顯微鏡下確定薊馬的種類(主要依據(jù)成蟲的形態(tài)特征)[26],并統(tǒng)計西花薊馬的發(fā)生區(qū)域、寄主植物種類以及種群數(shù)量。然后以西花薊馬數(shù)量超過10頭的寄主植物種群為材料,包括蜀葵Althaea rosea (Linn.) Cavan.、月季Rosa chinensis Jacq.、美人蕉Canna indica Linn.、辣椒Capsicum annuum Linn.、茄子Solanum melongena Linn.等,進行種群遺傳多樣性分析。檢測對象為雌性個體,具體的地域范圍為東經(jīng)98°18′42.25″~106°16′07.29″,北緯37°56′25.55″~39°43′23.21″,相關(guān)信息詳見表1。

    1.2 西花薊馬種群基因組DNA提取

    以5種寄主植物上的西花薊馬種群為材料,參照喬瑋娜等的方法[27],提取單頭薊馬的DNA。提取的DNA以20 μL超純水復(fù)溶,濃度為100~250 ng/μL (超微量分光光度計,NanoPhotometerTM P330, Implen, Munich, Germany),-20℃保存?zhèn)溆谩?/p>

    1.3 引物設(shè)計與PCR擴增及序列測定

    參照NCBI數(shù)據(jù)庫中西花薊馬線粒體基因組DNA的全長序列(GenBank 登錄號:JN835456.1),以CO Ⅱ區(qū)域為靶標(biāo),利用Oligo 7軟件自行設(shè)計引物FOCOIIZTF1 (5′TTT TGT ACG AAT CTC AAC A3′)和FOCOIIZTR1 (5′CTA GGT CGA AAC TAA ATG C3′)。引物由生工生物工程(上海)股份有限公司合成;擴增片段大小約為650 bp。CO Ⅰ基因序列的擴增使用標(biāo)準(zhǔn)的DNA條形碼引物L(fēng)epF1 (5′ATT CAA CCA ATC ATA AAG ATA TTG G3′)和LepR1 (5′TAA ACT TCT GGA TGT CCA AAA AAT CA3′),擴增片段大小約為700 bp[28]。PCR反應(yīng)體系為25 μL,其中,模板DNA 2 μL,10× Buffer 2.5 μL (含Mg2+),dNTPs 0.5 μL (0.2 mmol/L),上游引物和下游引物各0.5 μL(5 pmol/L),Taq DNA聚合酶(5 U/μL)0.25 μL,超純水18.75 μL。擴增程序:94℃預(yù)變性5 min;35個循環(huán):94℃ 30 s,51℃ (CO Ⅰ) / 55℃ (CO Ⅱ) 30 s,72℃ 60 s;最后72℃延伸5 min。

    取3 μL PCR擴增產(chǎn)物,加1 μL上樣緩沖液 (025%溴酚藍,40%蔗糖水溶液),在1% (W∶V)瓊脂糖凝膠(含0.5 μg/mL的Goldview Ⅱ)上以100 V電泳(電泳緩沖液為0.5×TAE)分離50 min后,以GelDoc Universal Hood Ⅱ型凝膠成像系統(tǒng)分析電泳結(jié)果,將含有目的片段的PCR產(chǎn)物送生工生物工程(上海)股份有限公司進行雙向測序。每一寄主種群檢測的數(shù)量,依據(jù)寄主植物的田間種植情況、西花薊馬的發(fā)生和采集情況以及測序質(zhì)量而定,各種群的檢測數(shù)量為14~63頭;5種寄主種群共計檢測獲得CO Ⅰ基因序列174條,CO Ⅱ基因序列177條(表1)。

    1.4 序列比對分析

    1.4.1 序列檢驗與分析

    以MEGA 7.0軟件中的TraceEditor程序檢驗序列的峰圖質(zhì)量,并進行序列閱讀和人工校閱,對峰圖質(zhì)量較差或出現(xiàn)套峰的序列直接棄用,并重新擴增、測序和檢驗。然后,以Clustal W軟件進行多重序列比對分析,并利用NCBI數(shù)據(jù)庫進行BLAST比對,以確保所得序列來自靶標(biāo)基因。序列特征分析,包括多態(tài)位點分析、轉(zhuǎn)換/顛換偏倚率、堿基組成等也在MEGA7.0軟件環(huán)境下完成。

    1.4.2 遺傳學(xué)參數(shù)計算與分析

    以DnaSP 5.10軟件計算核苷酸多樣性(π)、單倍型多樣性(haplotype diversity, Hd),以及種群間遺傳分化程度評價參數(shù)FST(固定指數(shù))、基因流及核苷酸平均差異數(shù)(K)等分子遺傳學(xué)參數(shù)[29],并分別對mtDNA CO Ⅰ和CO Ⅱ基因片段進行Tajimas D和Fus Fs中性檢驗[3031]。同時,利用Arlequin 3.5.2.2軟件進行分子變異分析(analysis of molecular variance, AMOVA)[32]。

    1.4.3 網(wǎng)狀樹和系統(tǒng)發(fā)育樹構(gòu)建

    以不同寄主植物上的西花薊馬種群單倍型為靶標(biāo),應(yīng)用Network 5.0程序中的中接網(wǎng)絡(luò)(medianjoining networks, MJ)法構(gòu)建單倍型網(wǎng)狀樹[33]。同時,以花薊馬Frankliniella intonsa (GenBank:JQ917403)為外群,采用MEGA 7.0中的鄰接法(neighborjoining, NJ)和最大似然法(maximum likelihood, ML)分別構(gòu)建系統(tǒng)發(fā)育樹,估算西花薊馬各單倍型之間的系統(tǒng)發(fā)育關(guān)系。其中,NJ樹主要通過MEGA 7.0中的Kimura雙參數(shù)模型(Kimura2Parameter)進行構(gòu)建[34],并對進化樹的各分支進行Bootstrap檢測,重復(fù)1 000次;而ML樹的構(gòu)建,首先由MEGA 7.0計算獲得最優(yōu)核苷酸替代模型(Tamura3parameter model),然后再在此基礎(chǔ)上建樹,各分支同樣進行Bootstrap檢測,重復(fù)1 000次。此外,以CO Ⅰ基因為靶標(biāo)構(gòu)建進化樹時,分別以數(shù)據(jù)庫中已公開的西花薊馬溫室品系(GenBank:JN790696)和羽扇豆品系(GenBank:JN790699)的相應(yīng)堿基序列作為陽性參照[11]。

    2 結(jié)果與分析

    2.1 西花薊馬種群mtDNA CO Ⅰ和CO Ⅱ基因片段的序列特征

    CO Ⅰ基因經(jīng)雙向拼接和Clustal W序列比對分析后,得到的片段長度為656 bp,共檢測到30個多態(tài)位點,包括6個自裔位點和24個簡約信息位點;其中,發(fā)生轉(zhuǎn)換的位點數(shù)為7個,發(fā)生顛換的位點數(shù)為1個,轉(zhuǎn)換/顛換偏倚率R為5.3;核苷酸A、T、C和G的平均含量分別為31.1%、39.0%、15.5% 和14.5%,堿基組成具有明顯的AT偏好性,(A+T)含量為70.1%。CO Ⅱ基因經(jīng)雙向拼接和Clustal W序列比對后,得到的片段長度為608 bp,共檢測到26個多態(tài)位點,均為簡約信息位點;其中,發(fā)生轉(zhuǎn)換的位點數(shù)為8個,發(fā)生顛換的位點數(shù)為0,轉(zhuǎn)換/顛換偏倚率R為27.3;核苷酸A、T、C、G的平均含量依次為34.8%、38.7%、12.9%和13.6%,堿基組成具有明顯的AT偏好性,(A+T)含量為73.5%。

    2.2 不同寄主植物上的西花薊馬種群mtDNA單倍型分布

    當(dāng)以CO Ⅰ基因為靶標(biāo)時,在測序成功的174條序列中共檢測到13種單倍型(haplotype, H)(H_1~H_13)。其中,H_1和H_2為優(yōu)勢單倍型,占單倍型的71.26% (表2)。此外,辣椒和茄子上的西花薊馬種群只有3種單倍型,且都為溫室品系單倍型,其中嘉峪關(guān)辣椒(JYGLJ)上的種群單倍型最少,只有2種;而蜀葵、月季和美人蕉上西花薊馬種群的單倍型多達7~9種,涉及溫室品系和羽扇豆品系的所有單倍型,其中張掖蜀葵(ZYSK)上的種群單倍型最少,僅有3種(表2)。當(dāng)以CO Ⅱ基因為靶標(biāo)時,在測序成功的177條序列中,共檢測到12種單倍型(H_1~H_12);其中,H_1和H_2為優(yōu)勢單倍型,占所有單倍型的70.62% (表3)。此外,辣椒和茄子上的種群均只有2種單倍型,都為西花薊馬溫室品系;而蜀葵、月季和美人蕉上西花薊馬種群的單倍型數(shù)量為8~9種,同樣包括了溫室品系和羽扇豆品系的所有單倍型,其中嘉峪關(guān)蜀葵(JYGSK)上的種群單倍型最少,只有2種(表3)。

    綜合分析,辣椒和茄子上的種群單倍型種類相對較少,主要以兩種最常見的單倍型H_1和H_2為優(yōu)勢單倍型,分別占96.96%和86.67% (CO Ⅰ基因)以及100%和100% (CO Ⅱ基因);而蜀葵、月季和美人蕉上的種群單倍型種類相對較多,幾乎所有的單倍型都有發(fā)現(xiàn),其中以代表溫室品系的單倍型H_1和H_2以及代表羽扇豆品系的單倍型H_6 (CO Ⅰ基因)或H_5 (CO Ⅱ基因)為優(yōu)勢單倍型(表2~3)。

    2.3 不同寄主植物上的西花薊馬種群mtDNA CO Ⅰ和CO Ⅱ基因遺傳多樣性分析

    遺傳多樣性分析結(jié)果顯示,總體上,無論以mtDNA CO Ⅰ基因還是CO Ⅱ基因為靶標(biāo),辣椒和茄子上的西花薊馬種群的單倍型多樣性、核苷酸多樣性以及核苷酸平均差異數(shù)等遺傳多樣性指數(shù)相對較低;而蜀葵、月季及美人蕉上的種群單倍型多樣性、核苷酸多樣性以及核苷酸平均差異數(shù)等遺傳多樣性指數(shù)則相對較高,但以CO Ⅰ基因為靶標(biāo)時的張掖蜀葵(ZYSK)上的種群和以CO Ⅱ基因為靶標(biāo)時的嘉峪關(guān)蜀葵(JYGSK)上的種群的遺傳多樣性指數(shù)偏低(表4)。CO Ⅰ和CO Ⅱ基因的總體寄主植物種群的Fus Fs中性檢驗結(jié)果分別為2.36和4.06,且無顯著性,表明西花薊馬5種寄主植物上的種群歷史上未經(jīng)歷明顯的擴張,其種群大小整體保持穩(wěn)定。此外,所有寄主植物上的種群的Tajimas D和Fus Fs中性檢驗結(jié)果均未達到顯著水平(表4)。

    2.4 不同寄主植物上的西花薊馬種群遺傳分化與分子變異分析

    當(dāng)以CO Ⅰ基因為靶標(biāo)時,西花薊馬總寄主種群固定指數(shù)FST為0.068,變異范圍為-0.046~0.262;總基因流Nm為6.16,變異范圍為1.41~26.03。當(dāng)以CO Ⅱ基因為靶標(biāo)時,西花薊馬總寄主種群固定指數(shù)FST為0.067,變異范圍為-0.051~0.106;總基因流Nm為8.24,變異范圍為4.21~40.74 (表5)。綜合分析顯示,本研究中5種常見寄主植物上的西花薊馬種群之間基因交流充分,尚未發(fā)生明顯的遺傳分化。但兩種基因反映的各種群間遺傳分化水平并不完全一致,如基于CO Ⅰ基因的分析顯示,張掖茄子上的種群(ZYQZ)與嘉峪關(guān)辣椒上的種群(JYGLJ)和蜀葵上的種群(JYGSK)以及銀川美人蕉上的種群(YCMRJ)之間,張掖辣椒上的種群(ZYLJ)與張掖蜀葵上的種群(ZYSK)、銀川美人蕉上的種群(YCMRJ)之間存在明顯的遺傳分化;而基于CO Ⅱ基因的分析則顯示,它們之間的遺傳分化并不顯著(表5)。進一步的AMOVA分子變異分析表明,基于CO Ⅰ基因的各種群內(nèi)的遺傳變異(FST)差異性顯著(P < 0.05),而基于CO Ⅱ基因上則無顯著性差異(P > 0.05),但兩基因的分析均顯示種群遺傳變異主要源自種群內(nèi)部;此外,不同寄主植物間(FCT)以及同一寄主不同地理種群間(FSC)均表現(xiàn)出較低的遺傳變異水平(P > 0.05),各寄主植物種群間的遺傳變異不明顯(表6)。

    2.5 單倍型系統(tǒng)發(fā)育及品系分析

    以花薊馬F.intonsa (GenBank:JN835456.1)為外群,采用ML法和NJ法分別構(gòu)建各寄主植物上西花薊馬種群不同單倍型的系統(tǒng)發(fā)育樹,并對系統(tǒng)發(fā)育樹顯示的拓?fù)潢P(guān)系進行分析。結(jié)果顯示,當(dāng)以CO Ⅰ基因為靶標(biāo)時,兩種方法構(gòu)建的系統(tǒng)進化樹的聚類結(jié)果基本一致,且所有的單倍型均與外群明顯區(qū)分。整體上,進化樹分為明顯的2支。其中,第一支包含H_1~H_5, H_7~H_11以及H_13等共計11種單倍型,對應(yīng)西花薊馬的溫室品系;第二支包含單倍型H_6和H_12,對應(yīng)西花薊馬的羽扇豆品系(圖1a, 2a)。進一步的單倍型寄主植物對應(yīng)關(guān)系分析表明,代表羽扇豆品系分支的單倍型H_6和H_12來自寄主植物蜀葵、月季和美人蕉,而代表溫室品系分支的單倍型則在所有的寄主植物上均有發(fā)生(表2; 圖1a, 2a)。當(dāng)以CO Ⅱ基因為靶標(biāo)時,兩種方法構(gòu)建的系統(tǒng)發(fā)育樹的聚類結(jié)果基本一致。整體顯示,系統(tǒng)發(fā)育樹同樣分為2支,其中,第一支包含H_1~H_4,H_6~H_7和H_9~H_12等10種單倍型;第二支包含2種單倍型,即單倍型H_5和H_8;兩個分支分別代表西花薊馬的溫室品系和羽扇豆品系(圖1b, 2b)。并且,單倍型寄主植物對應(yīng)關(guān)系的分析結(jié)果亦與CO Ⅰ基因的一致;亦即,代表羽扇豆品系分支的單倍型H_5和H_8來自寄主植物蜀葵、月季和美人蕉,而溫室品系的各單倍型則在所有寄主植物上都有發(fā)生(表3,圖1b, 2b)。

    此外,以中接法構(gòu)建的單倍型網(wǎng)狀樹顯示的聚類分析結(jié)果,也與系統(tǒng)發(fā)育樹的單倍型分析結(jié)果基本一致,亦即,西花薊馬所有的種群個體聚為2支,并分別對應(yīng)西花薊馬的溫室品系和羽扇豆品系(圖3)。當(dāng)以CO Ⅰ基因為靶標(biāo)時,兩個種群的比例分別為91.38%和8.62% (表2, 圖3a);當(dāng)以CO Ⅱ基因為靶標(biāo)時,兩個種群的比例分別為90.96%和9.04% (表3, 圖3b)。

    3 討論

    遺傳多樣性是生物多樣性的重要組成部分,也是維持物種進化的根本保證[35]。就入侵生物而言,其遺傳多樣性極易受到“奠基者效應(yīng)”、“瓶頸效應(yīng)”、遺傳漂變等因素的影響,并進而使入侵種群的遺傳多樣性降低[3637]。西花薊馬是本世紀(jì)初傳入我國的一種重要的外來入侵生物,為闡明其傳播擴散機制我國學(xué)者相繼開展了基因變異、群體遺傳結(jié)構(gòu)分化等方面的研究,并證實入侵我國的西花薊馬種群的遺傳多樣性較低,種群歷史上可能經(jīng)歷過“瓶頸效應(yīng)”[11, 13]。本研究以mtDNA CO Ⅰ和CO Ⅱ基因為靶標(biāo),以我國寧夏和甘肅田間常見的5種寄主植物上的西花薊馬為對象,分析了其種群遺傳多樣性和遺傳結(jié)構(gòu),結(jié)果顯示,辣椒和茄子上的西花薊馬種群其單倍型多樣性、核苷酸多樣性以及核苷酸平均差異數(shù)等遺傳參數(shù)總體上均較蜀葵、月季和美人蕉等寄主種群低(表4),結(jié)合采樣地特征分析其原因可能與辣椒和茄子受人為干擾較多有關(guān);進一步的種群間固定指數(shù)和基因流分析顯示,5種常見寄主植物上西花薊馬種群之間基因交流充分,未發(fā)生明顯的遺傳分化,但各種群間兩種基因反映的遺傳分化水平并非完全一致(表5)。有研究表明不同線粒體基因以及蛋白質(zhì)編碼基因的不同編碼位點的進化速率存在差異[3839],因此推測這種不一致可能與西花薊馬mtDNA CO Ⅰ和CO Ⅱ基因的突變速率不同有關(guān)。此外,研究還發(fā)現(xiàn),不同寄主植物上的西花薊馬種群的多態(tài)性位點數(shù)差異不明顯,且沒有明顯的規(guī)律性;同時,檢測出的CO Ⅰ和CO Ⅱ基因的單倍型數(shù)量均較少,分別只有13種和12種,低于原產(chǎn)地的30種(CO Ⅰ基因)[25],表明該區(qū)域西花薊馬種群可能經(jīng)歷過“瓶頸效應(yīng)”,或是由于初始定殖種群數(shù)量較少而顯示了“奠基者效應(yīng)”。

    Tajimas D和Fus Fs模型是目前兩種應(yīng)用最為廣泛的種群內(nèi)等位基因頻率中性檢驗?zāi)P?,該兩種模型可用于推算種群的歷史變化動態(tài)[40]。當(dāng)中性檢驗結(jié)果顯著小于0時,可推斷種群在歷史上曾發(fā)生過規(guī)模擴張和定向選擇;當(dāng)中性檢驗結(jié)果顯著大于0時,可認(rèn)為種群大小在歷史上曾發(fā)生萎縮和平衡選擇;而當(dāng)D值不顯著背離0時,中性零假說則不能被排除,表明種群大小歷史上保持相對穩(wěn)定[3031]。本研究的Tajimas D和Fus Fs中性檢驗結(jié)果均不具有顯著性,表明我國西北地區(qū)的寧夏和甘肅田間5種寄主植物上的西花薊馬種群歷史上并未經(jīng)歷明顯的擴張,究其原因可能與其入侵時間比較短,尚未積累到足夠的遺傳突變而偏離中性進化有關(guān)(表4)。進一步的系統(tǒng)發(fā)育分析顯示,與溫室品系不同,代表羽扇豆品系的單倍型主要來自多年生的蜀葵、月季和美人蕉等3種植物,而辣椒和茄子并無羽扇豆品系的發(fā)生,表明野外不同品系的西花薊馬其寄主偏好性可能不一樣,羽扇豆品系可能更嗜食多年生植物。此結(jié)論在以往的研究中亦得到了佐證,如,在RugmanJones等的研究中[41],除未確定種名的寄主植物外,羽扇豆品系均采自桃Prunus persica、歐洲李Prunus domestica、柳葉石楠Heteromeles arbutifolia、巨杉Sequoiadendron giganteum、加州葡萄Vitis californica、灌木羽扇豆Lupinus arboreus、花菱草Eschscholzia californica、全緣葉美洲茶Ceanothus integerrimus等多年生植物;在我國的相關(guān)研究中,羽扇豆品系亦多采自月季、三葉草Trifolium repens、繡球花Hydrangea macrophylla、洋薊Cynara scolymus、油菜Brassica campestris,以及多花薔薇Rosa multiflora、大麗花Dahlia pinnata和雞冠花Celosia cristata等植物[1214],這些寄主植物亦多為二年生、多年生或越年生。故此推測,多年生/越年生植物可能更有利于西花薊馬羽扇豆品系種群的發(fā)生。

    本研究對象為西北地區(qū)寧夏、甘肅范圍內(nèi)常見的5種寄主植物上的西花薊馬種群,寄主種類相對較少,而西花薊馬屬雜食性害蟲,寄主范圍廣泛,因此有關(guān)其寄主專化性以及種群遺傳多樣性尚需進一步研究,未來應(yīng)在更廣泛的地理范圍開展樣本采集,并需要結(jié)合生物學(xué)、生態(tài)學(xué)等特性予以佐證。

    參考文獻

    [1] KIRK W D J. The pest and vector from the west: Frankliniella occidentalis [C]∥Thrips and Tospoviruses: Proceedings of the 7th International Symposium on Thysanoptera, Canberra, Australia:Australian National Insect Collection,2002:3342.

    [2] 呂要斌, 貝亞維, 林文彩, 等. 西花薊馬的生物學(xué)特性、寄主范圍及危害特點[J]. 浙江農(nóng)業(yè)學(xué)報, 2004, 16(5): 317320.

    [3] KIRK W D J, TERRY L I. The spread of the western flower thrips Frankliniella occidentalis (Pergande)[J]. Agricultural and Forest Entomology, 2003, 5(4): 301310.

    [4] 張友軍,吳青君,徐寶云,等.危險性外來入侵生物——西花薊馬在北京發(fā)生危害[J].植物保護,2003,29(4):5859.

    [5] 吳青君,徐寶云,張治軍,等.京、浙、滇地區(qū)植物薊馬種類及其分布調(diào)查[J].中國植保導(dǎo)刊,2007,27(1):3234.

    [6] 鄭長英,劉云虹,張乃芹,等.山東省發(fā)現(xiàn)外來入侵有害生物——西花薊馬[J].青島農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版),2007,24(3):172174.

    [7] 袁成明,郅軍銳,李景柱,等.貴州省蔬菜薊馬種類調(diào)查研究[J].中國植保導(dǎo)刊,2008,28(7):810.

    [8] 劉佳,張林,盧焰梅,等.湖南外來入侵害蟲西花薊馬初步調(diào)查[J].安徽農(nóng)業(yè)科學(xué),2010,38(25):1380013801.

    [9] 楊華,崔元玗,張升,等.外來入侵害蟲西花薊馬在新疆的發(fā)生為害[J].新疆農(nóng)業(yè)科學(xué),2010,47(11):22522253.

    [10]呂要斌,張治軍,吳青君,等.外來入侵害蟲西花薊馬防控技術(shù)研究與示范[J].應(yīng)用昆蟲學(xué)報,2011,48(3):488496.

    [11]YANG Xianming, SUN Jingtao, XUE Xiaofeng, et al. Invasion genetics of the western flower thrips in China: evidence for genetic bottleneck, hybridization and bridgehead effect [J]. PLoS ONE, 2012, 7: e34567.

    [12]DUAN Huisheng, YU Yi, ZHANG Ansheng, et al. Genetic diversity and inferences on potential source areas of adventive Frankliniella occidentalis (Thysanoptera: Thripidae) in Shandong, China based on mitochondrial and microsatellite markers[J]. Florida Entomologist, 2013, 96(3): 933940.

    [13]沈登榮,宋文菲,袁盛勇,等.基于mtDNACOⅠ的云南西花薊馬的遺傳分析[J].植物保護,2014,40(5):7579.

    [14]張桂芬,喬瑋娜,古君伶,等.我國西花薊馬線粒體DNACOⅠ基因變異及群體遺傳結(jié)構(gòu)分析[J].生物安全學(xué)報.2014,23(3):196209.

    [15]LOOMANS A J M, LENTEREN J C V. Biological control of thrips pests: a review on thrips parasitoids [J]. Wageningen Agricultural University Papers, 1995, 95(1): 89201.

    [16]MORITZ G. The biology of thrips is not the biology of their adults: a developmental view [C]∥Thrips and Tospoviruses: Proceedings of the 7th International Symposium on Thysanoptera, Canberra, Australia: Australian National Insect Collection, 2002: 259267.

    [17]BRDSGAARD H F. Insecticide resistance in European and African strains of western flower thrips (Thysanoptera: Thripidae) tested in a new residueonglass test [J]. Journal of Economic Entomology, 1994, 87(5): 11411146.

    [18]NIELSEN M C, TEULON D A J, CHAPMAN R B, et al. Comparison of life history parameters of two Frankliniella occidentalis (Thysanoptera: Thripidae) strains in New Zealand [J]. Environmental Entomology, 2010, 39(2): 303311.

    [19]JANZ N, NYBLOM K, NYLIN S. Evolutionary dynamics of hostplant specialization: a case study of the tribe Nymphalini [J]. Evolution, 2001, 55(4): 783796.

    [20]NOVOTNY V, DROZD P, MILLER S E, et al. Why are there so many species of herbivorous insects in tropical rainforests?[J]. Science, 2006, 313(5790): 11151118.

    [21]BRUNNER P C, CHATZIVASSILIOU E K, KATIS N I, et al. Hostassociated genetic differentiation in Thrips tabaci (Insecta; Thysanoptera), as determined from mtDNA sequence data [J]. Heredity, 2004, 93(4): 364370.

    [22]JURADORIVERA J A, VOGLER A P, REID C A, et al. DNA barcoding insecthost plant associations [J]. Proceedings of the Royal Society of London B: Biological Sciences, 2009, 276(1657): 639648.

    [23]MOPPER S.Adaptive genetic structure in phytophagous insect populations [J]. Trends in Ecology and Evolution, 1996, 11(6): 235238.

    [24]PYRY J, LUOTO M, PAUKKUNEN J, et al. Different responses of plants and herbivore insects to a gradient of vegetation height: an indicator of the vertebrate grazing intensity and successional age [J]. OIKOS, 2006, 115(3): 401412.

    [25]BRUNNER P C, FREY J E.Habitatspecific population structure in native western flower thrips Frankliniella occidentalis (Insecta, Thysanoptera)[J]. Journal of Evolutionary Biology, 2010, 23(4): 797804.

    [26]FUNDERBURK J, DIFFIE S, SHARMA J, et al. Thrips of ornamentals in the southeastern US [J]. Entomology and Nematology Department, Florida Cooperative Extension Service, Institute of Food and Agricultural Sciences, University of Florida, 2007, ENY845 (IN754).

    [27]喬瑋娜, 萬方浩, 張愛兵, 等. DNA條形碼技術(shù)在田間常見薊馬種類識別中的應(yīng)用[J]. 昆蟲學(xué)報, 2012, 55(3): 344356.

    [28]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 butterfly Astraptes fulgerator [J]. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(41): 1481214817.

    [29]LIBRADO P, ROZAS J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J]. Bioinformatics, 2009, 25(11): 14511452.

    [30]TAJIMA F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism[J]. Genetics, 1989, 123(3): 585595.

    [31]FU Yunxin. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection[J]. Genetics, 1997, 147(2): 915925.

    [32]EXCOFFIER L, LISCHER H E.Arlequin suite ver3.5: a new series of programs to perform population genetics analyses under Linux and Windows [J]. Molecular Ecology Resources, 2010, 10(3): 564567.

    [33]BANDELT H J, FORSTER P, RHL A. Medianjoining networks for inferring intraspecific phylogenies [J]. Molecular Biology and Evolution, 1999, 16(1): 3748.

    [34]KUMAR S, STECHER G, TAMURA K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets [J]. Molecular Biology and Evolution, 2016, 33(7): 18701874.

    [35]GIENAPP P, TEPLITSKY C, ALHO J S, et al. Climate change and evolution: disentangling environmental and genetic responses [J]. Molecular Ecology, 2008, 17(1): 167178.

    [36]AUSTERLITZ F, JUNGMULLER B, GODELLE B, et al. Evolution of coalescence times, genetic diversity and structure during colonization [J]. Theoretical Population Biology, 1997, 51(2): 148164.

    [37]TSUCHIDA K, KUD K, ISHIGURO N. Genetic structure of an introduced paper wasp, Polistes chinensis antennalis (Hymenoptera, Vespidae) in New Zealand [J]. Molecular Ecology, 2014, 23(16): 40184034.

    [38]HUA Jimeng, LI Ming, DONG Pengzhi, et al. Comparative and phylogenomic studies on the mitochondrial genomes of Pentatomomorpha (Insecta: Hemiptera: Heteroptera)[J]. BMC Genomics, 2008, 9(1): 610.

    [39]NEGRISOLO E, BABBUCCI M, PATARNELLO T. The mitochondrial genome of the ascalaphid owlfly Libelloides macaronius and comparative evolutionary mitochondriomics of neuropterid insects [J]. BMC Genomics, 2011, 12(1): 221.

    [40]NIELSEN R. Molecular signatures of natural selection[J]. Annual Review of Genetics, 2005, 39: 197218.

    [41]RUGMANJONES P F, HODDLE M S, STOUTHAMER R. Nuclearmitochondrial barcoding exposes the global pest western flower thrips (Thysanoptera: Thripidae) as two sympatric cryptic species in its native California [J]. Journal of Economic Entomology, 2010, 103(3): 877886.

    (責(zé)任編輯:田 喆)

    猜你喜歡
    薊馬靶標(biāo)品系
    薊馬防治難的原因及綜防措施
    靶標(biāo)評改,讓習(xí)作評改有序更有效
    昆明地區(qū)危害蘋果的薊馬發(fā)生初步調(diào)查
    甬優(yōu)系列不同組合在江蘇省淮北地區(qū)試種表現(xiàn)初報
    用DTOPSIS法綜合評價陸地棉23個品系
    水稻品種(苗頭品系)對比試驗
    如何防治大棚茄子蔥薊馬
    新殺菌劑靶標(biāo)和先導(dǎo)化合物的探索
    對捕食螨在薊馬生物防治中的應(yīng)用研究
    低鹽脅迫下條斑紫菜不同品系體細(xì)胞與殼孢子的耐受性差異
    毛片一级片免费看久久久久| av福利片在线观看| 亚洲欧美精品专区久久| 狂野欧美激情性xxxx在线观看| 国产在视频线精品| 中文在线观看免费www的网站| 黄色一级大片看看| 免费播放大片免费观看视频在线观看| 国产黄片视频在线免费观看| 97精品久久久久久久久久精品| 国产色婷婷99| 国产精品久久久久久久电影| 99久久中文字幕三级久久日本| 免费观看无遮挡的男女| 欧美成人午夜免费资源| 色网站视频免费| 看非洲黑人一级黄片| 免费大片18禁| 99久国产av精品| 内射极品少妇av片p| 亚洲精品日韩在线中文字幕| or卡值多少钱| 免费看av在线观看网站| 天天一区二区日本电影三级| 日韩中字成人| 午夜福利在线观看吧| 国产91av在线免费观看| av福利片在线观看| 亚洲av免费在线观看| 一区二区三区乱码不卡18| 高清毛片免费看| 少妇熟女aⅴ在线视频| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 成人av在线播放网站| 蜜桃亚洲精品一区二区三区| 中文字幕亚洲精品专区| 精品亚洲乱码少妇综合久久| 国产高清三级在线| 成人午夜精彩视频在线观看| 日韩成人伦理影院| 国产黄片美女视频| 两个人的视频大全免费| 街头女战士在线观看网站| 国产精品av视频在线免费观看| 亚洲国产欧美在线一区| 激情 狠狠 欧美| 91久久精品国产一区二区成人| 午夜福利视频1000在线观看| 国产乱来视频区| 伊人久久国产一区二区| 80岁老熟妇乱子伦牲交| 国产精品国产三级专区第一集| 男女边吃奶边做爰视频| 国产精品美女特级片免费视频播放器| 成人午夜精彩视频在线观看| 国产真实伦视频高清在线观看| av在线天堂中文字幕| 日韩一区二区三区影片| 特级一级黄色大片| 97人妻精品一区二区三区麻豆| 午夜亚洲福利在线播放| 丝瓜视频免费看黄片| 三级经典国产精品| 国产精品人妻久久久影院| eeuss影院久久| 视频中文字幕在线观看| 亚洲成人一二三区av| 国产成人一区二区在线| 97人妻精品一区二区三区麻豆| 国产av国产精品国产| 极品教师在线视频| 精品国产一区二区三区久久久樱花 | 九色成人免费人妻av| 在现免费观看毛片| 欧美一级a爱片免费观看看| 精品熟女少妇av免费看| 国产男人的电影天堂91| 久久久久性生活片| 麻豆成人av视频| 国产真实伦视频高清在线观看| 最近中文字幕高清免费大全6| 99久国产av精品| 少妇熟女aⅴ在线视频| 在线观看一区二区三区| 亚洲四区av| 99热这里只有是精品50| 亚洲精品亚洲一区二区| av在线播放精品| 蜜桃久久精品国产亚洲av| 国产成人精品福利久久| 亚洲国产精品专区欧美| 欧美日韩国产mv在线观看视频 | 超碰97精品在线观看| 亚洲在久久综合| 3wmmmm亚洲av在线观看| 春色校园在线视频观看| 国产又色又爽无遮挡免| 免费高清在线观看视频在线观看| 亚洲国产精品成人综合色| 伊人久久国产一区二区| 精品一区二区三区人妻视频| 久久久精品欧美日韩精品| 国产单亲对白刺激| 大陆偷拍与自拍| 亚洲一级一片aⅴ在线观看| 日本欧美国产在线视频| 人妻一区二区av| 精华霜和精华液先用哪个| 国产免费视频播放在线视频 | 有码 亚洲区| 在线免费观看不下载黄p国产| 日韩在线高清观看一区二区三区| 久久这里只有精品中国| 国产片特级美女逼逼视频| 久久久精品免费免费高清| 日韩精品青青久久久久久| 欧美日韩视频高清一区二区三区二| 欧美 日韩 精品 国产| 日本熟妇午夜| 高清日韩中文字幕在线| 男人爽女人下面视频在线观看| 国产又色又爽无遮挡免| 男女下面进入的视频免费午夜| 国模一区二区三区四区视频| 国产成人精品一,二区| 日韩人妻高清精品专区| 成人午夜高清在线视频| 老女人水多毛片| 国产免费一级a男人的天堂| 丰满乱子伦码专区| 熟妇人妻不卡中文字幕| 我的老师免费观看完整版| or卡值多少钱| 中文字幕制服av| 人体艺术视频欧美日本| 国国产精品蜜臀av免费| 中文字幕亚洲精品专区| 国产一级毛片七仙女欲春2| 精品一区在线观看国产| 18禁动态无遮挡网站| 国产男女超爽视频在线观看| 有码 亚洲区| 免费看av在线观看网站| 精品国产露脸久久av麻豆 | 国产大屁股一区二区在线视频| 大话2 男鬼变身卡| 国产av国产精品国产| videos熟女内射| 天天躁夜夜躁狠狠久久av| 国产综合精华液| 直男gayav资源| 人妻夜夜爽99麻豆av| 国产精品爽爽va在线观看网站| 亚洲精品中文字幕在线视频 | 一夜夜www| 天堂影院成人在线观看| 男的添女的下面高潮视频| 在线观看美女被高潮喷水网站| 亚洲不卡免费看| 精品熟女少妇av免费看| 天美传媒精品一区二区| 十八禁国产超污无遮挡网站| 男女边摸边吃奶| 久久久精品欧美日韩精品| 在线 av 中文字幕| 国产乱来视频区| 在线免费观看不下载黄p国产| 亚洲国产成人一精品久久久| 极品教师在线视频| 人妻一区二区av| 午夜福利在线在线| 麻豆国产97在线/欧美| 免费黄色在线免费观看| 一级毛片aaaaaa免费看小| 肉色欧美久久久久久久蜜桃 | xxx大片免费视频| 高清欧美精品videossex| videos熟女内射| 在线a可以看的网站| 亚洲欧美成人综合另类久久久| 色播亚洲综合网| 特大巨黑吊av在线直播| www.av在线官网国产| 亚州av有码| 国产视频内射| 国产免费福利视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产大屁股一区二区在线视频| 亚洲精品亚洲一区二区| 国产亚洲av嫩草精品影院| kizo精华| 国产欧美另类精品又又久久亚洲欧美| 中文字幕亚洲精品专区| 成人毛片a级毛片在线播放| 久久97久久精品| 午夜精品一区二区三区免费看| 男女边摸边吃奶| 日韩大片免费观看网站| 国产v大片淫在线免费观看| 国产一级毛片七仙女欲春2| 综合色丁香网| 黄色欧美视频在线观看| 日本wwww免费看| 狂野欧美白嫩少妇大欣赏| 黄色配什么色好看| 一级a做视频免费观看| 久久久久久久大尺度免费视频| 人妻夜夜爽99麻豆av| 日日干狠狠操夜夜爽| 日韩视频在线欧美| 久久久久精品久久久久真实原创| 国产精品嫩草影院av在线观看| 亚洲性久久影院| 精品久久久久久电影网| 成人二区视频| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 亚洲无线观看免费| 麻豆国产97在线/欧美| 内地一区二区视频在线| 久久99蜜桃精品久久| 丰满少妇做爰视频| 成人高潮视频无遮挡免费网站| 少妇的逼好多水| 亚洲av不卡在线观看| 欧美精品一区二区大全| 亚洲国产最新在线播放| 久久久久国产网址| 中国国产av一级| 中文在线观看免费www的网站| 自拍偷自拍亚洲精品老妇| 18禁裸乳无遮挡免费网站照片| 亚洲国产日韩欧美精品在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人亚洲精品av一区二区| 一区二区三区高清视频在线| 精品久久久久久久末码| 亚洲欧美成人精品一区二区| 美女国产视频在线观看| 国产91av在线免费观看| 国产真实伦视频高清在线观看| 五月天丁香电影| 白带黄色成豆腐渣| 少妇裸体淫交视频免费看高清| 一级爰片在线观看| 精品久久久久久电影网| 亚洲av电影不卡..在线观看| 黄片无遮挡物在线观看| 日韩 亚洲 欧美在线| 国产黄片美女视频| 一级二级三级毛片免费看| 亚洲在线自拍视频| 久久久成人免费电影| 久久精品国产亚洲网站| 欧美精品一区二区大全| 午夜精品在线福利| 高清午夜精品一区二区三区| 久久韩国三级中文字幕| 免费看光身美女| 国产亚洲午夜精品一区二区久久 | 在线a可以看的网站| 一级毛片久久久久久久久女| 大又大粗又爽又黄少妇毛片口| 日韩不卡一区二区三区视频在线| 又爽又黄无遮挡网站| 中文天堂在线官网| 好男人在线观看高清免费视频| 伦理电影大哥的女人| 狂野欧美激情性xxxx在线观看| 国产精品麻豆人妻色哟哟久久 | 国产探花在线观看一区二区| 美女黄网站色视频| 高清在线视频一区二区三区| 成人国产麻豆网| 韩国av在线不卡| 国产一级毛片七仙女欲春2| 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| 淫秽高清视频在线观看| 国产在线男女| 国产成人福利小说| 午夜精品一区二区三区免费看| 99久久精品热视频| 免费黄网站久久成人精品| 国产亚洲5aaaaa淫片| 午夜精品在线福利| 天堂中文最新版在线下载 | 乱系列少妇在线播放| 插逼视频在线观看| 好男人视频免费观看在线| 国产精品99久久久久久久久| av又黄又爽大尺度在线免费看| 久久久久久久大尺度免费视频| 亚洲真实伦在线观看| 国产老妇伦熟女老妇高清| 波野结衣二区三区在线| 日韩av在线大香蕉| 国产精品一区二区三区四区免费观看| 日日干狠狠操夜夜爽| 亚洲欧美成人精品一区二区| or卡值多少钱| 日韩欧美一区视频在线观看 | 亚洲精品一二三| 777米奇影视久久| 日韩欧美国产在线观看| 日韩制服骚丝袜av| 久久久久久久久久久丰满| 亚洲成人中文字幕在线播放| 蜜桃亚洲精品一区二区三区| 中文欧美无线码| 午夜激情福利司机影院| 午夜日本视频在线| 国产精品美女特级片免费视频播放器| ponron亚洲| 2021少妇久久久久久久久久久| av天堂中文字幕网| 欧美日韩国产mv在线观看视频 | 久久久久久久久久成人| 欧美性猛交╳xxx乱大交人| 日韩一区二区视频免费看| 精品久久久久久电影网| 好男人在线观看高清免费视频| 国产激情偷乱视频一区二区| 麻豆久久精品国产亚洲av| 3wmmmm亚洲av在线观看| 三级毛片av免费| 国产午夜福利久久久久久| 噜噜噜噜噜久久久久久91| 国产在线一区二区三区精| 3wmmmm亚洲av在线观看| 成年av动漫网址| 国产黄色免费在线视频| 日韩一区二区三区影片| 人人妻人人澡欧美一区二区| 免费看av在线观看网站| 亚洲成人中文字幕在线播放| 午夜精品国产一区二区电影 | 久久久久久久久久黄片| 亚洲在线自拍视频| 久久久国产一区二区| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 国产亚洲av片在线观看秒播厂 | 国产精品久久久久久久久免| 哪个播放器可以免费观看大片| 久久99热6这里只有精品| 国产黄色免费在线视频| 99热这里只有是精品在线观看| 亚洲精品日韩在线中文字幕| 建设人人有责人人尽责人人享有的 | 男人狂女人下面高潮的视频| 国产精品无大码| 亚洲一区高清亚洲精品| 哪个播放器可以免费观看大片| 干丝袜人妻中文字幕| 久久99蜜桃精品久久| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 久久99热6这里只有精品| 久久99精品国语久久久| 亚洲一区高清亚洲精品| 婷婷色综合大香蕉| 国产91av在线免费观看| 免费av不卡在线播放| 色播亚洲综合网| 成年版毛片免费区| 卡戴珊不雅视频在线播放| 欧美不卡视频在线免费观看| 一区二区三区免费毛片| 亚洲国产高清在线一区二区三| 亚洲精品国产av成人精品| 精品国产一区二区三区久久久樱花 | 午夜福利在线观看免费完整高清在| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久人人人人人人| 99久久精品国产国产毛片| 亚洲电影在线观看av| 亚洲va在线va天堂va国产| 精品国产露脸久久av麻豆 | 国产黄片美女视频| 一本久久精品| 成人毛片60女人毛片免费| 一本一本综合久久| ponron亚洲| 最近的中文字幕免费完整| 国产乱人视频| 久久精品久久久久久久性| 亚洲av国产av综合av卡| 啦啦啦啦在线视频资源| 日韩伦理黄色片| 欧美成人午夜免费资源| 最近最新中文字幕免费大全7| 又黄又爽又刺激的免费视频.| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品久久久久久久人妻蜜臀av| 国产精品日韩av在线免费观看| 少妇的逼水好多| 噜噜噜噜噜久久久久久91| 91狼人影院| 国产黄片视频在线免费观看| 免费av不卡在线播放| 亚洲va在线va天堂va国产| 日本-黄色视频高清免费观看| 午夜精品在线福利| 日韩精品青青久久久久久| 麻豆精品久久久久久蜜桃| 国产精品无大码| 亚洲国产成人一精品久久久| 美女国产视频在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲欧洲日产国产| 国产成人一区二区在线| 91精品伊人久久大香线蕉| 最近中文字幕2019免费版| 免费观看的影片在线观看| 国产免费福利视频在线观看| 免费黄色在线免费观看| 日本午夜av视频| 欧美高清性xxxxhd video| 我的女老师完整版在线观看| 国产成人a区在线观看| 国产成人福利小说| 国产伦在线观看视频一区| 日本爱情动作片www.在线观看| 少妇高潮的动态图| 亚洲综合精品二区| 在线观看av片永久免费下载| 又爽又黄无遮挡网站| 日韩精品有码人妻一区| 在线观看人妻少妇| 欧美另类一区| 99re6热这里在线精品视频| 欧美高清成人免费视频www| 狂野欧美激情性xxxx在线观看| 欧美激情国产日韩精品一区| 六月丁香七月| 亚洲欧美一区二区三区黑人 | 国产精品av视频在线免费观看| 国产综合懂色| 麻豆国产97在线/欧美| 尾随美女入室| 六月丁香七月| 一区二区三区高清视频在线| 干丝袜人妻中文字幕| 26uuu在线亚洲综合色| 又爽又黄a免费视频| 久久人人爽人人爽人人片va| 人妻一区二区av| 色播亚洲综合网| av.在线天堂| 熟女人妻精品中文字幕| 亚洲国产精品sss在线观看| 别揉我奶头 嗯啊视频| av又黄又爽大尺度在线免费看| 久久久精品免费免费高清| 在线观看美女被高潮喷水网站| 网址你懂的国产日韩在线| 国产亚洲av嫩草精品影院| 日韩av不卡免费在线播放| 亚洲成人一二三区av| 日韩不卡一区二区三区视频在线| 日韩视频在线欧美| 超碰av人人做人人爽久久| 伦理电影大哥的女人| 国产亚洲精品久久久com| 日韩精品青青久久久久久| 欧美 日韩 精品 国产| 亚洲av.av天堂| 久久这里有精品视频免费| 青青草视频在线视频观看| 女人久久www免费人成看片| 亚洲四区av| 国产精品国产三级专区第一集| 亚洲国产成人一精品久久久| 亚洲自偷自拍三级| 国产成人一区二区在线| 一级毛片我不卡| 国产爱豆传媒在线观看| 99久国产av精品| 国产精品国产三级国产av玫瑰| 少妇的逼水好多| 女人十人毛片免费观看3o分钟| 美女脱内裤让男人舔精品视频| 国产美女午夜福利| 精品久久久久久久久亚洲| 又黄又爽又刺激的免费视频.| 综合色av麻豆| 三级经典国产精品| 亚洲欧美中文字幕日韩二区| 亚洲成人中文字幕在线播放| 国产成人精品一,二区| 亚洲av电影在线观看一区二区三区 | 天堂中文最新版在线下载 | 精品人妻一区二区三区麻豆| 丰满乱子伦码专区| h日本视频在线播放| av女优亚洲男人天堂| 国产 一区精品| 国产在视频线在精品| 久久久久久久午夜电影| 亚洲色图av天堂| 国产淫片久久久久久久久| 精品人妻视频免费看| 国产在视频线在精品| 国产亚洲最大av| 看十八女毛片水多多多| 国产人妻一区二区三区在| 久久久午夜欧美精品| 天堂俺去俺来也www色官网 | 嫩草影院精品99| 中文字幕制服av| 国产女主播在线喷水免费视频网站 | 成人一区二区视频在线观看| 九色成人免费人妻av| 久久亚洲国产成人精品v| 一区二区三区免费毛片| 日日啪夜夜撸| 午夜久久久久精精品| 激情 狠狠 欧美| av福利片在线观看| 麻豆乱淫一区二区| 七月丁香在线播放| 91狼人影院| 美女cb高潮喷水在线观看| 日本与韩国留学比较| 免费黄频网站在线观看国产| 欧美性感艳星| 免费观看在线日韩| 欧美3d第一页| 亚洲熟妇中文字幕五十中出| 久久鲁丝午夜福利片| 中文乱码字字幕精品一区二区三区 | 久99久视频精品免费| 在线观看一区二区三区| 国产激情偷乱视频一区二区| 美女内射精品一级片tv| 久久久亚洲精品成人影院| 亚洲av电影不卡..在线观看| 精品久久国产蜜桃| 亚洲av二区三区四区| 日韩 亚洲 欧美在线| 久久这里只有精品中国| 久久久久久久久中文| 亚洲国产最新在线播放| 2021天堂中文幕一二区在线观| 最近最新中文字幕大全电影3| 美女xxoo啪啪120秒动态图| 欧美日韩一区二区视频在线观看视频在线 | 大香蕉97超碰在线| 哪个播放器可以免费观看大片| 欧美97在线视频| 日日摸夜夜添夜夜爱| 在线观看av片永久免费下载| 色尼玛亚洲综合影院| 我要看日韩黄色一级片| av.在线天堂| 嘟嘟电影网在线观看| 久久人人爽人人爽人人片va| 国产精品一区二区三区四区免费观看| 最近2019中文字幕mv第一页| 草草在线视频免费看| 亚洲精华国产精华液的使用体验| 久久久久九九精品影院| 91aial.com中文字幕在线观看| 中文在线观看免费www的网站| 一级毛片久久久久久久久女| 久久久久久久大尺度免费视频| 有码 亚洲区| av天堂中文字幕网| 国产69精品久久久久777片| 久久久久久久久久黄片| 久久鲁丝午夜福利片| 精品人妻熟女av久视频| 亚洲伊人久久精品综合| 在线播放无遮挡| 日本欧美国产在线视频| 在线观看一区二区三区| 免费看不卡的av| 禁无遮挡网站| 欧美bdsm另类| 免费看不卡的av| 毛片女人毛片| 日韩中字成人| 18+在线观看网站| 毛片女人毛片| 日韩欧美精品v在线| 免费观看精品视频网站| 大片免费播放器 马上看| 只有这里有精品99| 国产成人精品一,二区| 在线观看免费高清a一片| 免费高清在线观看视频在线观看| 亚洲在线自拍视频| 日韩av在线大香蕉| 只有这里有精品99| 国产探花在线观看一区二区| 久久午夜福利片| 亚洲成人精品中文字幕电影| 大又大粗又爽又黄少妇毛片口| 大片免费播放器 马上看| 街头女战士在线观看网站| 欧美潮喷喷水| 国语对白做爰xxxⅹ性视频网站| 91av网一区二区| 国产午夜精品一二区理论片| 国产亚洲精品久久久com| 国内揄拍国产精品人妻在线| 欧美高清成人免费视频www| 亚洲av免费高清在线观看| 亚洲av成人精品一二三区| 麻豆成人av视频| 国产精品美女特级片免费视频播放器| 一级毛片 在线播放| 免费大片黄手机在线观看|