潘哲超， 徐 進(jìn)， 顧 鋼， 吳 畏，許景升， 陳順輝， 馮 潔＊

（1．中國農(nóng)業(yè)科學(xué)院植物保護(hù)研究所，植物病蟲害生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京 100193；2．云南省農(nóng)業(yè)科學(xué)院經(jīng)濟(jì)作物研究所，昆明 650205； 3．福建省煙草農(nóng)業(yè)科學(xué)研究所，福州 350003；4．西南大學(xué)植物保護(hù)學(xué)院，重慶 400715）

煙草青枯病是由茄科雷爾氏菌［Ralstonia sol anacear u m （Smit h）Yabuuchi et al．］引起的煙草生產(chǎn)上最具毀滅性的病害之一。該病1880年發(fā)現(xiàn)于美國的北卡羅來納州，1940年前后，該病于美國和印度尼西亞猖獗危害，此后，又逐漸成為日本、澳大利亞、韓國等許多產(chǎn)煙國家煙草上重要的細(xì)菌病害［1］。煙草青枯病在我國長江流域及其以南煙區(qū)普遍發(fā)生，其中以廣東、福建、湖南、四川及貴州煙區(qū)危害最為嚴(yán)重，個(gè)別年份的暴發(fā)流行可造成毀滅性損失［2］。近些年來，該病逐漸有擴(kuò)展蔓延至北方煙區(qū)的趨勢，山東、河南、陜西及遼寧等省均有該病發(fā)生的報(bào)道，且局部地區(qū)危害較重［3］。

傳統(tǒng)的分類框架根據(jù)青枯菌的寄主范圍或?qū)?種雙糖和3種己醇的氧化利用情況，將其劃分為5個(gè)不同的生理小種（race）或5個(gè)生化變種（biovar）［4－7］。煙草青枯菌包括青枯菌1號(hào)生理小種的生化變種Ⅰ、Ⅲ和Ⅳ（即R1／Bv1、R1／Bv3和R1／Bv4菌株）。

2005年，F(xiàn)egan和Prior［8］提出了青枯菌演化型分類框架（phylotype classification shemes），青枯菌被劃分為與地理起源密切相關(guān)的4個(gè)演化型—即演化型Ⅰ、Ⅱ、Ⅲ和Ⅳ型，在地理起源上依次對(duì)應(yīng)于亞洲分支、美洲分支、非洲分支以及印尼分支；每個(gè)演化型又由諸多序列變種（sequevar）構(gòu)成。序列變種定義為：以內(nèi)切葡聚糖酶基因（endogl ucanase，egl）進(jìn)行系統(tǒng)發(fā)育分析時(shí)，一組序列高度保守的菌株集合體。演化型分類框架是建立在青枯菌基因組長期緩慢累積而形成的遺傳變異的基礎(chǔ)上，與傳統(tǒng)的小種及生化變種分類方法相比，演化型的分類框架可以更精確地反映青枯菌這一復(fù)合種的地理起源及種內(nèi)的遺傳多樣性。國際上的研究表明無毒基因avr A是控制青枯菌對(duì)煙草致病能力的決定因素之一，其于青枯菌群體中的分布也存在較大的差異性［9］，亞洲分支青枯菌無毒基因avr A存在RS1000類型和GMI1000類型的兩種類型［10］。

國內(nèi)目前對(duì)煙草青枯病菌系的研究大都遲滯于基于傳統(tǒng)方法的分類框架上［11－17］，僅鄭向華等［18］采用RAPD對(duì)廣東省煙草青枯菌的遺傳多樣性進(jìn)行了研究。本研究采用國際上最新的演化型分類框架，對(duì)福建及貴州部分地區(qū)采集的煙草青枯病菌株進(jìn)行系統(tǒng)發(fā)育分析并進(jìn)行序列變種鑒定，并對(duì)各序列變種的無毒基因avr A進(jìn)行了序列類型的分析，為我國煙草青枯病的防控工作及煙草的抗病育種提供一定的基礎(chǔ)。

1 材料與方法

1．1 供試菌株

