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

    基于小RNA深度測(cè)序技術(shù)鑒定侵染我國(guó)部分省市草莓種苗的病毒種類

    2022-06-07 10:41:29王佳白瑩雪祝寧齊長(zhǎng)紅任俊達(dá)韓成貴尚巧霞
    植物保護(hù) 2022年3期

    王佳 白瑩雪 祝寧 齊長(zhǎng)紅 任俊達(dá) 韓成貴 尚巧霞

    摘要 隨著草莓保護(hù)地栽培面積的增加和無(wú)性繁殖種苗的繁殖與調(diào)運(yùn),草莓病毒病的發(fā)生與流行日益嚴(yán)重。為明確侵染我國(guó)部分省市草莓種苗的病毒種類,應(yīng)用小RNA深度測(cè)序技術(shù)進(jìn)行檢測(cè),并利用RT-PCR技術(shù)對(duì)結(jié)果進(jìn)行驗(yàn)證及序列分析。結(jié)果表明,從來(lái)自我國(guó)7省市的41株具有典型病毒病癥狀的草莓種苗樣品中檢測(cè)到草莓斑駁病毒strawberry mottle virus (SMoV)、草莓鑲脈病毒strawberry vein banding virus(SVBV)和草莓輕型黃邊病毒strawberry mild yellow edge virus (SMYEV)3種。SMoV、SVBV和SMYEV的檢出率分別為34.1%、24.4%和2.4%。選取不同產(chǎn)地草莓種苗上檢出的不同病毒進(jìn)行部分序列測(cè)定和分析,獲得了3個(gè)SMoV分離物(四川分離物schhy13、遼寧分離物lnhy23和河北分離物hbhy28)的部分RNA1 3′端非編碼區(qū)606 bp核苷酸序列,其一致性為98.12%~99.34%。測(cè)定并獲得了5個(gè)SVBV分離物(遼寧分離物lnhy15、lnhy17、lnhy24、河北分離物hbhy28和陜西分離物shaxhy35)的ORF6的部分基因及其5′端非編碼區(qū)的核苷酸序列,其中部分ORF6基因序列的一致性為97.98%~99.8%。僅1個(gè)遼寧樣品中檢測(cè)到SMYEV,獲得了遼寧分離物lnhy26,其CP基因序列與中國(guó)甘肅分離物sy02的一致性為99.02%。病毒侵染草莓植株產(chǎn)生的典型癥狀有葉片畸形、變色和植株矮化。存在2種病毒和3種病毒復(fù)合侵染的情況,占所有陽(yáng)性樣品的41.2%。遼寧和河北的‘紅顏種苗中檢出SMoV和SVBV或SMoV、SVBV和SMYEV的復(fù)合侵染情況;四川和陜西的‘紅顏種苗僅檢出SMoV或SVBV,而其他省市和品種均未檢出病毒。研究結(jié)果可以為草莓病毒的檢測(cè)以及草莓病毒病的科學(xué)防治提供依據(jù)。

    關(guān)鍵詞 草莓病毒病;小RNA深度測(cè)序技術(shù);RT-PCR;病毒檢測(cè)

    中圖分類號(hào): S436.639

    文獻(xiàn)標(biāo)識(shí)碼: A

    DOI: 10.16688/j.zwbh.2021079

    Abstract With the increase of protected cultivation area, propagation and transportation of strawberry asexual reproductive seedlings between different regions are increasingly frequent, and the occurrence and prevalence of strawberry viral disease are also becoming more and more serious. To identify the viruses infecting strawberry plants from different provinces of China, small RNA deep sequencing assay was used to detect the viruses in strawberry seedling samples, and then the detected viruses were further verified by RT-PCR and sequence analysis. The results showed that strawberry mottle virus (SMoV), strawberry vein banding virus (SVBV) and strawberry mild yellow edge virus (SMYEV) were detected in 41 strawberry seedlings with typical symptoms from seven provinces of China. The positive rates of SMoV, SVBV and SMYEV detection were 34.1%, 24.4% and 2.4%, respectively. Strawberry seedling samples infected by different viruses from different regions were selected, from which the partial viral nucleotide sequences of the detected viruses were sequenced and analyzed. The partial 3′ UTR nucleotide sequences of RNA1 of three samples infected by SMoV from Sichuan, Liaoning and Hebei were sequenced, which was 606 bp in length. The partial SMoV nucleotide sequences among Sichuan isolate schhy13, Liaoning isolate lnhy23 and Hebei isolate hbhy28 shared identity from 98.12% to 99.34%. The sequences of partial ORF6 gene and its 5′ UTR of SVBV Liaoning isolates lnhy15, lnhy17, lnhy24, Hebei isolate hbhy28 and Shaanxi isolate shaxhy35 were obtained. The partial ORF6 gene nucleotide sequences of five SVBV samples shared identity from 97.98% to 99.8%. Only one sample with SMYEV was detected from Liaoning, and the CP gene sequence of Liaoning isolate lnhy26 shared identity of 99.02% with that of Gansu isolate sy02. Symptoms of strawberry plants infected by virus included leaf deformity, leaf discoloration and dwarf. The detection results showed that combined infection by two or three viruses existed, and the detection rate accounted for 41.2% of virus-positive samples. ‘Benihoppe seedlings from Liaoning and Hebei were simultaneously infected by SMoV and SVBV or SMoV, SVBV and SMYEV. ‘Benihoppe seedlings from Sichuan and Shaanxi were only infected by SMoV or SVBV, but no virus was detected in strawberry seedlings sampled from other provinces or in other varieties. This study provided important information and techniques for virus detection and strawberry viral disease management.1617EC4E-91A2-4AC0-96D3-3A8AA5458DFA

    Key words strawberry viral disease;small RNA deep sequencing technology;RT-PCR;virus detection

    草莓因具有種植周期短、投資少、經(jīng)濟(jì)收益高等優(yōu)點(diǎn)成為重要的經(jīng)濟(jì)作物,在全球范圍內(nèi)廣泛種植[1]。我國(guó)的草莓產(chǎn)業(yè)經(jīng)過(guò)近40年的迅猛發(fā)展,生產(chǎn)面積和產(chǎn)量已躍居世界首位[2]。

    在我國(guó),草莓以冬春保護(hù)地栽培和依賴匍匐莖無(wú)性繁殖種苗為主,高溫高濕的生長(zhǎng)環(huán)境和無(wú)性繁殖的種植方式導(dǎo)致草莓斑駁病毒strawberry mottle virus (SMoV)、草莓鑲脈病毒strawberry vein banding virus (SVBV)、草莓輕型黃邊病毒strawberry mild yellow edge virus (SMYEV)和草莓皺縮病毒strawberry crinkle virus(SCV)引起的病毒病發(fā)生日益嚴(yán)重,同時(shí)伴有這些病毒的新株系和病毒種類產(chǎn)生,嚴(yán)重影響草莓的品質(zhì)和產(chǎn)量,對(duì)草莓產(chǎn)業(yè)構(gòu)成威脅[1,3-5]。2008年楊洪一等[6]從起源自我國(guó)甘肅的‘五葉草莓和自育的‘長(zhǎng)虹2號(hào)中分離到1個(gè)SMYEV的新株系類型。2014年Chen等[7]首次報(bào)道黃瓜花葉病毒cucumber mosaic virus (CMV)對(duì)我國(guó)的草莓產(chǎn)業(yè)構(gòu)成威脅,2017年Ding等[8]首次發(fā)現(xiàn)福建的草莓感染草莓白化病毒strawberry pallidosis-associated virus (SPaV)。

    小RNA深度測(cè)序技術(shù)相對(duì)于ELISA和RT-PCR等檢測(cè)方法,其結(jié)果更全面,也更易發(fā)現(xiàn)未知的病毒[9]。我國(guó)學(xué)者已利用該技術(shù)從福建和湖南等地的草莓上發(fā)現(xiàn)了一些新的病毒種類[8,10]。因此,筆者采用小RNA深度測(cè)序技術(shù)對(duì)來(lái)自北京、河北、遼寧、內(nèi)蒙古、陜西、四川和云南的草莓種苗進(jìn)行檢測(cè),并利用RT-PCR對(duì)檢測(cè)結(jié)果進(jìn)行驗(yàn)證,對(duì)所得序列進(jìn)行比對(duì)分析,以明確我國(guó)不同省市草莓種苗病毒病的發(fā)生情況和病毒種類,為進(jìn)一步繁育草莓脫毒種苗和病毒病的科學(xué)防控提供理論依據(jù)。

    1 材料與方法

    1.1 樣品采集

    2019年11月27日,從北京市昌平區(qū)南邵鎮(zhèn)草莓種植園采集41株具有明顯畸形、矮化、斑駁和皺縮等典型病毒病癥狀的草莓種苗葉片樣品,種苗的來(lái)源和品種包括產(chǎn)自北京的‘白雪公主(1株)和‘隋珠(4株),河北的‘光點(diǎn)(2株)、‘紅顏(4株)和‘圣誕紅(2株),遼寧的‘紅顏(10株),內(nèi)蒙古的‘紅顏(2株),陜西的‘紅顏(2株),四川的‘紅顏(10株)和云南的‘圣誕紅(4株)。

    1.2 主要試劑

    EASYspin植物RNA快速提取試劑盒、2×Taq PCR Mix購(gòu)自北京愛(ài)博森生物科技有限公司;dNTPs、5×M-MLV逆轉(zhuǎn)錄反應(yīng)緩沖液、M-MLV反轉(zhuǎn)錄酶、RNase抑制劑購(gòu)自TaKaRa公司;DEPC水購(gòu)自生工生物工程(上海)股份有限公司。

    1.3 方法

    1.3.1 具有典型病毒病癥狀的草莓種苗葉片樣品的小RNA深度測(cè)序

    文庫(kù)構(gòu)建、測(cè)序及原始數(shù)據(jù)整理:由北京優(yōu)吉科技有限公司完成。提取41株具有典型病毒病癥狀的草莓種苗葉片樣品的總RNA,并使用Agilent 2100和NanoDrop測(cè)定總RNA的濃度、完整性和純度。質(zhì)量合格的RNA樣品采用NEBNext Multiplex Small RNA Library Prep Set for Illumina構(gòu)建小RNA文庫(kù)。先后在3′端和5′端加上接頭,再反轉(zhuǎn)錄合成cDNA。然后用Qubit 2.0和實(shí)時(shí)熒光定量PCR對(duì)文庫(kù)的濃度進(jìn)行定量,并用Agilent 2100對(duì)文庫(kù)的插入片段長(zhǎng)度(insert size)進(jìn)行檢測(cè)。最后,使用Illumina Hiseq 2000測(cè)序平臺(tái)對(duì)檢測(cè)合格的小RNA文庫(kù)進(jìn)行測(cè)序。從獲得的序列中篩除低質(zhì)量、長(zhǎng)度小于18 nt或大于40 nt、帶polyA或polyN等條件的小RNA序列,獲得有效數(shù)據(jù)(clean reads)[11]。使用PFOR軟件和Velvet軟件對(duì)clean reads中的序列進(jìn)行拼接得到重疊群(contigs),與NCBI中的病毒數(shù)據(jù)庫(kù)進(jìn)行BLASTn和BLASTx比對(duì)分析并注釋,獲得可能存在的病毒序列信息。

    對(duì)獲得的病毒序列信息進(jìn)行整理分析,初步確定感染病毒的草莓葉片樣品,以及可能存在的病毒種類。

    1.3.2 RT-PCR驗(yàn)證小RNA深度測(cè)序結(jié)果

    1.3.2.1 草莓組織總RNA的提取

    參照陳柳等[12]的方法從初步確定感染病毒的草莓葉片樣品中提取總RNA。

    1.3.2.2 RT-PCR

    上述提取的總RNA進(jìn)行反轉(zhuǎn)錄合成cDNA。以1.5 μL的總RNA為模板,加入隨機(jī)引物 (9mer)(50 μmol/L)和Oligo(dT)18 Primer(50 μmol/L)各0.5 μL,dNTPs(10 mmol/L)2 μL,DEPC水將體系補(bǔ)至15.5 μL。65℃水浴5 min,冰浴5 min。向體系中加入5×M-MLV逆轉(zhuǎn)錄反應(yīng)緩沖液4 μL,M-MLV反轉(zhuǎn)錄酶(100 U/μL)和RNase抑制劑(40 U/μL)各0.5 μL。42℃ 1 h,70℃ 15 min。合成的cDNA于-20℃保存。

    根據(jù)小RNA深度測(cè)序確定的病毒種類,利用病毒的特異性引物(表1)對(duì)cDNA進(jìn)行PCR擴(kuò)增。擴(kuò)增體系: cDNA模板2 μL,正向和反向引物(10 mmol/L)各1 μL,2×Taq PCR mix 12.5 μL,DEPC水將體系補(bǔ)至25 μL。PCR擴(kuò)增程序: 94℃預(yù)變性 2 min;94℃變性 30 s,退火40 s(退火溫度見(jiàn)表1),72℃延伸 60 s,共35個(gè)循環(huán);最后72℃延伸5 min。取5 μL PCR產(chǎn)物用1%瓊脂糖凝膠電泳檢測(cè)。1617EC4E-91A2-4AC0-96D3-3A8AA5458DFA

    SMoV[13]、SVBV[14]和SMYEV[15]的特異性引物由生工生物工程(上海)股份有限公司合成,分別能擴(kuò)增出SMoV RNA1的3′端部分非編碼區(qū)核苷酸序列、SVBV ORF6部分基因及其5′ 端非編碼區(qū)核苷酸序列和SMYEV完整的外殼蛋白(coat protein, CP)基因及該基因兩端部分非編碼區(qū)核苷酸序列。

    1.3.3 序列測(cè)定和分析比較

    取目的條帶明亮的PCR產(chǎn)物送交北京博邁德基因技術(shù)有限公司進(jìn)行序列測(cè)定。測(cè)序結(jié)果在NCBI上與國(guó)內(nèi)外已報(bào)道的相關(guān)序列進(jìn)行同源性比對(duì),并利用MEGA 7.0的鄰接法構(gòu)建系統(tǒng)發(fā)育樹。

    1.4 病毒檢出率計(jì)算

    統(tǒng)計(jì)檢測(cè)樣品中各種病毒的陽(yáng)性樣品數(shù),并計(jì)算病毒檢出率。病毒檢出率=陽(yáng)性樣品數(shù)/總檢測(cè)樣品數(shù)×100%[16]。

    2 結(jié)果與分析

    2.1 小RNA深度測(cè)序檢測(cè)結(jié)果

    對(duì)41株具有典型病毒病癥狀的草莓種苗樣品進(jìn)行小RNA深度測(cè)序,經(jīng)質(zhì)量控制獲得142 945 074條clean reads,其中長(zhǎng)度為21、22 nt和24 nt的小RNA序列最多,分別占clean reads數(shù)量的50.35%、14.17%和12.13%。與NCBI中的病毒數(shù)據(jù)庫(kù)比對(duì)結(jié)果顯示,共有488 276條reads來(lái)自病毒,占clean reads數(shù)量的0.34%。

    將長(zhǎng)度為18~40 nt的序列拼接后獲得16 202個(gè)contigs,經(jīng)BLASTn比對(duì)分析表明,有261個(gè)contigs與SMoV、SVBV或SMYEV的核苷酸序列的一致性較高,未比對(duì)到其他植物病毒。其中51個(gè)contigs與SMoV核苷酸序列的一致性為92.59%~100%,184個(gè)contigs與SVBV核苷酸序列的一致性為96.67%~100%,26個(gè)contigs與SMYEV核苷酸序列的一致性為84.08%~100%。還有31個(gè)contigs利用BLASTn比對(duì)不到結(jié)果,使用BLASTx比對(duì)后顯示,31個(gè)contigs與SMoV、SVBV或SMYEV的氨基酸序列具有較高的一致性,也未比對(duì)到其他的植物病毒。其中4個(gè)contigs與SMoV氨基酸序列的一致性為95%~100%,1個(gè)contigs與SVBV氨基酸序列的一致性為100%,26個(gè)contigs與SMYEV氨基酸序列的一致性為90%~100%(表2)。

    2.2 RT-PCR檢測(cè)結(jié)果

    利用SMoV、SVBV和SMYEV的特異性引物對(duì)41株草莓種苗樣品進(jìn)行RT-PCR檢測(cè),并用DNAMAN和MEGA 7.0對(duì)獲得的3種病毒分離物進(jìn)一步分析,以驗(yàn)證小RNA深度測(cè)序結(jié)果的可靠性。結(jié)果表明,17份樣品中檢測(cè)到這3種病毒,總檢出率為41.5%,SMoV、SVBV和SMYEV的檢出率分別為34.1%、24.4%和2.4%(表3)。

    2.2.1 SMoV不同分離物的序列分析

    利用SMoV特異性引物SM6126F1/SM6732R1進(jìn)行RT-PCR檢測(cè),可以擴(kuò)增出約606 bp的目的片段,包含RNA1 3′端非編碼區(qū)部分核苷酸序列。根據(jù)草莓種苗產(chǎn)地不同,選取不同陽(yáng)性樣品進(jìn)行序列測(cè)定,獲得SMoV四川分離物schhy13、遼寧分離物lnhy23和河北分離物hbhy28的序列,GenBank登錄號(hào)分別為MW556623、MW556622和MW556621。利用DNAMAN進(jìn)行序列比對(duì)分析,3個(gè)SMoV分離物的核苷酸序列一致性為98.12%~99.34%。與GenBank上已報(bào)道的核苷酸序列進(jìn)行比對(duì)和系統(tǒng)進(jìn)化分析表明,四川分離物schhy13、遼寧分離物lnhy23 和河北分離物hbhy28與加拿大分離物NB926序列(GenBank登錄號(hào):KU200453.1)的一致性較高,分別為97.68%、99.14%和98.19%;四川分離物schhy13和河北分離物hbhy28親緣關(guān)系較近,與遼寧分離物lnhy23親緣關(guān)系較遠(yuǎn),形成兩個(gè)獨(dú)立的小分支(圖1)。

    2.2.2 SVBV不同分離物的序列分析

    利用SVBV特異性引物SV5508F/SV6606R進(jìn)行RT-PCR檢測(cè),可以擴(kuò)增出SVBV約1 098 bp的目的片段,包含ORF6部分基因及其5′端非編碼區(qū)的核苷酸序列。根據(jù)草莓種苗產(chǎn)地不同,選取不同地區(qū)草莓種苗陽(yáng)性樣品進(jìn)行序列擴(kuò)增,獲得SVBV遼寧分離物lnhy15、lnhy17、lnhy24、河北分離物hbhy28和陜西分離物shaxhy35,GenBank登錄號(hào)分別為MW556618、MW556617、MW556616、MW556619和MW556615。利用DNAMAN進(jìn)行序列比對(duì)分析,5個(gè)SVBV分離物的部分ORF6基因序列一致性為97.98%~99.8%。與GenBank上已報(bào)道的基因序列比對(duì)分析,遼寧分離物lnhy15和lnhy17與北京分離物BJ1(GenBank登錄號(hào): KT250632.1)的核苷酸序列一致性分別為99.01%和99.42%;遼寧分離物lnhy24和河北分離物hbhy28與北京分離物BJHXTX1(GenBank登錄號(hào): MF197916.1)的核苷酸序列一致性分別為99.36%和99.44%;陜西分離物shaxhy35與貴州分離物GZ(GenBank登錄號(hào): MH894295.1)的核苷酸序列一致性為99.52%。系統(tǒng)進(jìn)化樹分析結(jié)果如圖2所示,遼寧分離物lnhy15、lnhy17、lnhy24和河北分離物hbhy28親緣關(guān)系較近,與陜西分離物shaxhy35親緣關(guān)系較遠(yuǎn),均與中國(guó)不同省市的分離物有較近的親緣關(guān)系,而與美國(guó)和加拿大分離物親緣關(guān)系較遠(yuǎn)。

    2.2.3 SMYEV不同分離物的序列分析

    利用SMYEV特異性引物YT1/Y2進(jìn)行RT-PCR檢測(cè),可以擴(kuò)增出SMYEV約861 bp的目的片段,包含完整的CP基因及該基因兩端部分非編碼區(qū)核苷酸序列[17]。僅獲得1個(gè)遼寧分離物lnhy26的序列(GenBank登錄號(hào): MW556620),與GenBank的SMYEV CP基因進(jìn)行序列比對(duì)和系統(tǒng)進(jìn)化樹分析,該分離物與中國(guó)甘肅分離物sy02(GenBank登錄號(hào): EU107084.1)的核苷酸序列一致性為99.02%,且親緣關(guān)系較近(圖3)。1617EC4E-91A2-4AC0-96D3-3A8AA5458DFA

    2.3 草莓病毒病典型癥狀

    經(jīng)小RNA深度測(cè)序和RT-PCR檢測(cè)不同地區(qū)草莓種苗的病毒侵染情況,結(jié)合表4中的田間癥狀調(diào)查結(jié)果表明,病毒單獨(dú)侵染或2~3種病毒復(fù)合侵染時(shí),草莓植株可表現(xiàn)出單一癥狀,也可同時(shí)表現(xiàn)出多種癥狀。被同一種病毒單獨(dú)侵染的草莓植株表現(xiàn)出的癥狀也略有不同。主要癥狀如下:

    葉片畸形:葉片皺縮;三出復(fù)葉的1~2片小葉生長(zhǎng)較?。▓D4b-d)。

    葉片變色:葉片呈不同程度的斑駁、鑲脈或褪綠;葉脈變黑(圖4b-f)。

    植株矮化:株高較矮、葉片較小;葉柄細(xì);長(zhǎng)勢(shì)弱(圖4f)。

    2.4 不同產(chǎn)地草莓種苗感染病毒的情況

    對(duì)北京、河北、內(nèi)蒙古、四川、遼寧、云南和陜西共7個(gè)?。ㄊ?、自治區(qū))生產(chǎn)的草莓種苗病毒病進(jìn)行調(diào)查和檢測(cè),結(jié)果如表3所示。遼寧草莓種苗受病毒危害最為嚴(yán)重,病毒檢出率高達(dá)100.0%,陜西、四川和河北種苗的檢出率分別為50.0%、40.0%和25.0%,而北京、內(nèi)蒙古和云南的草莓種苗未檢測(cè)到病毒。

    不同省市的草莓種苗檢出病毒的種類也略有不同。調(diào)查結(jié)果顯示,四川、遼寧和河北的草莓種苗中檢出SMoV,遼寧、河北和陜西的種苗中檢出SVBV,僅遼寧的種苗中檢出SMYEV。病毒復(fù)合侵染發(fā)生較為嚴(yán)重,占總檢出率的41.2%。遼寧和河北的草莓種苗中檢出SMoV和SVBV兩種病毒復(fù)合侵染,僅在1株遼寧草莓種苗中檢出3種病毒復(fù)合侵染。

    對(duì)不同省市的‘紅顏‘白雪公主‘隋珠‘光點(diǎn)和‘圣誕紅5個(gè)草莓品種調(diào)查發(fā)現(xiàn),17個(gè)病毒陽(yáng)性樣品的品種均為‘紅顏,而‘白雪公主 ‘隋珠‘光點(diǎn)和‘圣誕紅均未檢出病毒。

    3 討論

    小RNA深度測(cè)序技術(shù)的出現(xiàn)提高了植物病毒病病原學(xué)的研究能力,該技術(shù)的原理是將植物被病毒侵染產(chǎn)生的源于病毒的特異小RNA分子序列進(jìn)行拼接,以此來(lái)鑒定病毒種類[18]。此技術(shù)可快速、靈敏、全面檢測(cè)植物體內(nèi)已知和未知病毒的核酸分子,被廣泛用于農(nóng)作物、藥用、園林和園藝植物病毒病的檢測(cè)和鑒定[11,18-22]。但是,它用于單個(gè)樣品的檢測(cè)成本較高,也無(wú)法確定混合樣品中感染病毒的具體樣品,還會(huì)因較短的小RNA分子序列使比對(duì)結(jié)果產(chǎn)生一定程度的誤判。因此,仍需要結(jié)合ELISA和RT-PCR等檢測(cè)方法進(jìn)行使用[23]。

    本研究對(duì)來(lái)源于北京、河北、遼寧、內(nèi)蒙古、陜西、四川和云南共7個(gè)省(市、自治區(qū))的草莓種苗進(jìn)行小RNA深度測(cè)序和生物信息學(xué)分析,并利用病毒特異性引物對(duì)檢測(cè)結(jié)果進(jìn)行RT-PCR驗(yàn)證。結(jié)果表明,侵染我國(guó)草莓種苗的病毒主要是SMoV、SVBV和SMYEV,且存在2種病毒和3種病毒復(fù)合侵染的情況。SMoV不同分離物的大外殼蛋白(large coat protein,CPL)基因和依賴RNA的RNA聚合酶(RNA-dependent RNA polymerase,RdRP)基因序列具有較大的差異,而3′非編碼區(qū)的核苷酸序列較保守[24]。因此,測(cè)定了SMoV的RNA1 3′端非編碼區(qū)部分核苷酸序列。經(jīng)系統(tǒng)進(jìn)化樹分析,本文獲得的SMoV四川分離物schhy13和河北分離物hbhy28與遼寧分離物lnhy23形成了兩個(gè)不同的分支。SVBV中國(guó)不同分離物的核苷酸序列保守性相對(duì)較高,且ORF6編碼的蛋白與寄主范圍和癥狀表現(xiàn)有關(guān)[5,14]。因此,測(cè)定了不同省市SVBV分離物ORF6的部分基因序列。本文獲得的5個(gè)SVBV分離物與GenBank中已登錄的中國(guó)分離物形成一個(gè)獨(dú)立的組群,與美國(guó)和加拿大分離物親緣關(guān)系較遠(yuǎn)。上述SMoV和SVBV的分組差異可能與地理起源有關(guān)[24]。另外,僅在一個(gè)來(lái)自遼寧的樣品中檢測(cè)到SMYEV,通過(guò)測(cè)定其完整的CP基因及該基因兩端部分非編碼區(qū)核苷酸序列,SMYEV遼寧分離物lnhy26 的CP基因序列可能是與中國(guó)甘肅分離物sy02同屬一個(gè)組群的新株系類型[6]。

    被病毒單獨(dú)侵染或2~3種病毒復(fù)合侵染的草莓植株可單獨(dú)或同時(shí)表現(xiàn)出葉片畸形、葉片變色和植株矮化等癥狀,調(diào)查表明田間草莓植株表現(xiàn)出的癥狀與感染病毒的種類關(guān)系不明確。此現(xiàn)象的產(chǎn)生可能是病毒不同株系侵染后的癥狀略有不同;也可能是RT-PCR無(wú)法檢出含量過(guò)低的病毒,使得檢測(cè)結(jié)果不夠準(zhǔn)確;或是與水肥、藥害、機(jī)械損傷等非病毒因素共同危害產(chǎn)生的[1,25]。在后續(xù)的研究中,可結(jié)合實(shí)時(shí)熒光定量PCR檢測(cè)技術(shù)提高檢測(cè)效率,使檢測(cè)結(jié)果更加準(zhǔn)確[1,26-27]。

    不同省市的草莓種苗病毒檢出情況差異顯著,在四川、遼寧、河北和陜西種苗中檢出病毒,而北京、內(nèi)蒙古和云南種苗中未檢出。其中,遼寧種苗的病毒檢出率最高、病毒種類最多,可能與遼寧發(fā)展草莓種植業(yè)的時(shí)間早、規(guī)模大有關(guān)[28]。對(duì)不同品種進(jìn)行病毒檢測(cè),僅在‘紅顏中檢出病毒,而‘白雪公主 ‘隋珠‘光點(diǎn)和‘圣誕紅均未檢出。可能是‘紅顏?zhàn)鳛橹匾耐茝V品種,因種植時(shí)間長(zhǎng)、種植范圍廣導(dǎo)致病毒病發(fā)生嚴(yán)重[29]。本研究對(duì)不同產(chǎn)地草莓種苗檢測(cè)的品種種類和數(shù)量不同,也可能對(duì)檢測(cè)結(jié)果產(chǎn)生影響,今后會(huì)增加檢測(cè)數(shù)量,進(jìn)一步明確我國(guó)不同省市草莓種苗病毒病發(fā)生情況。

    參考文獻(xiàn)

    [1] 王佳, 周宇, 祝寧, 等. 我國(guó)部分省市草莓種苗中黃瓜花葉病毒的檢測(cè)分析[J]. 植物病理學(xué)報(bào), 2020, 50(6): 649-656.

    [2] 聯(lián)合國(guó)糧食及農(nóng)業(yè)組織. 糧食和農(nóng)業(yè)數(shù)據(jù)[EB/OL]. [2021-01-30]http:∥www.fao.org/faostat/zh/#data/QC.

    [3] 楊肖芳, 苗立祥, 張?jiān)コ?等. 浙江省紅顏草莓鑲脈病毒(SVBV)的調(diào)查與莖尖脫毒技術(shù)研究[J]. 核農(nóng)學(xué)報(bào), 2015, 29(9): 1694-1700.

    [4] 王麗紅, 馬俊文, 肖紅, 等. 北方保護(hù)地草莓病害防治措施[J]. 中國(guó)果菜, 2019, 39(12): 95-98.1617EC4E-91A2-4AC0-96D3-3A8AA5458DFA

    [5] 尚巧霞. 介體傳播的草莓病毒研究進(jìn)展[J]. 植物保護(hù)學(xué)報(bào), 2015, 42(4): 488-496.

    [6] 楊洪一, 代紅艷, 李麗麗, 等. 草莓輕型黃邊病毒3′末端序列多態(tài)性研究[J]. 植物病理學(xué)報(bào), 2008, 38(1): 96-99.

    [7] CHEN Liu, SHANG Qiaoxia, CHEN Xiaoyu, et al. First report on the occurrence of cucumber mosaic virus on Fragaria ananassa in China [J]. Plant Disease, 2014, 98(7): 1015.

    [8] DING Xinlun, CHEN Dao, WANG Aiming, et al. First report of strawberry pallidosis-associated virus on strawberry in China [J]. Plant Disease, 2017, 101(12): 2154.

    [9] 湯亞飛, 裴凡, 李正剛, 等. 基于小RNA深度測(cè)序技術(shù)鑒定侵染廣東辣椒的病毒種類[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2019, 52(13): 2256-2267.

    [10]謝雪花. 小RNA深度測(cè)序鑒定草莓病毒的研究及草莓快繁體系的優(yōu)化[D]. 長(zhǎng)沙: 湖南農(nóng)業(yè)大學(xué), 2016.

    [11]LIU Qili, LI Yongqiang, SONG Puwen, et al. Diversity of viruses infecting sweet potato in Beijing based on small RNA deep sequencing and PCR or RT-PCR detection [J]. Journal of General Plant Pathology, 2020, 86: 283-289.

    [12]陳柳, 邢冬梅, 劉正雄, 等. 京郊設(shè)施蔬菜黃瓜花葉病毒的檢測(cè)[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2015, 31(20): 82-85.

    [13]褚明昕, 魏然, 席昕, 等. 北京地區(qū)不同草莓品種5種主要病毒的檢測(cè)[C]∥彭友良, 李向東. 中國(guó)植物病理學(xué)會(huì)2017年學(xué)術(shù)年會(huì)論文集. 北京: 中國(guó)農(nóng)業(yè)科技出版社, 2017: 217.

    [14]冉策. 京郊草莓病毒病調(diào)查及病毒檢測(cè)技術(shù)應(yīng)用[D]. 北京: 北京農(nóng)學(xué)院, 2015.

    [15]THOMPSON J R, JELKMANN W. Strain diversity and conserved genome elements in strawberry mild yellow edge virus [J]. Archives of Virology, 2004, 149(10): 1897-1909.

    [16]王佳, 崔高峰, 祝寧, 等. 北京地區(qū)設(shè)施草莓8種病毒的酶聯(lián)免疫吸附測(cè)定檢測(cè)[J]. 微生物學(xué)通報(bào), 2021, 48(3): 1048-1056.

    [17]陳柳, 冉策, 尚巧霞, 等. 草莓輕型黃邊病毒北京分離物的CP基因序列分析[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2015, 31(3): 229-233.

    [18]范旭東, 董雅鳳, 張尊平, 等. ‘陽(yáng)光玫瑰葡萄病毒小RNA測(cè)序鑒定及RT-PCR檢測(cè)[J]. 植物病理學(xué)報(bào), 2019, 49(6): 749-755.

    [19]錢亞娟, 徐毅, 周琦, 等. 利用深度測(cè)序技術(shù)發(fā)掘植物病毒資源[J]. 中國(guó)科學(xué): 生命科學(xué), 2014, 44(4): 351-363.

    [20]崔麗艷, 龐小靜, 齊永紅, 等. 利用小RNA深度測(cè)序技術(shù)鑒定半夏病毒病病原[J]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào), 2018, 34(12): 1334-1341.

    [21]蘇秀, 徐毅, 陳莎, 等. 利用小RNA深度測(cè)序和組裝技術(shù)鑒定紫藤花葉病病原[J]. 植物病理學(xué)報(bào), 2015, 45(1): 88-92.

    [22]HE Yan, YANG Zuokun, HONG Ni, et al. Deep sequencing reveals a novel closterovirus associated with wild rose leaf rosette disease [J]. Molecular Plant Pathology, 2015, 16(5): 449-458.

    [23]劉曉璇, 李威霖, 張彤, 等. 番茄斑駁花葉病毒海南分離物全基因組序列分析[J]. 熱帶作物學(xué)報(bào), 2019, 40(1): 139-144.

    [24]楊洪一, 李麗麗, 代紅艷, 等. 草莓斑駁病毒分子變異及PCR檢測(cè)技術(shù)研究[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41(7): 1983-1988.

    [25]MRAZ I, PETRZIK K, FRNOV-HONETSLEGROV J, et al. Detection of strawberry vein banding virus by polymerase chain reaction and dot blot hybridization [J]. Acta Virologica, 1997, 41(4): 241-242.

    [26]苗立祥, 榮寧寧, 張?jiān)コ?等. 實(shí)時(shí)熒光定量PCR法檢測(cè)草莓鑲脈病毒[J]. 浙江農(nóng)業(yè)學(xué)報(bào), 2016, 28(6): 1030-1036.

    [27]王建輝, 劉建軍, 陳克玲, 等. 草莓輕型黃邊病毒檢測(cè)研究[J]. 西南農(nóng)業(yè)學(xué)報(bào), 2016, 29(12): 2866-2870.

    [28]王玉坤, 張放, 祝庭耀. 國(guó)內(nèi)草莓生產(chǎn)現(xiàn)狀與發(fā)展趨勢(shì)[J]. 北方園藝, 2003(6): 6-7.

    [29]常琳琳, 董靜, 鐘傳飛, 等. 中國(guó)育成草莓品種的系譜分析[J]. 果樹學(xué)報(bào), 2018, 35(2): 158-167.

    (責(zé)任編輯:田 喆)1617EC4E-91A2-4AC0-96D3-3A8AA5458DFA

    久久久精品大字幕| 一边亲一边摸免费视频| 亚洲一级一片aⅴ在线观看| 久久这里只有精品中国| 欧美一区二区国产精品久久精品| 午夜福利在线观看吧| 热99re8久久精品国产| 色视频www国产| 亚洲图色成人| 不卡视频在线观看欧美| 亚洲精品自拍成人| 国产私拍福利视频在线观看| 99国产精品一区二区蜜桃av| 麻豆av噜噜一区二区三区| 99久久成人亚洲精品观看| 成人美女网站在线观看视频| 国语自产精品视频在线第100页| 99视频精品全部免费 在线| 欧美又色又爽又黄视频| 久久国产乱子免费精品| 久久久久久久久大av| 一级二级三级毛片免费看| 99久久精品国产国产毛片| 女人被狂操c到高潮| 级片在线观看| 国产精品国产高清国产av| 91久久精品电影网| 精品午夜福利在线看| 欧美激情在线99| 精品无人区乱码1区二区| 啦啦啦观看免费观看视频高清| 国产高清不卡午夜福利| 久久久成人免费电影| 精品免费久久久久久久清纯| 免费看av在线观看网站| 中国美白少妇内射xxxbb| 中文字幕av成人在线电影| 五月玫瑰六月丁香| 久久久色成人| 国产探花在线观看一区二区| 日本黄大片高清| 精品久久久久久电影网 | 中文字幕免费在线视频6| 在线播放无遮挡| 午夜福利网站1000一区二区三区| 成年女人永久免费观看视频| 国产成人aa在线观看| 国产精品麻豆人妻色哟哟久久 | 精品久久久噜噜| h日本视频在线播放| 午夜福利在线观看免费完整高清在| 最新中文字幕久久久久| 免费大片18禁| 欧美成人免费av一区二区三区| 亚洲国产日韩欧美精品在线观看| 蜜臀久久99精品久久宅男| 变态另类丝袜制服| 日韩 亚洲 欧美在线| 国产精品永久免费网站| 精品国产一区二区三区久久久樱花 | 精品酒店卫生间| 国产极品天堂在线| 成人特级av手机在线观看| 国产老妇伦熟女老妇高清| 十八禁国产超污无遮挡网站| 欧美三级亚洲精品| 中文精品一卡2卡3卡4更新| 亚洲精品影视一区二区三区av| 亚洲国产成人一精品久久久| 国产精品1区2区在线观看.| av免费观看日本| 精品免费久久久久久久清纯| 亚洲人成网站在线播| 久久鲁丝午夜福利片| 村上凉子中文字幕在线| 91在线精品国自产拍蜜月| 秋霞在线观看毛片| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日韩精品成人综合77777| 欧美性感艳星| 欧美极品一区二区三区四区| 老司机影院毛片| 久久亚洲精品不卡| 久久久久九九精品影院| 黄色欧美视频在线观看| 精品人妻偷拍中文字幕| 国产精品精品国产色婷婷| 一级黄片播放器| 日韩欧美在线乱码| 日本色播在线视频| av国产久精品久网站免费入址| 亚洲久久久久久中文字幕| www.av在线官网国产| 3wmmmm亚洲av在线观看| 欧美一区二区国产精品久久精品| 国产精品一区二区三区四区久久| 欧美日本亚洲视频在线播放| 精品少妇黑人巨大在线播放 | 亚洲成人精品中文字幕电影| 日韩在线高清观看一区二区三区| 国产 一区 欧美 日韩| 亚洲丝袜综合中文字幕| 国产午夜精品论理片| 国产精品一区二区在线观看99 | 日韩在线高清观看一区二区三区| 女人久久www免费人成看片 | 午夜福利视频1000在线观看| 三级毛片av免费| 亚洲国产成人一精品久久久| 久久人妻av系列| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩无卡精品| 天天躁日日操中文字幕| 天堂网av新在线| 插逼视频在线观看| 高清av免费在线| 成人性生交大片免费视频hd| 18+在线观看网站| 日本与韩国留学比较| 18禁在线播放成人免费| 国产乱人偷精品视频| 黄色日韩在线| 天堂√8在线中文| 人妻少妇偷人精品九色| 熟妇人妻久久中文字幕3abv| 亚洲伊人久久精品综合 | 特级一级黄色大片| 免费av毛片视频| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站高清观看| av在线观看视频网站免费| 少妇裸体淫交视频免费看高清| 亚洲av不卡在线观看| 亚洲色图av天堂| 精品一区二区免费观看| 99在线视频只有这里精品首页| 久久这里有精品视频免费| 日产精品乱码卡一卡2卡三| 午夜亚洲福利在线播放| videos熟女内射| 精品一区二区三区人妻视频| 一夜夜www| 欧美高清性xxxxhd video| 最近中文字幕2019免费版| 成人毛片a级毛片在线播放| 婷婷六月久久综合丁香| 99久久精品国产国产毛片| 秋霞伦理黄片| www日本黄色视频网| 毛片一级片免费看久久久久| 免费人成在线观看视频色| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区三区| 青春草国产在线视频| 亚洲欧美日韩东京热| 久久亚洲精品不卡| 51国产日韩欧美| 人人妻人人看人人澡| 男人舔奶头视频| 国产视频首页在线观看| 高清在线视频一区二区三区 | 在线观看美女被高潮喷水网站| 高清av免费在线| 国产美女午夜福利| 国产精品国产高清国产av| 亚洲国产精品专区欧美| 青春草国产在线视频| 国产精品乱码一区二三区的特点| 国产成年人精品一区二区| 少妇人妻精品综合一区二区| 亚洲成色77777| 久久精品久久久久久噜噜老黄 | 人妻夜夜爽99麻豆av| 国产成人福利小说| 精品久久久久久成人av| 日韩av不卡免费在线播放| 成人漫画全彩无遮挡| 99热这里只有精品一区| 国产乱人视频| 日韩欧美精品免费久久| 国产午夜精品论理片| 热99在线观看视频| 麻豆乱淫一区二区| 美女脱内裤让男人舔精品视频| 亚洲国产欧美人成| 九九爱精品视频在线观看| 国产一级毛片七仙女欲春2| 狂野欧美激情性xxxx在线观看| 一级av片app| 夫妻性生交免费视频一级片| 淫秽高清视频在线观看| 尤物成人国产欧美一区二区三区| 神马国产精品三级电影在线观看| 一本久久精品| 伊人久久精品亚洲午夜| 精品久久久久久久久av| 直男gayav资源| 国产精品综合久久久久久久免费| 丝袜美腿在线中文| 亚洲av中文字字幕乱码综合| 国产成人免费观看mmmm| 久久99热6这里只有精品| 天天一区二区日本电影三级| 有码 亚洲区| 午夜福利成人在线免费观看| 中文字幕av成人在线电影| 两个人视频免费观看高清| 国产色爽女视频免费观看| 汤姆久久久久久久影院中文字幕 | 久久国内精品自在自线图片| 看免费成人av毛片| 丰满乱子伦码专区| 日本三级黄在线观看| 你懂的网址亚洲精品在线观看 | 网址你懂的国产日韩在线| 亚洲av中文av极速乱| 欧美一区二区国产精品久久精品| 亚洲天堂国产精品一区在线| 午夜福利在线观看免费完整高清在| 在线免费观看的www视频| 亚洲欧美中文字幕日韩二区| 男女那种视频在线观看| 久久精品熟女亚洲av麻豆精品 | 18禁动态无遮挡网站| 久久久久久久亚洲中文字幕| 国产精品久久久久久久电影| 我的女老师完整版在线观看| 一区二区三区四区激情视频| 人妻系列 视频| 中文字幕av在线有码专区| 日韩欧美 国产精品| 亚州av有码| 国产精品久久视频播放| 高清毛片免费看| 国产毛片a区久久久久| 少妇被粗大猛烈的视频| 天美传媒精品一区二区| 干丝袜人妻中文字幕| 内射极品少妇av片p| 日本与韩国留学比较| 免费观看人在逋| 亚洲欧美一区二区三区国产| 黄色一级大片看看| 99九九线精品视频在线观看视频| 联通29元200g的流量卡| 最近最新中文字幕免费大全7| 韩国av在线不卡| 欧美区成人在线视频| 亚洲av电影在线观看一区二区三区 | 亚洲精品国产成人久久av| 国产精品永久免费网站| av国产免费在线观看| 国产精华一区二区三区| 欧美最新免费一区二区三区| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| 国产精品av视频在线免费观看| 午夜激情福利司机影院| 永久免费av网站大全| av视频在线观看入口| 中文字幕久久专区| 搡女人真爽免费视频火全软件| 熟女人妻精品中文字幕| 欧美三级亚洲精品| 一边摸一边抽搐一进一小说| 亚洲欧美日韩东京热| 男女边吃奶边做爰视频| 老女人水多毛片| 欧美日韩在线观看h| 日韩亚洲欧美综合| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 久久久精品94久久精品| 春色校园在线视频观看| 国产又黄又爽又无遮挡在线| 成人综合一区亚洲| h日本视频在线播放| 日韩av在线免费看完整版不卡| 国产一区亚洲一区在线观看| 麻豆国产97在线/欧美| 国内精品宾馆在线| 又黄又爽又刺激的免费视频.| 青春草视频在线免费观看| 在线观看66精品国产| 国产精品久久久久久精品电影| 最后的刺客免费高清国语| 十八禁国产超污无遮挡网站| 成年女人看的毛片在线观看| 99在线人妻在线中文字幕| 69av精品久久久久久| 18禁裸乳无遮挡免费网站照片| 三级国产精品片| 91精品一卡2卡3卡4卡| 婷婷色麻豆天堂久久 | 99热这里只有精品一区| 亚洲av二区三区四区| 直男gayav资源| 亚洲电影在线观看av| 我的女老师完整版在线观看| 精品欧美国产一区二区三| 看片在线看免费视频| 日产精品乱码卡一卡2卡三| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清专用| 亚洲成色77777| 青青草视频在线视频观看| 国产淫语在线视频| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 婷婷色麻豆天堂久久 | 中文欧美无线码| 国产成人aa在线观看| 欧美三级亚洲精品| 亚洲精品国产av成人精品| 女人久久www免费人成看片 | 午夜福利在线观看免费完整高清在| 国产成人a∨麻豆精品| 国产伦一二天堂av在线观看| 亚洲最大成人av| 久久久精品大字幕| av又黄又爽大尺度在线免费看 | 蜜臀久久99精品久久宅男| av女优亚洲男人天堂| 日本免费在线观看一区| 亚洲高清免费不卡视频| 三级经典国产精品| 两个人视频免费观看高清| 色5月婷婷丁香| 国产成人福利小说| av专区在线播放| 五月玫瑰六月丁香| 天堂网av新在线| 国产私拍福利视频在线观看| 国产亚洲av片在线观看秒播厂 | 五月玫瑰六月丁香| 午夜精品一区二区三区免费看| 亚洲国产精品专区欧美| 又爽又黄a免费视频| 精华霜和精华液先用哪个| 色视频www国产| 亚洲精品日韩在线中文字幕| 久久精品国产亚洲网站| 人人妻人人澡欧美一区二区| 一级黄片播放器| 男女啪啪激烈高潮av片| 国产麻豆成人av免费视频| 午夜福利在线观看免费完整高清在| 好男人视频免费观看在线| 99久久九九国产精品国产免费| 搡老妇女老女人老熟妇| 欧美日本视频| 乱人视频在线观看| 国产伦在线观看视频一区| 97在线视频观看| 啦啦啦观看免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 国产又黄又爽又无遮挡在线| 亚州av有码| 观看免费一级毛片| 少妇被粗大猛烈的视频| av女优亚洲男人天堂| 高清在线视频一区二区三区 | 少妇人妻精品综合一区二区| 免费观看的影片在线观看| 免费av观看视频| 亚洲欧美日韩无卡精品| 男人的好看免费观看在线视频| 赤兔流量卡办理| 干丝袜人妻中文字幕| 熟妇人妻久久中文字幕3abv| 哪个播放器可以免费观看大片| 最近最新中文字幕免费大全7| 午夜福利在线观看吧| 免费av观看视频| 久久久成人免费电影| 99热这里只有是精品在线观看| 国产伦精品一区二区三区四那| 高清日韩中文字幕在线| 成人二区视频| 国产女主播在线喷水免费视频网站 | 波野结衣二区三区在线| 婷婷六月久久综合丁香| 久久久国产成人免费| 成年版毛片免费区| 久久人人爽人人爽人人片va| 亚洲久久久久久中文字幕| 国产爱豆传媒在线观看| 丰满少妇做爰视频| 日本一本二区三区精品| 男女啪啪激烈高潮av片| av国产久精品久网站免费入址| 精品久久久久久久人妻蜜臀av| 最新中文字幕久久久久| 亚洲国产精品成人久久小说| 岛国毛片在线播放| 国产黄a三级三级三级人| 免费看日本二区| 亚洲中文字幕一区二区三区有码在线看| 最近视频中文字幕2019在线8| 国产久久久一区二区三区| 91午夜精品亚洲一区二区三区| 久久午夜福利片| 国产精品电影一区二区三区| 色噜噜av男人的天堂激情| 国产精品久久久久久精品电影小说 | 黄色欧美视频在线观看| 亚洲精品日韩av片在线观看| 国产亚洲一区二区精品| 日韩欧美精品v在线| 国产一级毛片在线| 男人舔女人下体高潮全视频| 日韩欧美国产在线观看| 久久人人爽人人爽人人片va| 中文字幕av在线有码专区| 中国国产av一级| 中文资源天堂在线| 免费观看性生交大片5| 不卡视频在线观看欧美| 18禁在线无遮挡免费观看视频| 国产高清不卡午夜福利| 日韩亚洲欧美综合| 日韩精品青青久久久久久| 亚洲最大成人手机在线| 熟妇人妻久久中文字幕3abv| av又黄又爽大尺度在线免费看 | 国产免费又黄又爽又色| 人人妻人人澡欧美一区二区| 久久久色成人| 波野结衣二区三区在线| 青青草视频在线视频观看| av黄色大香蕉| 国产午夜精品久久久久久一区二区三区| 美女cb高潮喷水在线观看| av国产久精品久网站免费入址| 综合色丁香网| 18禁动态无遮挡网站| 免费观看a级毛片全部| 亚洲aⅴ乱码一区二区在线播放| 亚洲美女搞黄在线观看| 国产av在哪里看| 国产亚洲精品久久久com| 黄色日韩在线| 日韩成人av中文字幕在线观看| 级片在线观看| 久久亚洲精品不卡| 免费黄网站久久成人精品| 精品久久久久久久久亚洲| 亚洲一区高清亚洲精品| 色网站视频免费| 精品久久久久久久末码| 国产白丝娇喘喷水9色精品| 午夜亚洲福利在线播放| 久久精品国产鲁丝片午夜精品| www日本黄色视频网| 国产视频内射| 日韩人妻高清精品专区| 精品欧美国产一区二区三| h日本视频在线播放| 69人妻影院| 欧美日本视频| 成人特级av手机在线观看| 日韩视频在线欧美| 在线免费观看的www视频| 久久久久久久久久久丰满| 看非洲黑人一级黄片| 建设人人有责人人尽责人人享有的 | 建设人人有责人人尽责人人享有的 | 欧美另类亚洲清纯唯美| 国产精品一区二区性色av| 永久网站在线| 欧美性感艳星| 国产在线男女| 国产精品国产三级国产av玫瑰| 国产视频内射| 日韩人妻高清精品专区| 亚洲国产色片| 国产精品一区二区在线观看99 | 午夜福利网站1000一区二区三区| 2021少妇久久久久久久久久久| 久久久久久久久久久免费av| 国产黄色小视频在线观看| 在线播放无遮挡| 久久久精品94久久精品| 国产探花在线观看一区二区| 国产亚洲91精品色在线| 久久久久久久亚洲中文字幕| 亚洲av.av天堂| 久久亚洲精品不卡| 一二三四中文在线观看免费高清| 久久亚洲国产成人精品v| 国产一区二区亚洲精品在线观看| 真实男女啪啪啪动态图| 亚洲国产精品专区欧美| 精品久久久久久成人av| 亚洲欧美精品专区久久| 一级黄片播放器| 亚洲精品乱码久久久久久按摩| 国产老妇伦熟女老妇高清| 亚洲国产高清在线一区二区三| 亚洲国产最新在线播放| av天堂中文字幕网| 一二三四中文在线观看免费高清| 国内精品一区二区在线观看| 成人午夜高清在线视频| 国产探花极品一区二区| or卡值多少钱| 国产白丝娇喘喷水9色精品| 欧美色视频一区免费| 亚洲真实伦在线观看| 直男gayav资源| 亚洲国产欧美人成| 国内精品一区二区在线观看| 99久国产av精品国产电影| 亚洲精品久久久久久婷婷小说 | 国产精品1区2区在线观看.| 亚洲国产色片| 日本五十路高清| 三级毛片av免费| 亚洲三级黄色毛片| 欧美一区二区精品小视频在线| 2021少妇久久久久久久久久久| 精品国产露脸久久av麻豆 | 日本免费一区二区三区高清不卡| 禁无遮挡网站| 亚洲国产高清在线一区二区三| 99久久精品国产国产毛片| av在线蜜桃| 一级黄片播放器| 精品人妻视频免费看| 成人亚洲精品av一区二区| 国产真实伦视频高清在线观看| 99热全是精品| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看 | 国产淫片久久久久久久久| 国产精品一区二区性色av| 国产老妇伦熟女老妇高清| 纵有疾风起免费观看全集完整版 | 日韩av在线大香蕉| 亚洲欧美精品自产自拍| 99热这里只有是精品在线观看| 亚洲av男天堂| 尾随美女入室| 美女cb高潮喷水在线观看| 亚洲欧美精品专区久久| 国产色爽女视频免费观看| 偷拍熟女少妇极品色| 亚洲中文字幕日韩| 欧美3d第一页| 日韩三级伦理在线观看| 亚洲在久久综合| 国产亚洲av嫩草精品影院| 2021天堂中文幕一二区在线观| 国产在线一区二区三区精 | 国产精品麻豆人妻色哟哟久久 | 如何舔出高潮| 九色成人免费人妻av| 国内精品一区二区在线观看| 国产人妻一区二区三区在| 国内揄拍国产精品人妻在线| 联通29元200g的流量卡| 亚洲国产精品成人综合色| 国产高潮美女av| 蜜桃久久精品国产亚洲av| 久久久久久久午夜电影| 国产精品.久久久| 在现免费观看毛片| 国产高潮美女av| 韩国av在线不卡| 久久久久久久久久久丰满| 村上凉子中文字幕在线| 欧美丝袜亚洲另类| 精品人妻视频免费看| 婷婷色av中文字幕| 91aial.com中文字幕在线观看| 婷婷色av中文字幕| 免费播放大片免费观看视频在线观看 | av女优亚洲男人天堂| 中文字幕av成人在线电影| 青春草国产在线视频| 亚洲欧美日韩无卡精品| 国产精品麻豆人妻色哟哟久久 | 中文字幕制服av| 成年女人永久免费观看视频| 一个人看的www免费观看视频| 亚洲人与动物交配视频| 国产精品一及| 不卡视频在线观看欧美| 熟女人妻精品中文字幕| 99热网站在线观看| 一个人看视频在线观看www免费| 一本一本综合久久| 午夜a级毛片| 欧美变态另类bdsm刘玥| 一级毛片我不卡| 哪个播放器可以免费观看大片| 国产淫语在线视频| 好男人视频免费观看在线| 村上凉子中文字幕在线| 啦啦啦韩国在线观看视频| 卡戴珊不雅视频在线播放| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久久久久久久| 国产v大片淫在线免费观看| 国产高清视频在线观看网站| 国产精品.久久久| 麻豆乱淫一区二区| 免费不卡的大黄色大毛片视频在线观看 | 精品酒店卫生间| 韩国av在线不卡| 中文精品一卡2卡3卡4更新| h日本视频在线播放| 久久精品久久久久久久性| 丰满人妻一区二区三区视频av| 少妇熟女aⅴ在线视频|