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/).

Copper-based antimicrobial compounds are widely used to control plant bacterial pathogens. Pathogens have adapted in response to this selective pressure. Xanthomonas citri pv. citri, a major citrus pathogen causing Asiatic citrus canker, was first reported to carry plasmid-encoded copper resistance in Argentina. This phenotype was conferred by the copLAB gene system. The emergence of resistant strains has since been reported in Réunion and Martinique. Using microsatellite-based genotyping and copLAB PCR, we demonstrated that the genetic structure of the copper-resistant strains from these three regions was made up of two distant clusters and varied for the detection of copLAB amplicons. In order to investigate this pattern more closely, we sequenced six copper-resistant X. citri pv. citri strains from Argentina, Martinique and Réunion, together with reference copper-resistant Xanthomonas and Stenotrophomonas strains using long-read sequencing technology. Genes involved in copper resistance were found to be strain dependent with the novel identification in X. citri pv. citri of copABCD and a cus heavy metal efflux resistance-nodulation-division system. The genes providing the adaptive trait were part of a mobile genetic element similar to Tn3-like transposons and included in a conjugative plasmid. This indicates the system's great versatility. The mining of all available bacterial genomes suggested that, within the bacterial community, the spread of copper resistance associated with mobile elements and their plasmid environments was primarily restricted to the Xanthomonadaceae family.

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.