煙草青枯菌株由本實(shí)驗(yàn)室采集分離。供試62個(gè)煙草青枯菌株和15個(gè)參考菌株見表1。

表1 供試煙草青枯菌株和參考菌株1）

續(xù)表1

1．2 供試菌株egl基因的擴(kuò)增

用 于 擴(kuò) 增egl 的 引 物 為：Endo－F （5′－ATGCATGCCGCTGGTCGCCGC－3′） 和 Endo－R （5′－GCGTTGCCCGGCACGAACACC－3′）［19］。 所 用 引物由上海生工生物公司合成。

PCR反應(yīng)體系組分（25μL）包含1．5 mmol／L MgCl2，250μmol／L d NTPs，50 mmol／L KCl，10 mmol／L Tris－HCl及1．25 U Taq DNA聚合酶，引物各4 p mol，50 ng模板DNA。

PCR反應(yīng)程序：96℃預(yù)變性9 min；95℃變性1 min，64℃復(fù)性1 min，72℃延伸2 min，30個(gè)循環(huán)；最后72℃延伸10 min，4℃保存［19］。PCR擴(kuò)增產(chǎn)物回收后交由上海生工公司測序。

1．3 avr A基因的擴(kuò)增

用于擴(kuò)增avr A的引物寡核苷酸序列為：avr AA1（5′－TGTAAAACGACGGCCAGTATGAGAAGAATYG GMAA－3′）和 avr AB1（5′－CAGGAAACAGCTATGACCTCGCTRTCGCTATCGYT－3′）［20］。

PCR反應(yīng)程序：94℃預(yù)變性3 min，94℃變性1 min，58℃復(fù)性1．5 min，72℃延伸1 min，30個(gè)循環(huán)；最后72℃延伸10 min，4℃保存［20］。

1．4 數(shù)據(jù)分析方法

測序結(jié)果經(jīng)Bio Edit［21］軟件編輯后應(yīng)用Clustal W 軟件［22］對(duì)egl序列比對(duì)，用 MEGA4．1軟件［23］進(jìn)行系統(tǒng)發(fā)育分析，采用模型鄰接法（NJ）構(gòu)建系統(tǒng)發(fā)育樹［24］，自舉分析為1 000次重復(fù)抽樣檢驗(yàn)得到的自舉置信度。從Gen Bank選取12個(gè)已確定序列變種歸屬的參考菌株egl序列作為參考，以確定參試菌株在系統(tǒng)發(fā)育樹中的位置；同時(shí)選取各序列變種代表性菌株及5個(gè)已確定avr A序列的參考菌株進(jìn)行avr A基因氨基酸序列比對(duì)。

2 結(jié)果與分析

2．1 egl基因系統(tǒng)發(fā)育學(xué)的分析結(jié)果

參試煙草青枯病菌株的egl基因序列分布于亞洲分支以下的4個(gè)亞組內(nèi)（圖1）。采用已知序列變種的標(biāo)準(zhǔn)菌株序列作為參照，4個(gè)亞組分別被鑒定為序列變種15、17、34和44。在該4個(gè)序列變種中，序列變種15所包含的菌株數(shù)在全部供試菌株中所占比例最高，為38．7%。序列變種17略次之，為37．1%。序列變種34再次之，占19．4%。序列變種44最低，僅占4．9%。

圖1 基于egl基因的系統(tǒng)發(fā)育分析

序列變種15的標(biāo)準(zhǔn)菌株為我國臺(tái)灣的番茄青枯病菌株P(guān)SS358，該序列變種包括了24個(gè)來源于福建和廣西的分離菌株，寄主品種分別為‘K326’、‘翠碧1號(hào)’、‘云煙85’以及‘云煙87’；以來源于廣西花生青枯病菌株P(guān)11作為標(biāo)準(zhǔn)菌株的序列變種17，包括23個(gè)分離自福建、廣西和貴州的菌株，分離寄主品種分別為‘K326’、‘翠碧1號(hào)’和‘紅花大金元’。其中分離自貴州的煙草青枯病菌株全部集中于該序列變種內(nèi)。

