Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Methods

New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations

Jérémy Guinard, Anne Latreille, Fabien Guérin, Stéphane Poussier, Emmanuel Wicker
Harold L. Drake, Editor
Jérémy Guinard
aUniversité de la Réunion, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Saint Denis, Réunion, France
bCIRAD, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Pôle de Protection des Plantes, Saint Pierre, Réunion, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anne Latreille
bCIRAD, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Pôle de Protection des Plantes, Saint Pierre, Réunion, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fabien Guérin
aUniversité de la Réunion, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Saint Denis, Réunion, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stéphane Poussier
aUniversité de la Réunion, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Saint Denis, Réunion, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emmanuel Wicker
bCIRAD, UMR 53 Peuplements Végétaux et Bioagresseurs en Milieu Tropical (PVBMT), Pôle de Protection des Plantes, Saint Pierre, Réunion, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Harold L. Drake
University of Bayreuth
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.03095-16
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is considered one of the most harmful plant diseases in the world. Special attention should be paid to R. pseudosolanacearum phylotype I due to its large host range, its worldwide distribution, and its high evolutionary potential. So far, the molecular epidemiology and population genetics of this bacterium are poorly understood. Until now, the genetic structure of the RSSC has been analyzed on the worldwide and regional scales. Emerging questions regarding evolutionary forces in RSSC adaptation to hosts now require genetic markers that are able to monitor RSSC field populations. In this study, we aimed to evaluate the multilocus variable-number tandem-repeat analysis (MLVA) approach for its ability to discriminate genetically close phylotype I strains and for population genetics studies. We developed a new MLVA scheme (MLVA-7) allowing us to genotype 580 R. pseudosolanacearum phylotype I strains extracted from susceptible and resistant hosts and from different habitats (stem, soil, and rhizosphere). Based on specificity, polymorphism, and the amplification success rate, we selected seven fast-evolving variable-number tandem-repeat (VNTR) markers. The newly developed MLVA-7 scheme showed higher discriminatory power than the previously published MLVA-13 scheme when applied to collections sampled from the same location on different dates and to collections from different locations on very small scales. Our study provides a valuable tool for fine-scale monitoring and microevolution-related study of R. pseudosolanacearum phylotype I populations.

IMPORTANCE Understanding the evolutionary dynamics of adaptation of plant pathogens to new hosts or ecological niches has become a key point for the development of innovative disease management strategies, including durable resistance. Whereas the molecular mechanisms underlying virulence or pathogenicity changes have been studied thoroughly, the population genetics of plant pathogen adaptation remains an open, unexplored field, especially for plant-pathogenic bacteria. MLVA has become increasingly popular for epidemiosurveillance and molecular epidemiology studies of plant pathogens. However, this method has been used mostly for genotyping and identification on a regional or global scale. In this study, we developed a new MLVA scheme, targeting phylotype I of the soilborne Ralstonia solanacearum species complex (RSSC), specifically to address the bacterial population genetics on the field scale. Such a MLVA scheme, based on fast-evolving loci, may be a tool of choice for field experimental evolution and spatial genetics studies.

INTRODUCTION

Intensification of agriculture during the 20th century led to global changes in cropping practices and to intensification of rapid and large-scale exchanges. Thus, it increased risks for food security by contributing to significant changes in plant pathogen populations, such as emergence and/or reemergence of pathogens and selection of virulent pathogen populations (1). Among the different strategies to control plant pathogens, deploying resistant plant varieties appears to be the most promising solution (2). The implementation of efficient and durable disease management strategies requires deciphering of the ecological and genetic mechanisms that lead to the adaptation of a pathogen population to a new host or a new resistant cultivar (3). Thorough investigation of these mechanisms on the field scale requires the use of efficient ad hoc molecular markers.

To date, several markers have been developed through the use of DNA fingerprinting methods (restriction fragment length polymorphism [RFLP], amplified fragment length polymorphism [AFLP], and random amplified polymorphic DNA [RAPD] analyses) and housekeeping gene sequencing (multilocus sequence typing [MLST] and multilocus sequence analysis [MLSA]). Nevertheless, variable-number tandem-repeat (VNTR) markers, with their fast molecular clock (4), remain the most suitable tool for molecular epidemiology studies (5). VNTR are tandem repeats of sequences ranging from 5 to 100 bp; they have been used more and more for genotyping due to their portability between different laboratories, their reproducibility, and their high level of polymorphism on monomorphic bacteria (6). Multilocus VNTR analysis (MLVA) has been used widely in the medical field, especially to monitor epidemics of monomorphic bacteria, such as Bacillus anthracis (7), Yersinia pestis (8, 9), Mycobacterium tuberculosis (10, 11), Burkholderia pseudomallei (12, 13), and methicillin-resistant Staphylococcus aureus (14), and to assess their evolutionary history on different time and space scales (15, 16). Although there is a growing interest in the use of MLVA on plant-associated bacteria, only a few studies have been published, focusing on their population biology. These studies aimed at tracing sources of water contamination (Ralstonia solanacearum [17]), assessing the host-driven pathogen population structure (Pseudomonas syringae pv. maculicola and P. syringae pv. tomato [18], Xanthomonas axonopodis pv. manihotis [19], and R. solanacearum [20]), reconstructing the emergence history (Xanthomonas citri pv. citri [21–23]), evaluating the structure and evolution of pathogen populations (X. citri pv. citri [22] and Erwinia amylovora [24]), and global surveillance (Xanthomonas oryzae pv. oryzae and X. oryzae pv. oryzicola [25]). Studies investigating population structure on a small geographical scale have seldom been realized for plant-associated bacteria or for plant-associated oomycetes and fungi (22, 26–28).

Bacterial wilt caused by the R. solanacearum species complex (RSSC) (29, 30), a soilborne proteobacterium, is one of the most economically important bacterial plant diseases worldwide (31), mainly due to its wide geographical distribution and its large host range (250 host species in 54 botanical families [32]). The RSSC is considered a rather primitive vascular plant pathogen (33) and is composed of four phylotypes that correlate with their geographical origins (29). These phylotypes display specific evolutionary dynamics as estimated by their recombination rates and demographic histories (34). A splitting of the RSSC was proposed recently, reclassifying (i) phylotypes I and III into Ralstonia pseudosolanacearum, (ii) phylotype II into R. solanacearum, and (iii) phylotype IV into Ralstonia syzygii (35). The R. pseudosolanacearum phylotype I population structure presents molecular signatures of a recent and rapid demographic expansion (34), suggesting a high dissemination capacity. Moreover, phylotype I is highly recombinogenic (34, 36) and distributed worldwide (32), affects the largest number of hosts (37, 38), and contains all six pathoprofiles as described by Lebeau et al. (virulence patterns on a core collection of tomato, eggplant, and pepper plants [39–41]). Altogether, these features clearly indicate that phylotype I possesses a high evolutionary potential sensu McDonald and Linde (42). This phylotype is thus an interesting model for studying the evolutionary dynamics of plant-bacterium interactions. Although the molecular basis of its interactions with plants and the variability of its genome have been investigated widely (43–49), the population genetics and dynamics of this fascinating plant pathogen remain poorly understood and considered. Until now, only two studies have addressed the population genetic structure on a local scale. The first one was conducted in India by using RAPD markers (50), and the second was performed on Trinidad Island by using repetitive sequence-based PCR (REP-PCR) (51). Unfortunately, the low repeatability and reproducibility with RAPD markers make the results difficult to interpret, and REP-PCR portability between laboratories remains problematic (6). To date, three MLVA schemes have already been developed, to decipher RSSC population genetics on different scales (41), to investigate sources of water contamination by Ralstonia solanacearum (phylotype IIB) (17), and to characterize African R. pseudosolanacearum phylotype III strains on local and global scales (52). The MLVA phylotype I-specific scheme of N′Guessan et al. (20) (13 loci; MLVA-13) and the MLVA phylotype III-specific scheme of Ravelomanantsoa et al. (52) gave promising results in discriminating populations on scales ranging from cropping areas to continents (41).

In the present study, we first evaluated the MLVA-13 scheme for its discriminatory power on strains from the same field. Given that this scheme was poorly polymorphic on the field scale, we developed a new MLVA phylotype I-specific scheme and tested its use for conducting population genetics and microevolution studies on the field scale. More specifically, our objectives were (i) to identify new VNTR loci displaying potentially higher allelic diversity within phylotype I and III genomes, (ii) to evaluate their discriminatory power on collections sampled from the same location at different dates and collections sampled from different locations but on very small scales (single bulk soils and different plant locations), and (iii) to verify their suitability for assessing population differentiation and population evolution.

RESULTS

The MLVA-13 scheme is not adapted for analysis of field populations.We first used the previously developed MLVA-13 scheme (20) on the Vallon 2012 field collection, composed of 319 strains of R. pseudosolanacearum phylotype I. Within this collection, egl-mutS typing identified a single haplotype (egl-ST043, mutS-ST022, sequevar 31), identical to the sequevar 31 strain RUN0471, previously reported from Réunion Island (41). MLVA-13 identified only three haplotypes, which were most highly related to the reference haplotype H123, a singleton also corresponding to RUN0471. Only one locus (RS3L20) was polymorphic, with only three alleles detected. The genetic diversity detected was too low (Table 1) to properly perform population genetic analyses, leading us to develop a more discriminative MLVA scheme.

View this table:
  • View inline
  • View popup
TABLE 1

Comparison of genetic diversity revealed by the MLVA-13 and MLVA-7 schemes for descending scales, i.e., worldwide, regional (Ivory Coast), and field (Vallon 2012) scalesb

