劉穎，唐元滿，張淑婷，郭兵，谷紀濤，丁偉,2

1 西南大學植物保護學院，重慶市北碚區(qū)天生路2號 400715；2 重慶煙草科學研究所，重慶市北碚區(qū)天生路2號 400715

生物技術

重慶煙區(qū)青枯雷爾氏菌生化變種及序列變種的特性研究

劉穎1，唐元滿1，張淑婷1，郭兵1，谷紀濤1，丁偉1,2

1 西南大學植物保護學院，重慶市北碚區(qū)天生路2號 400715；2 重慶煙草科學研究所，重慶市北碚區(qū)天生路2號 400715

為進一步明確重慶煙區(qū)煙草青枯菌的遺傳及生化特性，采用生化變種分類標準及演化型復合PCR、系統(tǒng)發(fā)育進化分析對重慶10個煙草青枯病發(fā)生區(qū)縣青枯雷爾氏菌的種下分化情況進行了研究。結果表明：46株煙草青枯雷爾氏菌均屬于生化變種3和亞洲分支菌株（演化型Ⅰ）。內切葡聚糖酶基因（egl）和DNA修復蛋白基因（mutS）部分序列的系統(tǒng)發(fā)育學分析表明：46株供試菌株可聚類到青枯雷爾氏菌演化型Ⅰ的3個序列變種，分別為序列變種15、17和44，其中序列變種17和44為優(yōu)勢菌株。表明重慶煙草青枯雷爾氏菌具有一定的遺傳分化和地理區(qū)域特性。

煙草青枯菌；生化變種；演化型；系統(tǒng)進化分析；內切葡聚糖酶基因；DNA修復蛋白基因

煙草青枯病是由青枯雷爾氏菌（Ralstonia solanacearum；以下簡稱青枯菌）引起的一種細菌性病害[1]。青枯菌廣泛分布于熱帶、亞熱帶及溫帶地區(qū)，能侵染54個科450余種植物[2]。煙草青枯病在我國許多煙區(qū)，如福建、湖南、云南、貴州、四川、湖北、重慶等地普遍發(fā)生[3-4]。近年來，由于環(huán)境、氣候、種植結構調整等方面的原因，該病的危害范圍有逐漸向北和冷涼地區(qū)擴展蔓延之勢[4]。

青枯菌因其寄主范圍、地理分布、致病性以及生理特性的多樣性和復雜性而被公認為是多變的復合種（species complex）[5]。根據(jù)傳統(tǒng)的分類框架青枯菌可分為不同的生理小種（race）和生化變種（biovar），生理小種是依據(jù)不同菌株對不同植物種類的致病性差異來劃分的，而生化變種則是根據(jù)其利用3種雙糖和3種己醇的能力差異來分類的[6-9]。隨著分子生物學技術的廣泛應用，F(xiàn)egan和Prior（2005年）共同提出了以演化型分類框架來描述青枯菌種以下的差異，該分類框架依次將青枯菌劃分為種(Species)、演化型(Phylotype)、序列變種(Sequevar)以及克隆(Clone)4種不同水平的分類單元，并分別建立了與之相對應的鑒定方法[10]。

青枯菌的種水平鑒定采用以lpxC基因為靶標的引物對759/760；演化型(Phylotype)為亞種水平，分為演化型Ⅰ、Ⅱ、Ⅲ、Ⅳ，可由基于ITS區(qū)域的演化型特異性復合PCR(Phylotype speci fi c multiplex PCR，Pmx-PCR)鑒定；每個演化型則又是由多個序列變種(Sequevar)組成的，一個序列變種可作為一組在特定基因范圍內序列高度保守的菌株集合體，目前用于界定序列變種的基因包括內切葡聚糖酶基因(Endoglucanase,egl)、Ⅲ型分泌系統(tǒng)調節(jié)基因(hrpB)和一個編碼錯配糾正蛋白的基因(mutS)等[11]。演化型分類框架與傳統(tǒng)的生理小種及生化變種分類方法相比可以更精確地反映青枯菌這一復合種的地理起源及種內的遺傳多樣性。