序列變種34的標(biāo)準(zhǔn)菌株為我國臺(tái)灣番茄青枯病菌株P(guān)SS219，該序列變種包括的參試菌株數(shù)目較少，且12個(gè)菌株全部分離自福建，分離品種為‘翠碧1號(hào)’和‘云煙87’；序列變種44僅包括3個(gè)分離菌株，其標(biāo)準(zhǔn)菌株為廣東油橄欖青枯病菌株O3。

2．2 avr A基因氨基酸序列比對(duì)結(jié)果

對(duì)各序列變種代表性菌株Tb2K91（序列變種34）、Tb2 K94（序列變種15）、Tb2 K99（序列變種44）和Tb2K918（序列變種17）及5個(gè)已確定avr A序列的參考菌株GMI1000、IPO1609、UW551、MOLK2和RS1000進(jìn)行了氨基酸比對(duì)，結(jié)果表明4個(gè)序列變種代表菌株與亞洲菌株RS1000的序列一致，且4個(gè)序列變種代表菌株的avr A氨基酸序列并不存在差異（圖2）。

圖2 基于avr A基因的氨基酸序列比對(duì)

3 討論

青枯菌寄主范圍極為廣泛，可侵染54個(gè)科的450余種植物［25］；不同地理起源的青枯菌在與其寄主長期協(xié)同進(jìn)化的過程中，演化出明顯的生理分化或菌系多樣性。以致美國科學(xué)家Buddenhagen在文章中如是寫道：“有如此之多的青枯病，和如此之多的青枯菌”［26］。

Xu等基于演化型分類框架，對(duì)國內(nèi)13個(gè)省份、17種不同寄主植物上的286個(gè)青枯菌株材料的研究結(jié)果表明：中國青枯菌種內(nèi)具有豐富的遺傳多樣性，存在亞洲分支的10個(gè)序列變種及1個(gè)美洲分支的序列變種［27］。在此基礎(chǔ)上，本研究采用該框架對(duì)福建、貴州等地區(qū)62個(gè)煙草青枯病菌株進(jìn)行了系統(tǒng)發(fā)育分析，結(jié)果表明：62個(gè)參試煙草青枯病菌株全部歸屬于亞洲分支的4個(gè)已知序列變種，即序列變種15、17、34和44；未發(fā)現(xiàn)歸屬于美洲或非洲分支的煙草青枯病菌株（即傳統(tǒng)上的R1／Bv1菌株）。其中序列變種15和17為優(yōu)勢菌系，占全部參試菌株數(shù)的75．8%。值得注意的是，作為新鑒定出的序列變種，以我國臺(tái)灣番茄青枯病菌株P(guān)SS219為標(biāo)準(zhǔn)菌株的序列變種34雖包含的參試菌株數(shù)目較少，但無一例外均系分離自福建省的菌株；同時(shí)，分離于貴州的菌株全部歸屬于序列變種17。但由于本研究采集的菌株數(shù)目有限，這種序列變種與國內(nèi)省份或地區(qū)之間的對(duì)應(yīng)關(guān)系尚需進(jìn)一步研究驗(yàn)證。

某些狹窄寄主范圍的菌株歸屬于特定的序列變種。如來源于香蕉的2號(hào)小種菌株屬于序列變種3、4、6、24；來源于馬鈴薯的3號(hào)小種菌株屬于序列變種1、2，來源于姜的4號(hào)小種菌株屬于序列變種16，來源于桑的5號(hào)小種的菌株屬于序列變種12、44、48。煙草青枯菌雖屬于廣泛寄主的生理小種1，但其在致病性上具有一定的特殊性。如1975年，煙草青枯病成為美國南、北卡羅來納州煙草生產(chǎn)上最主要的病害；但毗鄰的佐治亞和佛羅里達(dá)州卻鮮有該病發(fā)生危害的報(bào)道，盡管在這些地區(qū)番茄青枯病極為普遍［20］。何禮遠(yuǎn)對(duì)中國青枯菌致病性測定的研究結(jié)果表明：絕大多數(shù)非煙草寄主的青枯菌株不能對(duì)煙草致病；但煙草菌株Tb2不僅對(duì)煙草致病，還對(duì)茄科、花生和辣椒等其他寄主植物具不同程度的致病能力［28］。本研究表明福建、貴州等地區(qū)煙草青枯菌存在于序列變種15、17、34和44，但是否僅存于這4個(gè)序列變種尚需今后更多的煙草菌株來驗(yàn)證。

