From local surveys to global surveillance: three high-throughput genotyping methods for epidemiological monitoring of Xanthomonas citri pv. citri pathotypes

Genetic Structure and TALome Analysis Highlight a High Level of Diversity in Burkinabe Xanthomonas Oryzae pv. oryzae Populations

Bacterial Leaf Blight of rice (BLB) caused by Xanthomonas oryzae pv. oryzae (Xoo) is a major threat for food security in many rice growing countries including Burkina Faso, where the disease was first reported in the 1980's. In line with the intensification of rice cultivation in West-Africa, BLB incidence has been rising for the last 15 years. West-African strains of Xoo differ from their Asian counterparts as they (i) are genetically distant, (ii) belong to new races and, (iii) contain reduced repertoires of Transcription Activator Like (TAL) effector genes. In order to investigate the evolutionary dynamics of Xoo populations in Burkina Faso, 177 strains were collected from 2003 to 2018 in three regions where BLB is occurring. Multilocus VNTR Analysis (MLVA-14) targeting 10 polymorphic loci discriminated 24 haplotypes and showed that Xoo populations were structured according to their geographical localization and year of collection. Considering their major role in Xoo pathogenicity, we assessed the TAL effector repertoires of the 177 strains upon RFLP-based profiling. Surprisingly, an important diversity was revealed with up to eight different RFLP patterns. Finally, comparing neutral vs. tal effector gene diversity allowed to suggest scenarios underlying the evolutionary dynamics of Xoo populations in Burkina Faso, which is key to rationally guide the deployment of durably resistant rice varieties against BLB in the country.

Citrus Canker Pathogen, Its Mechanism of Infection, Eradication, and Impacts

Citrus canker is a ravaging bacterial disease threatening citrus crops. Its major types are Asiatic Canker, Cancrosis B, and Cancrosis C, caused by Xanthomonas citri pv. citri (Xcc), Xanthomonas citri pv. aurantifolii pathotype-B (XauB), and pathotype-C (XauC), respectively. The bacterium enters its host through stomata and wounds, from which it invades the intercellular spaces in the apoplast. It produces erumpent corky necrotic lesions often surrounded by a chlorotic halo on the leaves, young stems, and fruits, which causes dark spots, defoliation, reduced photosynthetic rate, rupture of leaf epidermis, dieback, and premature fruit drop in severe cases. Its main pathogenicity determinant gene is pthA, whose variants are present in all citrus canker-causing pathogens. Countries where citrus canker is not endemic adopt different methods to prevent the introduction of the pathogen into the region, eradicate the pathogen, and minimize its dissemination, whereas endemic regions require an integrated management program to control the disease. The main aim of the present manuscript is to shed light on the pathogen profile, its mechanism of infection, and fruitful strategies for disease management. Although an adequate method to completely eradicate citrus canker has not been introduced so far, many new methods are under research to abate the disease.

Molecular Epidemiology of Xanthomonas euvesicatoria Strains from the Balkan Peninsula Revealed by a New Multiple-Locus Variable-Number Tandem-Repeat Analysis Scheme

Bacterial spot of pepper and tomato is caused by at least three species of Xanthomonas, among them two pathovars of Xanthomonas euvesicatoria, which are responsible for significant yield losses on all continents. In order to trace back the spread of bacterial spot pathogens within and among countries, we developed the first multilocus variable number of tandem repeat analyses (MLVA) scheme for pepper- and tomato-pathogenic strains of X. euvesicatoria. In this work, we assessed the repeat numbers by DNA sequencing of 16 tandem repeat loci and applied this new tool to analyse a representative set of 88 X. euvesicatoria pepper strains from Bulgaria and North Macedonia. The MLVA-16 scheme resulted in a Hunter–Gaston Discriminatory Index (HGDI) score of 0.944 and allowed to resolve 36 MLVA haplotypes (MTs), thus demonstrating its suitability for high-resolution molecular typing. Strains from the different regions of Bulgaria and North Macedonia were found to be widespread in genetically distant clonal complexes or singletons. Sequence types of the variable number of tandem repeats (VNTR) amplicons revealed cases of size homoplasy and suggested the coexistence of different populations and different introduction events. The large geographical distribution of MTs and the existence of epidemiologically closely related strains in different regions and countries suggest long dispersal of strains on pepper in this area.

Contrasting genetic diversity and structure among Malagasy Ralstonia pseudosolanacearum phylotype I populations inferred from an optimized Multilocus Variable Number of Tandem Repeat Analysis scheme

The Ralstonia solanacearum species complex (RSSC), composed of three species and four phylotypes, are globally distributed soil-borne bacteria with a very broad host range. In 2009, a devastating potato bacterial wilt outbreak was declared in the central highlands of Madagascar, which reduced the production of vegetable crops including potato, eggplant, tomato and pepper. A molecular epidemiology study of Malagasy RSSC strains carried out between 2013 and 2017 identified R. pseudosolanacearum (phylotypes I and III) and R. solanacearum (phylotype II). A previously published population biology analysis of phylotypes II and III using two MultiLocus Variable Number of Tandem Repeats Analysis (MLVA) schemes revealed an emergent epidemic phylotype II (sequevar 1) group and endemic phylotype III isolates. We developed an optimized MLVA scheme (RS1-MLVA14) to characterize phylotype I strains in Madagascar to understand their genetic diversity and structure. The collection included isolates from 16 fields of different Solanaceae species sampled in Analamanga and Itasy regions (highlands) in 2013 (123 strains) and in Atsinanana region (lowlands) in 2006 (25 strains). Thirty-one haplotypes were identified, two of them being particularly prevalent: MT007 (30.14%) and MT004 (16.44%) (sequevar 18). Genetic diversity analysis revealed a significant contrasting level of diversity according to elevation and sampling region. More diverse at low altitude than at high altitude, the Malagasy phylotype I isolates were structured in two clusters, probably resulting from different historical introductions. Interestingly, the most prevalent Malagasy phylotype I isolates were genetically distant from regional and worldwide isolates. In this work, we demonstrated that the RS1-MLVA14 scheme can resolve differences from regional to field scales and is thus suited for deciphering the epidemiology of phylotype I populations.

