Next Article in Journal
Re-Emergence of Salmonellosis in Hog Farms: Outbreak and Bacteriological Characterization
Next Article in Special Issue
Bacteriophage-Mediated Control of Phytopathogenic Xanthomonads: A Promising Green Solution for the Future
Previous Article in Journal
Favipiravir Does Not Inhibit Chikungunya Virus Replication in Mosquito Cells and Aedes aegypti Mosquitoes
Previous Article in Special Issue
Trends in Molecular Diagnosis and Diversity Studies for Phytosanitary Regulated Xanthomonas
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 5
10.3390/microorganisms9050945
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity and Geographical Structure of Xanthomonas citri pv. citri on Citrus in the South West Indian Ocean Region

by
Olivier Pruvost
1,*,
Damien Richard
1,2,3,
Karine Boyer
1,
Stéphanie Javegny
1,
Claudine Boyer
1,
Frédéric Chiroleu
1,
Pierre Grygiel
1,
Evelyne Parvedy
1,
Isabelle Robène
1,
Véronique Maillot-Lebon
1,
Azali Hamza
4,
Kanta Kumar Lobin
5,
Marc Naiken
6 and
Christian Vernière
1,7
1
CIRAD, UMR PVBMT, F-97410 Saint Pierre, La Réunion, France
2
ANSES, Plant Health Laboratory, F-97410 St Pierre, La Réunion, France
3
UFR Sciences et Technologies, Université de la Réunion, UMR PVBMT, F-97490 St Denis, La Réunion, France
4
INRAPE, Moroni BP 289, Comoros
5
FAREI, Le Réduit 80835, Mauritius
6
National Biosecurity Agency, Victoria P.O Box 464, Mahé, Seychelles
7
PHIM Plant Health Institute, CIRAD, INRAE, Institut Agro, IRD, University Montpellier, F-34398 Montpellier, France
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(5), 945; https://doi.org/10.3390/microorganisms9050945
Submission received: 24 March 2021 / Revised: 16 April 2021 / Accepted: 26 April 2021 / Published: 27 April 2021
(This article belongs to the Special Issue Integrating Science on Xanthomonas and Xylella for Integrated Plant Disease Management)
Browse Figures
Versions Notes

Abstract

:
A thorough knowledge of genotypic and phenotypic variations (e.g., virulence, resistance to antimicrobial compounds) in bacteria causing plant disease outbreaks is key for optimizing disease surveillance and management. Using a comprehensive strain collection, tandem repeat-based genotyping techniques and pathogenicity assays, we characterized the diversity of X. citri pv. citri from the South West Indian Ocean (SWIO) region. Most strains belonged to the prevalent lineage 1 pathotype A that has a wide host range among rutaceous species. We report the first occurrence of genetically unrelated, nonepidemic lineage 4 pathotype A* (strains with a host range restricted to Mexican lime and related species) in Mauritius, Moheli and Réunion. Microsatellite data revealed that strains from the Seychelles were diverse, grouped in three different clusters not detected in the Comoros and the Mascarenes. Pathogenicity data suggested a higher aggressiveness of strains of one of these clusters on citron (Citrus medica). With the noticeable exception of the Comoros, there was no sign of recent interisland movement of the pathogen. Consistent with this finding, the copL gene, a marker for the plasmid-borne copLAB copper resistance that was recently identified in Réunion, was not detected in 568 strains from any islands in the SWIO region apart from Réunion.

1. Introduction

Emerging and invasive plant pests and pathogens are a major threat to food security and the sustainability of agricultural cash crop industry [1,2]. Emerging diseases are diseases that occur in new areas (i.e., geographical expansion of preadapted pests and pathogens), on new hosts (i.e., host jumps or changes in host range) or that present new traits (e.g., acquired resistance to antimicrobials) [3]. Geographical expansion is a major factor in the emergence of plant pathogens [4]. The rapid increase in the geographical expansion of preadapted pests and pathogens over the last century is undoubtedly linked to the globalization of international trade and the increase of international travel of humans [5,6].
Collectively, citrus is the most important fruit crop in terms of world production. Several major bacterial, fungal and viral pathogens, which usually figure on quarantine or action lists, represent a major challenge to citrus industries. Preventing their establishment or delaying their spread is a priority. Xanthomonas citri pv. citri (synonym = X. citri subsp. citri) causes Asiatic citrus canker (ACC). This bacterium originated from Asia (i.e., the center of origin of Citrus) [7] and has spread globally due to the transport of citrus plant material over long distances [8]. X. citri pv. citri induces canker-like lesions on the external tissues of leaf, fruit and stem. Severe infections cause premature fruit drop, defoliation and twig dieback, resulting in yield loss and export market restrictions. Infection occurs through natural openings (i.e., stomata) or wounds. Stomatal infections only occur on actively growing aerial organs of citrus, whereas wounds are susceptible to infection over much larger time frames [9,10].
Strains of X. citri pv. citri differ in host range and defense reaction types on nonhost species, which determines the pathotype classification. The host range of pathotypes A* and Aw is limited to certain acid citrus types. They have primarily been sampled from natural infections on Mexican lime (Citrus x aurantiifolia) [11,12]. In contrast, pathotype A strains infect nearly all Citrus species, as well as other rutaceous genera [8]. Molecular diagnostics and genotyping techniques are essential for efficient disease management in agroecosystems [13]. Minisatellite and microsatellite genotyping, respectively, are well suited for global and small-scale studies on the molecular epidemiology of X. citri pv. citri [14,15,16,17]. In general, they are congruent with SNP-based phylogenies [15,18]. Pathotype A strains are split into two lineages (lineage 1 and 2), while pathotype Aw and A* were each assigned to a single lineage (lineage 3 and 4, respectively) [15]. Pathotypes Aw and A* globally have a more restricted geographical distribution and have less impact on agriculture because of their narrow host range, although severe outbreaks have been reported in countries where their host species are grown widely [19,20,21]. Lineage 1 strains were largely responsible for the major geographical expansion of X. citri pv. citri in the 20th century and for most worldwide outbreaks [15]. Recently reported, the agricultural significance of lineage 2 strains has not been widely documented, although they are present in the pathogen’s area of origin and are expanding in West Africa [15,16].
ACC control is primarily based on eradication or integrated management. The latter involves a combination of measures. These include the use of disease-free, partially-resistant nursery plants; efficient windbreaks; the use of drip irrigation instead of overhead irrigation; the extensive application of copper compounds (or use of antibiotics, which is less frequent) when plant organs are most susceptible to infection and pruning of diseased shoots during grove maintenance operations [8,22]. The impact of ACC is host-, environment- and bacterial lineage-dependent. Bacterial spread and infection are primarily driven by the use of X. citri pv. citri-contaminated nursery plants for grove establishment but also by extreme weather events [14,16]. Another factor promoting infection is resistance to antibacterial copper compounds, which reduces the effectiveness of ACC control [23,24]. To date, copper resistance is mostly restricted to lineage 1 strains found in three areas: South America (Argentina), the Caribbean region (Martinique) and the South West Indian Ocean (SWIO) region (Réunion). Copper-resistant X. citri pv. citri strains emerged in Réunion since 2010 [23,25] but the presence of the phenotype in other SWIO islands has not been formally assessed yet. ACC was reported in the SWIO region in the 20th century. ACC currently prevents the cultivation of susceptible citrus cultivars in this region, where the regular occurrence of hurricanes greatly exacerbates disease impact. Disease incidence on fruit of 30–50% (altering fruit quality and inducing early fruit drop) are commonly observed on susceptible cultivars. The earliest report was in Mauritius in 1917, and the latest was in Anjouan and Grande Comore in 2014 [26,27]. ACC was listed as present in Madagascar (EPPO global database, CABI disease map 011), but there is some doubt about disease occurrence, which has not been confirmed by laboratory analyses during grove inspections.
In the present study, we built on the high discriminatory power of microsatellite markers to assess the genetic structure of X. citri pv. citri in the SWIO region and to estimate the importance of strain circulation between the islands. A subset of strains representative of the genetic diversity was assigned to X. citri pv. citri lineages using minisatellite data. More specifically, we set out to answer the following questions: (i) Do all islands in the SWIO region host X. citri pv. citri genetic lineage 1, previously reported in Réunion, and to what extent are they related? (ii) Does this region host previously unreported lineages, and how prevalent are they? (iii) Is copper resistance present on SWIO islands other than Réunion? To answer these questions, we assembled and characterized a comprehensive strain collection, sampled from citrus in the Comoros (Anjouan, Grande Comore, Moheli and Mayotte), the Mascarenes (Mauritius, Rodrigues) and the Seychelles (Mahé). We compared them to the previously analysed X. citri pv. citri strains from Réunion [23]. Strains representative of the observed genetic diversity were further compared for pathogenicity to several Citrus species.

