A Genomic redefinition of Pseudomonas avellanae species

Genetic and phenotypical diversity of Pseudomonas syringae population in the Russian Federation

Abstract Proteobacteria comprising species of Pseudomonas syringae group cause diseases of many plants around the world. The phytopathogen has a complex taxonomic structure, which is constantly being revised due to the emergence of new molecular and biochemical diagnostic methods. Here for the first time, we describe the genetic and phenotypic diversity of 57 strains of Pseudomonas syringae isolated from affected soybeans, cereals, sunflowers, and other plants in the Russian Federation from 1950 to 2019. Genetic diversity was assessed by Multi Locus Sequence Analysis (MLSA) using fragments of the genes of glyceraldehyde-3-phosphate dehydrogenase (gapdh), the DNA-directed RNA polymerase subunit D (rpoD), gyrase (topoisomerase) B subunit (gyrB), and citrate synthase I (gltA). The synthesis of syringomycin and coronatine by bacteria was assessed by the reaction of susceptible yeast culture, seedlings of barley, tomato, and sunflower, and by presence of toxin genes confirmed by PCR test. The pathogenicity of the strains was confirmed on seedlings of dicotyledonous and monocotyledonous plants of peas, soybean, sunflowers, barley and wheat, as the most affected crops. The sensitivity of bacteria to 10 antibiotics of the main mechanisms of activity and two bactericidal commercial products was tested by standard disc method. The obtained results showed a high genetic homogeneity of the Russian population of P. syringae, which infects various agricultural crops, and an increase in the proportion of antibiotic-resistant strains over the years.

Comparative Genomics of Prunus-Associated Members of the Pseudomonas syringae Species Complex Reveals Traits Supporting Co-evolution and Host Adaptation

Members of the Pseudomonas syringae species complex cause symptoms that are ranging from leaf spots to cankers on a multitude of plant species, including some of the genus Prunus. To date, a total of two species of the P. syringae species complex and six different pathovars have been associated with diseases on Prunus spp., which were shown to belong to different phylogenetic units (phylogroups, PG) based on sequence similarity of housekeeping genes or whole genomes, suggesting that virulence to Prunus spp. may be the result of convergent pathoadaptation. In this study, a comparative genomics approach was used to determine genes significantly associated with strains isolated from Prunus spp. across a phylogeny of 97 strains belonging to the P. syringae species complex. Our study revealed the presence of a set of orthologous proteins which were significantly associated with strains isolated from Prunus spp. than in strains isolated from other hosts or from non-agricultural environments. Among them, the type III effector HopAY predicted to encode for a C58 cysteine protease was found to be highly associated with strains isolated from Prunus spp. and revealed patterns supporting co-evolution and host adaptation.

Two New Bacterial Pathogens of Peanut, Causing Early Seedling Decline Disease, Identified in the Texas Panhandle

Peanut (Arachis hypogaea L.) is cultivated in tropical and subtropical regions of the world as an important source of oil and protein. Until now, bacterial wilt, caused by Ralstonia solanacearum, was the only known bacterial disease of peanut. In 2020, widespread incidence of poor stand establishment was observed in multiple production fields planted to the Spanish-type peanut varieties in the Texas Panhandle. The observed symptoms included seed rot, pre- and postemergence damping-off, poor seedling vigor, poorly developed root systems with little or no nodule formation, and death. Subsequent diagnosis of symptomatic seedlings recovered two bacterial species identified by BLAST using 676- and 661-bp 16S rRNA fragments as a Ralstonia sp. and a Pantoea sp., respectively. To investigate a possible causative role of these bacteria in the observed peanut disease, the pathogenicity of the two isolates was evaluated under greenhouse conditions relying on Koch's postulates. Cell suspensions of the two bacteria, separately and in combination, were used to inoculate seeds of a Valencia-type peanut variety with no history of the disease and found to be pathogenic on the resultant seedling plants. Symptoms that developed on the inoculated plants were similar to the symptoms initially observed in the field, including seed rot, pre- and postemergence damping-off, poor seedling vigor, and root establishment. The two bacteria were also successfully recovered from inoculated and symptomatic plants, thus satisfying Koch's postulates. Given the early onset of symptom development on affected seeds and seedlings, a seedborne origin of the disease, described here as early-decline bacterial disease of peanut, was investigated in the same batches of peanut seeds that were planted, as well as seeds later harvested in some of the affected fields. Identical bacterial species, on the basis of 16S rRNA identity, were recovered from all of the seeds evaluated indicating that the bacteria are both seedborne and seed-transmissible. Multilocus sequence analysis involving six genes (dnaK, fumC, gyrB, murG, trpB, and tuf) showed that these new strains are most closely related to R. pickettii and Pantoea dispersa, but also phylogenetically distinct. The two bacteria were designated Ralstonia sp. strain B265 and Pantoea sp. strain B270. Losses from the disease in affected fields in 2020 averaged 50% (US$1.12 million) from a total of nine production fields. Findings from this study provide evidence for two new bacterial pathogens of peanuts capable of infecting Spanish and Valencia peanut varieties.