Development and comparative validation of genomic-driven PCR-based assays to detect Xanthomonas citri pv. citri in citrus plants

Background

Asiatic Citrus Canker, caused by Xanthomonas citri pv. citri, severely impacts citrus production worldwide and hampers international trade. Considerable regulatory procedures have been implemented to prevent the introduction and establishment of X. citri pv. citri into areas where it is not present. The effectiveness of this surveillance largely relies on the availability of specific and sensitive detection protocols. Although several PCR- or real-time PCR-based methods are available, most of them showed analytical specificity issues. Therefore, we developed new conventional and real-time quantitative PCR assays, which target a region identified by comparative genomic analyses, and compared them to existing protocols.

Results

Our assays target the X. citri pv. citri XAC1051 gene that encodes for a putative transmembrane protein. The real-time PCR assay includes an internal plant control (5.8S rDNA) for validating the assay in the absence of target amplification. A receiver-operating characteristic approach was used in order to determine a reliable cycle cut-off for providing accurate qualitative results. Repeatability, reproducibility and transferability between real-time devices were demonstrated for this duplex qPCR assay (XAC1051-2qPCR). When challenged with an extensive collection of target and non-target strains, both assays displayed a high analytical sensitivity and specificity performance: LOD_95% = 754 CFU ml^− 1 (15 cells per reaction), 100% inclusivity, 97.2% exclusivity for XAC1051-2qPCR; LOD_95% = 5234 CFU ml^− 1 (105 cells per reaction), 100% exclusivity and inclusivity for the conventional PCR. Both assays can detect the target from naturally infected citrus fruit. Interestingly, XAC1051-2qPCR detected X. citri pv. citri from herbarium citrus samples. The new PCR-based assays displayed enhanced analytical sensitivity and specificity when compared with previously published PCR and real-time qPCR assays.

Conclusions

We developed new valuable detection assays useful for routine diagnostics and surveillance of X. citri pv. citri in citrus material. Their reliability was evidenced through numerous trials on a wide range of bacterial strains and plant samples. Successful detection of the pathogen was achieved from both artificially and naturally infected plants, as well as from citrus herbarium samples, suggesting that these assays will have positive impact both for future applied and academic research on this bacterium.

CRISPR elements provide a new framework for the genealogy of the citrus canker pathogen Xanthomonas citri pv. citri

Background

Xanthomonads are an important clade of Gram-negative bacteria infecting a plethora of economically important host plants, including citrus. Knowledge about the pathogen’s diversity and population structure are prerequisite for epidemiological surveillance and efficient disease management. Rapidly evolving genetic loci, such as Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), are of special interest to develop new molecular typing tools.

Results

We analyzed CRISPR loci of 56 Xanthomonas citri pv. citri strains of world-wide origin, a regulated pathogen causing Asiatic citrus canker in several regions of the world. With one exception, 23 unique sequences built up the repertoire of spacers, suggesting that this set of strains originated from a common ancestor that already harbored these 23 spacers. One isolate originating from Pakistan contained a string of 14 additional, probably more recently acquired spacers indicating that this genetic lineage has or had until recently the capacity to acquire new spacers. Comparison of CRISPR arrays with previously obtained molecular typing data, such as amplified fragment length polymorphisms (AFLP), variable-number of tandem-repeats (VNTR) and genome-wide single-nucleotide polymorphisms (SNP), demonstrated that these methods reveal similar evolutionary trajectories. Notably, genome analyses allowed to generate a model for CRISPR array evolution in X. citri pv. citri, which provides a new framework for the genealogy of the citrus canker pathogen.

Conclusions

CRISPR-based typing will further improve the accuracy of the genetic identification of X. citri pv. citri outbreak strains in molecular epidemiology analyses, especially when used concomitantly with another genotyping method.

Detection and identification of Xanthomonas pathotypes associated with citrus diseases using comparative genomics and multiplex PCR

Background

In Citrus cultures, three species of Xanthomonas are known to cause distinct diseases. X. citri subsp. citri patothype A, X. fuscans subsp. aurantifolii pathotypes B and C, and X. alfalfae subsp. citrumelonis, are the causative agents of cancrosis A, B, C, and citrus bacterial spots, respectively. Although these species exhibit different levels of virulence and aggressiveness, only limited alternatives are currently available for proper and early detection of these diseases in the fields. The present study aimed to develop a new molecular diagnostic method based on genomic sequences derived from the four species of Xanthomonas.

Results

Using comparative genomics approaches, primers were synthesized for the identification of the four causative agents of citrus diseases. These primers were validated for their specificity to their target DNA by both conventional and multiplex PCR. Upon evaluation, their sensitivity was found to be 0.02 ng/µl in vitro and 1.5 × 10⁴ CFU ml⁻¹ in infected leaves. Additionally, none of the primers were able to generate amplicons in 19 other genomes of Xanthomonas not associated with Citrus and one species of Xylella, the causal agent of citrus variegated chlorosis (CVC). This denotes strong specificity of the primers for the different species of Xanthomonas investigated in this study.

Conclusions

We demonstrated that these markers can be used as potential candidates for performing in vivo molecular diagnosis exclusively for citrus-associated Xanthomonas. The bioinformatics pipeline developed in this study to design specific genomic regions is capable of generating specific primers. It is freely available and can be utilized for any other model organism.

