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

    基于全基因序列的中國北方3省(區(qū))馬鈴薯Y病毒遺傳多樣性分析

    2020-08-25 10:03:25馬俊豐李小宇張春雨周雪平王永志
    植物保護 2020年4期
    關鍵詞:遺傳多樣性

    馬俊豐 李小宇 張春雨 周雪平 王永志

    摘要 本研究以馬鈴薯Y病毒(PVY)全基因組為基礎,分析吉林、黑龍江和內蒙古3?。▍^(qū))PVY群體遺傳多樣性和群體分化,并評估突變、重組、選擇等遺傳力所起的作用。根據已報道的PVY全基因序列保守區(qū)設計4對引物,采用片段重疊法對來自內蒙古和吉林的24個PVY分離物全基因序列進行測定,并聯(lián)合NCBI中已登錄的9個黑龍江分離物全基因組序列進行遺傳多樣性參數(shù)評估、群體分化檢驗和分子變異等分析。結果顯示,我國北方3?。▍^(qū))PVY群體遺傳多樣性高,其中內蒙古和黑龍江PVY群體遺傳多樣性高于吉林群體,并且3個群體之間呈現(xiàn)一定程度的遺傳分化。分子變異分析發(fā)現(xiàn)在PVY基因組中存在1 786個變異位點,表明我國北方3?。▍^(qū))PVY群體變異程度較高,并且這種高變異度有85.54%來自各個馬鈴薯種植區(qū)內PVY個體的遺傳變異。重組分析和系統(tǒng)發(fā)育分析發(fā)現(xiàn),我國北方3?。▍^(qū))PVY群體中重組株系占比高達90.3%,并具有明顯的株系多樣性,表明PVY重組株系已成為我國北方3?。▍^(qū))馬鈴薯種植區(qū)的流行株系。選擇壓力分析顯示,使用FEL和IFEL法分別檢測出501個和315個凈化壓力選擇位點,這表明3?。▍^(qū))PVY群體受凈化選擇壓力為主。以上結果表明,中國北方3?。▍^(qū))PVY群體遺傳多樣性高,突變、重組和自然選擇都對遺傳多樣性和群體分化存在一定影響。

    關鍵詞 馬鈴薯Y病毒(PVY); 遺傳多樣性; 基因組序列

    中圖分類號:

    S 432.41

    文獻標識碼: A

    DOI: 10.16688/j.zwbh.2019571

    Genetic diversity analysis of Potato virus Y in three provinces of north China based on the full genomic sequences

    MA Junfeng1, LI Xiaoyu1, ZHANG Chunyu1, ZHOU Xueping2, WANG Yongzhi1, 2*

    (1. Institute of Plant Protection, Jilin Academy of Agricultural Sciences, Gongzhuling 136100, China;

    2. Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China)

    Abstract

    Complete genome sequences of Potato virus Y (PVY) from three provinces (Jilin, Heilongjiang and Inner Mongolia) of north China were used to investigate its genetic diversity and population differentiation. Meanwhile, the effects of mutations, recombination and selection on population genetic dynamics were evaluated. The complete genomes of 24 PVY isolates from Jilin and Inner Mongolia were amplified and sequenced using 4 pairs of primers designed according to the conserved regions of the reported PVY isolates, and were then combined with complete genome sequences of 9 Heilongjiang isolates for estimation of genetic variation, population differentiation, molecular variation, etc. The genetic diversity of the PVY populations in the three provinces was high. Furthermore, the genetic diversities of PVY isolates from Heilongjiang and Inner Mongolia were higher than that of PVY from Jilin and there was a certain degree of genetic differentiation among the three populations. Analysis of molecular variation revealed that there were 1 786 mutation sites in PVY genomes, which indicated that the degree of variation in PVY populations in the three provinces was high; 85.54% of the variation was from individuals of different potato planting areas. Recombination and phylogenetic analysis showed that 90.3% strains were recombinant and there was obvious strain diversity in the PVY populations, indicating that the recombinant strains had become the epidemic strain in the potato cropping regions of the three provinces. Totally 501 and 315 purifying selection sites were detected by using FEL and IFEL methods, respectively, suggesting that the PVY populations in the three provinces were mainly affected by purifying selection pressure. In summary, the PVY populations in the three provinces of northern China had high genetic diversity, and mutations, recombination and selection all contributed to the genetic diversity and population differentiation.

    Key words

    Potato virus Y(PVY); genetic diversity; genome sequence

    馬鈴薯病毒是危害馬鈴薯的主要病原物之一,已有報道顯示,世界范圍內能夠侵染馬鈴薯的病毒至少有53種,其中10種已在我國被發(fā)現(xiàn)[1]。在我國已經發(fā)現(xiàn)的10種病毒中,馬鈴薯Y病毒Potato virus Y(PVY)發(fā)生范圍最廣,造成經濟損失最為嚴重。近些年在我國多個馬鈴薯種植區(qū)檢測到PVY[2-7]。PVY侵染馬鈴薯可引起葉脈壞死、葉面黃化褪綠、塊莖環(huán)斑壞死等癥狀,導致馬鈴薯種質退化,產量降低,給馬鈴薯的生產造成嚴重損失[8-10]。

    PVY是馬鈴薯Y病毒科Potyviridae馬鈴薯Y病毒屬Potyvirus代表成員,其病毒粒體呈彎曲線狀,每個病毒粒體包含一套完整的基因組。該基因組由正義單鏈RNA構成,長約9.7 kb,包含一個大的開放閱讀框(open reading frame, ORF),兩端含有非編碼區(qū)(untranslated regions,UTRs)?;蚪M翻譯時,開放閱讀框先編碼一個多聚蛋白(polyprotein),隨后再通過自身編碼的3個蛋白酶將多聚蛋白切割加工為P1、HC-Pro、P3等10個成熟蛋白[11]。除此之外,PVY基因組內還有一個通過移碼翻譯產生的蛋白——PIPO,該蛋白與P3蛋白的N′端以P3N-PIPO融合形式存在[12]。最終,PVY基因組通過這兩種方式生成11個成熟的多功能蛋白,使整個基因組呈現(xiàn)不同程度的遺傳多樣性。

    近些年,PVY呈高變異、多株系發(fā)展趨勢,因此,PVY研究已不局限在分子檢測、株系鑒定、蛋白質功能驗證等方面[13-17],遺傳變異及分子進化已成為新的研究方向。已有工作主要是選擇其中1~2個功能基因研究PVY的遺傳多樣性[18-20],不能全面地反映PVY進化和遺傳機制。因此,本研究通過擴增、克隆等手段獲得吉林(JL)和內蒙古(NMG)PVY分離物全基因組序列24條,并聯(lián)合GenBank中已提交的9條黑龍江(HLJ)分離物的全基因組序列進行序列比對分析、遺傳多樣性參數(shù)評估及群體分化檢驗等一系列分析,旨在基于全基因組序列更為準確地探明我國北方3省(區(qū))PVY群體遺傳多樣性及進化機制。該研究結果將有助于了解我國北方地區(qū)PVY的發(fā)生、流行及變異趨勢,為其有效防控提供理論依據。

    1 材料與方法

    1.1 毒源

    2017年—2018年從吉林、內蒙古馬鈴薯種植區(qū),用隨機抽樣法采集帶有花葉、褪綠和壞死等典型病毒癥狀的馬鈴薯葉片。先采用酶聯(lián)免疫吸附法(enzyme linked immunosorbent assay, ELISA)檢驗病毒,之后對病毒CP基因進行擴增,將確定為PVY分離物的葉片密封保存在-80℃超低溫冰箱中。

    1.2 主要試劑

    RNeasy Plant Mini Kit,QIAGEN公司;ReverAid First Strand cDNA Synthesis Kit,賽默飛世爾科技公司(Thermo Fisher Scientific);TaKaRa LA Taq、pMD18-T Vector Cloning Kit、DNA Marker,寶生物工程有限公司;AxyPrep DNA Gel Extraction Kit,愛思進生物技術公司(Axygen)。

    1.3 試驗方法

    1.3.1 RNA提取及基因克隆

    按照RNeasy Plant Mini Kit說明書,從感染PVY的馬鈴薯葉片中提取總RNA。以提取的葉片總RNA為模板,以OligodT為引物,按照反轉錄試劑盒說明書合成cDNA第一鏈。根據GenBank中已提交的PVY不同分離物基因組序列的保守區(qū),設計用于擴增PVY全基因序列的簡并引物和特異性引物(表1)。引物委托吉林省庫美生物科技有限公司合成。

    PCR擴增采用50 μL反應體系:TaKaRa LATaq (5 U/L)0.5 μL、10×LA PCR buffer Ⅱ(Mg2+Plus)5 μL、dNTP mixture 5 μL、cDNA 1 μL、正、反向引物(10 μmol/L)各 1 μL、超純水36.5 μL。PCR反應條件:94℃預變性5 min;94℃變性30 s,退火30 s(退火溫度視基因片段而定),72℃延伸若干分鐘(延伸時間按1 kb/min計算),共30個循環(huán);最后72℃再延伸10 min。

    取4 μL PCR產物用1.5%瓊脂糖凝膠電泳進行檢測。檢測后切取目的條帶,利用膠回收試劑盒進行純化,純化產物分別與pMD18-T載體連接,連接產物轉化到大腸桿菌DH5α感受態(tài)細胞中。經菌液PCR驗證后篩選出陽性克隆,隨機選取3~5個陽性克隆,委托吉林省庫美生物科技有限公司進行測序,并通過測序峰圖及序列比對排除PCR過程中可能引起的突變。

    1.3.2 序列處理

    測序獲得的吉林和內蒙古PVY分離物各序列片段使用DNAMAN 6.0軟件進行拼接,并通過PVY參考序列(GenBank登錄號為NC_001616)進行校正。校正時先使用MEGA 6.0軟件中的Muscle程序對PVY全基因組序列進行多重比對,之后使用Muscle(codon)程序對PVY基因組編碼區(qū)進行比對,最后將編碼區(qū)翻譯為相應的氨基酸序列,進行人工校正,以確?;蚪M序列的準確性。試驗所用的黑龍江PVY分離物基因組序列均為NCBI中已登錄的序列(序列名稱及登錄號見圖1)。最終用于本試驗研究的PVY全基因組序列共33條。

    1.3.3 遺傳多樣性參數(shù)評估及群體分化檢驗

    使用DnaSP 5.10軟件計算遺傳多樣性參數(shù):多態(tài)性位點數(shù)目(segregating sits, S)、單倍型數(shù)目(number of haplotypes, h)、單倍型多樣性(haplotype diversity, Hd)、核苷酸差異平均數(shù)(average number of nucleotide differences, k)、核苷酸多樣性(nucleotide diversity, Pi)。使用Arlequin 3.5軟件計算群體分化系數(shù)(fixation index, Fst)并根據Fst值的大小判斷群體間分化程度。當00.25時,群體嚴重分化。群體間遺傳分化顯著水平采用Kst、Z、Snn[21-22]3種獨立統(tǒng)計方法進行檢驗,并通過排列檢驗法(1 000 permutations)進行判斷。

    1.3.4 PVY基因組的分子變異

    以PVY全基因組序列為分子標記,利用Arlequin 3.5軟件進行不同群體的分子變異分析(analysis of molecular variance, AMOVA)。使用DnaSP 5.10軟件進行變異位點分析。在網站http:∥blast.ncbi.nlm.nih.gov/上進行序列一致性分析。

    1.3.5 PVY基因組的重組

    以PVYO株系分離物Oz(EF026074)和PVYN株系分離物Mont(AY884983)作為親本,使用SimPlot 3.5軟件進行重組分析。同時在網站(http:∥www.datamonkey.org)使用遺傳算法中的重組檢測法(genetic algorithm recombination detection,GARD)[23]進行重組位點檢測及可靠性評估。

    1.3.6 壓力分析

    使用在線網站(http:∥www.datamonkey.org),基于密碼子,通過IFEL(internal branches fixed-effects likelihood)、FEL(fixed effects likelihood)、MEME(mixed effects model of evolution)3種不同方法,對PVY基因組進行壓力檢測及統(tǒng)計檢驗。壓力選擇方向根據非同義核苷酸替換速率(dN)與同義核苷酸替換速率(dS)比值(ω=dN/dS)大小進行判斷。當ω>1且群體差異顯著時,說明基因組受到正向選擇(positive selection)作用;當0<ω<1時,說明基因組受到凈化選擇(purifying selection)作用;當ω=1時,基因組不受壓力選擇[24]。

    1.3.7 系統(tǒng)發(fā)育分析

    使用最大似然法(maximum likelihood, ML)基于全基因組序列編碼區(qū)(coding sequence, CDs)核苷酸序列構建PVY系統(tǒng)發(fā)育樹。建樹前,使用MAFFT軟件對建樹的33條PVY編碼區(qū)核苷酸序列進行多重比對,并使用DAMBE軟件檢測核苷酸替換是否飽和,之后利用MEGA 6.0軟件選擇最優(yōu)化的核苷酸替換模型,并參照BIC(Bayesian information criterion)標準設置相應參數(shù),最后通過1 000次自舉(bootstrap)抽樣來對各分支節(jié)點的置信值進行評估。

    2 結果與分析

    2.1 全基因組序列特征

    本研究共獲得24個PVY分離物全基因組序列(9 698 nt或9 699 nt),5′非編碼區(qū)長度為184 nt或185 nt,3′非編碼區(qū)長度為328 nt,編碼區(qū)長度為9 186 nt;其中在編碼區(qū)部分,大開放閱讀框對應的多聚蛋白編碼序列長度為9 186 nt,編碼一個包含3 061個氨基酸殘基的多聚蛋白;小開放閱讀框對應的蛋白PIPO編碼序列長度為195 nt或228 nt,并通過移碼翻譯編碼一個包含65或76個氨基酸殘基的蛋白PIPO(表2)。

    2.2 遺傳多樣性參數(shù)評估及群體分化檢驗

    PVY基因組遺傳多樣性參數(shù)如表3所示:3個群體的多態(tài)性位點數(shù)目(S)表現(xiàn)出一定差異,其中黑龍江和內蒙古群體S值相對較大,多態(tài)性位點數(shù)目分別為1 232個和1 229個,吉林群體相對較小,多態(tài)性位點數(shù)目為807個;核苷酸差異平均數(shù)(k)方面,黑龍江和內蒙古群體k值較大,分別為491.361和526.933,吉林群體k值較小為249.659;單倍型多樣性(Hd)方面,3個群體Hd值均為1,大于0.5,且33個分離物中共發(fā)現(xiàn)33種單倍型(h),未出現(xiàn)兩條完全一致的序列,各個群體均表現(xiàn)出較高的單倍型多樣性;核苷酸多樣性方面(Pi),3個群體的Pi值分別為0.050 61、0.025 72、0.054 29,均大于0.005,表明3個群體均具有較高的核苷酸多樣性。以上數(shù)據表明我國北方3省(區(qū))PVY群體遺傳多樣性高,且黑龍江和內蒙古PVY群體遺傳多樣性高于吉林群體。

    以PVY基因組為分子標記,黑龍江和吉林PVY群體之間Fst值為0.051 71,Z統(tǒng)計顯著,Kst和Snn值統(tǒng)計不顯著,表明在黑龍江和吉林種植區(qū)的PVY群體間呈現(xiàn)中度分化但差異不顯著;黑龍江和內蒙古,群體之間Fst值為0.072 52,Z、Kst和Snn統(tǒng)計均不顯著,表明黑龍江和內蒙古種植區(qū)的PVY群體呈現(xiàn)中度分化,但差異不顯著;吉林和內蒙古群體之間Fst值為0.233 52,Kst統(tǒng)計顯著,Z統(tǒng)計極顯著,Snn統(tǒng)計不顯著,表明吉林和內蒙古種植區(qū)的PVY群體呈現(xiàn)高度分化且差異顯著(表4)。

    2.3 分子變異分析

    全基因組序列一致性比對分析結果顯示,本研究獲得的24個分離物之間核苷酸一致性為90.3%~99.0%,氨基酸一致性為99.0%~100%,進一步與NCBI中已登錄的9個黑龍江PVY分離物比對發(fā)現(xiàn),核苷酸一致性在89.3%~99.5%之間,氨基酸一致性在98.6%~100%之間。綜上所述,本次研究3個群體之間核苷酸一致性為89.3%~99.5%,氨基酸一致性為98.6%~100%。

    突變位點分析發(fā)現(xiàn),3個群體PVY全基因組多態(tài)性位點(polymorphic sites)共1 786個,其中單一變異位點(singleton variable sites)450個,簡約信息位點(parsimony informative sites)1 336個,表明PVY基因組存在較高程度的變異。

    以PVY基因組為基礎,AMOVA分析結果如表5所示,14.46%的變異來自群體間,85.54%的變異來自群體內,因此不同種植區(qū)PVY群體之間的遺傳變異主要來自各個種植區(qū)內PVY個體的遺傳變異。

    2.4 重組分析

    使用Simplot軟件分析發(fā)現(xiàn),33個PVY分離物中有30個檢測到潛在重組信號,約占分離物總數(shù)的90.3%,并且所有檢測到的重組分離物均為N×O重組類型。吉林和黑龍江PVY群體分離物均存在重組,而內蒙古PVY群體存在3個O株系非重組分離物。進一步通過GARD程序確認,3個群體存在4個或5個重組位點,每個重組位點都具有較高的模型支持率,且KH檢驗均達到顯著水平(表6)。

    2.5 選擇壓力分析

    PVY全基因組選擇壓力分析顯示(表7),使用FEL、IFEL、MEME 3種檢測方法均檢測到大量選擇壓力位點,在3個群體中利用FEL和IFEL法分別檢測出501個和315個凈化選擇位點,但僅檢測出9個和12個正向選擇位點;使用MEME法檢測出43個正向選擇位點。而以上3種方法共同檢測出6個正向選擇位點,且這些位點均顯著(表8)。以上結果表明盡管PVY基因組受凈化選擇壓力為主,但基因組中也存在個別密碼子位點受到強烈的正向選擇壓力。

    2.6 系統(tǒng)發(fā)育分析

    建樹序列多重比對后長度為9 186 nt,根據MEGA 6.0軟件的BIC標準,建樹序列最合適的核苷酸替換模型為GTR+I+G。ML法構建的系統(tǒng)發(fā)育樹如圖1所示。從圖中可以看出3?。▍^(qū))PVY分離物并未以地區(qū)為標準相聚成簇,因此北方3省(區(qū))PVY分離物在系統(tǒng)發(fā)育關系上與地理分布無顯著聯(lián)系。根據Chang等[25]的研究結果,以PVY基因組編碼區(qū)構建系統(tǒng)發(fā)育樹,相同株系的分離物會相聚成簇,表現(xiàn)出明顯的株系特異性,因此,本研究中33個PVY分離物是以不同株系類型為標準相聚成簇,由此可以看出北方3?。▍^(qū))的PVY群體具有較為明顯的株系多樣性。

    3 討論

    遺傳多樣性是體現(xiàn)病原物群體對環(huán)境適應力的一項重要指標,是突變、基因重組、基因流(基因遷移)、隨機漂變和自然選擇5種遺傳力長期相互作用的結果。遺傳多樣性直接影響著病原物群體的進化潛力。高遺傳多樣性的病原物具有更強的生存和進化優(yōu)勢,能更快適應新的寄主和生存環(huán)境[26]。開展病原物遺傳多樣性及群體分化研究有助于了解植物病原物的發(fā)生規(guī)律,流行機制,并為相關病害的監(jiān)測及防控提供理論依據。

    中國北方3省(區(qū))PVY遺傳多樣性參數(shù)結果表明,黑龍江、吉林和內蒙古3?。▍^(qū))PVY群體遺傳多樣性高,與前人研究結果一致,符合中國PVY群體遺傳多樣性高的特征[27-28]。3?。▍^(qū))的馬鈴薯種植區(qū)存在一定程度的種薯調運,這會導致不同種植區(qū)間PVY群體的二次接觸,從而造成高遺傳多樣性。此外,種群內這種高遺傳多樣性的形成是一個長期過程,本次研究的3個?。▍^(qū))PVY群體基因組單倍型多樣性(Hd)均大于0.5,核苷酸多樣性(Pi)均大于0.005,且兩兩種植區(qū)之間PVY群體呈現(xiàn)一定程度的群體分化,表明3?。▍^(qū))PVY群體是由一個相對較大且穩(wěn)定的群體經過長期演化所形成的。

    PVY具有高度變異性,而突變在PVY的遺傳變異過程中起到重要作用。PVY是RNA病毒,其依賴于RNA的RNA聚合酶(RNA-dependent RNA polymerase, RdRp)缺乏校正功能,會在復制過程中產生大量的點突變[29],在所產生的突變中有利突變會被逐漸積累,而有害突變則在自然選擇過程中被逐漸淘汰,由此突變常被視作推動PVY進化的原始動力。本試驗基于PVY全基因序列發(fā)現(xiàn)了大量變異位點,表明我國北方3?。▍^(qū))PVY群體存在較高程度的變異。AMOVA分析結果顯示3?。▍^(qū))PVY群體之間的遺傳變異主要來自各個種植區(qū)內PVY個體的遺傳變異。在序列一致性方面,3個PVY群體之間的核苷酸一致性為89.3%~99.5%,符合同一種病毒全基因序列一致性不低于85.0%的標準[11]。

    基因重組同樣是PVY遺傳變異的主要驅動力,在PVY群體中發(fā)揮著重要的作用。已經確認的重組位點主要集中在P1、HC-Pro/P3、VPg、NIb和CP基因區(qū)域[30-32]。本試驗3個PVY群體中均檢測出大量的重組分離物,并且3個PVY群體中都存在著多個重組位點,表明相較于基因突變,基因重組可能是PVY進化過程中更為重要的機制,因為重組產生的遺傳變異遠比僅由突變造成的變異更快、更頻繁[33]。常飛等在進行PVY群體遺傳分析時并未檢測到重組信號[28],原因是選用的P3和PIPO基因位于PVY基因組的非重組區(qū)域,因此基于全基因序列的重組分析可以更準確的了解PVY群體的重組狀況。此外PVY N×O重組分離物的大量發(fā)現(xiàn)也表明此類型重組株系已成為北方3省(區(qū))馬鈴薯種植區(qū)的主流株系。而PVY重組株系的頻繁發(fā)生可以歸納于以下兩點:一是重組分離物相較于非重組分離物有著更強的致病性;二是重組分離物兼具非重組分離物的適應性優(yōu)勢,可以快速適應環(huán)境變化,并發(fā)展為優(yōu)勢株系[10, 34]。因此,在生產上需密切關注重組株系的流行趨勢。

    選擇壓力分析中,中國北方3省(區(qū))PVY群體基因組凈化選擇位點遠高于正向選擇位點,表明PVY基因組中多數(shù)位點的突變是有害的,在進化過程中這些有害突變會被剔除,從而保證PVY各個基因發(fā)揮正常功能。而采用3種方法共同檢測到6個顯著的正向選擇壓力位點,表明PVY基因組中存在個別位點的突變是有利于該病毒生存競爭的,在遺傳進化過程中這些有利突變會被不斷積累。

    CP和PIPO基因由于高度保守的特點,在進行Potyvirus病毒系統(tǒng)發(fā)育分析時,常被用以構建系統(tǒng)發(fā)育樹[12],但對于PVY等變異程度高,重組頻繁的病毒來說,單基因建樹不足以準確地反映同種病毒不同株系分離物之間的系統(tǒng)發(fā)育關系,因此,以全基因組構建系統(tǒng)發(fā)育樹成為解決此類問題的有效方法[25,35]。本研究以我國北方3?。▍^(qū))PVY分離物全基因序列進行系統(tǒng)發(fā)育分析,準確地反映出了PVY分離物的系統(tǒng)發(fā)育關系,為我國北方3?。▍^(qū))PVY群體遺傳多樣性研究提供了準確數(shù)據。

    參考文獻

    [1] 路秉翰, 湯亞飛, 于琳, 等. 基于小RNA深度測序和RT-PCR檢測侵染廣東省冬種馬鈴薯的病毒[J]. 植物保護學報, 2019, 46(3): 604-610.

    [2] 劉洪義, 張洪祥, 李明福, 等. 黑龍江省馬鈴薯病毒病的普查及鑒定[J]. 東北農業(yè)大學學報, 2006, 37(3): 307-310.

    [3] 鐘婷婷, 蒲志剛, 何俊蓉, 等. 四川省馬鈴薯主產區(qū)最新病毒病普查及血清學鑒定[J]. 西南農業(yè)學報, 2008, 21(1): 96-99.

    [4] 胡新喜, 雷艷, 何長征, 等. 湖南省馬鈴薯主產區(qū)馬鈴薯病毒種類及流行分析[J]. 中國馬鈴薯, 2012, 26(6): 358-361.

    [5] 黃丹, 余琨, 陳建斌, 等. 馬鈴薯病毒PVY、PVS和PLRV多重RT-PCR檢測[J]. 云南農業(yè)大學學報, 2015, 30(4): 535-540.

    [6] 苑智華. 烏蘭察布市馬鈴薯病毒病調查分析[J]. 江蘇農業(yè)科學, 2015, 43(11): 189-191.

    [7] 王炳森, 王麗, 陳惠蘭. 山東省馬鈴薯病害調查[J]. 中國馬鈴薯, 2017, 31(2): 104-112.

    [8] NOLTE P, WHITWORTH J L, THORNTON M K, et al. Effect of seedborne Potato virus Y on performance of Russet Burbank, Russet Norkotah, and Shepody potato [J]. Plant Disease, 2004, 88(3): 248-252.

    [9] RAHMAN M S, AKANDA A K. Performance of seed potato produced from sprout cutting, stem cutting and conventional tuber against PVY and PLRV [J]. Bangladesh Journal of Agricultural Research, 2009, 34(4): 609-622.

    [10]CHIKH A M, MAOKA T, NATSUAKI T, et al. PVYNTN-NW, a novel recombinant strain of Potato virus Y predominating in potato fields in Syria [J]. Plant Pathology, 2010, 59(1): 31-41.

    [11]KING A M, LEFKOWITZ E, ADAMS M J, et al. Virus taxonomy, ninth report of the international committee on taxonomy of viruses [M]. Amsterdam: Elsevier Academic Press, 2012.

    [12]高芳鑾, 沈建國, 史鳳陽, 等. 馬鈴薯Y病毒pipo基因的分子變異及結構特征分析[J]. 遺傳, 2013, 35(9): 1125-1134.

    [13]LORENZEN J H, MEACHAM T, BERGER P H, et al. Whole genome characterization of Potato virus Y isolates collected in the western USA and their comparison to isolates from Europe and Canada [J]. Archives of Virology, 2006, 151(6): 1055-1074.

    [14]CHIKH A M, GRAY S M, KARASEV A V. An improved multiplex IC-RT-PCR assay distinguishes nine strains of Potato virus Y [J]. Plant Disease, 2013, 97(10): 1370-1374.

    [15]GAO Fangluan, CHANG Fei, SHEN Jianguo, et al. Complete genome analysis of a novel recombinant isolate of Potato virus Y from China [J]. Archives of Virology, 2014, 159(12): 3439-3442.

    [16]GLAIS L, FAUREZ F, TRIBODET M, et al. The amino acid 419 in HC-Pro is involved in the ability of PVY isolate N605 to induce necrotic symptoms on potato tubers [J]. Virus Research, 2015, 208: 110-119.

    [17]KIM J, CHA D J, KWON M, et al. Potato virus Y (PVY) detection in a single aphid by one-step RT-PCR with boiling technique [J]. Entomological Research, 2016, 46(4): 278-285.

    [18]TIAN Yanping, LIU Jinliang, ZHANG Chengling, et al. Genetic diversity of Potato virus Y infecting tobacco crops in China [J]. Phytopathology, 2011, 101(3): 377-387.

    [19]GAO Fangluan, ZOU Wenchao, XIE Lianhui, et al. Adaptive evolution and demographic history contribute to the divergent population genetic structure of Potato virus Y between China and Japan [J]. Evolutionary Applications, 2017, 10(4): 379-390.

    [20]MAO Yanzhi, SUN Xuhong, SHEN Jianguo, et al. Molecular evolutionary analysis of Potato Virus Y infecting potato based on the VPg gene [J/OL]. Frontiers in Microbiology, 2019, 10: 1708. DOI:10.3389/fmicb.2019.01708.

    [21]HUDSON R R, BOOS D D, KAPLAN N L. A statistical test for detecting geographic subdivision [J]. Molecular Biology and Evolution, 1992, 9(1): 138-151.

    [22]HUDSON R R, RICHARD R. A new statistic for detecting genetic differentiation [J]. Genetics, 2000, 155(4):2011-2014.

    [23]KOSAKOVSKY P S L, POSADA D, GRAVENOR M B, et al. GARD: a genetic algorithm for recombination detection [J]. Bioinformatics, 2006, 22(24): 3096-3098.

    [24]NIELSEN R. Molecular signatures of natural selection [J]. Annual Review of Genetics, 2005, 39(1): 197-218.

    [25]CHANG Fei, GAO Fangluan, SHEN Jianguo, et al. Complete genome analysis of a PVYN-Wi recombinant isolate from Solanum tuberosum in China [J]. Potato Research, 2015, 58(4): 377-389.

    [26]祝雯, 詹家綏. 植物病原物的群體遺傳學[J]. 遺傳, 2012, 34(2): 157-166.

    [27]史鳳陽. 中國4個省份馬鈴薯Y病毒(PVY)P1基因的群體遺傳結構[D]. 福州: 福建農林大學, 2014.

    [28]常飛. 馬鈴薯Y病毒pipo基因的遺傳結構分析[D]. 福州: 福建農林大學, 2015.

    [29]謝聯(lián)輝, 林奇英. 植物病毒學[M]. 北京: 中國農業(yè)出版社, 2011.

    [30]LORENZEN J H, PICHE L M, GUDMESTAD N C, et al. A multiplex PCR assay to characterize Potato virus Y isolates and identify strain mixtures [J]. Plant Disease, 2006, 90(7): 935-940.

    [31]VISSER J C, BELLSTEDT D U, PIRIE M D. The recent recombinant evolution of a major crop pathogen, Potato virus Y [J/OL]. PLoS ONE, 2012, 7(11): 50631. DOI: 10.1371/journal. pone. 0050631.

    [32]QUENOUILLE J, VASSILAKOS N, MOURY B. Potato virus Y: a major crop pathogen that has provided major insights into the evolution of viral pathogenicity [J]. Molecular Plant Pathology, 2013, 14(5): 439-452.

    [33]CHAO Lin. Fitness of RNA virus decreased by Mullers ratchet [J]. Nature, 1990, 348(6300): 454-455.

    [34]STEINHAUER D A, HOLLAND J J. Rapid evolution of RNA viruses [J]. Annual Review of Microbiology, 1987, 41(1): 409-433.

    [35]ADAMS M J, ANTONIW J F, FAUQUET C M. Molecular criteria for genus and species discrimination within the family Potyviridae [J]. Archives of Virology, 2005, 150(3): 459-479.

    (責任編輯:田 喆)

    猜你喜歡
    遺傳多樣性
    從葉綠體DNA角度分析云南省砂梨地方品種遺傳多樣性
    寧夏外引水稻種質資源表型性狀遺傳多樣性分析
    寧夏外引水稻種質資源表型性狀遺傳多樣性分析
    茄子種質資源農藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    金魚起源及遺傳多樣性研究進展
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質資源遺傳多樣性研究進展
    金銀花SSR指紋圖譜的構建及遺傳多樣性分析
    赤兔流量卡办理| 特大巨黑吊av在线直播| 亚洲欧洲日产国产| 亚洲欧美精品专区久久| 精品国产露脸久久av麻豆| 在线观看人妻少妇| 国产精品久久久久久精品古装| 国产精品熟女久久久久浪| 免费观看无遮挡的男女| 久久久久久久久大av| 亚洲精品aⅴ在线观看| 欧美人与善性xxx| 国产爽快片一区二区三区| 老司机影院成人| 国产黄片视频在线免费观看| 亚洲av成人精品一二三区| 免费av不卡在线播放| av在线app专区| 亚洲欧美一区二区三区国产| 免费观看性生交大片5| 高清日韩中文字幕在线| 久久久久久九九精品二区国产| 国产探花极品一区二区| 成年女人在线观看亚洲视频| 女的被弄到高潮叫床怎么办| 色婷婷久久久亚洲欧美| 成人漫画全彩无遮挡| 十分钟在线观看高清视频www | 亚洲欧美一区二区三区国产| 国产综合精华液| 亚洲欧洲国产日韩| 男的添女的下面高潮视频| 免费观看在线日韩| 日韩成人伦理影院| a 毛片基地| 久久久久久久精品精品| 国产精品秋霞免费鲁丝片| 在线看a的网站| 国产成人午夜福利电影在线观看| 又爽又黄a免费视频| 王馨瑶露胸无遮挡在线观看| 伦理电影免费视频| 久久青草综合色| 大片免费播放器 马上看| 春色校园在线视频观看| 五月天丁香电影| 国产中年淑女户外野战色| 久久久久久久大尺度免费视频| 99热6这里只有精品| av女优亚洲男人天堂| 插阴视频在线观看视频| 日本av手机在线免费观看| 精品一区在线观看国产| 日韩三级伦理在线观看| 深夜a级毛片| 男女啪啪激烈高潮av片| 久久鲁丝午夜福利片| 又爽又黄a免费视频| 中文字幕精品免费在线观看视频 | 国产精品成人在线| 亚洲精品视频女| 男的添女的下面高潮视频| 三级经典国产精品| 亚洲成色77777| 欧美最新免费一区二区三区| 亚洲最大成人中文| 国产av一区二区精品久久 | av视频免费观看在线观看| 91精品国产九色| 亚洲av免费高清在线观看| 欧美成人a在线观看| 狂野欧美激情性xxxx在线观看| 亚洲自偷自拍三级| 久久精品久久久久久噜噜老黄| 亚洲精品日韩av片在线观看| 精品人妻偷拍中文字幕| 欧美精品一区二区大全| 寂寞人妻少妇视频99o| 国产伦在线观看视频一区| av免费观看日本| 蜜臀久久99精品久久宅男| 男女边摸边吃奶| 国产精品国产三级专区第一集| 观看av在线不卡| 亚洲天堂av无毛| 日本av免费视频播放| 亚洲电影在线观看av| 亚洲最大成人中文| 亚洲无线观看免费| h视频一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲av免费高清在线观看| 六月丁香七月| 国产高清有码在线观看视频| 男女边吃奶边做爰视频| 99久久中文字幕三级久久日本| 亚洲成人中文字幕在线播放| 日韩av不卡免费在线播放| 99视频精品全部免费 在线| 十八禁网站网址无遮挡 | 如何舔出高潮| 最近2019中文字幕mv第一页| 亚洲欧美清纯卡通| 偷拍熟女少妇极品色| 一级毛片电影观看| 五月玫瑰六月丁香| 亚洲国产精品专区欧美| 看非洲黑人一级黄片| 国产乱人视频| 国产女主播在线喷水免费视频网站| 亚洲国产高清在线一区二区三| 亚洲成色77777| 少妇高潮的动态图| 久久精品国产亚洲av涩爱| 久久久a久久爽久久v久久| 乱码一卡2卡4卡精品| 国产精品人妻久久久久久| 久久鲁丝午夜福利片| 高清视频免费观看一区二区| 一级毛片aaaaaa免费看小| 精品久久久噜噜| 精品久久久噜噜| 精品久久久噜噜| 国产午夜精品一二区理论片| 成人午夜精彩视频在线观看| 亚洲av成人精品一区久久| 老女人水多毛片| 国产成人a区在线观看| 黄色欧美视频在线观看| 一个人免费看片子| 亚洲av成人精品一区久久| 欧美一级a爱片免费观看看| 国产高清三级在线| 欧美一级a爱片免费观看看| 国产色婷婷99| 欧美97在线视频| 另类亚洲欧美激情| 97超碰精品成人国产| 小蜜桃在线观看免费完整版高清| 一个人免费看片子| 尤物成人国产欧美一区二区三区| 国产在线一区二区三区精| 纯流量卡能插随身wifi吗| 欧美成人精品欧美一级黄| 国产精品一二三区在线看| 免费黄色在线免费观看| 国产片特级美女逼逼视频| 热99国产精品久久久久久7| 建设人人有责人人尽责人人享有的 | 在线观看三级黄色| 美女高潮的动态| 久久久久性生活片| 国产欧美日韩一区二区三区在线 | 亚洲精品久久久久久婷婷小说| freevideosex欧美| 日产精品乱码卡一卡2卡三| 亚洲精品亚洲一区二区| 中文字幕av成人在线电影| 久久久久久九九精品二区国产| 日韩一本色道免费dvd| 日韩av不卡免费在线播放| 天天躁夜夜躁狠狠久久av| 精品人妻偷拍中文字幕| 亚洲av国产av综合av卡| 黄片无遮挡物在线观看| 久久国产乱子免费精品| 最黄视频免费看| www.色视频.com| 国产一区亚洲一区在线观看| 精品午夜福利在线看| 一级a做视频免费观看| 免费av中文字幕在线| 成人特级av手机在线观看| 亚洲天堂av无毛| 国产极品天堂在线| 人体艺术视频欧美日本| 亚洲欧美一区二区三区黑人 | 在线免费十八禁| 51国产日韩欧美| 建设人人有责人人尽责人人享有的 | 男人狂女人下面高潮的视频| 精品99又大又爽又粗少妇毛片| 国产日韩欧美亚洲二区| 黄色欧美视频在线观看| 亚洲av中文字字幕乱码综合| 久久久欧美国产精品| 尾随美女入室| 亚洲国产日韩一区二区| 波野结衣二区三区在线| 亚洲第一区二区三区不卡| 男女啪啪激烈高潮av片| 最黄视频免费看| 日韩三级伦理在线观看| 亚洲婷婷狠狠爱综合网| 亚洲四区av| 51国产日韩欧美| 国产精品久久久久久精品电影小说 | 又爽又黄a免费视频| 久热这里只有精品99| 久久久久久久久久成人| 精品国产一区二区三区久久久樱花 | 久久久久久久久久久免费av| 高清午夜精品一区二区三区| 国产永久视频网站| 欧美bdsm另类| 色吧在线观看| 国产精品一区二区在线不卡| 亚洲成人手机| 亚洲av综合色区一区| av线在线观看网站| 成人免费观看视频高清| 亚洲精品亚洲一区二区| 精华霜和精华液先用哪个| 十分钟在线观看高清视频www | 最近最新中文字幕大全电影3| 丝袜喷水一区| 不卡视频在线观看欧美| 亚洲欧美精品专区久久| 精品久久久久久久久av| 日韩强制内射视频| 久久6这里有精品| 国产在视频线精品| 国产 精品1| 久热久热在线精品观看| 亚洲不卡免费看| 亚洲国产毛片av蜜桃av| 欧美变态另类bdsm刘玥| 日本黄色日本黄色录像| 在线天堂最新版资源| 男女边摸边吃奶| 亚洲天堂av无毛| 97超视频在线观看视频| 亚洲精品国产色婷婷电影| 最近中文字幕高清免费大全6| 国模一区二区三区四区视频| 国产片特级美女逼逼视频| 亚洲成人一二三区av| 老熟女久久久| av专区在线播放| 在线播放无遮挡| freevideosex欧美| 日本色播在线视频| 最近最新中文字幕免费大全7| 麻豆国产97在线/欧美| 色视频在线一区二区三区| 欧美日韩国产mv在线观看视频 | 在线观看免费视频网站a站| 久久精品国产鲁丝片午夜精品| 欧美高清成人免费视频www| 国产在线免费精品| 亚洲成人av在线免费| 99精国产麻豆久久婷婷| 大香蕉久久网| 国产高清不卡午夜福利| 老司机影院毛片| 国产在线免费精品| 国产高清不卡午夜福利| 成年美女黄网站色视频大全免费 | 欧美精品国产亚洲| 日韩三级伦理在线观看| av又黄又爽大尺度在线免费看| 男女边吃奶边做爰视频| 国产亚洲91精品色在线| 自拍欧美九色日韩亚洲蝌蚪91 | 少妇人妻一区二区三区视频| 国产一区二区三区av在线| 亚洲精品乱久久久久久| 91在线精品国自产拍蜜月| 免费看日本二区| 美女中出高潮动态图| 国产探花极品一区二区| 黑丝袜美女国产一区| 成人毛片a级毛片在线播放| 国产亚洲5aaaaa淫片| 99久久精品热视频| 男女边吃奶边做爰视频| 欧美成人一区二区免费高清观看| 晚上一个人看的免费电影| 国产成人精品福利久久| 亚洲综合色惰| 99久久综合免费| 直男gayav资源| 五月开心婷婷网| 26uuu在线亚洲综合色| 久久久久久久精品精品| 日韩一区二区三区影片| 国产日韩欧美在线精品| 国产免费福利视频在线观看| 有码 亚洲区| 久久久色成人| 中文在线观看免费www的网站| 国产探花极品一区二区| 老司机影院毛片| 欧美日韩精品成人综合77777| 成人影院久久| 男人爽女人下面视频在线观看| 尤物成人国产欧美一区二区三区| 亚洲av二区三区四区| 日本午夜av视频| 麻豆国产97在线/欧美| 久久青草综合色| 一级黄片播放器| 日韩电影二区| 欧美日韩综合久久久久久| 欧美精品一区二区免费开放| 国产视频首页在线观看| 一级黄片播放器| 男女无遮挡免费网站观看| 永久网站在线| 18禁裸乳无遮挡动漫免费视频| 久久韩国三级中文字幕| 18禁动态无遮挡网站| 日本黄色日本黄色录像| 在线观看一区二区三区| 欧美成人a在线观看| 国产乱人偷精品视频| 国产av精品麻豆| 少妇人妻久久综合中文| 国精品久久久久久国模美| 美女高潮的动态| 91狼人影院| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 精品亚洲成a人片在线观看 | 亚洲欧美成人综合另类久久久| 欧美一区二区亚洲| 欧美xxxx性猛交bbbb| 成人综合一区亚洲| 国产熟女欧美一区二区| 欧美bdsm另类| 啦啦啦中文免费视频观看日本| 人妻系列 视频| 秋霞伦理黄片| av网站免费在线观看视频| 亚洲av福利一区| 美女cb高潮喷水在线观看| 丝瓜视频免费看黄片| 黄色配什么色好看| 亚洲欧美日韩无卡精品| 精品一区在线观看国产| 国产精品99久久99久久久不卡 | 精品亚洲成a人片在线观看 | 偷拍熟女少妇极品色| 秋霞在线观看毛片| 国产伦精品一区二区三区四那| 大话2 男鬼变身卡| 爱豆传媒免费全集在线观看| 最后的刺客免费高清国语| 久久久久久久久久久丰满| 麻豆成人av视频| 国产黄频视频在线观看| av一本久久久久| 韩国av在线不卡| 中国美白少妇内射xxxbb| 精品人妻偷拍中文字幕| 少妇的逼好多水| 人妻少妇偷人精品九色| 成人一区二区视频在线观看| 91久久精品电影网| 亚洲精品国产色婷婷电影| 一级黄片播放器| 99九九线精品视频在线观看视频| 简卡轻食公司| av专区在线播放| 在线观看av片永久免费下载| 中文精品一卡2卡3卡4更新| 最近手机中文字幕大全| 国产日韩欧美亚洲二区| 熟女电影av网| 97超视频在线观看视频| 中文字幕精品免费在线观看视频 | 中文字幕制服av| 黄色日韩在线| 七月丁香在线播放| 国产亚洲精品久久久com| 美女中出高潮动态图| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级 | 亚洲四区av| 男女无遮挡免费网站观看| 国产永久视频网站| 搡老乐熟女国产| 妹子高潮喷水视频| 中文资源天堂在线| 成人黄色视频免费在线看| 久久韩国三级中文字幕| 亚洲国产精品国产精品| 全区人妻精品视频| 2018国产大陆天天弄谢| 51国产日韩欧美| 国产伦理片在线播放av一区| 亚洲天堂av无毛| 国产高清有码在线观看视频| 亚洲三级黄色毛片| 日本与韩国留学比较| 天美传媒精品一区二区| 18禁动态无遮挡网站| 三级经典国产精品| 欧美国产精品一级二级三级 | 天堂俺去俺来也www色官网| 岛国毛片在线播放| 亚洲精品乱久久久久久| 简卡轻食公司| 十八禁网站网址无遮挡 | 在线观看一区二区三区| 能在线免费看毛片的网站| 赤兔流量卡办理| 亚洲天堂av无毛| 中文字幕精品免费在线观看视频 | 欧美97在线视频| 亚洲精品久久久久久婷婷小说| 一区二区三区四区激情视频| 亚洲精品一区蜜桃| 国产高潮美女av| 狂野欧美激情性xxxx在线观看| 亚洲av不卡在线观看| 久久人妻熟女aⅴ| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久久免| 久久久久久久久久人人人人人人| 啦啦啦在线观看免费高清www| 国产精品国产三级国产专区5o| 欧美激情极品国产一区二区三区 | 男女边摸边吃奶| 99九九线精品视频在线观看视频| 丝瓜视频免费看黄片| 在线看a的网站| 久久久久久久精品精品| 亚洲av在线观看美女高潮| 欧美精品人与动牲交sv欧美| 久久综合国产亚洲精品| 老师上课跳d突然被开到最大视频| 精品人妻偷拍中文字幕| 欧美日韩视频高清一区二区三区二| 久久久久久久大尺度免费视频| 最近中文字幕2019免费版| 九色成人免费人妻av| 美女视频免费永久观看网站| 一个人看的www免费观看视频| 国产亚洲欧美精品永久| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 中文乱码字字幕精品一区二区三区| 国产伦理片在线播放av一区| 日本vs欧美在线观看视频 | 久久6这里有精品| videos熟女内射| 国内少妇人妻偷人精品xxx网站| 色视频在线一区二区三区| 亚洲不卡免费看| 欧美精品人与动牲交sv欧美| av在线播放精品| 欧美97在线视频| 国产一区二区在线观看日韩| 久久久久网色| 不卡视频在线观看欧美| av又黄又爽大尺度在线免费看| 干丝袜人妻中文字幕| 黄色一级大片看看| 国产一区二区在线观看日韩| 麻豆成人午夜福利视频| 伊人久久精品亚洲午夜| 一级片'在线观看视频| 久热久热在线精品观看| 日本黄大片高清| 午夜免费鲁丝| 大又大粗又爽又黄少妇毛片口| 观看av在线不卡| 日本黄色片子视频| 成人毛片a级毛片在线播放| 国精品久久久久久国模美| 人妻 亚洲 视频| 亚洲av成人精品一二三区| a 毛片基地| 亚洲性久久影院| 日本wwww免费看| 男人爽女人下面视频在线观看| 国产精品人妻久久久久久| 视频区图区小说| 简卡轻食公司| 女性生殖器流出的白浆| 一级毛片电影观看| 亚洲av在线观看美女高潮| 免费少妇av软件| 午夜福利高清视频| 嫩草影院入口| 欧美激情国产日韩精品一区| 黑人高潮一二区| 亚洲精品国产av蜜桃| 午夜福利高清视频| 18禁裸乳无遮挡免费网站照片| 色5月婷婷丁香| 九色成人免费人妻av| 欧美高清性xxxxhd video| 久久热精品热| 18禁裸乳无遮挡免费网站照片| 色5月婷婷丁香| 少妇的逼好多水| 丝袜喷水一区| 一区二区av电影网| 中文字幕亚洲精品专区| 久久99蜜桃精品久久| 丝袜脚勾引网站| 夜夜看夜夜爽夜夜摸| 久久6这里有精品| 内射极品少妇av片p| 夫妻午夜视频| 国产日韩欧美亚洲二区| 成人黄色视频免费在线看| 狂野欧美白嫩少妇大欣赏| 欧美日韩一区二区视频在线观看视频在线| 日本av手机在线免费观看| 成人一区二区视频在线观看| 99久久精品热视频| 亚洲成人一二三区av| 中文字幕av成人在线电影| 国语对白做爰xxxⅹ性视频网站| 欧美+日韩+精品| 午夜激情久久久久久久| 国产在线一区二区三区精| 国产午夜精品久久久久久一区二区三区| 老司机影院成人| 亚洲成色77777| 久久99热这里只频精品6学生| 午夜免费男女啪啪视频观看| 国产欧美另类精品又又久久亚洲欧美| 欧美极品一区二区三区四区| 国产精品国产三级专区第一集| 在现免费观看毛片| 欧美最新免费一区二区三区| 天美传媒精品一区二区| 国产av一区二区精品久久 | 大码成人一级视频| 纯流量卡能插随身wifi吗| 欧美亚洲 丝袜 人妻 在线| 日韩亚洲欧美综合| 久久久久久久久久人人人人人人| 日本av手机在线免费观看| 亚洲在久久综合| 久久99热这里只频精品6学生| 国产伦理片在线播放av一区| 日本猛色少妇xxxxx猛交久久| 最近最新中文字幕大全电影3| 男女下面进入的视频免费午夜| 热re99久久精品国产66热6| 欧美精品人与动牲交sv欧美| 最近手机中文字幕大全| 女人久久www免费人成看片| 成人漫画全彩无遮挡| 欧美变态另类bdsm刘玥| 国产亚洲午夜精品一区二区久久| av黄色大香蕉| 少妇人妻 视频| av又黄又爽大尺度在线免费看| 久久精品国产亚洲av天美| 国产真实伦视频高清在线观看| 97精品久久久久久久久久精品| 边亲边吃奶的免费视频| 我要看黄色一级片免费的| 国产男女超爽视频在线观看| 欧美zozozo另类| 伊人久久精品亚洲午夜| 舔av片在线| 国产成人午夜福利电影在线观看| 麻豆国产97在线/欧美| 中文字幕久久专区| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产亚洲av天美| 3wmmmm亚洲av在线观看| 亚洲av综合色区一区| av播播在线观看一区| 在线免费十八禁| 我的老师免费观看完整版| 男人和女人高潮做爰伦理| 高清午夜精品一区二区三区| 久久午夜福利片| h日本视频在线播放| 久久久久久久精品精品| 99re6热这里在线精品视频| 一级毛片久久久久久久久女| 在线播放无遮挡| 久久久久人妻精品一区果冻| 国产成人午夜福利电影在线观看| 97在线视频观看| 小蜜桃在线观看免费完整版高清| 色哟哟·www| 午夜视频国产福利| 丰满乱子伦码专区| 另类亚洲欧美激情| 激情 狠狠 欧美| 在线播放无遮挡| 成人18禁高潮啪啪吃奶动态图 | 亚洲电影在线观看av| 国产淫片久久久久久久久| 日日啪夜夜撸| 18禁动态无遮挡网站| 亚洲伊人久久精品综合| videos熟女内射| 爱豆传媒免费全集在线观看| 又黄又爽又刺激的免费视频.| 免费不卡的大黄色大毛片视频在线观看| 中文在线观看免费www的网站| 99热网站在线观看| 久久久a久久爽久久v久久| 日韩大片免费观看网站| 最近最新中文字幕大全电影3| 久久久久视频综合| 免费人成在线观看视频色| 精品久久久久久久久av| 80岁老熟妇乱子伦牲交| 精品99又大又爽又粗少妇毛片| 欧美日韩精品成人综合77777| 成年av动漫网址| 亚洲国产精品专区欧美| 免费黄频网站在线观看国产|