# Multilocus variable number of tandem repeat analysis reveals multiple introductions in spain of xanthomonas arboricola pv. pruni, the causal agent of bacterial spot disease of stone fruits and almond

Pablo Lopez-Soriano, Karine Boyer, Denis Cesbron, Maria Clara Morente, Javier Penalver, Ana Palacio-Bielsa, Christian Vernière, Maria M. Lopez, Olivier Pruvost

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Pablo Lopez-Soriano, Karine Boyer, Denis Cesbron, Maria Clara Morente, Javier Penalver, et al. Multilocus variable number of tandem repeat analysis reveals multiple introductions in spain of xanthomonas arboricola pv. pruni, the causal agent of bacterial spot disease of stone fruits and almond. PLoS ONE, 2016, 11 (9), e0163729. doi:10.1371/journal.pone.0163729

HAL Id: hal-01601970

https://hal.science/hal-01601970

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

CrossMark
click for updates

RESEARCH ARTICLE

# Multilocus Variable Number of Tandem Repeat Analysis Reveals Multiple Introductions in Spain of *Xanthomonas arboricola* pv. *pruni*, the Causal Agent of Bacterial Spot Disease of Stone Fruits and Almond

Pablo López-Soriano¹, Karine Boyer², Sophie Cesbron³, María Clara Morente¹, Javier Peñalver¹, Ana Palacio-Bielsa⁴, Christian Vernière²,⁵, María M. López¹*, Olivier Pruvost²

1 Instituto Valenciano de Investigaciones Agrarias, Moncada, Valencia, Spain, 2 UMR Peuplement Végétaux et Bioagresseurs en Milieu Tropical, Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Saint-Pierre, La Réunion, France, 3 INRA, UMR1345 IRHS Institut de Recherche en Horticulture et Semences, Beaucouzé, France, 4 Centro de Investigación y Tecnología Agroalimentaria de Aragón. Instituto Agroalimentario de Aragón, IA2 (CITA-Universidad de Zaragoza), Zaragoza, Spain, 5 UMR Biologie et Génétique des Interactions Plante-Parasite, Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Montpellier, France

* mlopez@ivia.es

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** López-Soriano P, Boyer K, Cesbron S, Morente MC, Peñalver J, Palacio-Bielsa A, et al. (2016) Multilocus Variable Number of Tandem Repeat Analysis Reveals Multiple Introductions in Spain of *Xanthomonas arboricola* pv. *pruni*, the Causal Agent of Bacterial Spot Disease of Stone Fruits and Almond. PLoS ONE 11(9): e0163729. doi:10.1371/journal.pone.0163729

**Editor:** Boris Alexander Vinatzer, Virginia Tech, UNITED STATES

**Received:** May 30, 2016

**Accepted:** September 13, 2016

**Published:** September 26, 2016