Deciphering how plant pathogenic bacteria disperse and meet: Molecular epidemiology of Xanthomonas citri pv. citri at microgeographic scales in a tropical area of Asiatic citrus canker endemicity

Although some plant pathogenic bacteria represent a significant threat to agriculture, the determinants of their ecological success and evolutionary potential are still poorly understood. Refining our understanding of bacterial strain circulation at small spatial scales and the biological significance and evolutionary consequences of co-infections are key questions. The study of bacterial population biology can be challenging, because it requires high-resolution markers that can be genotyped with a high throughput. Here, we overcame this difficulty for Xanthomonas citri pv. citri, a genetically monomorphic bacterium causing Asiatic citrus canker (ACC). Using a genotyping method that did not require cultivating the bacterium or purifying DNA, we deciphered the pathogen's spatial genetic structure at several microgeographic scales, down to single lesion, in a situation of ACC endemicity. In a grove where copper was recurrently applied for ACC management, copper-susceptible and copper-resistant X. citri pv. citri coexisted and the bacterial population structured as three genetic clusters, suggesting a polyclonal contamination. The range of spatial dependency, estimated for the two largest clusters, was four times greater for the cluster predominantly composed of copper-resistant bacteria. Consistently, the evenness value calculated for this cluster was indicative of increased transmission. Linkage disequilibrium was high even at a tree scale, probably due to a combination of clonality and admixture. Approximately 1% of samples exhibited within-lesion multilocus polymorphism, explained at least in part by polyclonal infections. Canker lesions, which are of major biological significance as an inoculum source, may also represent a preferred niche for horizontal gene transfer. This study points out the potential of genotyping data for estimating the range of spatial dependency of plant bacterial pathogens, an important parameter for guiding disease management strategies.

An Optimized Microsatellite Scheme for Assessing Populations of Xanthomonas phaseoli pv. manihotis

Diverse molecular markers have been used to analyze the genetic diversity of plant pathogens. Compared with traditional fingerprinting methods, multiple loci variable number of tandem repeat analyses (MLVAs) have gained importance recently due to their reproducibility, high discriminatory power, ease of performance, low cost, and throughput potential. These characteristics are desirable for continuous pathogen monitoring, especially for pathogens with relatively low genetic diversity, and for disease epidemiology studies. Genetic diversity studies of Xanthomonas phaseoli pv. manihotis, which is the causal agent of cassava bacterial blight, have shown variability and changes in the bacterial population over time. Thus, an easy and fast method needs to be developed to type populations of this pathogen in different countries of the world, especially on small scales. In this study, we developed an MLVA scheme to analyze X. phaseoli pv. manihotis variability on a local scale. The MLVA-15 scheme comprises 15 variable number of tandem repeat loci grouped into four multiplex polymerase chain reaction pools. We showed that the MLVA-15 scheme had slightly higher discriminatory ability at the locality level when compared with amplified fragment length polymorphisms. The MLVA-15 scheme allowed for an accurate determination of the number of genotypes in the sample and showed reproducibility and portability. Additionally, this scheme could be used to analyze numerous strains in a reasonable timeframe. The MLVA-15 scheme was highly specific to X. phaseoli but up to eight tandem repeat loci could be amplified from other Xanthomonas spp. Finally, we assessed the utility of the scheme for analyses of X. phaseoli pv. manihotis genetic variability in the Colombian Caribbean region. MLVA-15 distinguished 88.9% of the haplotypes in our sample. Strains originating from the same field and isolated at the same time could be discriminated. In this study, the advantages of the MLVA-15 scheme targeting 6- or 7-bp repeats were demonstrated. Moreover, this scheme was a fast method that was appropriate for routine monitoring of X. phaseoli pv. manihotis populations on a local scale and, thus, was useful for addressing epidemiological questions.

A new Multi Locus Variable Number of Tandem Repeat Analysis Scheme for epidemiological surveillance of Xanthomonas vasicola pv. musacearum, the plant pathogen causing bacterial wilt on banana and enset

Xanthomonas vasicola pv. musacearum (Xvm) which causes Xanthomonas wilt (XW) on banana (Musa accuminata x balbisiana) and enset (Ensete ventricosum), is closely related to the species Xanthomonas vasicola that contains the pathovars vasculorum (Xvv) and holcicola (Xvh), respectively pathogenic to sugarcane and sorghum. Xvm is considered a monomorphic bacterium whose intra-pathovar diversity remains poorly understood. With the sudden emergence of Xvm within east and central Africa coupled with the unknown origin of one of the two sublineages suggested for Xvm, attention has shifted to adapting technologies that focus on identifying the origin and distribution of the genetic diversity within this pathogen. Although microbiological and conventional molecular diagnostics have been useful in pathogen identification. Recent advances have ushered in an era of genomic epidemiology that aids in characterizing monomorphic pathogens. To unravel the origin and pathways of the recent emergence of XW in Eastern and Central Africa, there was a need for a genotyping tool adapted for molecular epidemiology. Multi-Locus Variable Number of Tandem Repeat Analysis (MLVA) is able to resolve the evolutionary patterns and invasion routes of a pathogen. In this study, we identified microsatellite loci from nine published Xvm genome sequences. Of the 36 detected microsatellite loci, 21 were selected for primer design and 19 determined to be highly typeable, specific, reproducible and polymorphic with two- to four- alleles per locus on a sub-collection. The 19 markers were multiplexed and applied to genotype 335 Xvm strains isolated from seven countries over several years. The microsatellite markers grouped the Xvm collection into three clusters; with two similar to the SNP-based sublineages 1 and 2 and a new cluster 3, revealing an unknown diversity in Ethiopia. Five of the 19 markers had alleles present in both Xvm and Xanthomonas vasicola pathovars holcicola and vasculorum, supporting the phylogenetic closeliness of these three pathovars. Thank to the public availability of the haplotypes on the MLVABank database, this highly reliable and polymorphic genotyping tool can be further used in a transnational surveillance network to monitor the spread and evolution of XW throughout Africa.. It will inform and guide management of Xvm both in banana-based and enset-based cropping systems. Due to the suitability of MLVA-19 markers for population genetic analyses, this genotyping tool will also be used in future microevolution studies.

