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

    多主棒孢SdhB-H278R突變位點AS-real-time PCR定量檢測體系的建立

    2018-12-27 02:06:22孫炳學(xué)石延霞朱發(fā)娣謝學(xué)文柴阿麗李寶聚
    中國農(nóng)業(yè)科學(xué) 2018年24期
    關(guān)鍵詞:檢測

    孫炳學(xué),石延霞,朱發(fā)娣,謝學(xué)文,柴阿麗,李寶聚

    ?

    多主棒孢SdhB-H278R突變位點AS-real-time PCR定量檢測體系的建立

    孫炳學(xué),石延霞,朱發(fā)娣,謝學(xué)文,柴阿麗,李寶聚

    (中國農(nóng)業(yè)科學(xué)院蔬菜花卉研究所,北京 100081)

    【目的】建立一種快速、高效、定量檢測黃瓜多主棒孢()琥珀酸脫氫酶B亞基()H278R突變的實時熒光定量PCR(AS-real-time PCR)檢測方法,并進行效果驗證。【方法】從北京市大興區(qū)采集、分離純化病原菌,獲得24株多主棒孢的單孢菌株,采用菌絲生長速率法測定其對啶酰菌胺的EC50,隨機選取4株敏感(S)和8株抗性(R)菌株測其菌絲生長速率、產(chǎn)孢量和致病性。采用引物對Cc-SdhB-F/R測定多主棒孢基因序列,檢測的堿基變化?;诙嘀靼翩叩南嗤怂嵝蛄泻蚐NP位點設(shè)計內(nèi)參引物B-H278R-TY-F/R和特異性引物B-H278R-2F/2R14,建立并優(yōu)化AS-real-time PCR定量檢測體系,并對引物的特異性、熔解曲線和靈敏度進行評價。利用該體系分別檢測含有H278R不同突變比例的DNA和孢子懸浮液?!窘Y(jié)果】測定的24個菌株中,敏感菌株占66.7%,EC50值為0.057—0.563 μg·mL-1;抗性菌株占33.3%,EC50值為5.395—11.710 μg·mL-1??剐院兔舾芯陜H在菌絲生長速率方面存在顯著性差異,且菌絲生長速率與EC50之間存在顯著負相關(guān)。在產(chǎn)孢量和致病性方面不存在顯著性差異。測序分析發(fā)現(xiàn)抗性菌株均攜帶SdhB-H278R突變。本研究設(shè)計的引物B-H278R-2F和B-H278R-2F2R14特異性強,僅對多主棒孢SdhB-H278R突變菌株的DNA有擴增條帶,對其余供試菌株均無條帶擴增。普通AS-PCR檢測的靈敏度為91 pg·μL-1, 而AS-real-time PCR的靈敏度可達9.1 pg·μL-1,靈敏度為普通AS-PCR的10倍。以基因組DNA為標(biāo)準(zhǔn)品,構(gòu)建的AS-real-time PCR標(biāo)準(zhǔn)曲線ΔCT值與相對模板濃度的對數(shù)有良好的線性關(guān)系,相關(guān)系數(shù)為0.9857, 擴增效率為92.59%。熔解曲線吸收峰單一,內(nèi)參引物B-H278R-TY-F/R和特異性引物B-H278R-2F/2R14分別在87.81℃和91.62℃處出現(xiàn)單一特異性峰。用DNA和孢子懸浮液對標(biāo)準(zhǔn)曲線進行驗證,結(jié)果表明隨著相對模板濃度逐漸降低,該檢測體系的準(zhǔn)確性逐漸升高且檢測下限為5%,同時預(yù)期值與試驗值呈線性相關(guān)(2=0.9998和2=0.9922)?!窘Y(jié)論】建立了一種高效、定量、快速的AS-real-time PCR檢測體系用于SdhB-H278R突變位點的檢測,可為殺菌劑的抗性治理提供理論依據(jù)。

    黃瓜;多主棒孢;AS-real-time PCR;抗藥性;琥珀酸脫氫酶B亞基;啶酰菌胺

    0 引言

    【研究意義】多主棒孢()隸屬于格孢菌目(Pleosporales)棒孢科(Corynesporascaceae)棒孢屬(),由其引起的黃瓜棒孢葉斑病發(fā)生普遍,該病害在溫室、大棚和露地中均不斷加重,已成為黃瓜栽培中主要的病害之一[1]。目前大多采用化學(xué)藥劑防治黃瓜棒孢葉斑病,由于多主棒孢繁殖速度快、遺傳變異大,加之多年連續(xù)使用單一藥劑防治,導(dǎo)致其對苯并咪唑類、二羧酰亞胺類、-苯氨基甲酸酯類和QoI類的大部分殺菌劑產(chǎn)生了抗藥性[2-3]。啶酰菌胺是德國巴斯夫公司1996年開發(fā)的新型煙酰胺類內(nèi)吸性殺菌劑,該藥劑能夠阻斷電子在復(fù)合體II琥珀酸脫氫酶[3Fe-4S]結(jié)構(gòu)域與泛醌之間的傳遞,從而影響病原真菌的電子傳遞,阻礙病原菌的能量代謝[3-4]。因該藥劑作用位點單一,F(xiàn)RAC已將其歸為中度抗性風(fēng)險殺菌劑。自2008年啶酰菌胺在中國登記防治灰霉病并在黃瓜栽培中大量使用,間接導(dǎo)致多主棒孢對啶酰菌胺產(chǎn)生抗藥性,出現(xiàn)大量的抗性突變株,因此,建立一種快速、準(zhǔn)確、定量的多主棒孢SdhB-H278R突變檢測方法,對監(jiān)測啶酰菌胺田間的抗性頻率和病害防治具有重要意義。【前人研究進展】隨著啶酰菌胺的廣泛使用,在不同蔬菜和作物上出現(xiàn)了對啶酰菌胺有抗性的病原菌。2007年,首次在開心果上發(fā)現(xiàn)了對啶酰菌胺產(chǎn)生抗性的鏈格孢()[5],之后在灰葡萄孢()、小雙胞腔菌()和核盤菌()上也發(fā)現(xiàn)了對啶酰菌胺產(chǎn)生抗性的菌株[6-8]。2009年,Miyamoto等[9]在黃瓜上發(fā)現(xiàn)了對啶酰菌胺產(chǎn)生抗性的多主棒孢并在2010[10]年發(fā)現(xiàn)了SdhB-H278R/Y、SdhC-C73P、SdhD-S89P和SdhD-G109V等多種抗性突變類型,且SdhB-H278R基因型在田間的頻率為0.64%。傳統(tǒng)的抗藥性檢測方法包括離體測定和活體測定法,其中離體測定法主要有菌落生長速率法和孢子萌發(fā)法,均是通過測定藥劑對植物病原菌的EC50值以區(qū)分敏感和抗性菌株,比較費時費力。近年來,隨著分子生物學(xué)技術(shù)的發(fā)展,PCR-RFLP(restriction fragment length polymorphism PCR)、AS-PCR(allele specific PCR)、AS-real-time PCR、HRM(high resolution melting)和多重等位基因PCR等分子鑒定手段廣泛應(yīng)用于植物病原菌對藥劑的抗藥性檢測。Furuya等[11]和Aoki等[12]在霜霉病菌()上分別建立了對QoI殺菌劑抗性Cytb-G143A和對CAA殺菌劑PvCesA3-G1105S突變位點PCR-RFLP的檢測體系。基于等位基因特異性核苷酸引物,AS-PCR和AS-real-time PCR抗藥性檢測技術(shù)已在核盤菌、茄鏈格孢()、大麥云紋病菌()和灰葡萄孢等病原菌上應(yīng)用[13-16]。最近,新興的HRM抗藥性檢測技術(shù)因具有高通量、快速的特點并可用于SNP位點的篩選,已開始在核盤菌對甲基硫菌靈和灰葡萄孢對環(huán)酰菌胺以及啶酰菌胺的抗藥性檢測上應(yīng)用[17-18],準(zhǔn)確性高達100%,但該方法只能定性檢測?!颈狙芯壳腥朦c】目前為止尚未建立一種檢測多主棒孢抗性突變位點的檢測體系,基于筆者課題組前期發(fā)現(xiàn)的一種抗啶酰菌胺的黃瓜多主棒孢SdhB-H278R抗性突變位點,使用熒光定量PCR儀,建立一種定量、高效、快速檢測該突變位點的AS-real-time PCR體系,并對此體系的特異性、可重復(fù)性和靈敏度進行評價?!緮M解決的關(guān)鍵問題】建立一種定量、高效、快速檢測黃瓜多主棒孢SdhB-H278R突變位點的AS-real-time PCR體系,以期為該病原菌對啶酰菌胺的田間抗性檢測和病害防控提供更好的技術(shù)支持。

    1 材料與方法

    試驗于2017年3月至2018年4月在中國農(nóng)業(yè)科學(xué)院蔬菜花卉研究所完成。

    1.1 試驗材料

    1.1.1 菌株 2017年3月從北京市大興區(qū)采集黃瓜多主棒孢葉斑病病樣,在室內(nèi)進行分離培養(yǎng),單孢純化后得到24株多主棒孢菌株,其他5株多主棒孢菌株采自河北和遼寧。所有菌株由中國農(nóng)業(yè)科學(xué)院蔬菜花卉研究所蔬菜病害綜合防治課題組提供。

    1.1.2 試劑和培養(yǎng)基 96%啶酰菌胺原藥(陜西美邦農(nóng)藥有限公司)、植物組織DNA提取試劑盒、SuperReal PreMix Plus(天根生化科技(北京)有限公司),2×Taq MasterMix (北京博邁德生物技術(shù)有限公司)、PDA培養(yǎng)基、YBA培養(yǎng)基。

    1.1.3 主要儀器設(shè)備 NanoDrop2000超微量分光光度計(Thermo Fisher Scientific Inc.,USA)、3-18K型離心機(Sigma,德國)、7500 Real time PCR system(ABI,美國)、S1000 Thermal Cycle PCR儀(BIO- RAD,美國)。

    1.1.4 引物合成測序 熒光定量PCR和普通PCR所用引物及基因序列由北京博邁德基因技術(shù)有限公司合成并測定。

    1.2 方法

    1.2.1 多主棒孢對啶酰菌胺敏感性的測定 參照Miyamoto等[9]方法,采用菌絲生長速率法測定北京大興24株及河北、遼寧的5株多主棒孢對啶酰菌胺的敏感性??顾幮詼y定設(shè)置7個濃度梯度(0、0.01、0.1、1、5、10、30 μg·mL-1),以不加藥劑為空白對照;在新鮮的菌落邊緣打取直徑5 mm的菌餅,接種于含藥培養(yǎng)基的平板上,在26℃黑暗培養(yǎng)5 d,測量菌落直徑,每個處理重復(fù)3次。比較試驗組和對照組的菌落直徑,并計算EC50值、毒力回歸方程和相關(guān)系數(shù)。根據(jù)Miyamoto等[9]的敏感基線(0.04—0.59 μg·mL-1)判定多主棒孢對啶酰菌胺的敏感性。

    1.2.2 抗性和敏感菌株的生物學(xué)特性測定 菌絲生長速率的測定:從24株多主棒孢菌中隨機選取4株敏感菌株和8株抗性菌株作為測定生物學(xué)特性的菌株。在多主棒孢菌生長5 d的PDA平板上制備5 mm的菌餅,采用滅菌牙簽將菌碟菌絲面朝下貼于PDA平板中央,密封后置于26℃培養(yǎng)。分別于第2—11天,采用十字交叉法測量菌落生長直徑。每個菌株3次重復(fù),試驗重復(fù)2次。

    產(chǎn)孢能力的測定:待多主棒孢菌生長至第11天,用直徑5 mm的打孔器在菌落邊緣打取5個菌餅置于1 mL的無菌水中,加入5個直徑2 mm的研磨珠用渦旋振蕩機振蕩2 min洗脫分生孢子,用血球計數(shù)板計數(shù)計算分生孢子的濃度。每個菌株3次重復(fù),試驗重復(fù)2次。

    致病性的測定:采集離體組織測定法,采集葉齡一致,形狀、大小相似的健康黃瓜葉片,用滅菌的接種針刺破形成微傷口。將制備好直徑5 mm的菌餅分別接種于黃瓜葉片,接種后放置于保鮮盒中,置于相對濕度85%,自然光照的26℃溫室中,培養(yǎng)5 d后采用十字交叉法測量病斑直徑。每個菌株測6片葉片,試驗重復(fù)2次。

    1.2.3 DNA提取 將多主棒孢接種到PDA平板培養(yǎng)5 d,刮取菌絲作為DNA提取樣本。采取植物組織DNA提取試劑盒,根據(jù)試劑盒中的說明書提取DNA,提取后用NanoDrop 2000c測定DNA濃度,用ddH2O稀釋至40 ng·μL-1,4℃保存?zhèn)溆谩?/p>

    1.2.4 AS-real-time PCR定量檢測體系的建立序列測定:根據(jù)GenBank多主棒孢序列(AB548739.1), 應(yīng)用Primer Premier 5.0設(shè)計引物Cc-SdhB-F和Cc-SdhB-R(表1),以表2中所有多主棒孢的DNA為模板,用Cc-SdhB-F和Cc-SdhB-R引物擴增基因片段,PCR產(chǎn)物進行測序。50 μL PCR體系:2×Taq PCR MasterMix 25 μL,Cc-SdhB-F和Cc-SdhB-R引物2 μL,模板DNA 2 μL,ddH2O補足50 μL。PCR反應(yīng)程序:94℃ 5 min,94℃ 45 s,58℃ 45 s,72℃ 90 s,步驟2—4進行30個循環(huán),72℃ 10 min。

    引物設(shè)計:根據(jù)的測序信息,針對SNP位點和相同的序列信息,用Primer Premier 5.0設(shè)計特異性引物B-H278R-2F/2R14和內(nèi)參引物B-H278R- TY-F/R(表1),引物設(shè)計策略和引物信息如圖1和表1所示。為增加引物的特異性,引物3′端與突變堿基G配對,在靠近3′端的第2、3位引入錯配堿基[13-16,19]。

    通過對全球自動駕駛技術(shù)領(lǐng)域?qū)@袌鲞M行研究,發(fā)現(xiàn)當(dāng)前自動駕駛技術(shù)專利申請數(shù)量逐年上升,研發(fā)熱情逐年高漲。但從技術(shù)角度出發(fā),目前從L3有條件自動駕駛階段邁向未來L5完全無人駕駛階段仍需要很長的過程和時間。本文從專利時間分布、專利空間分布、專利內(nèi)容分布這三大角度出發(fā),全面揭示了全球自動駕駛技術(shù)的發(fā)展趨勢、技術(shù)公開國家、技術(shù)來源國家、技術(shù)研發(fā)機構(gòu)、技術(shù)研究方向、研究熱點等情況,最終總結(jié)出以下6點結(jié)論,期望為自動駕駛技術(shù)的發(fā)展提供決策參考。

    圖1 引物設(shè)計策略

    表1 引物信息

    單下劃線Single underline:突變位點Mutation site;雙下劃線double underline:人為錯配堿基additional mismatched base

    引物特異性檢測:以含有H278R突變、其他點突變菌株以及常見的8種植物病原菌DNA為模板,以ddH2O為陰性對照,用引物B-H278R-2F和B-H278R-2R14進行PCR擴增,檢測引物特異性。20 μL PCR反應(yīng)體系:2×Taq PCR MasterMix 10 μL,引物B-H278R-2F和B-H278R-2R14 0.4 μL,模板DNA 1.6 μL,ddH2O補足10 μL。PCR反應(yīng)程序:94℃ 5 min,94℃ 45 s,65℃ 45 s,72℃ 30 s,步驟2—4進行35個循環(huán),72℃ 10 min。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測后,凝膠成像系統(tǒng)分析結(jié)果。

    標(biāo)準(zhǔn)曲線的建立:為了定量檢測多主棒孢SdhB- H278R突變株在田間的比例,以含有H278R突變菌株(R278)和不含有H278R突變菌株(H278)的混合DNA作為模板,以R278 DNA摩爾質(zhì)量分別占總量DNA的6.25%、12.5%、25%、50%、75%和100%作為不同的處理,3次重復(fù),分別用特異性引物B-H278R-2F/2R14和內(nèi)參引物B-H278R-TY-F/R在7500 real-time PCR system儀中進行AS-real-time PCR,計算ΔCT(ΔCT=CT特異性引物-CT內(nèi)參引物),以ΔCT值為縱坐標(biāo),以R278 DNA相對摩爾質(zhì)量的對數(shù)為橫坐標(biāo),制作標(biāo)準(zhǔn)曲線。

    20 μL反應(yīng)體系:2×SuperRealPreMix Plus 10 μL,50×ROX Reference Dye 0.4 μL,上游引物和下游引物各1 μL,模板DNA 40 pg,RNase-free ddH2O補足20 μL,AS-real-time PCR反應(yīng)程序:95℃ 15 min;95℃ 10 s;60℃ 32 s,40個循環(huán);95℃ 15 s;60℃ 1 min;60℃開始,每個循環(huán)增加0.35℃,100個循環(huán),60℃ 15 s。

    普通PCR和AS-real-time PCR靈敏度檢測:將定量的R278基因組DNA按濃度進行10倍梯度稀釋,用引物B-H278R-2F和B-H278R-2R14進行普通PCR和AS-real-time PCR擴增,檢測引物的靈敏度。

    1.2.5 AS-real-time PCR定量檢測體系的驗證 DNA驗證:以R278和H278混合DNA作為模板,以R278 DNA摩爾質(zhì)量分別占總量DNA的80%、40%、20%、10%和5%作為不同的處理,3次重復(fù)。按照已建立的AS-real-time PCR體系驗證標(biāo)準(zhǔn)曲線的可行性。用Excel計算標(biāo)準(zhǔn)差以及與預(yù)期值和檢測值之間的線性相關(guān)性。

    孢子懸浮液驗證:將R278和H278菌株26℃下培養(yǎng)11 d,用0.1%的吐溫20刷下孢子懸浮液并用三層紗布過濾得到孢子懸浮液。用血球計數(shù)板調(diào)節(jié)孢子濃度為1×105個/mL,配置混合孢子懸浮液,其中以H278孢子數(shù)量分別占混合孢子數(shù)的80%、40%、20%、10%和5%作為不同的處理,12 000 r/min離心,收集沉淀,用植物組織DNA提取試劑盒提取不同處理的DNA。按照已建立的AS-real-time PCR體系驗證標(biāo)準(zhǔn)曲線的可行性。用Excel計算標(biāo)準(zhǔn)差以及與預(yù)期值和檢測值之間的線性相關(guān)性。

    2 結(jié)果

    2.1 多主棒孢對啶酰菌胺的敏感性測定

    采用菌絲生長速率法測定采集的多主棒孢對啶酰菌胺的敏感性,結(jié)果表明,采自北京的24株多主棒孢分為敏感和抗性2個群體,其中敏感菌株16株,占66.7%,其EC50值為0.057—0.563 μg·mL-1,平均值為0.368 μg·mL-1??剐跃?株,占33.3%,其EC50值為5.395—11.710 μg·mL-1,平均值為7.586 μg·mL-1(圖2、表2)。

    2.2 抗性和敏感菌株的生物學(xué)特性

    隨機選取的12株多主棒孢之間的平均菌絲生長速率存在差異,敏感菌株菌絲生長速率為0.49—0.61 cm·d-1,抗性菌株菌絲生長速率為0.30—0.56 cm·d-1(表3)。與敏感菌株相比,抗性菌株除第9天的菌絲生長速率相等外,其他時間段均比敏感菌株生長慢(圖3),且抗性菌株的平均生長速率與敏感菌株差異極顯著(<0.01)。菌絲生長速率與EC50之間呈顯著負相關(guān)(表4)。

    由表3可見,敏感菌株的產(chǎn)孢量在0.33×105—7.00×105個/mL,抗性菌株的產(chǎn)孢量在1.33×105—7.46×105個/mL,抗性菌株和敏感菌株在產(chǎn)孢量上不存在顯著差異(0.949)。產(chǎn)孢量與生長速率和致病性存在低度相關(guān),但不顯著(表4)。

    敏感菌株的病斑直徑在0.38—0.84 cm,抗性菌株的病斑直徑在0.33—0.91 cm??剐院兔舾芯暝谥虏⌒陨喜淮嬖陲@著差異(表3)。致病性與產(chǎn)孢量和生長速率之間存在低度相關(guān)性,但不顯著(表4)。

    圖2 24株多主棒孢菌的敏感性分布

    圖3 抗性和敏感菌株菌絲生長速率比較

    表2 供試菌株敏感性測定

    +:有SdhB-H278R突變Carried SdhB-H278R mutation;-:無SdhB-H278R突變No SdhB-H278R mutation;R:抗性菌株resistant strain;S:敏感菌株sensitive strain

    2.3 AS-real-time PCR定量檢測體系的建立

    應(yīng)用引物Cc-SdhB-F和Cc-SdhB-R測定的序列,結(jié)果如圖4 所示,抗性菌株(R)在第956位的A突變到G,導(dǎo)致在密碼子的第278組氨酸(H)突變成精氨酸(R),進一步發(fā)現(xiàn)含有H278R突變的菌株對啶酰菌胺具有抗性。

    以攜帶不同點突變的多主棒孢和8株常見的病原菌基因組DNA為模板,采用引物B-H278R-2F和B-H278R-2R14進行特異性驗證(圖5),結(jié)果表明僅有攜帶SdhB-H278R突變的多主棒孢基因組DNA有擴增條帶,大小為244 bp,攜帶其他點突變的多主棒孢以及8種病原菌DNA均未檢測到擴增條帶,可用于不同菌種間多主棒孢H278R基因型的檢測。

    表3 抗性和敏感菌株的生物學(xué)特性

    表中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)差。表5同Data aremean±standard deviation. The same as Table 5根據(jù)Duncan’s最小顯著性差異測定,同列數(shù)值后不同字母表示差異顯著(<0.05)According to a Duncan’s least-signi?cant difference, different lowercases after data in the same column indicate significant difference (<0.05)

    圖4 在SdhB上抗性和敏感菌株的核酸序列和氨基酸序列信息

    M:分子量標(biāo)準(zhǔn)BM5000 Molecular weight marker BM5000;1:SdhB-H278R突變SdhB-H278R mutation;2:無突變no mutation;3:SdhB-H278Y突變SdhB-H278Y mutation;4:SdhB-I280V突變SdhB-I280V mutation;5:SdhB-P199S突變SdhB-P199S mutation;6:SdhC-S73P突變SdhC-S73P mutation;7:黃瓜葉片cucumber leaf;8:灰葡萄孢B. cinerea;9:瓜枝孢Cladosporium cucnerinu;10:立枯絲核菌Rhizoctonia solani;11:茄鏈格孢A. solani;12:茄葡柄霉Stemphylium solani;13:尖鐮孢Fusarium oxysporum;14:核盤菌S. sclerotiorum;15:古巴假霜霉菌Pseudoperonospora cubensis;16:ddH2O

    利用10倍梯度稀釋的基因組DNA對引物B-H278R-2F和B-H278R-2R14進行普通PCR和AS-real-time PCR靈敏度檢測,結(jié)果表明普通PCR檢測的靈敏度為91 pg·μL-1(圖6-A),AS-real-time PCR檢測的靈敏度為9.1 pg·μL-1(圖6-B)。AS-real-time PCR的靈敏度為普通PCR的10倍。

    表4 EC50與生物學(xué)特性的相關(guān)性分析

    表中數(shù)據(jù)為Pearson相關(guān)系數(shù)Data in the table are Pearson correlation coefficient values;*:<0.05

    A:普通PCR靈敏度檢測sensitivity detection with ordinary PCR;B:AS-real-time PCR靈敏度檢測sensitivity detection with AS-real-time PCR;M:分子量標(biāo)準(zhǔn)BM5000 Molecular weight marker BM5000;1—8:DNA稀釋為9.1、9.1×10-1、9.1×10-2、9.1×10-3、9.1×10-4、9.1×10-5、9.1×10-6、9.1×10-7 ng·μL-1 DNA diluted to 9.1, 9.1×10-1, 9.1×10-2, 9.1×10-3, 9.1×10-4, 9.1×10-5, 9.1×10-6, 9.1×10-7 ng·μL-1;9:ddH2O

    利用優(yōu)化的條件進行AS-real-time PCR,特異性引物B-H278R-2F/2R14擴增產(chǎn)物在熔解溫度為91.62℃處出現(xiàn)單一的特異性峰(圖7-A),內(nèi)參引物B-H278R- TY-F/R擴增產(chǎn)物在熔解溫度為87.81℃處出現(xiàn)單一的特異性峰(圖7-B)。

    以R278 DNA摩爾質(zhì)量分別占總量DNA的6.25%、12.5%、25%、50%、75%和100%作為不同的處理,用引物B-H278R-2F/2R14、B-H278R-TY/R進行AS- real-time PCR檢測,建立了標(biāo)準(zhǔn)曲線=-3.513+9.807,斜率為-3.513,截距為9.807,相關(guān)系數(shù)2=0.9857,擴增效率為92.59%。R278 DNA占總DNA的6.25%—100%的比例中,呈現(xiàn)出良好的線性關(guān)系(圖8)。

    2.4 AS-real-time PCR定量檢測體系的驗證

    以R278菌株DNA摩爾質(zhì)量和孢子數(shù)量分別占總量的80%、40%、20%、10%和5%作為不同的處理,進行AS-real-time PCR檢測并驗證標(biāo)準(zhǔn)曲線的可行性。結(jié)果顯示預(yù)期值為80%、40%、20%、10%和5%時,相應(yīng)的DNA檢測值依次為85.95%、44.64%、22.24%、11.86%和6.40%,孢子檢測值依次為116.45%、48.39%、26.38%、11.86%和6.83%(表5)。隨著預(yù)期值逐漸降低檢測值的標(biāo)準(zhǔn)差也依次下降,且預(yù)期值在5%時具有最小的標(biāo)準(zhǔn)差。同時預(yù)期值和檢測值具有很高的相關(guān)性(2=0.9998和2=0.9922),說明此檢測體系具有可行性,在檢測SdhB-H278R突變比例為5%時具有更高的準(zhǔn)確性。

    圖8 AS-real-time PCR建立檢測SdhB-H278R突變的標(biāo)準(zhǔn)曲線

    3 討論

    多主棒孢是一種重要的植物病原菌,是棒孢屬內(nèi)發(fā)現(xiàn)最早、寄主范圍最廣的種,可侵染葫蘆科、茄科和豆科等大部分蔬菜,園藝花卉,橡膠等經(jīng)濟作物以及大豆等糧食作物,約380屬530多種植物[20]。黃瓜棒孢葉斑病在我國北京、山東、甘肅、河北、廣東和黑龍江等多個省份[21-23]以及在日本和韓國[9,10,24]均大范圍發(fā)生,一般田間發(fā)病率為10%—25%,嚴重時可達60%—70%,甚至100%。因此由多主棒孢引起的黃瓜棒孢葉斑病在黃瓜上已成為一種重要的病害。啶酰菌胺是德國巴斯夫公司開發(fā)的一種SDHIs類殺菌劑,自2006年開始在日本登記防治黃瓜棒孢葉斑病,隨后很快出現(xiàn)了對啶酰菌胺的抗藥性[9-10]。但在國內(nèi)該藥劑并未登記在防治黃瓜棒孢葉斑病上,只登記用于防治黃瓜灰霉病,隨著啶酰菌胺在防治黃瓜灰霉病中的廣泛使用,灰葡萄孢對啶酰菌胺產(chǎn)生了抗藥性[25],間接地多主棒孢對啶酰菌胺也相應(yīng)產(chǎn)生了抗藥性,導(dǎo)致黃瓜棒孢葉斑病大量發(fā)生。目前多主棒孢上的抗性突變包括SdhB-P199S、H278Y/R、I280V、SdhC-S73P、SdhD-S89P、D95E、H108R、G109V[10],而H278R為首次在國內(nèi)報道,因此建立一種AS-real-time PCR檢測方法非常必要,可用于指導(dǎo)后續(xù)的抗性治理。若H278R抗性頻率升高,應(yīng)使該藥劑減少甚至停止使用,并采用不同作用機制的殺菌劑輪換使用或具有負交互抗藥性藥劑復(fù)配使用,可延長啶酰菌胺的使用壽命,減少黃瓜棒孢葉斑病的發(fā)生。

    表5 AS-real-time PCR驗證

    通過測定在北京市大興區(qū)采集的黃瓜多主棒孢對啶酰菌胺的敏感性并參考Miyamoto等[9]的研究結(jié)果,將多主棒孢菌株分為敏感和抗性菌株兩個群體,通過對黃瓜多主棒孢序列的測定,抗性菌株均為SdhB-H278R突變。Miyamoto等[9]檢測到SdhB-H278R共4株且劃分為高等抗性菌株,EC50范圍為8.9—10.7 μg·mL-1,而本文檢測到SdhB-H278R突變共8株且EC50范圍更廣(5.395—11.710 μg·mL-1),推測可能與檢測到H278R基因型的數(shù)量有關(guān)。通過對抗性和敏感菌株的生物學(xué)特性分析,發(fā)現(xiàn)抗性和敏感菌株除在菌絲生長速率方面存在顯著差異外,在產(chǎn)孢量和致病性方面差異不顯著,說明H278R突變株與敏感菌株之間的田間適合度差異不大,在藥劑選擇壓力下易發(fā)展成優(yōu)勢群體導(dǎo)致黃瓜葉斑病的大量發(fā)生。-H272Y/R/L、H277Y/R/Y和H278R/Y幾種點突變類型已經(jīng)在小雙胞腔菌、鏈格孢和灰葡萄孢等多種病原菌中被報道[7,14,16,26-28]。這幾種點突變均發(fā)生在泛醌結(jié)合口袋及其附近[29],且氨基酸序列具有高度的保守性,可以將這幾種突變歸為一類。LALèVE等[30]研究發(fā)現(xiàn),灰葡萄孢H272R/L/Y突變在以葡萄糖作為碳源時,前期菌絲生長速率較慢,后期一致;另外H272Y突變不影響SDH和呼吸活性。推測這類點突變在田間的適應(yīng)性較強,極易導(dǎo)致田間病害的發(fā)生。結(jié)合AS-real-time PCR檢測體系,預(yù)測含有H278R突變位點菌株的發(fā)展趨勢,可指導(dǎo)藥劑的合理使用及制定合理的田間管理策略。

    本研究建立的AS-real-time PCR檢測方法可以快速定量檢測田間多主棒孢SdhB-H278R突變位點的頻率,經(jīng)驗證AS-real-time PCR檢測方法預(yù)期值和試驗值的線性相關(guān)性為2=0.9998和2=0.9922。此方法可以檢測到5%的SdhB-H278R抗性突變,靈敏度是普通PCR的10倍。相比高葦?shù)萚31-32]建立的檢測黃瓜多主棒孢實時熒光定量PCR檢測技術(shù)和普通PCR檢測技術(shù),此檢測體系可以定量檢測SdhB-H278R抗性突變位點以及抗性突變在田間頻率的變化,克服了PCR-RFLP檢測體系中變異系數(shù)大的難題[33-35]。與常規(guī)PCR檢測方法相比較,該檢測體系具有快速、無污染和靈敏度高的特點,并可根據(jù)檢測結(jié)果指導(dǎo)科學(xué)用藥以及降低用藥成本和減少環(huán)境污染。

    4 結(jié)論

    建立了快速檢測黃瓜多主棒孢菌SdhB-H278R突變位點的AS-real-time PCR技術(shù)體系,該體系自動化程度高,準(zhǔn)確性強、靈敏度高,整個過程可在一個封閉的體系中進行,無污染、操作簡單,可用于檢測該突變位點在田間的頻率。

    [1] 李寶聚, 高葦, 石延霞, 謝學(xué)文. 多主棒孢和棒孢葉斑病的研究進展. 植物保護學(xué)報, 2012, 39(2): 171-176.LI B J, GAO W, SHI Y X, XIE X W. Progress in researches onleaf spot., 2012, 39(2): 171-176. (in chinese)

    [2] 于淑晶, 王滿意, 田芳, 趙衛(wèi)光, 邊強, 李寶聚. 黃瓜棒孢葉斑病的防治及抗藥性研究進展. 農(nóng)藥, 2014, 53(1): 7-11.YU S J, WANG M Y, TIAN F, ZHAO W G, BIAN Q, LI B J. Progress in research on control of cucumberleaf spot and fungicide resistance., 2014, 53(1): 7-11. (in chinese)

    [3] 李良孔, 袁善奎, 潘洪玉, 王巖. 琥珀酸脫氫酶抑制劑類(SDHIs)殺菌劑及其抗性研究進展. 農(nóng)藥, 2011, 50(3): 165-169.LI L K, YUAN S K, PAN H Y, WANG Y. Progress in research on SDHIs fungicides and its resistance., 2011, 50(3): 165-169. (in chinese)

    [4] 顏范勇, 劉冬青, 司馬利鋒, 石恒, 胡欣. 新型煙酰胺類殺菌劑——啶酰菌胺. 農(nóng)藥, 2008, 47(2): 132-135.YAN F Y, LIU D Q, SIMA L F, SHI H, HU X. Boscalid, a novel carboxamide aka anilide class of fungicides., 2008, 47(2): 132-135. (in chinese)

    [5] AVENOT H, Michailides T J. Resistance to boscalid fungicide inisolates from pistachio in California., 2007, 91(10): 1345-1350.

    [6] YIN Y N, KIM Y K, XIAO C L. Molecular characterization of boscalid resistance in field isolates offrom apple., 2011, 101(8): 986-995.

    [7] AVENOT H F, THOMAS A, GITAITIS R D, JR LANGSTON D B, STEVENSON K L. Molecular characterization of boscalid- and penthiopyrad-resistant isolates ofand assessment of their sensitivity to fluopyram., 2012, 68(4): 645-651.

    [8] WANG Y, DUAN Y B, WANG J X, ZHOU M G. A new point mutation in the iron-sulfur subunit of succinate dehydrogenase confers resistance to boscalid in., 2015, 16(7): 653-661.

    [9] MIYAMOTO T, ISHII H, SEKO T, KOBORI S, TOMITA Y. Occurrence ofisolates resistant to boscalid on cucumber in Ibaraki Prefecture, Japan., 2009, 58(6): 1144-1151.

    [10] MIYAMOTO T, ISHII H, STAMMLER G, KOCH A, OGAWARA T, TOMITA Y, FOUNTAINE J M, USHIO S, SEKO T, KOBORI S. Distribution and molecular characterization ofisolates resistant to boscalid., 2010, 59(5): 873-881.

    [11] FURUYA S, SUZUKI S, KOBAYASHI H, SAITO S, TAKAYANAGI T. Rapid method for detecting resistance to a QoI fungicide inpopulations., 2010, 65(8): 840-843.

    [12] AOKI Y, FURUYA S, SUZUKI S. Method for rapid detection of thegene allele conferring resistance to mandipropamid, a carboxylic acid amide fungicide, inpopulations., 2011, 67(12): 1557-1561.

    [13] 李紅霞, 周明國. 用等位基因特異性寡核苷酸(ASO)-PCR快速檢測抗多菌靈的油菜菌核病菌. 中國農(nóng)業(yè)科學(xué), 2004, 37(9): 1396-1399.LI H X, ZHOU M G. Rapid identification of carbendazim resistant strains ofusing allele-specific oligonucleotide (ASO)-PCR., 2004, 37(9): 1396-1399. (in chinese)

    [14] MALLIK I, ARABIAT S, PASCHE J S, BOLTON M D, PATEL J S, GUDMESTAD N C. Molecular characterization and detection of mutations associated with resistance to succinate dehydrogenase- inhibiting fungicides in., 2014, 104(1): 40-49.

    [15] WHEELER I, KENDALL S, BUTTERS J, HOLLOMON D. Detection of benzimidazole resistance inusing allele-specific oligonucleotide probes./, 1995, 25(1/2): 113-116.

    [16] DE MICCOLIS ANGELINI R M, MASIELLO M, ROTOLO C, POLLASTRO S, FARETRA F. Molecular characterization and detection of resistance to succinate dehydrogenase inhibitor fungicides in()., 2014, 70(12): 1884-1893.

    [17] LEHNER M S, JúNIOR T J P, SILVA R A, VIEIRA R F, SCHNABEL G, MIZUBUTI E S G. Fungicide sensitivity of: A thorough assessment using discriminatory dose, EC50, high-resolution melting analysis, and description of new point mutation associated with thiophanate-methyl resistance., 2015, 99(11): 1537-1543.

    [18] SAMARAS A, MADESIS P, KARAOGLANIDIS G S. Detection ofB gene mutations in SDHI-resistant isolates ofusing high resolution melting (HRM) analysis., 2016, 7: Article 1815.

    [19] 聶燕釵, 王斌, 趙子琴, 周懷谷. 等位基因特異性PCR技術(shù)及其法醫(yī)學(xué)應(yīng)用. 法醫(yī)學(xué)雜志, 2014, 30(4): 282-287.NIE Y C, WANG B, ZHAO Z Q, ZHOU H G. Allele-specific PCR and its application in forensic science., 2014, 30(4): 282-287. (in chinese)

    [20] DIXON L J, SCHLUB R L, PERNEZNY K, DATNOFF L E. Host specialization and phylogenetic diversity of, 2009, 99(9): 1015-1027.

    [21] 楊苗.我國蔬菜棒孢葉斑病病原菌多樣性研究[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2013.YANG M. Diversity of pathogen of corynespora leaf spot on vegetables in china[D]. Beijing: Chinese Academy of Agricultural Sciences, 2013. (in chinese)

    [22] 高葦, 李寶聚, 石延霞, 謝學(xué)文. 河北青縣黃瓜棒孢葉斑病病原菌種群分化的研究. 華北農(nóng)學(xué)報, 2011, 26(5): 9-15.GAO W, LI B J, SHI Y X, XIE X W. Population differentiation ofin qing county, Hebei province., 2011, 26(5): 9-15. (in chinese)

    [23] 李長松, 張眉, 李林, 李凡, 齊軍山, 徐作珽, 張博. 山東省黃瓜棒孢葉斑病(褐斑病)病原菌鑒定和防治. 中國蔬菜, 2009(18): 29-33.Li C S, ZHANG M, Li L, LI F, QI J S, XU Z T, ZHANG B. Identification of cucumber target leaf spot (brown spot) pathogen and its control., 2009(18): 29-33. (in chinese)

    [24] KWON M K, KANG B R, CHO B H, KIM Y C. Occurrence of target leaf spot disease caused byon cucumber in Korea., 2010, 52(3): 424.

    [25] 余玲, 劉慧平, 韓巨才, 張寶俊. 山西省灰霉菌對啶酰菌胺的敏感性測定. 山西農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2012, 32(3): 232-234.YU L, LIU H P, HAN J C, ZHANG B J. Sensitivity offrom Shanxi province to boscalid., 2012, 32(3): 232-234. (in chinese)

    [26] VELOUKAS T, LEROCH M, HAHN M, KARAOGLANIDIS G S. Detection and molecular characterization of boscalid-resistantisolates from strawberry., 2011, 95(10): 1302-1307.

    [27] 史曉晶. 番茄早疫病菌對啶酰菌胺的抗性檢測及抗性機理初探[D]. 太谷: 山西農(nóng)業(yè)大學(xué), 2015. SHI X J. Study on sensitivity ofto boscalid and resistant mechanisms[D]. Taigu: Shanxi Agricultural University, 2015. (in chinese)

    [28] YANG J H, BRANNEN P M, SCHNABEL G. Resistance into SDHI fungicides causes rare disease outbreak in peach orchards., 2015, 99(1): 65-70.

    [29] HORSEFIELD R, YANKOVSKAYA V, SEXTON G, WHITTINGHAM W, SHIOMI K, OMURA S, BYRNE B, CECCHINI G, IWATA S. Structural and computational analysis of the quinone-binding site of complex II (succinate-ubiquinone oxidoreductase): a mechanism of electron transfer and proton conduction during ubiquinone reduction., 2006, 281(11): 7309-7316.

    [30] LALèVE A, GAMET S, WALKER A, DEBIEU D, TOQUIN V, FILLINGER S. Site-directed mutagenesis of the P225, N230 and H272 residues of succinate dehydrogenase subunit B fromhighlights different roles in enzyme activity and inhibitor binding., 2014, 16(7): 2253-2266.

    [31] 高葦, 李寶聚, 王萬立, 郝永娟, 石延霞. 土壤中黃瓜棒孢葉斑病病原菌實時熒光定量PCR檢測技術(shù)研究. 華北農(nóng)學(xué)報, 2014, 29(2): 71-74.GAO W, LI B J, WANG W L, HAO Y J, SHI Y X. Detection ofin soil with real-time quantitative PCR., 2014, 29(2): 71-74. (in chinese)

    [32] 陳璐. 黃瓜細菌性角斑病菌和多主棒孢菌PCR檢測技術(shù)的建立[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2014.CHEN L. PCR-based specific detection ofpv.and[D]. Beijing: Chinese Academy of Agricultural Sciences, 2014. (in chinese)

    [33] MICHALECKA M, MALINOWSKI T, BRONIAREKNIEMIEC A, BIELENIN A. Real-time PCR assay with SNP-specific primers for the detection of a G143A mutation level infield populations., 2011, 159(7/8): 569-578.

    [34] LUO Y, MA Z H, MICHAILIDES T J. Quantification of allele E198A in beta-tubulin conferring benzimidazole resistance inusing real-time PCR., 2007, 63(12): 1178-1184.

    [35] KIANIANMOMENI A, SCHWARZ G, FELSENSTEIN F G, WENZEL G. Validation of a real-time PCR for the quantitative estimation of a G143A mutation in the cytochrome1gene of., 2010, 63(3): 219-224.

    Establishment of AS-real-time PCR for Quantitatively Detecting the H278R Allele in theAssociated within Cucumber

    SUN BingXue, SHI YanXia, ZHU FaDI, XIE XueWen, CHAI ALi, LI BaoJu

    (Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081)

    【Objective】The objective of this study is to establish a rapid and efficient allele specific real-time PCR (AS-real-time PCR) method for quantitative detection of the H278R mutation in theassociated within cucumber, and to verify the effect. 【Method】A total of 24 single conidial strains ofwere isolated from Daxing, Beijing.The method of mycelial growth inhibition was used to determine the EC50value to boscalid. Mycelium growth rate, spore outputs and pathogenicity of 4 susceptible (S) and 8 resistant (R) strains were measured. Primer pair Cc-SdhB-F/R was used to detect the sequence ofand the base change ofBased on the sequencing results inby Cc-SdhB-F/R primers, specific primer pair B-H278R-2F/2R14 and internal primers B-H278R-TY-F/R were designed. An AS-real-time PCR reaction system was established and optimized. The specificity, sensitivity and melt curve of the system were also evaluated.DNA and spore suspension containing different proportions of H278R mutation were detected with the optimized reaction system, respectively. 【Result】The sensitive frequency ofto boscalid was 66.67% in 24 strains and the EC50values ranged from 0.057 to 0.563 μg·mL-1. The resistance frequency ofto boscalid was 33.33% and the EC50values ranged from 5.395 to 11.710 μg·mL-1. Resistant and susceptible strains only significantly differed in terms of mycelium growth rate, and there was a significant negative correlation between the rate of mycelial growth and EC50values, but no significant difference in spore outputs and pathogenicity. Sequencing analysis showed that all the resistant strains carried SdhB-H278R mutation. The primers B-H278R-2F/2R14 were of great specificity, the specific fragment was amplified from the DNA of H278R strains, but not from the DNA of other fungal strains. The sensitivity of ordinary AS-PCR was 91 pg·μL-1, while that of AS-real-time PCR was 9.1 pg·μL-1, which was 10 times higher than that of ordinary AS-PCR. The standard curve established by AS-real-time PCR showed a fine linear relationship between ΔCT value and lg of H278R frequency, the correlation coefficient of the standard curve was 0.9857 and with high amplification efficiency (92.59%). The absorption peak of melting curve was single. The internal reference primers B-H278R-TY-F/R and the specific primer B-H278R-2F/2R14 had a single specific peak at 87.81℃and 91.62℃, respectively. The result of mixtures DNA and mixtures spore suspension to verify the standard curve showed a fine linear relationship between expected percentage and detected percentage (2=0.9998 and2=0.9922). As the proportions of H278R mutation inDNA and spore suspension decreased, the accuracy of the system gradually increased, and the detection limit was 5%. 【Conclusion】An efficient, quantitative, and rapid AS-real-time PCR detection system was established for the detection of SdhB-H278R mutation site, which can provide a theoretical basis for the SDHIs resistance management.

    cucumber;; AS-real-time PCR; resistance;; boscalid

    2018-06-22;

    2018-09-15

    國家重點研發(fā)計劃(2018YFD0201300)、中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程(CAAS-ASTIP-IVFCAAS)、農(nóng)業(yè)部園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點實驗室開放課題(IVF2017)

    孫炳學(xué),E-mail:sunbingxuechina@163.com。石延霞,E-mail:shiyanxia@caas.cn。孫炳學(xué)和石延霞為同等貢獻作者。

    李寶聚,E-mail:libaoju@caas.cn

    10.3864/j.issn.0578-1752.2018.24.006

    (責(zé)任編輯 岳梅)

    猜你喜歡
    檢測
    QC 檢測
    “不等式”檢測題
    “一元一次不等式”檢測題
    “一元一次不等式組”檢測題
    “幾何圖形”檢測題
    “角”檢測題
    “有理數(shù)的乘除法”檢測題
    “有理數(shù)”檢測題
    “角”檢測題
    “幾何圖形”檢測題
    亚洲一卡2卡3卡4卡5卡精品中文| 亚洲全国av大片| 日韩精品免费视频一区二区三区| 久久精品综合一区二区三区| 免费在线观看黄色视频的| 久久伊人香网站| 欧美精品啪啪一区二区三区| 久久久久亚洲av毛片大全| 亚洲av成人av| 给我免费播放毛片高清在线观看| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区免费观看 | 国产久久久一区二区三区| 一本大道久久a久久精品| 日韩欧美 国产精品| 国产91精品成人一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 哪里可以看免费的av片| 1024手机看黄色片| 窝窝影院91人妻| 亚洲熟女毛片儿| 国内揄拍国产精品人妻在线| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区| 国产亚洲精品综合一区在线观看 | 狂野欧美激情性xxxx| 欧美日本视频| 国产在线观看jvid| 大型av网站在线播放| 久久久国产成人精品二区| 国产亚洲精品综合一区在线观看 | 欧美日韩亚洲国产一区二区在线观看| 手机成人av网站| 精华霜和精华液先用哪个| 听说在线观看完整版免费高清| 国产精品美女特级片免费视频播放器 | 亚洲国产欧美网| 日本免费一区二区三区高清不卡| 国产av又大| 免费看日本二区| 亚洲一区二区三区不卡视频| 美女大奶头视频| 亚洲性夜色夜夜综合| 美女 人体艺术 gogo| 国产av麻豆久久久久久久| 三级国产精品欧美在线观看 | 亚洲成av人片免费观看| 国产成人一区二区三区免费视频网站| 精品日产1卡2卡| 成人18禁在线播放| 久久久国产成人精品二区| 成人av在线播放网站| 给我免费播放毛片高清在线观看| 亚洲成a人片在线一区二区| 欧美国产日韩亚洲一区| 日日爽夜夜爽网站| 大型黄色视频在线免费观看| 久久久国产精品麻豆| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| 亚洲美女视频黄频| 国模一区二区三区四区视频 | 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 国产精品精品国产色婷婷| 99久久精品国产亚洲精品| 99久久99久久久精品蜜桃| 欧美+亚洲+日韩+国产| АⅤ资源中文在线天堂| 99久久国产精品久久久| 欧美乱色亚洲激情| 久久精品aⅴ一区二区三区四区| 夜夜爽天天搞| 国产成人啪精品午夜网站| 老熟妇乱子伦视频在线观看| 黄色女人牲交| 香蕉国产在线看| 欧美最黄视频在线播放免费| 久久婷婷成人综合色麻豆| 变态另类丝袜制服| 亚洲色图 男人天堂 中文字幕| 国产成人av激情在线播放| 国产一级毛片七仙女欲春2| 国产成年人精品一区二区| 精品欧美国产一区二区三| 日韩大尺度精品在线看网址| 国产精品影院久久| 18禁美女被吸乳视频| 亚洲av电影不卡..在线观看| 在线观看www视频免费| 亚洲熟妇中文字幕五十中出| 亚洲欧美激情综合另类| 老司机在亚洲福利影院| 一二三四社区在线视频社区8| 欧美日韩福利视频一区二区| 99国产精品一区二区三区| 国产又色又爽无遮挡免费看| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 国产高清视频在线播放一区| a级毛片a级免费在线| 国产日本99.免费观看| 一本综合久久免费| 免费av毛片视频| 18禁裸乳无遮挡免费网站照片| 国产在线精品亚洲第一网站| 欧美又色又爽又黄视频| 母亲3免费完整高清在线观看| 亚洲第一电影网av| 久久久久久人人人人人| 黄色视频,在线免费观看| 性欧美人与动物交配| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 熟女电影av网| 久久久精品欧美日韩精品| 欧美中文日本在线观看视频| 久久香蕉精品热| 波多野结衣高清作品| 精品熟女少妇八av免费久了| 国产伦在线观看视频一区| av欧美777| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 午夜视频精品福利| 美女大奶头视频| 女生性感内裤真人,穿戴方法视频| 叶爱在线成人免费视频播放| 欧美乱妇无乱码| 在线播放国产精品三级| 99热这里只有是精品50| 一区二区三区国产精品乱码| 白带黄色成豆腐渣| 亚洲 国产 在线| 亚洲,欧美精品.| 黄频高清免费视频| 女警被强在线播放| 欧美大码av| 国语自产精品视频在线第100页| 精品人妻1区二区| 国产欧美日韩精品亚洲av| 母亲3免费完整高清在线观看| 精品日产1卡2卡| 中文资源天堂在线| 日韩精品青青久久久久久| 99在线人妻在线中文字幕| 欧美乱妇无乱码| 国产亚洲精品综合一区在线观看 | 最近最新中文字幕大全电影3| 黄色 视频免费看| 中文字幕熟女人妻在线| 中文在线观看免费www的网站 | 亚洲第一欧美日韩一区二区三区| 动漫黄色视频在线观看| 舔av片在线| www.精华液| 老汉色∧v一级毛片| 一个人免费在线观看的高清视频| 欧美一级毛片孕妇| 亚洲电影在线观看av| 国产一区二区在线av高清观看| 女人被狂操c到高潮| 三级男女做爰猛烈吃奶摸视频| 久久精品影院6| 亚洲成a人片在线一区二区| 国产精品99久久99久久久不卡| 老汉色av国产亚洲站长工具| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久久电影 | 深夜精品福利| 成年女人毛片免费观看观看9| 亚洲色图av天堂| 啦啦啦免费观看视频1| 久久这里只有精品中国| 国产精品电影一区二区三区| 成人特级黄色片久久久久久久| 在线免费观看的www视频| 制服诱惑二区| 午夜久久久久精精品| 亚洲专区字幕在线| 国产av在哪里看| 亚洲精品色激情综合| 亚洲五月天丁香| 高潮久久久久久久久久久不卡| 久久久国产欧美日韩av| 亚洲欧美日韩高清专用| www日本在线高清视频| 久久中文看片网| 久久久久性生活片| 国产伦人伦偷精品视频| 国产熟女午夜一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 国产主播在线观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 757午夜福利合集在线观看| 欧美黄色片欧美黄色片| 成人三级做爰电影| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 正在播放国产对白刺激| 又粗又爽又猛毛片免费看| а√天堂www在线а√下载| 成人18禁高潮啪啪吃奶动态图| 波多野结衣高清无吗| 欧美黑人巨大hd| 免费观看精品视频网站| 香蕉av资源在线| aaaaa片日本免费| 国产精品98久久久久久宅男小说| 制服丝袜大香蕉在线| 制服人妻中文乱码| 欧美在线黄色| xxxwww97欧美| 精品国产亚洲在线| 欧美成人性av电影在线观看| 久久精品国产综合久久久| 香蕉av资源在线| 最新在线观看一区二区三区| 国产精品综合久久久久久久免费| 久久久久久久精品吃奶| 国产熟女xx| 久久天堂一区二区三区四区| 亚洲成a人片在线一区二区| 中文亚洲av片在线观看爽| 露出奶头的视频| 欧美国产日韩亚洲一区| 精华霜和精华液先用哪个| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品免费视频内射| 久久精品国产亚洲av高清一级| 淫妇啪啪啪对白视频| 97人妻精品一区二区三区麻豆| 毛片女人毛片| 精品福利观看| 国产伦人伦偷精品视频| 99久久99久久久精品蜜桃| 美女黄网站色视频| 亚洲国产精品久久男人天堂| 18禁黄网站禁片免费观看直播| 久久香蕉激情| 日本黄色视频三级网站网址| 最近最新中文字幕大全免费视频| 国产午夜精品论理片| 宅男免费午夜| √禁漫天堂资源中文www| 精品人妻1区二区| 成人手机av| 亚洲va日本ⅴa欧美va伊人久久| 国产午夜精品久久久久久| 99热这里只有精品一区 | 一区二区三区激情视频| 免费人成视频x8x8入口观看| 香蕉久久夜色| 舔av片在线| 久久久久久久久免费视频了| 久久婷婷人人爽人人干人人爱| 亚洲成a人片在线一区二区| 亚洲一区二区三区色噜噜| 岛国在线免费视频观看| 女人爽到高潮嗷嗷叫在线视频| www.999成人在线观看| 色av中文字幕| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 手机成人av网站| 亚洲性夜色夜夜综合| 中文字幕久久专区| 亚洲精品色激情综合| 亚洲精品国产一区二区精华液| 成人国产一区最新在线观看| 欧美日韩一级在线毛片| 久久婷婷成人综合色麻豆| 日本 欧美在线| 国产精品亚洲美女久久久| 亚洲精品粉嫩美女一区| 黑人欧美特级aaaaaa片| 免费av毛片视频| 99久久综合精品五月天人人| 午夜福利18| 亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品99久久久久| 757午夜福利合集在线观看| 丁香六月欧美| 国产精品1区2区在线观看.| 亚洲欧美日韩无卡精品| 亚洲九九香蕉| 蜜桃久久精品国产亚洲av| 岛国在线观看网站| 亚洲精品久久成人aⅴ小说| 久久精品91无色码中文字幕| 夜夜躁狠狠躁天天躁| 国产一区在线观看成人免费| 可以免费在线观看a视频的电影网站| 好男人电影高清在线观看| 免费看十八禁软件| 久热爱精品视频在线9| 亚洲精品中文字幕一二三四区| 在线国产一区二区在线| 91成年电影在线观看| 91老司机精品| 在线免费观看的www视频| 亚洲男人的天堂狠狠| av中文乱码字幕在线| 可以在线观看的亚洲视频| 999久久久精品免费观看国产| 久久精品91无色码中文字幕| 一级黄色大片毛片| 又大又爽又粗| 女同久久另类99精品国产91| 日韩av在线大香蕉| 色老头精品视频在线观看| 国产精品av久久久久免费| 中文字幕高清在线视频| 国产一区二区激情短视频| 亚洲人成网站高清观看| 免费在线观看完整版高清| 国产v大片淫在线免费观看| 成人国产综合亚洲| 美女黄网站色视频| 国内精品久久久久久久电影| 怎么达到女性高潮| 最好的美女福利视频网| 两个人视频免费观看高清| 国产片内射在线| 国产又色又爽无遮挡免费看| 中出人妻视频一区二区| 日韩欧美精品v在线| 很黄的视频免费| 欧美+亚洲+日韩+国产| 久久香蕉激情| 好男人电影高清在线观看| 99精品久久久久人妻精品| 日本黄大片高清| 亚洲五月婷婷丁香| 亚洲熟女毛片儿| 欧美成人一区二区免费高清观看 | 国产精品久久久久久久电影 | 18禁美女被吸乳视频| 久久久久免费精品人妻一区二区| 最近视频中文字幕2019在线8| 久久 成人 亚洲| 免费高清视频大片| 一级毛片女人18水好多| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久精免费| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 99精品欧美一区二区三区四区| 国产片内射在线| 少妇被粗大的猛进出69影院| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| 午夜成年电影在线免费观看| 久久99热这里只有精品18| 特级一级黄色大片| 久久久精品欧美日韩精品| 亚洲国产欧美网| 欧美一区二区国产精品久久精品 | 香蕉av资源在线| 精品久久久久久久久久久久久| 亚洲片人在线观看| 日韩欧美在线二视频| 国产高清激情床上av| 99热这里只有精品一区 | 长腿黑丝高跟| 99国产精品99久久久久| 欧美性猛交╳xxx乱大交人| 韩国av一区二区三区四区| av在线天堂中文字幕| 91在线观看av| 久久精品国产亚洲av高清一级| 欧美日韩瑟瑟在线播放| 久久久久久久精品吃奶| 午夜激情av网站| 69av精品久久久久久| 757午夜福利合集在线观看| 亚洲一码二码三码区别大吗| 欧美日韩精品网址| 这个男人来自地球电影免费观看| 亚洲电影在线观看av| 香蕉国产在线看| 悠悠久久av| 国产成人aa在线观看| 18禁观看日本| 在线看三级毛片| 村上凉子中文字幕在线| 亚洲avbb在线观看| 亚洲成人国产一区在线观看| 中文字幕久久专区| 亚洲人成77777在线视频| 亚洲国产日韩欧美精品在线观看 | 国产精品日韩av在线免费观看| 一级毛片高清免费大全| 制服丝袜大香蕉在线| 又黄又爽又免费观看的视频| 又大又爽又粗| 少妇粗大呻吟视频| aaaaa片日本免费| 村上凉子中文字幕在线| 99热这里只有精品一区 | 久久久久久大精品| 色综合婷婷激情| 欧美日韩亚洲综合一区二区三区_| 亚洲自偷自拍图片 自拍| 亚洲一码二码三码区别大吗| 国产精品亚洲一级av第二区| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻1区二区| 欧美中文日本在线观看视频| 国产又色又爽无遮挡免费看| 日韩欧美国产在线观看| 99在线人妻在线中文字幕| 波多野结衣巨乳人妻| av在线天堂中文字幕| 99国产精品一区二区三区| 久久亚洲真实| 国产一区二区三区视频了| 操出白浆在线播放| 免费看a级黄色片| 精品欧美一区二区三区在线| 久久亚洲精品不卡| 久久久久久久久免费视频了| 90打野战视频偷拍视频| 亚洲专区字幕在线| 精品久久久久久久末码| 亚洲成人免费电影在线观看| 女警被强在线播放| 搡老岳熟女国产| 日本在线视频免费播放| 免费高清视频大片| 搡老熟女国产l中国老女人| а√天堂www在线а√下载| 黄色成人免费大全| 欧美人与性动交α欧美精品济南到| www.自偷自拍.com| 99久久综合精品五月天人人| 宅男免费午夜| 99国产精品99久久久久| 亚洲精品中文字幕一二三四区| 婷婷精品国产亚洲av| 给我免费播放毛片高清在线观看| 日韩精品青青久久久久久| 少妇的丰满在线观看| 国产av一区二区精品久久| 欧美久久黑人一区二区| 国产精品九九99| 亚洲av片天天在线观看| 国产精品久久久久久人妻精品电影| 在线国产一区二区在线| 日本免费a在线| av免费在线观看网站| 久久午夜亚洲精品久久| 国产av在哪里看| 亚洲avbb在线观看| 天堂动漫精品| 国产伦在线观看视频一区| 老鸭窝网址在线观看| 国产99白浆流出| 亚洲人成伊人成综合网2020| 日韩精品免费视频一区二区三区| 午夜精品一区二区三区免费看| 露出奶头的视频| 久久久久国产精品人妻aⅴ院| 中亚洲国语对白在线视频| 身体一侧抽搐| 又爽又黄无遮挡网站| 国产精品久久久人人做人人爽| 91老司机精品| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 国产熟女午夜一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 最新美女视频免费是黄的| 国产精品日韩av在线免费观看| 欧美绝顶高潮抽搐喷水| 叶爱在线成人免费视频播放| 韩国av一区二区三区四区| 日韩三级视频一区二区三区| 制服人妻中文乱码| 国产蜜桃级精品一区二区三区| 欧美一级毛片孕妇| av国产免费在线观看| 少妇的丰满在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 一级作爱视频免费观看| 日本黄色视频三级网站网址| 久久久久久久精品吃奶| 国产午夜精品论理片| 久久久精品欧美日韩精品| 特级一级黄色大片| 男插女下体视频免费在线播放| 国产三级中文精品| 国产精品 欧美亚洲| 国产av不卡久久| 久久精品亚洲精品国产色婷小说| 99riav亚洲国产免费| 男人舔奶头视频| www.精华液| 变态另类丝袜制服| 色老头精品视频在线观看| 制服诱惑二区| 亚洲成人中文字幕在线播放| 精品国产乱子伦一区二区三区| 丝袜美腿诱惑在线| 亚洲中文字幕日韩| 国内精品久久久久久久电影| 色综合亚洲欧美另类图片| 一级毛片高清免费大全| 亚洲国产欧美一区二区综合| 老熟妇仑乱视频hdxx| 特大巨黑吊av在线直播| 天堂√8在线中文| 国产精品乱码一区二三区的特点| 一本一本综合久久| 精华霜和精华液先用哪个| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久5区| 亚洲 欧美 日韩 在线 免费| 欧美丝袜亚洲另类 | 久久久久国内视频| 午夜免费成人在线视频| 91大片在线观看| 黄色视频,在线免费观看| 日本成人三级电影网站| 国产日本99.免费观看| 精品国产美女av久久久久小说| 久久久久久久久中文| 中亚洲国语对白在线视频| 国产三级中文精品| 级片在线观看| 国产精品av久久久久免费| 全区人妻精品视频| 精品无人区乱码1区二区| 色综合亚洲欧美另类图片| 亚洲av成人一区二区三| 国产精品久久视频播放| 午夜免费成人在线视频| 最近最新中文字幕大全电影3| 久久中文字幕人妻熟女| 亚洲熟女毛片儿| 久久久久久亚洲精品国产蜜桃av| 88av欧美| 成人高潮视频无遮挡免费网站| 最近在线观看免费完整版| 美女午夜性视频免费| 日日夜夜操网爽| 精品免费久久久久久久清纯| 悠悠久久av| 亚洲精品国产一区二区精华液| 久久香蕉国产精品| 亚洲专区国产一区二区| 免费在线观看亚洲国产| 黄片小视频在线播放| 久久精品亚洲精品国产色婷小说| 精品欧美国产一区二区三| 99re在线观看精品视频| 黄色毛片三级朝国网站| 1024视频免费在线观看| 国产av在哪里看| 国产精品一区二区免费欧美| 狠狠狠狠99中文字幕| ponron亚洲| 亚洲国产精品久久男人天堂| 舔av片在线| 精品国产超薄肉色丝袜足j| 国产激情欧美一区二区| 日本 欧美在线| 亚洲真实伦在线观看| 久久午夜综合久久蜜桃| 一夜夜www| 国内少妇人妻偷人精品xxx网站 | 全区人妻精品视频| 免费看日本二区| 亚洲av美国av| 又黄又粗又硬又大视频| 久9热在线精品视频| 国产高清激情床上av| 欧美日本视频| 精品欧美国产一区二区三| av欧美777| 黑人操中国人逼视频| 亚洲成av人片免费观看| 午夜a级毛片| netflix在线观看网站| 国产亚洲av嫩草精品影院| 午夜福利欧美成人| 日韩欧美在线乱码| 成年人黄色毛片网站| 在线观看一区二区三区| 中文字幕高清在线视频| 青草久久国产| 五月伊人婷婷丁香| 国产成人精品久久二区二区免费| 亚洲精品久久国产高清桃花| 美女大奶头视频| 一区二区三区国产精品乱码| 男人的好看免费观看在线视频 | 午夜精品一区二区三区免费看| 亚洲黑人精品在线| 波多野结衣高清无吗| 淫妇啪啪啪对白视频| 母亲3免费完整高清在线观看| 2021天堂中文幕一二区在线观| 日韩大尺度精品在线看网址| 精品久久久久久久末码| 亚洲国产精品sss在线观看| 麻豆一二三区av精品| 久久性视频一级片| 一区福利在线观看| 亚洲成av人片在线播放无| 亚洲色图 男人天堂 中文字幕| 草草在线视频免费看| 女人爽到高潮嗷嗷叫在线视频| 一级作爱视频免费观看|