• <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.

    欧美日韩国产mv在线观看视频| 超碰97精品在线观看| 啦啦啦免费观看视频1| 十八禁网站免费在线| 满18在线观看网站| 女人精品久久久久毛片| 国产在视频线精品| 天天躁狠狠躁夜夜躁狠狠躁| 视频区欧美日本亚洲| 在线观看午夜福利视频| av天堂久久9| 久久精品亚洲精品国产色婷小说| 12—13女人毛片做爰片一| 国产一卡二卡三卡精品| 一边摸一边抽搐一进一小说 | 亚洲成人免费电影在线观看| 国产淫语在线视频| 日韩欧美免费精品| 亚洲成国产人片在线观看| 亚洲成人手机| 涩涩av久久男人的天堂| 成在线人永久免费视频| 99热网站在线观看| 99久久99久久久精品蜜桃| 午夜视频精品福利| 国产午夜精品久久久久久| 桃红色精品国产亚洲av| 精品第一国产精品| 国产精品免费视频内射| 在线免费观看的www视频| 老鸭窝网址在线观看| 亚洲一区高清亚洲精品| 国产亚洲精品一区二区www | 99久久综合精品五月天人人| 美国免费a级毛片| 亚洲 国产 在线| 一二三四社区在线视频社区8| 丝袜人妻中文字幕| 精品国产乱子伦一区二区三区| 中文亚洲av片在线观看爽 | 亚洲欧洲精品一区二区精品久久久| 又黄又爽又免费观看的视频| 51午夜福利影视在线观看| 婷婷精品国产亚洲av在线 | 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 国产一区二区三区在线臀色熟女 | 丝袜美足系列| 新久久久久国产一级毛片| 午夜福利在线免费观看网站| 欧美一级毛片孕妇| 久久国产精品影院| 亚洲,欧美精品.| 国产成人影院久久av| 中文字幕最新亚洲高清| 欧美日韩亚洲高清精品| 国产欧美日韩精品亚洲av| 欧美精品一区二区免费开放| 91精品国产国语对白视频| 久久久精品区二区三区| 色在线成人网| 在线观看免费视频网站a站| x7x7x7水蜜桃| 亚洲精品国产色婷婷电影| 看免费av毛片| 一级黄色大片毛片| tube8黄色片| 99精品欧美一区二区三区四区| 国产有黄有色有爽视频| 久久九九热精品免费| bbb黄色大片| 久久久精品国产亚洲av高清涩受| 国产高清videossex| 啦啦啦在线免费观看视频4| 搡老熟女国产l中国老女人| 中文字幕高清在线视频| www日本在线高清视频| 麻豆成人av在线观看| 亚洲欧美日韩高清在线视频| 女人爽到高潮嗷嗷叫在线视频| 超碰97精品在线观看| 国产野战对白在线观看| 午夜视频精品福利| 老汉色av国产亚洲站长工具| 51午夜福利影视在线观看| 男女免费视频国产| 男人操女人黄网站| 欧美国产精品一级二级三级| 一级毛片精品| 亚洲精品成人av观看孕妇| 搡老岳熟女国产| 十八禁网站免费在线| 露出奶头的视频| 亚洲九九香蕉| 交换朋友夫妻互换小说| 国产色视频综合| 久久九九热精品免费| 狂野欧美激情性xxxx| xxx96com| 在线观看舔阴道视频| 欧美不卡视频在线免费观看 | 国产精品久久视频播放| 午夜免费鲁丝| 久久精品国产a三级三级三级| 亚洲精品国产精品久久久不卡| 一区在线观看完整版| 男男h啪啪无遮挡| 亚洲少妇的诱惑av| 国产成人精品无人区| 色尼玛亚洲综合影院| 欧美成人午夜精品| 成人免费观看视频高清| 人妻久久中文字幕网| 国产精品一区二区精品视频观看| 亚洲熟妇熟女久久| 国产成人欧美在线观看 | 男人操女人黄网站| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看 | 欧美不卡视频在线免费观看 | 亚洲色图 男人天堂 中文字幕| av国产精品久久久久影院| 99久久人妻综合| 中文字幕另类日韩欧美亚洲嫩草| 搡老乐熟女国产| 身体一侧抽搐| 美女高潮到喷水免费观看| 亚洲一区二区三区欧美精品| 91精品国产国语对白视频| 亚洲美女黄片视频| 亚洲一区二区三区欧美精品| 午夜日韩欧美国产| 丁香欧美五月| 日韩大码丰满熟妇| 99久久99久久久精品蜜桃| 精品一区二区三卡| 欧美日韩视频精品一区| 又黄又粗又硬又大视频| 午夜福利,免费看| bbb黄色大片| 美国免费a级毛片| 亚洲国产看品久久| 成年人黄色毛片网站| 国产激情欧美一区二区| 亚洲美女黄片视频| 欧美乱妇无乱码| 真人做人爱边吃奶动态| avwww免费| 日韩人妻精品一区2区三区| av天堂久久9| 悠悠久久av| 成在线人永久免费视频| 欧美另类亚洲清纯唯美| 男人舔女人的私密视频| 一级片免费观看大全| 久久影院123| 欧美 日韩 精品 国产| 91成年电影在线观看| 欧美成人免费av一区二区三区 | 亚洲一码二码三码区别大吗| 日日摸夜夜添夜夜添小说| 人人妻,人人澡人人爽秒播| 亚洲七黄色美女视频| 伦理电影免费视频| 久久久国产欧美日韩av| 香蕉国产在线看| 亚洲欧美一区二区三区久久| 欧美激情 高清一区二区三区| 午夜精品久久久久久毛片777| 天天躁日日躁夜夜躁夜夜| 日韩欧美一区视频在线观看| 精品久久久久久电影网| 午夜两性在线视频| 18禁美女被吸乳视频| 日韩欧美在线二视频 | 十分钟在线观看高清视频www| 一二三四在线观看免费中文在| 亚洲熟妇熟女久久| 丰满人妻熟妇乱又伦精品不卡| 国产高清视频在线播放一区| 一级片免费观看大全| 日日夜夜操网爽| 啦啦啦免费观看视频1| 人成视频在线观看免费观看| 女人被狂操c到高潮| 亚洲国产欧美网| 久久久精品区二区三区| 国产又爽黄色视频| 19禁男女啪啪无遮挡网站| 亚洲黑人精品在线| 国产成人一区二区三区免费视频网站| 精品一区二区三区av网在线观看| 欧美日韩精品网址| 黑人巨大精品欧美一区二区蜜桃| 国产精品二区激情视频| 男女免费视频国产| 亚洲免费av在线视频| 亚洲中文日韩欧美视频| 一级a爱片免费观看的视频| a在线观看视频网站| 9191精品国产免费久久| 国产精品二区激情视频| 69精品国产乱码久久久| xxxhd国产人妻xxx| 又黄又爽又免费观看的视频| 新久久久久国产一级毛片| 午夜福利免费观看在线| 久9热在线精品视频| 久久午夜综合久久蜜桃| www日本在线高清视频| 丁香欧美五月| 亚洲av电影在线进入| 日本wwww免费看| 啦啦啦在线免费观看视频4| av免费在线观看网站| 他把我摸到了高潮在线观看| 一级,二级,三级黄色视频| 老司机福利观看| av不卡在线播放| 久久久国产欧美日韩av| 99国产综合亚洲精品| 国产成人av教育| 看片在线看免费视频| 老熟妇仑乱视频hdxx| 国产精品av久久久久免费| 欧美精品av麻豆av| 精品一区二区三卡| 99热国产这里只有精品6| 精品福利观看| 免费女性裸体啪啪无遮挡网站| xxx96com| 丰满的人妻完整版| e午夜精品久久久久久久| 日韩中文字幕欧美一区二区| 叶爱在线成人免费视频播放| 在线观看免费高清a一片| 国产高清激情床上av| 人人妻人人澡人人看| 少妇 在线观看| 成在线人永久免费视频| 久久国产亚洲av麻豆专区| 国产精品久久久人人做人人爽| 国产麻豆69| 一区二区三区激情视频| 老鸭窝网址在线观看| 99久久99久久久精品蜜桃| 99热只有精品国产| 欧美一级毛片孕妇| 国产精品98久久久久久宅男小说| 久久人妻av系列| 啦啦啦在线免费观看视频4| 亚洲一区中文字幕在线| 免费在线观看完整版高清| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品一区二区www | 中文字幕另类日韩欧美亚洲嫩草| √禁漫天堂资源中文www| a级毛片黄视频| a级片在线免费高清观看视频| 久久天躁狠狠躁夜夜2o2o| 多毛熟女@视频| 999久久久精品免费观看国产| 精品视频人人做人人爽| 国产黄色免费在线视频| 午夜成年电影在线免费观看| 成年人午夜在线观看视频| 无遮挡黄片免费观看| 一个人免费在线观看的高清视频| 欧美色视频一区免费| 精品国产乱子伦一区二区三区| 中文字幕制服av| 国产午夜精品久久久久久| 99精品久久久久人妻精品| 久久 成人 亚洲| 久久久久久久国产电影| 国产男女超爽视频在线观看| 国产成人av教育| 女人爽到高潮嗷嗷叫在线视频| 日韩 欧美 亚洲 中文字幕| 亚洲色图 男人天堂 中文字幕| 97人妻天天添夜夜摸| 国产精品av久久久久免费| 99精品在免费线老司机午夜| 国产欧美日韩一区二区三| 老司机福利观看| av中文乱码字幕在线| 午夜精品久久久久久毛片777| 精品少妇一区二区三区视频日本电影| 最近最新免费中文字幕在线| 一本大道久久a久久精品| 亚洲中文av在线| 久久午夜综合久久蜜桃| 91麻豆精品激情在线观看国产 | 色婷婷久久久亚洲欧美| 欧美日韩乱码在线| 美国免费a级毛片| 成人免费观看视频高清| 水蜜桃什么品种好| 久久久久久久精品吃奶| 亚洲一区中文字幕在线| 男女下面插进去视频免费观看| 天天躁日日躁夜夜躁夜夜| 国产一区二区激情短视频| xxx96com| 黄色女人牲交| 午夜成年电影在线免费观看| 精品人妻在线不人妻| 精品久久蜜臀av无| av中文乱码字幕在线| 操出白浆在线播放| 免费在线观看影片大全网站| 国产av又大| 在线观看免费视频日本深夜| 交换朋友夫妻互换小说| 岛国在线观看网站| 亚洲aⅴ乱码一区二区在线播放 | 两个人免费观看高清视频| 777米奇影视久久| 日本五十路高清| 免费在线观看黄色视频的| 久久久久久免费高清国产稀缺| a级毛片在线看网站| 老司机深夜福利视频在线观看| 成人免费观看视频高清| 在线国产一区二区在线| 亚洲熟妇熟女久久| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 国产人伦9x9x在线观看| 久久中文看片网| av一本久久久久| 成年动漫av网址| 欧美成人免费av一区二区三区 | 午夜老司机福利片| 亚洲av成人不卡在线观看播放网| 国产成人一区二区三区免费视频网站| 99久久人妻综合| 视频区图区小说| 午夜精品国产一区二区电影| 久久久精品国产亚洲av高清涩受| 十八禁高潮呻吟视频| 无人区码免费观看不卡| 亚洲av成人不卡在线观看播放网| 日本撒尿小便嘘嘘汇集6| 自线自在国产av| 国产极品粉嫩免费观看在线| 欧美+亚洲+日韩+国产| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产区一区二| 别揉我奶头~嗯~啊~动态视频| 国产一区二区三区视频了| 男女下面插进去视频免费观看| 一级毛片精品| 日韩欧美一区视频在线观看| 十八禁人妻一区二区| 亚洲成人手机| 午夜免费观看网址| 老司机午夜十八禁免费视频| 久久人人爽av亚洲精品天堂| 变态另类成人亚洲欧美熟女 | netflix在线观看网站| 日韩欧美一区视频在线观看| 国产单亲对白刺激| 一进一出抽搐gif免费好疼 | 欧美日韩一级在线毛片| 无限看片的www在线观看| 国产精品久久久人人做人人爽| 日韩一卡2卡3卡4卡2021年| 欧美日韩av久久| 亚洲熟妇中文字幕五十中出 | 精品久久蜜臀av无| 久久久精品区二区三区| 久久人妻福利社区极品人妻图片| 久久ye,这里只有精品| 大型黄色视频在线免费观看| 香蕉丝袜av| 亚洲情色 制服丝袜| 新久久久久国产一级毛片| 青草久久国产| 岛国在线观看网站| 岛国毛片在线播放| 亚洲avbb在线观看| 国产高清视频在线播放一区| 亚洲av第一区精品v没综合| 制服人妻中文乱码| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 看片在线看免费视频| 美女高潮到喷水免费观看| 99国产精品99久久久久| 成在线人永久免费视频| 高清欧美精品videossex| 性色av乱码一区二区三区2| 国产免费男女视频| 人成视频在线观看免费观看| 国产亚洲精品久久久久久毛片 | 人人妻,人人澡人人爽秒播| 亚洲成a人片在线一区二区| 久久国产乱子伦精品免费另类| 国产精品亚洲av一区麻豆| 免费少妇av软件| 人妻久久中文字幕网| 欧美黑人精品巨大| 黄片播放在线免费| 久久香蕉国产精品| 香蕉久久夜色| 久久久久久久精品吃奶| 国产精品亚洲一级av第二区| tocl精华| 国产精品综合久久久久久久免费 | 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 69av精品久久久久久| av免费在线观看网站| 亚洲人成电影观看| 午夜视频精品福利| 女人精品久久久久毛片| 亚洲国产看品久久| 最新美女视频免费是黄的| 亚洲精品美女久久久久99蜜臀| 一级毛片精品| 少妇裸体淫交视频免费看高清 | 中文字幕制服av| 亚洲人成伊人成综合网2020| 国产真人三级小视频在线观看| 午夜成年电影在线免费观看| 男女床上黄色一级片免费看| 亚洲色图综合在线观看| 亚洲性夜色夜夜综合| 久久精品亚洲熟妇少妇任你| 亚洲五月天丁香| 在线观看www视频免费| 69精品国产乱码久久久| videos熟女内射| av网站在线播放免费| 欧美午夜高清在线| 怎么达到女性高潮| 老司机影院毛片| 欧美国产精品va在线观看不卡| 交换朋友夫妻互换小说| 精品一区二区三区视频在线观看免费 | 欧美黄色片欧美黄色片| 99国产精品99久久久久| 在线看a的网站| 精品少妇一区二区三区视频日本电影| 亚洲一卡2卡3卡4卡5卡精品中文| 女人高潮潮喷娇喘18禁视频| 美女扒开内裤让男人捅视频| 免费在线观看完整版高清| 国产成人免费无遮挡视频| 久久国产精品大桥未久av| 亚洲国产精品合色在线| 丝袜人妻中文字幕| 大型黄色视频在线免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩一区二区三| 精品高清国产在线一区| 亚洲精品美女久久av网站| 一区二区三区国产精品乱码| 精品免费久久久久久久清纯 | 久久精品成人免费网站| 国产无遮挡羞羞视频在线观看| 热99re8久久精品国产| 久久精品国产a三级三级三级| 亚洲成a人片在线一区二区| 精品视频人人做人人爽| 黄色怎么调成土黄色| 亚洲性夜色夜夜综合| 亚洲午夜精品一区,二区,三区| 一区二区三区精品91| 久久人人97超碰香蕉20202| 老司机午夜福利在线观看视频| av视频免费观看在线观看| 国产成人av激情在线播放| 亚洲人成77777在线视频| 亚洲国产精品合色在线| 如日韩欧美国产精品一区二区三区| 自线自在国产av| 亚洲综合色网址| 人人妻人人爽人人添夜夜欢视频| 王馨瑶露胸无遮挡在线观看| 夜夜爽天天搞| 亚洲国产欧美日韩在线播放| 少妇的丰满在线观看| 少妇被粗大的猛进出69影院| 少妇猛男粗大的猛烈进出视频| 少妇裸体淫交视频免费看高清 | av中文乱码字幕在线| 午夜久久久在线观看| 一二三四社区在线视频社区8| 香蕉久久夜色| 国产成人影院久久av| 日韩制服丝袜自拍偷拍| 日韩欧美在线二视频 | 国产在视频线精品| a在线观看视频网站| 丝瓜视频免费看黄片| 婷婷丁香在线五月| 99精品欧美一区二区三区四区| 韩国av一区二区三区四区| 色精品久久人妻99蜜桃| 伊人久久大香线蕉亚洲五| 久久久国产一区二区| 人妻一区二区av| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品自拍成人| 乱人伦中国视频| 免费av中文字幕在线| 国产三级黄色录像| 9热在线视频观看99| 久9热在线精品视频| 精品久久久久久,| 丰满迷人的少妇在线观看| 亚洲性夜色夜夜综合| 亚洲欧美日韩另类电影网站| x7x7x7水蜜桃| 亚洲视频免费观看视频| 黑人操中国人逼视频| 亚洲五月天丁香| 自线自在国产av| √禁漫天堂资源中文www| 美女高潮到喷水免费观看| 美女国产高潮福利片在线看| 在线观看免费午夜福利视频| 国产精品香港三级国产av潘金莲| 欧美最黄视频在线播放免费 | 亚洲色图 男人天堂 中文字幕| 久久久久久久久免费视频了| 精品国产乱码久久久久久男人| 亚洲中文字幕日韩| 亚洲一区高清亚洲精品| 深夜精品福利| 超碰成人久久| 亚洲一卡2卡3卡4卡5卡精品中文| 精品视频人人做人人爽| 国产成人免费无遮挡视频| 免费人成视频x8x8入口观看| 91字幕亚洲| 免费久久久久久久精品成人欧美视频| 国产亚洲欧美在线一区二区| 老熟妇乱子伦视频在线观看| 国产成人av教育| 国产免费av片在线观看野外av| 精品卡一卡二卡四卡免费| 视频在线观看一区二区三区| 99久久人妻综合| 欧美精品av麻豆av| 香蕉丝袜av| 国产精品久久久久久精品古装| 久久精品国产99精品国产亚洲性色 | 亚洲色图综合在线观看| 亚洲黑人精品在线| 韩国精品一区二区三区| 亚洲欧洲精品一区二区精品久久久| 亚洲熟妇熟女久久| 国产乱人伦免费视频| 日韩有码中文字幕| 国产淫语在线视频| 免费一级毛片在线播放高清视频 | 少妇裸体淫交视频免费看高清 | 亚洲中文日韩欧美视频| 视频在线观看一区二区三区| 新久久久久国产一级毛片| xxx96com| 久久99一区二区三区| 纯流量卡能插随身wifi吗| 国产精品免费大片| 午夜福利在线免费观看网站| 国产亚洲一区二区精品| 日韩一卡2卡3卡4卡2021年| 在线免费观看的www视频| 欧美 亚洲 国产 日韩一| 久久香蕉激情| 国产免费av片在线观看野外av| 韩国av一区二区三区四区| 欧美成狂野欧美在线观看| 一级a爱片免费观看的视频| 日日爽夜夜爽网站| 69av精品久久久久久| 丝瓜视频免费看黄片| 丰满迷人的少妇在线观看| 日本a在线网址| 亚洲全国av大片| av一本久久久久| 欧美日韩乱码在线| 国产欧美日韩一区二区三区在线| 高清在线国产一区| 香蕉丝袜av| 老熟妇仑乱视频hdxx| 国产片内射在线| 黄频高清免费视频| 午夜视频精品福利| 伦理电影免费视频| 亚洲av第一区精品v没综合| 男人舔女人的私密视频| 久久影院123| 国产成人精品久久二区二区91| 91字幕亚洲| 亚洲情色 制服丝袜| 99国产精品一区二区蜜桃av | 欧美中文综合在线视频| 99精国产麻豆久久婷婷| 热re99久久精品国产66热6| 欧美另类亚洲清纯唯美| 黑丝袜美女国产一区| 丝袜在线中文字幕| 国内久久婷婷六月综合欲色啪| 国产有黄有色有爽视频| ponron亚洲| 成人影院久久| 亚洲七黄色美女视频| 嫩草影视91久久| 黄片播放在线免费| av超薄肉色丝袜交足视频| 精品久久久久久,|