Population genomics reveals additive and replacing horizontal gene transfers in the emerging pathogen Dickeya solani

Investigating Additive and Replacing Horizontal Gene Transfers Using Phylogenies and Whole Genomes

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.

Clonality and Diversity in the Soft Rot Dickeya solani Phytopathogen

Bacterial diversity analyses often suffer from a bias due to sampling only from a limited number of hosts or narrow geographic locations. This was the case for the phytopathogenic species Dickeya solani, whose members were mainly isolated from a few hosts-potato and ornamentals-and from the same geographical area-Europe and Israel, which are connected by seed trade. Most D. solani members were clonal with the notable exception of the potato isolate RNS05.1.2A and two related strains that are clearly distinct from other D. solani genomes. To investigate if D. solani genomic diversity might be broadened by analysis of strains isolated from other environments, we analysed new strains isolated from ornamentals and from river water as well as strain CFBP 5647 isolated from tomato in the Caribbean island Guadeloupe. While water strains were clonal to RNS05.1.2A, the Caribbean tomato strain formed a third clade. The genomes of the three clades are highly syntenic; they shared almost 3900 protein families, and clade-specific genes were mainly included in genomic islands of extrachromosomal origin. Our study thus revealed both broader D. solani diversity with the characterisation of a third clade isolated in Latin America and a very high genomic conservation between clade members.

Phyloproteomic study by MALDI-TOF MS in view of intraspecies variation in a significant homogenous phytopathogen Dickeya solani

Dickeya solani is an economically significant pectinolytic phytopathogen belonging to the Pectobacteriaceae family, which causes soft rot and blackleg diseases. Despite its notable impact on global potato production, there are no effective methods to control this pest. Here, we undertook a phyloproteomic study on 20 D. solani strains, of various origin and year of isolation, with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) supported by an in-depth characterization of the strains in terms of the virulence-associated phenotype. In spite of high homogeneity in this species, we herein revealed for the first time intraspecies variation in the MALDI-TOF MS protein profiles among the studied D. solani isolates. Finally, representative mass spectra for the four delineated clades are presented. A majority of the analysed D. solani strains showed high virulence potential, while two strains stood out in their growth dynamics, virulence factors production and ability to macerate plant tissue. Nonetheless, the metabolic profiles of D. solani strains turned out to be uniform, except for gelatinase activity. Given that all D. solani isolates distinctly grouped from the other Dickeya species in the MALDI-TOF MS analysis, there is strong evidence supporting the potential routine use of this method for fast and reliable to-species identification of D. solani isolates of environmental origin.

Insight into biodiversity of the recently rearranged genus Dickeya

The genus Dickeya includes plant pathogenic bacteria attacking a wide range of crops and ornamentals as well as a few environmental isolates from water. Defined on the basis of six species in 2005, this genus now includes 12 recognized species. Despite the description of several new species in recent years, the diversity of the genus Dickeya is not yet fully explored. Many strains have been analyzed for species causing diseases on economically important crops, such as for the potato pathogens D. dianthicola and D. solani. In contrast, only a few strains have been characterized for species of environmental origin or isolated from plants in understudied countries. To gain insights in the Dickeya diversity, recent extensive analyzes were performed on environmental isolates and poorly characterized strains from old collections. Phylogenetic and phenotypic analyzes led to the reclassification of D. paradisiaca (containing strains from tropical or subtropical regions) in the new genus, Musicola, the identification of three water species D. aquatica, D. lacustris and D. undicola, the description of a new species D. poaceaphila including Australian strains isolated from grasses, and the characterization of the new species D. oryzae and D. parazeae, resulting from the subdivision of the species D. zeae. Traits distinguishing each new species were identified from genomic and phenotypic comparisons. The high heterogeneity observed in some species, notably for D. zeae, indicates that additional species still need to be defined. The objective of this study was to clarify the present taxonomy of the genus Dickeya and to reassign the correct species to several Dickeya strains isolated before the current classification.

Dissimilar gene repertoires of Dickeya solani involved in the colonization of lesions and roots of Solanum tuberosum

