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

Characterization of Pseudomonas syringae Strains Associated with Shoot Blight of Raspberry and Blackberry in Serbia

During May 2016, severe blight symptoms were observed in several raspberry and blackberry fields in Serbia. In total, 22 strains were isolated: 16 from symptomatic raspberry shoots, 2 from asymptomatic raspberry leaves, and 4 from symptomatic blackberry shoots. Additionally, eight raspberry strains, isolated earlier from two similar outbreaks, were included in the study. Pathogenicity of the strains was confirmed on detached raspberry and blackberry shoots by reproducing the symptoms of natural infection. The strains were Gram-negative, fluorescent on King's medium B, ice nucleation positive, and utilized glucose oxidatively. All strains were levan positive, oxidase negative, nonpectolytic, arginine dihydrolase negative, and induced hypersensitivity in tobacco leaves (LOPAT + - - - +, Pseudomonas group Ia). Furthermore, all strains liquefied gelatin and hydrolyzed aesculin but did not show tyrosinase activity or utilize tartrate (GATTa + + - -). Tentative identification using morphology, LOPAT, GATTa, and ice-nucleating ability tests suggested that isolated strains belong to Pseudomonas syringae. The syrB gene associated with syringomycin production was detected in all strains. DNA fingerprints with REP, ERIC, and BOX primers generated identical profiles for 29 strains, except for strain KBI 222, which showed a unique genomic fingerprint. In all, 9 of 10 selected strains exhibited identical sequences of four housekeeping genes: gyrB, rpoD, gapA, and gltA. Five nucleotide polymorphisms were found in strain KBI 222 at the rpoD gene locus only. In the phylogenetic tree based on a concatenated sequence of all four housekeeping genes, strains clustered within phylogroup 2 (i.e., genomospecies 1) of the P. syringae species complex, with pathotype strains of P. syringae pv. aceris and P. syringae pv. solidagae as their closest relatives. There was no correlation between genotype and geographic origin, particular outbreak, host, or cultivar.

Bacteriocin Production Correlates with Epidemiological Prevalence of Phylotype I Sequevar 18 Ralstonia pseudosolanacearum in Madagascar

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is a major threat to vegetable crops in Madagascar. For more effective disease management, surveys were carried out in the main vegetable production areas of the country, leading to the collection of 401 new RSSC isolates. Phylogenetic assignment of the isolates revealed a high prevalence of phylotype I sequevar 18. This result contrasts sharply with the epidemiological pattern of RSSC in neighboring islands, including Reunion Island, Comoros, Mayotte, Mauritius, Rodrigues, and the Seychelles, where phylotype I sequevar 31 is widespread. Molecular typing characterization of the Malagasy isolates allowed the identification of 96 haplotypes. Some are found in various plots located in different provinces, which suggests that they were probably disseminated via infected plant material. To find out a potential explanation for the observed epidemiological pattern, we examined the capacity of the Malagasy strains to produce bacteriocin. Interestingly, the highly prevalent genetic lineages I-18 produce bacteriocins that are active against all the genetic lineages present in the country. This work sheds light on the potential impact of bacteriocins in the epidemiology of Malagasy RSSC. IMPORTANCE Knowledge of the epidemiology of a plant pathogen is essential to develop effective control strategies. This study focuses on the epidemiological pattern of Ralstonia pseudosolanacearum phylotype I populations responsible for bacterial wilt in Madagascar. We identified, with the newly collected isolates in three provinces, four genetic lineages probably propagated via infected plant material in Madagascar. We revealed that the epidemiological situation in Madagascar contrasts with that of neighboring Indian Ocean islands. Interestingly, our study on the bacteriocin-producing capacity of Malagasy isolates revealed a correlation between the inhibitory activity of the producing strains and the observed epidemiology. These results suggested that the epidemiology of plant pathogens may be impacted by bacteriocin production.

Three semi-selective media for Pseudomonas syringae pv. maculicola and P. cannabina pv. alisalensis