重慶是我國主要的植煙區(qū)，煙草種植分布在武陵山區(qū)和大巴山區(qū)的黔江、石柱、彭水、酉陽、武隆、南川、涪陵、豐都、萬州、巫溪、巫山、奉節(jié)等12個區(qū)縣，除巫山和奉節(jié)外，其余10個區(qū)縣均有青枯病發(fā)生。本研究分別采用生化變種和演化型分類框架對重慶煙草青枯菌群體進行種下分類研究，旨在了解重慶煙區(qū)煙草青枯菌的遺傳及生化特性，解釋該地區(qū)的青枯菌生態(tài)多樣性，揭示煙草青枯菌在地理環(huán)境、生態(tài)條件以及種植結構影響下的多樣性分化情況，為有針對性的控制病害提供參考依據(jù)。

1 材料與方法

1.1 青枯菌株的分離及DNA的提取

供試煙草青枯病菌株由西南大學天然產物農藥研究實驗室2013年以及2015年在重慶地區(qū)采集分離，分別采自黔江、石柱、彭水、酉陽、武隆、南川、涪陵、豐都、萬州、巫溪，共46株，供試菌株信息見表1。應用細菌基因組DNA提取試劑盒（TIANGEN,DP302）提取菌株DNA，-20℃保存?zhèn)溆谩?/p>

1.2 生化變種鑒定

供試菌株生化變種測定參照Hayward的方法[7]?；九囵B(yǎng)基包括：NH4H2PO4， 1.0 g·L-1； KCl， 0.2 g·L-1；MgSO4·7H2O， 0.2 g·L-1；蛋白胨， 1.0 g·L-1；溴百里酚藍， 0.03 g·L-1；瓊脂， 1.5 g·L-1（最終用1M NaOH調節(jié)pH值到7.1）。分別將3種雙糖（乳糖、麥芽糖和纖維二糖）和3種己醇（甘露醇、山梨醇和甜醇）配成10%的溶液，滅菌（3種雙糖用過濾法滅菌，3種己醇經121℃蒸氣滅菌）后分別加入基本培養(yǎng)基中，使其終濃度為1%。按實驗設計將含不同碳水化合物的培養(yǎng)基垂直分布于96孔板上前6孔，不含碳水化合物的培養(yǎng)基作為對照放于垂直的第7、8孔[12]。取200 μL培養(yǎng)基注入每個孔中，每孔接種青枯菌菌懸液5 μL（菌懸液濃度約2 × 109CFU·mL-1），第8孔作為對照不接種，置30℃的恒溫培養(yǎng)箱中培養(yǎng)21 d。根據(jù)供試菌株對3種雙糖和3種己醇的利用情況，即培養(yǎng)基的顏色變化，確定各菌株的生化變種分類歸屬。

表1 供試煙草菌株相關信息Tab.1 Tobacco Ralstonia solanacearum strains used in the study

續(xù)表1

1.3 演化型鑒定

根據(jù)演化型分類框架所設計的復合PCR引物對青枯菌菌株進行了演化型鑒定[9]。PCR擴增采用25 μL反應體系：2×Taq PCR MasterMix（TIANGEN,KT201）、 引 物 Nmult:21:1F，Nmult:21:2F和 Nmult:22:InF各 6 pmoles，Nmult:23:AF和Nmult:22:RR 各 18 pmoles，759/760各 4 pmoles， 青枯菌DNA，ddH2O。反應程序為96℃預變性5 min；94℃變性15 s，59℃退火30 s和72℃延伸30 s，30個循環(huán)；最后72℃延伸10 min ，4℃保存。取5 μL PCR產物于1%的瓊脂糖凝膠中電泳，電壓4 V/cm，通過Biorad凝膠成像儀觀察結果。實驗所用引物均由華大基因科技股份有限公司合成，引物序列見表2。

表2 用于遺傳多樣性分析的引物序列Tab.2 The primers used for genetic diversity analysis

1.4 部分egl基因擴增

用引物Endo-F和Endo-R（表2）擴增內源葡聚糖酶基因（egl）內部的750 bp片段。egl基因擴增PCR反應體系（25 μl）：2×Taq PCR MasterMix（TIANGEN, KT201），青枯菌DNA 50 ng，引物Endo-F和Endo-R各10 pmoles，ddH2O補足。egl擴增的PCR反應程序為96℃預變性9 min；95℃變性1min，64℃退火1 min，72℃延伸2 min，30個循環(huán)；最后72℃延伸10 min ，4℃保存。將PCR產物送華大基因科技股份有限公司測序。

