• <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指紋圖譜的構建及遺傳多樣性分析
    日韩欧美 国产精品| 亚洲真实伦在线观看| 另类亚洲欧美激情| 极品教师在线视频| 最近中文字幕高清免费大全6| 中文字幕av电影在线播放| 亚洲真实伦在线观看| 精品国产露脸久久av麻豆| 久久这里有精品视频免费| 国产一级毛片在线| 亚洲av电影在线观看一区二区三区| 五月伊人婷婷丁香| 欧美日本中文国产一区发布| 又爽又黄a免费视频| 简卡轻食公司| 午夜福利在线观看免费完整高清在| 新久久久久国产一级毛片| 高清不卡的av网站| 日日撸夜夜添| 亚洲国产精品一区二区三区在线| 亚洲精品一二三| 国产深夜福利视频在线观看| a 毛片基地| av天堂久久9| 免费人妻精品一区二区三区视频| 九九爱精品视频在线观看| 五月天丁香电影| 国产精品欧美亚洲77777| 久久女婷五月综合色啪小说| 免费少妇av软件| 丝袜在线中文字幕| 亚洲欧美中文字幕日韩二区| 噜噜噜噜噜久久久久久91| 国内少妇人妻偷人精品xxx网站| 日韩,欧美,国产一区二区三区| 青春草视频在线免费观看| 在线亚洲精品国产二区图片欧美 | 伊人久久国产一区二区| 精品卡一卡二卡四卡免费| 最近最新中文字幕免费大全7| 亚洲一区二区三区欧美精品| 大陆偷拍与自拍| 精品视频人人做人人爽| 国产精品国产三级国产专区5o| 国产成人a∨麻豆精品| 成人特级av手机在线观看| 少妇丰满av| 日韩一本色道免费dvd| 国产白丝娇喘喷水9色精品| 亚洲,欧美,日韩| 乱码一卡2卡4卡精品| 亚洲欧美日韩另类电影网站| 欧美日韩av久久| 一边亲一边摸免费视频| 26uuu在线亚洲综合色| 3wmmmm亚洲av在线观看| 国产黄色免费在线视频| freevideosex欧美| 久久久久久久久久久久大奶| 亚洲精品一二三| 国产精品不卡视频一区二区| 看免费成人av毛片| 午夜福利,免费看| 日韩视频在线欧美| 一级av片app| 建设人人有责人人尽责人人享有的| 老司机影院成人| 亚洲欧美清纯卡通| 午夜91福利影院| 色吧在线观看| 免费播放大片免费观看视频在线观看| 亚洲欧美精品专区久久| av天堂中文字幕网| av视频免费观看在线观看| 哪个播放器可以免费观看大片| 亚洲电影在线观看av| 99精国产麻豆久久婷婷| 国产熟女欧美一区二区| 精品少妇黑人巨大在线播放| 国产午夜精品一二区理论片| 亚洲av中文av极速乱| 你懂的网址亚洲精品在线观看| 三级经典国产精品| 日韩av不卡免费在线播放| 一区二区三区精品91| 少妇人妻 视频| 国产精品偷伦视频观看了| 国产精品一区二区性色av| 久久人妻熟女aⅴ| 美女大奶头黄色视频| 久久女婷五月综合色啪小说| 精品国产一区二区三区久久久樱花| 成人亚洲精品一区在线观看| 熟女电影av网| 亚洲图色成人| 女的被弄到高潮叫床怎么办| 久热久热在线精品观看| 热re99久久精品国产66热6| 欧美日韩亚洲高清精品| 青春草国产在线视频| 三级国产精品片| 精品人妻熟女av久视频| 亚洲人与动物交配视频| 亚洲av免费高清在线观看| 一级二级三级毛片免费看| 精品酒店卫生间| 国产欧美亚洲国产| 亚洲成人av在线免费| 一区二区三区四区激情视频| 国产又色又爽无遮挡免| 午夜影院在线不卡| 日日啪夜夜爽| 国产深夜福利视频在线观看| 日韩一区二区三区影片| 欧美 亚洲 国产 日韩一| 春色校园在线视频观看| 人人澡人人妻人| 免费av中文字幕在线| 噜噜噜噜噜久久久久久91| 亚洲情色 制服丝袜| 高清av免费在线| 黑人猛操日本美女一级片| 色婷婷av一区二区三区视频| 亚洲怡红院男人天堂| 视频区图区小说| 极品人妻少妇av视频| 有码 亚洲区| 人人妻人人添人人爽欧美一区卜| 精品少妇久久久久久888优播| 视频中文字幕在线观看| 美女xxoo啪啪120秒动态图| 欧美日韩视频高清一区二区三区二| 99久久中文字幕三级久久日本| 日韩大片免费观看网站| 女人精品久久久久毛片| 色婷婷av一区二区三区视频| 国产亚洲91精品色在线| 国产高清国产精品国产三级| 亚洲精品国产色婷婷电影| 久久午夜综合久久蜜桃| 亚洲自偷自拍三级| 搡老乐熟女国产| 亚洲自偷自拍三级| 成人国产麻豆网| 亚洲国产精品专区欧美| 99久久中文字幕三级久久日本| 五月玫瑰六月丁香| 国产男女内射视频| 日韩大片免费观看网站| 少妇人妻久久综合中文| 日本猛色少妇xxxxx猛交久久| 黄色毛片三级朝国网站 | 日韩视频在线欧美| av在线app专区| 中国美白少妇内射xxxbb| 寂寞人妻少妇视频99o| 在线观看免费日韩欧美大片 | 最近最新中文字幕免费大全7| 国产精品一区二区在线不卡| 免费观看无遮挡的男女| 丝袜脚勾引网站| 特大巨黑吊av在线直播| 秋霞伦理黄片| 黄色日韩在线| 欧美日韩在线观看h| 美女国产视频在线观看| 久久青草综合色| 水蜜桃什么品种好| 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 日韩av在线免费看完整版不卡| 亚洲一级一片aⅴ在线观看| 久久97久久精品| 国产一级毛片在线| 国产av码专区亚洲av| 天天操日日干夜夜撸| 精品国产乱码久久久久久小说| 一级片'在线观看视频| 9色porny在线观看| 亚洲综合精品二区| 中文字幕免费在线视频6| 亚洲欧美成人精品一区二区| 久久人妻熟女aⅴ| 日日撸夜夜添| 亚洲精品久久午夜乱码| 国产精品久久久久久精品电影小说| 九九爱精品视频在线观看| 97超碰精品成人国产| 久久人妻熟女aⅴ| 一区在线观看完整版| 女人精品久久久久毛片| 91精品伊人久久大香线蕉| 极品教师在线视频| 男女无遮挡免费网站观看| 夫妻性生交免费视频一级片| 国产探花极品一区二区| 亚洲性久久影院| 纵有疾风起免费观看全集完整版| 久久精品国产自在天天线| 九色成人免费人妻av| 中国国产av一级| 最后的刺客免费高清国语| 免费av中文字幕在线| a级毛片免费高清观看在线播放| 国精品久久久久久国模美| 一本一本综合久久| 黄色视频在线播放观看不卡| 亚洲人成网站在线播| 久久狼人影院| 国产伦在线观看视频一区| 永久免费av网站大全| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 国产综合精华液| 国产视频内射| 久久国产精品男人的天堂亚洲 | 最近手机中文字幕大全| 一级av片app| 久久久久国产网址| 中文欧美无线码| 亚洲国产欧美在线一区| 国模一区二区三区四区视频| 亚洲美女搞黄在线观看| 欧美日本中文国产一区发布| 欧美精品人与动牲交sv欧美| 男女啪啪激烈高潮av片| 美女国产视频在线观看| 免费av中文字幕在线| 人人妻人人添人人爽欧美一区卜| 国产精品国产三级国产av玫瑰| 精品一区二区免费观看| 22中文网久久字幕| 精品久久久久久电影网| 日本爱情动作片www.在线观看| 午夜福利影视在线免费观看| 狂野欧美白嫩少妇大欣赏| 久久精品国产a三级三级三级| 亚洲无线观看免费| 久久国产精品大桥未久av | 亚洲内射少妇av| 国产精品女同一区二区软件| 国产免费一区二区三区四区乱码| 下体分泌物呈黄色| 乱人伦中国视频| 亚洲国产精品999| 中文字幕人妻熟人妻熟丝袜美| 成人无遮挡网站| 在线免费观看不下载黄p国产| 亚洲图色成人| 日日摸夜夜添夜夜添av毛片| 肉色欧美久久久久久久蜜桃| 日本av免费视频播放| 日韩中文字幕视频在线看片| 日韩成人av中文字幕在线观看| 日韩亚洲欧美综合| 美女视频免费永久观看网站| 国产成人精品一,二区| 欧美日本中文国产一区发布| 欧美日韩综合久久久久久| 波野结衣二区三区在线| 天天操日日干夜夜撸| 成人国产av品久久久| 免费大片18禁| 中文字幕精品免费在线观看视频 | 久久 成人 亚洲| 伦理电影免费视频| 亚洲,一卡二卡三卡| 免费黄色在线免费观看| 亚洲av综合色区一区| 乱码一卡2卡4卡精品| 精品酒店卫生间| 亚洲一区二区三区欧美精品| 久久人人爽人人爽人人片va| 亚洲精品乱久久久久久| 亚洲精品久久午夜乱码| 观看av在线不卡| 亚洲欧洲精品一区二区精品久久久 | 国产精品国产三级国产专区5o| 久久精品国产a三级三级三级| 国产毛片在线视频| 男女免费视频国产| 国产真实伦视频高清在线观看| 久久久久人妻精品一区果冻| 精品一区二区三区视频在线| 亚洲精品第二区| 51国产日韩欧美| 欧美xxxx性猛交bbbb| 国产淫语在线视频| 人体艺术视频欧美日本| 国产精品福利在线免费观看| 3wmmmm亚洲av在线观看| 久久ye,这里只有精品| 国产日韩欧美视频二区| 一级爰片在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲av欧美aⅴ国产| 街头女战士在线观看网站| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 亚洲成人一二三区av| 成年女人在线观看亚洲视频| av福利片在线观看| 男女国产视频网站| 精品久久久久久久久av| 久久国产亚洲av麻豆专区| 麻豆乱淫一区二区| 国产成人精品一,二区| 久久久久久久久久人人人人人人| 日韩精品有码人妻一区| 亚洲欧洲精品一区二区精品久久久 | 亚洲色图综合在线观看| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 国产亚洲5aaaaa淫片| 一区二区三区精品91| 成人免费观看视频高清| .国产精品久久| 一本—道久久a久久精品蜜桃钙片| 久久久久精品性色| kizo精华| 3wmmmm亚洲av在线观看| 亚洲va在线va天堂va国产| 日韩精品免费视频一区二区三区 | 亚洲欧美成人精品一区二区| 久久99精品国语久久久| 好男人视频免费观看在线| 少妇精品久久久久久久| 男女边摸边吃奶| 国产成人免费无遮挡视频| 老女人水多毛片| 日韩一本色道免费dvd| 嫩草影院新地址| 涩涩av久久男人的天堂| 天天躁夜夜躁狠狠久久av| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| av网站免费在线观看视频| 日韩强制内射视频| 国产精品久久久久久精品古装| 久久久久久久久久久免费av| 街头女战士在线观看网站| 色视频在线一区二区三区| 免费看av在线观看网站| 中国美白少妇内射xxxbb| 国产真实伦视频高清在线观看| 午夜91福利影院| 在线观看免费视频网站a站| 3wmmmm亚洲av在线观看| 99久久精品国产国产毛片| 午夜免费男女啪啪视频观看| 久久精品国产自在天天线| 亚洲人与动物交配视频| 日本黄色日本黄色录像| 日韩大片免费观看网站| 久久精品国产亚洲网站| 久久国产精品大桥未久av | tube8黄色片| 高清午夜精品一区二区三区| 免费av中文字幕在线| 中文字幕免费在线视频6| 交换朋友夫妻互换小说| 嫩草影院新地址| 精品一品国产午夜福利视频| 99热国产这里只有精品6| 97精品久久久久久久久久精品| 伊人久久国产一区二区| 99热全是精品| 亚洲三级黄色毛片| 少妇人妻一区二区三区视频| 久久国内精品自在自线图片| 一个人看视频在线观看www免费| 亚洲一级一片aⅴ在线观看| a级毛色黄片| 亚洲成人av在线免费| 亚洲欧洲精品一区二区精品久久久 | 99久久精品一区二区三区| 成人亚洲精品一区在线观看| 爱豆传媒免费全集在线观看| 丰满饥渴人妻一区二区三| 你懂的网址亚洲精品在线观看| 97在线视频观看| 亚洲不卡免费看| 久久99一区二区三区| 欧美国产精品一级二级三级 | 亚洲美女黄色视频免费看| 能在线免费看毛片的网站| 国产午夜精品一二区理论片| 最新中文字幕久久久久| 人妻少妇偷人精品九色| 三级经典国产精品| 色哟哟·www| 成年av动漫网址| 日韩人妻高清精品专区| 黄色毛片三级朝国网站 | √禁漫天堂资源中文www| 亚洲国产成人一精品久久久| av黄色大香蕉| 久久女婷五月综合色啪小说| 少妇的逼好多水| 亚洲人成网站在线播| 亚洲精品视频女| 久久6这里有精品| 日韩精品免费视频一区二区三区 | 在现免费观看毛片| 少妇高潮的动态图| 国产亚洲91精品色在线| 美女视频免费永久观看网站| 日韩欧美 国产精品| 国产精品熟女久久久久浪| 国产成人午夜福利电影在线观看| 日韩制服骚丝袜av| 最近的中文字幕免费完整| 99视频精品全部免费 在线| 有码 亚洲区| 亚洲精品一区蜜桃| 男女啪啪激烈高潮av片| av卡一久久| 少妇的逼水好多| 一本—道久久a久久精品蜜桃钙片| 超碰97精品在线观看| 国产精品久久久久久久电影| 天美传媒精品一区二区| 亚洲国产精品国产精品| 亚洲一区二区三区欧美精品| 亚洲欧美清纯卡通| 91精品一卡2卡3卡4卡| 黑人高潮一二区| 三级国产精品片| 亚洲电影在线观看av| 91精品伊人久久大香线蕉| 欧美老熟妇乱子伦牲交| 久久久久久伊人网av| 极品教师在线视频| 精品人妻偷拍中文字幕| 美女福利国产在线| 免费人成在线观看视频色| 九九久久精品国产亚洲av麻豆| 在线免费观看不下载黄p国产| 国产精品久久久久久久电影| 亚洲精品久久午夜乱码| 国产淫语在线视频| 亚洲怡红院男人天堂| 热re99久久精品国产66热6| 22中文网久久字幕| 春色校园在线视频观看| 黄色配什么色好看| 国产视频首页在线观看| 国产成人精品一,二区| 日韩精品有码人妻一区| 天堂俺去俺来也www色官网| 亚洲精品aⅴ在线观看| 午夜91福利影院| 免费看光身美女| 一级a做视频免费观看| 久久精品国产亚洲av天美| 国产精品无大码| 最近2019中文字幕mv第一页| 国产 精品1| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 最近手机中文字幕大全| 99国产精品免费福利视频| 97超视频在线观看视频| 中文字幕免费在线视频6| freevideosex欧美| 色婷婷久久久亚洲欧美| 女性生殖器流出的白浆| 青春草亚洲视频在线观看| 精品一品国产午夜福利视频| 丰满人妻一区二区三区视频av| 一级黄片播放器| 一级毛片黄色毛片免费观看视频| 啦啦啦视频在线资源免费观看| 国产在线视频一区二区| 久久久久人妻精品一区果冻| 免费少妇av软件| 国产伦精品一区二区三区视频9| 丝袜在线中文字幕| 在线观看av片永久免费下载| 一级毛片我不卡| 成人亚洲精品一区在线观看| 午夜免费鲁丝| 亚洲人成网站在线观看播放| 爱豆传媒免费全集在线观看| 黑人巨大精品欧美一区二区蜜桃 | 精品99又大又爽又粗少妇毛片| av在线观看视频网站免费| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在| 99九九在线精品视频 | 国产一区二区在线观看日韩| 久久人人爽人人片av| 日本免费在线观看一区| 久久6这里有精品| 国产毛片在线视频| 免费久久久久久久精品成人欧美视频 | 极品教师在线视频| 伊人亚洲综合成人网| 一区二区三区乱码不卡18| 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 久久久精品免费免费高清| 极品人妻少妇av视频| 欧美成人精品欧美一级黄| 寂寞人妻少妇视频99o| 日韩免费高清中文字幕av| 久久女婷五月综合色啪小说| av在线播放精品| 九草在线视频观看| 人人妻人人添人人爽欧美一区卜| 亚洲无线观看免费| 国产亚洲最大av| 国产精品人妻久久久久久| 国产男女内射视频| 精品一区二区三卡| 亚洲国产欧美在线一区| 中文字幕久久专区| 一级片'在线观看视频| 成人美女网站在线观看视频| 亚洲美女视频黄频| 成年av动漫网址| 搡女人真爽免费视频火全软件| 久久久久久久国产电影| 少妇丰满av| 欧美一级a爱片免费观看看| 又大又黄又爽视频免费| 亚洲国产成人一精品久久久| 欧美日韩在线观看h| 国产精品久久久久久精品古装| 91aial.com中文字幕在线观看| 国产精品久久久久久精品电影小说| 亚洲av日韩在线播放| 国产淫片久久久久久久久| 国产成人精品无人区| 成年人午夜在线观看视频| 永久网站在线| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 一级毛片我不卡| 全区人妻精品视频| 国产在视频线精品| 全区人妻精品视频| 日本欧美国产在线视频| 国产免费一级a男人的天堂| 亚洲精品日本国产第一区| 午夜福利网站1000一区二区三区| 亚洲欧美中文字幕日韩二区| av网站免费在线观看视频| 一本一本综合久久| 国产精品久久久久久久电影| 亚洲欧美清纯卡通| 久久人妻熟女aⅴ| 国产片特级美女逼逼视频| 中文天堂在线官网| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 久久久久久久久久成人| 偷拍熟女少妇极品色| 尾随美女入室| 日本av手机在线免费观看| 亚洲欧美日韩东京热| 精品少妇内射三级| 波野结衣二区三区在线| 美女内射精品一级片tv| 精品国产一区二区三区久久久樱花| 插阴视频在线观看视频| 制服丝袜香蕉在线| 成年人午夜在线观看视频| 欧美精品一区二区大全| 少妇 在线观看| 中文欧美无线码| 亚洲精品国产av蜜桃| 亚洲美女视频黄频| 99热这里只有是精品50| 久久人人爽人人爽人人片va| 午夜免费男女啪啪视频观看| 亚洲自偷自拍三级| 久热这里只有精品99| 久久影院123| 看免费成人av毛片| 午夜福利视频精品| 久久久久久久久久久免费av| 亚洲综合精品二区| av网站免费在线观看视频| a级毛色黄片| 亚洲精品日韩av片在线观看| 黑人巨大精品欧美一区二区蜜桃 | 精品亚洲乱码少妇综合久久| 日日撸夜夜添| 久久国产精品男人的天堂亚洲 | 亚洲av二区三区四区| 夜夜骑夜夜射夜夜干| 成年美女黄网站色视频大全免费 | 久久精品国产亚洲网站| 深夜a级毛片| 亚洲av在线观看美女高潮| 在线观看三级黄色| 日韩欧美一区视频在线观看 | 亚洲精品第二区| 街头女战士在线观看网站| 久久久久久人妻| av免费在线看不卡| 日本欧美视频一区| 91久久精品国产一区二区三区| 国产精品99久久久久久久久| 在线观看免费高清a一片| 成年美女黄网站色视频大全免费 | 熟女人妻精品中文字幕| 一区二区三区精品91| 自线自在国产av| 色吧在线观看| 免费人成在线观看视频色| 一级毛片久久久久久久久女| 三级经典国产精品| 日韩亚洲欧美综合|