Dickeya and Pectobacterium species are necrotrophic pathogens that macerate stems (blackleg disease) and tubers (soft rot disease) of Solanum tuberosum. They proliferate by exploiting plant cell remains. They also colonize roots, even if no symptoms are observed. The genes involved in pre-symptomatic root colonization are poorly understood. Here, transposon-sequencing (Tn-seq) analysis of Dickeya solani living in macerated tissues revealed 126 genes important for competitive colonization of tuber lesions and 207 for stem lesions, including 96 genes common to both conditions. Common genes included acr genes involved in the detoxification of plant defense phytoalexins and kduD, kduI, eda (=kdgA), gudD, garK, garL, and garR genes involved in the assimilation of pectin and galactarate. In root colonization, Tn-seq highlighted 83 genes, all different from those in stem and tuber lesion conditions. They encode the exploitation of organic and mineral nutrients (dpp, ddp, dctA, and pst) including glucuronate (kdgK and yeiQ) and synthesis of metabolites: cellulose (celY and bcs), aryl polyene (ape), and oocydin (ooc). We constructed in-frame deletion mutants of bcsA, ddpA, apeH, and pstA genes. All mutants were virulent in stem infection assays, but they were impaired in the competitive colonization of roots. In addition, the ΔpstA mutant was impaired in its capacity to colonize progeny tubers. Overall, this work distinguished two metabolic networks supporting either an oligotrophic lifestyle on roots or a copiotrophic lifestyle in lesions. This work revealed novel traits and pathways important for understanding how the D. solani pathogen efficiently survives on roots, persists in the environment, and colonizes progeny tubers.

Solanimycin: Biosynthesis and Distribution of a New Antifungal Antibiotic Regulated by Two Quorum-Sensing Systems

The increasing emergence of drug-resistant fungal infections has necessitated a search for new compounds capable of combating fungal pathogens of plants, animals, and humans. Microorganisms represent the main source of antibiotics with applicability in agriculture and in the clinic, but many aspects of their metabolic potential remain to be explored. This report describes the discovery and characterization of a new antifungal compound, solanimycin, produced by a hybrid polyketide/nonribosomal peptide (PKS/NRPS) system in Dickeya solani, the enterobacterial pathogen of potato. Solanimycin was active against a broad range of plant-pathogenic fungi of global economic concern and the human pathogen Candida albicans. The genomic cluster responsible for solanimycin production was defined and analyzed to identify the corresponding biosynthetic proteins, which include four multimodular PKS/NRPS proteins and several tailoring enzymes. Antifungal production in D. solani was enhanced in response to experimental conditions found in infected potato tubers and high-density fungal cultures. Solanimycin biosynthesis was cell density dependent in D. solani and was controlled by both the ExpIR acyl-homoserine lactone and Vfm quorum-sensing systems of the bacterial phytopathogen. The expression of the solanimycin cluster was also regulated at the post-transcriptional level, with the regulator RsmA playing a major role. The solanimycin biosynthetic cluster was conserved across phylogenetically distant bacterial genera, and multiple pieces of evidence support that the corresponding gene clusters were acquired by horizontal gene transfer. Given its potent broad-range antifungal properties, this study suggests that solanimycin and related molecules may have potential utility for agricultural and clinical exploitation. IMPORTANCE Fungal infections represent a major clinical, agricultural, and food security threat worldwide, which is accentuated due to the difficult treatment of these infections. Microorganisms represent a prolific source of antibiotics, and current data support that this enormous biosynthetic potential has been scarcely explored. To improve the performance in the discovery of novel antimicrobials, there is a need to diversify the isolation niches for new antibiotic-producing microorganisms as well as to scrutinize novel phylogenetic positions. With the identification of the antifungal antibiotic solanimycin in a broad diversity of phytopathogenic Dickeya spp., we provide further support for the potential of plant-associated bacteria for the biosynthesis of novel antimicrobials. The complex regulatory networks involved in solanimycin production reflect the high metabolic cost of bacterial secondary metabolism. This metabolic regulatory control makes many antibiotics cryptic under standard laboratory conditions, and mimicking environmental conditions, as shown here, is a strategy to activate cryptic antibiotic clusters.

Patterns of Genomic Variations in the Plant Pathogen Dickeya solani

The plant pathogen Dickeya solani causes soft rot and blackleg diseases in several crops including Solanum tuberosum. Unveiling the patterns of its diversity contributes to understanding the emergence and virulence of this pathogen in potato agro-systems. In this study, we analyzed the genome of several D. solani strains exhibiting an atypically high number of genetic variations. Variant calling and phylogenomics support the evidence that the strains RNS10-105-1A, A623S-20A-17 and RNS05.1.2A belong to a divergent sub-group of D. solani for which we proposed RNS05.1.2A as a reference strain. In addition, we showed that the variations (1253 to 1278 snp/indels) in strains RNS13-30-1A, RNS13-31-1A and RNS13-48-1A were caused by a horizontal gene transfer event from a donor belonging to the D. solani RNS05.1.2A subgroup. The overall results highlight the patterns driving the diversification in D. solani species. This work contributes to understanding patterns and causes of diversity in the emerging pathogen D. solani.

BuscoPhylo: a webserver for Busco-based phylogenomic analysis for non-specialists

Here we present the BuscoPhylo tool that enables both students and established scientists to easily perform Busco-based phylogenomic analysis starting from a set of genomes sequences. BuscoPhylo is an efficient and user-friendly web server freely accessible at https://buscophylo.inra.org.ma/ . The source code, along with documentation, is freely available under an MIT license at https://github.com/alaesahbou/BuscoPhylo .

Influence of glucose on swarming and quorum sensing of Dickeya solani

Dickeya solani is a pathogen most frequently responsible for infecting potato plants in Europe. As in the case of most plant pathogens, its ability to colonize and invade the host depends on chemotaxis and motility. The coordinated movement of Dickeya over solid surfaces is governed by a quorum sensing mechanism. In D. solani motility is regulated by ExpI-ExpR proteins, homologous to luxI-luxR system from Vibrio fisheri, in which N-acyl-homoserine lactones (AHLs) serve as signaling molecules. Moreover, in many Gram-negative bacteria motility is coupled with central metabolism via carbon catabolite repression. This enables them to reach more nutrient-efficient niches. The aim of this study was to analyze the swarming motility of D. solani depending on the volume of the medium in the cultivation plate and glucose content. We show that the ability of this bacterium to move is strictly dependent on both these factors. Moreover, we analyze the production of AHLs and show that the quorum sensing mechanism in D. solani is also influenced by the availability of glucose in the medium and that the distribution of these signaling molecules are different depending on the volume of the medium in the plate.

Pathogenicity and Relative Abundance of Dickeya and Pectobacterium Species in Switzerland: An Epidemiological Dichotomy

Pectobacterium and Dickeya species are the causal agents of blackleg and soft rot diseases in potatoes. The main pathogenic species identified so far on potatoes are Dickeya dianthicola, Dickeya solani, Pectobacterium atrosepticum, Pectobacterium brasiliense, Pectobacterium carotovorum, and Pectobacterium parmentieri. Ten years ago, the most prevalent Soft Rot Pectobacteriaceae in Europe were the Dickeya species, P. atrosepticum and P. carotovorum, with some variations among countries. Since then, a drastic increase in the abundance of P. brasiliense has been observed in most European countries. This shift is difficult to explain without comparing the pathogenicity of all Dickeya and Pectobacterium species. The pathogenicity of all the above-mentioned bacterial species was assessed in field trials and in vitro tuber slice trials in Switzerland. Two isolates of each species were inoculated by soaking tubers of cv. Desiree in a suspension of 10⁵ CFU/mL, before planting in the field. For all trials, the Dickeya species were the most virulent ones, but long-term strain surveys performed in Switzerland indicate that P. brasiliense is currently the most frequent species detected. Our results show that the pathogenicity of the species is not the main factor explaining the high prevalence of P. brasiliense and P. parmentieri in the Swiss potato fields.

Early Emergence of Dickeya solani Revealed by Analysis of Dickeya Diversity of Potato Blackleg and Soft Rot Causing Pathogens in Switzerland

Blackleg and soft rot in potato caused by Pectobacterium and Dickeya enterobacteral genera are among the most destructive bacterial diseases in this crop worldwide. In Europe, over the last century, Pectobacterium spp. were the predominant causal agents of these diseases. As for Dickeya, before the large outbreak caused by D. solani in the 2000s, only D. dianthicola was isolated in Europe. The population dynamics of potato blackleg causing soft rot Pectobacteriaceae was, however, different in Switzerland as compared to that in other European countries with a high incidence (60 up to 90%) of Dickeya species (at the time called Erwinia chrysanthemi) already in the 1980s. To pinpoint what may underlie this Swiss peculiarity, we analysed the diversity present in the E. chrysanthemi Agroscope collection gathering potato isolates from 1985 to 2000s. Like elsewhere in Europe during this period, the majority of Swiss isolates belonged to D. dianthicola. However, we also identified a few isolates, such as D. chrysanthemi and D. oryzeae, two species that have not yet been reported in potatoes in Europe. Interestingly, this study allowed the characterisation of two "early" D. solani isolated in the 1990s. Genomic comparison between these early D. solani strains and strains isolated later during the large outbreak in the 2000s in Europe revealed only a few SNP and gene content differences, none of them affecting genes known to be important for virulence.

European Population of Pectobacterium punjabense: Genomic Diversity, Tuber Maceration Capacity and a Detection Tool for This Rarely Occurring Potato Pathogen

Enterobacteria belonging to the Pectobacterium and Dickeya genera are responsible for soft rot and blackleg diseases occurring in many crops around the world. Since 2016, the number of described species has more than doubled. However, some new species, such as Pectobacterium punjabense, are often poorly characterized, and little is known about their genomic and phenotypic variation. Here, we explored several European culture collections and identified seven strains of P. punjabense. All were collected from potato blackleg symptoms, sometimes from a long time ago, i.e., the IFB5596 strain isolated almost 25 years ago. We showed that this species remains rare, with less than 0.24% of P. punjabense strains identified among pectinolytic bacteria present in the surveyed collections. The analysis of the genomic diversity revealed the non-clonal character of P. punjabense species. Furthermore, the strains showed aggressiveness differences. Finally, a qPCR Taqman assay was developed for rapid and specific strain characterization and for use in diagnostic programs.

Diversity of limestone bacteriophages infecting Dickeya solani isolated in the Czech Republic

Seven novel tailed lytic viruses (Ds3CZ, Ds5CZ, Ds9CZ, Ds16CZ, Ds20CZ, Ds23CZ, Ds25CZ) infecting the bacterium Dickeya solani were isolated in the Czech Republic. Genomes of these viruses are dsDNA, 149,364 to 155,285 bp in length, and the genome arrangement is very similar to that of the type virus Dickeya virus LIMEstone 1. All but the Ds25CZ virus should be regarded as strains of a single species. Most of the sequence differences are due to the presence or absence of homing endonuclease (HE) genes, with 23 HEs found in Ds3CZ, Ds5CZ, and Ds20CZ, 22 in Ds9CZ, 19 in Ds16CZ, 18 in Ds25CZ, and 15 in Ds23CZ.

Pattern and causes of the establishment of the invasive bacterial potato pathogen Dickeya solani and of the maintenance of the resident pathogen D. dianthicola

Invasive pathogens can be a threat when they affect human health, food production or ecosystem services, by displacing resident species, and we need to understand the cause of their establishment. We studied the patterns and causes of the establishment of the pathogen Dickeya solani that recently invaded potato agrosystems in Europe by assessing its invasion dynamics and its competitive ability against the closely related resident D. dianthicola species. Epidemiological records over one decade in France revealed the establishment of D. solani and the maintenance of the resident D. dianthicola in potato fields exhibiting blackleg symptoms. Using experimentations, we showed that D. dianthicola caused a higher symptom incidence on aerial parts of potato plants than D. solani, while D. solani was more aggressive on tubers (i.e. with more severe symptoms). In co-infection assays, D. dianthicola outcompeted D. solani in aerial parts, while the two species co-existed in tubers. A comparison of 76 D. solani genomes (56 of which have been sequenced here) revealed balanced frequencies of two previously uncharacterized alleles, VfmB_Pro and VfmB_Ser , at the vfmB virulence gene. Experimental inoculations showed that the VfmB_Ser population was more aggressive on tubers, while the VfmB_Pro population outcompeted the VfmB_Ser population in stem lesions, suggesting an important role of the vfmB virulence gene in the ecology of the pathogens. This study thus brings novel insights allowing a better understanding of the pattern and causes of the D.solani invasion into potato production agrosystems, and the reasons why the endemic D. dianthicola nevertheless persisted.

Genome-Wide Identification of Dickeya solani Transcriptional Units Up-Regulated in Response to Plant Tissues From a Crop-Host Solanum tuberosum and a Weed-Host Solanum dulcamara

Dickeya solani is a Gram-negative bacterium able to cause disease symptoms on a variety of crop and ornamental plants worldwide. Weeds including Solanum dulcamara (bittersweet nightshade) growing near agricultural fields have been reported to support populations of soft rot bacteria in natural settings. However, little is known about the specific interaction of D. solani with such weed plants that may contribute to its success as an agricultural pathogen. The aim of this work was to assess the interaction of D. solani with its crop plant (Solanum tuberosum) and an alternative (S. dulcamara) host plant. From a collection of 10,000 Tn5 transposon mutants of D. solani IPO2222 carrying an inducible, promotorless gusA reporter gene, 210 were identified that exhibited plant tissue-dependent expression of the gene/operon into which the Tn5 insertion had occurred. Thirteen Tn5 mutants exhibiting the greatest plant tissue induction of such transcriptional units in S. tuberosum or S. dulcamara as measured by qRT-PCR were assessed for plant host colonization, virulence, and ability to macerate plant tissue, as well as phenotypes likely to contribute to the ecological fitness of D. solani, including growth rate, carbon and nitrogen source utilization, motility, chemotaxis toward plant extracts, biofilm formation, growth under anaerobic conditions and quorum sensing. These 13 transcriptional units encode proteins involved in bacterial interactions with plants, with functions linked to cell envelope structure, chemotaxis and carbon metabolism. The selected 13 genes/operons were differentially expressed in, and thus contributed preferentially to D. solani fitness in potato and/or S. dulcamara stem, leaf, and root tissues.

Comparative genomics and pangenome-oriented studies reveal high homogeneity of the agronomically relevant enterobacterial plant pathogen Dickeya solani

Background

Dickeya solani is an important plant pathogenic bacterium causing severe losses in European potato production. This species draws a lot of attention due to its remarkable virulence, great devastating potential and easier spread in contrast to other Dickeya spp. In view of a high need for extensive studies on economically important soft rot Pectobacteriaceae, we performed a comparative genomics analysis on D. solani strains to search for genetic foundations that would explain the differences in the observed virulence levels within the D. solani population.

Results

High quality assemblies of 8 de novo sequenced D. solani genomes have been obtained. Whole-sequence comparison, ANIb, ANIm, Tetra and pangenome-oriented analyses performed on these genomes and the sequences of 14 additional strains revealed an exceptionally high level of homogeneity among the studied genetic material of D. solani strains. With the use of 22 genomes, the pangenome of D. solani, comprising 84.7% core, 7.2% accessory and 8.1% unique genes, has been almost completely determined, suggesting the presence of a nearly closed pangenome structure. Attribution of the genes included in the D. solani pangenome fractions to functional COG categories showed that higher percentages of accessory and unique pangenome parts in contrast to the core section are encountered in phage/mobile elements- and transcription- associated groups with the genome of RNS 05.1.2A strain having the most significant impact. Also, the first D. solani large-scale genome-wide phylogeny computed on concatenated core gene alignments is herein reported.

Conclusions

The almost closed status of D. solani pangenome achieved in this work points to the fact that the unique gene pool of this species should no longer expand. Such a feature is characteristic of taxa whose representatives either occupy isolated ecological niches or lack efficient mechanisms for gene exchange and recombination, which seems rational concerning a strictly pathogenic species with clonal population structure. Finally, no obvious correlations between the geographical origin of D. solani strains and their phylogeny were found, which might reflect the specificity of the international seed potato market.

Diversity of Pectobacteriaceae Species in Potato Growing Regions in Northern Morocco

Dickeya and Pectobacterium pathogens are causative agents of several diseases that affect many crops worldwide. This work investigated the species diversity of these pathogens in Morocco, where Dickeya pathogens have only been isolated from potato fields recently. To this end, samplings were conducted in three major potato growing areas over a three-year period (2015-2017). Pathogens were characterized by sequence determination of both the gapA gene marker and genomes using Illumina and Oxford Nanopore technologies. We isolated 119 pathogens belonging to P. versatile (19%), P. carotovorum (3%), P. polaris (5%), P. brasiliense (56%) and D. dianthicola (17%). Their taxonomic assignation was confirmed by draft genome analyses of 10 representative strains of the collected species. D. dianthicola were isolated from a unique area where a wide species diversity of pectinolytic pathogens was observed. In tuber rotting assays, D. dianthicola isolates were more aggressive than Pectobacterium isolates. The complete genome sequence of D. dianthicola LAR.16.03.LID was obtained and compared with other D. dianthicola genomes from public databases. Overall, this study highlighted the ecological context from which some Dickeya and Pectobacterium species emerged in Morocco, and reported the first complete genome of a D. dianthicola strain isolated in Morocco that will be suitable for further epidemiological studies.

Species Diversity of Dickeya and Pectobacterium Causing Potato Blackleg Disease in Pakistan

Potato blackleg is caused by a diverse species of pectinolytic bacteria. In Pakistan, approximately 90% of the pathogens involved belong to Pectobacterium atrosepticum. Survey (2014 to 2017), sampling, and isolation from different potato growing areas of Punjab, Pakistan depicted an overall disease incidence of approximately 15%. Thirty-six pectinolytic strains confirmed through biochemical and pathogenicity testing were characterized via gapA gene to identify them at the species level. To further validate the identification, one strain from each species SS26 (P. atrosepticum), SS28 (Pectobacterium polaris), SS70 (Dickeya dianthicola), SS90 (Pectobacterium parmentieri), SS95 (Pectobacterium punjabense), and SS96 (Pectobacterium versatile) were selected for draft genome sequencing and multilocus sequence analysis of 13 housekeeping genes (fusA, rpoD, acnA, purA, gyrB, recA, mdh, mtlD, groEL, secY, glyA, gapA, and rplB). Phylogenetic analysis revealed considerable genetic diversity in the genus Pectobacterium. In silico DNA-DNA hybridization and average nucleotide identity values of the strains selected for genome sequencing were determined with other reference Pectobacterium and Dickeya strains. Moreover, all six representative strains were also phenotypically characterized on the basis of metabolism of different carbon sources. Overall, on the basis of genotypic and phenotypic characteristics, these 36 isolates were grouped into six species: P. atrosepticum, P. versatile, P. parmentieri, P. polaris, P. punjabense, and D. dianthicola.

Acetylome analysis of lysine acetylation in the plant pathogenic bacterium Brenneria nigrifluens

Protein lysine acetylation, a dynamic and reversible posttranslational modification, plays a crucial role in several cellular processes, including cell cycle regulation, metabolism, enzymatic activities, and protein interactions. Brenneria nigrifluens is a pathogen of walnut trees with shallow bark canker and can cause serious disease in walnut trees. Until now, a little has been known about the roles of lysine acetylation in plant pathogenic bacteria. In the present study, the lysine acetylome of B. nigrifluens was determined by high‐resolution LC‐MS/MS analysis. In total, we identified 1,866 lysine acetylation sites distributed in 737 acetylated proteins. Bioinformatics results indicated that acetylated proteins participate in many different biological functions in B. nigrifluens. Four conserved motifs, namely, LK^ac, K^ac*F, I*K^ac, and L*K^ac, were identified in this bacterium. Protein interaction network analysis indicated that all kinds of interactions are modulated by protein lysine acetylation. Overall, 12 acetylated proteins were related to the virulence of B. nigrifluens.

The phytopathogenic nature of Dickeya aquatica 174/2 and the dynamic early evolution of Dickeya pathogenicity

Dickeya is a genus of phytopathogenic enterobacterales causing soft rot in a variety of plants (e.g. potato, chicory, maize). Among the species affiliated to this genus, Dickeya aquatica, described in 2014, remained particularly mysterious because it had no known host. Furthermore, while D. aquatica was proposed to represent a deep-branching species among Dickeya genus, its precise phylogenetic position remained elusive. Here, we report the complete genome sequence of the D. aquatica type strain 174/2. We demonstrate the affinity of D. aquatica strain 174/2 for acidic fruits such as tomato and cucumber and show that exposure of this bacterium to acidic pH induces twitching motility. An in-depth phylogenomic analysis of all available Dickeya proteomes pinpoints D. aquatica as the second deepest branching lineage within this genus and reclassifies two lineages that likely correspond to new genomospecies (gs.): Dickeya gs. poaceaephila (Dickeya sp NCPPB 569) and Dickeya gs. undicola (Dickeya sp 2B12), together with a new putative genus, tentatively named Prodigiosinella. Finally, from comparative analyses of Dickeya proteomes, we infer the complex evolutionary history of this genus, paving the way to study the adaptive patterns and processes of Dickeya to different environmental niches and hosts. In particular, we hypothesize that the lack of xylanases and xylose degradation pathways in D. aquatica could reflect adaptation to aquatic charophyte hosts which, in contrast to land plants, do not contain xyloglucans.

Genomic characterisation of the new Dickeya fangzhongdai species regrouping plant pathogens and environmental isolates

BACKGROUND

The Dickeya genus is part of the Pectobacteriaceae family that is included in the newly described enterobacterales order. It comprises a group of aggressive soft rot pathogens with wide geographic distribution and host range. Among them, the new Dickeya fangzhongdai species groups causative agents of maceration-associated diseases that impact a wide variety of crops and ornamentals. It affects mainly monocot plants, but D. fangzhongdai strains have also been isolated from pear trees and water sources. Here, we analysed which genetic novelty exists in this new species, what are the D. fangzhongdai-specific traits and what is the intra-specific diversity.

RESULTS

The genomes of eight D. fangzhongdai strains isolated from diverse environments were compared to 31 genomes of strains belonging to other Dickeya species. The D. fangzhongdai core genome regroups approximately 3500 common genes, including most genes that encode virulence factors and regulators characterised in the D. dadantii 3937 model strain. Only 38 genes are present in D. fangzhongdai and absent in all other Dickeyas. One of them encodes a pectate lyase of the PL10 family of polysaccharide lyases that is found only in a few bacteria from the plant environment, soil or human gut. Other D. fangzhongdai-specific genes with a known or predicted function are involved in regulation or metabolism. The intra-species diversity analysis revealed that seven of the studied D. fangzhongdai strains were grouped into two distinct clades. Each clade possesses a pool of 100-150 genes that are shared by the clade members, but absent from the other D. fangzhongdai strains and several of these genes are clustered into genomic regions. At the strain level, diversity resides mainly in the arsenal of T5SS- and T6SS-related toxin-antitoxin systems and in secondary metabolite biogenesis pathways.

CONCLUSION

This study identified the genome-specific traits of the new D. fangzhongdai species and highlighted the intra-species diversity of this species. This diversity encompasses secondary metabolites biosynthetic pathways and toxins or the repertoire of genes of extrachromosomal origin. We however didn't find any relationship between gene content and phenotypic differences or sharing of environmental habitats.

Pectobacterium punjabense sp. nov., isolated from blackleg symptoms of potato plants in Pakistan

Pectobacterium isolates SS95^T, SS54 and SS56 were collected from a potato field in the Chiniot district in the plains of the Punjab province, Pakistan. Sequencing of the gapA barcode revealed that these strains belong to a novel phylogenetic group separated from P.ectobacterium wasabiae and Pectobacterium parmentieri species. Furthermore, multilocus sequence analyses of 13 housekeeping genes (fusA, rpoD, acnA, purA, gyrB, recA, mdh, mtlD, groEL, secY, glyA, gapA and rplB) clearly distinguished the type strain, SS95^T, from its closest relatives, i.e. P. parmentieri RNS 08-42-1A^T and P. wasabiae CFBP3304^T, as well as from all the other known Pectobacteriumspecies. In silico DNA-DNA hybridization (<44.1 %) and average nucleotide identity (<90.75 %) values of strain SS95^T compared with other Pectobacterium type strains supported the delineation of a new species. Genomic and phenotypic comparisons permitted the identification of additional traits that distinguished the Pakistani isolates from all other known Pectobacterium type strains. The name Pectobacterium punjabense sp. nov. is proposed for this taxon with the type strain SS95^T (=CFBP 8604^T=LMG 30622^T).

Comparison of Highly and Weakly Virulent Dickeya solani Strains, With a View on the Pangenome and Panregulon of This Species

Bacteria belonging to the genera Dickeya and Pectobacterium are responsible for significant economic losses in a wide variety of crops and ornamentals. During last years, increasing losses in potato production have been attributed to the appearance of Dickeya solani. The D. solani strains investigated so far share genetic homogeneity, although different virulence levels were observed among strains of various origins. The purpose of this study was to investigate the genetic traits possibly related to the diverse virulence levels by means of comparative genomics. First, we developed a new genome assembly pipeline which allowed us to complete the D. solani genomes. Four de novo sequenced and ten publicly available genomes were used to identify the structure of the D. solani pangenome, in which 74.8 and 25.2% of genes were grouped into the core and dispensable genome, respectively. For D. solani panregulon analysis, we performed a binding site prediction for four transcription factors, namely CRP, KdgR, PecS and Fur, to detect the regulons of these virulence regulators. Most of the D. solani potential virulence factors were predicted to belong to the accessory regulons of CRP, KdgR, and PecS. Thus, some differences in gene expression could exist between D. solani strains. The comparison between a highly and a low virulent strain, IFB0099 and IFB0223, respectively, disclosed only small differences between their genomes but significant differences in the production of virulence factors like pectinases, cellulases and proteases, and in their mobility. The D. solani strains also diverge in the number and size of prophages present in their genomes. Another relevant difference is the disruption of the adhesin gene fhaB2 in the highly virulent strain. Strain IFB0223, which has a complete adhesin gene, is less mobile and less aggressive than IFB0099. This suggests that in this case, mobility rather than adherence is needed in order to trigger disease symptoms. This study highlights the utility of comparative genomics in predicting D. solani traits involved in the aggressiveness of this emerging plant pathogen.

Interplay of classic Exp and specific Vfm quorum sensing systems on the phenotypic features of Dickeya solani strains exhibiting different virulence levels

Bacteria from the genus Dickeya cause severe symptoms on numerous economically important plants. Dickeya solani is the Dickeya species most frequently found on infected potato plants in Europe. D. solani strains from different countries show high genetic homogeneity, but significant differences in their virulence level. Dickeya species possess two quorum sensing (QS) mechanisms: the Exp system based on classic N-acyl-homoserine lactone (AHL) signals and a specific system depending on the production and perception of a molecule of unknown structure, Virulence Factor Modulating (VFM). To study the interplay between these two QS systems, five D. solani strains exhibiting different virulence levels were selected. Mutants were constructed by inactivating genes coding for each QS system. Double mutants were obtained by simultaneous inactivation of genes coding for both QS systems. Most of the D. solani mutants showed an attenuation of chicory maceration and a decreased production of plant cell wall-degrading enzymes (PCWDEs) and motility, but to different degrees depending on the strain. The VFM-QS system seems to regulate virulence in both D. solani and Dickeya dadantii, but the AHL-QS system has greater effects in D. solani than in D. dadantii. The inactivation of both QS systems in D. solani did not reveal any additive effect on the tested features. The inactivation of vfm genes generally has a more dominant effect relative to that of exp genes. Thus, VFM- and AHL-QS systems do not work in synergy to modulate the production of diverse virulence factors and the ability to macerate plant tissue.

Complete Chromosome and Plasmid Sequences of Two Plant Pathogens, Dickeya solani Strains D s0432-1 and PPO 9019

Dickeya solani species are emerging bacterial pathogens of Solanum tuberosum. Here, we announce the complete genome sequences of two strains, Dickeya solani D s0432-1 and PPO 9019. Strain PPO 9019 represents the first described member of the genus Dickeya with an extrachromosomal genetic element.

The uniform structure of O-polysaccharides isolated from Dickeya solani strains of different origin

O-polysaccharides were isolated from lipopolysaccharides obtained from four different strains of plant pathogenic bacteria belonging to the species Dickeya solani: two of them were isolated in Poland (IFB0099 and IFB0158), the third in Germany (IFB0223) and the last one, D. solani Type Strain IPO2222, originated from the Netherlands. In addition, the O-polysaccharide of a closely related species D. dadantii strain 3937 was isolated. The purified polysaccharides of the five strains were analyzed using NMR spectroscopy and chemical methods. Sugar and methylation analyses, including absolute configuration assignment, together with NMR data revealed that all O-polysaccharides tested are homopolymers of 6-deoxy-d-altrose (d-6dAlt) the following structure: →2)-β-d-6dAltp-(1→.

Complete genome anatomy of the emerging potato pathogen Dickeya solani type strain IPO 2222T

Several species of the genus Dickeya provoke soft rot and blackleg diseases on a wide range of plants and crops. Dickeya solani has been identified as the causative agent of diseases outbreaks on potato culture in Europe for the last decade. Here, we report the complete genome of the D. solani IPO 2222^T. Using PacBio and Illumina technologies, a unique circular chromosome of 4,919,833 bp was assembled. The G + C content reaches 56% and the genomic sequence contains 4,059 predicted proteins. The ANI values calculated for D. solani IPO 2222^T vs. other available D. solani genomes was over 99.9% indicating a high genetic homogeneity within D. solani species.