1.5 部分mutS基因擴增

用引物mutS-RsF.1570和mutS-RsR.1926（表2）擴增DNA蛋白修復基因mutS，其擴增PCR反應體系（25 μL）：2×Taq PCR MasterMix （TIANGEN,KT201），青枯菌DNA 50 ng，引物mutS-RsF.1570和 mutS-RsR.1926各10 pmoles，ddH2O補 足。PCR反應程序為96℃預變性5 min；94℃變性1 min，66℃退火1 min，72℃延伸90 s，35個循環(huán)；最后72℃延伸5 min，4℃保存。將PCR產物送華大基因科技股份有限公司測序。

1.6 序列分析

將所測得的供試青枯菌株egl和mutS基因序列與GenBank上已登錄的青枯菌核酸序列進行比對（參考序列信息見表3），構建系統(tǒng)發(fā)育樹。698 bp的egl基因和682 bp的mutS基因的測序結果用Clustal x軟件比對后，再用MEGA 5軟件進行系統(tǒng)發(fā)育分析，采用Jukes and Cantor模型鄰接法(neighbor-joining，NJ)構建系統(tǒng)發(fā)育樹，1000次重復取樣作出Bootstrap樹狀圖。序列之間的遺傳距離用Kimura-2-parameter方法計算。并將所有的序列提交到GenBank數(shù)據(jù)庫中。

表3 參考序列信息Tab.3 Reference strains used for sequence analysis

續(xù)表3

2 結果與分析

2.1 生化變種鑒定結果

測定了46株供試菌株對乳糖、麥芽糖、纖維二糖和甘露醇、山梨醇、甜醇的利用能力（圖1）。結果表明，所有的供試菌株均能利用3種雙糖和3種己醇，參照生化變種劃分標準，重慶地區(qū)的煙草青枯菌均屬于生化變種3（表1）。

圖1 部分菌株生化變種結果展示Fig.1 Results of biovar test of R. solanacearum

2.2 Pmx-PCR的演化型鑒定結果

供試菌株演化型特異性復合PCR結果表明，所有菌株均可同時擴增得到144 bp和281 bp的2條特異性片段，其中281 bp的片段為青枯菌種特異性擴增條帶，144 bp的片段為演化型Ⅰ的特異性擴增條帶。從而在種和演化型分類單元水平上，揭示出供試菌株屬于青枯菌演化型Ⅰ型，即亞洲分支菌株。

2.3 egl基因分析

對供試菌株的egl基因部分序列（698 bp）進行系統(tǒng)發(fā)育分析，并從GenBank 上選取32個菌株的egl基因的序列（表3）來確定供試菌株在系統(tǒng)發(fā)育結構中的位置。由系統(tǒng)發(fā)育分析知，青枯菌被分為4個分支，即4個演化型分支，而我國煙草青枯菌均被聚類到演化型Ⅰ型分支中（圖2），這與Pmx-PCR的鑒定結果一致。

供試46株菌株被聚類在演化型Ⅰ型之下的3個分支中，它們分別與序列變種15、17和44對應（圖2）。其中，序列變種15的參考菌株為中國番茄菌株PSS358，僅包含了2株菌，分別是酉陽的CQYY1和萬州的CQWZ2。序列變種17以分離自中國花生的P11為參考菌株，包含18株供試菌，分別為彭水的CQPS1、CQPS2、CQPS3、CQPS4，酉陽的 CQYY2， 武 隆 的 CQWL2、CQWL3、CQWL4、CQWL6，豐都的CQFD1、CQFD2，黔江的CQQJ4、CQQJ9、CQQJ10、CQQJ11， 石 柱 的 CQSZ-3、CQSZ-4，以及萬州的CQWZ1。而序列變種44包含了2個小分支，1個分支以中國煙草青枯菌Tb28和Tb43為參考菌株，共9株菌，其中酉陽2株（CQYY5、CQYY6）、黔江 7株（CQQJ1、CQQJ2、CQQJ3、CQQJ6、CQQJ7、CQQJ8、CQQJ12）；1個分支以中國木槿菌株Bd11為參考菌株，共有17株，其中彭水2株（CQPS5、CQPS6）、酉陽2株（CQYY3、CQYY4）、 武 隆 3株（CQWL1、CQWL5、CQWL7）、 南 川 2株（CQNC2、CQNC2）、 涪陵 2株（CQFL1、CQFL2）、黔江3株（CQQJ5、CQQJ13、CQQJ14）、石柱2株（CQSZ-1、CQSZ-2）以及巫溪1株（CQWX-1）。