Tandem repeat screening leads to development of a new MLVA scheme composed of seven polymorphic loci (MLVA-7).TR searches using the Tandem Repeats Database (TRDB) and the Microorganisms Tandem Repeats Database (MTRDB) led us to select 24 TR loci: 11 localized on the chromosome (four from strain GMI1000, three from strain Y45, and four from strain CMR15) and 13 localized on the megaplasmid (five from GMI1000, three from Y45, and five from CMR15). Four of the 11 chromosomal loci were intergenic, whereas all 13 megaplasmid loci were intergenic. No significant (NS) difference in polymorphism was observed either between chromosomal and megaplasmid loci (P = 1.000) or between intergenic and intragenic loci (P = 0.625). Prescreening of these 24 loci on a 19-strain collection subset by high-resolution MetaPhor gel electrophoresis revealed only six loci to be polymorphic. These were between 6 and 9 bp long and repeated 5 to 21 times in the genome of origin (Table 2). We thus combined these six loci with the RS3L20 locus (from the MLVA-13 scheme) to define the new MLVA scheme, called MLVA-7. The interexperimental (between plates) variation of fragment size calling by capillary electrophoresis reached 1 bp when the GeneScan-500 LIZ size standard was used and allowed unambiguous assignation of alleles. The number of expected TR (in silico prediction) and the number of observed TR for reference strains GMI1000 (phylotype I) and CMR15 (phylotype III) were statistically similar to each other (chi-square test P value = 0.976 [NS] and 0.988 [NS], respectively).

View this table:
  • View inline
  • View popup
TABLE 2

Features of the VNTR loci composing the newly developed MLVA-7 scheme

The MLVA-7 scheme is specific to the RSSC and to phylotype I strains.All MLVA-7 loci were amplified only from RSSC strains, not from the other Ralstonia species (Ralstoniamannitolilytica, Ralstoniainsidiosa, and Ralstoniapickettii) strains, indicating that they were RSSC specific (Table 3). Amplification of the MLVA-7 loci showed clearly variable specificity for the different RSSC phylotypes of the worldwide collection. All seven loci were amplified from all R. pseudosolanacearum phylotype I strains. GMch3461 was the only phylotype I-specific locus, whereas all six other loci were amplified from the phylotype III strains, with rates ranging from 66.66% to 100%. The GMch0133 locus was also amplified from R. solanacearum (phylotype II) and R. syzygii (phylotype IV) strains (87.5% and 75%, respectively), whereas CMmp0233 was amplified from 25% of the phylotype II strains (Table 3).

View this table:
  • View inline
  • View popup
TABLE 3

Typeability of each MLVA-7 locus in the RSSC and other Ralstonia species, as estimated by the number (percentage) of strains from which the locus could be amplified within each RSSC phylotype and the Ralstonia outgroupa

MLVA-7 loci reveal genetic diversity from the field scale to the bulk soil scale.All MLVA-7 loci showed allelic diversity in field collections, ranging from 2 to 10 alleles per locus (see Table S2 in the supplemental material). The lowest diversity was found in the Pirogue 2007 soil sample (Martinique), as expected, with one or two alleles per locus (Table S2) and five multilocus haplotypes (Table S3). The highest diversity was found in the Vallon 2012 field collection (Table 4), with two to nine alleles per locus (Table S2) and 32 haplotypes (Table S3). Genetic diversity was thus highest in the Vallon 2012 collection, with decreasing levels in the Vallon 2009, Vallon 1999, and Pirogue 2007 collections (Table 4; see also Table S3). However, it is striking that the ratios between numbers of haplotypes and numbers of strains (G/N ratios) were similar for the Vallon 2012 and Pirogue 2007 collections (Table 4), despite the huge difference in sampling efforts.

MLVA-7 is more discriminant than MLVA-13 on both global and field scales.To compare the discriminatory power between MLVA-13 and the newly developed MLVA-7 scheme on the global (worldwide), regional (countrywide), and local (field-wide) scales, we considered the following different subsets of strains: (i) the worldwide phylotype I collection (n = 34), among which was defined (ii) a subset composed of representative phylotype I strains from the Ivory Coast (referred to here as the “country-scale” collection; n = 18), and (iii) the Vallon 2012 field collection (n = 319) (Table 1).

Considering the worldwide collection, the unbiased diversity index for the MLVA-7 scheme was almost 2 times higher than that for MLVA-13, though the G/N ratios were the same (0.91) (Table 1). Considering the country-scale collection (Ivory Coast strains [41]), a similar trend was observed, with MLVA-7 revealing a genetic diversity 2 times higher than that with the MLVA-13 scheme, as well as a higher G/N ratio (0.94 versus 0.83). For the Vallon 2012 field collection, the MLVA-7 scheme revealed a genetic diversity 23 times higher than that with the MLVA-13 scheme. The G/N ratio followed the same trend (Table 1).

MLVA-7 is congruent with MLVA-13 on the global scale, giving geographically relevant clusterings.We evaluated the congruence of both MLVA schemes toward the worldwide collection of R. pseudosolanacearum phylotype I strains (n = 34) by using two approaches: Mantel tests on distance matrices and comparison of clusterings obtained via discriminant analysis of principal components (DAPC).

Distance matrices (Euclidean and Manhattan distances) were highly correlated, with respective Mantel correlation coefficients of 0.666 and 0.640 (P = 0.001 for both Euclidean and Manhattan matrices), indicating that MLVA-13 and MLVA-7 data were significantly congruent.

DAPC was performed on the worldwide collection of phylotype I strains (n = 34) genotyped by both MLVA-13 and MLVA-7. DAPC clusters were named by giving the MLVA scheme (“DAPC13.” for data from the MLVA-13 typing scheme and “DAPC7.” for data from the MLVA-7 scheme) followed by the number of the cluster (cl01 to cl14). In both DAPC analyses, all strains were assigned to a single cluster (Fig. S1), with one exception (RUN1533 with MLVA-7). Among the 14 clusters identified in both data sets, 10 DAPC clusters obtained with the MLVA-7 scheme each perfectly matched a single MLVA-13 cluster (Table 5). DAPC7.03 was split between DAPC13.13 and DAPC13.14, DAPC7.09 was split between DAPC13.9 and DAPC13.10, and DAPC7.13 was split between DAPC13.03 and DAPC13.11 (Fig. 1). Strains RUN0044 (Australia), RUN0215 and RUN0119 (Cameroon), and RUN1533 and RUN1878 (Ivory Coast) were thus ranked differently by MLVA-13 and MLVA-7 (Fig. S1). We also checked for correspondence between the clonal complexes (CC) previously identified with MLVA-13 (20) and the DAPC clusters. Most of the DAPC13 clusters corresponded to a single CC, except for DAPC13.03, which grouped strains of CC02 and CC05. The correspondence between CC and DAPC7 clusters was a little lower, since two clusters corresponded to more than one CC: (i) DAPC7.09 (CC02 and CC10) and (ii) DAPC07.13 (CC02 and CC03). CC02 was split into three clusters with both MLVA-13 and MLVA-7 (Fig. 1); CC07 and CC10 were each split into two DAPC7 clusters, whereas CC11 was split into two DAPC13 clusters. DAPC clusters were globally correlated with geographical origins for both data sets, even though 4 of 14 clusters gathered strains from two or three countries (Fig. S3).

View this table:
  • View inline
  • View popup
TABLE 5

Correspondence between clusters identified by DAPC within the MLVA-13 data set and within the MLVA-7 data seta

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Correspondence between clonal complexes (defined in reference 41) and DAPC clusters inferred by MLVA-13 (right) and MLVA-7 (left). CC, clonal complex; S, singleton. The figure was created using the table.value function in the R adegenet package (87).