**Copyright:** © 2016 López-Soriano et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by Instituto Nacional de Tecnología Agraria y Alimentaria. Project RTA-2011-00140-C03-01 (http://www.inia.es) to PLS MCM JP MML, Euphresco project (266505-FP7-ERANET Euphresco II) to PLS MCM JP APB MML, European Social Fund to PLS KB MCM JP APB CV MML OP, European Regional

## Abstract

*Xanthomonas arboricola* pv. *pruni* is the causal agent of the bacterial spot disease of stone fruits, almond and some ornamental *Prunus* species. In Spain it was first detected in 2002 and since then, several outbreaks have occurred in different regions affecting mainly Japanese plum, peach and almond, both in commercial orchards and nurseries. As the origin of the introduction(s) was unknown, we have assessed the genetic diversity of 239 *X. arboricola* pv. *pruni* strains collected from 11 Spanish provinces from 2002 to 2013 and 25 reference strains from international collections. We have developed an optimized multilocus variable number of tandem repeat analysis (MLVA) scheme targeting 18 microsatellites and five minisatellites. A high discriminatory power was achieved since almost 50% of the Spanish strains were distinguishable, confirming the usefulness of this genotyping technique at small spatio-temporal scales. Spanish strains grouped in 18 genetic clusters (conservatively delineated so that each cluster contained haplotype networks linked by up to quadruple-locus variations). Furthermore, pairwise comparisons among populations from different provinces showed a strong genetic differentiation. Our results suggest multiple introductions of this pathogen in Spain and redistribution through contaminated nursery propagative plant material.
Development Fund to PLS KB SC MCM JP APB CV MML OP, Conseil Régional de La Réunion (http://www.regionreunion.com/fr) to KB CV OP, Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD) (http://www.cirad.fr) to KB CV OP, FPI-INIA grant (http://www.inia.es) Instituto Nacional de Tecnología Agraria y Alimentaria to PLS, AIP Bioressources 2012 taxonomic project to SC. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

Geographical expansions of pathogens into new areas are one of the most important factors involved in plant disease emergence [1]. Such introductions have become more frequent over the last decades due to major increases in plant trade and transport [2], which has favored the exchange of contaminated plant material. The application of efficient molecular typing tools to the study of bacterial populations becomes essential for a better understanding of their genetic diversity and may allow the reconstruction of the invasion routes in order to develop preven- tive control strategies.

*Xantomonas arboricola* pv. *pruni* [3] (synonym, *Xanthomonas campestris* pv. *pruni* Smith) is the causal agent of bacterial spot disease of stone fruits, almond and some ornamental species of *Prunus*. It is considered one of the most important bacterial pathogens affecting plants of the genus *Prunus* due to its high economic impact since most of the commercial cultivars of peach, Japanese plum, apricot and almond are susceptible [4, 5]. Crop losses may exceed 10,000 € per ha in plum orchards in favorable years to the development of the disease [4] and, in Spain, decrease production in susceptible almond cultivars may reach 47% [6]. For these and other reasons, this bacterium is listed as a quarantine organism in the European Union (EU) phytosanitary legislation (EU Council Directive 2000/29/EC) and in the European and Mediterranean Plant Protection Organization (EPPO) list (EPPO A2 list).

The disease was first described in North America in 1903 on Japanese plum [7] and thereafter *X. arboricola* pv. *pruni* has been reported in the main stone fruit producing areas from the five continents [8]. In Europe this pathogen is present in Italy, where it is considered endemic [9] and it has emerged in Belgium, France, Germany, the Netherlands, Switzerland, Spain and some Eastern Europe countries [10] (EPPO-PQR, 2015). In the Americas, it is endemic in many areas of Argentina, Brazil, United States and Uruguay. Symptoms are observed on leaves and fruits, and may cause dieback [8, 11]. The severity of disease outbreaks was found strongly related to weather conditions, cultivar susceptibility and agricultural practices. Bacteria survive during winter in cankers, leaf scars and buds in the absence of symptoms [9, 12, 13]. Shepard and Zehr [14] also showed the epiphytic persistence of this pathogen in different plant organs in population ranging from $10^2$ to $10^6$ CFU per gram of fresh weight in asymptomatic plant material. As there is no effective chemical control available, and to be part of a sustainable management, integrated control strategies are advised. In this context, it is essential to control the sanitary status of plant material in nurseries in order to propagate it free from the disease [4]. Currently, in the EU legislation, a visual inspection in nurseries for symptoms is the only requirement to certify plants free of *X. arboricola* pv. *pruni*, but latent infections and/or epi- phytic populations in asymptomatic propagative material have been found [13, 15, 16].

Despite its economic importance, little is known about the population biology and epidemi- ology of this quarantine bacterium. Boudon et al. [17] used fluorescent amplified fragment length polymorphism (FAFLP) analysis and multilocus sequence analysis (MLSA) to study the diversity of 64 *X. arboricola* pv. *pruni* strains isolated from different countries worldwide. No polymorphism was observed by MLSA based on four gene portions indicating that *X. arbori- cola* pv. *pruni* is a monomorphic bacteria [18]. Subsequent studies based on seven gene por- tions highlighted minor polymorphism but further confirmed its monomorphic nature [19]. Similarly, FAFLP globally yielded little polymorphism, slightly larger in America than in Europe and consistent with the history of the reported outbreaks. Similar results were obtained by Barionovi and Scortichini [20] when screening the variability of 47 *X. arboricola* pv. *pruni* strains collected from Italy, Spain and Australia using an integron gene cassette array and BOX-PCR. Kawaguchi [21] also observed low genetic diversity in *X. arboricola* pv. *pruni* Japa- nese strains using inter-simple sequence repeat PCR and repetitive sequence-based (rep)-PCR.
Thus, all these methods lack resolution for deeper genetic diversity analyses at intrapathovar level, especially at small spatio-temporal (e.g. outbreak investigation) scales.

Multilocus variable number of tandem repeat analysis (MLVA) is a reliable genotyping method that proved efficient for assessing the genetic diversity of monomorphic bacterial pathogens [22] and it is based on the variation detected in the size of short repetitive DNA sequences, tandem repeats (TR), which are highly variable regions in bacterial genomes [23]. Variation in TR number is produced mainly by slipped strand mispairing during DNA replication [24]. Most TRs evolve following the stepwise mutation model where new alleles are created by the addition or deletion of a single repeat unit per mutation event [25], or a generalized two-phase model where multiple-repeat mutation frequencies follow a geometric distribution [26]. The first MLVA scheme targeting a bacterial plant pathogen was performed for *Xylella fastidiosa* [27]. In *Xanthomonas* species, MLVA schemes have been developed first for *X. citri* pv. *citri* [28, 29] and later on for *X. oryzae* [30, 31] and *X. axonopodis* pv. *manihotis* [32]. A comprehensive MLVA scheme targeting 26 TR loci was developed to discriminate 61 *X. arboricola* strains of different pathovars (*celebensis*, *corylina*, *fragariae*, *juglandis*, *poinsetticola*, *populi* and *pruni*) [33]. Essakhi et al. [34] used MLVA for assessing the genetic relatedness between *X. arboricola* pv. *juglandis* (a walnut pathogen) and non-pathogenic *X. arboricola* strains also isolated from walnut. Six TR loci were used to subtype 25 *X. arboricola* pv. *pruni* isolates from cherry laurel (*P. laurocerasus*) in the Netherlands [35].

Here, we study the genetic relatedness among strains causing bacterial spot disease in Spain, where *X. arboricola* pv. *pruni* was first detected in 2002 on Japanese plum, and in subsequent years several outbreaks have occurred in the main stone fruit producing regions [5, 8]. Our objective was to assess the genetic diversity in a collection of 239 Spanish strains of *X. arboricola* pv. *pruni* as well as in 25 reference strains. To achieve this goal, we used a MLVA scheme including two sets of molecular markers combining microsatellite and minisatellite loci. This is the first extensive study targeting the intra-pathovar genetic diversity of a large collection of *X. arboricola* pv. *pruni* strains in a single country as an attempt to provide clues for better understanding introduction and dispersal pathways.

# Materials and Methods

## Bacterial strains and DNA extraction

A total of 264 strains of *X. arboricola* pv. *pruni* were used in this study, 239 isolated from 11 Spanish provinces (Alicante, Badajoz, Huelva, Huesca, Lleida, Mallorca, Navarra, Tarragona, Teruel, Valencia and Zaragoza) as shown in Fig 1 and S1 Table. Spanish strains were isolated from commercial orchards, nurseries and experimental plots since the first detection of this pathogen in Spain in 2002 to 2013. The strains used in this study were obtained from samples received as Laboratory of Reference of the Spanish Ministry of Agriculture (IVIA). The Ministry has been consulted about publishing the results and as official laboratory we do not require specific permission for handling these strains.

Strains were routinely cultured on YPGA medium [36] (yeast extract 5 g l⁻¹ [Difco], bacteriological peptone 5 g l⁻¹ [Difco], glucose 10 g l⁻¹ and agar 20 g l⁻¹ [pH 7.1]) at 25°C for 72 h. DNA from bacterial pure cultures was obtained using a simple isopropanol-based extraction method, as previously described [37]. DNA concentration and purity were determined by using a ND1000 spectrophotometer (Thermo Fisher, Alcobendas, Spain) and the preparations were stored at -20°C. Most of the strains were isolated from Japanese plum (*Prunus salicina*) (n = 103), almond (*P. dulcis*) (n = 77), peach (*P. persica*) (n = 39), a few from apricot (*P. armeniaca*), nectarine (*P. persica* var. *nectarine*), flat peach (*P. persica* var. *platycarpa*), and from the rootstocks Cadaman® (*P. persica* x *P. davidiana*), Barrier (*P. davidiana* x *P. persica*), Garnem
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Chronological detection of outbreaks of *Xanthomonas arboricola* pv. *pruni* in Spain. 1, Badajoz, 2002; 2, Valencia, 2003; 3, Zaragoza, 2004; 4, Alicante, 2006; 5, Huesca, 2008; 6, Lleida, 2008; 7, Teruel, 2009; 8, Navarra, 2009; 9, Mallorca, 2010; 10, Tarragona, 2011; 11, Huelva, 2012. In brackets, number of isolates per province.</figcaption></figure>

doi:10.1371/journal.pone.0163729.g001

(P. dulcis x P. persica) and Monegro (P. dulcis x P. persica) (S1 Table). Spanish strains were also compared to 25 representative strains from Argentina, Australia, Brazil, Canada, France, Italy, New Zealand and United States conserved in international collections (S2 Table). A panel test comprising 16 strains representative of the geographical and genetic diversity of the pathogen was used for preliminary primers screening (S1 and S2 Tables).

## Genomic DNA isolation, sequencing and annotation

Genomic DNA from *X. arboricola* pv. *pruni* pathotype strain CFBP 3894 was isolated and purified using the Qiagen genome DNA isolation kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The Genomic DNA quality and quantity were assessed on an aga-rose gel and using a NanoDrop ND-1000 spectrophotometer (the NanoDrop Technologies, Wilmington, DE). Libraries with an average insert size of 350 bp were sequenced using the Illumina HiSeq 2000 platform (Genoscreen, France). Paired-end reads were assembled in contigs using SOAPdenovo 1.05 [38] and Velvet 1.2.02 [39]. The assemblies had a total length of 5.06 MB generating 77 contigs. The G+C content of the sequence was of 65.42%. Annotation was performed using EuGene-PP using similarities with known protein sequences [40]. Annotation
of the genome sequences revealed 4337 putative protein-coding sequences, 55 tRNA, and three rRNA.

## TR selection

In a first step 26 TR loci previously selected by Cesbron et al. [33] were tested with a test panel of 16 Spanish and reference strains from different geographical origins. TR primers were assessed by PCR as follows: 5 to 10 ng of DNA was used as template in a 15 µl reaction mix containing 0.3 µM of each primer, 1X Terra Buffer (Ozyme, Saint Quentin en Yvelines, France), 0.5X Q-solution and 0.375 U Terra Polymerase Mix (Ozyme). All reactions were performed in a Veriti Thermal Cycler (Applied Biosystems, Villebon sur Yvette, France) using the following PCR conditions: initial denaturation at 98°C for 2 min, followed by 35 cycles consisting of denaturation at 98°C for 10 s, annealing at temperatures ranging from 58 to 68°C for 15 s, and extension at 68°C for 1 min, with a final extension step at 68°C for 30 min. PCR products were visualized under UV light in agarose gels (2%) stained with ethidium bromide.

In a second step, the draft genome sequence of the *X. arboricola* pv. *pruni* pathotype strain CFBP 3894 was screened for minisatellites using the Tandem Repeat Finder (TRF) (http://minisatellites.u-psud.fr/; http://tandem.bu.edu/trf/tr.html) [41, 42]. Parameters selected were total length in a range of 50–1000 bp and tandem repeats length ≥ 10 bp. Other parameters were set as default. Twenty-three single TR loci were retained from the genome of the CFBP 3894. Primer pairs targeting single locus alleles were designed in the TR flanking regions using Primer3, included in Geneious 7 software. Primers were tested with the test panel strains by PCR with the same conditions described above. Amplicons were visualized as mentioned above.

## MLVA scheme

A total of 23 primer pairs targeting single locus alleles were used in a multiplex PCR format. Primers were grouped in multiplex pools according to their annealing temperature. Each primer in the PCR mix was 5'-labeled with one of the following fluorescent dyes: 6-FAM, NED, PET and VIC (Applied Biosystems) (Table 1). Each PCR reaction contained 5 to 10 ng of genomic DNA as template in 15 µl mix containing 0.3 to 1.2 µM of each primer, 1X Terra Buffer (Ozyme), 0.5X Q-solution and 0.375 U Terra polymerase (Ozyme). PCR amplifications were performed in a Veriti Thermal Cycler (Applied Biosystems) using the following conditions: 2 min at 98°C for polymerase activation, followed by 25 cycles at 98°C for 10 s, annealing temperatures ranging from 64 to 68°C for 15 s (Table 1), and 68°C for 1 min (except in one pool where elongation time was 2 min 30 s), and a final extension step at 68°C for 30 min. 1 µl of diluted amplicons was mixed with 0.1 µl of GeneScan-500 LIZ or 0.5 µl of GeneScan-1200 LIZ internal size standard (Applied Biosystems) and 10.9 µl or 10.5 µl of Hi-Di formamide (Sigma-Aldrich) (for GeneScan-500 LIZ and GeneScan-1200 LIZ, respectively). Capillary electrophoresis was performed in an ABI PRISM 3130xl Genetic Analyzer and results were analyzed with GeneMapper 4.0 (Applied Biosystems).

## Data scoring and analysis

Fragment sizes obtained for each TR locus using GeneMapper 4.0 (Applied Biosystems) were transformed to tandem repeat numbers and used as input data. According to Pourcel and Vergnaud [43] when a TR array was truncated, the TR number was rounded up to the higher integer.

Strains sharing the same MLVA profile were grouped in haplotypes. Phylogenetic relationships between strains were studied generating a minimum spanning tree (MST) with
<table><caption>Table 1. TR markers tested on Spanish strains of X. arboricola pv. pruni (n = 239), primers, PCR conditions, number of alleles and Nei's genetic diversity (H<sub>T</sub>).</caption><thead><tr><th>Name</th><th>TR length (bp)</th><th>Primers</th><th>Annealing T (°C)</th><th>Primer concentration (µM)</th><th>PCR pool</th><th>Range of repeat numbers</th><th>Number of alleles (H<sub>T</sub>)<sup>a</sup></th></tr></thead><tbody><tr><td>TR51I</td><td>7</td><td>• 5' FAM-CATGGCAGTCAGGTGGATC 3' • 5' CTGCAACTCCCGATTCCCGA 3'</td><td>68</td><td>0.3</td><td>1</td><td>4-6</td><td>3 (0.154)</td></tr><tr><td>TR37I</td><td>7</td><td>• 5' VIC-ATGGAGGATGCGGTTGCGGCT 3' • 5' CCAACAGAACCCCGCACCA 3'</td><td>68</td><td>0.3</td><td>1</td><td>4-6</td><td>3 (0.384)</td></tr><tr><td>TR37II</td><td>6</td><td>• 5' PET-CGTCATGGACGCCCTGGTCAG 3' • 5' CATTGGCATCGGCACGGCTACT 3'</td><td>68</td><td>0.3</td><td>1</td><td>9-11</td><td>3 (0.017)</td></tr><tr><td>TR05-06</td><td>7</td><td>• 5' NED-GTCGACGGGTTCGCGGAAGGT 3' • 5' GTGCAGCACCAGCCAAAGGCA 3'</td><td>68</td><td>0.3</td><td>1</td><td>15-21</td><td>7 (0.795)</td></tr><tr><td>TR40I</td><td>7</td><td>• 5' FAM-TATCAGGCAGCGCACCACT 3' • 5' TGGAATGTGGAGGCTGTTCG 3'</td><td>68</td><td>0.3</td><td>2</td><td>7-8</td><td>2 (0.168)</td></tr><tr><td>TR36I</td><td>7</td><td>• 5' VIC-GCAGGAGAAGGAAAGCGCCAG 3' • 5' CGATCGCATCTGTGTGGGTTAG 3'</td><td>68</td><td>0.3</td><td>2</td><td>5-6</td><td>2 (0.132)</td></tr><tr><td>TR03I</td><td>7</td><td>• 5' PET-GACATTCGCCGGAGTCAG 3' • 5' GGTTGCTTGGTCGTTGATCG 3'</td><td>68</td><td>0.3</td><td>2</td><td>7-8</td><td>2 (0.008)</td></tr><tr><td>TR68I</td><td>7</td><td>• 5' NED-CTTGCGGTACTGGCTGTTC 3' • 5' AAATCATCGGCGCCTGAAAC 3'</td><td>68</td><td>0.3</td><td>2</td><td>9-15</td><td>7 (0.632)</td></tr><tr><td>TR30II</td><td>7</td><td>• 5' FAM-TTCTGCCGTCTTTTCAGGGCTGG 3' • 5' CATCAGTGCGAGGCCACGAAC 3'</td><td>66</td><td>0.3</td><td>3</td><td>7-11</td><td>4 (0.371)</td></tr><tr><td>TR28II</td><td>9</td><td>• 5' VIC-GCGGCATGTTCCGACTGCACC 3' • 5' GGGTGGATGAGGGTCTGCATG 3'</td><td>66</td><td>0.6</td><td>3</td><td>3-4</td><td>2 (0.008)</td></tr><tr><td>TR38II</td><td>6</td><td>• 5' PET-TCTCGGTATCGATGTGGGTGC 3' • 5' CCCGTAGCTGTATCAGTGCCT 3'</td><td>66</td><td>0.9</td><td>3</td><td>6-12</td><td>7 (0.691)</td></tr><tr><td>TR58II</td><td>7</td><td>• 5' NED-GGAAGAGTACCCGGCAATTCT 3' • 5' TCTGATCGGTGCTGAGCGTCT 3'</td><td>66</td><td>0.6</td><td>3</td><td>9-17</td><td>8 (0.667)</td></tr><tr><td>TR33I</td><td>7</td><td>• 5' FAM-CGAGTGGATGTTATGGCGTGG 3' • 5' CTCGCAAAACCCTTGCCATC 3'</td><td>66</td><td>0.3</td><td>4</td><td>5-8</td><td>4 (0.207)</td></tr><tr><td>TR58I</td><td>7</td><td>• 5' VIC-ATCTGTTGCTGGCCGAGAGC 3' • 5' ACCAACACCGAGCTTGCCTC 3'</td><td>66</td><td>0.3</td><td>4</td><td>5-9</td><td>5 (0.355)</td></tr><tr><td>TR79I</td><td>7</td><td>• 5' NED-GCTGATCCTTCGTGGGCTTG 3' • 5' GGTGTGAATTCGTCGGTGAC 3'</td><td>66</td><td>0.3</td><td>4</td><td>6-10</td><td>5 (0.393)</td></tr><tr><td>TR50I</td><td>7</td><td>• 5' VIC-GTTGCGAGATCGGGCGCTTC 3' • 5' CGTGCATCAGACGCTTGCGT 3'</td><td>68</td><td>0.3</td><td>5</td><td>5-9</td><td>5 (0.727)</td></tr><tr><td>TR66I</td><td>8</td><td>• 5' PET-TGCAGTTGTGGTCTTCGGCA 3' • 5' CACGCATCAAGTTCGACATGGTGC 3'</td><td>68</td><td>0.3</td><td>5</td><td>9-12</td><td>4 (0.719)</td></tr><tr><td>TR15I</td><td>7</td><td>• 5' NED-GCCATGTCGCCGGGAAACGA 3' • 5' TCGAGCGGTTCCTGCGGTTGT 3'</td><td>68</td><td>0.3</td><td>5</td><td>6-7</td><td>2 (0.285)</td></tr><tr><td>TR67II</td><td>7</td><td>• 5' VIC-AGATACAAGGCGAACGCGAT 3' • 5' CAGGACAGGAACGGCAACC 3'</td><td>64</td><td>0.3</td><td>6</td><td>10-20</td><td>11 (0.779)</td></tr><tr><td>Xap4790</td><td>112</td><td>• 5' CGCGTATTGCAGGAATCCAC 3' • 5' FAM-CGATTGGAGATCCGGACCAG3'</td><td>68</td><td>0.3</td><td>A</td><td>2-4</td><td>2 (0.041)</td></tr><tr><td>Xap4422<sup>b</sup></td><td>116</td><td>• 5' GGTCTGATCCGCTTCTCACC • 5' -NED-ATCCGCGCCAACTACAAGAA</td><td>68</td><td>0.3</td><td>A</td><td></td><td></td></tr><tr><td>Xap0897</td><td>15</td><td>• 5' ATTACTTCTTCCCTGCGGC 3' • 5' FAM-TCTACAACACCAAGACCGGC 3'</td><td>64</td><td>1.2</td><td>B</td><td>1-3</td><td>2 (0.008)</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>Name</th><th>TR length (bp)</th><th>Primers</th><th>Annealing T (°C)</th><th>Primer concentration (µM)</th><th>PCR pool</th><th>Range of repeat numbers</th><th>Number of alleles (H<sub>T</sub>)<sup>a</sup></th></tr></thead><tbody><tr><td>Xap2280</td><td>20</td><td>• 5' AGAGCCTACACGGACGTACT 3'<br/>• 5' VIC-TGAACGGGATGGTGCAAGTT 3'</td><td>64</td><td>0.3</td><td>B</td><td>7–13</td><td>6 (0.508)</td></tr><tr><td>TR10II</td><td>12</td><td>• 5' TGGTTGCGCCCTTGCCTTCTC 3'<br/>• 5' PET-TCGCCGGCATCAACATGGCCG 3'</td><td>64</td><td>0.9</td><td>B</td><td>5–7</td><td>3 (0.292)</td></tr></tbody></table>

doi:10.1371/journal.pone.0163729.t001

PHYLOViZ v1.0 [44]. Clonal complexes (CC) were defined as networks of single locus variants (SLVs, i.e. groups of strains differing at a single TR locus), and were obtained and analyzed with eBURST v3 [45]. The founder was defined as the haplotype with the largest number of SLVs.

The population structure of the collection of the Spanish strains was assessed using PHILOViZ v1.0. A categorical MST was generated using the algorithm recommended for MLVA data combining global optimal eBURST (goeBURST) and Euclidean distances. Genetic clusters (GCs) were delineated so that each cluster contains a network of haplotypes linked by up to quadruple-locus variations.

## Genetic diversity of Spanish strains

Only Spanish strains were used to calculate genetic diversity indices and to perform analysis of population structure at the province level, although only those provinces where more than 12 strains were isolated were selected (i.e. 224 strains from seven provinces) which represent 94% of total strain collection studied. Nei's unbiased estimates of genetic diversity ($H_E$) were calculated using ARLEQUIN version 3.01 [46]. Allelic richness (A), the mean number of alleles per locus per population, was calculated using the rarefaction procedure with HP-RARE version 1.0 [47]. Genetic differentiation among provinces was evaluated with Wright's fixation index ($F_{ST}$) and Slatkin's $F_{ST}$ analogue ($R_{ST}$), which were calculated by computing distance matrices for haplotypic data. Significance was tested with 1000 permutations using ARLEQUIN version 3.1. [46].

Metric multidimensional scaling (MDS), based on a Manhattan distance matrix, was used to represent distances among the Spanish and the worldwide strains. MDS transforms a distance matrix (which cannot be analyzed by eigen-decomposition) into a cross-product matrix, then solving the eigen-vector problem to find the coordinates of individuals so that distortions to the distance matrix are minimized. As in principal component analysis, individuals are projected into *n* dimensions [48]. MDS was performed using the cmd-scale function in the R software.

## Results

### Selection of TR markers

The panel test of 16 *X. arboricola* pv. *pruni* strains was first tested by the 26 TR loci previously selected by Cesbron et al. [33]. Primers were redesigned when original primers were found to hybridize within the TR array (e.g. TR18II and TR67II) or when mispriming was observed (TR21II). Four TR loci were not further considered (TR19I, TR54I, TR39II and TR18II). In the

<footnote>a Nei's genetic diversity calculated using ARLEQUIN version 3.01.</footnote>

<footnote>b This marker was monomorphic within Spanish strains and was only used to compare them to the worldwide collection.</footnote>
context of xanthomonads, which often contain compound microsatellites in their genome (i.e. loci for which TRs are composed of two or more sequence types), we considered TRs to be different when they had less than 70% sequence identity. TR19I was found to be a combination of two microsatellites sharing only 58% identity. TR54I showed a combination of two microsatellites (showing only 29% identity) and one minisatellite. TR39II and TR18II showed very different profiles in closely related strains (i.e. single- or double-locus variants), in some cases differences were up to 12–14 repeats in the same locus between strains isolated in the same orchard and the same year. Sanger sequencing confirmed that the observed size polymorphism concerned the TRs and not the flanking regions. These differences may be attributed to recombination events and this suggested that these two loci might not follow a stepwise or two-phase mutation model. Therefore, they were not further considered in the analysis.

Moreover, from the genome sequence of strain CFBP 3894, markers TR05II and TR06II corresponded to a same and single tandem repeat locus. Furthermore, TR05II forward primer and TR06II reverse primer were found to prime within the repeated sequence. Therefore, we used the forward primer of TR06II and the reverse primer of TR05II (new marker named TR05-06), which did not prime in the tandem repeat array (S1 Fig).

The panel test of strains was also used to test primer pairs targeting 23 minisatellites identified with TRF screening the draft genome sequence of strain CFBP 3894. Four of them were polymorphic (Xap4790, Xap0897, Xap4422 and Xap2280) and were included in the MLVA scheme, whereas the other markers were monomorphic. Xap4422 was monomorphic among Spanish strains and was only retained to compare Spanish strains to the worldwide collection of strains.

Consequently we eventually retained 23 TR loci (18 microsatellites (TR < 9bp) and 5 minisatellites (TR ≥ 9 bp) [43] for assessing the diversity of the Spanish strains collection and 24 loci for the worldwide collection.

## MLVA genotyping of the Spanish strains

The set of 239 strains of *X. arboricola* pv. *pruni* isolated in Spain was successfully genotyped, and all of them yielded reproducible amplicons for all 23 TR loci. Nei's genetic diversity at each locus ranged from 0.008 to 0.795 (Table 1). Three TR loci (TR03I, TR28II and Xap0897) showed very low levels of polymorphism and differentiated a single isolate in the whole collection. TR37II and Xap4790 were also weakly polymorphic and only differentiated two and four isolates, respectively. Allele number per locus ranged from 2 to 11 (Table 1).

Out of 239 Spanish strains, 119 MLVA haplotypes were identified. The genetic diversity of the whole collection was $H_T = 0.363$. Genetic diversity indices for each province are shown in Table 2. In order to assess the population structure of the Spanish strains, a categorical minimum spanning tree was built. Eighteen GCs were obtained grouping haplotypes differing by up to four TR loci (Fig 2).

The largest GC (i.e. GC1) contained 48 haplotypes and 111 strains, which represent 40 and 46% of the whole collection, respectively. Strains from nine provinces were found in GC1, including all the strains from Valencia, Mallorca and Tarragona. Only those originating from Alicante and Huelva were not linked to this group. Other GCs were much smaller ($\le 16$ haplotypes), and six of them were formed by a single haplotype (GC13, GC14, GC15, GC16, GC17 and GC18). Strains from the province of Badajoz, where the disease was first reported in Spain, and Zaragoza split in six and five different GCs, respectively (Fig 2). All but three GCs (GC1, GC2 and GC7) were composed of a small number of strains originating from a single province with five of them containing strains originating from Badajoz. For example, all strains isolated in Alicante were included in GC4. This was correlated with the allelic richness (A) calculated
<table><caption>Table 2. Genetic diversity estimated from MLVA data of *X. arboricola* pv. *pruni* for the strains from 11 sampled Spanish provinces.</caption><thead><tr><th>Province<sup>a</sup></th><th>N<sup>b</sup></th><th>Na<sup>c</sup></th><th>A<sup>d</sup></th><th>N<sub>H</sub><sup>e</sup></th><th>H<sub>E</sub><sup>f</sup></th></tr></thead><tbody><tr><td>Alicante</td><td>22</td><td>33</td><td>1.30</td><td>8</td><td>0.059</td></tr><tr><td>Badajoz</td><td>57</td><td>62</td><td>2.37</td><td>24</td><td>0.344</td></tr><tr><td>Huesca</td><td>38</td><td>70</td><td>2.56</td><td>30</td><td>0.311</td></tr><tr><td>Mallorca</td><td>12</td><td>26</td><td>1.04</td><td>2</td><td>0.290</td></tr><tr><td>Teruel</td><td>21</td><td>52</td><td>2.09</td><td>14</td><td>0.230</td></tr><tr><td>Valencia</td><td>32</td><td>34</td><td>1.28</td><td>9</td><td>0.053</td></tr><tr><td>Zaragoza</td><td>42</td><td>68</td><td>2.48</td><td>27</td><td>0.281</td></tr><tr><td>Huelva</td><td>6</td><td></td><td></td><td>3</td><td></td></tr><tr><td>Lleida</td><td>3</td><td></td><td></td><td>3</td><td></td></tr><tr><td>Navarra</td><td>2</td><td></td><td></td><td>2</td><td></td></tr><tr><td>Tarragona</td><td>1</td><td></td><td></td><td>1</td><td></td></tr><tr><td>Unknown</td><td>3</td><td></td><td></td><td>1</td><td></td></tr><tr><td>Total</td><td>239</td><td></td><td>4.30</td><td>119</td><td>0.363<sup>g</sup></td></tr></tbody></table>

doi:10.1371/journal.pone.0163729.t002

for each province (Table 2), where Alicante, Mallorca and Valencia showed the lowest values (A ≤ 1.3), while Badajoz, Huesca and Zaragoza displayed the highest values (A > 2.3). Also consistent with the delineation of GCs, pairwise comparisons among populations sampled in seven provinces (n ≥ 12) showed highly significant F<sub>ST</sub> and R<sub>ST</sub> values (Table 3) supporting a strong genetic differentiation between strains originating from different provinces. Only strains from Mallorca and Zaragoza showed a R<sub>ST</sub> value slightly significant (P < 0.05).

Haplotypes grouped in 16 CC and 34 singletons (i.e. strains with no SLV). The correspondence between GCs and CCs is highlighted in Fig 2. Single-repeat variants represented 70.7% of the total number of SLVs, suggesting a predominant stepwise mutation model. Double-repeats variants were observed (25.6% of the total SLVs) only with TR05-06, TR33I, TR38II, TR50I and TR58II loci, which are among the most polymorphic ones. A plausible explanation may be that some intermediate allelic states were missing in our strain collection, or perhaps these loci do not strictly follow a strict stepwise mutation model but a two-phase model [26]. Only three triple-repeat variants were found in TR05-06 locus. Variations involving more than three repeats, which may be attributed to recombination events, were not observed in our dataset.

The largest clonal complex (GC1-A) was formed by 23 haplotypes and 70 strains isolated from six different provinces (Huesca, Lleida, Mallorca, Teruel, Valencia and Zaragoza) between 2003 and 2013 (Fig 2). Two other clonal complexes, GC2-A and GC2-B, were constituted by strains from different provinces (Huesca and Zaragoza). There were provinces with all their strains included in a single clonal complex, e.g. all the strains from Mallorca and Valencia that grouped in GC1-A. Most clonal complexes (13 out of 16) grouped only strains from the same province or the same location.

The primary founder in GC1-A (haplotype 29), the haplotype with the highest number of SLVs (n = 10) corresponded to the most frequent haplotype (n = 18), with strains isolated from

<footnote><sup>a</sup> Only provinces with n ≥ 12 were selected to calculate diversity indices (Na, A and H<sub>E</sub>).</footnote>

<footnote><sup>b</sup> N, number of isolates per province.</footnote>

<footnote><sup>c</sup> Na, number of alleles.</footnote>

<footnote><sup>d</sup> A, allelic richness calculated by rarefaction method (n = 12).</footnote>

<footnote><sup>e</sup> N<sub>H</sub>, number of haplotypes in each province.</footnote>

<footnote><sup>f</sup> H<sub>E</sub>, Nei's genetic diversity within provinces.</footnote>

<footnote><sup>g</sup> H<sub>T</sub>, total Nei's genetic diversity.</footnote>
<figure><img src="image_6.png" /></figure>

**Fig 2.** Categorical minimum spanning tree from MLVA data (239 strains and 119 haplotypes) representing the genetic diversity of the Spanish strains of X. arboricola pv. pruni in relation with its province of origin. Dot diameter represents the number of strains per haplotype. Numbers represent the 18 genetic clusters identified. Haplotypes in the same genetic cluster are up to quadruple locus variants. Thick links indicate single locus variants and thin links indicates double locus variants. Shaded areas show different clonal complexes in genetic clusters 1 and 2 and are identified with letters (A-E). Genetic clusters not enclosed in dashed lines are formed by a unique clonal complex or singleton. *Indicates the primary founder haplotype.

doi:10.1371/journal.pone.0163729.g002

<table><caption>**Table 3.** Genetic differentiation of Spanish strains of X. arboricola pv. pruni from seven provinces (n ≥ 12) estimated by FST (above the diagonal) and RST (below the diagonal) pairwise comparisons based on MLVA data.</caption><thead><tr><th></th><th></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th></tr></thead><tbody><tr><td>1</td><td>Alicante</td><td></td><td>0.535***</td><td>0.634***</td><td>0.907***</td><td>0.703***</td><td>0.885***</td><td>0.643***</td></tr><tr><td>2</td><td>Badajoz</td><td>0.367***<sup>a</sup></td><td></td><td>0.206***</td><td>0.379***</td><td>0.263***</td><td>0.421***</td><td>0.238***</td></tr><tr><td>3</td><td>Huesca</td><td>0.648***</td><td>0.231***</td><td></td><td>0.341***</td><td>0.266***</td><td>0.345***</td><td>0.106***</td></tr><tr><td>4</td><td>Mallorca</td><td>0.961***</td><td>0.341***</td><td>0.307***</td><td></td><td>0.447***</td><td>0.573***</td><td>0.241***</td></tr><tr><td>5</td><td>Teruel</td><td>0.821***</td><td>0.325***</td><td>0.331***</td><td>0.458***</td><td></td><td>0.507***</td><td>0.220***</td></tr><tr><td>6</td><td>Valencia</td><td>0.954***</td><td>0.401***</td><td>0.448***</td><td>0.559***</td><td>0.414***</td><td></td><td>0.176***</td></tr><tr><td>7</td><td>Zaragoza</td><td>0.679***</td><td>0.268***</td><td>0.134***</td><td>0.091*</td><td>0.228***</td><td>0.175***</td><td></td></tr></tbody></table>

<sup>a</sup> Significance level of FST and RST pairwise comparisons

*P* < 0.05;

**P* < 0.01;

***P* < 0.001;

NSP* > 0.005 (non-significant).

doi:10.1371/journal.pone.0163729.t003
4 different provinces (Huesca n = 1, Lleida n = 1, Valencia n = 12 and Zaragoza n = 4) (Fig 2). In addition, two other haplotypes sharing strains isolated from different provinces were identified in this CC: haplotype 30 composed by eight strains isolated from Mallorca (n = 2) and Valencia (n = 6) and haplotype 32 formed by three strains from Huesca (n = 1) and Zaragoza (n = 2).

Sixty-three strains isolated in seven Spanish nurseries located in five different provinces were distributed in five different GCs (S1 Table). Among them, 32 strains isolated from nursery plants in Valencia in 2003 and seven strains isolated from nursery plants in Zaragoza between 2004 and 2008 grouped into the major clonal complex (GC1-A). Among those, 13 strains (Valencia n = 12, Zaragoza n = 1) were assigned to the primary founder haplotype, suggesting a probable transmission of contaminated plant material by human activities. Fourteen strains were isolated in a single nursery located in Badajoz in 2002 showing very different MLVA profiles. Indeed, two of them clustered in GC1 and the other ones in GC6, suggesting a possible polyclonal or multiple introduction(s) through the same nursery. Other strains from nursery plants were clearly separated with different allelic profiles, such as strains from Zaragoza clustering in GC8 or strains from Huesca and Lleida nurseries in GC13 and GC16, respectively.

Six haplotypes grouped strains isolated both in nurseries and orchards. For three of them the corresponding nurseries and the orchards were located in the same province (two in Badajoz and the other one in Zaragoza). For the three other cases, the nurseries were located in different provinces from those of the orchards.

## MLVA genotyping on a worldwide collection

The 25 *X. arboricola* pv. *pruni* strains from our world collection were classified in 23 haplotypes, which means that almost every strain formed a unique haplotype. No haplotype was shared by strains from different countries. MLVA profiles clearly separated the 239 Spanish strains from the world strain collection, showing in most of the cases differences in a range of 4 to 15 loci, even within strains originated from a same country (haplotypes from Australia resulted from 5 to 9-loci variants) (Fig 3). Only two strains, CFBP 3894 from New Zealand and CFBP 5530 from Italy, were SLV (TR79I) and double locus variants (DLV) (TR33I and TR67II) of a Spanish haplotype, respectively, as shown by their close genetic distance on the MDS (Fig 3).

## Discussion

The understanding of the invasion routes, biology and epidemiology of plant pathogens is of striking importance in order to elucidate which are the main factors involved in the success of invasions and to develop strategies aiming at preventing new introductions. Traditionally, two types of methods are employed to retrace invasion routes: direct methods, based on recorded data by quarantine services or agricultural agencies and indirect methods, based on genetic data obtained by the use of molecular markers [49]. In this study we have tried to combine both of them, using all the available information from the 239 Spanish strains analyzed.

## MLVA is sufficiently resolutive for molecular epidemiology analyses of *X. arboricola* pv. *pruni* over small spatial or temporal scales

The analysis of literature data suggested that *X. arboricola* pv. *pruni* is a monomorphic bacterium and all available genotyping techniques but MLVA are not suitable for deciphering the genetic structure of strains during outbreaks. MLVA had been used to distinguish among *X. arboricola* pathovars [33] and for assessing the genetic diversity among Dutch strains of *X. arboricola* pv. *pruni* isolated from of *P. lauroceraus* [35]. Using a comprehensive strain
<figure><img src="image_7.png" /><figcaption>**Fig 3.** MDS representation of the distances among 264 strains of *X. arboricola* pv. *pruni*. Countries of origin are represented by different symbols. MDS axes 1–2, 1–3 and 1–4 described 59.3%, 47.4% and 44.4% of the total variation, respectively.</figcaption></figure>

doi:10.1371/journal.pone.0163729.g003

collection, we implemented and evaluated a MLVA scheme derived and improved from Cesbron et al. [33] in order to assess the genetic structure of *X. arboricola* pv. *pruni*, eleven years after its emergence in Spain [8].

Around 50% of the Spanish strains were discriminated as different haplotypes, which confirms the appropriateness of MLVA as a genotyping tool for small-scale analyses of this monomorphic pathogen in a relatively small geographical area compared to other techniques such as MLST [17]. A set of 11 Spanish strains shown genetically diverse based on MLVA (i.e. assigned
to different GCs) constituted a single MLST haplotype using four housekeeping genes (Garita-Cambronero et al., personal communication).

Conversely and consistent with earlier data on *X. citri* pv. *citri* [29, 50] and *X. oryzae* [31], our MLVA scheme, which is mainly based on microsatellites, was not appropriate for understanding deep phylogenies and therefore strain relatedness at a global scale.

Our study did not aim at analyzing the role of the hosts in structuring the genetic diversity of *X. arboricola* pv. *pruni* because the strain collection available could not allow to extensively address this question. However, based on partial data that have been analyzed for this purpose, no strong or clear structure by the host was observed (data not shown) and several hosts from different provinces shared populations that were not statistically differentiated.

## MLVA supports multiple introductions of *X. arboricola* pv. *pruni* in Spain

Since the original description of bacterial spot disease of stone fruits and almond in Badajoz in 2002, the outbreaks caused by *X. arboricola* pv. *pruni* in Spain are likely the result of multiple and/or polyclonal introductions. A scenario involving multiple introductions of a pest is a common phenomenon in invasions [51], and it has already been described in some plant pathogenic bacteria such as *Erwinia amylovora* [52, 53].

This hypothesis was suggested by several independent pieces of evidence. Firstly, the 239 Spanish strains were assigned to 18 conservatively delineated GCs (i.e. strains in a network linked by up to quadruple-locus variations were assigned to a single GC). Secondly, pairwise comparisons among populations from Spanish provinces showed highly significant values of FST and RST among seven provinces, indicating a strong genetic differentiation. Thirdly, the assignation of isolates collected in 2002 in the Badajoz province further suggested that the initial establishment of the pathogen in Spain was caused by at least four distinct GCs (GC1, GC6, GC9 and GC12).

Few studies focus on the genetic diversity of local emerging strains of a *Xanthomonas* species within the first years following its detection. The situation observed in Spain for *X. arboricola* pv. *pruni* was strikingly different from the one reported for *X. citri* pv. *citri* strains having emerged in Mali and Burkina Faso [50]. In the latter study, strains from two genetic lineages were identified but only one was detected in commercial citrus nurseries making it markedly prevalent likely through the massive spread of diseased propagative plant material. As a result, a large majority of emerging strains clustered as a single large CC and all epidemic strains would have been assigned to a single GC using the same criteria for GC delineation as herein. On the contrary, the genetic structure revealed in the present study was markedly different, suggestive of more or less concomitant outbreaks having a distinct origin.

Some GCs (i.e. GC3, GC4 and GC10) were formed by a few strains that were isolated in the same orchard and in the same year. These groups of strains showed MLVA profiles completely different to the rest of Spanish strains suggesting that these outbreaks were spatially restricted. These groups are probably the result of local dispersal as a result of natural dispersion by rain and wind or mechanical transmission during agricultural practices within the orchard. Interestingly, strains from these GCs have not been detected from nurseries.

## MLVA together with records from plant protection agencies support the role of nurseries in the introduction and spread of *X. arboricola* pv. *pruni* in Spain

Nursery strains analyzed herein clustered in five GCs. Several strains from nurseries in different provinces grouped in the major CC (found in cluster GC1). Interestingly, the primary
founder in this CC shared strains from different provinces that were isolated both from two nurseries and several orchards. Field investigations indicated that there had been an exchange of plant material between these two nurseries located in Valencia and Zaragoza, from which strains isolated in 2003 and 2004 respectively, were assigned to this primary founder (haplotype 29). Consistent with these results, several interceptions of contaminated plant material have been reported in Spanish nurseries and most of them consisted of symptomless plants with latent infections [5]. X. arboricola pv. pruni has been detected in seedlings of nectarine, Japanese plum and different stone fruit rootstocks, some of them originating from Italy where the disease is endemic in some regions [9]. A striking example of an interception was documented in GC8, which is composed of four strains from the rootstocks Garnem and Cadaman®. All of them were isolated in a nursery from Zaragoza in 2010 but the plants were produced in an Italian nursery (unpublished data). Such strains presented an allelic profile very different from the rest of the Spanish strains.

The largest CC identified herein was formed by strains originated from six different provinces. Some of them, as Huesca, Lleida, Teruel and Zaragoza belong to one of the greatest stone fruit producing areas in Spain, and some of the orchards where the strains were collected are separated by distances ranging from 5 to 30 km. Although natural spread of X. arboricola pv. pruni may occur within orchards during rain showers and hail storms [9, 54, 55], its spread has not been documented over such large distances. To our knowledge, the sole documented case involving extreme weather events (not reported in Spain) was described for X. citri pv. citri in a subtropical environment by Irey et al. [56]. Contaminated nursery plant material may have been associated with such long distance dissemination.

X. arboricola pv. pruni is a pathogen which can survive epiphytically or endophytically in different organs of the plants, such as leaves, twigs, buds, flowers and fruits at low concentration, without causing symptoms of bacterial spot disease [14, 57, 58]. These populations are a potential inoculum source for the development of the disease when weather conditions are favorable and may be transported between nurseries or from them to orchards inadvertently. Our results suggest that nursery plant material has contributed to spread X. arboricola pv. pruni in Spain and that the current inspection procedure consisting of visual observations of typical symptoms in plants in nurseries is insufficient to ascertain their sanitary status. Consequently, survey and sampling protocols combined with high accuracy detection techniques should be applied [59, 60]. There is a need to produce new research data combining microbial ecology and genotyping to improve the knowledge on the biological significance of X. arboricola pv. pruni epiphytic and endophytic populations as a primary inoculum source. Assessing mutation rates of TR markers in future experimental evolution studies would also be valuable.

**Supporting Information**

S1 Fig. Nucleotide sequence in *Xanthomonas arboricola* pv. *pruni* strain CFBP 3894 of TR loci that were excluded from the analysis or for which new primer combinations were designed.
(PDF)

S1 Table. *Xanthomonas arboricola* pv. *pruni* strains isolated from 11 provinces in Spain.
(PDF)

S2 Table. X. arboricola pv. pruni strains from international collections.
(PDF)
Acknowledgments

We especially thank Spanish Diagnostic Laboratories of Badajoz, Barcelona, Huelva, Mallorca, Valencia and Zaragoza for sending samples and information, José Miguel Quesada and Isabel M. Berruete for their excellent technical work and Jerson Garita-Cambronero and Jaime Cubero for sharing unpublished data. We thank Géraldine Taghouti for genomic DNA extraction and Martial Briand, Jerome Gouzy and Sébastien Carrère for performing assembly and automatic annotation of the genome.

Author Contributions

Conceptualization: PLS CV MML OP.

Data curation: PLS KB CV MML OP.

Formal analysis: CV OP.

Funding acquisition: CV MML OP.

Investigation: PLS KB SC MCM JP.

Methodology: PLS KB CV MML OP.

Project administration: CV MML OP.

Resources: SC APB CV MML OP.

Supervision: CV MML OP.

Validation: PLS CV MML OP.

Visualization: PLS MML.

Writing – original draft: PLS.

Writing – review & editing: APB CV MML OP.

References

1. Anderson PK, Cunningham AA, Patel NG, Morales FJ, Epstein PR, Daszak P. Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers. *Trends Ecol Evol.* 2004; 19: 535–544. PMID: 16701319

2. Hulme PE. Trade, transport and trouble: managing invasive species pathways in an era of globalization. *J Appl Ecol.* 2009; 46: 10–18.

3. Vauterin L, Hoste B, Kersters K, Swings J. Reclassification of *Xanthomonas*. *Int J Syst Bacteriol.* 1995; 45: 472–489.

4. Stefani E. Economic significance and control of bacterial spot/canker of stone fruits caused by *Xanthomonas arboricola* pv. *pruni*. *J Plant Pathol.* 2010; 92(Supplement 1): S1.99–S1.103.

5. Palacio-Bielsa A, Cambra MA, Cubero J, Garita-Cambronero J, Roselló M, Lopez MM. La mancha bacteriana de los frutales de hueso y del almendro (*Xanthomonas arboricola* pv. *pruni*), una grave enfermedad emergente en España. *Phytoma-España*. 2014; 259: 36–42.

6. Palacio-Bielsa A, Lopez MM, Peñalver J, Morente C, Cubero J, Garita-Cambronero J, et al. La mancha bacteriana de los frutales de hueso y del almendro (*Xanthomonas arboricola* pv. *pruni*) en España y Sudamérica. *Phytoma-España*. 2015; 271: 21–28.

7. Smith E. Observations on a hitherto unreported bacterial disease, the cause of which enters the plant through ordinary stomata. *Science*. 1903; 17: 456–457.

8. Roselló M, Santiago R, Palacio-Bielsa A, García-Figueres F, Montón C, Cambra MA, et al. Current status of bacterial spot of stone fruits and almond caused by *Xanthomonas arboricola* pv. *pruni* in Spain. *J Plant Pathol.* 2012; 94(Supplement 1): S1.15–S1.21.
9. Battilani P, Rossi V, Saccardi A. Development of *Xanthomonas arboricola* pv. *pruni* epidemics on peaches. J Plant Pathol. 1999; 81: 161–171.

10. López MM, Peñalver J, Morente MC, Quesada JM, Navarro I, López-Soriano P, et al. Evaluation of the efficiency of a conventional PCR protocol for the diagnosis of bacterial spot disease caused by *Xanthomonas arboricola* pv. *pruni*. in stone fruits and almond. J Plant Pathol. 2012; 94(Supplement 1): S1.75–S1.82.

11. Anonymous. EPPO standards PM 7/64. Diagnostics *Xanthomonas arboricola* pv. *pruni*. Bull OEPP/EPPO Bull. 2006; 36: 129–133.

12. Dhavantari BN. Overwintering sources of inoculum of bacterial spot of peach (*Xanthomonas pruni*) in southwestern Ontario. Proc Can Phytopathol Soc. 1971; 37: 21–30.

13. Zaccardelli M, Malaguti S, Bazzi C. Biological and epidemiological aspects of *Xanthomonas arboricola* pv. *pruni* on peach in Italy. J Plant Pathol. 1998; 80: 125–132.

14. Shepard DP, Zehr EI. Epiphytic persistence of *Xanthomonas campestris* pv. *pruni* on peach and plum. Plant Dis. 1994; 78: 627–629.

15. Goodman CA, Hattingh MJ. Transmission of *Xanthomonas campestris* pv. *pruni* in plum and apricot nursery trees by budding. HortSci. 1986; 21: 995–996.

16. López MM, Roselló M, Palacio-Bielsa A. Diagnosis and detection of the main bacterial pathogens of stone fruit and almond. J Plant Pathol. 2010; 92(Supplement 1): S1.57–S1.66.

17. Boudon S, Manceau C, Nottéghem JL. Structure and origin of *Xanthomonas arboricola* pv. *pruni* populations causing bacterial spot of stone fruit trees in Western Europe. Phytopathology. 2005; 95: 1081–1088. doi: 10.1094/PHYTO-95-1081 PMID: 18943306

18. Achtman M. Evolution, population structure, and phylogeography of genetically monomorphic bacterial pathogens. Annu Rev Microbiol. 2008; 62: 53–70. doi: 10.1146/annurev.micro.62.081307.162832 PMID: 18785837

19. Fischer-Le Saux M, Bonneau S, Essakhi S, Manceau C, Jacques MA. Aggressive emerging pathovars of *Xanthomonas arboricola* represent widespread epidemic clones that are distinct from poorly pathogenic strains, as revealed by multilocus sequence typing. Appl Environ Microbiol. 2015; 81: 4651–4668. doi: 10.1128/AEM.00050-15 PMID: 25934623

20. Barionovi D, Scortichini M. Integron variability in *Xanthomonas arboricola* pv. *juglandis* and *Xanthomonas arboricola* pv. *pruni* strains. FEMS Microbiol Lett. 2008; 288: 19–24. doi: 10.1111/j.1574-6968.2008.01315.x PMID: 18785932

21. Kawaguchi A. Genetic diversity of *Xanthomonas arboricola* pv. *pruni* strains in Japan revealed by DNA fingerprinting. J Gen Plant Pathol. 2014; 80: 366–369.

22. Lindstedt BA. Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis. 2005; 26: 2567–2582. PMID: 15937984

23. Van Belkum A, Scherer S, Van Alphen L, Verbrugh H. Short-sequence DNA repeats in prokaryotic genomes. Microbiol Mol Biol Rev. 1998; 62: 275–293. PMID: 9618442

24. Levinson G, Gutman GA. Slipped-strand mispairing: A major mechanism for DNA sequence evolution. Mol Biol Evol. 1987; 4: 203–221. PMID: 3328815

25. Balloux F, Lugon-Moulin N. The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002; 11: 155–165. PMID: 11856418

26. Vogler AJ, Keys C, Nemoto Y, Colman RE, Jay Z, Keim P. Effect of repeat copy number on Variable-Number Tandem Repeat mutations in *Escherichia coli* O157:H7. J Bacteriol. 2006; 188: 4253–4263. PMID: 16740932

27. Coletta-Filho HD, Takita MA, Alves de Souza A, Aguilar-Vildoso CI, Machado MA. Differentiation of strains of *Xylella fastidiosa* by a variable number of tandem repeat analysis. Appl Environ Microbiol. 2001; 67: 4091–4095. PMID: 11526010

28. Bui Thi Ngoc L, Vernière C, Jarne P, Brisse S, Guérin F, Boutry S, et al. From local surveys to global surveillance: three high-throughput genotyping methods for epidemiological monitoring of *Xanthomonas citri* pv. *citri* pathotypes. Appl Environ Microbiol. 2009; 75: 1173–1184. doi: 10.1128/AEM.02245-08 PMID: 19088309

29. Pruvost O, Magne M, Boyer K, Leduc A, Tourterel C, Drevet C, et al. 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 One. 2014; 9(6): e98129. doi: 10.1371/journal.pone.0098129 PMID: 24897119

30. Zhao S, Poulin L, Rodriguez-R LM, Serna NF, Liu SY, Wonni I, et al. Development of a variable number of tandem repeats typing scheme for the bacterial rice pathogen *Xanthomonas oryzae* pv. *oryzicola*. Phytopathology. 2012; 102: 948–956. doi: 10.1094/PHYTO-04-12-0078-R PMID: 22957820
31. Poulin L, Grygiel P, Magne M, Gagnevin L, Rodriguez-R LM, Forero Serna N, et al. 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 Microbiol*. 2015; 81: 688–698. doi: 10.1128/AEM.02768-14 PMID: 25398857

32. Arrieta-Ortiz ML, Rodríguez-R LM, Pérez-Quintero AL, Poulin L, Díaz AC, Arias Rojas N, et al. Genomic survey of pathogenicity determinants and VNTR markers in the cassava bacterial pathogen *Xanthomonas axonopodis* pv. *manihotis* strain CIO151. PLoS One. 2013; 8(11): e79704. doi: 10.1371/journal.pone.0079704 PMID: 24278159

33. Cesbron S, Pothier J, Gironde S, Jacques MA, Manceau C. Development of multilocus variable-number tandem repeat analysis (MLVA) for *Xanthomonas arboricola* pathovars. J Microbiol Meth. 2014; 100: 84–90.

34. Essakhi S, Cesbron S, Fischer-Le Saux M, Bonneau S, Jacques MA, Manceau C. Phylogenetic and variable-number tandem-repeat analyses identify nonpathogenic *Xanthomonas arboricola* lineages lacking the canonical type III secretion system. Appl Environ Microbiol. 2015; 81: 5395–5410. doi: 10.1128/AEM.00835-15 PMID: 26048944

35. Bergsma-Vlami M, Martin W, Koenraadt H, Teunissen H, Pothier JF, Duffy B, et al. Molecular typing of Dutch isolates of *Xanthomonas arboricola* pv. *pruni* isolated from ornamental cherry laurel. J Plant Pathol. 2012; 94(Supplement 1): S1.29–S1.35.

36. Ridé M. Bactéries phytopathogènes et maladies bactériennes des végétaux. In: Bourgin CVM editor. Les bactérioses et les viroses des arbres fruitiers. Ponsot, Paris, France; 1969.

37. Llop P, Caruso P, Cubero J, Morente C, López MM. A simple extraction procedure for efficient routine detection of pathogenic bacteria in plant material by polymerase chain reaction. J Microbiol Meth. 1999; 37: 23–31.

38. Li Y, Hu Y, Bolund L, Wang J. State of the art de novo assembly of human genomes from massively parallel sequencing data. Hum Genomics. 2010; 4: 271–277. PMID: 20511140

39. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using De Bruijn graphs. Genome Res. 2008; 18: 821–829. doi: 10.1101/gr.074492.107 PMID: 18349386

40. Sallet E, Gouzy J, Schiex T. EuGene-PP: a next generation automated annotation pipeline for prokaryotic genomes. Bioinformatics. 2014; 30: 2659–2661. doi: 10.1093/bioinformatics/btu366 PMID: 24880686

41. Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999; 27: 573–580. PMID: 9862982

42. Denoeud F, Vergnaud G. Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource. BMC Bioinformatics. 2004; 5: 1–12.

43. Pourcel C, Vergnaud G. Strain typing using Multiple “Variable Number of Tandem Repeat” Analysis and genetic element CRISPR. In: Persing DH, Tenover FC, Tang YW, Nolte FS, Hayden RT, et al., editors. Molecular microbiology: Diagnostic principles and practice. 2nd ed. ASM Press, Washington, DC; 2011. pp. 179–197.

44. Francisco AP, Vaz C, Monteiro PT, Melo-Cristino J, Ramirez M, Carriço JA. PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods. BMC Bioinformatics. 2012; 13: 87. doi: 10.1186/1471-2105-13-87 PMID: 22568821

45. Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004; 186: 1518–1530. PMID: 14973027

46. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinf. 2005; 1: 47–50.

47. Kalinowski ST. HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes. 2005; 5: 187–189.

48. Abdi H. Metric multidimensional scaling (MDS): Analyzing distance matrices. In: Salkind N, editor. Encyclopedia of measurement and statistics. Sage, Thousand Oaks, CA; 2007. pp. 1–13.

49. Estoup A, Guillemaud T. Reconstructing routes of invasion using genetic data: why, how and so what? Mol Ecol. 2010; 19: 4113–4130. doi: 10.1111/j.1365-294X.2010.04773.x PMID: 20723048

50. Leduc A, Traoré YN, Boyer K, Magne M, Grygiel P, Juhasz CC, et al. Bridgehead invasion of a monomorphic plant pathogenic bacterium: *Xanthomonas citri* pv. *citri*, an emerging citrus pathogen in Mali and Burkina Faso. Environ Microbiol. 2015; 17: 4429–4442. doi: 10.1111/1462-2920.12876 PMID: 25866121

51. Duglosch KM, Parker IM. Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol. 2008; 17: 431–449. PMID: 17908213
52. Donat V, Biosca EG, Peñalver J, López MM. Exploring diversity among Spanish strains of *Erwinia amylovora* and possible infection sources. J Appl Microbiol. 2007; 103: 1639–1649. PMID: 17953575

53. Bülhmann A, Dreo T, Rezzonico F, Pothier JF, Smits THM, Ravnikar M, et al. Phylogeography and population structure of the biologically invasive phytopathogen *Erwinia amylovora* inferred using minisatellites. Environ Microbiol. 2014; 16: 2112–2125. doi: 10.1111/1462-2920.12289 PMID: 24112873

54. Ritchie DF. Bacterial spot. In: Ogawa JM, Zehr EI, Bird GW, Ritchie DF, Uriu K, Uyemoto JK, editors. Compendium of stone fruit diseases. American Phytopathological Society, Saint Paul, MN, USA; 1995. pp. 50–52.

55. Zehr EI, Shepard DP. Bacterial spot of peach as influenced by water congestion, leaf wetness duration, and temperature. Plant Dis. 1996; 80: 339–341.

56. Irey M, Gottwald TR, Graham JH, Riley TD, Carlton G. Post-hurricane analysis of citrus canker spread and progress towards the development of a predictive model to estimate disease spread due to catastrophic weather events. Plant Health Prog. 2006; doi: 10.1094/PHP-2006-0822-01-RS

57. Gasperini C, Bazzi C, Mazzucchi U. Autumn inoculation of *Xanthomonas campestris* pv. *pruni* through leaf scars in plum trees in the Po valley. Phytopathol Mediterr. 1984; 23: 60–62.

58. Loreti S, Pucci N, Perez G, Catara V, Scortichini M, Bella P, et al. Detection and identification of *Xanthomonas arboricola* pv. *pruni* from symptomless plant material: results of an Italian test performance study. Bull OEPP/EPPO Bull. 2015; 45: 41–51.

59. Palacio-Bielsa A, Cubero J, Cambra MA, Collados R, Berruete IM, López MM. Development of an efficient real-time quantitative PCR protocol for detection of *Xanthomonas arboricola* pv. *pruni* in *Prunus* species. Appl Environ Microbiol. 2011; 77: 89–97. doi: 10.1128/AEM.01593-10 PMID: 21037298

60. Palacio-Bielsa A, Pothier JF, Roselló M, Duffy B, López MM. Detection and identification methods and new tests as developed and used in the framework of COST 873 for bacteria pathogenic to stone fruits and nuts. J Plant Pathol. 2012; 94(Supplement 1): S1.135–S1.146.