本研究進(jìn)一步對(duì)采集的煙草青枯病菌4個(gè)序列變種的代表性菌株進(jìn)行了氨基酸比對(duì)分析。結(jié)果表明，4個(gè)序列變種代表菌株的avr A氨基酸序列并不存在差異，都屬于日本菌株RS1000類型而非GMI1000類型。但因?yàn)檫@兩種類型的無毒基因avr A讀碼框是一致的，所以無毒基因avr A類型與青枯菌株對(duì)煙草致病力之間是否存在相關(guān)性尚需明確。

作為青枯病防控策略的首選措施之一，抗病育種一直以來備受各國研究者的青睞。但由于青枯菌群體所特有的高度適應(yīng)性及變異性，又成為長期困擾抗病育種研究人員的難題。故而，成功培育抗性品種的先決條件是揭示出青枯菌這一復(fù)雜群體內(nèi)部豐富的遺傳多樣性［29］。國際上以美國為代表的煙草抗青枯病育種工作，主要是針對(duì)演化型Ⅱ型（美洲分支）煙草青枯病菌株開展的。本研究的結(jié)果揭示出國內(nèi)對(duì)煙草致病的菌株均系演化型Ⅰ型菌株（亞洲分支），因此亟需評(píng)價(jià)國內(nèi)煙草青枯病菌4個(gè)序列變種對(duì)現(xiàn)有煙草抗病品種的毒性水平差異，藉以指導(dǎo)現(xiàn)有品種的合理布局；同時(shí)，根據(jù)我國煙草青枯菌的遺傳背景選育新的抗病品種。

［1］ 朱賢朝，王彥亭，王智發(fā)．中國煙草病害［M］．北京：中國農(nóng)業(yè)出版社，2002．

［2］ 陳瑞泰，朱賢朝，王智發(fā)，等．全國16個(gè)主產(chǎn)煙?。▍^(qū)）煙草侵染性病害調(diào)研報(bào)告［J］．中國煙草科學(xué)，1997（4）：1－7．

［3］ 孔凡玉．煙草青枯病的綜合防治［J］．煙草科技，2003（4）：42－48．

［4］ Buddenhagen I，Sequeira L，Kel man A．Designation of races in Pseudomonas sol anacear u m ［J］． Phytopat hology，1962，52：726．

［5］ He L Y，Sequeira L，Kel man A．Characteristics of strains of Pseudomonas sol anacear u m ［J］．Plant Disease，1983，67：1357－1361．

［6］ Pegg K G，Moffett M．Host range of the ginger strain of Pseudomonas sol anacear u m in Queensland［J］．Australian Jour nal of Experi mental Agriculture and Ani mal Husbandry，1971，11：696－698．

［7］ Hay ward A C．Characteristics of Pseudomonas sol anacear u m［J］．Jour nal of Applied Bacteriology，1964，27：265－277．

［8］ Fegan M，Prior P．How complex is the“Ralstonia solanacear um species complex”［M］∥Allen C，Prior P，Hay ward A C．Bacterial wilt disease and the Ralstonia solanacear um species complex．St．Paul：APS，2005：449－462．

［9］ Pouey miro M，Cunnac S，Barberis P，et al．Two typeⅢsecretion syste m effectors fro m Ralstonia sol anacear u m GMI1000 deter mine host－range specificity on tobacco［J］． Molecular Plant－Microbe Interact，2009，22（5）：538－550．

［10］Liu Y，Kanda A，Kiba A，et al．Distribution of avir ulence genes avr A and pop P1 in 22 Japanese phylotype I strains of Ralstonia sol anacear um［J］．Journal of General Plant Pathology，2009，75（5）：362－368．

［11］鄧正平，匡傳富，周志成，等．湖南煙草青枯病菌生理小種測定［J］．湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2004，30（1）：47－49．

