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

    我國茶毛蟲及其布尼亞病毒(EpBYV)的遺傳多樣性分析

    2024-11-08 00:00:00陳世春江宏燕廖姝然陳亭旭牛金志王曉慶
    茶葉科學(xué) 2024年5期

    摘要:茶毛蟲是一種危害茶葉生產(chǎn)的重要害蟲和人類皮膚致敏原,布尼亞病毒在茶毛蟲的不同地理種群中廣泛分布,研究茶毛蟲與茶毛蟲布尼亞病毒(Euproctis pseudoconspersa bunyavirus,EpBYV)遺傳背景可為更好地防治茶毛蟲和開發(fā)利用病毒資源提供理論支撐。采集了15個地理種群共148個茶毛蟲幼蟲樣本,測定了茶毛蟲的COI、ND5序列及EpBYV的RdRp序列,利用軟件DnaSP 6.12.03、Arlequin 3.5.2.2和MEGA 7.0.26等進行遺傳多樣性分析。合并COI、ND5序列分析發(fā)現(xiàn),茶毛蟲的15個地理種群具有較高的單倍型多樣性(h=0.880 68)和較低的核苷酸多樣性(π=0.003 17),99個種群對間遺傳分化顯著較高(FST>0.290,P<0.05),AMOVA分析顯示遺傳變異主要來自于種群之間(87.12%),且組間分化與我國地勢第二、三級階梯分界吻合,種群歷史動態(tài)分析推測茶毛蟲總?cè)后w在近期較為穩(wěn)定。EpBYV的RdRp序列在除重慶城口種群外的138個樣本中成功擴增,具有較高的單倍型多樣性(h=0.935 26)和相對較低的核苷酸多樣性(π=0.017 95),93個種群對間遺傳分化顯著較高(FST>0.257,P<0.05),遺傳變異主要來自于種群之間(62.13%),種群歷史動態(tài)分析支持EpBYV種群近期經(jīng)歷了種群擴張事件。綜合研究結(jié)果表明,我國茶毛蟲總?cè)后w在近期將較為穩(wěn)定,但在重慶城口和福建寧德有種群擴張風(fēng)險;EpBYV群體經(jīng)歷了種群擴張,EpBYV在茶毛蟲群體中的感染率和種群擴張能力都較高,具有應(yīng)用于茶毛蟲生物防治的潛力。

    關(guān)鍵詞:茶樹;茶毛蟲;布尼亞病毒;EpBYV;遺傳多樣性

    中圖分類號:S571.1;S435.711 文獻標識碼:A 文章編號:1000-369X(2024)05-793-14

    Genetic Diversity Analysis of Euproctis pseudoconspersa and Its Bunyavirus (EpBYV) in China

    CHEN Shichun1, JIANG Hongyan1, LIAO Shuran1, CHEN Tingxu1, NIU Jinzhi2, WANG Xiaoqing1*

    1. Tea Research Institute of Chongqing Academy of Agricultural Sciences, Chongqing 402160, China;

    2. College of Plant Protection, Southwest University, Chongqing 400715, China

    Abstract: Tea tussock moth, Euproctis pseudoconspersa, is an important pest which damages tea plants and causes human dermatitis. Euproctis pseudoconspersa bunyavirus, EpBYV, is a bunyavirus that widely distributed in different geographical populations of E. pseudoconspersa. In order to control the E. pseudoconspersa and utilize the virus resources, it is necessary to fully understand the genetic background of E. pseudoconspersa and EpBYV. In this study, 148 samples of E. pseudoconspersa larvae from 15 geographic populations were collected. COI and ND5 gene sequences of E. pseudoconspersa and RdRp sequences of EpBYV were determined. The genetic diversities of E. pseudoconspersa and EpBYV were analyzed by DnaSP 6.12.03, Arlequin 3.5.2.2 and MEGA 7.0.26. Sequence analysis of the combined sequences of COI and ND5 genes shows that 15 geographic populations have high haplotype diversity (h=0.880 68) and low nucleotide diversity (π=0.003 17). Significantly high genetic differentiation among 99 population pairs (FST>0.290, P<0.05) was identified. Molecular variance analysis (AMOVA) shows that the genetic differentiation of E. pseudoconspersa was mainly among populations (87.12%), and the differentiation among groups was consistent with the second and third ladder boundaries in China. Demographic history analysis suggests that the population of E. pseudoconspersa is relatively stable. RdRp sequences were successfully amplified in the 138 samples except CK population. RdRp sequence analysis reveals that the 14 geographic populations of EpBYV had high haplotype diversity (h=0.935 26) and relatively low nucleotide diversity (π=0.017 95). The 93 population pairs had significantly higher genetic differentiation (FST>0.257, P<0.05). AMOVA analysis shows that the genetic differentiation of EpBYV was mainly between populations (62.13%). Demographic history analysis reveals that EpBYV might have undergone population expansions in the past. Based on the analysis of this study, the population of E. pseudoconspersa in China is relatively stable, and there is a risk of population expansion in Chengkou, Chongqing and Ningde, Fujian. The population of EpBYV has experienced population expansion. The infection rate and population expansion ability of EpBYV in E. pseudoconspersa are high, which has good potential for biological control of E. pseudoconspersa.

    Keywords: tea plant, Euproctis pseudoconspersa, bunyavirus, EpBYV, genetic diversity

    茶毛蟲(Euproctis pseudoconspersa Strand),又名茶黃毒蛾,隸屬于鱗翅目(Lepidoptera)毒蛾科(Lymantriidae),是一種遍布我國各大茶區(qū)的重要害蟲,以幼蟲咀食葉片為害,嚴重時可致茶叢光禿,造成茶葉產(chǎn)量嚴重受損,且各蟲態(tài)的茶毛蟲毒毛觸及人體均易引起皮炎,嚴重影響茶園農(nóng)事活動及周邊居民健康[1-5]。昆蟲病毒廣泛應(yīng)用于鱗翅目害蟲的防治[6-7],目前茶毛蟲的病毒研究和應(yīng)用主要是茶毛蟲核型多角體病毒(Euproctis pseudoconspersa nuclear polyhedrosis virus,EpNPV)[8-11],其蘇云金桿菌(Bacillus thuringiensis,BT)混劑已實現(xiàn)了商品化[12-14]。EpNPV對茶毛蟲控制效果較好,但在田間調(diào)控存在滯后、株系間毒力差異[8],以及不同地理種群茶毛蟲對同株EpNPV敏感性不同等問題[15]。進一步收集并研究其他病毒,對茶毛蟲生物防治資源的儲備及生物防治效果的增強具有重要意義。茶毛蟲布尼亞病毒(Euproctis pseudoconspersa bunyavirus,EpBYV)屬于布尼亞病毒目(Bunyavirales)白纖病毒科(Phenuiviridae),是一種新鑒定的感染茶毛蟲的負鏈RNA病毒[16]。已有研究證實,布尼亞病毒目的費拉病毒(Ferak virus,F(xiàn)ERV)可通過蚊子產(chǎn)卵進行垂直傳播,并可能具有影響宿主傳播其他病毒的能力[17-19]。茶毛蟲感染EpBYV雖無明顯感染癥狀,但EpBYV可在茶毛蟲體內(nèi)復(fù)制并誘發(fā)體內(nèi)的RNAi抗病毒免疫反應(yīng)[15]。因此,研究EpBYV可為構(gòu)建基于RNA病毒的茶毛蟲綠色防控體系提供基礎(chǔ)數(shù)據(jù)。

    茶毛蟲作為一種全國茶區(qū)均有分布且在部分茶區(qū)屢有暴發(fā)的害蟲,近年來對其研究主要集中在田間防控方面[20-21]。對茶毛蟲開展遺傳多樣性研究,可更全面地從分子水平了解其種群遺傳結(jié)構(gòu),解析其生態(tài)適應(yīng)機制和成災(zāi)規(guī)律,為各茶區(qū)制定有效的防控策略提供基礎(chǔ)信息。此前,已有學(xué)者對茶小綠葉蟬[22]、茶網(wǎng)蝽[23]、茶棍薊馬[24]開展了相關(guān)研究。因此,本研究采集了我國四大茶區(qū)11個省份共15個地區(qū)茶園的茶毛蟲幼蟲樣本,擴增茶毛蟲的線粒體COI和ND5序列,以及EpBYV的RNA聚合酶(RNA-dependent RNA polymerase,RdRp)編碼基因序列,聯(lián)合分析茶毛蟲和EpBYV的種群遺傳多樣性和遺傳結(jié)構(gòu)特征,探討茶毛蟲種群歷史動態(tài),評估其近期在全國茶區(qū)的成災(zāi)風(fēng)險,明確EpBYV對茶毛蟲的感染率、遺傳變異和種群發(fā)展動態(tài)歷史,評估其未來用于生物防治的潛力。

    1 材料與方法

    1.1 試驗材料

    試驗所用茶毛蟲幼蟲樣本于2020年6月至2023年9月采集自我國四大茶區(qū)(江南茶區(qū)、江北茶區(qū)、西南茶區(qū)、華南茶區(qū))11個省份的15個茶園,具體信息如表1所示。

    1.2 試驗方法

    1.2.1 DNA、RNA提取及cDNA合成

    取PBS緩沖液清洗后的茶毛蟲幼蟲,縱切后取半頭經(jīng)液氮速凍、研磨,使用血液/細胞/組織基因組DNA提取試劑盒[天根生化科技(北京)有限公司]提取DNA,﹣20 ℃保存,用于茶毛蟲線粒體基因序列擴增;取另半頭茶毛蟲經(jīng)速凍研磨后采用TRIzol法提取總RNA,使用RQ1 RNase-Free DNase試劑盒(Promega)去除RNA樣品中的基因組DNA,再使用PrimeScriptTM 1st Strand cDNA Synthesis Kit試劑盒(TAKARA)反轉(zhuǎn)錄合成cDNA,﹣20 ℃保存,用于EpBYV的RdRp序列擴增。

    1.2.2 PCR擴增及測序

    茶毛蟲線粒體基因COI、ND5,以及EpBYV病毒中RdRp序列的PCR擴增引物如表2所示。使用Taq HS DNA聚合酶(TaKaRa)進行PCR擴增,總反應(yīng)體積為25 μL。擴增條件為95 ℃預(yù)變性3 min;95 ℃變性30 s,按表2中的退火溫度退火30 s,72 ℃延伸1 min,34

    個循環(huán);72 ℃終延伸10 min。擴增完成后取2 μL PCR產(chǎn)物進行1%的瓊脂糖凝膠電泳,檢測成功后將PCR產(chǎn)物送至北京六合華大基因有限責(zé)任公司進行雙向測序。測序得到的序列使用DNAMAN軟件拼接,序列翻譯后在NCBI數(shù)據(jù)庫中進行BLAST比對,驗證序列的準確性。

    1.2.3 數(shù)據(jù)整理與分析

    利用DnaSP 6.12.03軟件[27-28]計算茶毛蟲不同地理種群的單倍型多樣性(Haplotype diversity,h)、核苷酸多樣性(Nucleotide diversity,π)等遺傳多樣性指標。基因流(Nm)依據(jù)公式Nm=(1-FST)/4FST計算[29]。在Arlequin 3.5.2.2軟件[30]中使用Tajima’s D[31]和Fu’s Fs[32]統(tǒng)計進行中性檢驗,并進行種群的分子遺傳變異方差分析(Analysis of molecular variance,AMOVA),計算種群對間的遺傳分化指數(shù)FST和誤配等。基于Kimura 2-Parameter模型使用MEGA 7.0.26軟件[33]計算種群間的遺傳距離。采用Network 10.2軟件[34]繪制單倍型之間的網(wǎng)絡(luò)進化圖(Median-joining network),單倍型的最大似然法系統(tǒng)發(fā)育樹使用PhyML3.0(www.atgc-montpellier.fr/phyml)構(gòu)建[35]。

    2 結(jié)果與分析

    2.1 序列變異和遺傳多樣性

    本研究共采集到15個地理種群共148頭茶毛蟲幼蟲,均擴增獲得644 bp的COI序列和715 bp的ND5序列。將同一樣本的兩個序列片段分別合并,獲得148條長度為1 359 bp的序列,其中保守位點1 327個,變異位點32個,未發(fā)現(xiàn)堿基缺失或插入現(xiàn)象。在32個變異位點中自裔位點5個,簡約信息位點27個,均為兩堿基變異位點,共含有22個單倍型(表3),將M1中的COI和ND5序列上傳GenBank,獲得登錄號分別為PQ097955和PQ096005。遺傳多樣性分析顯示(表3),總?cè)后w的核苷酸多樣性指數(shù)為0.003 17,單倍型多樣性指數(shù)為0.880 68。15個地理種群中,永川、武隆、六盤水、無錫、恩施、杭州、桂林、清遠等地的8個種群不具有多態(tài)性位點,各項指數(shù)均為0,表明其種群內(nèi)部序列保守,遺傳多樣性極低;其余7個種群的核苷酸多樣性指數(shù)為0.000 34~0.001 96,單倍型多樣性指數(shù)為0.200 00~0.755 56,各種群間遺傳多樣性差異較為明顯,其中信陽和安康兩個種群的遺傳多樣性較高。

    EpBYV的RdRp序列擴增結(jié)果顯示,除重慶城口10個樣本外,其他138個樣本均成功獲得488 bp的RdRp序列,擴增成功率為93.24%。獲得的138條RdRp序列中均無堿基缺失或插入現(xiàn)象,其中保守位點有402個,變異位點86個,約占全長的17.62%,無堿基缺失或插入現(xiàn)象。86個變異位點中有自裔位點49個,簡約信息位點37個,其中兩堿基變異位點82個,三堿基變異位點4個,形成單倍型47個(表3),將R1序列上傳GenBank獲得登錄號為PQ096006。遺傳多樣性參數(shù)顯示(表3),總?cè)后w的核苷酸多樣性指數(shù)(0.017 95)和單倍型多樣性指數(shù)(0.935 26)均大于前文中茶毛蟲遺傳多樣性的相應(yīng)指標。相對于茶毛蟲線粒體序列無突變位點的8個種群,EpBYV的RdRp序列僅在桂林和恩施種群中無突變位點,表明EpBYV種群內(nèi)部的遺傳多樣性較高。13個EpBYV有突變位點的種群中,核苷酸多樣性指數(shù)為0.000 41~0.018 67,單倍型多樣性指數(shù)為0.200 00~1.000 00,其中杭州和永川兩個種群的單倍型多樣性指數(shù)均為1.000 00,遺傳多樣性最高,對應(yīng)的每個樣本單獨成為1個單倍型。

    2.2 單倍型及進化關(guān)系

    在茶毛蟲線粒體基因序列的22個單倍型中,種群間共享的單倍型有3個,占比較少(13.64%),表明茶毛蟲線粒體基因序列在種群間的遺傳分化程度較高。共享單倍型中M7的出現(xiàn)頻率最高(44次),占所有個體的29.73%,其次是M1(11次)和M11(9次)。從種群來看,永川、武隆、六盤水、無錫、恩施、杭州、桂林、清遠等地的8個種群內(nèi)部均無突變位點,各自僅有1種單倍型,其中六盤水、恩施、桂林和清遠4個地區(qū)的種群共享單倍型M7,此外信陽種群有5個單倍型,安康種群4個,遺傳分化程度相對較大。單倍型的進化網(wǎng)絡(luò)圖顯示,大部分來自相同種群的單倍型聚在一起(開州、城口、安康、信陽、安慶)。進化網(wǎng)絡(luò)圖在結(jié)構(gòu)上呈現(xiàn)明顯的兩個分枝(圖1A),上部一枝有9個單倍型,包含永川、開州、城口、安康4個種群的38個樣本;下部一枝有13個單倍型,由其他11個種群110個樣本組成。安康在傳統(tǒng)茶區(qū)劃分中被分在江北茶區(qū),但在地理位置上與上部一枝上的開州、城口共同位于秦巴山區(qū),而武隆雖隸屬于重慶市,但在地理位置上與恩施、銅仁共同位于武陵山區(qū),表明茶毛蟲線粒體基因序列的變異與種群的地理位置有一定的關(guān)系。采用最大似然法構(gòu)建22個單倍型的系統(tǒng)發(fā)育樹,顯示其分枝結(jié)構(gòu)與網(wǎng)絡(luò)進化圖相吻合,進化網(wǎng)絡(luò)圖中的上下兩枝在進化樹均形成單系群,并互為姊妹群(圖1B)。

    不同種群EpBYV的RdRp序列的47個單倍型中有7個共享單倍型,占比較?。?4.89%)。共享單倍型中R13出現(xiàn)的頻率最高(25次),占所有個體的18.12%,其次是R30(15次)。種群方面,桂林種群僅有1個單倍型,其他種群的單倍型數(shù)在2~10個,其中杭州和永川種群均是1個樣本為1個單倍型,與之對應(yīng)的是兩個種群的茶毛蟲線粒體基因序列在種群內(nèi)部無突變位點,由此可見EpBYV的RdRp序列比茶毛蟲線粒體序列的突變速率更快,種群內(nèi)遺傳分化更大。RdRp序列的單倍型數(shù)量較大,相對于茶毛蟲線粒體序列單倍型的進化網(wǎng)絡(luò)圖,其進化網(wǎng)絡(luò)圖結(jié)構(gòu)更為復(fù)雜,難以通過地理種群明確劃分枝系。大部分來自相同種群的單倍型較好地聚在一起,但永川、開州、安慶、六盤水種群的單倍型在圖中分布距離較遠(圖1C),表明這4個種群內(nèi)部出現(xiàn)較大的遺傳分化。此外,地理種群間的進化關(guān)系也更為復(fù)雜,在圖1A中,茶毛蟲的武隆種群與重慶的其他3個種群的進化距離較遠,在圖1C中,EpBYV的永川、開州種群的部分單倍型與武隆種群距離較遠,但另一些單倍型與武隆種群聚在一起(R4、R5、R9、R10、R12、R15)。47個單倍型的最大似然法系統(tǒng)發(fā)育樹與網(wǎng)絡(luò)進化圖結(jié)構(gòu)相似,在進化圖上聚在一起的單倍型也以較高的支持率在進化樹上聚為一枝(圖1D)。

    2.3 種群遺傳結(jié)構(gòu)分析

    多個水平的AMOVA分析均顯示茶毛蟲存在顯著的遺傳分化(表4),兩級水平分析顯示茶毛蟲具有顯著的遺傳結(jié)構(gòu)(ΦST=0.871 18,P<0.001),其中87.12%的變異來自于種群之間,僅12.88%的變異來自于種群內(nèi)部。105個種群對間的遺傳分化指數(shù)FST值有99個達到顯著水平,數(shù)值為0.290~1.000(表5),支持茶毛蟲多數(shù)地理種群間存在顯著的遺傳分化,與之對應(yīng)的是多數(shù)種群對間基因流Nm水平較低(表5)。三級水平AMOVA分析將地理種群按傳統(tǒng)茶區(qū)(表1)分成4個組時,茶毛蟲不具有顯著的遺傳結(jié)構(gòu)(ΦCT=0.024 09,P>0.05),遺傳變異主要來自于組內(nèi)種群間,即各地理種群之間的差異;但根據(jù)進化網(wǎng)絡(luò)結(jié)構(gòu)圖中上下兩枝分為2個組時,茶毛蟲存在顯著的遺傳結(jié)構(gòu)(ΦCT=0.421 23,P<0.01),茶毛蟲在兩個組之間顯著分化,組間變異占42.12%。對EpBYV二級水平AMOVA分析顯示其遺傳結(jié)構(gòu)顯著(ΦST=0.621 32,P<0.001),遺傳變異大部分來自于種群之間(62.13%),種群內(nèi)部的遺傳變異較?。?7.87%),但仍大于茶毛蟲各種群內(nèi)部的遺傳變異占比(12.88%)。與茶毛蟲類似,EpBYV種群對間的遺傳分化指數(shù)FST值大多數(shù)達到了顯著水平(表6),Nm值大多數(shù)小于1(表6),支持EpBYV的地理種群間具有顯著的遺傳分化,基因交流較少。三級水平AMOVA分析顯示,EpBYV無顯著的遺傳結(jié)構(gòu),兩種分組下ΦCT的P值均大于0.05,組間的遺傳變異占比均為負值。種群間遺傳距離與地理距離的相關(guān)性分析結(jié)果顯示,茶毛蟲(R2=0.010,P=0.305)和EpBYV(R2=0.012,P=0.296)的遺傳距離與地理距離的相關(guān)性不顯著(圖2),且多數(shù)EpBYV種群間的遺傳距離較茶毛蟲種群間更大。進一步計算發(fā)現(xiàn),EpBYV種群間的遺傳距離與茶毛蟲對應(yīng)種群

    間遺傳距離的相關(guān)性不顯著(R2=0.006,P=0.461),表明EpBYV進化方向不受茶毛蟲的影響。

    2.4 種群歷史動態(tài)分析

    采用中性檢驗和誤配分析進行茶毛蟲和EpBYV歷史動態(tài)分析。茶毛蟲7個有多態(tài)性位點的種群中,城口和寧德種群的Tajima’s D值顯著小于0,可能是種群擴張引起種群偏離了遺傳平衡狀態(tài),但Fu’s Fs值未達到顯著水平;其他種群的中性檢驗值均不顯著,說明群體處于較穩(wěn)定狀態(tài);總?cè)后w的中性檢驗雖均為負值,但未達到顯著水平(表3)。茶毛蟲總?cè)后w的誤配分布表現(xiàn)為單峰(圖3A),且誤配分布統(tǒng)計參數(shù)粗糙度指數(shù)(Raggedness index,rg)為0.036 97,期望與觀測錯配分布之間的平方和(SSD)為0.013 21,兩個值均較小,且未達到顯著水平,推測我國茶毛蟲群體可能經(jīng)歷了較弱的突然增長過程。EpBYV地理種群的中性檢驗顯示,六盤水種群的Tajima’s D值和Fu’s Fs值均顯著小于0,武隆、安慶種群的Tajima’s D值和杭州種群的Fu’s Fs值均顯著小于0,推測這4個種群可能經(jīng)歷了不同程度擴張。EpBYV總?cè)后w的Tajima’s D值和Fu’s Fs值均顯著小于0,其中Fu’s Fs值達到了極顯著水平,且誤配分布圖表現(xiàn)為明顯的單峰(圖3B),rg值為0.020 78(P=0.01),SSD值為0.011 51(P=0.07),表明EpBYV可能經(jīng)歷了突然增長過程。

    3 討論

    本研究采集并提取了茶毛蟲的DNA和RNA,以此為模板分別擴增茶毛蟲的線粒體COI、ND5和EpBYV的RdRp編碼序列,其中茶毛蟲線粒體兩個基因的擴增及比對成功率均為100%,表明所采集樣本均為茶毛蟲;而EpBYV的RdRp序列擴增成功率為93.24%,在城口種群10個樣本中未能成功擴增,原因可能是茶毛蟲城口種群未感染EpBYV或體內(nèi)病毒量極低。因此,148個茶毛蟲樣本的EpBYV感染率大于或等于93.24%,表明該病毒廣泛分布于茶毛蟲的田間種群中。

    基于線粒體COI和ND5合并序列分析茶毛蟲的種群遺傳結(jié)構(gòu),二級水平的AMOVA分析結(jié)果和遺傳分化指數(shù)ΦST值均支持茶毛蟲的遺

    傳變異主要來自于地理種群之間,種群之間遺傳分化極顯著。三級水平AMOVA分析顯示,我國傳統(tǒng)茶區(qū)劃分不能指示茶毛蟲地理種群間的遺傳分化程度,參考進化網(wǎng)絡(luò)結(jié)構(gòu)圖能較為正確地指示種群間的分化程度。將該分組置于我國地勢中可以發(fā)現(xiàn),圖1A中上部一枝的4個地理種群(永川、開州、城口、安康)位于我國地勢的第二階梯,而下部一枝的地理種群均位于第三階梯及第二、三階梯的分界上,推測茶毛蟲在兩個區(qū)域的顯著分化可能是兩個階梯海拔和氣候差異引起的。突變、自然選擇、基因流和隨機遺傳漂變(Random genetic drift)被公認為是物種進化的四大動力[29]。其中,突變、自然選擇和隨機遺傳漂變加大種群間的遺傳分化,基因流則可以阻止種群間的分化,因此茶毛蟲種群間低水平的基因流與高水平的遺傳分化相吻合,由此推測茶毛蟲擴散能力不強。同樣不具備長距離傳播能力的茶小綠葉蟬、茶網(wǎng)蝽和茶棍薊馬均顯示出較高的遺傳分化[22-24]。地理種群間的顯著遺傳分化也可能是造成不同地理種群茶毛蟲對同株EpNPV敏感性不同的原因。基于RdRp編碼序列分析EpBYV的種群遺傳結(jié)構(gòu),ΦST值顯示EpBYV種群之間遺傳分化顯著,二級水平AMOVA分析支持遺傳變異主要來自于地理種群之間,但種群內(nèi)部的遺傳變異占比(37.87%)明顯大于茶毛蟲種群內(nèi)部的遺傳變異占比(12.88%),遺傳多樣性參數(shù)(h=0.935 26、π=0.017 95)也支持EpBYV的種群具有更大的遺傳變異(表3),表明EpBYV與茶毛蟲的進化差異明顯,這與茶毛蟲成蟲飛翔能力較弱,對病毒水平傳播不足相關(guān)。病毒毒株的遺傳分化可能影響其對宿主的毒力水平[8],EpBYV的顯著遺傳分化和快速突變有利于進化出對茶毛蟲產(chǎn)生致病力的毒株。

    本研究的種群歷史動態(tài)分析采用中性檢驗和誤配分析相結(jié)合,當(dāng)中性檢驗Tajima’s D和Fu’s Fs值均為顯著負值時,說明所研究的種群可能經(jīng)歷了突然擴張事件[32],誤配分布表現(xiàn)為單峰表示種群在近期出現(xiàn)過群體擴張事件或持續(xù)增長模式[36-37]。茶毛蟲地理種群中,城口和寧德種群的Tajima’s D值顯著小于0,推測近期經(jīng)歷了群體擴張事件。茶毛蟲總?cè)后w的中性檢驗均為不顯著的負值,誤配分布表現(xiàn)為單峰,推測茶毛蟲總?cè)后w可能經(jīng)歷了較弱的群體擴張。EpBYV總?cè)后w的Tajima’s D值和Fu’s Fs值均顯著小于0,其中Fu’s Fs值達到了極顯著水平,誤配分布圖為明顯的單峰,表明EpBYV群體可能經(jīng)歷了突然增長過程,近似星狀分布的單倍型進化網(wǎng)絡(luò)圖也佐證EpBYV在近期經(jīng)歷了種群有效數(shù)量突然增長的過程。此外,RdRp序列具有較高水平的單倍型多樣性(0.935 26)和相對較低的核苷酸多樣性(0.017 95),也支持EpBYV種群經(jīng)歷了種群擴張事件[38]。EpBYV群體擴張事件主要出現(xiàn)在六盤水、武隆、安慶、杭州4個地理種群中。

    綜上所述,采自我國四大茶區(qū)15個地區(qū)的148個樣本均確認為茶毛蟲,且其群體相對穩(wěn)定,近期可能僅經(jīng)歷了較弱的群體擴張事件,僅城口和寧德種群經(jīng)歷了較明顯的群體擴張,表明茶毛蟲在部分茶區(qū)的適宜環(huán)境條件下也具有潛在暴發(fā)的風(fēng)險;在茶毛蟲總?cè)后w相對穩(wěn)定的背景下,EpBYV種群經(jīng)歷了明顯的種群擴張事件,說明EpBYV即使在經(jīng)歷種群萎縮后,仍可以迅速由少量個體建立龐大的種群。同時,EpBYV的種群擴張將增加病毒在復(fù)制時的突變機率,這有利于進化出對茶毛蟲產(chǎn)生致病力的毒株。EpBYV在茶毛蟲群體中的感染率和種群擴張能力均較高,用于茶毛蟲生物防治具有較好的潛力。

    參考文獻

    [1] 張漢鵠, 譚濟才. 中國茶樹害蟲及其無公害治理[M]. 合肥: 安徽科學(xué)技術(shù)出版社, 2004.

    Zhang H H, Tan J C. Chinese tea pests and their pollution-free control [M]. Hefei: Anhui Science and Technology Press, 2004.

    [2] 戈峰, 王永模, 劉向輝. 茶毛蟲性引誘劑誘殺效果分析[J]. 昆蟲知識, 2003, 40(3): 237-239.

    Ge F, Wang Y M, Liu X H. The effectiveness of synthetic sex pheromone by mass trapping as a control method for the tea tussock moth, Euproctis pseudoconspersa [J]. Entomological Knowledge, 2003, 40(3): 237-239.

    [3] 韓方岸, 胡云, 李世榮, 等. 茶毛蟲生態(tài)特性及人群皮炎暴發(fā)的觀察[J]. 中國寄生蟲病防治雜志, 2005, 18(5): 378-381.

    Han F A, Hu Y, Li S R, et al. Observation on the ecology of tea Euproctis pseudoconspersa strand and the outbreak of larvae dermatitis in crowd [J]. Chinese Journal of Parasitic Disease Control, 2005, 18(5): 378-381.

    [4] Dong W W, Dong S Y, Jiang G F, et al. Characterization of the complete mitochondrial genome of tea tussock moth, Euproctis pseudoconspersa (Lepidoptera: Lymantriidae) and its phylogenetic implications [J]. Gene, 2016, 577(1): 37-46.

    [5] 彭萍, 王曉慶, 李品武. 茶樹病蟲害測報與防治技術(shù)[M]. 北京: 中國農(nóng)業(yè)出版社, 2013.

    Peng P, Wang X Q, Li P W. Forecast and control techniques of tea diseases and pests [M]. Beijing: China Agriculture Press, 2013.

    [6] Hazarika L K, Bhuyan M, Hazarika B N. Insect pests of tea and their management [J]. Annual Review of Entomology, 2009, 54: 267-284.

    [7] Ye G Y, Xiao Q, Chen M, et al. Tea: biological control of insect and mite pests in China [J]. Biological Control, 2014, 68: 73-91.

    [8] 付建玉, 席羽, 唐美君, 等. 茶毛蟲核型多角體病毒不同分離株的毒力水平與遺傳結(jié)構(gòu)關(guān)系分析[J]. 茶葉科學(xué), 2011, 31(4): 289-294.

    Fu J Y, Xi Y, Tang M J, et al. Study on the relationship between virulence and genetic structure of four wild isolates of Euproctis pseudoconspersa nuclear polyhedrosis virus [J]. Journal of Tea Science, 2011, 31(4): 289-294.

    [9] 厲曉臘, 金軼偉, 柴一秋, 等. 茶毛蟲核型多角體病毒對茶毛蟲的致病性研究[J]. 茶葉科學(xué), 2006, 26(4): 265-269.

    Li X L, Jin Y W, Chai Y Q, et al. The pathogenicity of Euproctis pseudoconspersa nuclear polyhedrosis virus on the larvae of tea tussork (Euproctis pseudoconspersa Strand) [J]. Journal of Tea Science, 2006, 26(4): 265-269.

    [10] 王禮中, 肖強, 張傳溪. 茶毛蟲核型多角體病毒ph基因抗體的制備與利用[J]. 茶葉科學(xué), 2008, 28(4): 260-266.

    Wang L Z, Xiao Q, Zhang C X. Preparation of the antibody against EupsNPV polyhedrin and its utilization in detection of the viral pesticide [J]. Journal of Tea Science, 2008, 28(4): 260-266.

    [11] 李兆群. 我國茶園有害生物生物防治技術(shù)研究及應(yīng)用[J]. 中國茶葉, 2022, 44(5): 8-12.

    Li Z Q. Research and application of tea pest control technology in China [J]. China Tea, 2022, 44(5): 8-12.

    [12] 陳宗懋, 蔡曉明, 周利, 等. 中國茶園有害生物防控40年[J]. 中國茶葉, 2020, 42(1): 1-8.

    Chen Z M, Cai X M, Zhou L, et al. Developments on tea plant pest control in past 40 years in China [J]. China Tea, 2020, 42(1): 1-8.

    [13] 冷楊, 肖強, 殷坤山. 茶毛蟲核型多角體病毒Bt混劑的作用特性[J]. 植物保護學(xué)報, 2007, 34(2): 177-181.

    Leng Y, Xiao Q, Yin K S. The effect specialty of EpNPV-Bt preparation [J]. Acta Phytophylacica Sinica, 2007, 34(2): 177-181.

    [14] 唐美君, 殷坤山, 郭華偉, 等. 茶毛蟲核型多角體病毒和Bt混劑的配比篩選及藥效評價[J]. 植物保護, 2010, 36(5): 165-167.

    Tang M J, Yin K S, Guo H W, et al. Screening for ratio of EpNPV mixed with Bt and evaluation of its control effect [J]. Plant Protection, 2010, 36(5): 165-167.

    [15] 王曉慶. 茶毛蟲的發(fā)生規(guī)律及其感染的病毒分析[D]. 重慶: 西南大學(xué), 2021.

    Wang X Q. Study on the occurrence and viral infection of Euproctis pseudoconspersa [D]. Chongqing: Southwest University, 2021.

    [16] Wang X Q, Gu Q Y, Zhang W, et al. Prevalence of a novel bunyavirus in tea tussock moth Euproctis pseudoconspersa (Lepidoptera: Lymantriidae) [J]. Journal of Insect Science, 2021, 21(4): 5. doi: 10.1093/jisesa/ieab045.

    [17] Marklewitz M, Zirkel F, Kurth A, et al. Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family [J]. PNAS, 2015, 112(24): 7536-7541.

    [18] Bolling B G, Olea-Popelka F J, Eisen L, et al. Transmission dynamics of an insect-specific flavivirus in a naturally infected Culex pipiens laboratory colony and effects of co-infection on vector competence for West Nile virus [J]. Virology, 2012, 427(2): 90-97.

    [19] Newman C M, Cerutti F, Anderson T K, et al. Culex flavivirus and West Nile virus mosquito coinfection and positive ecological association in Chicago, United States [J]. Vector-Borne and Zoonotic Diseases, 2011, 11(8): 1099-1105.

    [20] 王星馳, 晏嫦妤. 茶毛蟲發(fā)生規(guī)律與生物防治研究進展[J]. 廣東茶業(yè), 2023(5): 16-22.

    Wang X C, Yan C Y. Research progress on occurrence regularity and biological control of Euproctis pseudoconspersa [J]. Guangdong Tea Industry, 2023(5): 16-22.

    [21] 馬蕊, 白阿先, 盧木樸. 4種生物農(nóng)藥防治云南綠春縣茶毛蟲的效果研究[J]. 昆明學(xué)院學(xué)報, 2022, 44(3): 53-56.

    Ma R, Bai A X, Lu M P. Effect of four kinds of biological pesticides on controlling Euproctis pseudoconspersa in Lüchun county of Yunnan province [J]. Journal of Kunming University, 2022, 44(3): 53-56.

    [22] 李金玉, 牛東升, 陳杰, 等. 茶園周邊景觀格局對茶小綠葉蟬種群遺傳結(jié)構(gòu)的影響[J]. 昆蟲學(xué)報, 2020, 63(10): 1242-1259.

    Li J Y, Niu D S, Chen J, et al. Effects of landscape pattern around tea plantation on the population genetic structure of the tea green leafhopper, Empoasca onukii (Hemiptera: Cicadellidae) [J]. Acta Entomologica Sinica, 2020, 63(10): 1242-1259.

    [23] 陳世春, 江宏燕, 廖姝然, 等. 基于COI基因解析我國茶網(wǎng)蝽種群遺傳多樣性和遺傳結(jié)構(gòu)[J]. 茶葉科學(xué), 2023, 43(6): 795-805.

    Chen S C, Jiang H Y, Liao S R, et al. Analysis of genetic diversity and genetic structure in geographic populations of Stephanitis chinensis from China based on mitochondrial DNA COI sequence [J]. Journal of Tea Science, 2023, 43(6): 795-805.

    [24] 羅林麗, 孟澤洪, 李帥, 等. 中國南方茶棍薊馬地理種群遺傳分化分析[J]. 昆蟲學(xué)報, 2022, 65(4): 500-511.

    Luo L L, Meng Z H, Li S, et al. Genetic differentiation analysis of geographical populations of Dendrothrips minowai (Thysanoptera: Thripidae) in South China [J]. Acta Entomologica Sinica, 2022, 65(4): 500-511.

    [25] Zhang D X, Hewitt G M. Assessment of the universality and utility of a set of conserved mitochondrial COI primers in insects [J]. Insect Molecular Biology, 1997, 6(2): 143-150.

    [26] Simon C, Buckley T R, Frati F, et al. Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA [J]. Annual Review of Ecology Evolution and Systematics, 2006, 37: 545-579.

    [27] Rozas J. DNA sequence polymorphism analysis using DnaSP [J]. Methods in Molecular Biology, 2009, 537: 337-350.

    [28] Rozas J, Ferrer-Mata A, Sanchez-DelBarrio J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets [J]. Molecular Biology and Evolution, 2017, 34(12): 3299-3302.

    [29] 魏丹丹. 書虱種群遺傳多樣性及線粒體基因組進化研究[D]. 重慶: 西南大學(xué), 2012.

    Wei D D. Population genetic diversity and mitochondrial genome analysis of Psocids (Psocoptera: Liposcelididae) [D]. Chongqing: Southwest University, 2012.

    [30] Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis [J]. Evolutionary Bioinformatics Online, 2005, 1: 47-50.

    [31] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J]. Genetics, 1989, 123(3): 585-595.

    [32] Fu Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection [J]. Genetics, 1997, 147(2): 915-925.

    [33] Chen J, Liu Q, Gao L W. Visual tea leaf disease recognition using a convolutional neural network model [J]. Symmetry, 2019, 11(3): 343. doi: 10.3390/sym11030343.

    [34] Bandelt H J, Forster P, R?hl A. Median-joining networks for inferring intraspecific phylogenies [J]. Molecular Biology and Evolution, 1999, 16(1): 37-48.

    [35] Guindon S, Dufayard J F, Lefort V, et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0 [J]. Systematic Biology, 2010, 59(3): 307-321.

    [36] Wei D D, Yuan M L, Wang B J, et al. Population genetics of two asexually and sexually reproducing psocids species inferred by the analysis of mitochondrial and nuclear DNA sequences [J]. Plos One, 2012, 7(3): e33883. doi: 10.1371/journal.pone.0033883.

    [37] Rogers A R, Harpending H. Population growth makes waves in the distribution of pairwise genetic differences [J]. Molecuar Biology and Evolution, 1992, 9(3): 552-569.

    [38] Slatkin M. Isolation by distance in equilibrium and non-equilibrium populations [J]. Evolution, 1993, 47(1): 264-279.

    最新美女视频免费是黄的| 久久人人精品亚洲av| 免费在线观看完整版高清| 亚洲一码二码三码区别大吗| 可以免费在线观看a视频的电影网站| 精品久久久久久久久久免费视频| 激情在线观看视频在线高清| 最好的美女福利视频网| 免费观看人在逋| 男人操女人黄网站| 窝窝影院91人妻| 动漫黄色视频在线观看| 欧美黑人欧美精品刺激| 久久午夜亚洲精品久久| 人妻丰满熟妇av一区二区三区| 亚洲自偷自拍图片 自拍| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| 欧美人与性动交α欧美精品济南到| 色尼玛亚洲综合影院| 精品福利观看| 搡老岳熟女国产| 99国产综合亚洲精品| 精品无人区乱码1区二区| 成人欧美大片| 亚洲成人久久性| 男女床上黄色一级片免费看| 午夜精品久久久久久毛片777| 久久影院123| av免费在线观看网站| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 免费不卡黄色视频| 欧美绝顶高潮抽搐喷水| 日韩 欧美 亚洲 中文字幕| 精品免费久久久久久久清纯| 免费久久久久久久精品成人欧美视频| 欧美激情久久久久久爽电影 | 麻豆av在线久日| 欧美另类亚洲清纯唯美| 国产精品免费视频内射| 一区二区日韩欧美中文字幕| 日韩av在线大香蕉| 色av中文字幕| 变态另类成人亚洲欧美熟女 | 亚洲 欧美一区二区三区| 亚洲熟女毛片儿| 在线永久观看黄色视频| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 精品国产美女av久久久久小说| 久久久久久久久久久久大奶| 亚洲国产中文字幕在线视频| 亚洲成人免费电影在线观看| 老汉色av国产亚洲站长工具| 日韩一卡2卡3卡4卡2021年| 真人做人爱边吃奶动态| 亚洲 国产 在线| 国产亚洲欧美在线一区二区| 国产区一区二久久| 日韩精品中文字幕看吧| 国产私拍福利视频在线观看| 在线观看免费日韩欧美大片| 日日摸夜夜添夜夜添小说| 久久久久久人人人人人| 99国产综合亚洲精品| 亚洲免费av在线视频| 欧美日本视频| av欧美777| 母亲3免费完整高清在线观看| 国产激情久久老熟女| 免费在线观看视频国产中文字幕亚洲| 在线观看免费视频网站a站| 亚洲黑人精品在线| 9色porny在线观看| 97碰自拍视频| а√天堂www在线а√下载| 亚洲第一欧美日韩一区二区三区| 侵犯人妻中文字幕一二三四区| 午夜福利一区二区在线看| 欧美成人一区二区免费高清观看 | 国产av一区二区精品久久| 亚洲五月天丁香| 黄色丝袜av网址大全| 女同久久另类99精品国产91| 国产黄a三级三级三级人| 操出白浆在线播放| 在线观看舔阴道视频| 色在线成人网| 亚洲国产欧美一区二区综合| 免费观看精品视频网站| 国产亚洲精品第一综合不卡| 男女下面插进去视频免费观看| 熟女少妇亚洲综合色aaa.| 手机成人av网站| 男女做爰动态图高潮gif福利片 | 一二三四在线观看免费中文在| 他把我摸到了高潮在线观看| 精品熟女少妇八av免费久了| 国产免费男女视频| 国产私拍福利视频在线观看| 美女 人体艺术 gogo| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 丝袜美腿诱惑在线| 久久久久精品国产欧美久久久| 亚洲熟女毛片儿| 国产欧美日韩精品亚洲av| 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 国内久久婷婷六月综合欲色啪| 亚洲人成网站在线播放欧美日韩| 欧美老熟妇乱子伦牲交| 国产色视频综合| 啦啦啦韩国在线观看视频| 伊人久久大香线蕉亚洲五| 俄罗斯特黄特色一大片| 精品日产1卡2卡| 精品少妇一区二区三区视频日本电影| 日韩欧美三级三区| 久久久国产成人免费| 在线国产一区二区在线| 成人三级黄色视频| 人成视频在线观看免费观看| 正在播放国产对白刺激| 亚洲中文字幕一区二区三区有码在线看 | 久久精品91无色码中文字幕| 久久热在线av| 午夜福利高清视频| 亚洲人成77777在线视频| 中出人妻视频一区二区| 精品福利观看| 人人妻人人爽人人添夜夜欢视频| 人人妻人人澡人人看| 国产精品亚洲美女久久久| avwww免费| 久久久久九九精品影院| 99精品在免费线老司机午夜| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 国产激情久久老熟女| 精品一区二区三区视频在线观看免费| 亚洲中文av在线| 最新在线观看一区二区三区| 日韩精品中文字幕看吧| 国产精品,欧美在线| 久久久国产成人免费| 午夜福利影视在线免费观看| 日韩精品免费视频一区二区三区| 久久性视频一级片| 麻豆一二三区av精品| 色在线成人网| 视频在线观看一区二区三区| 国产成人av激情在线播放| 神马国产精品三级电影在线观看 | 亚洲人成77777在线视频| 色av中文字幕| 激情视频va一区二区三区| 久久久久久国产a免费观看| 人人妻人人澡人人看| www日本在线高清视频| 欧美午夜高清在线| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 国产欧美日韩一区二区三| 久久久国产成人免费| 国产亚洲精品av在线| 婷婷六月久久综合丁香| 波多野结衣一区麻豆| 韩国av一区二区三区四区| 国产主播在线观看一区二区| 国产av又大| 久久精品亚洲精品国产色婷小说| 久久精品aⅴ一区二区三区四区| 后天国语完整版免费观看| 国产熟女xx| 两个人免费观看高清视频| 99精品久久久久人妻精品| 亚洲最大成人中文| 欧美老熟妇乱子伦牲交| 日本撒尿小便嘘嘘汇集6| 成年版毛片免费区| 变态另类丝袜制服| www国产在线视频色| 精品少妇一区二区三区视频日本电影| 久久久久久久午夜电影| 自拍欧美九色日韩亚洲蝌蚪91| 中文亚洲av片在线观看爽| 久久久久久大精品| 亚洲国产中文字幕在线视频| 纯流量卡能插随身wifi吗| 亚洲av第一区精品v没综合| 久久久久久久久免费视频了| av天堂在线播放| a级毛片在线看网站| 久久国产亚洲av麻豆专区| avwww免费| 操出白浆在线播放| 国语自产精品视频在线第100页| 99在线人妻在线中文字幕| 久久草成人影院| 国产99白浆流出| 精品国产乱码久久久久久男人| 国产精品永久免费网站| 亚洲专区中文字幕在线| 大型av网站在线播放| 啦啦啦观看免费观看视频高清 | 日韩欧美国产在线观看| 久久欧美精品欧美久久欧美| 18美女黄网站色大片免费观看| 丁香六月欧美| 国产高清激情床上av| 国产欧美日韩综合在线一区二区| 一边摸一边抽搐一进一出视频| 99国产精品免费福利视频| 中文字幕色久视频| 亚洲中文字幕一区二区三区有码在线看 | xxx96com| 不卡av一区二区三区| 午夜福利,免费看| 久久久久久亚洲精品国产蜜桃av| 午夜两性在线视频| 俄罗斯特黄特色一大片| 久久人人爽av亚洲精品天堂| 成人三级黄色视频| 99在线视频只有这里精品首页| 又黄又爽又免费观看的视频| 变态另类丝袜制服| www.精华液| 亚洲中文字幕一区二区三区有码在线看 | 欧美一级a爱片免费观看看 | 涩涩av久久男人的天堂| a在线观看视频网站| 搡老熟女国产l中国老女人| 国产一卡二卡三卡精品| 夜夜夜夜夜久久久久| 美女国产高潮福利片在线看| 国产精品日韩av在线免费观看 | 长腿黑丝高跟| 成人亚洲精品av一区二区| 91字幕亚洲| 色精品久久人妻99蜜桃| 91在线观看av| 欧美黑人欧美精品刺激| 成人亚洲精品av一区二区| 伊人久久大香线蕉亚洲五| 日韩大码丰满熟妇| 久久精品亚洲精品国产色婷小说| 十八禁网站免费在线| 91麻豆av在线| 亚洲 国产 在线| 亚洲人成电影观看| 中文字幕精品免费在线观看视频| 国产欧美日韩精品亚洲av| 亚洲国产中文字幕在线视频| 国产精品影院久久| 香蕉国产在线看| 激情在线观看视频在线高清| 免费高清在线观看日韩| 欧美大码av| 免费少妇av软件| 国产欧美日韩一区二区精品| 欧美日韩精品网址| 久久久精品欧美日韩精品| 精品一区二区三区av网在线观看| 久久狼人影院| 国产成人一区二区三区免费视频网站| 激情视频va一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 黄色毛片三级朝国网站| 国产91精品成人一区二区三区| 亚洲第一欧美日韩一区二区三区| 好男人电影高清在线观看| 91国产中文字幕| 亚洲一区中文字幕在线| 国产一区二区三区在线臀色熟女| 国产成人av激情在线播放| 精品人妻在线不人妻| 两人在一起打扑克的视频| 麻豆久久精品国产亚洲av| 1024视频免费在线观看| 在线十欧美十亚洲十日本专区| 90打野战视频偷拍视频| 夜夜躁狠狠躁天天躁| 在线观看一区二区三区| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 国产成人影院久久av| 欧美国产精品va在线观看不卡| 久久香蕉国产精品| 级片在线观看| 亚洲国产看品久久| 国产精品秋霞免费鲁丝片| 亚洲中文日韩欧美视频| 精品人妻在线不人妻| 操美女的视频在线观看| 欧美国产精品va在线观看不卡| 国产成年人精品一区二区| 国产精品永久免费网站| 精品国产一区二区三区四区第35| 日韩一卡2卡3卡4卡2021年| 少妇粗大呻吟视频| 一区二区三区精品91| 成年版毛片免费区| 亚洲av第一区精品v没综合| 亚洲午夜理论影院| 啦啦啦 在线观看视频| 两个人看的免费小视频| 久久亚洲精品不卡| 久久久久国产一级毛片高清牌| 波多野结衣高清无吗| 欧美色视频一区免费| 黄网站色视频无遮挡免费观看| 老鸭窝网址在线观看| 国产精品乱码一区二三区的特点 | 精品少妇一区二区三区视频日本电影| 国产高清有码在线观看视频 | 婷婷丁香在线五月| 男女午夜视频在线观看| 无人区码免费观看不卡| 国产成人欧美| 大香蕉久久成人网| 亚洲专区国产一区二区| 色播亚洲综合网| 亚洲最大成人中文| 制服丝袜大香蕉在线| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 国产精品一区二区免费欧美| 亚洲片人在线观看| 人妻久久中文字幕网| 精品久久久久久,| 欧美色欧美亚洲另类二区 | 亚洲avbb在线观看| 亚洲国产精品sss在线观看| 很黄的视频免费| 天天躁夜夜躁狠狠躁躁| 制服诱惑二区| 国产99白浆流出| 亚洲五月婷婷丁香| 精品不卡国产一区二区三区| 久久亚洲真实| 亚洲午夜精品一区,二区,三区| 免费观看人在逋| 午夜福利影视在线免费观看| 大陆偷拍与自拍| 亚洲av成人一区二区三| 亚洲av日韩精品久久久久久密| 日日夜夜操网爽| 国产精品一区二区三区四区久久 | 国产伦一二天堂av在线观看| 国产精品自产拍在线观看55亚洲| 精品久久久久久久人妻蜜臀av | 久久人妻av系列| 他把我摸到了高潮在线观看| 国产av一区在线观看免费| 久久婷婷人人爽人人干人人爱 | 啦啦啦 在线观看视频| 午夜免费鲁丝| 88av欧美| 精品一区二区三区av网在线观看| 人人妻人人爽人人添夜夜欢视频| xxx96com| 亚洲熟妇中文字幕五十中出| 黄片播放在线免费| 午夜福利,免费看| 国产精品一区二区免费欧美| 亚洲男人天堂网一区| 久久久久久久久免费视频了| 深夜精品福利| 淫妇啪啪啪对白视频| 99热只有精品国产| 久久国产精品人妻蜜桃| 欧美黑人欧美精品刺激| 人妻丰满熟妇av一区二区三区| 亚洲成av片中文字幕在线观看| 国产精品一区二区免费欧美| 99国产精品免费福利视频| 欧美一级a爱片免费观看看 | 国产又爽黄色视频| 99久久综合精品五月天人人| 国产真人三级小视频在线观看| 国产成人欧美在线观看| 久久香蕉精品热| 色精品久久人妻99蜜桃| 在线观看免费视频日本深夜| 丝袜人妻中文字幕| 动漫黄色视频在线观看| 好男人在线观看高清免费视频 | 国产亚洲欧美在线一区二区| 波多野结衣巨乳人妻| 久久久国产精品麻豆| 一区二区三区国产精品乱码| 亚洲av电影在线进入| 亚洲av成人一区二区三| 丝袜美腿诱惑在线| 在线十欧美十亚洲十日本专区| 69av精品久久久久久| 欧美日韩精品网址| www.www免费av| 精品国产超薄肉色丝袜足j| 91国产中文字幕| 老汉色∧v一级毛片| 男人操女人黄网站| 天天躁夜夜躁狠狠躁躁| 操美女的视频在线观看| 一级a爱片免费观看的视频| 一卡2卡三卡四卡精品乱码亚洲| 久久久久亚洲av毛片大全| 久久久久久亚洲精品国产蜜桃av| 巨乳人妻的诱惑在线观看| www.精华液| 1024香蕉在线观看| 怎么达到女性高潮| 美女高潮到喷水免费观看| 一级毛片高清免费大全| 成在线人永久免费视频| 欧美国产精品va在线观看不卡| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合一区二区三区| 国产精华一区二区三区| 少妇粗大呻吟视频| 性少妇av在线| 啦啦啦免费观看视频1| 久久伊人香网站| 国产精品免费一区二区三区在线| 怎么达到女性高潮| 亚洲狠狠婷婷综合久久图片| 超碰成人久久| 久久久久精品国产欧美久久久| 亚洲自拍偷在线| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 天堂影院成人在线观看| 天堂动漫精品| 久久 成人 亚洲| 黄色视频,在线免费观看| www日本在线高清视频| 久久天躁狠狠躁夜夜2o2o| 妹子高潮喷水视频| 久久这里只有精品19| 美女免费视频网站| 中文字幕最新亚洲高清| 女性生殖器流出的白浆| 欧美 亚洲 国产 日韩一| 亚洲国产精品成人综合色| 午夜两性在线视频| 黄色片一级片一级黄色片| 亚洲自拍偷在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲久久久国产精品| 波多野结衣巨乳人妻| 露出奶头的视频| www日本在线高清视频| 在线观看午夜福利视频| 欧美亚洲日本最大视频资源| 国内精品久久久久久久电影| 男人舔女人下体高潮全视频| xxx96com| 免费高清视频大片| 色综合婷婷激情| 欧美成人性av电影在线观看| 久久国产精品影院| 在线观看免费午夜福利视频| 波多野结衣av一区二区av| 日韩一卡2卡3卡4卡2021年| 女人高潮潮喷娇喘18禁视频| 欧美日韩福利视频一区二区| 激情视频va一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看66精品国产| 91av网站免费观看| 国产精品国产高清国产av| 男人的好看免费观看在线视频 | 麻豆成人av在线观看| 一级作爱视频免费观看| 免费少妇av软件| 九色国产91popny在线| 久久久国产精品麻豆| 岛国在线观看网站| 午夜精品在线福利| 午夜福利18| 久久久久国产一级毛片高清牌| 国产私拍福利视频在线观看| 亚洲av片天天在线观看| 亚洲黑人精品在线| 一进一出抽搐gif免费好疼| 女人精品久久久久毛片| 国产精品久久久久久精品电影 | 亚洲 国产 在线| 成人三级做爰电影| 久久精品影院6| 国产伦一二天堂av在线观看| www.www免费av| 97碰自拍视频| 国产精品 欧美亚洲| 在线天堂中文资源库| 99国产精品一区二区三区| 一区二区三区精品91| av中文乱码字幕在线| 亚洲国产中文字幕在线视频| 国产黄a三级三级三级人| 不卡av一区二区三区| 成人手机av| 国产高清视频在线播放一区| 国产区一区二久久| 亚洲第一青青草原| 欧美乱码精品一区二区三区| 久久精品成人免费网站| av天堂在线播放| 国产av一区在线观看免费| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧美激情在线| 欧美色欧美亚洲另类二区 | 一区二区三区国产精品乱码| 两个人视频免费观看高清| 精品少妇一区二区三区视频日本电影| 午夜精品在线福利| 97人妻精品一区二区三区麻豆 | 91九色精品人成在线观看| a级毛片在线看网站| 18禁观看日本| 日本三级黄在线观看| 亚洲av日韩精品久久久久久密| 老司机在亚洲福利影院| 欧美成狂野欧美在线观看| 精品人妻1区二区| 午夜久久久久精精品| 成在线人永久免费视频| 亚洲人成网站在线播放欧美日韩| av片东京热男人的天堂| 欧美午夜高清在线| 香蕉久久夜色| 亚洲五月婷婷丁香| 欧美av亚洲av综合av国产av| 精品一区二区三区av网在线观看| 亚洲av第一区精品v没综合| 亚洲精品中文字幕在线视频| 狂野欧美激情性xxxx| 中文字幕人妻丝袜一区二区| 少妇熟女aⅴ在线视频| tocl精华| 欧美在线黄色| 国产99久久九九免费精品| 久久精品国产亚洲av高清一级| 老司机靠b影院| 欧美成人免费av一区二区三区| 宅男免费午夜| 国产亚洲精品综合一区在线观看 | 黄频高清免费视频| 麻豆久久精品国产亚洲av| 午夜福利成人在线免费观看| 国产午夜精品久久久久久| 午夜激情av网站| 亚洲人成伊人成综合网2020| 亚洲,欧美精品.| 超碰成人久久| 精品熟女少妇八av免费久了| 两个人免费观看高清视频| 9191精品国产免费久久| 禁无遮挡网站| 国产xxxxx性猛交| 丁香六月欧美| 国产精品 国内视频| 婷婷精品国产亚洲av在线| 多毛熟女@视频| 免费高清在线观看日韩| 激情视频va一区二区三区| 99久久综合精品五月天人人| 女生性感内裤真人,穿戴方法视频| 国产精品98久久久久久宅男小说| 久久久精品国产亚洲av高清涩受| 18禁国产床啪视频网站| 精品久久久精品久久久| 好男人在线观看高清免费视频 | 欧美成人一区二区免费高清观看 | 久久影院123| 欧美激情极品国产一区二区三区| 亚洲专区国产一区二区| 性少妇av在线| 日本 av在线| 黄色a级毛片大全视频| 淫秽高清视频在线观看| 看片在线看免费视频| 国产99白浆流出| 亚洲激情在线av| 久久中文看片网| 欧美最黄视频在线播放免费| 嫩草影院精品99| 日韩精品中文字幕看吧| 日韩高清综合在线| 嫩草影视91久久| 午夜福利成人在线免费观看| 91成人精品电影| 久久久久国产一级毛片高清牌| 国产精品,欧美在线| 99精品在免费线老司机午夜| 18禁美女被吸乳视频| 在线观看www视频免费| 三级毛片av免费| 午夜福利一区二区在线看| 中文字幕高清在线视频| 人成视频在线观看免费观看| 久久久久久免费高清国产稀缺| 精品熟女少妇八av免费久了| 最好的美女福利视频网| 一卡2卡三卡四卡精品乱码亚洲| 国产xxxxx性猛交| 91大片在线观看| 啦啦啦 在线观看视频| 一边摸一边抽搐一进一出视频| 午夜免费激情av|