2.4 mutS基因分析

對供試菌株的mutS基因部分序列（682 bp）進行系統(tǒng)發(fā)育分析，并從GenBank 上選取27個菌株的mutS基因的序列（表2）作為參考序列進行比較。由系統(tǒng)發(fā)育分析可知，所有的供試菌株可聚到2個分支中，如圖3所示，第一分支與egl基因的系統(tǒng)發(fā)育圖中的序列變種44分支對應。第二分支則與egl基因的系統(tǒng)發(fā)育圖的序列變種15和17對應。

圖2 基于egl基因部分序列的系統(tǒng)發(fā)育分析Fig.2 Phylogenetic neighbor-joining tree based on partial sequences of the egl gene

圖3 基于mutS基因部分序列的系統(tǒng)發(fā)育分析Fig.3 Phylogenetic neighbor-joining tree based on partial sequences of the mutS gene

3 結論與討論

青枯菌受不同的環(huán)境、寄主植物和各種農藝操作的影響，自身不斷進化，形成了一個復合種[14]。Fegan和Prior所提出的演化型分類框架表明青枯菌的種下分化與地理起源密切相關[10]。在分析美國東南部的青枯菌多樣性時，分離自北卡羅來納州和佐治亞州的煙草青枯菌均屬于生化變種1和演化型Ⅱ（美洲分支）[15]。Wang等在研究中指出中國煙草青枯菌屬于生化變種3[16]。對我國青枯菌遺傳多樣性分析中9株分離自煙草的菌屬于演化型Ⅰ[17]。福建煙草青枯菌被鑒定出均屬于演化型Ⅰ[18]?？傮w來看，我國的煙草青枯菌屬于演化型Ⅰ（亞洲分支），與美國的煙草青枯菌（演化型Ⅱ即美洲分支）屬于不同的地理起源。Wicker等利用多位點序列分析（multilocus sequence analysis, MLSA）追溯了世界范圍內的青枯菌的進化歷史，發(fā)現(xiàn)該世系中存在21個重組事件，而演化型Ⅰ是多態(tài)性持續(xù)發(fā)展的分支之一[19]。重慶地區(qū)種下分化情況也屬于生化變種3及演化型Ⅰ,并在序列變種上存在分化。

重慶12個煙草種植區(qū)縣中，除巫溪、巫山和奉節(jié)地處大巴山區(qū)外，其余9個區(qū)縣均處于武陵山區(qū)。其中，巫山、奉節(jié)目前尚未發(fā)現(xiàn)青枯病，巫溪縣有零星發(fā)生，而另外的9個區(qū)縣均發(fā)生青枯病，且較嚴重。巫山和奉節(jié)所種植的煙草品種與重慶其他植煙地并無差異。同時，本研究的結果未發(fā)現(xiàn)青枯菌序列變種分化與煙草品種有直接關系，由此推測，煙草品種可能不是青枯菌種下分化的直接原因。而武陵山區(qū)與大巴山區(qū)的氣候、土壤條件等地理區(qū)域因子存在明顯差異。鄂西煙區(qū)也大多屬于武陵山區(qū)。劉海龍等對湖北地區(qū)煙草青枯菌系統(tǒng)發(fā)育分析的結果表明，湖北煙草青枯菌菌株屬于青枯菌亞洲分支（演化型Ⅰ）的3個序列變種，分別為序列變種15、17和44，其優(yōu)勢菌株也為序列變種17和44[20]。本研究結果與這個結果基本一致。福建和貴州的煙草青枯菌系統(tǒng)發(fā)育分析結果證實，福建的煙草青枯菌菌株包含序列變種15、17、34和44，而貴州的供試菌株僅被聚到序列變種17分支中[21]。比起重慶、湖北及靠近武陵山區(qū)的貴州，距離較遠的福建省所分離的煙草青枯菌菌株除序列變種15、17和44外還存在序列變種34。而其他地區(qū)是否存在異同還需進一步驗證，因此，對于全國煙草青枯菌分化情況的全面評估是很有必要的。青枯菌與寄主協(xié)同進化過程中，地域分布是否是其種內遺傳分化的主要驅動力，尚有待于大尺度范圍內進一步深入研究。