Molecular Epidemiology of Bacterial Wilt in the Madagascar Highlands Caused by Andean (Phylotype IIB-1) and African (Phylotype III) Brown Rot Strains of the Ralstonia solanacearum Species Complex

The Ralstonia solanacearum species complex (RSSC) is a highly diverse cluster of bacterial strains found worldwide, many of which are destructive and cause bacterial wilt (BW) in a wide range of host plants. In 2009, potato production in Madagascar was dramatically affected by several BW epidemics. Controlling this disease is critical for Malagasy potato producers. The first important step toward control is the characterization of strains and their putative origins. The genetic diversity and population structure of the RSSC were investigated in the major potato production areas of the Highlands. A large collection of strains (n = 1224) was assigned to RSSC phylotypes based on multiplex polymerase chain reaction (PCR). Phylotypes I and III have been present in Madagascar for a long time but rarely associated with major potato BW outbreaks. The marked increase of BW prevalence was found associated with phylotype IIB sequevar 1 (IIB-1) strains (n = 879). This is the first report of phylotype IIB-1 strains in Madagascar. In addition to reference strains, epidemic IIB-1 strains (n = 255) were genotyped using the existing MultiLocus Variable-Number Tandem Repeat Analysis (MLVA) scheme RS2-MLVA9, producing 31 haplotypes separated into two related clonal complexes (CCs). One major CC included most of the worldwide haplotypes distributed across wide areas. A regional-scale investigation suggested that phylotype IIB-1 strains were introduced and massively spread via latently infected potato seed tubers. Additionally, the genetic structure of phylotype IIB-1 likely resulted from a bottleneck/founder effect. The population structure of phylotype III, described here for the first time in Madagascar, exhibited a different pattern. Phylotype III strains (n = 217) were genotyped using the highly discriminatory MLVA scheme RS3-MLVA16. High genetic diversity was uncovered, with 117 haplotypes grouped into 11 CCs. Malagasy phylotype III strains were highly differentiated from continental African strains, suggesting no recent migration from the continent. Overall, population structure of phylotype III involves individual small CCs that correlate to restricted geographic areas in Madagascar. The evidence suggests, if at all, that African phylotype III strains are not efficiently transmitted through latently infected potato seed tubers.

Multiple Locus Variable number of tandem repeat Analysis: A molecular genotyping tool for Paenibacillus larvae

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P. larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P. larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated.

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

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.

Comparative genomics of 43 strains of Xanthomonas citri pv. citri reveals the evolutionary events giving rise to pathotypes with different host ranges

BACKGROUND

The identification of factors involved in the host range definition and evolution is a pivotal challenge in the goal to predict and prevent the emergence of plant bacterial disease. To trace the evolution and find molecular differences between three pathotypes of Xanthomonas citri pv. citri that may explain their distinctive host ranges, 42 strains of X. citri pv. citri and one outgroup strain, Xanthomonas citri pv. bilvae were sequenced and compared.

RESULTS

The strains from each pathotype form monophyletic clades, with a short branch shared by the A(w) and A pathotypes. Pathotype-specific recombination was detected in seven regions of the alignment. Using Ancestral Character Estimation, 426 SNPs were mapped to the four branches at the base of the A, A*, A(w) and A/A(w) clades. Several genes containing pathotype-specific nonsynonymous mutations have functions related to pathogenicity. The A pathotype is enriched for SNP-containing genes involved in defense mechanisms, while A* is significantly depleted for genes that are involved in transcription. The pathotypes differ by four gene islands that largely coincide with regions of recombination and include genes with a role in virulence. Both A* and A(w) are missing genes involved in defense mechanisms. In contrast to a recent study, we find that there are an extremely small number of pathotype-specific gene presences and absences.

CONCLUSIONS

The three pathotypes of X. citri pv. citri that differ in their host ranges largely show genomic differences related to recombination, horizontal gene transfer and single nucleotide polymorphism. We detail the phylogenetic relationship of the pathotypes and provide a set of candidate genes involved in pathotype-specific evolutionary events that could explain to the differences in host range and pathogenicity between them.

Genomics and transcriptomics of Xanthomonas campestris species challenge the concept of core type III effectome

BACKGROUND

The bacterial species Xanthomonas campestris infects a wide range of Brassicaceae. Specific pathovars of this species cause black rot (pv. campestris), bacterial blight of stock (pv. incanae) or bacterial leaf spot (pv. raphani).

RESULTS

In this study, we extended the genomic coverage of the species by sequencing and annotating the genomes of strains from pathovar incanae (CFBP 1606R and CFBP 2527R), pathovar raphani (CFBP 5828R) and a pathovar formerly named barbareae (CFBP 5825R). While comparative analyses identified a large core ORFeome at the species level, the core type III effectome was limited to only three putative type III effectors (XopP, XopF1 and XopAL1). In Xanthomonas, these effector proteins are injected inside the plant cells by the type III secretion system and contribute collectively to virulence. A deep and strand-specific RNA sequencing strategy was adopted in order to experimentally refine genome annotation for strain CFBP 5828R. This approach also allowed the experimental definition of novel ORFs and non-coding RNA transcripts. Using a constitutively active allele of hrpG, a master regulator of the type III secretion system, a HrpG-dependent regulon of 141 genes co-regulated with the type III secretion system was identified. Importantly, all these genes but seven are positively regulated by HrpG and 56 of those encode components of the Hrp type III secretion system and putative effector proteins.

