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

T6SS and T4SS Redundantly Secrete Effectors to Govern the Virulence and Bacterial Competition in Pectobacterium PccS1

Previous studies revealed that the type VI secretion system (T6SS) has an essential role in bacterial competition and virulence in many gram-negative bacteria. However, the role of T6SS in virulence in Pectobacterium atrosepticum remains controversial. We examined a closely related strain, PccS1, and discovered that its T6SS comprises a single-copy cluster of 17 core genes with a higher identity to homologs from P. atrosepticum. Through extensive phenotypic and functional analyses of over 220 derivatives of PccS1, we found that three of the five VgrGs could be classified into group I VgrGs. These VgrGs interacted with corresponding DUF4123 domain proteins, which were secreted outside of the membrane and were dependent on either the T6SS or type IV secretion system (T4SS). This interaction directly governed virulence and competition. Meanwhile, supernatant proteomic analyses with strains defective in the T6SS and/or T4SS confirmed that effectors, such as FhaB, were secreted redundantly to control the virulence and suppress host callose deposition in the course of infection. Notably, this redundant secretion mechanism between the T6SS and T4SS is believed to be the first of its kind in bacteria.

Antibiotic Resistance in Plant Pathogenic Bacteria: Recent Data and Environmental Impact of Unchecked Use and the Potential of Biocontrol Agents as an Eco-Friendly Alternative

The economic impact of phytopathogenic bacteria on agriculture is staggering, costing billions of US dollars globally. Pseudomonas syringae is the top most phytopathogenic bacteria, having more than 60 pathovars, which cause bacteria speck in tomatoes, halo blight in beans, and so on. Although antibiotics or a combination of antibiotics are used to manage infectious diseases in plants, they are employed far less in agriculture compared to human and animal populations. Moreover, the majority of antibiotics used in plants are immediately washed away, leading to environmental damage to ecosystems and food chains. Due to the serious risk of antibiotic resistance (AR) and the potential for environmental contamination with antibiotic residues and resistance genes, the use of unchecked antibiotics against phytopathogenic bacteria is not advisable. Despite the significant concern regarding AR in the world today, there are inadequate and outdated data on the AR of phytopathogenic bacteria. This review presents recent AR data on plant pathogenic bacteria (PPB), along with their environmental impact. In light of these findings, we suggest the use of biocontrol agents as a sustainable, eco-friendly, and effective alternative to controlling phytopathogenic bacteria.

Impact of Homologous Recombination on Core Genome Evolution and Host Adaptation of Pectobacterium parmentieri

Homologous recombination is a major force mechanism driving bacterial evolution, host adaptability, and acquisition of novel virulence traits. Pectobacterium parmentieri is a plant bacterial pathogen distributed worldwide, primarily affecting potatoes, by causing soft rot and blackleg diseases. The goal of this investigation was to understand the impact of homologous recombination on the genomic evolution of P. parmentieri. Analysis of P. parmentieri genomes using Roary revealed a dynamic pan-genome with 3,742 core genes and over 55% accessory genome variability. Bayesian population structure analysis identified 7 lineages, indicating species heterogeneity. ClonalFrameML analysis displayed 5,125 recombination events, with the lineage 4 exhibiting the highest events. fastGEAR analysis identified 486 ancestral and 941 recent recombination events ranging from 43 bp to 119 kb and 36 bp to 13.96 kb, respectively, suggesting ongoing adaptation. Notably, 11% (412 genes) of the core genome underwent recent recombination, with lineage 1 as the main donor. The prevalence of recent recombination (double compared to ancient) events implies continuous adaptation, possibly driven by global potato trade. Recombination events were found in genes involved in vital cellular processes (DNA replication, DNA repair, RNA processing, homeostasis, and metabolism), pathogenicity determinants (type secretion systems, cell-wall degrading enzymes, iron scavengers, lipopolysaccharides (LPS), flagellum, etc.), antimicrobial compounds (phenazine and colicin) and even CRISPR-Cas genes. Overall, these results emphasize the potential role of homologous recombination in P. parmentieri's evolutionary dynamics, influencing host colonization, pathogenicity, adaptive immunity, and ecological fitness.

Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum species complex