Three semi-selective media, DTarTA, SPbc, and SPamt, were developed and tested to isolate Pseudomonas syringae pv. maculicola (Psm) and P. cannabina pv. alisalensis (Pca) from Raphanus sativus seeds. DTarTA contained D-tartaric acid as a carbon source and potassium tellurite, ampicillin sodium, and methyl violet as antibiotics. DTarTA suppressed growth in 19 of the 24 pathovars from the P. syringae complex, whereas Psm and Pca grew and formed gray to black colonies. SPamt contained sucrose and peptone as nutrient sources and was supplemented with bromothymol blue and the same antibiotics present in DTarTA and Psm and Pca formed yellowish to dark brown colonies on the SPamt medium. SPbc contained sucrose and peptone and was supplemented with cephalexin and boric acid as antibiotics and Psm and Pca formed semi-translucent to white colonies on the SPbc medium. SPamt and SPbc suppressed the growth of several plant-associated bacteria (except the P. syringae complex). The growth of saprophytic bacteria in seeds on the different media was compared with that on King's B medium, using five types of commercially available Raphanus sativus seeds. The suppression rate of DTarTA was 85-99% and was lower for seeds with more saprophytic bacteria. The suppression rates of SPamt and SPbc were 90-99%. In detection tests using 10,000 seed samples mixed with Pca or Psm-contaminated seeds, it was possible to selectively isolate Psm and Pca using SPamt and SPbc, even when the colony numbers of the target bacterium constituted less than 10% of the total colonies. KEY POINTS: • Bacterial leaf spot and blight pathogens were selectively isolated from seeds. • DTarTA medium distinguishes these pathogens from P. syringae complex pathovars. • SPamp and SPbc media have different selectivity for plant-associated bacteria.

Characterization and Genetic Diversity of Pseudomonas syringae pv. syringae Isolates Associated with Rice Bacterial Leaf Spot in Heilongjiang, China

Simple Summary

In Northeast China, rice bacterial brown leaf spot caused by Pseudomonas syringae pv. syringae is among the most damaging rice diseases. This disease results in significant yield losses. This study focuses on a comprehensive analysis of the pathogen, population structure, and genetic diversity within the species based on biochemical tests and genetic characterizations. Our results indicate high genetic heterogeneity in Pseudomonas syringae pv. syringae isolates, and clustering of testing isolates and reference strains are related with the genomospecies 1. This work contributes to the physiological classification of the P. s. pv. syringae isolated from Heilongjiang Province, China and the results present new data concerning the phylogeny and genetic diversity. Such studies have not been reported about P. s. pv. syringae from this region yet.

Abstract

In China, rice is one of the most important cereal crops. Rice bacterial brown leaf spot caused by P. s. pv. syringae is among the most damaging rice diseases in the Heilongjiang Province of China and results in substantial yield losses. In this study, a comprehensive analysis of the pathogen, population structure, and genetic diversity within the species was performed. For this purpose, 176 bacterial isolates of P. s. pv. syringae collected from 15 locations were characterized by using biochemical tests such as the LOPAT test, and genetic characterizations such as multilocus sequence analysis (MLSA) and repetitive PCR, using BOX, REP and ERIC primers. Biochemical testing and detection of syrB genes confirm the presence of P. s. pv. syringae, genetic characterization by MLSA and genetic fingerprinting by repetitive PCR confirmed that high genetic heterogeneity exists in the P. s. pv. syringae isolates, and clustering of the tested isolates and reference strains are related with the same genomospecies 1. This work contributes to the physiological classification of the P. s. pv. syringae isolated from Heilongjiang Province, China, and the results present new data concerning the phylogeny and genetic diversity. This type of study about P. s. pv. syringae has been not reported from this region until now.

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.

Primers for specific detection and identification of Pseudomonas syringae pv. maculicola and P. cannabina pv. alisalensis

