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

Comparative genomic and transcriptome analyses of two Pectobacterium brasiliense strains revealed distinct virulence determinants and phenotypic features

Potato soft rot caused by Pectobacterium spp. are devastating diseases of potato which cause severe economic losses worldwide. Pectobacterium brasiliense is considered as one of the most virulent species. However, the virulence mechanisms and pathogenicity factors of this strain have not been fully elucidated. Here, through pathogenicity screening, we identified two Pectobacterium brasiliense isolates, SM and DQ, with distinct pathogenicity levels. SM exhibits higher virulence compared to DQ in inducing aerial stem rot, blackleg and tuber soft rot. Our genomic and transcriptomic analyses revealed that SM encodes strain specific genes with regard to plant cell wall degradation and express higher level of genes associated with bacterial motility and secretion systems. Our plate assays verified higher pectinase, cellulase, and protease activities, as well as fast swimming and swarming motility in SM. Importantly, a unique endoglucanase S specific to SM was identified. Expression of this cellulase in DQ greatly enhances its virulence compared to wild type strain. Our study sheds light on possible determinants causing different pathogenicity of Pectobacterium brasiliense species with close evolutionary distance and provides new insight into the direction of genome evolution in response to host variation and environmental stimuli.

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.

Genome characterization of Dickeya solani bacteriophage W2B

We have successfully characterized the complete genome sequence of the lytic Dickeya solani bacteriophage W2B, isolated from the Bunus Sewage Treatment Plant. The lytic phage from the Ningirsuvirus family has a 40,385-bp linear double-stranded DNA genome containing 51 coding sequences (CDSs).

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.

Pan-Genomics of Escherichia albertii for Antibiotic Resistance Profiling in Different Genome Fractions and Natural Product Mediated Intervention: In Silico Approach

Escherichia albertii is an emerging, enteric pathogen of significance. It was first isolated in 2003 from a pediatric diarrheal sample from Bangladesh. In this study, a comprehensive in silico strategy was followed to first list out antibiotic-resistant genes from core, accessory and unique genome fractions of 95 available genomes of E. albertii. Then, 56 drug targets were identified from the core essential genome. Finally, ZipA, an essential cell division protein that stabilizes the FtsZ protofilaments by cross-linking them and serves as a cytoplasmic membrane anchor for the Z ring, was selected for further downstream processing. It was computationally modeled using a threading approach, followed by virtual screening of two phytochemical libraries, Ayurvedic (n = 2103 compounds) and Traditional Chinese Medicine (n = 36,043 compounds). ADMET profiling, followed by PBPK modeling in the central body compartment, in a population of 250 non-diseased, 250 cirrhotic and 250 renally impaired people was attempted. ZINC85624912 from Chinese medicinal library showed the highest bioavailability and plasma retention. This is the first attempt to simulate the fate of natural products in the body through PBPK. Dynamics simulation of 20 ns for the top three compounds from both libraries was also performed to validate the stability of the compounds. The obtained information from the current study could aid wet-lab scientists to work on the scaffold of screened drug-like compounds from natural resources and could be useful in our quest for therapy against antibiotic-resistant E. albertii.

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.

Genomic Comparisons and Phenotypic Diversity of Dickeya zeae Strains Causing Bacterial Soft Rot of Banana in China

Bacterial soft rot of banana, caused by Dickeya zeae, is spreading rapidly in important banana growing areas in China and seriously threatens banana production. In this study, we sequenced the high-quality complete genomes of three typical banana strains, MS1 (size: 4,831,702-bp; genome coverages: 538x), MS_2014 (size: 4,740,000-bp; genome coverages: 586x) and MS_2018 (size: 4,787,201-bp; genome coverages: 583x), isolated in 2009, 2014, and 2018, respectively. To determine their genomic and phenotypic diversity with respect to their hosts of origin, they were compared with other D. zeae strains, including another representative banana strain MS2 from China. The sequenced strains were similar in utilization of carbon source and chemical substrates, and general genomic features of GC content, and tRNA and rRNA regions. They were also conserved in most virulence determinants, including gene-encoding secretion systems, plant cell wall degrading enzymes, and exopolysaccharides. We further explored their genomic diversity in the predicted genomic islands (GIs). These GIs were rich in integrases and transposases, where some genomic dissimilarity was observed in the flagellar gene cluster and several secondary metabolite gene clusters. Different constituents of core biosynthetic modules were found within the bacteriocin and aryl polyene (APE) pigment gene clusters, and the strains from banana showed different phenotypes with respect to antibiosis effects and colony pigmentation. Additionally, clustered regularly interspaced short palindromic repeat (CRISPR) and prophage elements, such as type I-F and III-A CRISPR arrays and an intact prophage of MS1-P5, contributed to bacterial diversity. Phylogenetic tree analysis and genome-genome nucleotide comparison confirmed the genomic divergence among the strains isolated from banana. Considering these characteristics, MS2 and MS_2014 probably diverged later than MS1, while MS_2018 was different and more similar to foreign strains isolated from other hosts in several characteristics. Strain MS_2018 caused severe symptoms on banana varieties previously considered moderately resistant or moderately susceptible, including varieties of Cavendish (Musa AAA) and Plantain (Musa ABB). Our study of genomic and phenotypic diversity raises public attention to the risk of spreading new pathogenic variants within banana growing regions and supports development of predictive strategies for disease control.

Genome-Wide Identification of Genes Important for Growth of Dickeya dadantii and Dickeya dianthicola in Potato (Solanum tuberosum) Tubers

Dickeya species are causal agents of soft rot diseases in many economically important crops, including soft rot disease of potato (Solanum tuberosum). Using random barcode transposon-site sequencing (RB-TnSeq), we generated genome-wide mutant fitness profiles of Dickeya dadantii 3937, Dickeya dianthicola ME23, and Dickeya dianthicola 67-19 isolates collected after passage through several in vitro and in vivo conditions. Though all three strains are pathogenic on potato, D. dadantii 3937 is a well-characterized model while D. dianthicola strains ME23 and 67-19 are recent isolates. Strain ME23 specifically was identified as a representative strain from a 2014 outbreak on potato. This study generated comparable gene fitness measurements across ecologically relevant conditions for both model and non-model strains. Tubers from the potato cultivars “Atlantic,” “Dark Red Norland,” and “Upstate Abundance” provided highly similar conditions for bacterial growth. Using the homolog detection software PyParanoid, we matched fitness values for orthologous genes in the three bacterial strains. Direct comparison of fitness among the strains highlighted shared and variable traits important for growth. Bacterial growth in minimal medium required many metabolic traits that were also essential for competitive growth in planta, such as amino acid, carbohydrate, and nucleotide biosynthesis. Growth in tubers specifically required the pectin degradation gene kduD. Disruption in three putative DNA-binding proteins had strain-specific effects on competitive fitness in tubers. Though the Soft Rot Pectobacteriaceae can cause disease with little host specificity, it remains to be seen the extent to which strain-level variation impacts virulence.

Pangenomic Analysis of Dickeya dianthicola Strains Related to the Outbreak of Blackleg and Soft Rot of Potato in the United States

Dickeya dianthicola has caused an outbreak of blackleg and soft rot of potato in the eastern half of the United States since 2015. To investigate genetic diversity of the pathogen, a comparative analysis was conducted on genomes of D. dianthicola strains. Whole genomes of 16 strains from the United States outbreak were assembled and compared with 16 previously sequenced genomes of D. dianthicola isolated from potato or carnation. Among the 32 strains, eight distinct clades were distinguished based on phylogenomic analysis. The outbreak strains were grouped into three clades, with the majority of the strains in clade I. Clade I strains were unique and homogeneous, suggesting a recent incursion of this strain into potato production from alternative hosts or environmental sources. The pangenome of the 32 strains contained 6,693 genes, 3,377 of which were core genes. By screening primary protein subunits associated with virulence from all U.S. strains, we found that many virulence-related gene clusters, such as plant cell wall degrading enzyme genes, flagellar and chemotaxis related genes, two-component regulatory genes, and type I/II/III secretion system genes, were highly conserved but that type IV and type VI secretion system genes varied. The clade I strains encoded two clusters of type IV secretion systems, whereas the clade II and III strains encoded only one cluster. Clade I and II strains encoded one more VgrG/PAAR spike protein than did clade III. Thus, we predicted that the presence of additional virulence-related genes may have enabled the unique clade I strain to become predominant in the U.S. outbreak.

Tailed Lytic Bacteriophages of Soft Rot Pectobacteriaceae

The study of the ecological and evolutionary traits of Soft Rot Pectobacteriaceae (SRP) comprising genera Pectobacterium and Dickeya often involves bacterial viruses (bacteriophages). Bacteriophages are considered to be a prospective tool for the ecologically safe and highly specific protection of plants and harvests from bacterial diseases. Information concerning bacteriophages has been growing rapidly in recent years, and this has included new genomics-based principles of taxonomic distribution. In this review, we summarise the data on phages infecting Pectobacterium and Dickeya that are available in publications and genomic databases. The analysis highlights not only major genomic properties that assign phages to taxonomic families and genera, but also the features that make them potentially suitable for phage control applications. Specifically, there is a discussion of the molecular mechanisms of receptor recognition by the phages and problems concerning the evolution of phage-resistant mutants.

Species of Dickeya and Pectobacterium Isolated during an Outbreak of Blackleg and Soft Rot of Potato in Northeastern and North Central United States

An outbreak of bacterial soft rot and blackleg of potato has occurred since 2014 with the epicenter being in the northeastern region of the United States. Multiple species of Pectobacterium and Dickeya are causal agents, resulting in losses to commercial and seed potato production over the past decade in the Northeastern and North Central United States. To clarify the pathogen present at the outset of the epidemic in 2015 and 2016, a phylogenetic study was made of 121 pectolytic soft rot bacteria isolated from symptomatic potato; also included were 27 type strains of Dickeya and Pectobacterium species, and 47 historic reference strains. Phylogenetic trees constructed based on multilocus sequence alignments of concatenated dnaJ, dnaX and gyrB fragments revealed the epidemic isolates to cluster with type strains of D. chrysanthemi, D. dianthicola, D. dadantii, P. atrosepticum, P. brasiliense, P. carotovorum, P. parmentieri, P. polaris, P. punjabense, and P. versatile. Genetic diversity within D. dianthicola strains was low, with one sequence type (ST1) identified in 17 of 19 strains. Pectobacterium parmentieri was more diverse, with ten sequence types detected among 37 of the 2015–2016 strains. This study can aid in monitoring future shifts in potato soft rot pathogens within the U.S. and inform strategies for disease management.

Interplay between Coumarin Accumulation, Iron Deficiency and Plant Resistance to Dickeya spp

Coumarins belong to a group of secondary metabolites well known for their high biological activities including antibacterial and antifungal properties. Recently, an important role of coumarins in plant resistance to pathogens and their release into the rhizosphere upon pathogen infection was discovered. It is also well documented that coumarins play a crucial role in the Arabidopsis thaliana growth under Fe-limited conditions. However, the mechanisms underlying interplay between plant resistance, accumulation of coumarins and Fe status, remain largely unknown. In this work, we investigated the effect of both mentioned factors on the disease severity using the model system of Arabidopsis/Dickeya spp. molecular interactions. We evaluated the disease symptoms in Arabidopsis plants, wild-type Col-0 and its mutants defective in coumarin accumulation, grown in hydroponic cultures with contrasting Fe regimes and in soil mixes. Under all tested conditions, Arabidopsis plants inoculated with Dickeya solani IFB0099 strain developed more severe disease symptoms compared to lines inoculated with Dickeya dadantii 3937. We also showed that the expression of genes encoding plant stress markers were strongly affected by D. solani IFB0099 infection. Interestingly, the response of plants to D. dadantii 3937 infection was genotype-dependent in Fe-deficient hydroponic solution.

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.

Multivariate Optimization of the FLC-dc-APGD-Based Reaction-Discharge System for Continuous Production of a Plasma-Activated Liquid of Defined Physicochemical and Anti-Phytopathogenic Properties

To the present day, no efficient plant protection method against economically important bacterial phytopathogens from the Pectobacteriaceae family has been implemented into agricultural practice. In this view, we have performed a multivariate optimization of the operating parameters of the reaction-discharge system, employing direct current atmospheric pressure glow discharge, generated in contact with a flowing liquid cathode (FLC-dc-APGD), for the production of a plasma-activated liquid (PAL) of defined physicochemical and anti-phytopathogenic properties. As a result, the effect of the operating parameters on the conductivity of PAL acquired under these conditions was assessed. The revealed optimal operating conditions, under which the PAL of the highest conductivity was obtained, were as follows: flow rate of the solution equaled 2.0 mL min^-1, the discharge current was 30 mA, and the inorganic salt concentration (ammonium nitrate, NH₄NO₃) in the solution turned out to be 0.50% (m/w). The developed PAL exhibited bacteriostatic and bactericidal properties toward Dickeya solani IFB0099 and Pectobacterium atrosepticum IFB5103 strains, with minimal inhibitory and minimal bactericidal concentrations equaling 25%. After 24 h exposure to 25% PAL, 100% (1-2 × 10⁶) of D. solani and P. atrosepticum cells lost viability. We attributed the antibacterial properties of PAL to the presence of deeply penetrating, reactive oxygen and nitrogen species (RONS), which were, in this case, OH, O, O₃, H₂O₂, HO₂, NH, N₂, N₂⁺, NO₂^-, NO₃^-, and NH₄⁺. Putatively, the generated low-cost, eco-friendly, easy-to-store, and transport PAL, exhibiting the required antibacterial and physicochemical properties, may find numerous applications in the plant protection sector.

Comparative Genomics, from the Annotated Genome to Valuable Biological Information: A Case Study

High availability of fast, cheap, and high-throughput next generation sequencing techniques resulted in acquisition of numerous de novo sequenced and assembled bacterial genomes. It rapidly became clear that digging out useful biological information from such a huge amount of data presents a considerable challenge. In this chapter we share our experience with utilization of several handy open source comparative genomic tools. All of them were applied in the studies focused on revealing inter- and intraspecies variation in pectinolytic plant pathogenic bacteria classified to Dickeya solani and Pectobacterium parmentieri. As the described software performed well on the species within the Pectobacteriaceae family, it presumably may be readily utilized on some closely related taxa from the Enterobacteriaceae family. First of all, implementation of various annotation software is discussed and compared. Then, tools computing whole genome comparisons including generation of circular juxtapositions of multiple sequences, revealing the order of synteny blocks or calculation of ANI or Tetra values are presented. Besides, web servers intended either for functional annotation of the genes of interest or for detection of genomic islands, plasmids, prophages, CRISPR/Cas are described. Last but not least, utilization of the software designed for pangenome studies and the further downstream analyses is explained. The presented work not only summarizes broad possibilities assured by the comparative genomic approach but also provides a user-friendly guide that might be easily followed by nonbioinformaticians interested in undertaking similar studies.

PacBio-Based Protocol for Bacterial Genome Assembly

Acquisition of high-quality bacterial genomes is fundamental, while having in mind investigation of subtitle intraspecies variation in addition to development of sensitive species-specific tools for detection and identification of the pathogens. In this view, Pacific Biosciences technology seems highly tempting taking into consideration over 10,000 bp length of the generated reads. In this work, we describe a bacterial genome assembly pipeline based on open-source software that might be handled also by non-bioinformaticians interested in transformation of sequencing data into reliable biological information. With the use of this method, we successfully closed six Dickeya solani genomes, while the assembly process was run just on a slightly improved desktop computer.

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.

High genomic variability in the plant pathogenic bacterium Pectobacterium parmentieri deciphered from de novo assembled complete genomes

BACKGROUND

Pectobacterium parmentieri is a newly established species within the plant pathogenic family Pectobacteriaceae. Bacteria belonging to this species are causative agents of diseases in economically important crops (e.g. potato) in a wide range of different environmental conditions, encountered in Europe, North America, Africa, and New Zealand. Severe disease symptoms result from the activity of P. parmentieri virulence factors, such as plant cell wall degrading enzymes. Interestingly, we observe significant phenotypic differences among P. parmentieri isolates regarding virulence factors production and the abilities to macerate plants. To establish the possible genomic basis of these differences, we sequenced 12 genomes of P. parmentieri strains (10 isolated in Poland, 2 in Belgium) with the combined use of Illumina and PacBio approaches. De novo genome assembly was performed with the use of SPAdes software, while annotation was conducted by NCBI Prokaryotic Genome Annotation Pipeline.

RESULTS

The pan-genome study was performed on 15 genomes (12 de novo assembled and three reference strains: P. parmentieri CFBP 8475T, P. parmentieri SCC3193, P. parmentieri WPP163). The pan-genome includes 3706 core genes, a high number of accessory (1468) genes, and numerous unique (1847) genes. We identified the presence of well-known genes encoding virulence factors in the core genome fraction, but some of them were located in the dispensable genome. A significant fraction of horizontally transferred genes, virulence-related gene duplications, as well as different CRISPR arrays were found, which can explain the observed phenotypic differences. Finally, we found also, for the first time, the presence of a plasmid in one of the tested P. parmentieri strains isolated in Poland.

CONCLUSIONS

We can hypothesize that a large number of the genes in the dispensable genome and significant genomic variation among P. parmentieri strains could be the basis of the potential wide host range and widespread diffusion of P. parmentieri. The obtained data on the structure and gene content of P. parmentieri strains enabled us to speculate on the importance of high genomic plasticity for P. parmentieri adaptation to different environments.