Multilocus Variable-Number-Tandem-Repeats Analysis (MLVA) distinguishes a clonal complex of Clavibacter michiganensis subsp. michiganensisstrains isolated from recent outbreaks of bacterial wilt and canker in Belgium
BMC Microbiology volume 13, Article number: 126 (2013)
Clavibacter michiganensis subsp. michiganensis (Cmm) causes bacterial wilt and canker in tomato. Cmm is present nearly in all European countries. During the last three years several local outbreaks were detected in Belgium. The lack of a convenient high-resolution strain-typing method has hampered the study of the routes of transmission of Cmm and epidemiology in tomato cultivation. In this study the genetic relatedness among a worldwide collection of Cmm strains and their relatives was approached by gyrB and dnaA gene sequencing. Further, we developed and applied a multilocus variable number of tandem repeats analysis (MLVA) scheme to discriminate among Cmm strains.
A phylogenetic analysis of gyrB and dnaA gene sequences of 56 Cmm strains demonstrated that Belgian Cmm strains from recent outbreaks of 2010–2012 form a genetically uniform group within the Cmm clade, and Cmm is phylogenetically distinct from other Clavibacter subspecies and from non-pathogenic Clavibacter-like strains. MLVA conducted with eight minisatellite loci detected 25 haplotypes within Cmm. All strains from Belgian outbreaks, isolated between 2010 and 2012, together with two French strains from 2010 seem to form one monomorphic group. Regardless of the isolation year, location or tomato cultivar, Belgian strains from recent outbreaks belonged to the same haplotype. On the contrary, strains from diverse geographical locations or isolated over longer periods of time formed mostly singletons.
We hypothesise that the introduction might have originated from one lot of seeds or contaminated tomato seedlings that was the source of the outbreak in 2010 and that these Cmm strains persisted and induced infection in 2011 and 2012. Our results demonstrate that MLVA is a promising typing technique for a local surveillance and outbreaks investigation in epidemiological studies of Cmm.
Clavibacter michiganensis subsp. michiganensis, a Gram positive bacterium, is the causative agent of bacterial canker and wilting, one of the most destructive bacterial diseases in tomato . Contaminated tomato seeds are considered to be the main source of infection. The bacterium survives for a long period of time in seeds, soil and plant debris [2, 3]. Every year, new or reoccurring outbreaks are detected causing substantial economic losses worldwide . Bacterial canker was described for the first time in 1905 in Michigan, USA, and since that moment it has been reported in nearly all tomato growing areas of the world . Difficulties in controlling the spread of the pathogen, the lack of resistant tomato varieties and severity of disease symptoms led to the classification of Cmm as quarantine organisms. Cmm is listed as an A2 quarantine pest by the European and Mediterranean Plant Protection Organization (EPPO)  in Europe and in many countries all over the world .
The epidemiology and the population structure of Cmm in areas where outbreaks of Cmm are common remains scantily investigated and poorly understood. Recent studies describing outbreaks of Cmm in Europe and Asia [5–8] have shed some light on this issue. In Italy a clonal population of Cmm was responsible for the outbreak in 2007 . A high homogeneity was also observed among strains isolated from 2002 to 2007 in Canary Islands suggesting a single introduction of the pathogen as a source of infection . Primary infections in many countries were attributed to the introductions of contaminated tomato seeds and/or seedlings [7, 10]. These findings indicate that seeds play an important role in long-distance spread of the pathogen. A direct link between tomato cultivar, year or place of isolation and Cmm type mostly could not be recognized [6, 8, 9] except the outbreak in 2001 in Turkey where bacterial canker was detected only on one tomato cultivar ‘Target’ . Interestingly, in Israel and Serbia Cmm strains showing the same haplotypes were repeatedly isolated from the same locations during several subsequent years [7, 10]. Reoccurring outbreaks suggest that despite intensified efforts for eradication, reliable control of this disease remains an unattainable goal. The limited progress in improving its management is mainly due to the sporadic nature of the disease outbreaks and to limited and scattered epidemiological data. Therefore, access to an accurate, efficient and cost-effective strain typing technique could be very useful.
Bacterial typing techniques are applied to quickly and reliably differentiate closely related strains in an epidemiological survey, to determinate the relatedness among the strains and to track their origin and pathways of spread. Over the past decades a variety of different typing methods have been developed to generate strain-specific patterns. They are also applied for comprehensive investigation of bacterial population structure and dynamics. A range of methods has already been applied to study the diversity of Clavibacter, particularly to investigate Cmm strains. Rep-PCR (repetitive-element-based PCR), a relatively easy and fast technique, was shown to be of moderate utility , mainly because of the lack of a database and the rather low discriminatory power needed to study closely related strains. Moreover, rep-PCR is mostly not portable between different laboratories . PFGE (pulsed-field gel electrophoresis of macro-restricted bacterial DNA), one of the oldest techniques used in epidemiology, is labor intensive and expensive but is still used as a gold standard in typing of some bacterial species [10, 13]. PFGE was applied to study the diversity of Cmm strains from outbreaks in Serbia  and in Israel  where the results of PFGE showed similar resolution of those obtained by gene sequence analysis and rep-PCR, respectively. Also, AFLP, a high resolution molecular typing method was applied by De Leon and coworkers to study genetic diversity of Cmm strains from Canary Islands . This technique generated more bands per strain and resulted in more reproducible and robust discriminatory clustering of the strains . Highly reproducible multilocus sequence typing (MLST) was used to analyze Cmm population from Serbia. Cmm strains were divided into seven groups and the results were confirmed by PFGE analysis .
MLVA (Multiple-Locus Variable number tandem repeat Analysis) is a PCR-based typing technique that has been widely applied in medical microbiology . It takes advantage of the inherent variability encountered in regions with a number of tandem repeats. The origin of the repetitive regions can be accounted to slipped strand mispairing events occurring during DNA duplication, in which repetitive regions are incorrectly copied resulting in deletion or insertion of one or several copies of the repeat . PCR primers designed to board different VNTR (Variable Number of Tandem Repeats) regions in the genome can be easily combined in a multiplex PCR in an MLVA scheme. The differences between strains are assessed by the different lengths of the repeats visualized by gel electrophoresis or automated fragment analysis on a sequencer. From these sizes, the number of repeat units at each locus can be deduced. The resulting information forms a strain-specific numerical code which can be easily compared to a reference database. The MLVA technique was introduced to bacterial typing as a promising alternative or a complement to already existing typing methods such as AFLP, MLST, rep-PCR or PFGE. The discriminatory power of MLVA is generally higher than other standard typing techniques . However, the final result is group dependent and can vary considerably between different bacterial species. VNTRs have been used to discriminate among individual strains within many food-borne pathogens with little genetic differences, including Escherichia coli O157:H7  and Vibrio cholerae and to study other important human pathogens, such as Neisseria gonorrhoeae, Streptococcus pneumoniae, and Mycobacterium tuberculosis. MLVA has been extensively used for tracking transmissions of important human and animal pathogens [22, 23] and for typing monomorphic bacterial pathogens including Bacillus anthracis and Yersinia pestis. To date, several MLVA schemes have been published on plant pathogens such as Xanthomonas citri pv. citri, X. oryzae pv. oryzicola, Pseudomonas syringae pv. maculicola and tomato, Xylella fastidiosa and on fungi e.g. Aspergillus flavus, but not for Clavibacter subspecies. In plant pathogens, such as Xanthomonas arbolicola pv. pruni, MLVA was proposed as a complementary molecular typing method to AFLP, BOX and ERIC-PCR . In the epidemiological study of pathotypes of Xanthomonas citri MLVA was compared to AFLP and insertion sequence ligation-mediated PCR (IS-LM-PCR) and was found the best method to describe the variations among strains originating from the same country or group of neighboring countries .
The objectives of this study were: 1) to characterize a Belgian population of Cmm strains by a newly developed MLVA scheme; 2) to compare its genetic variability with some strains of Cmm isolated in other countries; 3) to investigate whether the strains responsible for bacterial canker outbreaks in Belgium in 2010–2012 have one or several infection sources and 4) to assess the genetic relatedness of the Cmm strains from Belgium by gyrB and dnaA gene sequence analysis.
The bacterial strains used in this study are listed in Table 1. The strains were obtained from the BCCM/LMG Bacteria Collection (Ghent, Belgium), the GBBC (ILVO Plant Clinic, Merelbeke, Belgium) and the PD collection (Wageningen, The Netherlands). The Clavibacter strain subset consisted of five type strains Cmm LMG 7333T (species type strain), Clavibacter michiganensis subsp. nebraskensis (Cmn) LMG 5627T, Clavibacter michiganensis subsp. sepedonicus (Cms) LMG 2889T, Clavibacter michiganensis subsp. insidiosus (Cmi) LMG 3663T, Clavibacter michiganensis subsp. tessellarius (Cmt) LMG 7294T, two non-pathogenic Clavibacter-like strains and fifty five Cmm originating from Belgian outbreaks and other geographical locations. Twenty three Cmm strains were sampled from symptomatic tomato plants in fields and greenhouses in northeast Belgium. They were isolated from five different tomato cultivars and seven different locations, in the period February 2010 till February 2012 (Table 1). Clavibacter-like isolates from tomato seed are phenotypically similar to Cmm in the common diagnostic semi-selective media and are identified as Cmm in the standard tests but are non-pathogenic to tomato [32, 33]. They were isolated according to the current method for detection of Cmm in tomato seed recommended by International Seed Federation (ISF) . The strains were cultured aerobically on MTNA (mannitol, trimethoprim, nalidixic acid, amphotericin) medium without antibiotics  at 25°C for 24-48 h. Stock cultures were stored at −80°C in MicrobankTM beads (Pro-Lab Diagnostics, Canada).
DNA extraction, amplification and sequencing
Total genomic DNA was extracted according to the guanidium-thiocyanate-EDTA-sarkosyl method described by Pitcher et al.  which was adapted for Gram-positive bacteria by a pre-treatment with lysozyme (5 mg/μl lysozyme in TE buffer). Amplification and sequencing primers are listed in Table 2. The expected amplicons were generated with the Qiagen Taq DNA polymerase kit (supplemented with a Q-Solution) and GeneAmp® dNTP’s (Applied Biosystems, Belgium) according to the manufacturer specifications and with primers from Sigma Aldrich (Belgium). Amplicons were purified using the Nucleofast®96 PCR clean up membrane system (Macherey-Nagel, Germany). Sequencing PCR was performed in a total volume of 10 μl with 3 μl of a purified amplicon, 0.286 μl of BigDye™ mixture (Terminator Cycle Sequencing Kit version 3.1, Applied Biosystems), 1x sequencing buffer and 1.2 μM of each of the amplification primers listed in Table 2. The PCR program consisted of 30 cycles (96°C for 15 s, 35°C for 1 s, 60°C for 4 min). Subsequently, the sequencing products were purified using the BigDye XTerminator Kit (Applied Biosystems) and analyzed on a 3130xl Genetic Analyzer (Applied Biosystems).
In the frame of the European project QBOL (Quarantine Barcoding Of Life) we developed a gyrB barcode that was proven suitable to identify members of the genus Clavibacter at the subspecies level (http://www.q-bank.eu/) . Moreover, gyrB gene was used in MLST schemes developed to type Cmm strains [7, 33, 37]. DnaA sequence was shown a good taxonomic marker to identify and classify plant pathogenic bacteria such as Clavibacter, Xanthomonas and Ralstonia. The partial sequencing of dnaA was successfully used to study genetic diversity of non-pathogenic Clavibacter-like strains and to identify members of the genus Clavibacter (J. Zaluga, data unpublished). The gyrB and dnaA sequences were assembled with BioNumerics version 5.1 (Applied Maths, Belgium) and aligned using ClustalW . GyrB sequences and dnaA sequences were checked by amino acid translation with Transseq (http://www.ebi.ac.uk/Tools/emboss/transeq/) and presence of the GyrB and DnaA protein domain was confirmed with BlastP . DnaA and gyrB amplicons were 675 bp and 440 bp long (equal length was used for all strains), respectively. A phylogenetic tree was constructed on dnaA-gyrB concatenated sequence data with Molecular Evolutionary Genetics Analysis software (Mega 5.1) , using the Maximum Likelihood method with the Tamura-Nei model  and 1000 bootstrap replicates. The position of the sequenced gyrB and dnaA amplicons were checked by comparison to the reference Cmm genome sequence (AM711867). Newly generated gyrB and dnaA sequences have following accession numbers KC521547-521623 and have been deposited in NCBI database. Each unique sequence of a gene was assigned an allele number and the combination of allele numbers for each isolate defined the haplotype. Number of haplotypes, haplotype diversity and number of polymorphic sites were estimated for gyrB and dnaA genes using DnaSP version 5.0 . Percentages of polymorphic sites at the analyzed loci were calculated by dividing the number of polymorphic positions by the total length of the gene. The Discriminatory Power (D) was calculated using a discriminatory power calculator (http://insilico.ehu.es/mini_tools/discriminatory_power/index.php). The Discriminatory Power (D), as shown by Hunter can be expressed by the formula of Simpson’s index of diversity, which reads:
Where D is the index of discriminatory power, N the number of unrelated strains tested, S the number of different types, and xj the number of strains belonging to the jth type, assuming that strains will be classified into mutually exclusive categories. Thus, a D value of 1.0 would indicate that a typing method was able to distinguish each member of a strain population from all other members of that population. Conversely, an index of 0.0 would indicate that all members of a strain population were of an identical type. An index of 0.50 would mean that if one strain was chosen at random from a strain population, then there would be a 50% probability that the next strain chosen at random would be indistinguishable from the first .
Design of VNTR primers
The complete genome sequence of Clavibacter michiganensis subsp. michiganensis NCPPB 382 deposited under accession number AM711867 was screened for VNTR loci. Tandem Repeat Finder program (http://tandem.bu.edu)  was used to detect potential VNTR loci. Primer3 software  was used to design locus-specific amplifications and sequencing primers in regions flanking VNTR loci. Eight loci (Table 3) of 20 bp to 45 bp long tandem repeat (TR) units were selected. TRs longer than 20 bp were chosen to enable easier interpretation of results from an agarose gel. Primer pairs targeting single locus alleles were manually designed in the conserved regions to obtain amplicons of no more than 450 bp in length.
VNTR PCR amplification and sequencing
The PCR mixture had a total volume of 25 μl, containing 1 x PCR buffer (100 mM Tris–HCl, 15 mM MgCl2, 500 mM KCl [pH 8.3]) (Qiagen), dNTP’s 0.2 mM each, 0.6 μM of each primer, 0.5 U DNA Taq polymerase, and 50–60 ng template DNA. The PCR amplifications were performed under following conditions: 3 min denaturation step at 94˚C; 35 cycles of 94˚C for 1 min, annealing at 60˚C for 1 min, and extention at 72˚C for 1 min; and a final extension step at 72˚C for 10 min. Amplified products were run on a 2.5% Gel Pilot® Small Fragment Agarose (Qiagen) at 110 V for 2.5 hrs at 4°C using 25 bp size marker (Invitrogen), and visualized by ethidium bromide staining. PCR amplicons from one representative strain per different locus of a particular VNTR were sequenced using sequencing primers (Table 2) according to the sequencing protocol described above for gyrB and dnaA genes.
VNTR analysis and statistics
Product sizes were estimated and the exact number of repeats present was calculated using a derived allele-naming table, based on the number of repeats which could theoretically be present in a PCR product of a given size, allowing for extra flanking nucleotides and primer size. Theoretical number of repeats was confirmed subsequently by sequencing. Loci were named simply on the basis of the order in which they were found by the initial search. VNTR allele calls were analyzed in BioNumerics as ‘character’ data. Composite datasets were created for the eight Clav-VNTR loci. Distance trees were derived by clustering with the unweighted pair group method with arithmetic means (UPGMA), using ‘categorical’ character table values. All markers were given equal weight, irrespective of the number of repeats. The percentages in the dendrogram reflect the percentage of homology between the specific markers. Relatedness between the different haplotypes was investigated based on comparison of allelic profiles using the minimum spanning tree (MST) method from BioNumerics v 5.1. We used the classical criterium of one allelic mismatch to group haplotypes into clonal complexes. In order to assess the evolutionary relatedness between haplotypes the MLVA data was analyzed taking into account the number of repeat differences. The type strain LMG 7333T served as a reference and a starting point for calculations of the differences in other strains. For each VNTR locus the Hunter–Gaston and Simpson’s diversity indices were calculated using the VNTR diversity and confidence extractor software (V-DICE) available at the Health Protection Agency bioinformatics tools website (http://www.hpa-bioinformatics.org.uk/cgi-bin/DICI/DICI.pl) . Shannon-Wiener index of diversity was calculated using BioNumerics version 5.1.
Assessment of genetic diversity among Clavibacterstrains
In total, 62 strains representing the Clavibacter subspecies and non-pathogenic Clavibacter-like strains were included in this study. The identity of included Cmm strains was confirmed by analysis of the gyrB and dnaA gene sequences. The gene sequence analyses were performed on several related Clavibacter strains in order to study the genetic diversity in the genus Clavibacter. Phylogenetic analysis of two tested genes confirmed a clear separation of Clavibacter subspecies and a distinct position of non-pathogenic Clavibacter-like strains. Phylogenetic relationship between the Clavibacter subspecies and non-pathogenic Clavibacter-like strains was strongly supported by high bootstrap values (Figure 1). The number of polymorphic sites was 47 (10.7%) and 87 (12.9%), for gyrB and dnaA, respectively. It has to be noted that diversity among Cmm strains, especially among strains from recent Belgian outbreaks, was small which resulted in a limited number of clusters. Despite a low genetic diversity, a number of groups could be distinguished in a Cmm cluster (Figure 1). The largest cluster, containing Belgian strains from recent outbreaks and two French strains from 2010 (GBBC 1077 and GBBC 1078), was separated from the Cmm strains isolated previously in Belgium (Figure 1). Furthermore, strains originating from the same location mostly grouped together, such as French strains GBBC 1079, GBBC 1080 and PD 5719. However, based on the concatenated Maximum Likelihood tree of gyrB and dnaA no clear geographical separation among Cmm strains could be demonstrated. In gyrB and dnaA trees (data not shown) and in a concatenated tree Clavibacter subspecies are separated from each other and from non-pathogenic strains which suggests that they present the same phylogenetic information (Figure 1).
Development and implementation of MLVA
In parallel with the sequence analysis Cmm strains were investigated with MLVA. Fifty eight VNTR loci were identified in the genome of Cmm NCPPB 382. Thirty one of them were tested on a set of eight genetically diverse Cmm strains originating from geographically spread locations (Table 1). Subsequently, eight loci that were successfully amplified and showed to be polymorphic in the tested subset of strains were selected for further analysis. Successful amplification was obtained in all tested Cmm strains. Regarding the non-pathogenic, seed-borne Clavibacter-like strains the results varied from no amplification for Clav-VNTR5 or unspecific (more than one band, not expected product size) bands in Clav-VNTR26 (data not shown). Similar findings were observed for Clavibacter subspecies other than Cmm. In the cluster analysis, a total of 24 MLVA types were detected among 56 Cmm strains when the data from eight loci were combined, with allele numbers per locus ranging from two (Clav-VNTR22, Clav-VNTR26) to six (Clav-VNTR5) (Table 3, Figure 2). A large cluster, comprised of Cmm strains from recent Belgian outbreaks together with two French strains isolated in 2010, exhibited identical MLVA haplotypes. Strains from other countries formed mostly a separate branch or a cluster with two strains with an identical MLVA haplotype. No direct connection between strains from recent Belgian outbreaks of 2010–2012 and other Belgian strains included in this study could be observed. Remarkably, Belgian strains PD 5736 and GBBC 285, isolated in 1983 and 2008, respectively, showed the same MLVA haplotypes. In the concatenated tree of gyrB and dnaA these two Belgian strains clustered together among strains originating from other countries (Figure 1). Similar findings were observed for other two Belgian strains PD 1953 and GBBC 283, isolated in 1984 and 2002, respectively.
The discriminatory abilities of the MLVA technique was determined by calculating the discriminatory index (D) for 56 typed strains. MLVA differentiated 25 Cmm strains and showed a level of discrimination, with a D value of 0.8006. The discriminatory power of each VNTR was estimated by the number of alleles detected and the allele diversity. The number of different alleles ranged from two for Cmm-V22 and Cmm-V26 to six for Cmm-V5. Highest allelic diversities measured by Hunter–Gaston, Simpson’s and Shannon-Wiener diversity indices were 0.664; 0.652; 1.3377, respectively and were observed for the loci Clav-VNTR5 (Table 3). For the set under study, 27 different alleles of eight VNTR loci were observed. The relationship among the strains based on MLVA results is presented in a minimum spanning tree (MST) (Figure 3). The 56 Cmm strains were resolved into 24 types distributed into five complexes separating double locus variants (DLV). In addition, a large clonal group of Belgian strains from recent outbreaks (W), six singletons (S, T, Q, X, V, U) each represented by an isolate from a different country, and one separate group consisting of two strains (R) were detected (Table 1, Figure 3). Based on MLVA results, strains from Belgian outbreaks 2010–2012 were identical; no differences could be observed between strains originating from different years of isolation, tomato varieties or geographic locations in Belgium (Table 1, Figure 2, and Figure 3). To receive more information about evolutionary relatedness of strains from Belgium and France the MLVA data was analyzed taking into account the number of repeat differences (Additional file 1: Figure S1). Interestingly, Belgian strain PD 5737 and French strain PD 5749 clustered closer to ES2686.1 and CL01TF02 strains isolated in Spain during bacterial canker outbreak in 2002–2003. Moreover, these four strains showed to have a more similar MLVA haplotype to the group of strains from recent Belgian outbreaks 2010–2012.
Discussion and conclusion
Over the last few decades, bacterial canker has been frequently detected in tomato production areas, leading to substantial financial and economical losses. Only during the last three years several local outbreaks of Cmm were reported in Belgium. In some cases, reoccurring infections were detected in the primarily contaminated farms, suggesting a persistence of an initial infection source. Despite a quite frequent detection of tomato canker and wilting in Belgian tomato production areas there is little known about the genetic diversity of Cmm strains which hinders the correct conclusions about the probable sources of epidemics and transmission routes of Cmm.
This study is the first MLVA approach developed for efficient genotyping of Cmm strains. To date typing of Cmm strains was performed by RAPD-PCR , BOX-PCR [8, 48], AFLP , PFGE  and MLST . Despite the fact that some of these methods were found to have a good resolution most of them have limitations such as a poor interlaboratory portability or limited exchangeability of results that were generated on a specific machine or compared to an in-house database. Nowadays, fully sequenced genomes give a unique opportunity for a development of more robust and accurate typing methods such as MLVA. Its advantages, such as, high reproducibility, exchangeability of results and the possibility to add loci greatly facilitates epidemiological studies of economically important pathogens such as Cmm.
In this work, Clav-VNTR5 showed to be the most polymorphic loci with five different alleles and the highest HGDI of 0.664. Combined data from MLVA analysis of all eight investigated loci resulted in 25 different haplotypes and a discriminatory power of 0.8006. Cmm strains from the recent epidemics in Belgium in 2010–2012 showed identical MLVA haplotypes which suggests that a clonal population was responsible for these outbreaks. The presence of the same MLVA haplotypes of Cmm strains from 2011 and 2012 could mean that bacteria persisted in the used equipment, devices or soil and induced the outbreaks in the following years. Population of Belgian strains isolated from 2010–2011 is epidemiologically related to at least two French strains that exhibited the same MLVA haplotype. Moreover, based on minimum spanning tree, Belgian strains were found to be evolutionary related to the French strain PD 5749. When MLVA data was analyzed taking into account differences in the number of repeats it appeared that two French and two Spanish strains were found to have a similar MLVA haplotype to the group of Belgian strains from 2010–2012 suggesting that there might be a common origin of these strains (Additional file 1: Figure S1). It is worth mentioning that the strain ES 2686.1 isolated in Spain in 2002 was linked to outbreaks of Cmm in 2002–2007 in Canary Islands . Two French strains isolated in 2010 showed the same MLVA haplotype as strains from recent Belgian outbreaks which may imply that the contaminated material was spread also in France. Different MLVA patterns between strains from the recent Belgian outbreaks of 2010–2012 and Belgian strains isolated previously support our hypothesis about a novel introduction, presumably originating from a single lot of seeds or contaminated tomato seedlings. Remarkably, all Belgian Cmm strains from 2010–2012 (Table 1), were purchased from the same nursery.
In this study, VNTR loci were chosen to be longer than or equal to 20 bp to simplify the interpretation of the results from an agarose gel and to allow performing the analysis in standard laboratories not equipped in sophisticated tools (fragment analyzer or sequencer) required to analyze small (a few nucleotides) differences in an amplicon size. Shorter repeats are represented in a higher number of copies and are more likely to be polymorphic . However, many studies showed successful application of longer repeats which gave satisfactory resolution and discriminatory power [16, 50]. Moreover, in silico analysis of tandem repeats in the Cmm genome NCPPB 382 revealed only a few short repeats (6–8 bp) that had remarkably higher number of copies (around 10 copies).These microsatellite loci might be investigated in the future and combined with currently available MLVA scheme. MLVA can provide phylogenetic information even with a limited number of loci . MLVA assays are relatively robust [17, 52] but as any other technique they have their limitations. In MLVA, a need to develop a new set of loci for every species or serovar under investigation might be necessary. Moreover, some loci are ‘not stable’ and can ‘disappear’ from some strains or lineages what will result in an uninformative ‘zero’ allele .
VNTRs might possibly contribute to the genomic polymorphism and/or evolution. Comparative genomics of pathogenic Mycobacterium tuberculosis showed that a variation in size and number of repeats, located in coding regions, can result in a variable expression of surface-exposed proteins that play a role in pathogenicity . These changes could possibly help the pathogen to avoid the host immune response. Expansion or reduction of the number of tandem repeats can influence the expression, structure and activity of cellular proteins. Tandem repeats located within regulatory regions can result in a modification of gene expression at the transcriptional level . All tested Clav-VNTR loci were found in putative coding regions (Table 2). At least two of them were found within genes linked to processes taking place in a cell envelope (Clav-VNTR-13: putative NAD (FAD)-dependent dehydrogenase and Clav-VNTR 16: putative glycine/betaine ABC transporter). We could speculate that variability observed within these regions might possibly help bacteria to alternate the proteins of a cell envelope. However, more research has to be performed on the role of tandem repeat copy, and virulence in Cmm.
The genetic structure of the studied strains was assessed by the sequence analysis of two housekeeping genes, gyrB and dnaA, which were previously reported to be good molecular markers for studying populations of the genus Clavibacter[32, 38]. The phylogenetic position of Cmm strains was supported by high bootstrap values in a Maximum Likelihood tree. High similarity of Belgian strains from recent outbreaks was detected both, in a gene sequence analysis and by an MLVA typing method, supporting the hypothesis about their monomorphic nature. The percentages of polymorphic sites observed for the concatenated set of gyrB and dnaA genes (Table 4) was higher than the value obtained from five concatenated genes described in a recently published MLSA scheme of Clavibacter michiganensis subsp. michiganensis, (12 versus 8.8) . Based on these parameters the genes selected in this work can be applied in MLST studies to investigate highly similar Cmm populations.
In this study, MLVA was successfully applied to investigate a genetic relationship of Cmm strains from recent Belgian outbreaks. Its discriminatory power, measured by HGDI, was higher than these of each of the tested genes, gyrB and dnaA (Table 4). Our study has shown that MLVA analysis offers better discrimination of Cmm strains (HGDI = 0.8) than the typing method based on the concatenated tree of gyrB and dnaA (HGDI = 0.758) (Table 4). A significant advantage of the MLVA method is the excellent interlaboratory reproducibility  which makes this method well-suited for accurate and reproducible bacterial typing applicable in epidemiological studies of Clavibacter. MLVA, with its high discriminatory power to separate closely related strains, might be very useful for tracking sources of epidemic outbreaks as well as for investigating various haplotypes occurring during these outbreaks, as illustrated in the differentiation of Cmm strains. The technique is fast (results within one day), easy to perform, user-friendly, cost-effective compared to other typing techniques (e.g. AFLP) with an excellent reproducibility (intra- and interlaboratory). Additionally, data storage, comparison and exchange of the results are possible and easy. Moreover, the use of fluorescence-labeled primers enables multiplex PCR and subsequent analysis in a fragment analyzer. It is worth mentioning that the MLVA scheme, derived from in silico analysis of a complete genome sequence of Cmm, was experimentally confirmed to be accurate. It is consistent with previous findings demonstrated for Xanthomonas citri pv. citri and is advantageous over other experimentally tested techniques such as AFLP or IS-LM-PCR, where in vitro vs. in silico accuracy values of 75% and 87%, respectively, were reported .
The MLVA method, with eight novel VNTR loci identified within the genome of Cmm, demonstrated its applicability as a new tool for the molecular investigation of bacterial wilting and canker outbreaks.
In the future, additional VNTR loci and Clavibacter isolates might enable unraveling intrapopulation genetic variation and assessing the robustness of the method for investigating bacterial canker outbreaks on a global scale.
De León L, Siverio F, Lopez MM, Rodriguez A: Clavibacter michiganensis subsp. michiganensis, a seedborne tomato pathogen: healthy seeds are still the goal. Plant Dis. 2011, 95 (11): 1328-1338. 10.1094/PDIS-02-11-0091.
EPPO: Clavibacter michiganensis subsp. michiganensis Diagnostics, PM7/42. Bulletin OEPP/EPPO. 2005, 35: 273-283.
Strider DL: Bacterial canker of tomato caused by Corynebacterium michiganense: a literature review and bibliography. 1969, Technical Bulletin North Carolina Agricultural Experiment Station, 193-
Jahr H, Bahro R, Burger A, Ahlemeyer J, Eichenlaub R: Interactions between Clavibacter michiganensis and its host plants. Environ Microbiol. 1999, 1 (2): 113-118. 10.1046/j.1462-2920.1999.00011.x.
Lamichhane JR, Balestra GM, Varvaro L: Severe Outbreak of Bacterial Canker Caused by Clavibacter michiganensis subsp. michiganensis on Tomato in Central Italy. Plant Dis. 2011, 95 (2): 221-221.
De León L, Rodriguez A, Llop P, Lopez MM, Siverio F: Comparative study of genetic diversity of Clavibacter michiganensis subsp. michiganensis isolates from the Canary Islands by RAPD-PCR, BOX-PCR and AFLP. Plant Pathol. 2009, 58 (5): 862-871. 10.1111/j.1365-3059.2009.02117.x.
Milijašević-Marčić S, Gartemann KH, Frohwitter J, Eichenlaub R, Todorović B, Rekanović E, Potočnik I: Characterization of Clavibacter michiganensis subsp. michiganensis strains from recent outbreaks of bacterial wilt and canker in Serbia. Eur J Plant Pathol. 2012, 134: 697-711. 10.1007/s10658-012-0046-x.
Kawaguchi A, Tanina K, Inoue K: Molecular typing and spread of Clavibacter michiganensis subsp. michiganensis in greenhouses in Japan. Plant Pathol. 2010, 59 (1): 76-83. 10.1111/j.1365-3059.2009.02207.x.
Bella P, Ialacci G, Licciardello G, La Rosa R, Catara V: Characterization of atypical Clavibacter michiganensis subsp. michiganensis populations in greenhouse tomatoes in Italy. J Plant Pathol. 2012, 94 (3): 635-642.
Kleitman F, Barash I, Burger A, Iraki N, Falah Y, Sessa G, Weinthal D, Chalupowicz L, Gartemann K, Eichenlaub R: Characterization of a Clavibacter michiganensis subsp. michiganensis population in Israel. Eur J Plant Pathol. 2008, 121 (4): 463-475. 10.1007/s10658-007-9264-z.
Sahin F, Uslu H, Kotan R, Donmez M: Bacterial canker, caused by Clavibacter michiganensis subsp. michiganensis, on tomatoes in eastern Anatolia region of Turkey. Plant Pathol. 2002, 51 (3): 399-399. 10.1046/j.1365-3059.2002.00715.x.
Cangelosi GA, Freeman RJ, Lewis KN, Livingston-Rosanoff D, Shah KS, Milan SJ, Goldberg SV: Evaluation of a high-throughput repetitive-sequence-based PCR system for DNA fingerprinting of Mycobacterium tuberculosis and Mycobacterium avium complex strains. J Clin Microbiol. 2004, 42 (6): 2685-2693. 10.1128/JCM.42.6.2685-2693.2004.
Bosch T, de Neeling A, Schouls L, Zwaluw K, Kluytmans J, Grundmann H, Huijsdens X: PFGE diversity within the methicillin-resistant Staphylococcus aureus clonal lineage ST398. BMC Microbiol. 2010, 10 (1): 40-10.1186/1471-2180-10-40.
Van Belkum A: Tracing isolates of bacterial species by multilocus variable number of tandem repeat analysis (MLVA). FEMS Immunol Med Microbiol. 2006, 49 (1): 22-27.
Van Belkum A, Scherer S, Van Alphen L, Verbrugh H: Short-sequence DNA repeats in prokaryotic genomes. Microbiol Mol Biol Rev. 1998, 62 (2): 275-
Harth-Chu E, Espejo RT, Christen R, Guzman CA, Hofle MG: Multiple-locus variable-number tandem-repeat analysis for clonal identification of Vibrio parahaemolyticus isolates by using capillary electrophoresis. Appl Environ Microbiol. 2009, 75 (12): 4079-4088. 10.1128/AEM.02729-08.
Lindstedt BA, Heir E, Gjernes E, Vardund T, Kapperud G: DNA fingerprinting of Shiga-toxin producing Escherichia coli O157 based on multiple-locus variable-number tandem-repeats analysis (MLVA). Ann Clin Microbiol Antimicrob. 2003, 2 (1): 12-10.1186/1476-0711-2-12.
Danin-Poleg Y, Cohen LA, Gancz H, Broza YY, Goldshmidt H, Malul E, Valinsky L, Lerner L, Broza M, Kashi Y: Vibrio cholerae strain typing and phylogeny study based on simple sequence repeats. J Clin Microbiol. 2007, 45 (3): 736-746. 10.1128/JCM.01895-06.
Heymans R, Schouls LM, van der Heide HG: Schim van der Loeff MF, Bruisten SM: Multiple-locus variable-number tandem repeat analysis of Neisseria gonorrhoeae. J Clin Microbiol. 2011, 49 (1): 354-363. 10.1128/JCM.01059-10.
van Cuyck H, Pichon B, Leroy P, Granger-Farbos A, Underwood A, Soullie B, Koeck J-L: Multiple-locus variable-number tandem-repeat analysis of Streptococcus pneumoniae and comparison with multiple loci sequence typing. BMC Microbiol. 2012, 12 (1): 241-10.1186/1471-2180-12-241.
Skuce RA, McCorry TP, McCarroll JF, Roring SMM, Scott AN, Brittain D, Hughes SL, Hewinson RG, Neill SD: Discrimination of Mycobacterium tuberculosis complex bacteria using novel VNTR-PCR targets. Microbiology. 2002, 148 (2): 519-528.
Marsh JW, O’Leary MM, Shutt KA, Pasculle AW, Johnson S, Gerding DN, Muto CA, Harrison LH: Multilocus variable-number tandem-repeat analysis for investigation of Clostridium difficile transmission in hospitals. J Clin Microbiol. 2006, 44 (7): 2558-2566. 10.1128/JCM.02364-05.
Hidalgo Ă, Carvajal A, La T, Naharro G, Rubio P, Phillips ND, Hampson DJ: Multiple-locus variable-number tandem-repeat analysis of the swine dysentery pathogen, Brachyspira hyodysenteriae. J Clin Microbiol. 2010, 48 (8): 2859-2865. 10.1128/JCM.00348-10.
Le Fleche P, Hauck Y, Onteniente L, Prieur A, Denoeud F, Ramisse V, Sylvestre P, Benson G, Ramisse F, Vergnaud G: A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. 2001, 1 (1): 2-10.1186/1471-2180-1-2.
Li Y, Cui Y, Hauck Y, Platonov ME, Dai E, Song Y, Guo Z, Pourcel C, Dentovskaya SV, Anisimov AP: Genotyping and phylogenetic analysis of Yersinia pestis by MLVA: insights into the worldwide expansion of Central Asia plague foci. PLoS One. 2009, 4 (6): e6000-10.1371/journal.pone.0006000.
Zhao WJ, Chen HY, Zhu SF, Xia MX, Tan TW: One-step detection of Clavibacter michiganensis subsp. michiganensis in symptomless tomato seeds using a Taqman probe. J Plant Pathol. 2007, 89 (3): 349-351.
Gironde S, Manceau C: Housekeeping Gene Sequencing and Multilocus Variable-Number Tandem-Repeat Analysis To Identify Subpopulations within Pseudomonas syringae pv. maculicola and Pseudomonas syringae pv. tomato That Correlate with Host Specificity. Appl Environ Microbiol. 2012, 78 (9): 3266-3279. 10.1128/AEM.06655-11.
Coletta-Filho HD, Takita MA, De Souza AA, Aguilar-Vildoso CI, Machado MA: Differentiation of strains of Xylella fastidiosa by a variable number of tandem repeat analysis. Appl Environ Microbiol. 2001, 67 (9): 4091-4095. 10.1128/AEM.67.9.4091-4095.2001.
Wang DY, Hadj-Henni L, Thierry S, Arna P, Chermette R, Botterel F, Hadrich I, Makni F, Ayadi A, Ranque S: Simple and Highly Discriminatory VNTR-Based Multiplex PCR for Tracing Sources of Aspergillus flavus Isolates. PLoS One. 2012, 7 (9): e44204-10.1371/journal.pone.0044204.
Bergsma-Vlami M, Martin W, Koenraadt H, Teunissen H, Pothier J, Duffy B, van Doorn J: Molecular typing of Dutch isolates of Xanthomonas arboricola pv. pruni isolated from ornamental cherry laurel. J Plant Pathol. 2012, 94 (1): S1. 29-S21. 35.
Bui Thi Ngoc L, Vernire C, Jarne P, Brisse S, Guerin F, Boutry S, Gagnevin L, Pruvost O: From local surveys to global surveillance: three high-throughput genotyping methods for epidemiological monitoring of Xanthomonas citri pv. citri pathotypes. Appl Environ Microbiol. 2009, 75 (4): 1173-1184. 10.1128/AEM.02245-08.
Zaluga J, Heylen K, Van Hoorde K, Hoste B, Van Vaerenbergh J, Maes M, De Vos P: GyrB sequence analysis and MALDI-TOF MS as identification tools for plant pathogenic Clavibacter. Syst Appl Microbiol. 2011, 34 (6): 400-407. 10.1016/j.syapm.2011.05.001.
Jacques MA, Durand K, Orgeur G, Balidas S, Fricot C, Bonneau S, Quillévéré A, Audusseau C, Olivier V, Grimault V: Phylogenetic analysis and polyphasic characterization of Clavibacter michiganensis strains isolated from tomato seeds reveal that non-pathogenic strains are distinct from C. michiganensis subsp. michiganensis. Appl Environ Microbiol. 2012, 78 (23): 8388-8402. 10.1128/AEM.02158-12.
ISF: Methods for the detection of Clavibacter michiganensis ssp michiganensis on tomato seeds Version 4. 2011, http://www.worldseed.org/isf/ishi_vegetable.html,
Jansing H, Rudolph K: Physiological capabilities of Clavibacter michiganensis subsp. sepedonicus and development of a semi-selective medium. Zeitschrift Fur Pflanzenkrankheiten Und Pflanzenschutz-Journal of Plant Diseases and Protection. 1998, 105 (6): 590-601.
Pitcher D, Saunders N, Owen R: Rapid extraction of bacterial genomic DNA with guanidium thiocyanate. Lett Appl Microbiol. 1989, 8 (4): 151-156. 10.1111/j.1472-765X.1989.tb00262.x.
Waleron M, Waleron K, Kamasa J, Przewodowski W, Lojkowska E: Polymorphism analysis of housekeeping genes for identification and differentiation of Clavibacter michiganensis subspecies. Eur J Plant Pathol. 2011, 131 (2): 341-354. 10.1007/s10658-011-9812-4.
Schneider KL, Marrero G, Alvarez AM, Presting GG: Classification of plant associated bacteria using RIF, a computationally derived DNA marker. PLoS One. 2011, 6 (4): e18496-10.1371/journal.pone.0018496.
Thompson J, Higgins D, Gibson T: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-10.1093/nar/22.22.4673.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.
Hunter PR, Gaston MA: Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988, 26 (11): 2465-2466.
Gelfand Y, Rodriguez A, Benson G: TRDB - the tandem repeats database. Nucleic Acids Res. 2007, 35 (suppl 1): D80-D87.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132 (3): 365-386.
Simpson EH: Measurement of diversity. 1949, Nature: Nature
Nazari F, Niknam GR, Ghasemi A, Taghavi SM, Momeni H, Torabi S: An investigation on strains of Clavibacter michiganensis subsp. michiganensis in north and north west of Iran. J Phytopathol. 2007, 155 (9): 563-569. 10.1111/j.1439-0434.2007.01304.x.
Klevytska AM, Price LB, Schupp JM, Worsham PL, Wong J, Keim P: Identification and characterization of variable-number tandem repeats in the Yersinia pestis genome. J Clin Microbiol. 2001, 39 (9): 3179-3185. 10.1128/JCM.39.9.3179-3185.2001.
Sobral D, Schwarz S, Bergonier D, Brisabois A, Feßler AT, Gilbert FB, Kadlec K, Lebeau B, Loisy-Hamon F, Treilles M: High Throughput Multiple Locus Variable Number of Tandem Repeat Analysis (MLVA) of Staphylococcus aureus from Human. Animal and Food Sources. PLoS One. 2012, 7 (5): e33967-10.1371/journal.pone.0033967.
Call DR, Orfe L, Davis MA, Lafrentz S, Kang M-S: Impact of compounding error on strategies for subtyping pathogenic bacteria. Foodborne Pathog Dis. 2008, 5 (4): 505-516. 10.1089/fpd.2008.0097.
Gulati P, Varshney R, Virdi J: Multilocus variable number tandem repeat analysis as a tool to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A. J Appl Microbiol. 2009, 107 (3): 875-884. 10.1111/j.1365-2672.2009.04267.x.
Broschat S, Call D, Davis M, Meng D, Lockwood S, Ahmed R, Besser T: Improved identification of epidemiologically related strains of Salmonella enterica by use of a fusion algorithm based on pulsed-field gel electrophoresis and multiple-locus variable-number tandem-repeat analysis. J Clin Microbiol. 2010, 48 (11): 4072-4082. 10.1128/JCM.00659-10.
Domenech P, Barry C, Cole ST: Mycobacterium tuberculosis in the post-genomic age. Curr Opin Microbiol. 2001, 4 (1): 28-10.1016/S1369-5274(00)00160-0.
Pourcel C, Minandri F, Hauck Y, D’Arezzo S, Imperi F, Vergnaud G, Visca P: Identification of variable-number tandem-repeat (VNTR) sequences in Acinetobacter baumannii and interlaboratory validation of an optimized multiple-locus VNTR analysis typing scheme. J Clin Microbiol. 2011, 49 (2): 539-548. 10.1128/JCM.02003-10.
Kremer K, Arnold C, Cataldi A, Gutierrez MC, Haas WH, Panaiotov S, Skuce RA, Supply P, van der Zanden AGM, van Soolingen D: Discriminatory power and reproducibility of novel DNA typing methods for Mycobacterium tuberculosis complex strains. J Clin Microbiol. 2005, 43 (11): 5628-5638. 10.1128/JCM.43.11.5628-5638.2005.
We thank the PD, GBBC and BCCM/LMG collections and Ana Rodríguez Pérez (Spain) for providing necessary strains. This work was performed in the Seventh Framework Programme of project KBBE-2008-1-4-01 (QBOL) nr 226482 funded by the European Commission. Het Fonds Wetenschappelijk Onderzoek-Vlaanderen (FWO) is acknowledged for the postdoctoral fellowship of Pieter Stragier, and the Belgian NPPO (FAVV) for partially financing ILVO-research. We thank dr. Kim Heylen for her critical reading and valuable comments on the manuscript.
The other authors declare that they have no competing interests.
PS, MM, JVV and PDV conceived the study and participated in its design and coordination. JVV and PDV provided the bacterial culture collection for the study. JZ participated in the design of the study, carried out the molecular work, performed the data analysis and drafted the manuscript. PS coordinated the work and performed the statistical analysis. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Grouping of 56 Cmm strains using categorical values and the UPGMA (Unweighted-Pair Group Method with Arithmetic Mean) algorithm, generated with BioNumerics 5.1 software based on the number of repeats differences. Numbers in the Cmm-V2-26 columns indicate numbers of repeats differences. (DOCX 30 KB)
About this article
Cite this article
Zaluga, J., Stragier, P., Van Vaerenbergh, J. et al. Multilocus Variable-Number-Tandem-Repeats Analysis (MLVA) distinguishes a clonal complex of Clavibacter michiganensis subsp. michiganensisstrains isolated from recent outbreaks of bacterial wilt and canker in Belgium. BMC Microbiol 13, 126 (2013). https://doi.org/10.1186/1471-2180-13-126
- Bacterial typing
- Bacterial wilt and canker
- Plant pathogen