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

    欧美亚洲日本最大视频资源| 久久影院123| 亚洲欧美日韩高清在线视频 | 久久狼人影院| 亚洲av美国av| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区精品| 欧美一级毛片孕妇| 99精品在免费线老司机午夜| 侵犯人妻中文字幕一二三四区| 新久久久久国产一级毛片| 日本精品一区二区三区蜜桃| 亚洲黑人精品在线| 日韩一卡2卡3卡4卡2021年| 国产精品 欧美亚洲| 在线观看舔阴道视频| 欧美成人午夜精品| 少妇精品久久久久久久| 日韩三级视频一区二区三区| 国产黄色免费在线视频| 色婷婷av一区二区三区视频| 一本—道久久a久久精品蜜桃钙片| 最新美女视频免费是黄的| 不卡av一区二区三区| 另类精品久久| 午夜福利欧美成人| 精品久久久久久久毛片微露脸| 91国产中文字幕| 欧美中文综合在线视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品中文字幕在线视频| 国产高清激情床上av| 精品国产乱码久久久久久男人| 69精品国产乱码久久久| 成人18禁在线播放| 国产精品熟女久久久久浪| 亚洲 国产 在线| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美一区二区综合| 在线av久久热| 啦啦啦视频在线资源免费观看| 欧美成人午夜精品| 久久精品国产亚洲av香蕉五月 | 桃红色精品国产亚洲av| 成人国产av品久久久| 欧美在线黄色| 精品久久久久久久毛片微露脸| 欧美人与性动交α欧美精品济南到| 久久久久精品国产欧美久久久| 免费日韩欧美在线观看| 国产高清视频在线播放一区| 欧美老熟妇乱子伦牲交| 精品福利永久在线观看| 亚洲免费av在线视频| 国产精品电影一区二区三区 | 老熟女久久久| 亚洲视频免费观看视频| 久久久久久久久久久久大奶| 国产成人欧美在线观看 | 黄色成人免费大全| 午夜福利视频在线观看免费| 国产亚洲欧美在线一区二区| 国产欧美日韩一区二区精品| 成人18禁高潮啪啪吃奶动态图| 国产精品国产高清国产av | 中文字幕人妻丝袜一区二区| 中国美女看黄片| 欧美日韩亚洲国产一区二区在线观看 | 嫩草影视91久久| 另类亚洲欧美激情| 亚洲av欧美aⅴ国产| 在线av久久热| 国产在视频线精品| 老熟女久久久| 露出奶头的视频| 精品熟女少妇八av免费久了| 精品少妇久久久久久888优播| 日韩一区二区三区影片| 亚洲精品久久成人aⅴ小说| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频| 精品卡一卡二卡四卡免费| 亚洲国产成人一精品久久久| 久久久久精品人妻al黑| 日本av免费视频播放| 中文字幕最新亚洲高清| 在线观看66精品国产| 久热这里只有精品99| 99国产精品一区二区蜜桃av | av网站在线播放免费| 久久精品国产综合久久久| 欧美日韩亚洲高清精品| 亚洲精品美女久久久久99蜜臀| 汤姆久久久久久久影院中文字幕| 高清在线国产一区| 久久国产精品影院| 国产精品国产av在线观看| 人成视频在线观看免费观看| av又黄又爽大尺度在线免费看| 欧美性长视频在线观看| 一个人免费看片子| 丝袜美腿诱惑在线| 亚洲国产av新网站| 91精品国产国语对白视频| 久久婷婷成人综合色麻豆| 91精品三级在线观看| 久久午夜亚洲精品久久| 人人妻人人添人人爽欧美一区卜| 蜜桃国产av成人99| 久久久久精品人妻al黑| 欧美人与性动交α欧美精品济南到| 50天的宝宝边吃奶边哭怎么回事| 亚洲成人免费电影在线观看| 久久毛片免费看一区二区三区| a级毛片在线看网站| √禁漫天堂资源中文www| 岛国在线观看网站| 精品一区二区三卡| 天天操日日干夜夜撸| 成人三级做爰电影| 欧美精品啪啪一区二区三区| 巨乳人妻的诱惑在线观看| 中文字幕人妻丝袜制服| 久久久久视频综合| 757午夜福利合集在线观看| 啦啦啦免费观看视频1| 男女下面插进去视频免费观看| 精品久久蜜臀av无| 国产精品熟女久久久久浪| 啦啦啦免费观看视频1| 久久久久久久精品吃奶| 美女扒开内裤让男人捅视频| 下体分泌物呈黄色| 国产高清激情床上av| 亚洲精品自拍成人| 丰满迷人的少妇在线观看| 日本黄色视频三级网站网址 | 狠狠婷婷综合久久久久久88av| 亚洲伊人久久精品综合| 久久久久久久久免费视频了| 国产在视频线精品| 亚洲美女黄片视频| 国产伦人伦偷精品视频| 大片电影免费在线观看免费| 曰老女人黄片| 午夜福利一区二区在线看| 免费看十八禁软件| 在线播放国产精品三级| 日韩中文字幕欧美一区二区| 高清视频免费观看一区二区| 午夜精品久久久久久毛片777| 日日夜夜操网爽| 久久人人97超碰香蕉20202| 成人亚洲精品一区在线观看| 夜夜夜夜夜久久久久| 99国产精品一区二区三区| 国产不卡一卡二| 最近最新免费中文字幕在线| 正在播放国产对白刺激| 欧美亚洲日本最大视频资源| 亚洲av日韩精品久久久久久密| 亚洲国产av新网站| 两性午夜刺激爽爽歪歪视频在线观看 | 中国美女看黄片| 久久人人爽av亚洲精品天堂| 欧美日韩亚洲高清精品| 黑丝袜美女国产一区| 成年女人毛片免费观看观看9 | 欧美国产精品va在线观看不卡| 精品国产一区二区三区四区第35| 天堂中文最新版在线下载| 熟女少妇亚洲综合色aaa.| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区三| 后天国语完整版免费观看| 少妇猛男粗大的猛烈进出视频| 免费看a级黄色片| 精品免费久久久久久久清纯 | 王馨瑶露胸无遮挡在线观看| 一个人免费看片子| 久久精品91无色码中文字幕| 黑人欧美特级aaaaaa片| 欧美 日韩 精品 国产| 精品久久久久久久毛片微露脸| 久久精品国产99精品国产亚洲性色 | 搡老熟女国产l中国老女人| 国产成人系列免费观看| 满18在线观看网站| 久久国产精品人妻蜜桃| 少妇粗大呻吟视频| 大陆偷拍与自拍| 国产亚洲欧美精品永久| 在线观看66精品国产| 国产在线视频一区二区| av片东京热男人的天堂| 老司机在亚洲福利影院| 免费在线观看完整版高清| 人人妻人人澡人人看| 99久久99久久久精品蜜桃| 欧美在线黄色| 大香蕉久久成人网| 国产成人精品久久二区二区免费| 国产成人系列免费观看| 十八禁人妻一区二区| 亚洲精品av麻豆狂野| 丁香六月天网| 精品福利永久在线观看| 人人妻人人添人人爽欧美一区卜| 99久久人妻综合| 最近最新中文字幕大全免费视频| 91大片在线观看| 自线自在国产av| 最近最新中文字幕大全免费视频| 久久热在线av| 国产欧美日韩一区二区三区在线| 国产av精品麻豆| 亚洲avbb在线观看| 国产av国产精品国产| 国产日韩欧美亚洲二区| 交换朋友夫妻互换小说| 亚洲情色 制服丝袜| 韩国精品一区二区三区| 国产成人精品无人区| 国产高清国产精品国产三级| 老熟妇仑乱视频hdxx| 亚洲av第一区精品v没综合| av天堂久久9| 亚洲专区中文字幕在线| 狠狠狠狠99中文字幕| 成人18禁高潮啪啪吃奶动态图| 欧美+亚洲+日韩+国产| 满18在线观看网站| 国产在线精品亚洲第一网站| 亚洲少妇的诱惑av| 一级,二级,三级黄色视频| 日日摸夜夜添夜夜添小说| 少妇精品久久久久久久| 最近最新免费中文字幕在线| 日日爽夜夜爽网站| 日韩视频在线欧美| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 久久狼人影院| 99香蕉大伊视频| 久久精品亚洲精品国产色婷小说| 国产不卡av网站在线观看| 精品国产乱码久久久久久小说| 欧美国产精品va在线观看不卡| 一本一本久久a久久精品综合妖精| 99久久精品国产亚洲精品| 亚洲专区国产一区二区| 亚洲中文av在线| 婷婷丁香在线五月| 国产无遮挡羞羞视频在线观看| 国产精品自产拍在线观看55亚洲 | 国产成人精品久久二区二区免费| 日韩精品免费视频一区二区三区| 91国产中文字幕| 丰满少妇做爰视频| 国产91精品成人一区二区三区 | 啦啦啦免费观看视频1| 制服诱惑二区| 免费在线观看视频国产中文字幕亚洲| 黄片小视频在线播放| 国产区一区二久久| 在线观看www视频免费| 日韩视频一区二区在线观看| 国产精品一区二区在线不卡| 亚洲少妇的诱惑av| 嫁个100分男人电影在线观看| 午夜91福利影院| 91麻豆精品激情在线观看国产 | 亚洲av成人不卡在线观看播放网| 人人妻人人爽人人添夜夜欢视频| 69精品国产乱码久久久| 美国免费a级毛片| 国产精品久久久人人做人人爽| 久久午夜亚洲精品久久| 香蕉久久夜色| 91麻豆精品激情在线观看国产 | 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 成人三级做爰电影| 亚洲精品在线美女| 一个人免费在线观看的高清视频| 法律面前人人平等表现在哪些方面| 飞空精品影院首页| 亚洲欧美一区二区三区黑人| 欧美国产精品一级二级三级| 成人av一区二区三区在线看| 国产有黄有色有爽视频| 国产人伦9x9x在线观看| www.自偷自拍.com| 亚洲九九香蕉| 91大片在线观看| 国产在线一区二区三区精| 大型av网站在线播放| 午夜免费鲁丝| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| 国产精品久久电影中文字幕 | 我的亚洲天堂| 午夜免费成人在线视频| 超色免费av| 欧美性长视频在线观看| 啦啦啦 在线观看视频| 法律面前人人平等表现在哪些方面| 国产深夜福利视频在线观看| 大香蕉久久成人网| 女性被躁到高潮视频| 国产亚洲欧美在线一区二区| 黄色丝袜av网址大全| 亚洲成人免费av在线播放| 欧美精品一区二区免费开放| 曰老女人黄片| 一个人免费在线观看的高清视频| 国产aⅴ精品一区二区三区波| 亚洲欧美色中文字幕在线| 99热国产这里只有精品6| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 久久热在线av| 日韩精品免费视频一区二区三区| 一区二区av电影网| 少妇 在线观看| 少妇裸体淫交视频免费看高清 | 亚洲第一av免费看| 国产欧美日韩综合在线一区二区| 亚洲精品一卡2卡三卡4卡5卡| 精品国产一区二区三区久久久樱花| 午夜免费鲁丝| 少妇被粗大的猛进出69影院| 国产97色在线日韩免费| 女人被躁到高潮嗷嗷叫费观| 又大又爽又粗| 在线 av 中文字幕| 国产成人精品在线电影| 国产精品 欧美亚洲| 日韩成人在线观看一区二区三区| 精品熟女少妇八av免费久了| 午夜两性在线视频| 国产成人精品无人区| 咕卡用的链子| 日韩三级视频一区二区三区| 黄色丝袜av网址大全| 免费少妇av软件| 色婷婷av一区二区三区视频| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 国产在视频线精品| 国产亚洲欧美在线一区二区| 久久久久久久久免费视频了| 两个人免费观看高清视频| 欧美黄色淫秽网站| 男人舔女人的私密视频| 青草久久国产| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯 | 在线观看免费视频日本深夜| 极品教师在线免费播放| 精品国内亚洲2022精品成人 | 久久青草综合色| 亚洲视频免费观看视频| 老司机午夜福利在线观看视频 | 啪啪无遮挡十八禁网站| 午夜日韩欧美国产| 十八禁网站免费在线| 成人影院久久| 一区二区日韩欧美中文字幕| 成年女人毛片免费观看观看9 | 亚洲国产成人一精品久久久| a在线观看视频网站| 亚洲色图av天堂| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 天堂动漫精品| 中文字幕最新亚洲高清| 午夜福利乱码中文字幕| 老司机午夜福利在线观看视频 | 亚洲精华国产精华精| 老司机深夜福利视频在线观看| 露出奶头的视频| 啦啦啦在线免费观看视频4| 国产伦理片在线播放av一区| 高清在线国产一区| 午夜福利乱码中文字幕| 精品欧美一区二区三区在线| 人人妻,人人澡人人爽秒播| 久9热在线精品视频| 精品国产亚洲在线| 久久久久国内视频| 国产精品1区2区在线观看. | 香蕉丝袜av| 色综合婷婷激情| 美女高潮到喷水免费观看| 国产日韩一区二区三区精品不卡| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 午夜福利视频在线观看免费| 久久这里只有精品19| 亚洲 国产 在线| 99热国产这里只有精品6| 亚洲国产欧美日韩在线播放| 成人手机av| 亚洲熟女毛片儿| 亚洲一区二区三区欧美精品| 色老头精品视频在线观看| 精品免费久久久久久久清纯 | 亚洲国产欧美日韩在线播放| 中文欧美无线码| 亚洲欧美精品综合一区二区三区| 国产福利在线免费观看视频| 亚洲自偷自拍图片 自拍| 香蕉丝袜av| 女人爽到高潮嗷嗷叫在线视频| 国产精品一区二区在线观看99| 免费一级毛片在线播放高清视频 | 欧美日韩亚洲国产一区二区在线观看 | www.精华液| 日韩中文字幕欧美一区二区| 夜夜爽天天搞| 久久精品亚洲熟妇少妇任你| 久久人妻福利社区极品人妻图片| 精品卡一卡二卡四卡免费| 性高湖久久久久久久久免费观看| 久久 成人 亚洲| 国产精品九九99| 91精品国产国语对白视频| 国产欧美亚洲国产| 国产99久久九九免费精品| 男女床上黄色一级片免费看| 黑人欧美特级aaaaaa片| 国产免费av片在线观看野外av| 久久热在线av| 伊人久久大香线蕉亚洲五| 亚洲人成电影免费在线| 免费在线观看影片大全网站| 麻豆国产av国片精品| 老汉色∧v一级毛片| 亚洲天堂av无毛| 久久婷婷成人综合色麻豆| 后天国语完整版免费观看| 亚洲精品中文字幕在线视频| 制服诱惑二区| 免费在线观看影片大全网站| 麻豆国产av国片精品| 成人av一区二区三区在线看| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久久久精品古装| 青草久久国产| 午夜两性在线视频| 在线观看免费日韩欧美大片| 日本黄色视频三级网站网址 | 久久精品国产99精品国产亚洲性色 | 狠狠精品人妻久久久久久综合| 在线观看免费日韩欧美大片| 欧美日韩黄片免| 欧美人与性动交α欧美精品济南到| 不卡一级毛片| 午夜精品久久久久久毛片777| 成人av一区二区三区在线看| 在线观看免费午夜福利视频| 中文字幕制服av| 热99re8久久精品国产| 99国产极品粉嫩在线观看| 国产成人免费观看mmmm| 在线观看人妻少妇| 狂野欧美激情性xxxx| 日本欧美视频一区| 国产激情久久老熟女| 久久毛片免费看一区二区三区| 制服人妻中文乱码| 精品乱码久久久久久99久播| 日本av手机在线免费观看| 成年人免费黄色播放视频| 在线观看免费高清a一片| 国产亚洲欧美在线一区二区| 精品卡一卡二卡四卡免费| 国产有黄有色有爽视频| 五月天丁香电影| 久久久久久久久久久久大奶| 亚洲精品成人av观看孕妇| 国产色视频综合| 欧美 日韩 精品 国产| 三级毛片av免费| videos熟女内射| 亚洲国产欧美在线一区| 成年动漫av网址| 黑人猛操日本美女一级片| 黄片小视频在线播放| 亚洲欧美日韩高清在线视频 | 亚洲欧美一区二区三区久久| 极品教师在线免费播放| 免费日韩欧美在线观看| 国产成人精品久久二区二区91| 一级a爱视频在线免费观看| 免费人妻精品一区二区三区视频| 国产成人欧美| 国产一区二区在线观看av| 18禁观看日本| 夜夜夜夜夜久久久久| 日韩 欧美 亚洲 中文字幕| 伦理电影免费视频| 欧美黑人欧美精品刺激| 黄色毛片三级朝国网站| 国产精品自产拍在线观看55亚洲 | 在线亚洲精品国产二区图片欧美| bbb黄色大片| videos熟女内射| 亚洲国产欧美在线一区| 国产av一区二区精品久久| 国产男女内射视频| 国产精品国产av在线观看| 色视频在线一区二区三区| 国产在线视频一区二区| 久久国产亚洲av麻豆专区| 一二三四社区在线视频社区8| 中文字幕最新亚洲高清| www.999成人在线观看| 国产1区2区3区精品| 天天添夜夜摸| 成年人午夜在线观看视频| 色精品久久人妻99蜜桃| 精品国产亚洲在线| 18禁黄网站禁片午夜丰满| 亚洲美女黄片视频| 午夜久久久在线观看| 美女高潮喷水抽搐中文字幕| 国产成人免费无遮挡视频| 日韩欧美免费精品| 午夜激情av网站| 国产av又大| 老熟女久久久| 91麻豆av在线| 18禁国产床啪视频网站| 色视频在线一区二区三区| 亚洲人成77777在线视频| 757午夜福利合集在线观看| a级毛片在线看网站| 国产视频一区二区在线看| 中文字幕人妻熟女乱码| 国产日韩欧美亚洲二区| 午夜久久久在线观看| 亚洲久久久国产精品| 亚洲午夜理论影院| 亚洲av欧美aⅴ国产| 国产精品久久久久久精品古装| 亚洲国产av影院在线观看| 国产xxxxx性猛交| 最新美女视频免费是黄的| 国产精品九九99| av不卡在线播放| 麻豆av在线久日| 老熟妇仑乱视频hdxx| 国产精品亚洲av一区麻豆| 在线观看免费视频网站a站| 国产成人av激情在线播放| 久久国产精品大桥未久av| 曰老女人黄片| 国产黄色免费在线视频| 亚洲专区中文字幕在线| 免费不卡黄色视频| 国产高清激情床上av| 精品国内亚洲2022精品成人 | 久久精品人人爽人人爽视色| 在线十欧美十亚洲十日本专区| 久久精品成人免费网站| 亚洲专区中文字幕在线| 精品午夜福利视频在线观看一区 | 亚洲成a人片在线一区二区| 欧美黑人精品巨大| 国产av国产精品国产| 精品免费久久久久久久清纯 | 丝袜喷水一区| 国产精品欧美亚洲77777| 国产欧美亚洲国产| 国产精品熟女久久久久浪| 久久九九热精品免费| www.自偷自拍.com| 欧美一级毛片孕妇| 男女床上黄色一级片免费看| 精品久久久精品久久久| a级毛片在线看网站| 三上悠亚av全集在线观看| 色综合欧美亚洲国产小说| 99re6热这里在线精品视频| 美国免费a级毛片| 十八禁人妻一区二区| 99国产极品粉嫩在线观看| svipshipincom国产片| 欧美激情极品国产一区二区三区| 亚洲专区字幕在线| 亚洲一区二区三区欧美精品| 91字幕亚洲| √禁漫天堂资源中文www| 亚洲精品国产精品久久久不卡| 一区二区三区乱码不卡18| 久久人妻福利社区极品人妻图片| 曰老女人黄片| 国产欧美日韩一区二区精品| 亚洲美女黄片视频| 最新美女视频免费是黄的| 99精国产麻豆久久婷婷| 国产精品亚洲一级av第二区| 国产麻豆69| 久久久精品免费免费高清| 窝窝影院91人妻| 人成视频在线观看免费观看| 女人被躁到高潮嗷嗷叫费观| 欧美成人免费av一区二区三区 | 国产不卡一卡二| 99精品在免费线老司机午夜| 在线看a的网站| 飞空精品影院首页| 老司机福利观看|