CONCLUSIONS

This dataset is an important resource to mine for novel type III effector proteins as well as for bacterial genes which could contribute to pathogenicity of X. campestris.

Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA(w)) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA(w) and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in 'carbohydrate transport and metabolism' and 'DNA replication, recombination and repair' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.

Origin of the Outbreak in France of Pseudomonas syringae pv. actinidiae Biovar 3, the Causal Agent of Bacterial Canker of Kiwifruit, Revealed by a Multilocus Variable-Number Tandem-Repeat Analysis

The first outbreaks of bacterial canker of kiwifruit caused by Pseudomonas syringae pv. actinidiae biovar 3 were detected in France in 2010. P. syringae pv. actinidiae causes leaf spots, dieback, and canker that sometimes lead to the death of the vine. P. syringae pv. actinidifoliorum, which is pathogenic on kiwi as well, causes only leaf spots. In order to conduct an epidemiological study to track the spread of the epidemics of these two pathogens in France, we developed a multilocus variable-number tandem-repeat (VNTR) analysis (MLVA). MLVA was conducted on 340 strains of P. syringae pv. actinidiae biovar 3 isolated in Chile, China, France, Italy, and New Zealand and on 39 strains of P. syringae pv. actinidifoliorum isolated in Australia, France, and New Zealand. Eleven polymorphic VNTR loci were identified in the genomes of P. syringae pv. actinidiae biovar 3 ICMP 18744 and of P. syringae pv. actinidifoliorum ICMP 18807. MLVA enabled the structuring of P. syringae pv. actinidiae biovar 3 and P. syringae pv. actinidifoliorum strains in 55 and 16 haplotypes, respectively. MLVA and discriminant analysis of principal components revealed that strains isolated in Chile, China, and New Zealand are genetically distinct from P. syringae pv. actinidiae strains isolated in France and in Italy, which appear to be closely related at the genetic level. In contrast, no structuring was observed for P. syringae pv. actinidifoliorum. We developed an MLVA scheme to explore the diversity within P. syringae pv. actinidiae biovar 3 and to trace the dispersal routes of epidemic P. syringae pv. actinidiae biovar 3 in Europe. We suggest using this MLVA scheme to trace the dispersal routes of P. syringae pv. actinidiae at a global level.

Development of a Multiple Loci Variable Number of Tandem Repeats Analysis (MLVA) to Unravel the Intra-Pathovar Structure of Pseudomonas syringae pv. actinidiae Populations Worldwide

The bacterial canker of kiwifruit by Pseudomonas syringae pv. actinidiae is an emblematic example of a catastrophic disease of fruit crops. In 2008 a new, extremely virulent form of the pathogen emerged and rapidly devastated many Actinidia spp. orchards all over the world. In order to understand differences in populations within this pathovar and to elucidate their diffusion and movements on world scale, it is necessary to be able to quickly and on a routine basis compare new isolates with previous records. In this report a worldwide collection of 142 strains was analyzed by MLVA, chosen as investigative technique for its efficacy, reproducibility, simplicity and low cost. A panel of 13 Variable Number of Tandem Repeats (VNTR) loci was identified and used to describe the pathogen population. The MLVA clustering is highly congruent with the population structure as previously established by other molecular approaches including whole genome sequencing and correlates with geographic origin, time of isolation and virulence. For convenience, we divided the VNTR loci in two panels. Panel 1 assay, using six loci, recognizes 23 different haplotypes, clustered into ten complexes with highest congruence with previous classifications. Panel 2, with seven VNTR loci, provides discriminatory power. Using the total set of 13 VNTR loci, 58 haplotypes can be distinguished. The recent hypervirulent type shows very limited diversity and includes, beside the strains from Europe, New Zealand and Chile, a few strains from Shaanxi, China. A broad genetic variability is observed in China, but different types are also retrievable in Japan and Korea. The low virulent strains cluster together and are very different from the other MLVA genotypes. Data were used to generate a public database in MLVAbank. MLVA represents a very promising first-line assay for large-scale routine genotyping, prior to whole genome sequencing of only the most relevant samples.

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

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

Development of multilocus variable-number tandem repeat analysis (MLVA) for Xanthomonas arboricola pathovars

Xanthomonas arboricola is an important bacterial species, the pathovars of which are responsible for bacterial blight diseases on stone fruit, hazelnut, Persian walnut, poplar, strawberry, poinsettia and banana. In this study, we evaluated variable number tandem repeats (VNTR) as a molecular typing tool for assessing the genetic diversity within pathovars of X. arboricola. Screening of the X. arboricola pv. pruni genome sequence (CFBP5530 strain) predicted 51 candidate VNTR loci. Primer pairs for polymerase chain reaction (PCR) amplification of all 51 loci were designed, and their discriminatory power was initially evaluated with a core collection of 8 X. arboricola strains representative of the different pathovars. Next, the 26 polymorphic VNTR loci present in all strains were used for genotyping a collection of 61 strains. MLVA is a typing method that clearly differentiates X. arboricola strains. The MLVA scheme described in this study is a rapid and reliable molecular typing tool that can be used for further epidemiological studies of bacterial diseases caused by X. arboricola pathovars.

Highly polymorphic markers reveal the establishment of an invasive lineage of the citrus bacterial pathogen Xanthomonas citri pv. citri in its area of origin