A DAPC analysis performed on the MLVA-7 data set obtained from the Vallon 2012 (Réunion Island) (n = 319), Pirogue 2007 (Martinique) (n = 38), and worldwide (n = 37) collections clustered Réunion Island strains into several clusters, together with strains from different African countries (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Scatterplot representing DAPC clusters for the Vallon 2012, Pirogue 2007 (Martinique), and worldwide collections, inferred by use of the MLVA-7 scheme. Réunion Island-specific clusters and the Martinique cluster (box 4) are associated with negative and positive values, respectively, on the x axis. Some clusters are shared across countries, i.e., Ivory Coast and Réunion Island (boxes 2 and 3); Martinique and French Guiana (box 4); Cameroon, Ivory Coast, Réunion Island, and Taiwan (box 7); Cameroon, Ivory Coast, South Africa, and Réunion Island (box 8); and Cameroon, French Guiana, Ivory Coast, Réunion Island, and South Africa (box 9). Numbers on axes indicate respective inertia rates. The figure was created using the find.cluster and scatterplot functions in the R package adegenet (87).

Altogether, these results indicate that MLVA-7 gave clusterings similar to those obtained with MLVA-13 on the global scale, and the geographical relevancy of these clusters seemed equivalent between both MLVA schemes. Since several clusters had mixed origins (Fig. 2), the usefulness of such an MLVA scheme for biogeographic studies looks questionable.

The MLVA-7 scheme allows for identification of epidemiological relationships within and between field collections.The epidemiological relationships between the 78 haplotypes identified within the entire phylotype I collection (Table S3) were visualized using a minimum spanning tree (Fig. 3). Two clonal complexes were identified. The main clonal complex (Fig. 2), whose founder (H018) was specific to the Vallon 2012 collection, consisted mostly of Vallon 2012 haplotypes, Vallon 2009 haplotypes, Vallon 1999 haplotypes, and three worldwide haplotypes: H58 (Ivory Coast), H20 (Ivory Coast and Vallon 2012), and H21 (Kenya/Ivory Coast). The other small clonal complex gathered all five Pirogue collection-specific haplotypes (Martinique) (shown in red in Fig. 3) and two haplotypes (H035 and H039) that were specific to French Guiana. Double- and triple-haplotype clonal complexes and singletons were mostly composed of worldwide haplotypes, as well as six Vallon 2012 haplotypes and one Vallon 2009 haplotype.

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Haplotype network built from the MLVA-7 data sets for different collections. The Vallon 1999 collection is represented in blue, the Vallon 2009 collection is represented in green, the Vallon 2012 collection is represented in black, the Pirogue collection (Martinique) is represented in red, and the worldwide collection is represented in gray. Networks link single-locus variants (SLVs) in blue, double-locus variants (DLVs) in green, and triple-locus variants (TLVs) in red. Labels within circles indicate haplotype numbers. (A) Global view. (B) Detailed view of clonal complex 01 (CC01). The founder of CC01 is circled in yellow (H18), whereas the founders of subgroups within CC are circled in light green. This figure was built using the goeBURST algorithm implemented in PHYLOViZ 1.1.

The MLVA-7 scheme allows monitoring of the temporal evolution and genetic structure of a field population (“Vallon” collections).To investigate the temporal evolution of the genetic diversity and structure of a field population, we considered the three “Vallon” collections only. Ten DAPC clusters were inferred, among which two (cl07 and cl08) were common to all three collections (“universal” clusters) and six (cl01, cl02, cl03, cl04, cl09, and cl010) were specific to the 2012 sampling (Fig. 4). DAPC strongly differentiated the universal clusters from four of the Vallon 2012-specific clusters (cl01, cl02, cl04, and cl10) and from cluster 05, composed of Vallon 2009 and Vallon 2012 strains. Interestingly, cluster 05 was first detected within the Vallon 2009 collection and then increased in size within the Vallon 2012 collection. Since the Vallon 1999 and Vallon 2009 collections were sampled from plant stems only, we hypothesized that the Vallon 2012-specific clusters may have been isolated from a specific habitat (rhizosphere or soil). However, there was no correlation between DAPC clusters and the habitat of isolation (data not shown).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

DAPC clusters and plot inferred for the “Vallon” collections. Strains in the Vallon 1999 (n = 96), Vallon 2009 (n = 90), and Vallon 2012 (n = 319) collections were genotyped using the MLVA-7 scheme (left), and their distribution within the three collections was determined using the table.value tool (right). Among the 10 clusters, five (cl01, cl02, cl07, cl08, and cl09) are specific to the Vallon 2012 collection, whereas three (cl03, cl04, and cl10) are common to all three collections. Numbers on axes indicate respective inertia rates. The figure was obtained using the find.cluster, scatterplot, and table.value functions in the R package adegenet (89). Inf, DAPC clusters inferred from the data set mentioned.

MLVA-7 allows inference of the reproductive mode of R. pseudosolanacearum phylotype I within field collections as estimated by multilocus LD.The Pirogue 2007 (Martinique) collection, though collected from a very small amount of soil, displayed a significant linkage disequilibrium (LD), suggesting a clonal reproductive mode (Table 6). For the Réunion Island “Vallon” collections, the results contrasted: the Vallon 1999 collection appeared to be nearly recombining, while the Vallon 2009 and Vallon 2012 collections appeared to be highly clonal.

View this table:
  • View inline
  • View popup
TABLE 6

Estimates of multilocus linkage disequilibriumb

Estimation of a TR locus mutation model.In order to estimate the mutation model followed by our molecular markers, we looked at the locus variation along the haplotype network of the minimum spanning tree, specifically into the TR number variation across single-locus variants (SLV) (Fig. 3). Three loci (GMch0133, GMch3461, and RS3L20) clearly followed the stepwise mutation model (more than 80% of −1/+1 TR variation), while CMmp0131 gave a majority (55.6%) of single-TR variants (STRV) but also 44.44% 2-TR variants. Two other loci gave a majority of STRV, but at much lower rates. GMch0754 indeed gave 40% STRV, but it also gave 2-, 5-, and 6-TR variants, at a rate of 20% each. GMmp0266 gave 26.09% STRV, as well as 2-TR (8.70%), 4-TR (13.04%), 5-TR (13.04%), and 7-, 8-, 9-, 10-, and 15-TR variants (at rates of 4.35, 8.70, 4.35, 17.39, and 4.35%). Finally, the CMmp0233 locus gave only 3- and 4-TR variants (50% for both).

DISCUSSION

The proposed new MLVA scheme is the fourth one developed for the RSSC, after the initial work of N′Guessan et al. (including the phylotype I-specific MLVA-13 scheme) (20), the phylotype IIB-1-specific scheme of Parkinson et al. (17), and the phylotype III-specific scheme of Ravelomanantsoa et al. (52). These three schemes, used on regional scales, allowed the inference of possible dissemination routes of phylotype I within the Ivory Coast and Africa (41), addressed the global surveillance of phylotype III (52), and identified sources of contamination of the Thames by phylotype IIB-1 (17). These studies improved the comprehension of RSSC epidemiology, its spread, its structuration on different scales, and its colonization pattern. However, all these studies focused only on large scales and did not pay particular attention to how microevolution acted on the field scale. Until now, field-scale studies focusing on RSSC genetic diversity were achieved only with markers with poor interlaboratory portability (RAPD [50], REP-PCR [51, 53], and AFLP [54, 55]).

In this study, we decided to develop a new MLVA scheme, as VNTR loci have already proved their use in field-scale microevolution studies of plant-pathogenic bacterial epidemics, e.g., for local populations of Xanthomonas citri pv. citri from different citrus orchards in Vietnam (22, 27) and in Mali and Burkina Faso (23) as well as local populations of X. oryzae pv. oryzae and X. oryzae pv. oryzicola (25).

Our primary goal in developing this new MLVA scheme was to differentiate genetically closely related strains of R. pseudosolanacearum phylotype I from a unique field, and even from a unique bulk soil. We indeed demonstrated that the first phylotype I-specific MLVA-13 scheme (20) was not powerful enough on the field scale, with only one polymorphic locus (CMmp0985; also named RS3L20). We thus identified six new, polymorphic, fast-evolving loci of 6 to 9 nucleotides that were repeated 7 to 22 times in the genomes of origin, thus matching the rules previously enounced for fast-evolving loci (20). The search for loci that were both potentially fast evolving and polymorphic had a relatively low success rate, since only six loci were retained among 24 selected. Our VNTR loci were specifically present in RSSC genomes and were even specific to phylotype I as well as phylotype III, to a lesser extent. These findings, which confirmed earlier reports (17, 20), tend to support the proposed concept of three genomic species within the RSSC: the first (R. pseudosolanacearum) gathering phylotypes I and III, the second (R. solanacearum) corresponding to phylotype II, and the third (R. syzygii) corresponding to phylotype IV (35, 43, 56, 57).

Our second goal was to assess the discriminatory power of the newly developed MLVA scheme on the field scale. The new MLVA-7 scheme is clearly adapted to monitoring the evolution of strains sampled from the same location at different dates and from different fields on very small scales (either from single bulk soils or from different plant locations). On the field scale, the discriminatory power of MLVA-7 is much higher than that of MLVA-13 (20). Genetic diversity could be revealed on a scale as small as a single soil sample (5 g) collected in a banana fallow in Martinique; this point may open research prospects for study of bacterial microevolution in soil habitats. The MLVA-7 scheme is also discriminant enough to reveal evolutionary patterns of a bacterial population sampled from the same field at different dates (“Vallon” collections). Interestingly, our markers could distinguish between haplotypes prevalent over 14-year and 4-year periods (from 1999 to 2012 and from 2009 to 2012, respectively) and haplotypes specific to one or another sampling date. This point also paves the way for future experimental evolution assays performed on the field scale.

The third goal of this study was to evaluate the suitability of our markers for assessing population genetics and evolution studies. We demonstrated that clustering inferred by use of MLVA-7 on a worldwide collection of isolates was globally congruent with clustering inferred from the previously published MLVA-13 scheme (20). MLVA-7 loci also identified the specificity of African phylotype I populations toward Asian ones and the close relatedness between some Réunion Island strains and some Ivorian strains, previously established by multilocus sequence analysis (C. A. N′Guessan and E. Wicker, unpublished data). The geographical coherence of our markers was good, as we did not identify common allelic profiles between field collections and the worldwide collection (with one exception), though we realize that these markers are of limited value for biogeographic studies. This geographical coherence looked better than that of phylotype IIB-1 markers, which identified common allelic profiles in English and Swedish strains (17). Finally, we demonstrated that these markers could be used to decipher the evolution of field populations over time, even with different sampling efforts. We were able to distinguish different clusters in a single field and their evolution, as previously evoked. We could also address the reproductive mode of R. pseudosolanacearum phylotype I field populations, distinguishing some nearly recombining populations from more clonal populations. Such evolutionary patterns need to be investigated further in collections generated by similar sampling efforts and with complementary approaches, including spatial genetics.

The genetic diversity estimated for our R. pseudosolanacearum phylotype I field collections (0.101 < HE < 0.298) appears to be low compared to values reported for plant-pathogenic bacteria, such as E. amylovora (HE = 0.50) (24), X. citri pv. citri (HE = 0.77) (22, 27) and R. solanacearum (HE = 0.595) (41). It is important that all these studies analyzed the genetic diversity in areas that were larger than ours (continent, country, and regional scales, respectively). Interestingly, our values were higher than that reported for an R. pseudosolanacearum phylotype III field population sampled in Cameroon (HE = 0.03; G/N = 0.235) (52). An explanation for the low genetic diversity detected, apart from a possible sampling bias effect, might be that the RSSC is capable of surviving in plants and in soil in a viable but not cultivable (VBNC) state (58). Therefore, this means that the haplotypes identified might not reflect the total variability in the field and that we detected only haplotypes specifically able to grow on an agar medium.

The mutation model for our MLVA-7 loci has been determined. It has been reported that VNTR loci evolve following a stepwise mutation model (SMM) or another model, e.g., an infinite allele model (IAM) (59), a two-phase model (TPM) (60), or a generalized stepwise mutation model (GSM), or k-allele model (KAM) (61). We showed that within the MLVA-7 scheme, six loci strictly or predominantly follow SMM, whereas one locus does not follow this model. This finding is important information for future population genetics analyses, notably the assessment of population differentiation. It also illustrates that mutation dynamics of VNTR loci can be different depending on the loci and on the species studied (62, 63). SMM is the predominant mutational model for some plant-pathogenic bacteria (24, 27, 64), but data remain scarce on the actual mutation model at play and the possible variations around this model within the Proteobacteria phylum.

In conclusion, this study provides an adapted tool for characterizing the genetic diversity and structure of R. pseudosolanacearum phylotype I on fine scales, from fields to bulk soil. This MLVA-7 scheme opens research prospects toward an understanding of the evolutionary mechanisms underlying population structure and adaptation as well as plant selective pressure on the field scale. This will thus allow characterization of the impacts of resistant or susceptible plant accession and of crop successions on the spatial genetic structure of soilborne bacterial plant pathogens. Such tools will allow unraveling of the bacterial population dynamics on the plant and field scales, help us to understand plant resistance bypass processes, and pave the way to optimization of innovative disease control strategies, such as the use of resistant/susceptible varietal mixtures.

MATERIALS AND METHODS

Bacterial strains.A total of 601 strains—including 598 RSSC strains—were used in this study (Table 7; see Table S1 in the supplemental material). The species specificity of the MLVA markers was checked with DNAs from Ralstonia pickettii (LMG5942T), R. mannitolilytica (LMG6866T), and R. insidiosa (LMG21421T) (Table S1). The phylotype specificity and phylogenetic coherence of our markers were tested on 55 RSSC strains of different phylotypes representative of the worldwide diversity (Table S1). This collection is referred to here as the “worldwide collection.” The suitability of our markers for discriminating populations on the field and soil scales was assessed on (i) three collections sampled from the same R. pseudosolanacearum phylotype I-infected field on Réunion Island at different dates (“Vallon” collections) and (ii) a collection extracted from 5 g of soil sampled at one single location point (“Pirogue” collection) in another R. pseudosolanacearum phylotype I-infected field in Martinique (Table 1; Table S1).

View this table:
  • View inline
  • View popup
TABLE 7

Summary of the collections of RSSC strains

Strains of the Vallon 1999 collection were isolated as previously described by Poussier (54). For the other collections, bacterial extractions were performed on 5-cm stem segments following a previously described protocol (40). For bacterial extraction from rhizospheres, the complete root systems were soaked in 200 to 500 ml of Van Elsas and Smalla buffer (VES buffer) (65) and shaken overnight at 120 rpm and 28°C. For bacterial extraction from soil, 20 g of fresh soil (0- to 20-cm top layer) was soaked in 100 ml of VES buffer, rotated for 30 min at 100 rpm, and then left to settle for 5 min. Series of 1-ml aliquots were sampled and diluted down to 10−4. Afterwards, macerates (50 μl) were streaked either on Imazaki medium (rhizosphere and soil) (66) or on modified Sequeira medium (67). Plates were then incubated at 28°C for 2 to 3 days. Isolates presenting the typical RSSC morphology—mucoid, “bird's eye”-shaped beige colonies with formazan coloration—were finally plated on Kelman's triphenyltetrazolium chloride (TZC) medium (68). Finally, bacteria were stored on microbeads (Microbank ProLab Diagnostic, Neston, Wirral, United Kingdom) or in 20% glycerol at −80°C.

Molecular typing and DNA extraction.Molecular identification of the RSSC and determination of the phylotype were performed using multiplex PCR (mxPCR) on bacterial suspensions (∼108 cells ml−1 in high-pressure liquid chromatography [HPLC]-grade water) as previously described (69, 70). PCR products were then subjected to 2% agarose electrophoresis and visualized under UV light after ethidium bromide staining (5 μg ml−1). For isolates identified as belonging to the RSSC, bacterial suspensions (∼5 × 108 to 1 × 109 CFU ml−1) were washed with 1 ml of 0.5 M NaCl and centrifuged at 12,000 × g. The pellets were subjected to Wizard genomic DNA purification lysis buffer (Promega, Fitchburg, WI), and DNA was extracted following the protocol provided by the manufacturer. The DNA concentration and quality were checked using a NanoDrop ND-8000 device (NanoDrop Technologies Inc.). For genotyping, DNA was diluted to 2 ng μl−1 in HPLC-grade water and stored at −30°C.

Sequence-based molecular typing (egl and mutS).Multilocus sequence typing based on partial sequences of DNA mismatch repair (mutS) and virulence-associated endoglucanase (egl) genes has been a reference method for RSSC molecular diagnostics (40, 69, 71, 72). egl sequevars and mutS sequence types were determined using the Endo-F/Endo-R and mutS-RsF.1570/mutSRsR1926 primer pairs, respectively. The PCR conditions and the sequencing procedure were described previously (40). Sequences were assembled with Geneious v. 6.0.3 (Biomatters Ltd., Auckland, New Zealand), aligned with the algorithm ClustalW (73) implemented in MEGA5 (74), and then trimmed as explained previously (34). egl-based sequevars and STs (eSTs) and mutS-based STs (mSTs) were assigned using egl and mutS reference sequences (34) within DnaSP, using the “generate haplotype list” command (75).

Genotyping with the MLVA-13 scheme (20).To assess whether MLVA-13, the first phylotype I-specific MLVA scheme, developed by N′Guessan et al. (20), was suitable for analysis on the field scale, 319 phylotype I strains from one field (Vallon 2012 collection) were genotyped using this scheme. Multiplexed primer mixes and the PCR conditions were described previously (20). Capillary electrophoreses were performed in an ABI Prism 3130XL sequencer, and the electropherograms generated were analyzed with GeneMapper v. 4.0 (Applied Biosystems, Life Technologies, Carlsbad, CA).

Development of a new MLVA scheme (MLVA-7).Three fully sequenced and annotated genomes (phylotype I strains GMI1000 [46] and Y45 [49] and phylotype III strain CMR15 [43]) that are publicly available via the MaGe interface (RalstoniaScope [76]) were screened for detection of tandem repeats using the Microorganisms Tandem Repeats Database (MTRDB) (77, 78) and the Tandem Repeats Database (TRDB) (79). It has been reported that the number of alleles correlates with the number of TR in the genome and conversely correlates with the pattern size (80), but so far, no general rules have been published. From the literature, expert personal communications (O. Pruvost, personal communication), and the development of the first RSSC MLVA scheme (20), we observed that the highest level of polymorphism was observed for loci with more than 10 repeats and pattern sizes ranging from 5 to 9 bp (20). As a consequence, after a preliminary search, we set the following criteria for TR detection in the MTRDB: total length in the range of 50 to 600 bp, TR length in the range of 5 to 9 bp, number of copies in the range of 9 to 50, matching percentage above 80%, and the default setting for other parameters. TR detection in the TRDB was performed using default parameters (alignment weights at 2 [match], 7 [mismatch], and 7 [indels]) (see the detailed explanation at the TRDB website [79]). PCR primers (18 to 27 bp) were designed within the 500-bp flanking sequences of the TR by using Geneious software v. 6.0.3 (Biomatters Ltd.). Primer sequences were then subjected to BLAST searches with the BLASTN program (81), available on the MaGe website (76), to determine the locations of loci and to check the specificity of each pattern to a single locus per genome. This first step led to the identification of 24 loci (3 detected by both platforms, 5 detected only by the MTRDB, and 16 detected only by the TRDB). These loci were further screened for amplification ability, repeatability, and polymorphism on a 19-strain collection subset (Table S1) including (i) 11 Vallon 2012 strains, selected on the basis of their allelic profiles according to the MLVA-13 scheme, their locations in the field, their hosts (tomato or susceptible or resistant eggplant), and their sampling habitats (stem, rhizosphere, or soil), and (ii) 8 worldwide strains (41), selected on the basis of the clonal complex to which they belong and their geographic origin. Simplex PCRs were performed as previously detailed (20), followed by electrophoresis migration in high-resolution Metaphor agarose (3% Metaphor and 1% SeaKem agarose; Lonza Ltd., Basel, Switzerland). Based on these results, six loci and primer pairs were selected to develop the new MLVA scheme (Table 2). As we found that CMmp0985 (RS3L20) from the MLVA-13 scheme (20) was the only polymorphic locus in the Vallon 2012 collection, this locus was added to the six other loci to compose the new MLVA-7 scheme. Primers were multiplexed into two pools of four and three pairs (Table 2), labeled on the forward primers with different fluorescent dyes, i.e., 6-carboxyfluorescein (FAM) (blue), VIC (green), NED (black/yellow), and PET (red) (Applied Biosystems). Primer concentrations and annealing temperatures were optimized for amplification of loci on four reference strains according to the manufacturer's instructions (Applied Biosystems). Capillary electrophoresis was performed as mentioned above.

VNTR analysis and performance measures.Peak data were analyzed using GeneMapper v. 4.0 (Applied Biosystems) and automatically converted to fragment sizes (in base pairs). The analyses described below were realized by using the amplification sizes. When needed, amplicon sizes were converted to repeat numbers by use of the following formula: observed size = extra TR sequence length + (TR size × repeat number). The number of TR was rounded up to the nearest integer, as previously recommended (82). The reproducibility of the amplification of our markers was assessed using GMI1000, the reference phylotype I strain (46), as a positive control for each genotyping run, and also by twice amplifying the phylotype I worldwide collection from independent DNA extractions. Allelic profiles were identical across the different replicates.

Analysis of genetic data.Haplotypes (defined as allelic profiles), Nei's unbiased diversity index (HE or Hub) (83), the total number of alleles (Na), and their frequencies were determined using GenAlex software v. 6.5 (84) (Table 7). To assess the correlation between locus location (chromosome versus megaplasmid and intergenic versus intragenic) and locus polymorphism, allele numbers of loci were compared statistically by the Fisher exact test under the null hypothesis (H0) that a locus has equal probabilities of being polymorphic for any position on the gene.

View this table:
  • View inline
  • View popup
TABLE 4

Summary statistics for the different collections of R. pseudosolanacearum phylotype I strains, based on MLVA-7 locia

Haplotype networks were achieved with the globally optimized eBURST analysis (goeBURST) implemented in PHYLOViZ v. 1.0 (85). These analyses allowed us to highlight different clonal complexes, defined as groups of SLVs. It also allowed us to visualize allelic profiles diverging by one locus (SLVs) or several loci (double-locus variants [DLVs], three-locus variants [TLVs], etc.). Using the eBURST algorithm, STs are clustered based on their numbers of SLVs, DLVs, and TLVs and occurrence frequencies. The group founder is first connected to all its SLVs; then, each of those, following the ST ordering, is connected to its SLVs not yet present in the tree. This process is repeated until each element of the disjoint graph is present in the tree. The second step consists of a local optimization of the tree. The goeBURST algorithm provides a globally optimal solution and corrects links that do not strictly follow the eBURST rules.

To assess the congruence of MLVA-13 and MLVA-7 clusterings, distance matrices were calculated and compared by the Mantel test, using the CADM.post function of the R package ape 3.1-4, with 9,999 permutations (86). The Mantel correlation coefficients were computed on rank-transformed distances among the distance matrices.

The population structure of R. pseudosolanacearum phylotype I (the worldwide strain collection and the three “Vallon” collections) was assessed using DAPC with the adegenet 1.4.2. package (87) implemented in R software v. 3.3.1. This clustering method, free of any assumption related to a population genetics model (e.g., Hardy-Weinberg equilibrium or an absence of linkage disequilibrium) (6), provides an efficient description of genetic clusters by using a few synthetic variables (called the discriminant functions). DAPC yields linear combinations of the original variables (alleles) that maximize differences between groups while minimizing variation within clusters. Clusters were defined by independent k-means, using the Bayesian internal criterion (BIC). Twenty DAPC runs were performed in order to assess the stability of clusters. The xvalDAPC function was used to choose the optimal number of PCs for DAPC.

To detect signatures of recombination and population structure within field collections, multilocus linkage disequilibrium (LD) was estimated using the index of association calculated by the R package poppr (88). The poppr tool calculates both Brown's index (IA), which was found to increase with the number of loci, and rd, which corrects for this scaling and forces the index to lie between 0 (linkage equilibrium) and 1 (full disequilibrium). To test the significance of both indices, the package generated 999 resamplings of the data under the null hypothesis (H0) of free recombination, using the “permutation over alleles” and “parametric bootstrap” methods, and then calculated the P value.

ACKNOWLEDGMENTS

We thank S. Lebon, J. M. Baptiste, and E. Lallemand for their technical assistance in both the field and the laboratory and O. Pruvost and C. Vernière for enriching discussions.

This research was financially supported by the European Union (FEDER project), the Conseil Régional de la Réunion, and the French Ministry of Agriculture (project CASDAR C-2013-09 “RESAUBER”).

FOOTNOTES

    • Received 14 November 2016.
    • Accepted 13 December 2016.
    • Accepted manuscript posted online 21 December 2016.
  • Supplemental material for this article may be found at https://doi.org/10.1128/AEM.03095-16 .

  • Copyright © 2017 American Society for Microbiology.

All Rights Reserved .

REFERENCES

  1. 1.↵
    1. Grunwald NJ,
    2. Goss EM
    . 2011. Evolution and population genetics of exotic and re-emerging pathogens: novel tools and approaches. Annu Rev Phytopathol49:249–267. doi:10.1146/annurev-phyto-072910-095246.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Leach JE,
    2. Cruz CMV,
    3. Bai JF,
    4. Leung H
    . 2001. Pathogen fitness penalty as a predictor of durability of disease resistance genes. Annu Rev Phytopathol39:187–224. doi:10.1146/annurev.phyto.39.1.187.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. McDonald BA,
    2. Linde C
    . 2002. Pathogen population genetics, evolutionary potential, and durable resistance. Annu Rev Phytopathol40:349–379. doi:10.1146/annurev.phyto.40.120501.101443.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    1. van Belkum A
    . 1999. The role of short sequence repeats in epidemiologic typing. Curr Opin Microbiol2:306–311. doi:10.1016/S1369-5274(99)80053-8.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Vergnaud G,
    2. Pourcel C
    . 2006. Multiple locus VNTR (variable number of tandem repeat) analysis, p 83–104. InStackebrandt E (ed), Molecular identification, systematics, and population structure of prokaryotes. Springer-Verlag, Berlin, Germany.
  6. 6.↵
    1. Roumagnac P,
    2. Gagnevin L,
    3. Pruvost O,
    4. Achtman M
    . 2007. Insights into structure and evolution of bacterial species that are revealed by molecular methods, p 473–492. InTibayrenc M (ed), Encyclopedia of infectious diseases: modern methodologies. Wiley-Liss, Hoboken, NJ.
  7. 7.↵
    1. Keim P,
    2. Price LB,
    3. Klevytska AM,
    4. Smith KL,
    5. Schupp JM,
    6. Okinaka R,
    7. Jackson PJ,
    8. Hugh-Jones ME
    . 2000. Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis. J Bacteriol182:2928–2936. doi:10.1128/JB.182.10.2928-2936.2000.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Le Fleche P,
    2. Hauck Y,
    3. Onteniente L,
    4. Prieur A,
    5. Denoeud F,
    6. Ramisse V,
    7. Sylvestre P,
    8. Benson G,
    9. Ramisse F,
    10. Vergnaud G
    . 2001. A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol1:2. doi:10.1186/1471-2180-1-2.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Li Y,
    2. Cui Y,
    3. Hauck Y,
    4. Platonov ME,
    5. Dai E,
    6. Song Y,
    7. Guo Z,
    8. Pourcel C,
    9. Dentovskaya SV,
    10. Anisimov AP,
    11. Yang R,
    12. Vergnaud G
    . 2009. Genotyping and phylogenetic analysis of Yersinia pestis by MLVA: insights into the worldwide expansion of Central Asia plague foci. PLoS One4:e6000. doi:10.1371/journal.pone.0006000.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Jiang Y,
    2. Liu HC,
    3. Zheng HJ,
    4. Tang B,
    5. Dou XF,
    6. Zhao XQ,
    7. Zhu YQ,
    8. Lu B,
    9. Wang SY,
    10. Dong HY,
    11. Zhao GP,
    12. Zhang YY,
    13. Kan B,
    14. Wan KL
    . 2012. Evaluation of four candidate VNTR loci for genotyping 225 Chinese clinical Mycobacterium tuberculosis complex strains. Biomed Environ Sci25:82–90. doi:10.3967/0895-3988.2012.01.012.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Le Flèche P,
    2. Fabre M,
    3. Denoeud F,
    4. Koeck JL,
    5. Vergnaud G
    . 2002. High resolution, on-line identification of strains form the Mycobacterium tuberculosis complex based on tandem repeat typing. BMC Microbiol2:37. doi:10.1186/1471-2180-2-37.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Pearson T,
    2. U'Ren J,
    3. Schupp JM,
    4. Allan GJ,
    5. Foster PG,
    6. Mayo MJ,
    7. Gal D,
    8. Choy JL,
    9. Daugherty RL,
    10. Kachur S,
    11. Friedman CL,
    12. Leadem B,
    13. Georgia S,
    14. Hornstra H,
    15. Vogler AJ,
    16. Wagner DM,
    17. Keim P,
    18. Currie BJ
    . 2007. VNTR analysis of selected outbreaks of Burkholderia pseudomallei in Australia. Infect Genet Evol7:416–423. doi:10.1016/j.meegid.2006.12.002.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. U'Ren JM,
    2. Schupp JM,
    3. Pearson T,
    4. Hornstra H,
    5. Friedman CL,
    6. Smith KL,
    7. Daugherty RR,
    8. Rhoton SD,
    9. Leadem B,
    10. Georgia S,
    11. Cardon M,
    12. Huynh LY,
    13. DeShazer D,
    14. Harvey SP,
    15. Robison R,
    16. Gal D,
    17. Mayo MJ,
    18. Wagner DM,
    19. Currie BJ,
    20. Keim P
    . 2007. Tandem repeat regions within the Burkholderia pseudomallei genome and their application for high resolution genotyping. BMC Microbiol7:23. doi:10.1186/1471-2180-7-23.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Pourcel C,
    2. Hormigos K,
    3. Onteniente L,
    4. Sakwinska O,
    5. Deurenberg RH,
    6. Vergnaud G
    . 2009. Improved multiple-locus variable-number tandem-repeat assay for Staphylococcus aureus genotyping, providing a highly informative technique together with strong phylogenetic value. J Clin Microbiol47:3121–3128. doi:10.1128/JCM.00267-09.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Lindstedt BA
    . 2005. Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis26:2567–2582. doi:10.1002/elps.200500096.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Van Belkum A
    . 2007. Tracing isolates of bacterial species by multilocus variable number of tandem repeat analysis (MLVA). FEMS Immunol Med Microbiol49:22–27. doi:10.1111/j.1574-695X.2006.00173.x.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Parkinson N,
    2. Bryant R,
    3. Bew J,
    4. Conyers C,
    5. Stones R,
    6. Alcock M,
    7. Elphinstone J
    . 2013. Application of variable number tandem repeat (VNTR) typing to discriminate Ralstonia solanacearum strains associated with English watercourses and disease outbreaks. Appl Environ Microbiol79:6016–6022. doi:10.1128/AEM.01219-13.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Gironde S,
    2. Manceau C
    . 2012. 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 Microbiol78:3266–3279. doi:10.1128/AEM.06655-11.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Arrieta-Ortiz ML,
    2. Rodriguez RL,
    3. Perez-Quintero A,
    4. Poulin L,
    5. Diaz AC,
    6. Arias Rojas N,
    7. Trujillo C,
    8. Restrepo Benavides M,
    9. Bart R,
    10. Boch J,
    11. Boureau T,
    12. Darrasse A,
    13. David P,
    14. Duge de Bernonville T,
    15. Fontanilla P,
    16. Gagnevin L,
    17. Guerin F,
    18. Jacques MA,
    19. Lauber E,
    20. Lefeuvre P,
    21. Medina C,
    22. Medina E,
    23. Montenegro N,
    24. Munoz Bodnar A,
    25. Noel LD,
    26. Ortiz Quinones JF,
    27. Osorio D,
    28. Pardo C,
    29. Patil PB,
    30. Poussier S,
    31. Pruvost O,
    32. Robene-Soustrade I,
    33. Ryan RP,
    34. Tabima J,
    35. Urrego Morales OG,
    36. Verniere C,
    37. Carrere S,
    38. Verdier V,
    39. Szurek B,
    40. Restrepo S,
    41. Lopez C,
    42. Koebnik R,
    43. Bernal A
    . 2013. Genomic survey of pathogenicity determinants and VNTR markers in the cassava bacterial pathogen Xanthomonas axonopodis pv. manihotis strain CIO151. PLoS One8:e79704. doi:10.1371/journal.pone.0079704.
    OpenUrlCrossRef
  20. 20.↵
    1. N′Guessan CA,
    2. Brisse S,
    3. Le Roux-Nio A-C,
    4. Poussier S,
    5. Koné D,
    6. Wicker E
    . 2013. Development of variable number of tandem repeats typing schemes for Ralstonia solanacearum, the agent of bacterial wilt, banana Moko disease and potato brown rot. J Microbiol Methods92:366–374. doi:10.1016/j.mimet.2013.01.012.
    OpenUrlCrossRef
  21. 21.↵
    1. Bui Thi Ngoc L,
    2. Vernière C,
    3. Vital K,
    4. Guérin F,
    5. Gagnevin L,
    6. Brisse S,
    7. Ah-You N,
    8. Pruvost O
    . 2009. Development of fourteen minisatellite markers for the citrus canker bacterium, Xanthomonas citri pv. citri. Mol Ecol Res9:125–127. doi:10.1111/j.1755-0998.2008.02242.x.
    OpenUrlCrossRef
  22. 22.↵
    1. Bui Thi Ngoc L
    . 2009. Apport de la génétique évolutive dans la taxonomie et l'épidémiologie de Xanthomonas citri pv. citri, bactérie responsable du chancre asiatique des agrumes. PhD thesis. Université de la Réunion, Réunion, France.
  23. 23.↵
    1. Leduc A,
    2. Traore YN,
    3. Boyer K,
    4. Magne M,
    5. Grygiel P,
    6. Juhasz CC,
    7. Boyer C,
    8. Guerin F,
    9. Wonni I,
    10. Ouedraogo L,
    11. Verniere C,
    12. Ravigne V,
    13. Pruvost O
    . 2015. Bridgehead invasion of a monomorphic plant pathogenic bacterium: Xanthomonas citri pv. citri, an emerging citrus pathogen in Mali and Burkina Faso. Environ Microbiol17:4429–4442. doi:10.1111/1462-2920.12876.
    OpenUrlCrossRef
  24. 24.↵
    1. Buhlmann A,
    2. Dreo T,
    3. Rezzonico F,
    4. Pothier JF,
    5. Smits TH,
    6. Ravnikar M,
    7. Frey JE,
    8. Duffy B
    . 2014. Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites. Environ Microbiol16:2112–2125. doi:10.1111/1462-2920.12289.
    OpenUrlCrossRef
  25. 25.↵
    1. Poulin L,
    2. Grygiel P,
    3. Magne M,
    4. Gagnevin L,
    5. Rodriguez RL,
    6. Forero Serna N,
    7. Zhao S,
    8. El Rafii M,
    9. Dao S,
    10. Tekete C,
    11. Wonni I,
    12. Koita O,
    13. Pruvost O,
    14. Verdier V,
    15. Verniere C,
    16. Koebnik R
    . 2015. New multilocus variable-number tandem-repeat analysis tool for surveillance and local epidemiology of bacterial leaf blight and bacterial leaf streak of rice caused by Xanthomonas oryzae. Appl Environ Microbiol81:688–698. doi:10.1128/AEM.02768-14.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Daverdin G,
    2. Rouxel T,
    3. Gout L,
    4. Aubertot JN,
    5. Fudal I,
    6. Meyer M,
    7. Parlange F,
    8. Carpezat J,
    9. Balesdent MH
    . 2012. Genome structure and reproductive behaviour influence the evolutionary potential of a fungal phytopathogen. PLoS Pathog8:e1003020. doi:10.1371/journal.ppat.1003020.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Vernière C,
    2. Bui Thi Ngoc L,
    3. Jarne P,
    4. Ravigné V,
    5. Guérin F,
    6. Gagnevin L,
    7. Le Mai N,
    8. Chau NM,
    9. Pruvost O
    . 2014. Highly polymorphic markers reveal the establishment of an invasive lineage of the citrus bacterial pathogen Xanthomonas citri pv. citri in its area of origin. Environ Microbiol16:2226–2237. doi:10.1111/1462-2920.12369.
    OpenUrlCrossRef
  28. 28.↵
    1. Dunn AR,
    2. Bruening SR,
    3. Grunwald NJ,
    4. Smart CD
    . 2014. Evolution of an experimental population of Phytophthora capsici in the field. Phytopathology104:1107–1117. doi:10.1094/PHYTO-12-13-0346-R.
    OpenUrlCrossRef
  29. 29.↵
    1. Allen C,
    2. Prior P,
    3. Hayward AC
    (ed). 2005. Bacterial wilt disease and the Ralstonia solanacearum species complex. APS Press, Saint Paul, MN.
  30. 30.↵
    1. Gillings MR,
    2. Fahy P
    . 1994. Genomic fingerprinting: towards a unified view of the Pseudomonas solanacearum species complex, p 95–112. InHayward AC, Hartman GL (ed), Bacterial wilt: the disease and its causative agent, Pseudomonas solanacearum. CAB International, Wallingford, United Kingdom.
  31. 31.↵
    1. Elphinstone JG
    . 2005. The current bacterial wilt situation: a global overview, p 9–28. InAllen C, Prior P, Hayward AC (ed), Bacterial wilt disease and the Ralstonia solanacearum species complex. APS Press, Saint Paul, MN.
  32. 32.↵
    1. Hayward AC
    . 1991. Biology and epidemiology of bacterial wilt caused by Pseudomonas solanacearum. Annu Rev Phytopathol29:65–87. doi:10.1146/annurev.py.29.090191.000433.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    1. Vinatzer BA,
    2. Bull CT
    . 2009. The impact of genomic approaches on our understanding of diversity and taxonomy of plant pathogenic bacteria, p 37–62. InJackson RW (ed), Plant pathogenic bacteria—genomics and molecular biology. Caister Academic Press, Reading, United Kingdom.
  34. 34.↵
    1. Wicker E,
    2. Lefeuvre P,
    3. de Cambiaire JC,
    4. Lemaire C,
    5. Poussier S,
    6. Prior P
    . 2012. Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA. ISME J6:961–974. doi:10.1038/ismej.2011.160.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Safni I,
    2. Cleenwerck I,
    3. De Vos P,
    4. Fegan M,
    5. Sly L,
    6. Kappler U
    . 2014. Polyphasic taxonomic revision of the Ralstonia solanacearum species complex: proposal to emend the descriptions of Ralstonia solanacearum and Ralstonia syzygii and reclassify current R. syzygii strains as Ralstonia syzygii subsp. syzygii subsp. nov., R. solanacearum phylotype IV strains as Ralstonia syzygii subsp. indonesiensis subsp. nov., banana blood disease bacterium strains as Ralstonia syzygii subsp. celebesensis subsp. nov. and R. solanacearum phylotype I and III strains as Ralstonia pseudosolanacearum sp. nov. Int J Syst Evol Microbiol64:3087–3103. doi:10.1099/ijs.0.066712-0.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Coupat B,
    2. Chaumeille-Dole F,
    3. Fall S,
    4. Prior P,
    5. Simonet P,
    6. Nesme X,
    7. Bertolla F
    . 2008. Natural transformation in the Ralstonia solanacearum species complex: number and size of DNA that can be transferred. FEMS Microbiol Ecol66:14–24. doi:10.1111/j.1574-6941.2008.00552.x.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Hayward AC
    . 1994. The hosts of Pseudomonas solanacearum, p 9–24. InHayward AC, Hartman GL (ed), Bacterial wilt: the disease and its causative agent, Pseudomonas solanacearum. CAB International, Wallingford, United Kingdom.
  38. 38.↵
    1. Wicker E,
    2. Grassart L,
    3. Mian D,
    4. Coranson Beaudu R,
    5. Dufeal D,
    6. Guilbaud C,
    7. Prior P
    . 2002. Cucumis melo, Cucumis sativus, Cucurbita moschata, and Anthurium spp., new hosts of Ralstonia solanacearum in Martinique (French West Indies). Bacterial Wilt Newsl17:20–21.
    OpenUrl
  39. 39.↵
    1. Lebeau A,
    2. Daunay MC,
    3. Frary A,
    4. Palloix A,
    5. Wang JF,
    6. Dintinger J,
    7. Chiroleu F,
    8. Wicker E,
    9. Prior P
    . 2011. Bacterial wilt resistance in tomato, pepper, and eggplant: genetic resources respond to diverse strains in the Ralstonia solanacearum species complex. Phytopathology101:154–165. doi:10.1094/PHYTO-02-10-0048.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. N′Guessan CA,
    2. Abo K,
    3. Fondio L,
    4. Chiroleu F,
    5. Lebeau A,
    6. Poussier S,
    7. Wicker E,
    8. Kone D
    . 2012. So near and yet so far: the specific case of Ralstonia solanacearum populations from Cote d'Ivoire in Africa. Phytopathology102:733–740. doi:10.1094/PHYTO-11-11-0300.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. N′Guessan CA
    . 2013. Phylogénie, structure génétique et diversité de virulence de Ralstonia solanacearum, agent du flétrissement bactérien, en Côte d'Ivoire. PhD thesis. Université de Cocody-Abidjan, Abidjan, Ivory Coast.
  42. 42.↵
    1. McDonald BA,
    2. Linde C
    . 2002. The population genetics of plant pathogens and breeding strategies for durable resistance. Euphytica124:163–180. doi:10.1023/A:1015678432355.
    OpenUrlCrossRefWeb of Science
  43. 43.↵
    1. Remenant B,
    2. Coupat-Goutaland B,
    3. Guidot A,
    4. Cellier G,
    5. Wicker E,
    6. Allen C,
    7. Fegan M,
    8. Pruvost O,
    9. Elbaz M,
    10. Calteau A,
    11. Salvignol G,
    12. Mornico D,
    13. Mangenot S,
    14. Barbe V,
    15. Medigue C,
    16. Prior P
    . 2010. Genomes of three tomato pathogens within the Ralstonia solanacearum species complex reveal significant evolutionary divergence. BMC Genomics11:379. doi:10.1186/1471-2164-11-379.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Remenant B,
    2. de Cambiaire JC,
    3. Cellier G,
    4. Jacobs JM,
    5. Mangenot S,
    6. Barbe V,
    7. Lajus A,
    8. Vallenet D,
    9. Medigue C,
    10. Fegan M,
    11. Allen C,
    12. Prior P
    . 2011. Ralstonia syzygii, the blood disease bacterium and some Asian R. solanacearum strains form a single genomic species despite divergent lifestyles. PLoS One6:e24356. doi:10.1371/journal.pone.0024356.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Remenant B,
    2. Babujee L,
    3. Lajus A,
    4. Medigue C,
    5. Prior P,
    6. Allen C
    . 2012. Sequencing of K60, type strain of the major plant pathogen Ralstonia solanacearum. J Bacteriol194:2742–2743. doi:10.1128/JB.00249-12.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Salanoubat M,
    2. Genin S,
    3. Artiguenave F,
    4. Gouzy J,
    5. Mangenot S,
    6. Arlat M,
    7. Billault A,
    8. Brottier P,
    9. Camus JC,
    10. Cattolico L,
    11. Chandler M,
    12. Choisne N,
    13. Claudel-Renard C,
    14. Cunnac S,
    15. Demange N,
    16. Gaspin C,
    17. Lavie M,
    18. Moisan A,
    19. Robert C,
    20. Saurin W,
    21. Schiex T,
    22. Siguier P,
    23. Thebault P,
    24. Whalen M,
    25. Wincker P,
    26. Levy M,
    27. Weissenbach J,
    28. Boucher CA
    . 2002. Genome sequence of the plant pathogen Ralstonia solanacearum. Nature415:497–502. doi:10.1038/415497a.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Guidot A,
    2. Elbaz M,
    3. Carrère S,
    4. Siri MI,
    5. Pianzzola MJ,
    6. Prior P,
    7. Boucher C
    . 2009. Specific genes from the potato brown rot strains of Ralstonia solanacearum and their potential use for strain detection. Phytopathology99:1105–1112. doi:10.1094/PHYTO-99-9-1105.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Gabriel DW,
    2. Allen C,
    3. Schell M,
    4. Denny TP,
    5. Greenberg JT,
    6. Duan YP,
    7. Flores CZ,
    8. Huang Q,
    9. Clifford JM,
    10. Presting G,
    11. Gonzalez ET,
    12. Reddy J,
    13. Elphinstone J,
    14. Swanson J,
    15. Yao J,
    16. Mulholland V,
    17. Liu L,
    18. Farmerie W,
    19. Patnaikuni M,
    20. Balogh B,
    21. Norman D,
    22. Alvarez A,
    23. Castillo JA,
    24. Jones J,
    25. Saddler G,
    26. Walunas T,
    27. Zhukov A,
    28. Mikhailova N
    . 2006. Identification of open reading frames unique to a select agent: Ralstonia solanacearum race 3 biovar 2. Mol Plant Microb Interact19:69–79. doi:10.1094/MPMI-19-0069.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Li Z,
    2. Wu S,
    3. Bai X,
    4. Liu Y,
    5. Lu J,
    6. Liu Y,
    7. Xiao B,
    8. Lu X,
    9. Fan L
    . 2011. Genome sequence of the tobacco bacterial wilt pathogen Ralstonia solanacearum. J Bacteriol193:6088–6089. doi:10.1128/JB.06009-11.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Grover A,
    2. Azmi W,
    3. Gadewar AV,
    4. Pattanayak D,
    5. Naik PS,
    6. Shekhawat GS,
    7. Chakrabarti SK
    . 2006. Genotypic diversity in a localized population of Ralstonia solanacearum as revealed by random amplified polymorphic DNA markers. J Appl Microbiol101:798–806. doi:10.1111/j.1365-2672.2006.02974.x.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Ramsubhag A,
    2. Lawrence D,
    3. Cassie D,
    4. Fraser R,
    5. Umaharan P,
    6. Prior P,
    7. Wicker E
    . 2012. Wide genetic diversity of Ralstonia solanacearum strains affecting tomato in Trinidad, West Indies. Plant Pathol61:844–857. doi:10.1111/j.1365-3059.2011.02572.x.
    OpenUrlCrossRef
  52. 52.↵
    1. Ravelomanantsoa S,
    2. Robène I,
    3. Chiroleu F,
    4. Guérin F,
    5. Poussier S,
    6. Pruvost O,
    7. Prior P
    . 2016. A novel multilocus variable number tandem repeat analysis typing scheme for African phylotype III strains of the Ralstonia solanacearum species complex. PeerJ4:e1949. doi:10.7717/peerj.1949.
    OpenUrlCrossRef
  53. 53.↵
    1. Ivey MLL,
    2. Gardener BBM,
    3. Opina N,
    4. Miller SA
    . 2007. Diversity of Ralstonia solanacearum infecting eggplant in the Philippines. Phytopathology97:1467–1475. doi:10.1094/PHYTO-97-11-1467.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Poussier S
    . 2000. Exploration de la diversité génétique de Ralstonia solanacearum, agent du flétrissement bactérien. Détection et dynamique des populations dans les réservoirs d'inoculum. PhD thesis. Université de Rennes I, Rennes, France.
  55. 55.↵
    1. Cruz L,
    2. Eloy M,
    3. Quirino F,
    4. Oliveira H,
    5. Tenreiro R
    . 2012. Molecular epidemiology of Ralstonia solanacearum strains from plants and environmental sources in Portugal. Eur J Plant Pathol133:687–706. doi:10.1007/s10658-012-9947-y.
    OpenUrlCrossRef
  56. 56.↵
    1. Ailloud F
    . 2015. Le pouvoir pathogène chez Ralstonia solanacearum phylotype II: génomique intégrative et paysages transcriptomiques en relation avec l'adaptation à l'hôte. PhD thesis. Université de La Réunion, Saint Pierre, Réunion, France.
  57. 57.↵
    1. Prior P,
    2. Ailloud F,
    3. Dalsing BL,
    4. Remenant B,
    5. Sanchez B,
    6. Allen C
    . 2016. Genomic and proteomic evidence supporting the division of the plant pathogen Ralstonia solanacearum into three species. BMC Genomics17:90. doi:10.1186/s12864-016-2413-z.
    OpenUrlCrossRef
  58. 58.↵
    1. Grey B,
    2. Steck T
    . 2001. The viable but nonculturable state of Ralstonia solanacearum may be involved in long-term survival and plant infection. Appl Environ Microbiol67:3866–3872. doi:10.1128/AEM.67.9.3866-3872.2001.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Kimura M,
    2. Crow JF
    . 1964. The number of alleles that can be maintained in a finite population. Genetics49:725–738.
    OpenUrlFREE Full Text
  60. 60.↵
    1. Di Rienzo A,
    2. Peterson AC,
    3. Garza JC,
    4. Valdes AM,
    5. Slatkin M,
    6. Freimer NB
    . 1994. Mutational processes of simple-sequence repeat loci in human populations. Proc Natl Acad Sci U S A91:3166–3170. doi:10.1073/pnas.91.8.3166.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Crow JF,
    2. Kimura M
    . 1970. An introduction to population genetics theory. Harper and Row, New York, NY.
  62. 62.↵
    1. Ellegren H
    . 2000. Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet16:551–558. doi:10.1016/S0168-9525(00)02139-9.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    1. Schlotterer C
    . 2000. Evolutionary dynamics of microsatellite DNA. Chromosoma109:365–371. doi:10.1007/s004120000089.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Pruvost O,
    2. Magne M,
    3. Boyer K,
    4. Leduc A,
    5. Tourterel C,
    6. Drevet C,
    7. Ravigne V,
    8. Gagnevin L,
    9. Guerin F,
    10. Chiroleu F,
    11. Koebnik R,
    12. Verdier V,
    13. Verniere C
    . 2014. A MLVA genotyping scheme for global surveillance of the citrus pathogen Xanthomonas citri pv. citri suggests a worldwide geographical expansion of a single genetic lineage. PLoS One9:e98129. doi:10.1371/journal.pone.0098129.
    OpenUrlCrossRef
  65. 65.↵
    1. Van Elsas JD,
    2. Smalla K
    . 1997. Methods for sampling soil microbes, p 383–390. InHurst CJ, Knudsen GR, McInerney MJ, Stetzenbach CH, Walter MV (ed), Manual of environmental microbiology. American Society for Microbiology Press, Washington, DC.
  66. 66.↵
    1. Imazaki I,
    2. Nakaho K
    . 2010. Pyruvate-amended modified SMSA medium: improved sensitivity for detection of Ralstonia solanacearum. J Gen Plant Pathol76:52–61. doi:10.1007/s10327-009-0208-7.
    OpenUrlCrossRef
  67. 67.↵
    1. Poussier S,
    2. Cheron JJ,
    3. Couteau A,
    4. Luisetti J
    . 2002. Evaluation of procedures for reliable PCR detection of Ralstonia solanacearum in common natural substrates. J Microbiol Methods51:349–359. doi:10.1016/S0167-7012(02)00111-2.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Kelman A
    . 1954. The relationship of pathogenicity in Pseudomonas solanacearum to colony appearance on a tetrazolium medium. Phytopathology44:693–695.
    OpenUrlWeb of Science
  69. 69.↵
    1. Fegan M,
    2. Prior P
    . 2005. How complex is the “Ralstonia solanacearum species complex,” p 449–462. InAllen C, Prior P, Hayward AC (ed), Bacterial wilt disease and the Ralstonia solanacearum species complex. APS Press, Madison, WI.
  70. 70.↵
    1. Opina N,
    2. Tavner F,
    3. Hollway G,
    4. Wang J-F,
    5. Li T-H,
    6. Maghirang R,
    7. Fegan M,
    8. Hayward AC,
    9. Krishnapillai V,
    10. Hong WF,
    11. Holloway BW,
    12. Timmis J
    . 1997. A novel method for development of species and strain-specific DNA probes and PCR primers for identifying Burkholderia solanacearum (formerly Pseudomonas solanacearum). Asia Pac J Mol Biol Biotechnol5:19–30.
    OpenUrl
  71. 71.↵
    1. Wicker E,
    2. Grassart L,
    3. Coranson-Beaudu R,
    4. Mian D,
    5. Guilbaud C,
    6. Fegan M,
    7. Prior P
    . 2007. Ralstonia solanacearum strains from Martinique (French West Indies) exhibiting a new pathogenic potential. Appl Environ Microbiol73:6790–6801. doi:10.1128/AEM.00841-07.
    OpenUrlAbstract/FREE Full Text
  72. 72.↵
    1. Sanchez Perez A,
    2. Mejia L,
    3. Fegan M,
    4. Allen C
    . 2008. Diversity and distribution of Ralstonia solanacearum strains in Guatemala and rare occurrence of tomato fruit infection. Plant Pathol57:320–331. doi:10.1111/j.1365-3059.2007.01769.x.
    OpenUrlCrossRef
  73. 73.↵
    1. Thompson JD,
    2. Higgins DG,
    3. Gibson TJ
    . 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res22:4673–4680. doi:10.1093/nar/22.22.4673.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Tamura K,
    2. Peterson D,
    3. Peterson N,
    4. Stecher G,
    5. Nei M,
    6. Kumar S
    . 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol28:2731–2739. doi:10.1093/molbev/msr121.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Librado P,
    2. Rozas J
    . 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics25:1451–1452. doi:10.1093/bioinformatics/btp187.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Vallenet D,
    2. Belda E,
    3. Calteau A,
    4. Cruveiller S,
    5. Engelen S,
    6. Lajus A,
    7. Le Fevre F,
    8. Longin C,
    9. Mornico D,
    10. Roche D,
    11. Rouy Z,
    12. Salvignol G,
    13. Scarpelli C,
    14. Thil Smith AA,
    15. Weiman M,
    16. Medigue C
    . 2013. MicroScope: an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. Nucleic Acids Res41(Database Issue):D636–D647. doi:10.1093/nar/gks1194.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Grissa I,
    2. Bouchon P,
    3. Pourcel C,
    4. Vergnaud G
    . 2008. On-line resources for bacterial micro-evolution studies using MLVA or CRISPR typing. Biochimie90:660–668. doi:10.1016/j.biochi.2007.07.014.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Denoeud F,
    2. Vergnaud G
    . 2004. Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource. BMC Bioinformatics5:4. doi:10.1186/1471-2105-5-4.
    OpenUrlCrossRefPubMed
  79. 79.↵
    1. Gelfand Y,
    2. Rodriguez A,
    3. Benson G
    . 2007. TRDB—The Tandem Repeats Database. Nucleic Acids Res35(Database Issue):D80–D87. doi:10.1093/nar/gkl1013.
    OpenUrlCrossRefPubMedWeb of Science
  80. 80.↵
    1. Estoup A,
    2. Angers B
    . 1998. Microsatellites and minisatellites for molecular ecology: theoretical and empirical considerations. Adv Mol Ecol306:55–86.
    OpenUrl
  81. 81.↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    . 1990. Basic local alignment search tool. J Mol Biol215:403–410. doi:10.1016/S0022-2836(05)80360-2.
    OpenUrlCrossRefPubMedWeb of Science
  82. 82.↵
    1. Pourcel C,
    2. Vergnaud G
    . 2011. Strain typing using multiple “variable number of tandem repeat” analysis (MLVA) and genetic element CRISPR, p 179–197. InPersing D, Tenover F, Tang Y, Nolte F, Hayden R (ed), Molecular microbiology: diagnostic principles and practice, 2nd ed. ASM Press, Washington, DC.
  83. 83.↵
    1. Nei M
    . 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics89:583–590.
    OpenUrlAbstract/FREE Full Text
  84. 84.↵
    1. Peakall R,
    2. Smouse PE
    . 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes6:288–295.
    OpenUrlCrossRefWeb of Science
  85. 85.↵
    1. Francisco AP,
    2. Vaz C,
    3. Monteiro PT,
    4. Melo-Cristino J,
    5. Ramirez M,
    6. Carrico JA
    . 2012. PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods. BMC Bioinformatics13:87. doi:10.1186/1471-2105-13-87.
    OpenUrlCrossRefPubMed
  86. 86.↵
    1. Paradis E,
    2. Claude J,
    3. Strimmer K
    . 2004. APE: analyses of phylogenetics and evolution in R language. Bioinformatics20:289–290. doi:10.1093/bioinformatics/btg412.
    OpenUrlCrossRefPubMedWeb of Science
  87. 87.↵
    1. Jombart T,
    2. Devillard S,
    3. Balloux F
    . 2010. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet11:94. doi:10.1186/1471-2156-11-94.
    OpenUrlCrossRefPubMed
  88. 88.↵
    1. Kamvar ZN,
    2. Tabima JF,
    3. Grunwald NJ
    . 2014. Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ2:e281. doi:10.7717/peerj.281.
    OpenUrlCrossRefPubMed
  89. 89.↵
    R Development Core Team. 2013. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org.
PreviousNext
Back to top
Download PDF
Citation Tools
New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations
Jérémy Guinard, Anne Latreille, Fabien Guérin, Stéphane Poussier, Emmanuel Wicker
Applied and Environmental Microbiology Feb 2017, 83 (5) e03095-16; DOI: 10.1128/AEM.03095-16

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations
Jérémy Guinard, Anne Latreille, Fabien Guérin, Stéphane Poussier, Emmanuel Wicker
Applied and Environmental Microbiology Feb 2017, 83 (5) e03095-16; DOI: 10.1128/AEM.03095-16
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Evolution, Molecular
genotype
Minisatellite Repeats
phylogeny
Ralstonia solanacearum
VNTR
microevolution
epidemiosurveillance
molecular epidemiology
population biology
bacterial wilt

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336