MultiLocus Variable number of tandem repeat Analysis (MLVA) has been extensively used to examine epidemiological and evolutionary issues on monomorphic human pathogenic bacteria, but not on bacterial plant pathogens of agricultural importance albeit such tools would improve our understanding of their epidemiology, as well as of the history of epidemics on a global scale. Xanthomonas citri pv. citri is a quarantine organism in several countries and a major threat for the citrus industry worldwide. We screened the genomes of Xanthomonas citri pv. citri strain IAPAR 306 and of phylogenetically related xanthomonads for tandem repeats. From these in silico data, an optimized MLVA scheme was developed to assess the global diversity of this monomorphic bacterium. Thirty-one minisatellite loci (MLVA-31) were selected to assess the genetic structure of 129 strains representative of the worldwide pathological and genetic diversity of X. citri pv. citri. Based on Discriminant Analysis of Principal Components (DAPC), four pathotype-specific clusters were defined. DAPC cluster 1 comprised strains that were implicated in the major geographical expansion of X. citri pv. citri during the 20th century. A subset of 12 loci (MLVA-12) resolved 89% of the total diversity and matched the genetic structure revealed by MLVA-31. MLVA-12 is proposed for routine epidemiological identification of X. citri pv. citri, whereas MLVA-31 is proposed for phylogenetic and population genetics studies. MLVA-31 represents an opportunity for international X. citri pv. citri genotyping and data sharing. The MLVA-31 data generated in this study was deposited in the Xanthomonas citri genotyping database (http://www.biopred.net/MLVA/).

Molecular epidemiology studies further our understanding of migrations of phytopathogenic bacteria, the major determining factor in their emergence. Asiatic citrus canker, caused by Xanthomonas citri pv. citri, was recently reported in Mali and Burkina Faso, a region remote from other contaminated areas. To identify the origin and pathways of these emergences, we used two sets of markers, minisatellites and microsatellites, for investigating different evolutionary scales. Minisatellite typing suggested the introduction of two groups of strains in Mali (DAPC 1 and DAPC 2), consistent with microsatellite typing. DAPC 2 was restricted to Bamako district, whereas DAPC 1 strains were found much more invasive. The latter strains formed a major clonal complex based on microsatellite data with the primary and secondary founders detected in commercial citrus nurseries and orchards. This suggests that human activities played a major role in the spread of DAPC 1 strains via the movement of contaminated propagative material, further supported by the frequent lack of differentiation between populations from geographically distant nurseries and orchards. Approximate Bayesian Computation analyses supported the hypothesis that strains from Burkina Faso resulted from a bridgehead invasion from Mali. Multi-locus variable number of tandem repeat analysis and Approximate Bayesian Computation are useful for understanding invasion routes and pathways of monomorphic bacterial pathogens.

dMultilocus variable-number tandem-repeat analysis (MLVA) is efficient for routine typing and for investigating the genetic structures of natural microbial populations. Two distinct pathovars of Xanthomonas oryzae can cause significant crop losses in tropical and temperate rice-growing countries. Bacterial leaf streak is caused by X. oryzae pv. oryzicola, and bacterial leaf blight is caused by X. oryzae pv. oryzae. For the latter, two genetic lineages have been described in the literature. We developed a universal MLVA typing tool both for the identification of the three X. oryzae genetic lineages and for epidemiological analyses. Sixteen candidate variable-number tandem-repeat (VNTR) loci were selected according to their presence and polymorphism in 10 draft or complete genome sequences of the three X. oryzae lineages and by VNTR sequencing of a subset of loci of interest in 20 strains per lineage. The MLVA-16 scheme was then applied to 338 strains of X. oryzae representing different pathovars and geographical locations. Linkage disequilibrium between MLVA loci was calculated by index association on different scales, and the 16 loci showed linear Mantel correlation with MLSA data on 56 X. oryzae strains, suggesting that they provide a good phylogenetic signal. Furthermore, analyses of sets of strains for different lineages indicated the possibility of using the scheme for deeper epidemiological investigation on small spatial scales. M olecular typing of pathogen populations is essential to gain insight into their genetic diversity and population dynamics in order to elaborate efficient strategies for disease control (1, 2). In agricultural systems, pests are ideally controlled by integrated approaches, including eradication or treatment of diseased organisms and planting of resistant varieties. However, the durability of resistance can be challenged if pathogen diversity is significant. Importantly, gene flow between pathogen populations can facilitate the breakdown of resistance in crop plants (3). Hence, efficient and precise molecular-typing tools for identifying strains and differentiating among related bacterial isolates are essential for microevolutionary reconstruction as a population genetics approach for integrated plant protection.Rice, one of the major crops worldwide, is affected by two bacterial diseases that are caused by strains of Xanthomonas oryzae, bacterial leaf blight (BLB), caused by X. oryzae pv. oryzae, and bacterial leaf streak (BLS), caused by X. oryzae pv. oryzicola. Collectively, these two diseases cause significant yield losses in tropical and temperate rice-growing areas. X. oryzae pv. oryzae colonizes xylem vessels upon entry into the vascular system. X. oryzae pv. oryzicola infects the plant via natural openings and colonizes the mesophyll (4). Genomes of members of both pathovars have been sequenced; however, the determinants of tissue specificity are still largely unknown (5, 6). While X. oryzae pv. oryzicola has been shown to be seedborne and seed transmitted (7,8), the evidenc...