掌握各地煙草青枯菌的具體分化情況，對今后進一步深入研究青枯菌的系統(tǒng)演化及防治措施具有重要意義。本研究分析了重慶地區(qū)的煙草青枯菌生化變種、演化型及系統(tǒng)發(fā)育情況，證實重慶地區(qū)46株供試菌株均屬于生化變種3和演化型Ⅰ；系統(tǒng)發(fā)育進化分析結果表明重慶煙草青枯菌可聚類到3個序列變種（序列變種15、17和44）中。為進一步研究影響青枯菌分化的關鍵因子及其遺傳多樣性機制提供了相應的依據(jù)。

[1]Yabuuchi E, Kosako Y, Yano I, et al. Transfer of twoBurkholderiaand anAlcaligenesspecies toRalstoniagen. nov.:proposal ofRalstoniapickettii(Ralston, Palleroni and Doudoro ff 1973) comb. nov.,Ralstonia solanacearum(Smith 1896)comb. nov. andRalstoniaeutropha(Davis 1969) comb. nov[J].Microbiology Immunology, 1995, 39(11): 897-904.

[2]Wicker E, Grassart L, Coranson-Beaudu R, et al.Ralstonia solanacearumstrains from Martinique (French west indies)exhibiting a new pathogenic potential[J]. Applied and Environmental Microbiology, 2007, 73(21): 6790-6801.

[3]朱朝賢, 王彥亭, 王智發(fā). 中國煙草病害[M]. 北京: 中國農業(yè)出版社, 2002: 152-153.ZHU Caoxian, WANG Yanting, WANG Zhifa. Tobacco Diseases of China[M]. Beijing: China Agriculture Press, 2002: 152-153.

[4]孔凡玉. 煙草青枯病的綜合防治[J]. 煙草科技, 2003(4): 42-43.KONG Fanyu. Integrated Control of Tobacco Bacterial Wilt Disease[J]. Tobacco science & Technology, 2003(4): 42-43.

[5]Hayward A C. Biology and epidemiology of bacterial wilt caused byPseudomonasSolanacearum[J]. Annual Review of Phytopathology, 1991, 29(1): 65-87.

[6]Buddenhagen I W, Sequeira L, Kelman A. Designation of races inPseudomonasSolanacearum[J]. Phytopathology, 1962, 52(8):726.

[7]Hayward A C. Characteristics ofPseudomonassolanacearum[J].Journal of Applied Microbiology, 1964, 27(2): 265-277.

[8]He L Y, Sequeira L, Kelman A. Characteristics of strains ofPseudomonasSolanacearumfrom China[J]. Plant Disease, 1983,67(12): 1357-1361.

[9]Pegg K G, Moffett M L. Host range of the ginger strain ofPseudomonassolanacearumin Queensland[J]. Animal Production Science, 1971, 11(53): 696-698.

[10]Fegan M, Prior P. How complex is the“Ralstonia solanacearumspecies complex”?[M]// Allen C, Prior P, Hayward A C.Bacterial wilt disease and theRalstonia solanacearumspecies complex. St. Paul, MN: American Phytopathological Society Press, 2005: 449-461.

[11]Prior P, Fegan M. Recent developments in the phylogeny and classi fi cation ofRalstonia solanacearum[J]. Acta horticulturae,2005, 695, 127-136.

[12]Sagar V, Jeevalatha A, Mian S, et al. Potato bacterial wilt in India caused by strains of phylotype I, II and IV ofRalstonia solanacearum[J]. European Journal of Plant Pathology, 2014,138(1): 51-65.

[13]Poussier S, Prior P, Luisetti J, et al. Partial sequencing of the hrpB and endoglucanase genes con fi rms and expands the known diversity within theRalstonia solanacearumspecies complex[J].Systematic and applied microbiology, 2000, 23(4): 479-486.

[14]Genin S, Denny T P. Pathogenomics of theRalstonia solanacearumspecies complex[J]. Annual Review of Phytopathology, 2012, 50: 67-89.