Bacterial leaf spot and bacterial leaf blight are global threats to the cultivation of cruciferous vegetables, and it is necessary to develop methods to easily detect, identify, and distinguish the causative pathogens Pseudomonas syringae pv. maculicola (Psm) and P. cannabina pv. alisalensis (Pca). Here, we used the sequence specificity of the exchangeable effector loci flanking the hrp gene cluster to design primers that can help detect and discriminate between Psm and Pca. Primers common to both bacteria (hrpK_fw1 and hrpK_fw2) were designed within hrpK at the end of the hrp gene cluster. Psm-specific primers (MAC_rv1 and MAC_rv2) were designed in hopPtoB1 and Pca-specific primers (ALS_rv1 and ALS_rv2) were designed in hopX1 adjacent to hrpK. PCR using hrpK_fw1 and MAC_rv1 or hrpK_fw2 and MAC_rv2 amplified DNA fragments of only Psm, P. syringae pv. tomato (causal agent of tomato bacterial speck), and P. syringae pv. spinaciae (causal agent of spinach bacterial leaf spot), among 76 strains of phytopathogenic bacteria. PCR using hrpK_fw1 and ALS_rv1 or hrpK_2 and ALS_rv2 amplified DNA fragments of only Pca. Multiplex PCR with these primers could easily distinguish Psm and Pca from bacterial colonies isolated on growth media and detect the pathogen in symptomatic leaves. Multiplex nested PCR with the primers detected contamination in one Psm- and/or one Pca-infected seeds in 1000 seeds. These results suggest that these PCR primers could help detect and discriminate Psm and Pca. KEY POINTS: • We investigated Pseudomonas syringae pv. maculicola and P. cannabina pv. alisalensis. • Novel primers common to both bacteria were designed following genome comparison. • Multiplex PCR with new primers could discriminate Psm and Pca.

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.

Genetic Diversity of Pseudomonas syringae pv. actinidiae: Seasonal and Spatial Population Dynamics

Pseudomonassyringae pv. actinidiae (Psa) is a gram-negative bacterium responsible for the bacterial canker in Actinidia chinensis var. deliciosa and A. chinensis var. chinensis, a quarantine organism threatening the kiwifruit industry sustainability. The present study aimed to determine the genetic structure of the endophytic and epiphytic populations of Psa isolated from four different Portuguese orchards with distinct abiotic conditions in two consecutive seasons. The results identified several coexisting and highly heterogeneous Psa populations. Moreover, evident changes in population structure occurred between the epiphytic and endophytic populations, and between seasons with a notable decrease in Psa diversity in autumn. This work provided solid evidence that the initial clonal expansion of Psa in Europe was followed by a wide genomic diversification. This perspective is important for the understanding of kiwifruit bacterial canker disease occurrence and Psa evolution, namely when adopting strategies for management of epidemics.

Host Range Determinants of Pseudomonas savastanoi Pathovars of Woody Hosts Revealed by Comparative Genomics and Cross-Pathogenicity Tests

The study of host range determinants within the Pseudomonas syringae complex is gaining renewed attention due to its widespread distribution in non-agricultural environments, evidence of large variability in intra-pathovar host range, and the emergence of new epidemic diseases. This requires the establishment of appropriate model pathosystems facilitating integration of phenotypic, genomic and evolutionary data. Pseudomonas savastanoi pv. savastanoi is a model pathogen of the olive tree, and here we report a closed genome of strain NCPPB 3335, plus draft genome sequences of three strains isolated from oleander (pv. nerii), ash (pv. fraxini) and broom plants (pv. retacarpa). We then conducted a comparative genomic analysis of these four new genomes plus 16 publicly available genomes, representing 20 strains of these four P. savastanoi pathovars of woody hosts. Despite overlapping host ranges, cross-pathogenicity tests using four plant hosts clearly separated these pathovars and lead to pathovar reassignment of two strains. Critically, these functional assays were pivotal to reconcile phylogeny with host range and to define pathovar-specific genes repertoires. We report a pan-genome of 7,953 ortholog gene families and a total of 45 type III secretion system effector genes, including 24 core genes, four genes exclusive of pv. retacarpa and several genes encoding pathovar-specific truncations. Noticeably, the four pathovars corresponded with well-defined genetic lineages, with core genome phylogeny and hierarchical clustering of effector genes closely correlating with pathogenic specialization. Knot-inducing pathovars encode genes absent in the canker-inducing pv. fraxini, such as those related to indole acetic acid, cytokinins, rhizobitoxine, and a bacteriophytochrome. Other pathovar-exclusive genes encode type I, type II, type IV, and type VI secretion system proteins, the phytotoxine phevamine A, a siderophore, c-di-GMP-related proteins, methyl chemotaxis proteins, and a broad collection of transcriptional regulators and transporters of eight different superfamilies. Our combination of pathogenicity analyses and genomics tools allowed us to correctly assign strains to pathovars and to propose a repertoire of host range-related genes in the P. syringae complex.

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.

New Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Scheme for Fine-Scale Monitoring and Microevolution-Related Study of Ralstonia pseudosolanacearum Phylotype I Populations

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is considered one of the most harmful plant diseases in the world. Special attention should be paid to R. pseudosolanacearum phylotype I due to its large host range, its worldwide distribution, and its high evolutionary potential. So far, the molecular epidemiology and population genetics of this bacterium are poorly understood. Until now, the genetic structure of the RSSC has been analyzed on the worldwide and regional scales. Emerging questions regarding evolutionary forces in RSSC adaptation to hosts now require genetic markers that are able to monitor RSSC field populations. In this study, we aimed to evaluate the multilocus variable-number tandem-repeat analysis (MLVA) approach for its ability to discriminate genetically close phylotype I strains and for population genetics studies. We developed a new MLVA scheme (MLVA-7) allowing us to genotype 580 R. pseudosolanacearum phylotype I strains extracted from susceptible and resistant hosts and from different habitats (stem, soil, and rhizosphere). Based on specificity, polymorphism, and the amplification success rate, we selected seven fast-evolving variable-number tandem-repeat (VNTR) markers. The newly developed MLVA-7 scheme showed higher discriminatory power than the previously published MLVA-13 scheme when applied to collections sampled from the same location on different dates and to collections from different locations on very small scales. Our study provides a valuable tool for fine-scale monitoring and microevolution-related study of R. pseudosolanacearum phylotype I populations.IMPORTANCE Understanding the evolutionary dynamics of adaptation of plant pathogens to new hosts or ecological niches has become a key point for the development of innovative disease management strategies, including durable resistance. Whereas the molecular mechanisms underlying virulence or pathogenicity changes have been studied thoroughly, the population genetics of plant pathogen adaptation remains an open, unexplored field, especially for plant-pathogenic bacteria. MLVA has become increasingly popular for epidemiosurveillance and molecular epidemiology studies of plant pathogens. However, this method has been used mostly for genotyping and identification on a regional or global scale. In this study, we developed a new MLVA scheme, targeting phylotype I of the soilborne Ralstonia solanacearum species complex (RSSC), specifically to address the bacterial population genetics on the field scale. Such a MLVA scheme, based on fast-evolving loci, may be a tool of choice for field experimental evolution and spatial genetics studies.

Characterization of Pseudomonas syringae pv. syringae, Causal Agent of Citrus Blast of Mandarin in Montenegro

Citrus blast caused by bacterium Pseudomonas syringae is a very important disease of citrus occuring in many areas of the world, but with few data about genetic structure of the pathogen involved. Considering the above fact, this study reports genetic characterization of 43 P. syringae isolates obtained from plant tissue displaying citrus blast symptoms on mandarin (Citrus reticulata) in Montenegro, using multilocus sequence analysis of gyrB, rpoD, and gap1 gene sequences. Gene sequences from a collection of 54 reference pathotype strains of P. syringae from the Plant Associated and Environmental Microbes Database (PAMDB) was used to establish a genetic relationship with our isolates obtained from mandarin. Phylogenetic analyses of gyrB, rpoD, and gap1 gene sequences showed that P. syringae pv. syringae causes citrus blast in mandarin in Montenegro, and belongs to genomospecies 1. Genetic homogeneity of isolates suggested that the Montenegrian population might be clonal which indicates a possible common source of infection. These findings may assist in further epidemiological studies of this pathogen and for determining mandarin breeding strategies for P. syringae control.

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.