2. Materials and Methods

2.1. Bacterial Strains and Media

We analysed a total of 568 SWIO strains sampled from (i) the Comoros, i.e., Anjouan (n = 21), Grande Comore (n = 78), Moheli (n = 87) and Mayotte (n = 92); (ii) the Mascarenes, i.e., Mauritius (n = 135) and Rodrigues (n = 73); and (iii) the Seychelles, i.e., Mahe (n = 82) (Figure 1 and Table S1). While most strains of X. citri pv. citri were isolated in 2011–2017, the collection also included reference strains isolated between 1980 and 2010 (n = 57). These strains were compared to strains from Réunion that were analysed in a previous study [23]. Stock cultures (stored at −80°C) were initially produced from single colonies recovered after isolation from ACC symptoms, using KC semiselective medium. Subcultures of single colonies were prepared on YPGA plates [28]. Bacterial suspensions were produced in 0.01 M sterile Sigma 7–9 Tris buffer pH 7.2 (Sigma–Aldrich, Saint-Quentin Fallavier, France), unless otherwise stated. We determined copper phenotype (susceptible vs. resistant) using YPGA plates, with or without copper sulphate pentahydrate (470 mg L−1). Strain LH201 was used as the copper-resistant control [23].

2.2. Microsatellite Genotyping (MLVA-14)