[15]Hong, J C, Norman, D J, Reed, D L, et al. Diversity amongRalstonia solanacearumstrains isolated from the southeastern United States[J]. Phytopathology, 2012, 102(10): 924-936.

[16]Wang Han-cheng, Huang Yan-fei, Wang Jin, et al. Phenotypic Fingerprints ofRalstonia solanacearumBiovar 3 Strains from Tobacco and Tomato in China Assessed by Phenotype MicroArray Analysis[J]. Plant Pathology Journal, 2015, 14(1):38.

[17]Xu J, Pan Z C, Prior P, et al. Genetic diversity ofRalstonia solanacearumstrains from China[J]. European Journal of Plant Pathology, 2009, 125(4): 641-653.

[18]徐進, 顧鋼, 潘哲超, 等. 福建煙草青枯菌演化型及生化變種鑒定研究[J]. 中國煙草學報, 2010, 16(6): 66.XU Jin, GU Gang, PAN Zhechao, et al. Phylotype identi fi cation and biovar determination ofRalstonia solanacearumstrains isolated from Fujian province[J]. Acta Tobacaria Sinica, 2010,16(6): 66.

[19]Wicker E, Lefeuvre P, de Cambiaire J C, et al. Contrasting recombination patterns and demographic histories of the plant pathogenRalstonia solanacearuminferred from MLSA[J]. The ISME journal, 2012, 6(5): 961-974.

[20]劉海龍, 黎妍妍, 鄭露, 李錫宏, 黃俊斌. 湖北地區(qū)煙草青枯菌系統(tǒng)發(fā)育分析[J]. 中國煙草科學, 2015, 36(2): 81-86.LIU Hailong, LI Yanyan, ZHENG Lu, LI Xihong, HUANG Junbin. Phylogenetic Analysis of TobaccoRalstonia solanacearumstrains from Hubei Province[J]. Chinese Tobacco Science, 2015, 36(2): 6.

[21]潘哲超, 徐進, 顧鋼, 等. 福建及貴州等地煙草青枯菌系統(tǒng)發(fā)育分析[J]. 植物保護, 2012, 38(1): 18-23.PAN Zhechao, XU Jin, GU Gang, et al. Phylogeny of tobaccoRalstonia solanacearumstrains from Fujian and Guizhou provinces[J]. Plant Protection, 2012, 38(1): 18-23.

Biovar and sequevar ofRalstonia solanacearumin tobacco growing areas in Chongqing

LIU Ying1, TANG Yuanman1, ZHANG Shuting1, GUO Bing1, GU Jitao1, DING Wei1,2
1 College of Plant Protection, Southwest University, Chongqing 400715, China;2 Chongqing Tobacco Research Institute, Chongqing 400715, China

46 isolates were assessed to explore the diversity of tobaccoRalstonia solanacearumin Chongqing by using biovar, phylotype speci fi c multiplex PCR, and phylogenetic relationships (eglandmutS). Results showed that all tested isolates belonged to biovar 3 and phylotype I. The phylogenetic neighbor-joining tree based on partial sequences ofeglandmutSgenes indicated that the tested isolates were clustered into phylotype I, and were further divided into 3 sequevars, namely sequevars 15, 17, and 44. Sequevar 17 and 44 were predominant.R.solanacearumstrains on tobacco showed genetic differentiation to some extent, which provided basis for studying relationship between classi fi cation and geographic areas.

tobacco;Ralstoniasolanacearum; biovar; phylogenetic analysis;egl;mutS

劉穎，唐元滿，張淑婷，等. 重慶煙區(qū)青枯雷爾氏菌生化變種及序列變種的特性研究[J]. 中國煙草學報，2016,22（4）

中國煙草總公司重點項目(110201202002)；中國煙草總公司重慶市公司重點項目(NY20130501070005)

劉穎（1991—），碩士研究生，研究方向：天然產物農藥及植物與病原菌互作， Email: lyying1201@163.com

丁偉（1966—），博士，教授，研究方向：天然產物農藥和煙草有害生物系統(tǒng)控制，Email: dwing818@163.com

2015-12-11

:LIU Ying, TANG Yuanman, ZHANG Shuting, et al. Biovar and sequevar ofRalstoniasolanacearumin tobacco growing areas in Chongqing [J]. Acta Tabacaria Sinica, 2016,22(4)