［12］顧江濤，許大鳳，李英，等．安徽皖南煙區(qū)青枯病病原菌生化型研究［J］．中國煙草科學(xué)，2008，29（3）：60－61．

［13］黃福新，陳永惠，華靜月，等．廣西草青枯菌系及其主要生理特性研究［J］．植物保護(hù)學(xué)報(bào)，1998，25（3）：240－244．

［14］劉瓊光，陳澤鵬，董春，等．廣東煙草青枯病菌菌系研究［J］．中國煙草學(xué)報(bào)，1999，5（4）：25－27．

［15］張建華，鄭繼法，陳瑞泰，等．我國南方主要煙區(qū)煙草青枯病菌的生物型及產(chǎn)細(xì)菌特性的研究［J］．煙草學(xué)刊，1990（2）：85－93．

［16］鄒陽，肖崇剛，易龍．重慶地區(qū)煙草青枯病菌的生物型測定［J］．煙草科技，2007（6）：62－64．

［17］鄒陽，肖崇剛．重慶地區(qū)煙草青枯菌生理小種的初步鑒定［J］．煙草科技，2008（5）：60－64．

［18］鄭向華，鄧海濱，劉瓊光，等．廣東省煙草青枯菌的菌系和遺傳多樣性［J］．華中農(nóng)業(yè)大學(xué)學(xué)報(bào)，2007，26（4）：463－468．

［19］Poussier S，Prior P，Luisetti J，et al．Partial sequencing of the hrp B and endoglucanase genes confir ms and expands the known diversity within the Ralstonia solanacearum species complex［J］．Systematic and Applied Microbiology，2000，23：479－486．

［20］Robertson A E，Wechter W P，Denny T P，et al．Relationship bet ween avir ulence gene （avr A）diversity in Ralstonia sol anacearu m and bacterial wilt incidence［J］．Molecular Plant－Micr obe Interact，2004，17：1376－1384．

［21］Hall T A．Bio Edit：a user－friendly biological sequence alignment editor and analysis pr ogra m for Windows 95／98／NT［C］．Nucleic Acids Sy mposiu m Series 41，1999．

［22］Tho mpson J D，Higgins D G，Gibson T J．CL USTAL W：i mpr oving the sensitivity of progressive multiple sequence alignment t hr ough sequence weighting，positions－specific gap penalties and weight matrix choice［J］．Nucleic Acids Research，1994，22：4673－4680．

［23］Ku mar S，Ta mura K，Nei M．MEGA3：integrated soft ware for molecular evol utionary genetics analysis and sequence alignment［J］．Brief Bioinfor m，2004，5：150－163．

［24］Jukes T H，Cantor C R．Evol ution of pr otein molecules［M］∥Munro H N．Mammalian protein metabolism．New Yor k，USA：Academic Press，1969：121－132．

［25］Wicker E，Grassart L，Coranson－Beaudu R，et al．Ralstonia sol anacear u m Strains fro m Martinique（French West Indies）Exhibiting a New Pathogenic Potential［J］．Applied and Envir on mental Microbiology，2007，71：6790－6801．

［26］Buddenhagen I W．Bacterial wilt revisited［R］∥Persley G．Bacterial wilt disease in Asia and the South Pacific：proceedings of an international workshop held at PCARRD，Los Banos，Philippines，8－10 October 1985，Canberra：ACIAR，1986．

［27］Xu J，Pan Z C，Prior P，et al．Genetic diversity of Ralstonia sol anacear u m strains from China［J］．Eur opean Jour nal of Plant Pathology，2009，186：6186－6197．

［28］He L Y，Sequeira L，Kel man A．Characteristics of strains of Pseudomonas solanacear u m fro m China［J］．Plant Disease，1983，67：1357－1361．

［29］Gillings M R，F(xiàn)ahy P．Genomic fingerprinting：towards a unified view of the Pseudomonas sol anacear u m species co mplex［M］∥Hay ward A C，Hart man G L．Bacterial wilt：the disease and its causative agent，Pseudomonas solanacear u m，Wallingford，United Kingdom：CAB International，Wallingfor d，1994：95－112．