Phylogenetic and Variable-Number Tandem-Repeat Analyses Identify Nonpathogenic Xanthomonas arboricola Lineages Lacking the Canonical Type III Secretion System

Xanthomonas arboricola is conventionally known as a taxon of plant-pathogenic bacteria that includes seven pathovars. This study showed that X. arboricola also encompasses nonpathogenic bacteria that cause no apparent disease symptoms on their hosts. The aim of this study was to assess the X. arboricola population structure associated with walnut, including nonpathogenic strains, in order to gain a better understanding of the role of nonpathogenic xanthomonads in walnut microbiota. A multilocus sequence analysis (MLSA) was performed on a collection of 100 X. arboricola strains, including 27 nonpathogenic strains isolated from walnut. Nonpathogenic strains grouped outside clusters defined by pathovars and formed separate genetic lineages. A multilocus variable-number tandem-repeat analysis (MLVA) conducted on a collection of X. arboricola strains isolated from walnut showed that nonpathogenic strains clustered separately from clonal complexes containing Xanthomonas arboricola pv. juglandis strains. Some nonpathogenic strains of X. arboricola did not contain the canonical type III secretion system (T3SS) and harbored only one to three type III effector (T3E) genes. In the nonpathogenic strains CFBP 7640 and CFBP 7653, neither T3SS genes nor any of the analyzed T3E genes were detected. This finding raises a question about the origin of nonpathogenic strains and the evolution of plant pathogenicity in X. arboricola. T3E genes that were not detected in any nonpathogenic isolates studied represent excellent candidates to be those responsible for pathogenicity in X. arboricola.

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.

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

Knowing the identity of bacterial plant pathogens is essential to strategic and sustainable disease management in agricultural systems. This knowledge is critical for growers, diagnosticians, extension agents, and others dealing with crops. However, such identifications are linked to bacterial taxonomy, a complicated and changing discipline that depends on methods and information that are often not used by those who are diagnosing field problems. Modern molecular tools for fingerprinting and sequencing allow for pathogen identification in the absence of distinguishing or conveniently tested phenotypic characteristics. These methods are also useful in studying the etiology and epidemiology of phytopathogenic bacteria from epidemics, as was done in numerous studies conducted in California's Salinas Valley. Multilocus and whole-genome sequence analyses are becoming the cornerstones of studies of microbial diversity and bacterial taxonomy. Whole-genome sequence analysis needs to become adequately accessible, automated, and affordable in order to be used routinely for identification and epidemiology. The power of molecular tools in accurately identifying bacterial pathogenesis is therefore of value to the farmer, diagnostician, phytobacteriologist, and taxonomist.

Evaluating the impacts of stressors of Pseudomonas syringae pathovar tomato on the effectiveness of multi-locus variable number tandem repeat analysis and multi-locus sequence typing in microbial forensic investigations

BACKGROUND

Crops in the USA are vulnerable to natural and criminal threats because of their widespread cultivation and lack of surveillance, and because of implementation of growing practices such as monoculture. To prepare for investigation and attribution of such events, forensic assays, including determination of molecular profiles, are being adapted for use with plant pathogens. The use of multi-locus variable number tandem repeat (VNTR) analysis (MLVA) and multi-locus sequence typing (MLST) in investigations involving plant pathogens may be problematic because the long lag periods between pathogen introduction and discovery of associated disease may provide enough time for evolution to occur in the regions of the genome employed in each assay. Thus, more information on the stability of the loci employed in these methods is needed.

RESULTS

The MLVA fingerprints and MLST profiles were consistent throughout the experiment, indicating that, using a specific set of primers and conditions, MLVA and MLST typing systems reliably identify P.s. tomato DC3000. This information is essential to forensic investigators in interpreting comparisons between MLVA and MLST typing profiles observed in P.s. tomato isolates.

CONCLUSIONS

Our results indicate that MLVA and MLST typing systems, utilizing the specified primers and conditions, could be employed successfully in forensics investigations involving P.s. tomato. Similar experiments should be conducted in the field and with other high-consequence plant pathogens to ensure that the assays are reliable for pathogens infecting plants in their natural environment and for organisms that may display faster rates of mutation.

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.

Genomics-Based Exploration of Virulence Determinants and Host-Specific Adaptations of Pseudomonas syringae Strains Isolated from Grasses

The Pseudomonas syringae species complex has recently been named the number one plant pathogen, due to its economic and environmental impacts, as well as for its role in scientific research. The bacterium has been repeatedly reported to cause outbreaks on bean, cucumber, stone fruit, kiwi and olive tree, as well as on other crop and non-crop plants. It also serves as a model organism for research on the Type III secretion system (T3SS) and plant-pathogen interactions. While most of the current work on this pathogen is either carried out on one of three model strains found on dicot plants with completely sequenced genomes or on isolates obtained from recent outbreaks, not much is known about strains isolated from grasses (Poaceae). Here, we use comparative genomics in order to identify putative virulence-associated genes and other Poaceae-specific adaptations in several newly available genome sequences of strains isolated from grass species. All strains possess only a small number of known Type III effectors, therefore pointing to the importance of non-Type III secreted virulence factors. The implications of this finding are discussed.

Genotyping methods

Genotyping allows for the identification of bacterial isolates to the strain level and provides basic information about the evolutionary biology, population biology, taxonomy, ecology, and genetics of bacteria. Depending on the underlying question and available resources, Pseudomonas aeruginosa strains may be typed by anonymous fingerprinting techniques or electronically portable sequence-based typing methods such as multiple locus variable number tandem repeat (VNTR) analysis (MLVA), multilocus sequence typing, or oligonucleotide microarray. Macrorestriction fragment pattern analysis is a genotyping method that is globally applicable to all bacteria and hence has been and still is the reference method for strain typing in bacteriology. Agarose-embedded chromosomal DNA is cleaved with a rare-cutting restriction endonuclease and the generated 20-70 fragments are then separated by pulsed-field gel electrophoresis. The chapter provides a detailed step-by-step manual for SpeI genome fingerprinting of Pseudomonas chromosomes that has been optimized for SpeI fragment pattern analysis of P. aeruginosa.

Genomic survey of pathogenicity determinants and VNTR markers in the cassava bacterial pathogen Xanthomonas axonopodis pv. Manihotis strain CIO151

Xanthomonas axonopodis pv. manihotis (Xam) is the causal agent of bacterial blight of cassava, which is among the main components of human diet in Africa and South America. Current information about the molecular pathogenicity factors involved in the infection process of this organism is limited. Previous studies in other bacteria in this genus suggest that advanced draft genome sequences are valuable resources for molecular studies on their interaction with plants and could provide valuable tools for diagnostics and detection. Here we have generated the first manually annotated high-quality draft genome sequence of Xam strain CIO151. Its genomic structure is similar to that of other xanthomonads, especially Xanthomonas euvesicatoria and Xanthomonas citri pv. citri species. Several putative pathogenicity factors were identified, including type III effectors, cell wall-degrading enzymes and clusters encoding protein secretion systems. Specific characteristics in this genome include changes in the xanthomonadin cluster that could explain the lack of typical yellow color in all strains of this pathovar and the presence of 50 regions in the genome with atypical nucleotide composition. The genome sequence was used to predict and evaluate 22 variable number of tandem repeat (VNTR) loci that were subsequently demonstrated as polymorphic in representative Xam strains. Our results demonstrate that Xanthomonas axonopodis pv. manihotis strain CIO151 possesses ten clusters of pathogenicity factors conserved within the genus Xanthomonas. We report 126 genes that are potentially unique to Xam, as well as potential horizontal transfer events in the history of the genome. The relation of these regions with virulence and pathogenicity could explain several aspects of the biology of this pathogen, including its ability to colonize both vascular and non-vascular tissues of cassava plants. A set of 16 robust, polymorphic VNTR loci will be useful to develop a multi-locus VNTR analysis scheme for epidemiological surveillance of this disease.

A Pseudomonas syringae diversity survey reveals a differentiated phylotype of the pathovar syringae associated with the mango host and mangotoxin production

Pseudomonas syringae pv. syringae, the causal agent of bacterial apical necrosis (BAN) in mango crops, has been isolated in different mango-producing areas worldwide. An extensive collection of 87 P. syringae pv. syringae strains isolated from mango trees affected by BAN from different countries, but mainly from Southern Spain, were initially examined by repetitive sequence-based polymerase chain reaction (rep-PCR) to analyze the genetic diversity with an epidemiological aim. rep-PCR was powerful in assessing intrapathovar distribution and also allowing clustering of the P. syringae pv. syringae strains isolated from mango, depending on the isolation area. A clear pattern of clustering was observed for all the P. syringae pv. syringae strains isolated from mango distinct from strains from other hosts, including strains for the same geographical regions as the mango isolates. For this reason, a representative group of 51 P. syringae pv. syringae strains isolated from mango and other hosts, as well as some P. syringae strains from other pathovars, were further characterized to determine their possible genetic, phenotypic, and phylogenetic relationships. Similar to the rep-PCR results, the randomly amplified polymorphic DNA PCR (RAPD-PCR) and catabolic diversity analysis using the Biolog GN2 profile grouped 90% of the mango isolates together in a unique cluster. Interestingly, the majority of P. syringae pv. syringae strains isolated from mango produced mangotoxin. The analysis of the phylogenetic distribution using the multilocus sequence typing analysis strongly supports the existence of a differentiated phylotype of the pathovar syringae mainly associated with the mango host and characterized by the mangotoxin production.

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.

Application of variable-number tandem-repeat typing to discriminate Ralstonia solanacearum strains associated with English watercourses and disease outbreaks

Variable-number tandem-repeat (VNTR) analysis was used for high-resolution discrimination among Ralstonia solanacearum phylotype IIB sequevar 1 (PIIB-1) isolates and further evaluated for use in source tracing. Five tandem-repeat-containing loci (comprising six tandem repeats) discriminated 17 different VNTR profiles among 75 isolates from potato, geranium, bittersweet (Solanum dulcamara), tomato, and the environment. R. solanacearum isolates from crops at three unrelated outbreak sites where river water had been used for irrigation had distinct VNTR profiles that were shared with PIIB-1 isolates from infected bittersweet growing upriver of each site. The VNTR profiling results supported the implication that the source of R. solanacearum at each outbreak was contaminated river water. Analysis of 51 isolates from bittersweet growing in river water at different locations provided a means to evaluate the technique for studying the epidemiology of the pathogen in the environment. Ten different VNTR profiles were identified among bittersweet PIIB-1 isolates from the River Thames. Repeated findings of contiguous river stretches that produced isolates that shared single VNTR profiles supported the hypothesis that the pathogen had disseminated from infected bittersweet plants located upriver. VNTR profiles shared between bittersweet isolates from two widely separated Thames tributaries (River Ray and River Colne) suggested they were independently contaminated with the same clonal type. Some bittersweet isolates had VNTR profiles that were shared with potato isolates collected outside the United Kingdom. It was concluded that VNTR profiling could contribute to further understanding of R. solanacearum epidemiology and assist in control of future disease outbreaks.

Nonagricultural reservoirs contribute to emergence and evolution of Pseudomonas syringae crop pathogens

While the existence of environmental reservoirs of human pathogens is well established, less is known about the role of nonagricultural environments in emergence, evolution, and spread of crop pathogens. Here, we analyzed phylogeny, virulence genes, host range, and aggressiveness of Pseudomonas syringae strains closely related to the tomato pathogen P. syringae pv. tomato (Pto), including strains isolated from snowpack and streams. The population of Pto relatives in nonagricultural environments was estimated to be large and its diversity to be higher than that of the population of Pto and its relatives on crops. Ancestors of environmental strains, Pto, and other genetically monomorphic crop pathogens were inferred to have frequently recombined, suggesting an epidemic population structure for P. syringae. Some environmental strains have repertoires of type III-secreted effectors very similar to Pto, are almost as aggressive on tomato as Pto, but have a wider host range than typical Pto strains. We conclude that crop pathogens may have evolved through a small number of evolutionary events from a population of less aggressive ancestors with a wider host range present in nonagricultural environments.

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.

The mangotoxin biosynthetic operon (mbo) is specifically distributed within Pseudomonas syringae genomospecies 1 and was acquired only once during evolution

Mangotoxin production was first described in Pseudomonas syringae pv. syringae strains. A phenotypic characterization of 94 P. syringae strains was carried out to determine the genetic evolution of the mangotoxin biosynthetic operon (mbo). We designed a PCR primer pair specific for the mbo operon to examine its distribution within the P. syringae complex. These primers amplified a 692-bp DNA fragment from 52 mangotoxin-producing strains and from 7 non-mangotoxin-producing strains that harbor the mbo operon, whereas 35 non-mangotoxin-producing strains did not yield any amplification. This, together with the analysis of draft genomes, allowed the identification of the mbo operon in five pathovars (pathovars aptata, avellanae, japonica, pisi, and syringae), all of which belong to genomospecies 1, suggesting a limited distribution of the mbo genes in the P. syringae complex. Phylogenetic analyses using partial sequences from housekeeping genes differentiated three groups within genomospecies 1. All of the strains containing the mbo operon clustered in groups I and II, whereas those lacking the operon clustered in group III; however, the relative branching order of these three groups is dependent on the genes used to construct the phylogeny. The mbo operon maintains synteny and is inserted in the same genomic location, with high sequence conservation around the insertion point, for all the strains in groups I and II. These data support the idea that the mbo operon was acquired horizontally and only once by the ancestor of groups I and II from genomospecies 1 within the P. syringae complex.

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.

Development of a variable number of tandem repeats typing scheme for the bacterial rice pathogen Xanthomonas oryzae pv. oryzicola

Xanthomonas oryzae pv. oryzicola is an important bacterial pathogen responsible for outbreaks of bacterial leaf streak (BLS) on rice, mostly occurring in Asia and parts of Africa. To better monitor epidemics and assess population structures, efficient tools that allow the precise identification and diagnosis of pathogenic populations are needed. In this study, we explored variable numbers of tandem repeats (VNTR) as a fast, reliable, and cost-effective molecular typing tool. Screening of three X. oryzae pv. oryzicola genome sequences (Philippine strain BLS256, Chinese strain GX01, and Malian strain MAI10) predicted 28 candidate VNTR loci. Primer pairs for polymerase chain reaction (PCR) amplification of all 28 loci were designed and applied to a panel of 20 X. oryzae pv. oryzicola strains originating from Asia and Africa. Sequencing of PCR amplicons revealed 25 robust and polymorphic VNTR loci that are shared among Asian and African X. oryzae pv. oryzicola strains. A dendrogram constructed from 25 VNTR loci indicated that most Asian strains are clearly discriminated from African strains. However, in agreement with previous reports, one strain from Mali is related to Asian strains, pointing to a possible introduction of Asian strains to the African continent. The new VNTR-based tool described here is useful for studies of population structures and epidemiological monitoring of X. oryzae pv. oryzicola.