Investigating the population biology of plant pathogens in their native areas is essential to understand the factors that shape their population structure and favour their spread. Monomorphic pathogens dispatch extremely low genetic diversity in invaded areas, and native areas constitute a major reservoir for future emerging strains. One of these, the gammaproteobacterium Xanthomonas citri pv. citri, causes Asiatic canker and is a considerable threat to citrus worldwide. We studied its population genetic structure by genotyping 555 strains from 12 Vietnam provinces at 14 tandem repeat loci and insertion sequences. Discriminant analysis of principal components identified six clusters. Five of them were composed of endemic strains distributed heterogeneously across sampled provinces. A sixth cluster, VN6, displayed a much lower diversity and a clonal expansion structure, suggesting recent epidemic spread. No differences in aggressiveness on citrus or resistance to bactericides were detected between VN6 and other strains. VN6 likely represents a case of bioinvasion following introduction in a native area likely through contaminated plant propagative material. Highly polymorphic markers are useful for revealing migration patterns of recently introduced populations of a monomorphic bacterial plant pathogen.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

Multilocus variable-number-tandem-repeats analysis (MLVA) distinguishes a clonal complex of Clavibacter michiganensis subsp. michiganensis strains isolated from recent outbreaks of bacterial wilt and canker in Belgium

BACKGROUND

Clavibacter michiganensis subsp. michiganensis (Cmm) causes bacterial wilt and canker in tomato. Cmm is present nearly in all European countries. During the last three years several local outbreaks were detected in Belgium. The lack of a convenient high-resolution strain-typing method has hampered the study of the routes of transmission of Cmm and epidemiology in tomato cultivation. In this study the genetic relatedness among a worldwide collection of Cmm strains and their relatives was approached by gyrB and dnaA gene sequencing. Further, we developed and applied a multilocus variable number of tandem repeats analysis (MLVA) scheme to discriminate among Cmm strains.

RESULTS

A phylogenetic analysis of gyrB and dnaA gene sequences of 56 Cmm strains demonstrated that Belgian Cmm strains from recent outbreaks of 2010-2012 form a genetically uniform group within the Cmm clade, and Cmm is phylogenetically distinct from other Clavibacter subspecies and from non-pathogenic Clavibacter-like strains. MLVA conducted with eight minisatellite loci detected 25 haplotypes within Cmm. All strains from Belgian outbreaks, isolated between 2010 and 2012, together with two French strains from 2010 seem to form one monomorphic group. Regardless of the isolation year, location or tomato cultivar, Belgian strains from recent outbreaks belonged to the same haplotype. On the contrary, strains from diverse geographical locations or isolated over longer periods of time formed mostly singletons.

CONCLUSIONS

We hypothesise that the introduction might have originated from one lot of seeds or contaminated tomato seedlings that was the source of the outbreak in 2010 and that these Cmm strains persisted and induced infection in 2011 and 2012. Our results demonstrate that MLVA is a promising typing technique for a local surveillance and outbreaks investigation in epidemiological studies of Cmm.

Variations in type III effector repertoires, pathological phenotypes and host range of Xanthomonas citri pv. citri pathotypes

The mechanisms determining the host range of Xanthomonas are still undeciphered, despite much interest in their potential roles in the evolution and emergence of plant pathogenic bacteria. Xanthomonas citri pv. citri (Xci) is an interesting model of host specialization because of its pathogenic variants: pathotype A strains infect a wide range of Rutaceous species, whereas pathotype A*/A(W) strains have a host range restricted to Mexican lime (Citrus aurantifolia) and alemow (Citrus macrophylla). Based on a collection of 55 strains representative of Xci worldwide diversity assessed by amplified fragment length polymorphism (AFLP), we investigated the distribution of type III effectors (T3Es) in relation to host range. We examined the presence of 66 T3Es from xanthomonads in Xci and identified a repertoire of 28 effectors, 26 of which were shared by all Xci strains, whereas two (xopAG and xopC1) were present only in some A*/A(W) strains. We found that xopAG (=avrGf1) was present in all A(W) strains, but also in three A* strains genetically distant from A(W) , and that all xopAG-containing strains induced the hypersensitive response (HR) on grapefruit and sweet orange. The analysis of xopAD and xopAG suggested horizontal transfer between X. citri pv. bilvae, another citrus pathogen, and some Xci strains. A strains were genetically less diverse, induced identical phenotypic responses and possessed indistinguishable T3E repertoires. Conversely, A*/A(W) strains exhibited a wider genetic diversity in which clades correlated with geographical origin and T3E repertoire, but not with pathogenicity, according to T3E deletion experiments. Our data outline the importance of taking into account the heterogeneity of Xci A*/A(W) strains when analysing the mechanisms of host specialization.

Characterization of a variant of Xanthomonas citri subsp. citri that triggers a host-specific defense response

Citrus is an economically important fruit crop that is severely afflicted by Asiatic citrus bacterial canker (CBC), a disease caused by the phytopathogen Xanthomonas citri subsp. citri (X. citri). To gain insight into the molecular epidemiology of CBC, 42 Xanthomonas isolates were collected from a range of Citrus spp. across 17 different orchards in Tucumán, Argentina and subjected to molecular, biochemical, and pathogenicity tests. Analysis of genome-specific X. citri markers and DNA polymorphisms based on repetitive elements-based polymerase chain reaction showed that all 42 isolates belonged to X. citri. Interestingly, pathogenicity tests showed that one isolate, which shares >90% genetic similarity to the reference strain X. citri T, has host range specificity. This new variant of X. citri subsp. citri, named X. citri A(T), which is deficient in xanthan production, induces an atypical, noncankerous chlorotic phenotype in Citrus limon and C. paradisi and weak cankerous lesions in C. aurantifolia and C. clementina leaves. In C. limon, suppression of canker development is concomitant with an oxidative burst; xanthan is not implicated in the phenotype induced by this interaction, suggesting that other bacterial factors would be involved in triggering the defense response.

Revisiting the Specificity of PCR Primers for Diagnostics of Xanthomonas citri pv. citri by Experimental and In Silico Analyses