Genome Sequence and Antifungal Activity of Two Niche-Sharing Pseudomonas protegens Related Strains Isolated from Hydroponics

This work reports the comparison of the genome sequence and the ability to inhibit fungal growth of two Pseudomonas protegens related strains that were isolated from the same hydroponic culture of lamb's lettuce. The two strains were very similar in their core genome but one strain, Pf4, contained three gene clusters for the production of secondary metabolites, i.e., pyoluteorin (plt), pyrrolnitrin (prn), and rhizoxin (rzx), that were missing in the other strain, Pf11. The difference between the two strains was not due to simple insertion events, but to a relatively complex differentiation focused on the accessory genomes. In dual culture assays, both strains inhibited nearly all tested fungal strains, yet Pf4 exerted a significantly stronger fungal growth inhibition than Pf11. In addition to the differences in the secondary metabolite production associated genes abundance, the genome of Pf4 was more stable, smaller in size and with a lower number of transposons. The preservation of a dynamic equilibrium within natural populations of different strains comprised in the same species but differing in their secondary metabolite repertoire and in their genome stability may be functional to the adaptation to environmental changes.

Inference of Convergent Gene Acquisition Among Pseudomonas syringae Strains Isolated From Watermelon, Cantaloupe, and Squash

Pseudomonas syringae sensu stricto (phylogroup 2; referred to as P. syringae) consists of an environmentally ubiquitous bacterial population associated with diseases of numerous plant species. Recent studies using multilocus sequence analysis have indicated the clonal expansion of several P. syringae lineages, located in phylogroups 2a and 2b, in association with outbreaks of bacterial spot disease of watermelon, cantaloupe, and squash in the United States. To investigate the evolutionary processes that led to the emergence of these epidemic lineages, we sequenced the genomes of six P. syringae strains that were isolated from cucurbits grown in the United States, Europe, and China over a period of more than a decade, as well as eight strains that were isolated from watermelon and squash grown in six different Florida counties during the 2013 and 2014 seasons. These data were subjected to comparative analyses along with 42 previously sequenced genomes of P. syringae stains collected from diverse plant species and environments available from GenBank. Maximum likelihood reconstruction of the P. syringae core genome revealed the presence of a hybrid phylogenetic group, comprised of cucurbit strains collected in Florida, Italy, Serbia, and France, which emerged through genome-wide homologous recombination between phylogroups 2a and 2b. Functional analysis of the recombinant core genome showed that pathways involved in the ATP-dependent transport and metabolism of amino acids, bacterial motility, and secretion systems were enriched for recombination. A survey of described virulence factors indicated the convergent acquisition of several accessory type 3 secreted effectors (T3SEs) among phylogenetically distinct lineages through integrative and conjugative element and plasmid loci. Finally, pathogenicity assays on watermelon and squash showed qualitative differences in virulence between strains of the same clonal lineage, which correlated with T3SEs acquired through various mechanisms of horizontal gene transfer (HGT). This study provides novel insights into the interplay of homologous recombination and HGT toward pathogen emergence and highlights the dynamic nature of P. syringae sensu lato genomes.

Some strains that have converged to infect Prunus spp. trees are members of distinct Pseudomonas syringae genomospecies and ecotypes as revealed by in silico genomic comparison

A complementary taxonomic and population genetic study was performed to delineate genetically and ecologically distinct species within the Pseudomonas syringae complex by assessing 16 strains including pathovar strains that have converged to infect Prunus spp. trees, and two outgroups. Both average nucleotide identity and genome-to-genome distance comparison methods revealed the occurrence of distinct genomospecies, namely 1, 2, 3 and 8 (sensu Gardan et al.), with the latter two being closely related. Strains classified as P. s. pv. morsprunorum clustered into two distinct genomospecies, namely 2 and 8. Both the AdaptML and hierarchical Bayesian analysis of population structure methods highlighted the presence of three ecotypes, and the taxonomically related genomospecies 3 and 8 strains were members of the same ecotype. The distribution of pathogenic and virulence-associated genetic traits among Pseudomonas strains did not reveal any distinct type III secretion system effector or phytotoxin distribution pattern that characterized single genomospecies and strains that infect Prunus spp. The complete WHOP (Woody HOst and Pseudomonas spp.) genomic region and the entire β-ketoadipate gene cluster, including the catBCA operon, were found only in the members of genomospecies 2 and in the two P. s. pv. morsprunorum strains of genomospecies 8. A reduced gene flow between the three ecotypes suggested that point mutations played a larger role during the evolution of the strains than recombination. Our data support the idea that Prunus trees can be infected by different strains of distinct Pseudomonas genomospecies/ecotypes through diverse mechanisms of host colonization and infection. Such strains may represent particular lineages that emerged from environments other than that of the infected plant upon acquiring genetic traits that gave them the ability to cause plant diseases. The complementary assessment of bacterial strains using both taxonomic approaches and methods that reveal ecologically homogeneous populations has proven useful in confirming the cohesion of bacterial clusters.

Characterization of the pathogenicity of strains of Pseudomonas syringae towards cherry and plum

Bacterial canker is a major disease of Prunus avium (cherry), Prunus domestica (plum) and other stone fruits. It is caused by pathovars within the Pseudomonas syringae species complex including P. syringae pv. morsprunorum (Psm) race 1 (R1), Psm race 2 (R2) and P. syringae pv. syringae (Pss). Psm R1 and Psm R2 were originally designated as the same pathovar; however, phylogenetic analysis revealed them to be distantly related, falling into phylogroups 3 and 1, respectively. This study characterized the pathogenicity of 18 newly genome-sequenced P. syringae strains on cherry and plum, in the field and laboratory. The field experiment confirmed that the cherry cultivar Merton Glory exhibited a broad resistance to all clades. Psm R1 contained strains with differential specificity on cherry and plum. The ability of tractable laboratory-based assays to reproduce assessments on whole trees was examined. Good correlations were achieved with assays using cut shoots or leaves, although only the cut shoot assay was able to reliably discriminate cultivar differences seen in the field. Measuring bacterial multiplication in detached leaves differentiated pathogens from nonpathogens and was therefore suitable for routine testing. In cherry leaves, symptom appearance discriminated Psm races from nonpathogens, which triggered a hypersensitive reaction. Pathogenic strains of Pss rapidly induced disease lesions in all tissues and exhibited a more necrotrophic lifestyle than hemibiotrophic Psm. This in-depth study of pathogenic interactions, identification of host resistance and optimization of laboratory assays provides a framework for future genetic dissection of host-pathogen interactions in the canker disease.

Genomic Structural Variations Affecting Virulence During Clonal Expansion of Pseudomonas syringae pv. actinidiae Biovar 3 in Europe

Pseudomonas syringae pv. actinidiae (Psa) biovar 3 caused pandemic bacterial canker of Actinidia chinensis and Actinidia deliciosa since 2008. In Europe, the disease spread rapidly in the kiwifruit cultivation areas from a single introduction. In this study, we investigated the genomic diversity of Psa biovar 3 strains during the primary clonal expansion in Europe using single molecule real-time (SMRT), Illumina and Sanger sequencing technologies. We recorded evidences of frequent mobilization and loss of transposon Tn6212, large chromosome inversions, and ectopic integration of IS sequences (remarkably ISPsy31, ISPsy36, and ISPsy37). While no phenotype change associated with Tn6212 mobilization could be detected, strains CRAFRU 12.29 and CRAFRU 12.50 did not elicit the hypersensitivity response (HR) on tobacco and eggplant leaves and were limited in their growth in kiwifruit leaves due to insertion of ISPsy31 and ISPsy36 in the hrpS and hrpR genes, respectively, interrupting the hrp cluster. Both strains had been isolated from symptomatic plants, suggesting coexistence of variant strains with reduced virulence together with virulent strains in mixed populations. The structural differences caused by rearrangements of self-genetic elements within European and New Zealand strains were comparable in number and type to those occurring among the European strains, in contrast with the significant difference in terms of nucleotide polymorphisms. We hypothesize a relaxation, during clonal expansion, of the selection limiting the accumulation of deleterious mutations associated with genome structural variation due to transposition of mobile elements. This consideration may be relevant when evaluating strategies to be adopted for epidemics management.

Clarification of Taxonomic Status within the Pseudomonas syringae Species Group Based on a Phylogenomic Analysis

The Pseudomonas syringae phylogenetic group comprises 15 recognized bacterial species and more than 60 pathovars. The classification and identification of strains is relevant for practical reasons but also for understanding the epidemiology and ecology of this group of plant pathogenic bacteria. Genome-based taxonomic analyses have been introduced recently to clarify the taxonomy of the whole genus. A set of 139 draft and complete genome sequences of strains belonging to all species of the P. syringae group available in public databases were analyzed, together with the genomes of closely related species used as outgroups. Comparative genomics based on the genome sequences of the species type strains in the group allowed the delineation of phylogenomic species and demonstrated that a high proportion of strains included in the study are misclassified. Furthermore, representatives of at least 7 putative novel species were detected. It was also confirmed that P. ficuserectae, P. meliae, and P. savastanoi are later synonyms of P. amygdali and that “P. coronafaciens” should be revived as a nomenspecies.

Pseudomonas caspiana sp. nov., a citrus pathogen in the Pseudomonas syringae phylogenetic group

In a screening by multilocus sequence analysis of Pseudomonas strains isolated from diverse origins, 4 phylogenetically closely related strains (FBF58, FBF102^T, FBF103, and FBF122) formed a well-defined cluster in the Pseudomonas syringae phylogenetic group. The strains were isolated from citrus orchards in northern Iran with disease symptoms in the leaves and stems and its pathogenicity against citrus plants was demonstrated. The whole genome of the type strain of the proposed new species (FBF102^T=CECT 9164^T=CCUG 69273^T) was sequenced and characterized. Comparative genomics with the 14 known Pseudomonas species type strains of the P. syringae phylogenetic group demonstrated that this strain belonged to a new genomic species, different from the species described thus far. Genome analysis detected genes predicted to be involved in pathogenesis, such as an atypical type 3 secretion system and two type 6 secretion systems, together with effectors and virulence factors. A polyphasic taxonomic characterization demonstrated that the 4 plant pathogenic strains represented a new species, for which the name Pseudomonas caspiana sp. nov. is proposed.

WHOP, a Genomic Region Associated With Woody Hosts in the Pseudomonas syringae Complex Contributes to the Virulence and Fitness of Pseudomonas savastanoi pv. savastanoi in Olive Plants

Bacteria from the Pseudomonas syringae complex belonging to phylogroups 1 and 3 (PG1 and PG3, respectively) isolated from woody hosts share a genomic region herein referred to as WHOP (from woody host and Pseudomonas spp.), which is absent in strains infecting herbaceous organs. In this work, we show that this region is also encoded in P. syringae pv. actinidifoliorum (PG1) and six additional members of PG3, namely, Pseudomonas savastanoi pv. retacarpa, three P. syringae pathovars, Pseudomonas meliae, and Pseudomonas amygdali. Partial conservation of the WHOP occurs in only a few PG2 strains. In P. savastanoi pv. savastanoi NCPPB 3335, the WHOP region is organized into four operons and three independently transcribed genes. While the antABC and catBCA operons mediate the catabolism of anthranilate and catechol, respectively, the ipoABC operon confers oxygenase activity to aromatic compounds. The deletion of antABC, catBCA, or ipoABC in NCPPB 3335 caused reduced virulence in woody olive plants without affecting knot formation in nonwoody plants; catBCA, dhoAB, and PSA3335_3206 (encoding a putative aerotaxis receptor) were also required for the full fitness of this strain exclusively in woody olive plants. Overall, this study sheds light on the evolution and adaptation of bacteria from the P. syringae complex to woody hosts and highlights the enzymatic activities encoded within the WHOP region that are essential for this process.