Background

Ralstonia solanacearum species complex (RSSC) strains are destructive plant pathogenic bacteria and the causative agents of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, their virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape RSSC virulence, the global diversity, distribution, and wider functional gene content of RSSC prophages are unknown. In this study, prophages were identified in a diverse collection of 192 RSSC draft genome assemblies originating from six continents.

Results

Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distributions were characterised using metadata on RSSC strain geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 313 intact prophages were identified, forming ten genetically distinct clusters. These included six prophage clusters with similarity to the Inoviridae, Myoviridae, and Siphoviridae phage families, and four uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distributions, being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained many auxiliary genes involved in metabolism and virulence of both phage and bacteria.

Conclusions

Our results show that while RSSC prophages are highly diverse globally, they make lineage-specific contributions to the RSSC accessory genome, which could have resulted from shared coevolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08909-7.

Genomic characterization of a new phage BUCT541 against Klebsiella pneumoniae K1-ST23 and efficacy assessment in mouse and Galleria mellonella larvae

Over the past decades, the spread of multi-drug-resistant Klebsiella pneumoniae (MDR-KP) is becoming a new threat and new effective therapies against this pathogen are needed. Bacteriophage (phage) therapy is considered to be a promising alternative treatment for MDR-KP infections compared with antibacterial drug usage. Here, we reported a new phage BUCT541 which can lyse MDR-KP ST23. The genome of BUCT541 is a double-stranded linear 46,100-bp long DNA molecule with 48% GC content through the Next generation sequencing (NGS) data. A total of 81 open reading frames and no virulence or antimicrobial resistance genes are annotated in the BUCT541 genome. BUCT541 was able to lyse 7 of the 30 tested MDR-KP according to the host range analysis. And the seven sensitive strains belonged to the K. pneumoniae K1-ST23. BUCT541 exhibited high thermal stability (4–70°C) and broad pH tolerance (pH 3-11) in the stability test. The in vivo results showed that BUCT541 (4 × 10⁵ plaque-forming units (PFU)/each) significantly increased the survival rate of K. pneumoniae infected Galleria mellonella from 5.3% to 83.3% within 48 h. Moreover, in the mouse lung infection model, high doses of BUCT541 (2 × 10⁷ PFU/each) cured 100% of BALB/c mice that were infected with K. pneumoniae. After 30 h of treatment with phage BUCT541 of the multiplicity of infection (MOI) = 10, the K. pneumoniae in the lungs of mice was lower than 10⁴ CFU/mL, compared to the control group 10⁹ CFU/mL. Together, these findings indicate that phage BUCT541 holds great promise as an alternative therapy with excellent stability and a wide lysis range for the treatment of MDR-KP ST23 infection.

Characterization and Comparative Genomics Analysis of a New Bacteriophage BUCT610 against Klebsiella pneumoniae and Efficacy Assessment in Galleria mellonella Larvae

The spread of multidrug-resistant Klebsiella pneumoniae (MDR-KP) has become an emerging threat as a result of the overuse of antibiotics. Bacteriophage (phage) therapy is considered to be a promising alternative treatment for MDR-KP infection compared with antibiotic therapy. In this research, a lytic phage BUCT610 was isolated from hospital sewage. The assembled genome of BUCT610 was 46,774 bp in length, with a GC content of 48%. A total of 83 open reading frames (ORFs) and no virulence or antimicrobial resistance genes were annotated in the BUCT610 genome. Comparative genomics and phylogenetic analyses showed that BUCT610 was most closely linked with the Vibrio phage pYD38-A and shared 69% homology. In addition, bacteriophage BUCT610 exhibited excellent thermal stability (4–75 °C) and broad pH tolerance (pH 3–12) in the stability test. In vivo investigation results showed that BUCT610 significantly increased the survival rate of Klebsiella pneumonia-infected Galleria mellonella larvae from 13.33% to 83.33% within 72 h. In conclusion, these findings indicate that phage BUCT610 holds great promise as an alternative agent with excellent stability for the treatment of MDR-KP infection.