The MLVA-14 scheme targeting 14 microsatellites (6–7 bp-long) in a multiplex PCR format was used to decipher the genetic structure of all SWIO X. citri pv. citri strains. We proceeded as reported earlier [17], with one exception: in the multiplex PCR assay (PCR mix pool #1), we included a primer pair targeting the copL gene [29], required for copper resistance of X. citri pv. citri in Réunion. Briefly, one of each primer in each PCR mix was 5′-labelled with one of the following fluorescent dyes: 6-FAM, NED, PET and VIC (Applied Biosystems, Courtaboeuf, France). As a control in each experiment, we used the X. citri pv. citri strain LN002-01, a rifamycin-resistant derivative of IAPAR 306 [30] in which the copper resistance plasmid pLH201.1 [23] was introduced by conjugation.
The number of tandem repeats at each microsatellite locus was computed from fragment length for all assayed strains using GeneMapper 4.0 (Applied Biosystems, Courtaboeuf, France) and used as input data. Nei’s unbiased estimates of genetic diversity for MLVA-14 data were calculated using the poppr 2.8.3 package [31] in R version 3.6.1 (R Core Team, 2019. A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/ (accessed on 27 April 2021)). Allelic richness (A) was calculated using the rarefaction procedure for unequal sample sizes (subsample size n = 21) using the hierfstat 0.04–22 R package [32]. Private alleles were identified using poppr. Genetic clusters (GCs) were identified as networks of haplotypes that differed by more than five microsatellite loci, using the algorithm recommended for MLVA data, combining global optimal eBURST and Euclidean distances in PHYLOViZ v1.2 [33]. Relaxed clonal complexes (RCCs) were delineated as networks of MLVA-14 haplotypes linked by up to triple-locus variations. The population structure of X. citri pv. citri was further assessed using discriminant analysis of principal components (DAPC) [34]. Here, we used the adegenet V.2.1.1. R package [35], which involved optimizing the retained number of principal components from cross-validation analysis. DAPC was used, as it is free of any assumption linked to a population genetic model (e.g., Hardy–Weinberg equilibrium or linkage equilibrium), and, thus, it is suited for analysis of datasets produced from predominantly clonal bacteria.

2.3. Minisatellite Typing (MLVA-31)

The MLVA-31 scheme was used to analyse the genetic relatedness (assessed from the genotyping of 31 minisatellites) between a subcollection of SWIO strains selected on the basis of MLVA-14 data (Table S1) and a worldwide strain collection of X. citri pv. citri [15,16,21] (http://bioinfo-web.mpl.ird.fr/MLVA_bank/Genotyping/view.php (accessed on 27 April 2021)). Primer pairs targeting single-locus alleles were used in a multiplex PCR format (Clontech Terra PCR Direct Polymerase Mix). One of each primer in the PCR mix was 5′-labelled with one of the following fluorescent dyes: 6-FAM, NED, PET and VIC (Applied Biosystems). PCRs and capillary electrophoresis were performed with an initial injection of 23 s, as described earlier [17]. X. citri pv. citri strain IAPAR 306 [30] was used as a control in each experiment.
Input data were obtained as described above (Section 2.2). The MLVA-31 dataset generated in this study was deposited in the Xanthomonas citri genotyping database (http://bioinfo-web.mpl.ird.fr/MLVA_bank/Genotyping/view.php (accessed on 27 April 2021)) [15]. PHYLOViZ v1.2 was used to build a minimum-spanning tree (MST) from the MLVA-31 dataset [33]. We assigned SWIO strains to genetic lineages using discriminant analysis of principal components (DAPC), as previously reported [16,34]. Subclusters were identified as groups of strains linked by up to triple-locus variations using PHYLOViZ [15].

2.4. Typeability of SWIO Strains with XAC1051-qPCR

A subset of 42 strains (Table S1), including strains assigned to all lineage 1 GCs, as well as the four pathotype A* strains, were assayed using the X. citri pv. citri-specific real-time quantitative assay XAC1051-qPCR [36]. Amplifications were performed with the Quantstudio 5 (QS5) (Applied Biosystems) real-time PCR system, using the GoTaq® probe qPCR master mix kit (Promega, Charbonnières-les-Bains, France). Amplifications were carried out on 2 µL bacterial suspensions containing approximately 1 × 108 cells mL−1 (prepared as described above) using only the primers/MGB Taqman® probe specific to X. citri pv. citri (the amplification of the plant internal control was not relevant here).

2.5. Detached Leaf Assay-Based Pathogenicity Tests

Pathogenicity tests were conducted on eighty-six SWIO strains (all strains assigned to pathotype A* and a subcollection of lineage 1 strains, including all geographical origins selected on the basis of MLVA-14 data) (Table S1). Inoculations were performed on at least two Citrus species, namely, C. × aurantifolia (Mexican lime) plus one or more of the following: C. x sinensis (Washington and/or New Hall navel sweet orange), C. × sinensis × C. reticulata (Ortanique tangor) and C. x paradisi (Henderson, Marsh and/or Star Ruby grapefruit). We followed a method described previously [11]; the only difference here was that we used 5 µl droplets of bacterial suspensions containing approximately 1 x 108 cells mL−1 (prepared as described above) for inoculation. Sterile buffer and strain IAPAR 306 were used as negative and positive controls, respectively.

2.6. Attached Leaf Assay-Based Pathogenicity and in Planta Growth of Lineage 1 Strains

Assays were conducted on three different citrus lines: New Hall navel sweet orange (C. × sinensis), Zanzibar mandarin (C. reticulata) and Buddha’s hand citron (C. medica). We used spray inoculation, because it reveals subtle differences in virulence, which are not shown by infiltration inoculation [37]. We conducted assays on sixteen X. citri pv. citri strains, four of each main GC (Table S2). Six- to eight-month-old citrus plants (three plants per strain–host combination) were used for inoculations. Plants were placed in growth chambers at 28 ± 1 °C day and 26 ± 1 °C night and 80 ± 5% relative humidity with a photoperiod of 12 h. Newly expanded young flushes were spray inoculated until run-off with bacterial suspensions containing approximately 1 × 108 cells mL−1, prepared as described above. Plants were then covered with large clear polyethylene bags overnight. All plants were incubated for 25 days. All inoculated leaves were examined (daily during the first week after inoculation and then every other day) for lesion development. The following scale was used: 0 = no lesion; 1 = one to 10 lesions; 2 = 11 to 50 lesions; 3 ≥ 50 lesions. The susceptibility of leaves (from young citrus flushes) to stomatal infections by X. citri pv. citri is clearly age-dependent. Although stomata develop during the very early stages of leaf development, stomatal pores are not open at growth stages < 50% expansion, as shown by scanning electron microscopy [9,38]. When leaves are nearing their mature size, the fully developed cuticle acts as a protective layer against stomatal infections [8]. Therefore, we expected that, after spray inoculation, all the leaves along an actively growing flush would present different susceptibility to X. citri pv. citri. We expected to see partial resistance on older leaves. When we spray-inoculated very young leaves, we expected a delayed lesion development until the stomata became susceptible. Areas under disease progression stairs (AUDPS) were computed from disease severity data using the agricolae 1.3.3 R package [39]. AUDPS were converted into a ratio by dividing by the maximal area recorded and corrected to within ]0; 1[. A mixed generalized linear model with beta distribution was built for each host species with the strain or GC as the fixed effect and the plant as random effect using glmmTMB 1.0.2.1 R package [40]. The strain effect was analysed by deviance test using car 3.0.6 R package [41], and multiple comparisons, using emmeans 1.5.3 R package. In order to limit the heterogeneity in leaf susceptibility, the amount of censored data was estimated by maximizing the log-likelihood of the beta regression model for each host species. Binary disease status datasets were also produced from disease severity scores and were analysed using survival analysis [42] with the survival 3.1.8 [43] and survminer 0.4.6 R packages. The analysis estimates the probability for a leaf to remain healthy until a given point in time. Comparisons between strains or GCs were performed using chi-square-based test statistics [42].
Twenty-five days after inoculation (DAI), we enumerated X. citri pv. citri population sizes from single canker lesions. For each strain–cultivar combination, nine single canker lesions were collected using a sterile 8 mm punch. Each lesion was placed in a sterile microtube containing 1 mL of sterile buffer and gently agitated for 15 min. (50 rpm) at room temperature, mimicking the natural release of the bacterium from canker lesions in the presence of free water. Each lesion was then transferred to a fresh tube and homogenized for 30 s in 1 mL of sterile buffer with sterile ceramic beads in a FastPrep24 device (MP Biomedicals). Appropriate dilutions of washates and grindates were plated on semiselective KC agar plates using a Spiral plater (Interscience). Plates were incubated for 3–4 days at 28 °C, and Xanthomonas-like colonies were enumerated, as recommended by the manufacturer. Exuded X. citri pv. citri population sizes, expressed as log(ufc.lesion−1) were calculated from washate enumerations. Total X. citri pv. citri population sizes were calculated from washate + grindate enumerations. For each plant species, the strain effect on bacterial population size was analysed by ANOVA with Box–Cox power data transformation (residual variance homogenisation) using car and MASS 7.3–51.4 packages [44], followed by Tukey’s multiple pairwise comparisons using the multcomp 1.4–15 R package [45].

3. Results

3.1. Genetic Diversity among Strains from Different Islands in the SWIO Region Revealed by Microsatellite Typing and Cop-PCR

We used microsatellite typing to resolve the genetic structure among X. citri pv. citri strains at small evolutionary scales [14,16,17,21,23]. We identified 390 haplotypes among the 568 assayed strains. At the island scale, allelic richness (A) ranged from 3.1 to 5.6 (Table 1). The islands that have hosted the pathogen for long periods tended to display higher A values. The highest value of private allelic richness AP (i.e., alleles solely present in strains from a single island) was found in the Seychelles population. It was also high in Réunion. For other islands, AP was low to null (Anjouan and Grande Comore) (Table 1). The multiplex PCRs did not reveal copL PCR amplification from any assayed strain, although amplification was obtained for the copper-resistant control strain from Réunion. As expected, none of the 47 subsampled strains (originating from all sampled islands) grew on media supplemented with copper, confirming their copper-susceptible phenotype.
The genetic relatedness among SWIO strains was analysed from the microsatellite allelic profiles from our strain collection together with the ones previously obtained for Réunion strains [23]. The minimum spanning tree (not shown) suggested four main genetic clusters (GC) in the dataset. This structure was further supported by DAPC (Figure 2). Ninety-three percent of haplotypes identified in the Comoros (Anjouan, Grande Comore, Moheli and Mayotte) and the Mascarenes (Mauritius, Réunion and Rodrigues) clustered as GC1. The remaining haplotypes from these locations were structured as minor clusters (i.e., networks grouping ≤3 haplotypes) or singletons (i.e., haplotypes whose multilocus allelic profile differed from its closest relative at ≥6 microsatellite loci) (Table 2). Interestingly, all strains originating from the Seychelles were grouped into three distinct GCs, GC2 (n = 44), GC3 (n = 25) and GC4 (n = 13), which did not contain any haplotypes from any other island.
The structure of GC1 from the Comoros and Mascarene Archipelagos was further examined by delineating RCCs, which grouped together closely related strains, in order to identify putative epidemiological links between strains. A total of 75 RCCs, nine of which grouped ≥20 haplotypes, and 141 singletons were delineated. In the Mascarenes, very few strains originating from distinct islands grouped in the same RCC. Only two pairs of strains originating from Mauritius and Réunion (all isolated in the 1980s) shared identical haplotypes. In Réunion, RCC1, RCC2 and RCC8 included haplotypes that were distributed all over the island. In contrast, the distribution of RCC4, RCC5 and RCC7 was geographically restricted to one or two localities no more than 10 km apart. In Mauritius, most of the strains grouped in RCC6, which was distributed all over the island. In contrast, the structure of populations in Rodrigues generally had a large number of singletons and small RCCs, often restricted to a single or spatially close citrus blocks. The situation was markedly different in the Comoros, where most haplotypes originating from Grande Comore (n = 21), Anjouan (n = 10) and Mayotte (n = 46) clustered in RCC3, revealing a close genetic relatedness among strains from the three islands (Figure 3). Strains from Moheli appeared only distantly related to strains from other Comoros islands.
All 42 strains, representative of the genetic diversity among strains from different islands in the SWIO region, tested positive with the X. citri pv. citri-specific XAC1051-qPCR assay. The cycle threshold (Ct) values ranged from 10.8 to 17.5.

3.2. Minisatellite-Based Assignation of SWIO Strains to X. citri pv. citri Genetic Lineages

A subset of 215 strains from all GCs, minor clusters and singletons was submitted to minisatellite genotyping in order to assign them to the genetic lineages identified in X. citri pv. citri (Table 2) [15]. This complemented a dataset previously obtained from Réunion [23], yielding a total of 388 MLVA-31 genotyped SWIO strains.
Based on DAPC, most strains from the SWIO region were assigned to genetic lineage 1 (pathotype A). A network of 73 MLVA haplotypes differing by single- or double-locus variations was produced among these strains (Figure S1). Sixteen of these haplotypes were identified from at least two distinct SWIO islands. Haplotype #2, which was detected in five islands (Anjouan, Mayotte, Moheli, Mauritius and Réunion) and previously identified in different regions (Asia, Middle East, Oceania and South America) [16], was identified as the most frequent MLVA-31 haplotype in the SWIO region and as the putative founder haplotype. Although 13 minisatellite loci were shown to be polymorphic among SWIO lineage 1 strains, most (50/72) polymorphisms detected over the minimum spanning tree path occurred at three fast-evolving loci (Xcc3816, Xcc3993 and Xcc4748). Polymorphism at these three loci was previously identified among epidemiologically related strains [15]. An allelic state for locus Xcc2059 was found specific to GC2 strains (Seychelles), and the corresponding haplotypes (#168, 210 and 217) have not been identified from other countries to date [16].
Four strains isolated from Mauritius, Moheli and Réunion were assigned to genetic lineage 4, so far solely comprising pathotype A* strains (Table 2). These strains were genetically diverse (i.e., their respective allelic states differed at ≥8 minisatellite loci). The two Mauritian strains JJ238-41 and LM053-15 were single-locus variants of strains NCPPB 3615 (India) and NIGEB196-1 (Iran), respectively. Strain LK135-01 from Moheli shared the same allelic profile as the Thai strain JJ238-24, although the two strains were distinguishable using microsatellite data. The Réunion strain JN658 shared no close relatedness with other previously genotyped pathotype A* nor to any pathotype Aw strain.

3.3. Pathogenicity Assays Confirm the Prevalence of Pathotype A and the Limited Presence of Pathotype A* in the SWIO Region

A subset of strains from all GCs, minor clusters and singletons were submitted to pathogenicity tests on detached Mexican lime, as well as sweet orange, grapefruit and/or tangor leaves. Consistent with previous data, all SWIO strains that were assigned to genetic lineage 1 by minisatellite typing produced typical callus-like reactions on all assayed citrus species, similar to the positive control strain IAPAR 306. In contrast, all four strains assigned to genetic lineage 4 only produced callus-like reactions on Mexican lime, and no typical canker lesions were observed on the other citrus species assayed. This corresponds to the phenotype of the disease development index DDi III, reported earlier [46]. No lesion developed from negative controls, regardless of the citrus species assayed.

3.4. Relative Virulence of SWIO Strains from the Four Major Genetic Clusters on Three Citrus Species

We further assayed pathogenicity and in planta growth of selected strains from the four major GCs on attached leaves of four Citrus species. As expected, all 16 strains of X. citri pv. citri pathotype A lineage 1, whatever their GC assignment, started producing small canker-like lesions at 6 to 7 DAI on each assayed plant species. Based on survival analyses, the disease incidence progress of GC3 strains was greater than for the other GCs on citron and sweet orange (p < 0.0001). In contrast, survival analyses suggested a lack of GC-dependent response (p = 0.720) to inoculations on mandarin (Figures S2–S4). Differences among strains (p-values) are shown in Table S3.
Disease development patterns were investigated in more detail using AUDPS data. It further confirmed that GC3 was more aggressive on citron than GC1 (p = 0.004) and GC2 (p = 0.045), but not GC4 (p = 0.664). No significant differences among GCs were shown on mandarin and sweet orange (p > 0.38). The differences among strains (p-values) are shown in Table S4.
Large X. citri pv. citri population sizes were present in canker lesions, irrespective of the assayed bacterial strain and the citrus species at 25 DAI. Mean population sizes ranged from 4 × 106 to 3 × 107, 3 × 106 to 4 107 and 8 × 106 to 4 × 107 cfu lesion−1 on citron, mandarin and sweet orange, respectively (Table S5). When lesions were soaked in buffer, exudation of all strains was fast and massive for all assayed host species, consistent with earlier results [47,48]. Both exuded and total population sizes showed a significant strain * host interaction (ANOVA, p < 0.0001). Analyses were conducted again using split datasets. Differences among strains in total or exuded population sizes were ≤ 1 log unit and did not appear to be GC-related (Table 3).

4. Discussion

We conducted a comprehensive molecular and pathological characterization of 568 X. citri pv. citri strains from the SWIO region. The present study confirmed a high prevalence of ACC in the SWIO region (the Comoros, the Mascarenes and the Seychelles), where X. citri pv. citri causes severe repeated disease outbreaks in citrus groves. X. citri pv. citri populations in the SWIO region were genetically diverse. Importantly, no copper resistance was detected in islands other than Réunion.

4.1. Genetic Lineage 1 is Markedly Prevalent in the SWIO Region, Displays a Geographic Structure and Sporadically Cohabits with Nonepidemic Lineage 4 Strains

Consistent with its worldwide prevalence [15], the present study highlighted the marked predominance of pathotype A genetic lineage 1 strains in the SWIO region. Very recently, SNP data showed that genetic lineage 1 strains causing ACC in the Seychelles are genetically diverse and clearly differ from strains originating from the Comoros and the Mascarenes [49]. Although microsatellite data have recognized limitations for accurately assessing deep genetic relatedness among groups of individuals [50,51], the present study revealed large genetic differences between strains from the Seychelles vs. the Comoros and the Mascarenes. These results suggest that several distinct introduction events of X. citri pv. citri have occurred at both regional and island scale (e.g., Mahé in the Seychelles where strains isolated the same year from the same citrus species can be assigned to different GCs).
Here, we provide the first report of the presence of genetic lineage 4 pathotype A* strains in three islands in the SWIO region, namely Moheli, Mauritius and Réunion. These strains were found to be genetically diverse based on minisatellite and microsatellite data. Therefore, they are likely to represent distinct introductions in the SWIO region. The Moheli strain was distinct from but related to a Thai pathotype A* strain. The two Mauritian strains were distinct from but related to pathotype A* strains from India and Iran, respectively [21]. The pathotype A* strain from Réunion shared no close relatedness to any other pathotype A* genotyped to date. The citrus block where this strain was detected has been removed, and no other pathotype A* strains were detected despite extensive sampling [23]. Similarly, additional bacterial isolation attempts from the lime blocks in Mauritius and Moheli, where pathotype A* strains originated, confirmed the marked prevalence of pathotype A strains (data not shown). Epidemicity of pathotype A* was primarily reported from countries where lime cultivation is widespread (e.g., Iran, Ethiopia, Sudan). In the context of SWIO islands, these strains were of no agricultural significance. This is probably due to lime’s limited presence in the highly heterogeneous agricultural landscapes, which may prevent their epidemicity.

4.2. The Seychelles Host Genetically Diverse X. citri pv. citri Pathotype A Strains Differing in Pathogenicity

In the present study, we used a spray-inoculation technique to assess the pathogenicity of strains representative of the four GCs, which were identified in the SWIO region on three citrus species. A recent study evidenced that the structural and chemical properties of leaf surface components impact the susceptibility to X. citri pv. citri, i.e., pathogenicity differs among citrus lines. This was observed after spray inoculation but not after infiltration inoculation [52]. Interestingly, strains assigned to one of the unique groups detected in the Seychelles (GC3) differed to strains assigned to other GCs in terms of the dynamics of canker lesion development, i.e., disease incidence and AUDPS data on citron, an ancestral Citrus species from which some modern lines originated [7]. However, there was no marked difference in X. citri pv. citri population sizes, which were measured in canker lesions at 25 DAI, for all assayed GCs. This suggests that, compared to other GCs, GC3 strains can produce more canker lesions, which develop faster after spray inoculation, but do not survive in or exude from lesions at larger population sizes. Thus, in addition to a broad genetic polymorphism, the strains of X. citri pv. citri from the Seychelles showed greater polymorphism in terms of their pathogenicity. Given the Seychelles’ small size and relative isolation, the island of Mahé could be an interesting site for assessing gene flow patterns among admixed, genetically remote strains. This has been evidenced recently for Xanthomonas perforans using genomic and metagenomic approaches [53,54].

4.3. Microsatellite Typing Suggests that the Pathogen’s Inter- and Intra-Island Movements Differ

We found no sign of close genetic relatedness among strains originating from distinct archipelagos (i.e., no RCC grouped strains originating from distinct archipelagos). This suggests that the pathogen’s recent interisland movement through contaminated citrus material has been very limited. Consistent with this finding, no copper resistance associated with the copLAB system was identified in any strain originating from the SWIO islands, apart from Réunion, where it has been previously reported [23,25]. Indeed, besides a few exceptions, RCCs that were delineated among GC1 strains seemed largely island-specific.
In the Mascarenes, there was little evidence of high genetic relatedness among strains from distinct islands. The situation in the Comoros was more contrasted. Most haplotypes originating from Grande Comore, Anjouan and Mayotte displayed close genetic relatedness, suggesting that interisland dispersal of X. citri pv. citri, most likely through contaminated citrus propagative material, occurs in the Comoros. Conversely, none of the strains from Moheli appeared to be closely related to strains originating from other neighbour islands in the Comoros or, more broadly, from the SWIO region.
The intra-island distribution of haplotypes also revealed contrasting situations. Populations in Rodrigues and Moheli were structured primarily as a large number of singletons and small RCCs. The latter are frequently restricted to a single or spatially close citrus block. This structure might originate from the very low development of citrus industries. Indeed, local citrus growers often produce the plants they need for grove establishment. In addition, many strains from Rodrigues and Moheli were isolated from Mexican lime, which is generally grown for local consumption. The Mexican lime species is commonly propagated locally from seeds (which do not transmit the pathogen), thereby limiting exogenous contamination. In contrast, most of the strains from Mauritius or Réunion grouped in RCCs that were distributed all over the island, suggesting the diffusion of contaminated propagative citrus plant material from commercial nurseries. In Réunion, strains from these RCCs were also detected from commercial citrus nurseries. The trade of infected citrus nursery plants was repeatedly identified as a major pathway for X. citri pv. citri long-distance spread [14,16]. Overall, the present study emphasizes the importance of improving the control of ACC and other citrus diseases in the nurseries in the SWIO region. New legislation targeting ACC, Huanglongbing, tristeza (and their respective vectors) is currently being implemented in Réunion to improve the sanitary quality of citrus nursery plants. The XAC1051-qPCR assay, whose ability to detect all X. citri pv. citri variants identified in the SWIO region was herein confirmed, represents a valuable tool for ACC surveillance and management in the region.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9050945/s1, Figure S1: Minimum spanning tree representing the genetic relatedness based on minisatellite genotyping (MLVA-31) among lineage 1 strains of Xanthomonas citri pv. citri from the SWIO region. Dot diameter is representative of the number of strains per haplotype. Dot colour is representative of the archipelago/island of isolation (Comoros Archipelago: light green for Anjouan; medium green for Moheli; dark green for Grande Comore; khaki for Mayotte. Mascarene Archipelago: light blue for Rodrigues; medium blue for Mauritius; dark blue for Réunion. Seychelles Archipelago: red for Mahé). Single and double locus variations are represented as thick and thin solid lines joining haplotypes, respectively, Figure S2: Survival probability of citron (Citrus medica) leaf infection by Xanthomonas citri pv. citri over time, Figure S3: Survival probability of sweet orange (Citrus x sinensis) leaf infection by Xanthomonas citri pv. citri over time, Figure S4: Survival probability of mandarin (Citrus reticulata) leaf infection by Xanthomonas citri pv. citri over time, Table S1: Xanthomonas citri pv. citri strains isolated from citrus in the South West Indian Ocean (SWIO) region, Table S2: Xanthomonas citri pv. citri lineage 1 strains used for spray-inoculation experiments of three citrus species, Table S3: Pairwise strain comparisons (p-values) based on disease incidence data on citron, mandarin and sweet orange, Table S4 Pairwise strain comparisons (p-values) based on AUDPS data on citron, mandarin and sweet orange, Table S5: Comparisons among citrus species of exuded and total Xanthomonas citri pv. citri population sizes recorded from single canker lesions 25 days after spray-inoculation of 16 X. citri pv. citri strains.

Author Contributions

Conceptualization, O.P. and C.V.; Methodology, O.P., D.R., K.B., S.J., C.B., F.C., I.R., V.M.-L. and C.V.; Software: F.C.; Validation, O.P., D.R., K.B., F.C., I.R. and C.V.; Formal Analysis, O.P., D.R., K.B., F.C., I.R. and C.V.; Investigation: O.P., K.B., S.J., C.B., P.G., E.P., I.R., V.M.-L. and C.V.; Resources, O.P., A.H., K.K.L., M.N. and C.V.; Data Curation, O.P. and C.V.; Writing—Original Draft Preparation, O.P.; Writing—Review and Editing, D.R., K.B., S.J., C.B., F.C., P.G., E.P., I.R., K.K.L. and C.V.; Visualization, O.P., D.R. and C.V.; Supervision, O.P.; Project Administration, O.P.; Funding Acquisition, O.P., A.H. and M.N. All authors have read and agreed to the published version of the manuscript.

Funding

The European Union (ERDF contract GURDT I2016-1731-0006632, ERDF INTERREG IV and INTERREG V), Conseil Régional de La Réunion, the French Agropolis Foundation (Labex Agro—Montpellier, E-SPACE project number 1504-004), ANSES and CIRAD provided financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data produced in this study are presented in this publication. Minisatellite data are available in the Xanthomonas citri genotyping database (http://bioinfo-web.mpl.ird.fr/MLVA_bank/Genotyping/view.php (accessed on 27 April 2021)).

Acknowledgments

We would like to express our thanks to the plant protection platform (3P, IBISA) and to V. Ledoux and M. Panchbaya for their helpful contribution.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Pimentel, D.; McNair, S.; Janecka, J.; Wightman, J.; Simmonds, C.; O’Connell, C.; Wong, E.; Russel, L.; Zern, J.; Aquino, T.; et al. Economic and environmental threats of alien plant, animal, and microbe invasions. Agric. Ecosyst. Environ. 2001, 84, 1–20. [Google Scholar] [CrossRef]
  2. Vurro, M.; Bonciani, B.; Vannacci, G. Emerging infectious diseases of crop plants in developing countries: Impact on agriculture and socio-economic consequences. Food Secur. 2010, 2, 113–132. [Google Scholar] [CrossRef]
  3. Engering, A.; Hogerwerf, L.; Slingenbergh, J. Pathogen-host-environment interplay and disease emergence. Emerg. Microbes Infect. 2013, 2, e5. [Google Scholar] [CrossRef] [PubMed]
  4. Anderson, P.K.; Cunningham, A.A.; Patel, N.G.; Morales, F.J.; Epstein, P.R.; Daszak, P. Emerging infectious diseases of plants: Pathogen pollution, climate change and agrotechnology drivers. Trends Ecol. Evol. 2004, 19, 535–544. [Google Scholar] [CrossRef] [PubMed]
  5. Hulme, P.E. Trade, transport and trouble: Managing invasive species pathways in an era of globalization. J. Appl. Ecol. 2009, 46, 10–18. [Google Scholar] [CrossRef]
  6. Bebber, D.P.; Holmes, T.; Gurr, S.J. The global spread of crop pests and pathogens. Glob. Ecol. Biogeogr. 2014, 23, 1398–1407. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, G.A.; Terol, J.; Ibanez, V.; Lopez-Garcia, A.; Perez-Roman, E.; Borreda, C.; Domingo, C.; Tadeo, F.R.; Carbonell-Caballero, J.; Alonso, R.; et al. Genomics of the origin and evolution of Citrus. Nature 2018, 554, 311–316. [Google Scholar] [CrossRef] [Green Version]
  8. Graham, J.H.; Gottwald, T.R.; Cubero, J.; Achor, D.S. Xanthomonas axonopodis pvcitri: Factors affecting successful eradication of citrus canker. Mol. Plant Pathol. 2004, 5, 1–15. [Google Scholar] [CrossRef]
  9. Gottwald, T.R.; Graham, J.H. A device for precise and nondisruptive stomatal inoculation of leaf tissue with bacterial pathogens. Phytopathology 1992, 82, 930–935. [Google Scholar] [CrossRef]
  10. Vernière, C.; Gottwald, T.R.; Pruvost, O. Disease development and symptom expression of Xanthomonas axonopodis pv. citri in various citrus plant tissues. Phytopathology 2003, 93, 832–843. [Google Scholar] [CrossRef] [Green Version]
  11. Vernière, C.; Hartung, J.S.; Pruvost, O.P.; Civerolo, E.L.; Alvarez, A.M.; Maestri, P.; Luisetti, J. Characterization of phenotypically distinct strains of Xanthomonas axonopodis pv. citri from Southwest Asia. Eur. J. Plant Pathol. 1998, 104, 477–487. [Google Scholar] [CrossRef]
  12. Sun, X.A.; Stall, R.E.; Jones, J.B.; Cubero, J.; Gottwald, T.R.; Graham, J.H.; Dixon, W.N.; Schubert, T.S.; Chaloux, P.H.; Stromberg, V.K.; et al. Detection and characterization of a new strain of citrus canker bacteria from key Mexican lime and Alemow in South Florida. Plant Dis. 2004, 88, 1179–1188. [Google Scholar] [CrossRef] [Green Version]
  13. Bull, C.T.; Koike, S.T. Practical benefits of knowing the enemy: Modern molecular tools for diagnosing the etiology of bacterial diseases and understanding the taxonomy and diversity of plant-pathogenic bacteria. Annu. Rev. Phytopathol. 2015, 53, 157–180. [Google Scholar] [CrossRef]
  14. Vernière, C.; Bui Thi Ngoc, L.; Jarne, P.; Ravigné, V.; Guérin, F.; Gagnevin, L.; Le Mai, N.; Chau, N.M.; Pruvost, O. Highly polymorphic markers reveal the establishment of an invasive lineage of the citrus bacterial pathogen Xanthomonas citri pv. citri in its area of origin. Environ. Microbiol. 2014, 16, 2226–2237. [Google Scholar] [CrossRef]
  15. Pruvost, O.; Magne, M.; Boyer, K.; Leduc, A.; Tourterel, C.; Drevet, C.; Ravigné, V.; Gagnevin, L.; Guérin, F.; Chiroleu, F.; 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, e98129. [Google Scholar] [CrossRef]
  16. Leduc, A.; Traoré, Y.N.; Boyer, K.; Magne, M.; Grygiel, P.; Juhasz, C.; Boyer, C.; Guérin, F.; Wonni, I.; Ouedraogo, L.; 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. [Google Scholar] [CrossRef]
  17. Pruvost, O.; Boyer, K.; Ravigné, V.; Richard, D.; Vernière, C. 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. Evol. Appl. 2019, 12, 1523–1538. [Google Scholar] [CrossRef] [Green Version]
  18. Gordon, J.L.; Lefeuvre, P.; Escalon, A.; Barbe, V.; Cruveiller, S.; Gagnevin, L.; Pruvost, O. Comparative genomics of 43 strains of Xanthomonas citri pv. citri reveals the evolutionary events giving rise to pathotypes with different host ranges. BMC Genom. 2015, 16, 1098. [Google Scholar] [CrossRef] [Green Version]
  19. Bui Thi Ngoc, L.; Vernière, C.; Pruvost, O.; Kositcharoenkul, N.; Phawichit, S. First report in Thailand of Xanthomonas axonopodis pv. citri-A* causing citrus canker on lime. Plant Dis. 2007, 91, 771. [Google Scholar] [CrossRef]
  20. Derso, E.; Vernière, C.; Pruvost, O. First report of Xanthomonas citri pv. citri-A* causing citrus canker on lime in Ethiopia. Plant Dis. 2009, 93, 203. [Google Scholar] [CrossRef]
  21. Pruvost, O.; Goodarzi, T.; Boyer, K.; Soltaninejad, H.; Escalon, A.; Alavi, S.M.; Javegny, S.; Boyer, C.; Cottyn, B.; Gagnevin, L.; et al. Genetic structure analysis of strains causing citrus canker in Iran reveals the presence of two different lineages of Xanthomonas citri pv. citri pathotype A*. Plant Pathol. 2015, 64, 776–784. [Google Scholar] [CrossRef]
  22. Gottwald, T.R.; Graham, J.H.; Schubert, T.S. Citrus canker: The pathogen and its impact. Plant Health Prog. 2002, 3, 15. [Google Scholar] [CrossRef] [Green Version]
  23. Richard, D.; Ravigné, V.; Rieux, A.; Facon, B.; Boyer, C.; Boyer, K.; Grygiel, P.; Javegny, S.; Terville, M.; Canteros, B.I.; et al. Adaptation of genetically monomorphic bacteria: Evolution of copper resistance through multiple horizontal gene transfers of complex and versatile mobile genetic elements. Mol. Ecol. 2017, 26, 2131–2149. [Google Scholar] [CrossRef] [PubMed]
  24. Behlau, F.; Gochez, A.M.; Jones, J.B. Diversity and copper resistance of Xanthomonas affecting citrus. Trop. Plant Pathol. 2020, 45, 200–212. [Google Scholar] [CrossRef]
  25. Richard, D.; Tribot, N.; Boyer, C.; Terville, M.; Boyer, K.; Javegny, S.; Roux-Cuvelier, M.; Pruvost, O.; Moreau, A.; Chabirand, A.; et al. First report of copper-resistant Xanthomonas citri pv. citri pathotype A causing Asiatic citrus canker in Réunion, France. Plant Dis. 2017, 101, 503. [Google Scholar] [CrossRef]
  26. Aubert, B. Vergers de la Réunion et de l’Océan Indien. In Hommes et Fruits en Pays du Sud; CIRAD: Montpellier, France, 2014; pp. 111–166. [Google Scholar]
  27. Grygiel, P.; Seny-Couty, A.; Abdou Hassani, F.; Boyer, C.; Boyer, K.; Vernière, C.; Pruvost, O.; Hamza, A.A. First report of Xanthomonas citri pv. citri pathotype A causing Asiatic citrus canker in Grande Comore and Anjouan. Plant Dis. 2014, 98, 1739. [Google Scholar] [CrossRef] [Green Version]
  28. Pruvost, O.; Roumagnac, P.; Gaube, C.; Chiroleu, F.; Gagnevin, L. New media for the semi-selective isolation and enumeration of Xanthomonas campestris pv. mangiferaeindicae, the causal agent of mango bacterial black spot. J. Appl. Microbiol. 2005, 99, 803–815. [Google Scholar] [CrossRef]
  29. Behlau, F.; Canteros, B.I.; Jones, J.B.; Graham, J.H. Copper resistance genes from different xanthomonads and citrus epiphytic bacteria confer resistance to Xanthomonas citri subsp citri. Eur. J. Plant Pathol. 2012, 133, 949–963. [Google Scholar] [CrossRef]
  30. Da Silva, A.C.; Ferro, J.A.; Reinach, F.C.; Farah, C.S.; Furlan, L.R.; Quaggio, R.B.; Monteiro-Vitorello, C.B.; Van Sluys, M.A.; Almeida, N.F.; Alves, L.M.; et al. Comparison of the genomes of two Xanthomonas pathogens with differing host specificities. Nature 2002, 417, 459–463. [Google Scholar] [CrossRef]
  31. Kamvar, Z.N.; Tabima, J.F.; Grünwald, N.J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2014, 2, e281. [Google Scholar] [CrossRef] [Green Version]
  32. Goudet, J. Hierfstat, a package for R to compute and test hierarchical F-statistics. Mol. Ecol. Notes 2005, 5, 184–186. [Google Scholar] [CrossRef] [Green Version]
  33. Francisco, A.P.; Vaz, C.; Monteiro, P.T.; Melo-Cristino, J.; Ramirez, M.; Carriço, J.A. PHYLOViZ: Phylogenetic inference and data visualization for sequence based typing methods. BMC Bioinform. 2012, 13, 87. [Google Scholar] [CrossRef] [Green Version]
  34. Jombart, T.; Devillard, S.; Balloux, F. Discriminant analysis of principal components: A new method for the analysis of genetically structured populations. BMC Genet. 2010, 11, 94. [Google Scholar] [CrossRef] [Green Version]
  35. Jombart, T. adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics 2008, 24, 1403–1405. [Google Scholar] [CrossRef] [Green Version]
  36. Robene, I.; Maillot-Lebon, V.; Chabirand, A.; Moreau, A.; Becker, N.; Moumene, A.; Rieux, A.; Campos, P.; Gagnevin, L.; Gaudeul, M.; et al. Development and comparative validation of genomic-driven PCR-based assays to detect Xanthomonas citri pv. citri in citrus plants. BMC Microbiol. 2020, 20, 296. [Google Scholar] [CrossRef]
  37. Favaro, M.A.; Micheloud, N.G.; Roeschlin, R.A.; Chiesa, M.A.; Castagnaro, A.P.; Vojnov, A.A.; Gmitter, F.G., Jr.; Gadea, J.; Rista, L.M.; Gariglio, N.F.; et al. Surface barriers of mandarin ‘Okitsu’ leaves make a major contribution to canker disease resistance. Phytopathology 2014, 104, 970–976. [Google Scholar] [CrossRef] [Green Version]
  38. Reed, H.S.; HIrano, E. The density of stomata in citrus leaves. J. Agric. Res. 1931, 43, 209–222. [Google Scholar]
  39. Simko, I.; Piepho, H.P. The Area Under the Disease Progress Stairs: Calculation, advantage, and application. Phytopathology 2012, 102, 381–389. [Google Scholar] [CrossRef] [Green Version]
  40. Brooks, M.E.; Kristensen, K.; Van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Machler, M.; Bolker, B.M. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017, 9, 378–400. [Google Scholar] [CrossRef] [Green Version]
  41. Fox, J.; Weisberg, S. An R Companion to Applied Regression, 3rd ed.; SAGE Publications, Inc.: Thousand Oaks, CA, USA, 2019; p. 608. [Google Scholar]
  42. Scherm, H.; Ojiambo, P.S. Applications of survival analysis in botanical epidemiology. Phytopathology 2004, 94, 1022–1026. [Google Scholar] [CrossRef] [Green Version]
  43. Therneau, T.M.; Grambsch, P.M. Modeling Survival Data: Extending the Cox Model; Springer: New York, NY, USA, 2000. [Google Scholar]
  44. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: New York, NY, USA, 2002; p. 495. [Google Scholar]
  45. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous inference in general parametric models. Biom. J. 2008, 50, 346–363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Escalon, A.; Javegny, S.; Vernière, C.; Noël, L.D.; Vital, K.; Poussier, S.; Hajri, A.; Boureau, T.; Pruvost, O.; Arlat, M.; et al. Variations in type III effector repertoires, pathological phenotypes and host range of Xanthomonas citri pv. citri pathotypes. Mol. Plant Pathol. 2013, 14, 483–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Pruvost, O.; Boher, B.; Brocherieux, C.; Nicole, M.; Chiroleu, F. Survival of Xanthomonas axonopodis pv. citri in leaf lesions under tropical environmental conditions and simulated splash dispersal of inoculum. Phytopathology 2002, 92, 336–346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Timmer, L.W.; Gottwald, T.R.; Zitko, S.E. Bacterial exudation from lesions of asiatic citrus canker and citrus bacterial spot. Plant Dis. 1991, 75, 192–195. [Google Scholar] [CrossRef]
  49. Richard, D.; Rieux, A.; Lefeuvre, P.; Hamza, A.; Lobin, K.K.; Naiken, M.; Stravens, R.; Boyer, C.; Boyer, K.; Javegny, S.; et al. Draft genome sequences of 284 Xanthomonas citri pv. citri strains causing Asiatic citrus canker. Microbiol. Resour. Announc. 2021, 10, 1. [Google Scholar] [CrossRef]
  50. Estoup, A.; Jarne, P.; Cornuet, J.M. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol. Ecol. 2002, 11, 1591–1604. [Google Scholar] [CrossRef]
  51. Ellegren, H. Microsatellites: Simple sequences with complex evolution. Nat. Rev. Genet. 2004, 5, 435–445. [Google Scholar] [CrossRef]
  52. Favaro, M.A.; Molina, M.C.; Roeschlin, R.A.; Gadea, J.; Gariglio, N.; Marano, M.R. Different responses in mandarin cultivars uncover a role of cuticular waxes in the resistance to citrus canker. Phytopathology 2020, 110, 1791–1801. [Google Scholar] [CrossRef]
  53. Newberry, E.A.; Bhandari, R.; Minsavage, G.V.; Timilsina, S.; Jibrin, M.O.; Kemble, J.; Sikora, E.J.; Jones, J.B.; Potnis, N. Independent evolution with the gene flux originating from multiple Xanthomonas species explains genomic heterogeneity in Xanthomonas perforans. Appl. Environ. Microbiol. 2019, 85, e00885-19. [Google Scholar] [CrossRef]
  54. Newberry, E.; Bhandari, R.; Kemble, J.; Sikora, E.; Potnis, N. Genome-resolved metagenomics to study co-occurrence patterns and intraspecific heterogeneity among plant pathogen metapopulations. Environ. Microbiol. 2020, 22, 2693–2708. [Google Scholar] [CrossRef]
Microorganisms 09 00945 g001 550
Figure 1. Geographical location and sample numbers of the Xanthomonas citri pv. citri strain collection analysed in the present study.
Figure 1. Geographical location and sample numbers of the Xanthomonas citri pv. citri strain collection analysed in the present study.
Microorganisms 09 00945 g001
Microorganisms 09 00945 g002 550
Figure 2. Genetic structure of Xanthomonas citri pv. citri lineage 1 in the SWIO region based on the discriminant analysis of principal components (DAPC) of microsatellite data. Numbers and colours represent the four genetic clusters retained from Bayesian information criterion (BIC) values. Strains from the Comoros and the Mascarenes grouped in a single genetic cluster (GC1), whereas strains from the Seychelles grouped in three distinct clusters. (A) Scatterplot representing axes 1 and 2 of the DAPC. (B) Scatterplot representing axes 1 and 3 of the DAPC. A subset of haplotypes originating from the Comoros and the Mascarenes (GC1, minor clusters and singletons) were used to avoid too large within-cluster variance (large haplotype number disequilibrium among clusters).
Figure 2. Genetic structure of Xanthomonas citri pv. citri lineage 1 in the SWIO region based on the discriminant analysis of principal components (DAPC) of microsatellite data. Numbers and colours represent the four genetic clusters retained from Bayesian information criterion (BIC) values. Strains from the Comoros and the Mascarenes grouped in a single genetic cluster (GC1), whereas strains from the Seychelles grouped in three distinct clusters. (A) Scatterplot representing axes 1 and 2 of the DAPC. (B) Scatterplot representing axes 1 and 3 of the DAPC. A subset of haplotypes originating from the Comoros and the Mascarenes (GC1, minor clusters and singletons) were used to avoid too large within-cluster variance (large haplotype number disequilibrium among clusters).
Microorganisms 09 00945 g002
Microorganisms 09 00945 g003 550
Figure 3. Categorical minimum spanning tree from MLVA-14 data showing the genetic relatedness among strains of Xanthomonas citri pv. citri in the Comoros. Dots represent haplotypes. Dot diameter and colour are representative of the number of strains per haplotype and island of isolation, respectively (red = Anjouan; blue = Grande Comore; green = Mayotte). Single, double and triple locus variations are represented as thick solid, thin solid and dotted lines joining haplotypes, respectively.
Figure 3. Categorical minimum spanning tree from MLVA-14 data showing the genetic relatedness among strains of Xanthomonas citri pv. citri in the Comoros. Dots represent haplotypes. Dot diameter and colour are representative of the number of strains per haplotype and island of isolation, respectively (red = Anjouan; blue = Grande Comore; green = Mayotte). Single, double and triple locus variations are represented as thick solid, thin solid and dotted lines joining haplotypes, respectively.
Microorganisms 09 00945 g003
Table
Table 1. Genetic diversity parameters calculated from microsatellite data (MLVA-14) for X. citri pv. citri pathotype A in the SWIO region.
Table 1. Genetic diversity parameters calculated from microsatellite data (MLVA-14) for X. citri pv. citri pathotype A in the SWIO region.
Island N 1 NH 2 A 3 AP 4 HE 5 Cluster 6
The Comoros
Anjouan21193.7860.0000.430GC1 (95%)
Grande Comore78513.6040.0000.490GC1 (100%)
Mayotte92563.1170.1850.439GC1 (99%)
Moheli86785.6440.0820.666GC1 (74%)
The Mascarenes
Mauritius133734.0130.0830.468GC1 (95%)
Réunion 713207895.2890.4550.633GC1 (98 %)
Rodrigues73635.2650.0140.649GC1 (63%)
The Seychelles
Mahé82463.5120.5490.680GC2 (54%)
GC3 (30%)
GC4 (16%)
1 Number of isolates per island; 2 number of haplotypes per island; 3 allelic richness calculated by rarefaction method (n = 21); 4 private allelic richness; 5 Nei’s genetic diversity; 6 cluster number derived from analysis of microsatellite data; 7 data primarily from [23]. Values in brackets indicate the percentage of haplotypes assigned to GCs. Other strains were assigned to minor GCs or singletons (see Materials and Methods section for details). Strains assigned to pathotype A* (genetic lineage 4) were not considered.
Table
Table 2. Genetic cluster and genetic lineage assignation of strains of Xanthomonas citri pv. citri originating from the SWIO region based on microsatellite and minisatellite genotyping data, respectively.
Table 2. Genetic cluster and genetic lineage assignation of strains of Xanthomonas citri pv. citri originating from the SWIO region based on microsatellite and minisatellite genotyping data, respectively.
Cluster (Microsatellite Data)Country (Archipelago)SourcenGenetic Lineage (Minisatellite Data) 1
GC1Anjouan (The Comoros)This study161
GC1Grande Comore (The Comoros)This study51
GC1Mayotte (The Comoros)This study561
GC1Moheli (The Comoros)This study221
GC1Mauritius (The Mascarenes)This study251
GC1Rodrigues (The Mascarenes)This study111
GC1Réunion (The Mascarenes)This study251
GC1Réunion (The Mascarenes)[23]1621
GC2Mahé (The Seychelles)This study131
GC3Mahé (The Seychelles)This study61
GC4Mahé (The Seychelles)This study51
minor/singletonAnjouan (The Comoros)This study11
minor/singletonMayotte (The Comoros)This study11
minor/singletonMoheli (The Comoros)This study181
minor/singletonMoheli (The Comoros)This study14
minor/singletonMauritius (The Mascarenes)This study41
minor/singletonMauritius (The Mascarenes)This study24
minor/singletonRodrigues (The Mascarenes)This study241
minor/singletonRéunion (The Mascarenes)This study101
minor/singletonRéunion (The Mascarenes)This study14
minor/singletonRéunion (The Mascarenes)[23]51
1 Assignment to genetic lineages was achieved using Discriminant Analysis of Principal Components (DAPC), as reported earlier [16], lineage 1 = pathotype A, lineage 4 = pathotype A*.
Table
Table 3. Comparisons among strains of exuded or total Xanthomonas citri pv. citri population sizes enumerated from single canker lesions 25 days after spray inoculation of three host species.
Table 3. Comparisons among strains of exuded or total Xanthomonas citri pv. citri population sizes enumerated from single canker lesions 25 days after spray inoculation of three host species.
CitronMandarinSweet Orange
ExudedTotalExudedTotalExudedTotal
ANOVA probability0.0020.0040.006< 0.001< 0.0010.025
GCStrainLog-transformed cfu lesion−1 compared using Tukey’s tests
1LH2416.69 ab7.43 a6.20 ab7.33 abc6.97 ab7.43 ab
1LM089-417.00 a7.41 a6.63 a7.62 a6.88 abc7.27 ab
1LN005-46.65 ab7.08 ab6.20 ab6.99 abc6.91 abc7.58 a
1LN007-36.47 b7.30 ab6.30 ab7.21 abc6.79 abc7.47 ab
2JZ0926.45 b7.12 ab6.29 ab7.19 abc6.98 abc7.32 ab
2LB100-16.56 ab7.14 ab5.93 ab6.79 bc6.52 c6.90 b
2LJ0016.54 b6.92 ab6.50 ab7.49 ab6.72 abc7.07 ab
2LP029-156.80 ab7.09 ab5.65 ab6.60 c6.39 c7.19 ab
3JZ0946.16 b6.75 ab6.12 ab7.13 abc6.70 abc7.01 ab
3LP027-36.75 ab7.28 ab5.64 ab6.49 c6.88 abc7.38 ab
3LP027-56.32 b6.64 b6.49 ab7.07 abc6.60 bc7.23 ab
3LP027-136.58 ab7.06 ab5.85 ab6.89 abc6.52 c7.12 ab
4LP028-26.84 ab7.36 a5.60 b6.57 c6.46 bc7.17 ab
4LP028-36.73 ab7.15 ab5.70 ab6.64 c7.02 a7.39 ab
4LP028-56.66 ab7.00 ab5.92 ab6.63 bc6.63 abc7.34 ab
4LP028-66.80 ab7.16 ab6.11 ab6.82 c6.69 abc7.31 ab
Population sizes were determined on KC semiselective medium. Nine lesions per strain–host combination were used for assessing population size. Values (i.e., means) followed by the same letter(s) are not significantly different (p = 0.05) based on Tukey’s tests.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pruvost, O.; Richard, D.; Boyer, K.; Javegny, S.; Boyer, C.; Chiroleu, F.; Grygiel, P.; Parvedy, E.; Robène, I.; Maillot-Lebon, V.; et al. Diversity and Geographical Structure of Xanthomonas citri pv. citri on Citrus in the South West Indian Ocean Region. Microorganisms 2021, 9, 945. https://doi.org/10.3390/microorganisms9050945

AMA Style

Pruvost O, Richard D, Boyer K, Javegny S, Boyer C, Chiroleu F, Grygiel P, Parvedy E, Robène I, Maillot-Lebon V, et al. Diversity and Geographical Structure of Xanthomonas citri pv. citri on Citrus in the South West Indian Ocean Region. Microorganisms. 2021; 9(5):945. https://doi.org/10.3390/microorganisms9050945

Chicago/Turabian Style

Pruvost, Olivier, Damien Richard, Karine Boyer, Stéphanie Javegny, Claudine Boyer, Frédéric Chiroleu, Pierre Grygiel, Evelyne Parvedy, Isabelle Robène, Véronique Maillot-Lebon, and et al. 2021. "Diversity and Geographical Structure of Xanthomonas citri pv. citri on Citrus in the South West Indian Ocean Region" Microorganisms 9, no. 5: 945. https://doi.org/10.3390/microorganisms9050945

APA Style

Pruvost, O., Richard, D., Boyer, K., Javegny, S., Boyer, C., Chiroleu, F., Grygiel, P., Parvedy, E., Robène, I., Maillot-Lebon, V., Hamza, A., Lobin, K. K., Naiken, M., & Vernière, C. (2021). Diversity and Geographical Structure of Xanthomonas citri pv. citri on Citrus in the South West Indian Ocean Region. Microorganisms, 9(5), 945. https://doi.org/10.3390/microorganisms9050945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop