Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2009, p. 4539-4549, Vol. 75, No. 13
0099-2240/09/$08.00+0 doi:10.1128/AEM.01336-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Plant Pathology,1 Genome Center,2 Department of Genetics, University of Wisconsin—Madison, Madison, Wisconsin3
Received 3 July 2008/ Accepted 26 April 2009
|
|
|---|
|
|
|---|
Pectobacterium strains isolated from potato are diverse based on serology, genome structure, and fatty acid composition (5, 35). Previous epidemiological studies of pectolytic Enterobacteriaceae were complicated by the diversity of this group and the lack of tools capable of placing all isolates into clades. For example, Gross et al. (14) were unable to classify over 50% of Pectobacterium isolates obtained from potato, and Pitman et al. (23) were unable to type 13% of their isolates. Novel PCR-based methods potentially capable of classifying all Pectobacterium isolates have been described, but they were developed prior to the recognition of P. carotovorum subsp. brasiliensis (1, 34).
The main virulence determinants of Pectobacterium are the pectolytic enzymes secreted through the type II secretion system. Although these enzymes are required for development of symptoms, many other virulence genes have been shown to contribute to Pectobacterium pathogenicity, including the type III secretion system (T3SS) genes, the cfa gene cluster, and the type IV secretion system genes (3, 15, 19). Recent genomic analysis showed that some of these gene clusters, such as the cfa and type IV secretion system cluster genes, as well as genes important for interactions with insects, are present in only some Pectobacterium species (10). Thus, Pectobacterium species appear to use different genetic tools to overcome plant host barriers and to interact with insect vectors.
Many gram-negative pathogenic bacteria secrete virulence proteins, known as effectors, through the T3SS into host cells. Once inside host cells, the effectors manipulate host defenses and promote bacterial growth (13). Unlike many other gram-negative plant pathogens, Pectobacterium does not require the T3SS for pathogenicity. Rather, this secretion system makes a small, but measurable, contribution to the early stages of P. carotovorum growth in leaves of the model plant Arabidopsis thaliana (26) and contributes to the virulence of P. atrosepticum on potato (15). Recently, we isolated Pectobacterium strains that lack the T3SS from potatoes and also found P. wasabiae and P. carotovorum subsp. brasiliensis on potatoes in Wisconsin (35). The first goal of this study was to determine if P. wasabiae and P. carotovorum subsp. brasiliensis are common in agricultural fields or if soft rot disease is typically caused by P. carotovorum subsp. carotovorum and P. atrosepticum, which have been the focus of nearly all previous studies of potato soft rot, stem rot, and blackleg disease. Second, since we recently isolated a strain lacking the T3SS (35), we also aimed to determine if strains lacking the T3SS are common in infected potatoes and if these strains tend to be less virulent on potato stems and tubers than strains encoding a T3SS.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Strains and plasmids used in this study
|
DNA hybridization and PCR analysis.
Pectobacterium T3SS genes were detected by DNA hybridization and PCR amplification with T3SS-specific probes and primers for genes encoding regulatory, structural, and secreted proteins. For DNA hybridization assays, chromosomal DNA was digested with EcoRI, electrophoresed through a 0.8% agarose gel, and transferred to a nylon membrane (Millipore Co., Billerica, MA). Probe DNA was either PCR amplified from P. carotovorum WPP14 genomic DNA (hrpL, hrpN, hrpW-dspE/F, hecB, and a region upstream of hrcU) or amplified from plasmids containing DNA cloned from P. carotovorum WPP14 [a 5.8-kb PCR-amplified fragment containing hrpD through hrpN amplified from pThrpDN(Ecc) and a 3.5-kb PCR-amplified fragment containing hrpB through hrcC amplified from pThrpBC(Ecc)]. The primers and rregions used are described in Table 2. A PCR amplification program consisting of 30 cycles of 94°C for 30 s, 54°C for 30 s, and 72°C for 30 s to 5 min, depending on the size of the amplified fragment, was used. The probes were labeled with a Gene Images Alkphos direct labeling and detection system kit (GE Healthcare Bio-Sciences Corp., Piscataway, NJ). Hybridization results were recorded on X-ray film (Kodak, Los Angeles, CA). Additionally, the presence of hrcC, which encodes the T3SS porin, was detected with the primer set listed in Table 2. 16S-23S ITS-PCR and restriction fragment length polymorphism (RFLP) analyses were performed as described by Toth et al. (34). Briefly, the rRNA intergenic spacer region was amplified with primers G1 and L1, the amplified DNA was digested with RsaI (New England Biolabs, Ipswich, MA), and the DNA was analyzed by gel electrophoresis on 3% GenePure HiRes agarose gels (ISC Bioexpress, Kaysville, UT) in Tris-borate-EDTA buffer.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotides used in this study
|
-glutamylphosphate reductase), were PCR amplified from the type strains of Pectobacterium sp. and Dickeya chrysanthemi and the strains isolated in Wisconsin, including WPP2, WPP5, WPP12, WPP16, WPP20, WPP127, and WPP156 (Table 1), as described previously (20). The sequences of the PCR-amplified fragments were obtained with a BigDye Terminator kit (Perkin-Elmer). Sequence chromatogram output files were initially aligned and edited using SeqMan 5.08 (DNASTAR, Inc., Madison, WI). A phylogenetic analysis of the sequences of 51 bacterial strains was conducted using PAUP* 4.0b10 and the concatenated gene sequence data set. Yersinia sp. strains were used as outgroups in all reconstructions. We used ModelTest 3.7 to select the standard AIC model (TrN with a gamma correction for rate variation among sites [32]) as the best-fitting model of evolution (24, 25). The corresponding likelihood parameters were estimated and then applied in PAUP* 4.0b10. Two different algorithms were then used, maximum likelihood-corrected neighbor joining (NJ) and weighted maximum parsimony (MP) with substitutions weighted according to the instantaneous rate matrix or characters weighted according to their rescaled consistency index values. Bootstrap values were estimated in order to evaluate the support for each clade (1,000 replicates for maximum likelihood-corrected NJ and MP heuristic search) with sampling limited to parsimony-informative characters with a cutoff value of 50.
Pulsed-field gel electrophoresis (PFGE) analysis.
Agarose plugs containing chromosomal DNA of each bacterial strain were prepared and digested with I-CeuI (New England Biolabs, Inc.) as described previously (27, 31, 35). Briefly, a Pulsaphor Plus system with a hexagonal electrode array (Pharmacia, Uppsala, Sweden) was used for electrophoresis by following the manufacturer's instructions. The digested genomic DNA was separated by electrophoresis on a 1% agarose gel at 12°C for 20 h at 5 V/cm with the switch time ramping from 25 to 45 s. ProMega lambda ladder (Promega, Madison, WI) was used as the size marker, and the DNA was stained with ethidium bromide and visualized with a UV light transilluminator.
Virulence assays.
Virulence assays were performed with potato tubers and stems (Solanum tuberosum cv. Atlantic) in the laboratory and field, respectively. The relative virulence of 10 Pectobacterium strains in potato tubers (cv. Atlantic) was determined by measuring the amount of macerated potato tissue (35). Ten to 12 potato tubers were inoculated by placing 10 µl of a 108-CFU/ml bacterial suspension into 1-cm-deep holes poked into the tubers with a pipette tip. The tubers were placed in plastic bags and incubated for 2 days at 28°C. This inoculum level was chosen because it was the lowest level at which the tubers reliably developed symptoms within 2 days with the most virulent strains examined. After 2 days the inoculated tubers were cut open, and the macerated tissue was removed and weighed. Five assays with P. carotovorum subsp. carotovorum WPP14, P. carotovorum subsp. brasiliensis WPP17, P. wasabiae WPP163, and P. carotovorum subsp. brasiliensis WPP165 were conducted in the summers of 2005 and 2006 at the University of Wisconsin Hancock Agricultural Research Station using potato plants (S. tuberosum cv. Atlantic). S. tuberosum cv. Atlantic was chosen because it developed consistent and reliable symptoms in preliminary greenhouse assays. Bacterial treatments were arranged in a randomized complete block design with three blocks. Each plot was 4.5 m (16 rows) by 16.7 m, and every two rows were separated by two boarder rows in which S. tuberosum cv. Red Norland was planted. Twenty potato tubers were planted for each bacterial treatment on 23 April 2005 and 21 April 2006. In 2005 and 2006, 8 to 10 and 15 to 20 potato stems per block, respectively, were inoculated with WPP14, WPP17, WPP163, and WPP165. Prior to inoculation, the Pectobacterium strains were grown overnight on LB agar. The bacteria were scraped from the plates with sterile toothpicks, and then the plants were stab inoculated with the toothpicks. One stem per plant was inoculated with one strain, and a second stem was stabbed with a sterile toothpick as a negative control. The toothpicks were left in the stem wounds, and each wound was sealed with a strip of Parafilm, both to protect the bacteria from desiccation and to mark the inoculated and control stems. One week after inoculation, the stems were collected from the field. The stems were sliced in half, and the length of the lesion in each stem was measured.
Statistical analysis.
Statistical analysis of data was conducted using Statistical Analysis Systems (SAS Institute, Cary, NC). Analysis of variance was determined using the general linear model procedures, and means were separated with the least significant difference.
Nucleotide sequence accession numbers.
The GenBank accession numbers for the sequences of six housekeeping genes of WPP2, WPP5, WPP12, WPP16, WPP20, WPP127, and WPP156 are EU684027 to EU684051 and EU250362 to EU250373. The GenBank accession numbers for the sequences of six housekeeping genes of the five type strains used are FJ895836 to FJ895865.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Description of Pectobacterium strains isolated from seven fields in Wisconsin in 2001 and 2004 based on PFGE and elicitation of an HR on tobacco leaves (N. tabacum L. cv. Xanthi NN)
|
![]() View larger version (80K): [in a new window] |
FIG. 1. Sixteen of 45 Pectobacterium isolates obtained in 2004 and P. wasabiae SCRI488 were unable to elicit an HR on N. tabacum cv. Xanthi. Bacterial isolates were infiltrated into the leaves of 6- to 7-week-old N. tabacum cv. Xanthi, and results were recorded the following day. Leaf panels 1 through 8 were infiltrated with Pectobacterium sp. strains WPP14, WPP17, WPP161, WPP163, WPP165, and WPP172, P. wasabiae SCRI488, and a water control, respectively.
|
![]() View larger version (59K): [in a new window] |
FIG. 2. Pectobacterium isolates unable to elicit an HR did not contain the genes for the Pectobacterium T3SS. (A) The P. carotovorum WPP14 T3SS is indicated by open arrows, and the border genes, hecB, and ECA2075 are indicated by filled arrows. (B) Genomic DNA from isolates unable to elicit an HR did not hybridize to DNA amplified from the P. carotovorum WPP14 T3SS gene cluster, including hrpL, hrpN, hrpW, dspE, dspF, hrpD, hrpE, hrpF, hrpG, hrcC, hrpT, hrpV, hrpB, and hrcJ.
|
![]() View larger version (25K): [in a new window] |
FIG. 5. Character and substitution weight parsimony phylogeny generated by using concatenated sequences of six housekeeping genes of 51 isolates, including 4 Brenneria sp. isolates, 10 Dickeya sp. isolates, 4 Yersinia sp. isolates, and 33 Pectobacterium sp. isolates. The numbers for Pectobacterium sp. indicate the clades, as follows; 1, P. carotovorum subsp. brasiliensis; 2, P. carotovorum subsp. carotovorum; 3, P. atrosepticum; 4, P. betavasculorum; and 5, P. wasabiae. The strains whose designations begin with ATCC are the type strains of species. Primer sets were designed to amplify fragments of six conserved housekeeping genes: acnA, gapA, icdA, mdh, pgi, and proA. Bootstrap values (expressed as percentages for 1,000 replications) are indicated only for branches also retrieved by MP and NJ analyses. Yersinia sp. strains were used as an outgroup for these analyses. Details for strains used in a previous analysis were described by Ma et al. (20), and the new strains are described in Table 1.
|
![]() View larger version (30K): [in a new window] |
FIG. 3. Genomic DNA from WPP17 and WPP172 were unable to hybridize to two genes flanking the T3SS region, ECA2075 and hecB, respectively. Genomic DNA was digested with SacII and EcoRV and then hybridized with hecB or ECA2075 probes. The presence of T3SS and the flanking regions is indicated on the right. Lane 1, WPP14; lane 2, WPP17; lane 3, WPP161; lane 4, WPP163; lane 5, WPP165; lane 6, WPP168; lane 7, WPP172; lane 8, WPP178.
|
![]() View larger version (33K): [in a new window] |
FIG. 4. Presence of the Pectobacterium T3SS was not correlated with the ability to macerate potato tubers (A) or cause disease in stems (B). For the inoculated tubers, the bars indicate the amounts of tissue macerated, and the error bars indicate the standard errors for 30 to 36 replicates. Bars with the same letter are not significantly different according to a least significant difference test at a P value of 0.05. Pc, Pw, and Pb indicate that isolates were identified as P. carotovorum subsp. carotovorum, P. wasabiae, and P. carotovorum subsp. brasiliensis, respectively. For the inoculated stems, the bars indicate lesion lengths 7 days postinoculation, and the error bars indicate the standard errors for 14 to 20 replicates. Bars with the same letter were not significantly different according to a least significant difference test at a P value of 0.05. dH2O, distilled water; NC, control.
|
Multiple Pectobacterium species were present in fields and, occasionally, in individual diseased tubers.
Isolates obtained in 2001 and 2004 were used to determine if infected plants from individual fields were typically infected with one or multiple species. Two previously described DNA-based methods, multilocus sequence analysis (MLSA) (20) and 16S-23S ITS-PCR-restriction fragment length polymorphism (RFLP) analysis (34), were used to determine the species of the isolates (Fig. 5 and Fig. 6). Our MLSA relied on building phylogenetic trees from concatenated gene sequences of six housekeeping genes, acnA, gapA, icdA, mdh, pgi, and proA, and it provided robust support for clades corresponding to the four previously described Pectobacterium species (clades 2 through 5) and for P. carotovorum subsp. brasiliensis (clade 1) (20). Sequences of seven new Pectobacterium isolates and five type strains (four Pectobacterium strains and one Dickeya strain) were added to the MLSA developed by Ma et al. (20). The new isolates were chosen because they produced PFGE or 16S-23S ITS-PCR-RFLP patterns not produced by strains included in the previously described phylogenetic tree and therefore could not easily be identified to the species level. For example, WPP127 and WPP156, both of which were placed in the P. carotovorum subsp. carotovorum clade (Fig. 5), produced a novel 16S-23S ITS-PCR-RFLP pattern (Fig. 6), and their PFGE patterns were different from each other (Table 3 and Fig. 7). The tree in Fig. 5 lacks sequences of 14 Pectobacterium isolates used by Ma et al. (20); thus, approximately two-thirds of the Pectobacterium data are shared by the two trees.
![]() View larger version (49K): [in a new window] |
FIG. 6. 16S-23S ITS-PCR-RFLP patterns of Pectobacterium strains. The 16S-23S ITS region was amplified with primers L1 and G1, and the amplified DNA was digested with RsaI and then analyzed by gel electrophoresis. The strains used are indicated above the lanes, and strains are grouped into the clades defined by the MLSA. dH2O, distilled water.
|
![]() View larger version (53K): [in a new window] |
FIG. 7. Four representative I-CeuI patterns for the 42 Pectobacterium isolates obtained from diseased potato tubers in Wisconsin in 2004. WPP14 was isolated in 2001 from a diseased potato plant with aerial stem rot. Fragments more than 1,000 kb long are not shown. The fields from which the strains originated and the numbers of strains with each pattern are indicated at the bottom. Lane M contained markers.
|
We examined the isolates using the 16S-23S ITS-PCR-RFLP method as described by Toth et al. (34) with the hope that a method that is faster and simpler than MLSA could be used to reliably identify Pectobacterium species. We found two ITS-PCR-RFLP patterns for P. carotovorum subsp. carotovorum strains, one of which was described previously (34) (Fig. 6). Similarly, we found two patterns for P. carotovorum subsp. brasiliensis strains, one of which was described previously (34). The pattern for WPP17, which is at the base of the P. carotovorum subsp. brasiliensis clade and which lacks a T3SS, was yet a third pattern, which was identical to a pattern described by Toth et al. (34). We also found two ITS-PCR-RFLP patterns for P. wasabiae, one of which was described previously (34) and one of which was identical to the pattern obtained for some P. carotovorum subsp. carotovorum strains. WPP19, which is at the base of the P. wasabiae clade, had a pattern similar to that obtained for the majority of the P. wasabiae strains and also similar to that of some of the P. carotovorum subsp. carotovorum strains. Since P. carotovorum subsp. carotovorum strains, such as WPP14, and P. wasabiae strains, such as WPP161, have identical ITS-PCR-RFLP patterns, the 16S-23S ITS-PCR-RFLP method could not be used to reliably place Pectobacterium isolates into taxa.
To determine if Pectobacterium strains isolated from individual tubers were clonal and to determine if PFGE patterns could be correlated with Pectobacterium clades, genomic DNA was digested with endonuclease I-CeuI and separated by PFGE (Fig. 7). I-CeuI specifically targets a 26-bp sequence on 23S rRNA genes of bacterial rrn operons; thus, the number of digested fragments represents the number of rrn operons on the chromosome. PFGE profiles were obtained for 42 of the 45 isolates obtained in 2004 examined (Table 3 and Fig. 7). As expected, digestion of all Pectobacterium isolates yielded five or six fragments ranging from 50 to 1,000 kb long (Fig. 7) and a large fragment over 3,000 kb long (not shown) (35). A total of four distinct I-CeuI PFGE pulsotypes were observed, and all isolates from each individual tuber had identical patterns. Unfortunately, we were unable to obtain PFGE data for WPP168, which was a T3SS-deficient isolate obtained from a tuber also containing T3SS-encoding isolates. Strains with identical PFGE pulsotypes were also obtained from fields A and C; thus, isolates with the same pulsotype may be obtained from multiple locations. We gave the PFGE pulsotypes arbitrary numbers, and only pulsotypes 1 and 8 were seen in 2001 and 2004. P. carotovorum Ecc71, a model strain used for many years to study P. carotovorum genetics, was also a pulsotype 1 strain (35). The P. wasabiae PFGE pattern differed from the other patterns in that only five fragments smaller than 1,000 kb were apparent, suggesting that strains in the P. wasabiae clade contain six rRNA operons rather than the seven operons found in the other Pectobacterium species. WPP17, which is most closely related to P. carotovorum subsp. brasiliensis and which also lacks a T3SS, also appears to contain only six rRNA operons (35).
|
|
|---|
We also found that Pectobacterium strains lacking the T3SS can be isolated from diseased tubers and that the T3SS-deficient strains are still virulent in potato tubers and stems. A recent draft sequence of WPP163 shows that the T3SS is indeed not present in this strain (N. T. Perna, unpublished data; NCBI Genome Project ID 31293). Since we used high inoculum levels, we cannot rule out the possibility that the T3SS-encoding Pectobacterium strains may be able to better colonize plants or cause disease when bacteria are inoculated at lower concentrations, under environmental conditions, or into hosts or tissues that were not examined. It should not be assumed the T3SS-deficient P. wasabiae isolates are secondary invaders since P. wasabiae was the only species isolated from some samples, and P. wasabiae is capable of causing disease in both stems and tubers.
Recent genomic analyses of three Pectobacterium genomes revealed that gene clusters identified as clusters important for virulence in P. atrosepticum are not present in one or both sequenced P. carotovorum genomes (10). Similarly, our results show that Pectobacterium sp. strains lacking the T3SS, which has been reported to play a role in Pectobacterium virulence (15, 26), can be isolated from diseased tubers. The T3SS-deficient strains are likely to have lost a T3SS rather than to never have had one since closely related T3SS genes are present in the same relative chromosomal position in both Dickeya and Pectobacterium, suggesting that the T3SS was acquired by the common ancestor of these two soft rot genera. The mosaic of the virulence genes present in these Pectobacterium species suggests that different strains use different tools to overcome plant barriers. P. wasabiae WPP172, but not the other P. wasabiase strains, also lacks at least one gene, hecB, that flanks the T3SS gene cluster in Pectobacterium and Dickeya and that plays a role in Dickeya virulence (28). Finally, both the T3SS-deficient strains in the P. wasabiae clade and P. carotovorum WPP17 appear to contain only six rRNA operons, rather than the seven rRNA operons found in all other Pectobacterium strains analyzed. Spontaneous deletion of rrn operons has been reported for Yersinia species as well (6). It appears that deletion of an rrn operon, as well as deletion of the T3SS, occurred at least twice in the genus Pectobacterium.
Strains lacking a functional T3SS have been reported for the plant pathogens Pseudomonas syringae (22) and Erwinia pyrifoliae (18). P. wasabiae differs from these plant pathogens since it lacks a T3SS, but it is still virulent. Also, in contrast to T3SS-deficient P. syringae strains, which are closely related to each other and separate from pathogenic P. syringae (22), T3SS-deficient Pectobacterium strains fall into two species, P. carotovorum and P. wasabiae. The T3SS-deficient Pectobacterium strains may be most analogous to the animal pathogen Pseudomonas aeruginosa. Strains of this human pathogen isolated from patients with chronic lung infections are typically T3SS deficient, even though 90% of environmental P. aeruginosa isolates encode a T3SS (17).
The reported host range of Pectobacterium expands regularly, but it is challenging to identify the Pectobacterium species isolated from newly identified hosts, as well as from commonly studied hosts, such as potato. Of the three DNA-based methods that we used to characterize isolates, MLSA provided the most unambiguous results and allowed us to place isolates into defined clades with strong branch support (20), while the other two methods could not be used to unambiguously place isolates into taxa. All three methods tested suffer from the limitation that they can be used only with pure cultures. However, combination of the MLSA data with genome sequence data should aid in the development of sensitive assays capable of detecting and classifying Pectobacterium strains directly from field samples, thereby allowing tracking of these taxa in the environment that is more efficient than the tracking that can be done with currently available assays.
The MLSA method appears to be robust for classification of Pectobacterium. One-third of the Pectobacterium isolates in the MLSA tree differ from the isolates used previously (20), and this did not change the clades into which the remaining Pectobacterium isolates were placed. We defined clade 2 as a clade that includes both P. carotovorum subsp. carotovorum and P. carotovorum subsp. P. carotovorum subsp. odoriferum (SCRI482). P. carotovorum subsp. odoriferum was established based on DNA hybridization with three P. carotovorum strains (8), one of which was later determined not to be a P. carotovorum strain (9). Further review of P. carotovorum phylogeny is needed to determine if maintaining this subspecies designation is appropriate. As shown previously (20), Brenneria isolates did not form a monophyletic clade. Sequence data for additional Brenneria isolates are required to resolve how they are related to each other and to the soft rot enterobacteria.
This study was supported by National Science Foundation award 0412599 ("BE/GenEn: Genome-Enabled Analyses of Natural Populations of Pathogens on Natural Hosts") and funds from the Wisconsin Potato and Vegetable Growers.
Published ahead of print on 1 May 2009. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»