Comprehensive analysis of draft genomes of two closely related pseudomonas syringae phylogroup 2b strains infecting mono- and dicotyledon host plants

Background

In recent years, the damage caused by bacterial pathogens to major crops has been increasing worldwide. Pseudomonas syringae is a widespread bacterial species that infects almost all major crops. Different P. syringae strains use a wide range of biochemical mechanisms, including phytotoxins and effectors of the type III and type IV secretion systems, which determine the specific nature of the pathogen virulence.

Results

Strains 1845 (isolated from dicots) and 2507 (isolated from monocots) were selected for sequencing because they specialize on different groups of plants. We compared virulence factors in these and other available genomes of phylogroup 2 to find genes responsible for the specialization of bacteria. We showed that strain 1845 belongs to the clonal group that has been infecting monocots in Russia and USA for a long time (at least 50 years). Strain 1845 has relatively recently changed its host plant to dicots.

Conclusions

The results obtained by comparing the strain 1845 genome with the genomes of bacteria infecting monocots can help to identify the genes that define specific nature of the virulence of P. syringae strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3358-y) contains supplementary material, which is available to authorized users.

Comparative genomics reveals genes significantly associated with woody hosts in the plant pathogen Pseudomonas syringae

The diversification of lineages within Pseudomonas syringae has involved a number of adaptive shifts from herbaceous hosts onto various species of tree, resulting in the emergence of highly destructive diseases such as bacterial canker of kiwi and bleeding canker of horse chestnut. This diversification has involved a high level of gene gain and loss, and these processes are likely to play major roles in the adaptation of individual lineages onto their host plants. In order to better understand the evolution of P. syringae onto woody plants, we have generated de novo genome sequences for 26 strains from the P. syringae species complex that are pathogenic on a range of woody species, and have looked for statistically significant associations between gene presence and host type (i.e. woody or herbaceous) across a phylogeny of 64 strains. We have found evidence for a common set of genes associated with strains that are able to colonize woody plants, suggesting that divergent lineages have acquired similarities in genome composition that may form the genetic basis of their adaptation to woody hosts. We also describe in detail the gain, loss and rearrangement of specific loci that may be functionally important in facilitating this adaptive shift. Overall, our analyses allow for a greater understanding of how gene gain and loss may contribute to adaptation in P. syringae.

Draft Genome Sequence of a Pseudomonas sp. Strain Carrying blaIMP-25 and blaVIM-2 Carbapenemase Genes from Hospital Sewage

Pseudomonas strain WCHP16 recovered from hospital sewage in West China Hospital, Chengdu, China was found to carry two carbapenemase genes bla_IMP-25 and bla_VIM-2. Here, we report its 5.7-Mb draft genome sequence, comprising 141 contigs and an average 59.53% G+C content. The genome contained 5,504 coding sequences and 67 tRNA genes.

Genome-wide comparison and taxonomic relatedness of multiple Xylella fastidiosa strains reveal the occurrence of three subspecies and a new Xylella species

A total of 21 Xylella fastidiosa strains were assessed by comparing their genomes to infer their taxonomic relationships. The whole-genome-based average nucleotide identity and tetranucleotide frequency correlation coefficient analyses were performed. In addition, a consensus tree based on comparisons of 956 core gene families, and a genome-wide phylogenetic tree and a Neighbor-net network were constructed with 820,088 nucleotides (i.e., approximately 30-33 % of the entire X. fastidiosa genome). All approaches revealed the occurrence of three well-demarcated genetic clusters that represent X. fastidiosa subspecies fastidiosa, multiplex and pauca, with the latter appeared to diverge. We suggest that the proposed but never formally described subspecies 'sandyi' and 'morus' are instead members of the subspecies fastidiosa. These analyses support the view that the Xylella strain isolated from Pyrus pyrifolia in Taiwan is likely to be a new species. A widely used multilocus sequence typing analysis yielded conflicting results.

Genomic and proteomic evidence supporting the division of the plant pathogen Ralstonia solanacearum into three species

Background

The increased availability of genome sequences has advanced the development of genomic distance methods to describe bacterial diversity. Results of these fast-evolving methods are highly correlated with those of the historically standard DNA-DNA hybridization technique. However, these genomic-based methods can be done more rapidly and less expensively and are less prone to technical and human error. They are thus a technically accessible replacement for species delineation. Here, we use several genomic comparison methods, supported by our own proteomic analyses and metabolic characterization as well as previously published DNA-DNA hybridization analyses, to differentiate members of the Ralstonia solanacearum species complex into three species. This pathogen group consists of diverse and widespread strains that cause bacterial wilt disease on many different plants.

Results

We used three different methods to compare the complete genomes of 29 strains from the R. solanacearum species complex. In parallel we profiled the proteomes of 73 strains using Matrix-Assisted Laser Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF-MS). Proteomic profiles together with genomic sequence comparisons consistently and comprehensively described the diversity of the R. solanacearum species complex. In addition, genome-driven functional phenotypic assays excitingly supported an old hypothesis (Hayward et al. (J Appl Bacteriol 69:269–80, 1990)), that closely related members of the R. solanacearum could be identified through a simple assay of anaerobic nitrate metabolism. This assay allowed us to clearly and easily differentiate phylotype II and IV strains from phylotype I and III strains. Further, genomic dissection of the pathway distinguished between proposed subspecies within the current phylotype IV. The assay revealed large scale differences in energy production within the R. solanacearum species complex, indicating coarse evolutionary distance and further supporting a repartitioning of this group into separate species.

Conclusions

Together, the results of these studies support the proposed division of the R. solanacearum species complex into three species, consistent with recent literature, and demonstrate the utility of proteomic and genomic approaches to delineate bacterial species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2413-z) contains supplementary material, which is available to authorized users.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.

Diversity of Pseudomonas Genomes, Including Populus-Associated Isolates, as Revealed by Comparative Genome Analysis

The Pseudomonas genus contains a metabolically versatile group of organisms that are known to occupy numerous ecological niches, including the rhizosphere and endosphere of many plants. Their diversity influences the phylogenetic diversity and heterogeneity of these communities. On the basis of average amino acid identity, comparative genome analysis of >1,000 Pseudomonas genomes, including 21 Pseudomonas strains isolated from the roots of native Populus deltoides (eastern cottonwood) trees resulted in consistent and robust genomic clusters with phylogenetic homogeneity. All Pseudomonas aeruginosa genomes clustered together, and these were clearly distinct from other Pseudomonas species groups on the basis of pangenome and core genome analyses. In contrast, the genomes of Pseudomonas fluorescens were organized into 20 distinct genomic clusters, representing enormous diversity and heterogeneity. Most of our 21 Populus-associated isolates formed three distinct subgroups within the major P. fluorescens group, supported by pathway profile analysis, while two isolates were more closely related to Pseudomonas chlororaphis and Pseudomonas putida. Genes specific to Populus-associated subgroups were identified. Genes specific to subgroup 1 include several sensory systems that act in two-component signal transduction, a TonB-dependent receptor, and a phosphorelay sensor. Genes specific to subgroup 2 contain hypothetical genes, and genes specific to subgroup 3 were annotated with hydrolase activity. This study justifies the need to sequence multiple isolates, especially from P. fluorescens, which displays the most genetic variation, in order to study functional capabilities from a pangenomic perspective. This information will prove useful when choosing Pseudomonas strains for use to promote growth and increase disease resistance in plants.

Comparative Genomic Analyses of Multiple Pseudomonas Strains Infecting Corylus avellana Trees Reveal the Occurrence of Two Genetic Clusters with Both Common and Distinctive Virulence and Fitness Traits

The European hazelnut (Corylus avellana) is threatened in Europe by several pseudomonads which cause symptoms ranging from twig dieback to tree death. A comparison of the draft genomes of nine Pseudomonas strains isolated from symptomatic C. avellana trees was performed to identify common and distinctive genomic traits. The thorough assessment of genetic relationships among the strains revealed two clearly distinct clusters: P. avellanae and P. syringae. The latter including the pathovars avellanae, coryli and syringae. Between these two clusters, no recombination event was found. A genomic island of approximately 20 kb, containing the hrp/hrc type III secretion system gene cluster, was found to be present without any genomic difference in all nine pseudomonads. The type III secretion system effector repertoires were remarkably different in the two groups, with P. avellanae showing a higher number of effectors. Homologue genes of the antimetabolite mangotoxin and ice nucleation activity clusters were found solely in all P. syringae pathovar strains, whereas the siderophore yersiniabactin was only present in P. avellanae. All nine strains have genes coding for pectic enzymes and sucrose metabolism. By contrast, they do not have genes coding for indolacetic acid and anti-insect toxin. Collectively, this study reveals that genomically different Pseudomonas can converge on the same host plant by suppressing the host defence mechanisms with the use of different virulence weapons. The integration into their genomes of a horizontally acquired genomic island could play a fundamental role in their evolution, perhaps giving them the ability to exploit new ecological niches.

The bacterial pangenome as a new tool for analysing pathogenic bacteria

The bacterial pangenome was introduced in 2005 and, in recent years, has been the subject of many studies. Thanks to progress in next-generation sequencing methods, the pangenome can be divided into two parts, the core (common to the studied strains) and the accessory genome, offering a large panel of uses. In this review, we have presented the analysis methods, the pangenome composition and its application as a study of lifestyle. We have also shown that the pangenome may be used as a new tool for redefining the pathogenic species. We applied this to the Escherichia coli and Shigella species, which have been a subject of controversy regarding their taxonomic and pathogenic position.

•

Pangenome is a new way of studying pathogenic bacteria.

•

Pangenome can be used as a taxonomic tool.

•

This review describes pangenome in the world of pathogenic bacteria.

Insights from 20 years of bacterial genome sequencing

Since the first two complete bacterial genome sequences were published in 1995, the science of bacteria has dramatically changed. Using third-generation DNA sequencing, it is possible to completely sequence a bacterial genome in a few hours and identify some types of methylation sites along the genome as well. Sequencing of bacterial genome sequences is now a standard procedure, and the information from tens of thousands of bacterial genomes has had a major impact on our views of the bacterial world. In this review, we explore a series of questions to highlight some insights that comparative genomics has produced. To date, there are genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. However, the distribution is quite skewed towards a few phyla that contain model organisms. But the breadth is continuing to improve, with projects dedicated to filling in less characterized taxonomic groups. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system provides bacteria with immunity against viruses, which outnumber bacteria by tenfold. How fast can we go? Second-generation sequencing has produced a large number of draft genomes (close to 90 % of bacterial genomes in GenBank are currently not complete); third-generation sequencing can potentially produce a finished genome in a few hours, and at the same time provide methlylation sites along the entire chromosome. The diversity of bacterial communities is extensive as is evident from the genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. Genome sequencing can help in classifying an organism, and in the case where multiple genomes of the same species are available, it is possible to calculate the pan- and core genomes; comparison of more than 2000 Escherichia coli genomes finds an E. coli core genome of about 3100 gene families and a total of about 89,000 different gene families. Why do we care about bacterial genome sequencing? There are many practical applications, such as genome-scale metabolic modeling, biosurveillance, bioforensics, and infectious disease epidemiology. In the near future, high-throughput sequencing of patient metagenomic samples could revolutionize medicine in terms of speed and accuracy of finding pathogens and knowing how to treat them.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.

Definition of Plant-Pathogenic Pseudomonas Genomospecies of the Pseudomonas syringae Complex Through Multiple Comparative Approaches

A total of 34 phytopathogenic strain genomes belonging to the Pseudomonas syringae species complex and related species, including many pathotype strains, were assessed using average nucleotide identity (ANI) analysis. Their taxonomic relationships were consistently confirmed by the tetranucleotide frequency correlation coefficient (TETRA) values, multilocus sequence typing analysis (MLSA) performed with seven housekeeping genes, using both maximum likelihood and Bayesian methods, and split consensus network analyses. The ANI, MLSA, and split consensus analyses provided consistent and identical results. We confirmed the occurrence of the well-demarcated genomospecies inferred sensu Gardan et al. using DNA-DNA hybridization and ribotyping analyses. However, some P. syringae strains of the pathovars morsprunorum and lachrymans were placed in different genomospecies in our analyses. Genomospecies 1, 2, 4, 6, and 9 resulted well demarcated, whereas strains of genomospecies 3 and 8 had ANI values between 95 and 96% in some cases, confirming that this threshold reveals very closely related species that might represent cases of splitting entities or the convergence of different species to the same ecological niche. This study confirms the robustness of the combination of genomic and phylogenetic approaches in revealing taxonomic relationships among closely related bacterial strains and provides the basis for a further reliable demarcation of the phytopathogenic Pseudomonas species. Within each species, the pathovars might represent distinct ecological units. The possibility of performing extensive and standardized host range and phenotypic tests with many strains of different pathovars can assist phytobacteriologists for better determining the boundaries of these ecological units.