Asiatic citrus canker disease, caused by Xanthomonas citri pv. citri, seriously impacts citrus production worldwide. Two pathogenic variants, A and A*/A^w, have been described within this pathovar. Two additional pathovars of X. citri with a limited geographic distribution and reduced pathogenicity, namely X. citri pvs. aurantifolii and bilvae, are also pathogenic to citrus and some rutaceous species. Rapid and reliable identification is required for these citrus pathogens, which are classified as a quarantine organism in citrus-producing countries. The specificity of nine polymerase chain reaction primers previously designed for the identification of X. citri pv. citri or citrus bacterial canker strains (both pvs. citri and aurantifolii) was assayed on a large strain collection (n = 87), including the two pathotypes of X. citri pv. citri, other genetic related or unrelated pathogenic xanthomonads, and saprophytic xanthomonads. This study gave congruent results with the original articles when testing the same strains or pathovars but the use of a broad inclusivity and exclusivity panel of strains highlighted new findings. Particularly, primers 2/3, 4/7, and KingF/R failed to provide amplification for three strains from the pathotype A*/A^w. Moreover, all pairs of primers detected at least one non-target strain. These data were supported by in silico analysis of the DNA sequences available from National Center for Biotechnology Information databases.

Development of variable number of tandem repeats typing schemes for Ralstonia solanacearum, the agent of bacterial wilt, banana Moko disease and potato brown rot

Ralstonia solanacearum is an important soil borne bacterial plant pathogen causing bacterial wilt on many important crops. To better monitor epidemics, efficient tools that can identify and discriminate populations are needed. In this study, we assessed variable number of tandem repeats (VNTR) genotyping as a new tool for epidemiological surveillance of R. solanacearum phylotypes, and more specifically for the monitoring of the monomorphic ecotypes "Moko" (banana-pathogenic) and "brown rot" (potato-pathogenic under cool conditions). Screening of six R. solanacearum genome sequences lead to select 36 VNTR loci that were preliminarily amplified on 24 strains. From this step, 26 single-locus primer pairs were multiplexed, and applied to a worldwide collection of 337 strains encompassing the whole phylogenetic diversity, with revelation on a capillary-electrophoresis genotype. Four loci were monomorphic within all phylotypes and were not retained; the other loci were highly polymorphic but displayed a clear phylotype-specificity. Phylotype-specific MLVA schemes were thus defined, based on 13 loci for phylotype I, 12 loci for phylotype II, 11 loci for phylotype III and 6 for phylotype IV. MLVA typing was significantly more discriminative than egl-based sequevar typing, particularly on monomorphic "brown rot" ecotype (phylotype IIB/sequevar 1) and "Moko disease" clade 4 (Phylotype IIB/sequevar 4). Our results raise promising prospects for studies of population genetic structures and epidemiological monitoring.

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

Pseudomonas syringae pv. maculicola causes bacterial spot on Brassicaceae worldwide, and for the last 10 years severe outbreaks have been reported in the Loire Valley, France. P. syringae pv. maculicola resembles P. syringae pv. tomato in that it is also pathogenic for tomato and causes the same types of symptoms. We used a collection of 106 strains of P. syringae to characterize the relationships between P. syringae pv. maculicola and related pathovars, paying special attention to P. syringae pv. tomato. Phylogenetic analysis of gyrB and rpoD gene sequences showed that P. syringae pv. maculicola, which causes diseases in Brassicaceae, forms six genetic lineages within genomospecies 3 of P. syringae strains as defined by L. Gardan et al. (Int. J. Syst. Bacteriol. 49[Pt 2]:469-478, 1999), whereas P. syringae pv. tomato forms two distinct genetic lineages. A multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) conducted with eight minisatellite loci confirmed the genetic structure obtained with rpoD and gyrB sequence analyses. These results provide promising tools for fine-scale epidemiological studies on diseases caused by P. syringae pv. maculicola and P. syringae pv. tomato. The two pathovars had distinct host ranges; only P. syringae pv. maculicola strains were pathogenic for Brassicaceae. A subpopulation of P. syringae pv. maculicola strains that are pathogenic for Pto-expressing tomato plants were shown to lack avrPto1 and avrPtoB or to contain a disrupted avrPtoB homolog. Taking phylogenetic and pathological features into account, our data suggest that the DC3000 strain belongs to P. syringae pv. maculicola. This study shows that P. syringae pv. maculicola and P. syringae pv. tomato appear multiclonal, as they did not diverge from a single common ancestral group within the ancestral P. syringae genomospecies 3, and suggests that pathovar specificity within P. syringae may be due to independent genetic events.

Genetic Diversity of a Brazilian Strain Collection of Xanthomonas citri subsp. citri Based on the Type III Effector Protein Genes

Exclusion and eradication or management based on an integrated approach with less susceptible varieties, copper-based bactericides, and windbreaks are the two main strategies used to prevent or control citrus canker. Field tolerance or resistance to citrus canker is not found in the most important commercial sweet orange cultivars, and pathogen-derived resistance has been developed and applied in different crops to obtain resistant genotypes to plant pathogens. We describe the development of DNA primers and probes based on the type III effector genes avrXacE1, avrXacE2, avrXacE3, avrBs2, pthA4, hpaF, and XAC3090 (leucine rich protein), and their application in the evaluation of the genetic diversity of the pathogen. A total of 49 haplotypes were identified in 157 strains by Southern blot analysis. No genetic polymorphism was detected by BOX elements - and enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR) analysis, nor with the genes avrBs2, XAC3090, and hpaF. Nei's genetic diversity indexes varied from 0.65 to 0.96 for subcollections of the pathogen. One or few haplotypes were most frequent in the strain collection, but several haplotypes were represented by solely one or few strains. The PthA4 probe resulted in the higher number of haplotypes identified in the Brazilian subcollections. Greater variation in the frequency of haplotypes occurred within subcollections (93.7%) than among subcollections. Only some haplotypes were genetically distant from all others, especially those originated from Rio Grande do Sul and Santa Catarina states. These bacterial effectors are widely spread in the collections and are useful for a better understanding of the host-pathogen interaction and the search for resistance genes in host and nonhost plants.