Analysis of the Taxonomy, Synteny, and Virulence Factors for Soft Rot Pathogen Pectobacterium aroidearum in Amorphophallus konjac Using Comparative Genomics

Bacterial soft rot is a devastating disease for a wide range of crops, vegetables, and ornamental plants including konjac (Amorphophallus konjac). However, the pangenome and genomic plasticity of the konjac soft rot pathogens is little explored. In this study, we reported the complete genome sequences of 11 bacterial isolates that can cause typical soft rot symptoms in konjac by in vitro and in vivo pathogenicity tests. Based on in silico DNA-DNA hybridization, average nucleotide identity and phylogenomic analysis, all 11 isolates were determined to be Pectobacterium aroidearum. In addition, synteny analysis of these genomes revealed considerable chromosomal inversions, one of which is triggered by homologous recombination of ribose operon. Pangenome analysis and COG enrichment analysis showed that the pangenome of P. aroidearum is open and that accessory genes are enriched in replication, recombination, and repair. Variations in type IV secretion system and type VI secretion system were found, while plant cell wall degrading enzymes were conserved. Furthermore, sequence analyses also provided evidence for the presence of a type V secretion system in Pectobacterium. These findings advance our understanding of the pathogenicity determinants, genomic plasticity, and evolution of P. aroidearum.

Heterogenicity within the LPS Structure in Relation to the Chosen Genomic and Physiological Features of the Plant Pathogen Pectobacterium parmentieri

Pectobacterium parmentieri is a pectinolytic plant pathogenic bacterium causing high economic losses of cultivated plants. The highly devastating potential of this phytopathogen results from the efficient production of plant cell wall-degrading enzymes, i.e., pectinases, cellulases and proteases, in addition to the impact of accessory virulence factors such as motility, siderophores, biofilm and lipopolysaccharide (LPS). LPS belongs to pathogen-associated molecular patterns (PAMPs) and plays an important role in plant colonization and interaction with the defense systems of the host. Therefore, we decided to investigate the heterogeneity of O-polysaccharides (OPS) of LPS of different strains of P. parmentieri, in search of an association between the selected genomic and phenotypic features of the strains that share an identical structure of the OPS molecule. In the current study, OPS were isolated from the LPS of two P. parmentieri strains obtained either in Finland in the 1980s (SCC3193) or in Poland in 2013 (IFB5432). The purified polysaccharides were analyzed by utilizing 1D and 2D NMR spectroscopy (¹H, DQF-COSY, TOCSY, ROESY, HSQC, HSQC-TOCSY and HMBC) in addition to chemical methods. Sugar and methylation analyses of native polysaccharides, absolute configuration assignment of constituent monosaccharides and NMR spectroscopy data revealed that these two P. parmentieri strains isolated in different countries possess the same structure of OPS with a very rare residue of 5,7-diamino-3,5,7,9-tetradeoxy-l-glycero-l-manno-non-2-ulosonic acid (pseudaminic acid) substituted in the position C-8: →3)-β-d-Galf-(1→3)-α-d-Galp-(1→8)-β-Pse4Ac5Ac7Ac-(2→6)-α-d-Glcp-(1→6)-β-d-Glcp-(1→. The previous study indicated that three other P. parmentieri strains, namely IFB5427, IFB5408 and IFB5443, exhibit a different OPS molecule than SCC3193 and IFB5432. The conducted biodiversity-oriented assays revealed that the P. parmentieri IFB5427 and IFB5408 strains possessing the same OPS structure yielded the highest genome-wide similarity, according to average nucleotide identity analyses, in addition to the greatest ability to macerate chicory tissue among the studied P. parmentieri strains. The current research demonstrated a novel OPS structure, characteristic of at least two P. parmentieri strains (SCC3193 and IFB5432), and discussed the observed heterogenicity in the OPS of P. parmentieri in a broad genomic and phenotype-related context.

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.

Quantitative Real-Time PCR Assay for the Detection of Pectobacterium parmentieri, a Causal Agent of Potato Soft Rot

Pectobacterium parmentieri is a plant-pathogenic bacterium, recently attributed as a separate species, which infects potatoes, causing soft rot in tubers. The distribution of P. parmentieri seems to be global, although the bacterium tends to be accommodated to moderate climates. Fast and accurate detection systems for this pathogen are needed to study its biology and to identify latent infection in potatoes and other plant hosts. The current paper reports on the development of a specific and sensitive detection protocol based on a real-time PCR with a TaqMan probe for P. parmentieri, and its evaluation. In sensitivity assays, the detection threshold of this protocol was 10² cfu/mL on pure bacterial cultures and 10²–10³ cfu/mL on plant material. The specificity of the protocol was evaluated against P. parmentieri and more than 100 strains of potato-associated species of Pectobacterium and Dickeya. No cross-reaction with the non-target bacterial species, or loss of sensitivity, was observed. This specific and sensitive diagnostic tool may reveal a wider distribution and host range for P. parmentieri and will expand knowledge of the life cycle and environmental preferences of this pathogen.

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.

Analysis of the Taxonomy and Pathogenic Factors of Pectobacterium aroidearum L6 Using Whole-Genome Sequencing and Comparative Genomics

Soft rot pectobacteria are devastating plant pathogens with a global distribution and a broad host range. Pectobacterium aroidearum L6, previously isolated from leaves of Syngonium podophyllum, is a pectolytic bacterial pathogen that causes typical soft rot on S. podophyllum. There is a shortage for genome data of P. aroidearum, which seriously hinders research on classification and pathogenesis of Pectobacterium. We present here the complete genome sequence of P. aroidearum L6. The L6 strain carries a single 4,995,896-bp chromosome with 53.10% G + C content and harbors 4,306 predicted protein-coding genes. We estimated in silico DNA-DNA hybridization and average nucleotide identity values in combination with the whole-genome-based phylogeny from 19 Pectobacterium strains including P. aroidearum L6. The results showed that L6 and PC1 formed a population distinct from other populations of the Pectobacterium genus. Phylogenetic analysis based on 16S rRNA and genome sequences showed a close evolutionary relationship among Pectobacterium species. Overall, evolutionary analysis showed that L6 was in the same branch with PC1. In comparison with 18 Pectobacterium spp. reference pathogens, strain L6 had 2,712 gene families, among which 1,632 gene families were identified as orthologous to those strains, as well as 1 putative unique gene family. We discovered 478 genes, 10.4% of the total of predicted genes, that were potentially related to pathogenesis using the Virulence Factors of Pathogenic Bacteria database. A total of 25 genes were related to toxins, 35 encoded plant cell-wall degrading enzymes, and 122 were involved in secretion systems. This study provides a foundation for a better understanding of the genomic structure of P. aroidearum and particularly offers information for the discovery of potential pathogenic factors and the development of more effective strategies against this pathogen.

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.

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.

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.

Single gene enables plant pathogenic Pectobacterium to overcome host-specific chemical defence

Plants of the Brassicales order, including Arabidopsis and many common vegetables, produce toxic isothiocyanates to defend themselves against pathogens. Despite this defence, plant pathogenic microorganisms like Pectobacterium cause large yield losses in fields and during storage of crops. The bacterial gene saxA was previously found to encode isothiocyanate hydrolase that degrades isothiocyanates in vitro. Here we demonstrate in planta that saxA is a virulence factor that can overcome the chemical defence system of Brassicales plants. Analysis of the distribution of saxA genes in Pectobacterium suggests that saxA from three different phylogenetic origins are present within this genus. Deletion of saxA genes representing two of the most common classes from P. odoriferum and P. versatile resulted in significantly reduced virulence on Arabidopsis thaliana and Brassica oleracea. Furthermore, expressing saxA from a plasmid in a potato-specific P. parmentieri strain that does not naturally harbour this gene significantly increased the ability of the strain to macerate Arabidopsis. These findings suggest that a single gene may have a significant role in defining the host range of a plant pathogen.