Bioinformatic genome comparisons for taxonomic and phylogenetic assignments using Aeromonas as a test case

Prokaryotic taxonomy is the underpinning of microbiology, as it provides a framework for the proper identification and naming of organisms. The “gold standard” of bacterial species delineation is the overall genome similarity determined by DNA-DNA hybridization (DDH), a technically rigorous yet sometimes variable method that may produce inconsistent results. Improvements in next-generation sequencing have resulted in an upsurge of bacterial genome sequences and bioinformatic tools that compare genomic data, such as average nucleotide identity (ANI), correlation of tetranucleotide frequencies, and the genome-to-genome distance calculator, or in silico DDH (isDDH). Here, we evaluate ANI and isDDH in combination with phylogenetic studies using Aeromonas, a taxonomically challenging genus with many described species and several strains that were reassigned to different species as a test case. We generated improved, high-quality draft genome sequences for 33 Aeromonas strains and combined them with 23 publicly available genomes. ANI and isDDH distances were determined and compared to phylogenies from multilocus sequence analysis of housekeeping genes, ribosomal proteins, and expanded core genes. The expanded core phylogenetic analysis suggested relationships between distant Aeromonas clades that were inconsistent with studies using fewer genes. ANI values of ≥96% and isDDH values of ≥70% consistently grouped genomes originating from strains of the same species together. Our study confirmed known misidentifications, validated the recent revisions in the nomenclature, and revealed that a number of genomes deposited in GenBank are misnamed. In addition, two strains were identified that may represent novel Aeromonas species.

Improvements in DNA sequencing technologies have resulted in the ability to generate large numbers of high-quality draft genomes and led to a dramatic increase in the number of publically available genomes. This has allowed researchers to characterize microorganisms using genome data. Advantages of genome sequence-based classification include data and computing programs that can be readily shared, facilitating the standardization of taxonomic methodology and resolving conflicting identifications by providing greater uniformity in an overall analysis. Using Aeromonas as a test case, we compared and validated different approaches. Based on our analyses, we recommend cutoff values for distance measures for identifying species. Accurate species classification is critical not only to obviate the perpetuation of errors in public databases but also to ensure the validity of inferences made on the relationships among species within a genus and proper identification in clinical and veterinary diagnostic laboratories.

A user's guide to a data base of the diversity of Pseudomonas syringae and its application to classifying strains in this phylogenetic complex

The Pseudomonas syringae complex is composed of numerous genetic lineages of strains from both agricultural and environmental habitats including habitats closely linked to the water cycle. The new insights from the discovery of this bacterial species in habitats outside of agricultural contexts per se have led to the revelation of a wide diversity of strains in this complex beyond what was known from agricultural contexts. Here, through Multi Locus Sequence Typing (MLST) of 216 strains, we identified 23 clades within 13 phylogroups among which the seven previously described P. syringae phylogroups were included. The phylogeny of the core genome of 29 strains representing nine phylogroups was similar to the phylogeny obtained with MLST thereby confirming the robustness of MLST-phylogroups. We show that phenotypic traits rarely provide a satisfactory means for classification of strains even if some combinations are highly probable in some phylogroups. We demonstrate that the citrate synthase (cts) housekeeping gene can accurately predict the phylogenetic affiliation for more than 97% of strains tested. We propose a list of cts sequences to be used as a simple tool for quickly and precisely classifying new strains. Finally, our analysis leads to predictions about the diversity of P. syringae that is yet to be discovered. We present here an expandable framework mainly based on cts genetic analysis into which more diversity can be integrated.