Controlling rice bacterial blight in Africa: needs and prospects

Rice cultivation has drastically increased in Africa over the last decade. During this time, the region has also seen a rise in the incidence of rice bacterial blight caused by the pathogen Xanthomonas oryzae pv. oryzae. The disease is expanding to new rice production areas and threatens food security in the region. Yield losses caused by X. oryzae pv. oryzae range from 20 to 30% and can be as high as 50% in some areas. Employing resistant cultivars is the most economical and effective way to control this disease. To facilitate development and strategic deployment of rice cultivars with resistance to bacterial blight, biotechnology tools and approaches, including marker-assisted breeding, gene combinations for disease control, and multiplex-PCR for pathogen diagnosis, have been developed. Although these technologies are routinely used elsewhere, their application in Africa remains limited, usually due to high cost and advanced technical skills required. To combat this problem, developers of the technologies at research institutions need to work with farmers from an early stage to create and promote the integration of successful, low cost applications of research biotech products. Here, we review the current knowledge and biotechnologies available to improve bacterial blight control. We will also discuss how to facilitate their application in Africa and delivery to the field.

Insertion sequence- and tandem repeat-based genotyping techniques for Xanthomonas citri pv. mangiferaeindicae

Molecular fingerprinting techniques that have the potential to identify or subtype bacteria at the strain level are needed for improving diagnosis and understanding of the epidemiology of pathogens such as Xanthomonas citri pv. mangiferaeindicae, which causes mango bacterial canker disease. We developed a ligation-mediated polymerase chain reaction targeting the IS1595 insertion sequence as a means to differentiate pv. mangiferaeindicae from the closely related pv. anacardii (responsible for cashew bacterial spot), which has the potential to infect mango but not to cause significant disease. This technique produced weakly polymorphic fingerprints composed of ≈70 amplified fragments per strain for a worldwide collection of X. citri pv. mangiferaeindicae but produced no or very weak amplification for pv. anacardii strains. Together, 12 tandem repeat markers were able to subtype X. citri pv. mangiferaeindicae at the strain level, distinguishing 231 haplotypes from a worldwide collection of 299 strains. Multilocus variable number of tandem repeats analysis (MLVA), IS1595-ligation-mediated polymerase chain reaction, and amplified fragment length polymorphism showed differences in discriminatory power and were congruent in describing the diversity of this strain collection, suggesting low levels of recombination. The potential of the MLVA scheme for molecular epidemiology studies of X. citri pv. mangiferaeindicae is discussed.

Genetic and Pathological Diversity Among Xanthomonas Strains Responsible for Bacterial Spot on Tomato and Pepper in the Southwest Indian Ocean Region

Bacterial spot of tomato and pepper, a major problem in tropical climates, can be caused by several Xanthomonas genospecies. We examined the genetic and pathological diversity of a collection of 72 strains from the southwest Indian Ocean region as part of a regional research and development program to update inventories of agricultural pests and pathogens. Xanthomonas euvesicatoria, X. perforans, X. gardneri, and X. vesicatoria were identified in our strain collection. The identification of strains at the species level was consistently achieved by amplified fragment length polymorphism (AFLP) and multilocus sequence analysis (MLSA). Overall, X. euvesicatoria was the species recovered prevalently. MLSA data based on four housekeeping genes identified two to three sequence types per genospecies. It suggested that sequence variations primarily consisted of synonymous mutations, although a recombination event spanning several hundred nucleotides was detected for some strains of X. euvesicatoria on the atpD gene coding for the F1-F0-ATPase β subunit. The pathogenicity of strains was consistent with data found in the literature. Some pathological variations were primarily observed among strains identified as X. euvesicatoria. This study provides the first ever comprehensive description of the status of Xanthomonas species that cause bacterial spot of tomato and pepper in the southwest Indian Ocean region.

Amplified fragment length polymorphism and multilocus sequence analysis-based genotypic relatedness among pathogenic variants of Xanthomonas citri pv. citri and Xanthomonas campestris pv. bilvae

Three pathogenic variants (i.e. pathotypes) have been described within Xanthomonas citri pv. citri, the causal agent of Asiatic citrus canker. Pathotype A strains naturally infect a wide range of Citrus species and members of some related genera. In contrast, pathotypes A* and Aw have narrow host ranges within the genus Citrus and have been isolated from Mexican lime (Citrus aurantifolia L.) and from Mexican lime and alemow (Citrus macrophylla L.), respectively. We used amplified fragment length polymorphism (AFLP) and multilocus sequence analysis (MLSA) based on four partial housekeeping gene sequences (atpD, dnaK, efp and gyrB) for the genotypic classification of Xanthomonas citri pv. citri and the poorly characterized citrus pathogen Xanthomonas campestris pv. bilvae. A Mantel test showed that genetic distances derived from AFLP and MLSA were highly correlated. X. campestris pv. bilvae showed a close relatedness to the type strain of X. citri, indicating that this pathovar should be reclassified as X. citri pv. bilvae. All pathotype A* and Aw strains were most closely related to X. citri pv. citri strains with a wide host range (pathotype A), confirming previous DNA–DNA hybridization data. Pathotype Aw should be considered a junior synonym of pathotype A* on the basis of pathogenicity tests, AFLP, MLSA and PCR using pathovar-specific primers. Evolutionary genome divergences computed from AFLP data suggested that pathotype A* (including Aw strains) is a group of strains that shows a wider genetic diversity than pathotype A.