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

Gene expression analyses on Dickeya solani strains of diverse virulence levels unveil important pathogenicity factors for this species

Dickeya solani causes soft rot and blackleg mainly on potato crops. High pathogenicity of this species results from efficient production of plant cell wall-degrading enzymes, especially pectate lyases, potent root colonization, and fast vascular movement. Despite genomic homogeneity, variations in virulence-related phenotypes suggest differences in the gene expression patterns between diverse strains. Therefore, the methylomes and transcriptomes of two strains (virulent IFB0099 and low virulent IFB0223), differing in tissue maceration capacity and virulence factors production, have been studied. Methylation analysis revealed no significant differences. However, the analysis of transcriptomes, studied under both non-induced and induced by polygalacturonic acid conditions (in order to mimic diverse stages of plant infection process), unveiled higher expression of pectate lyases (pelD, pelE, pelL), pectin esterase (pemA), proteases (prtE, prtD) and Vfm-associated quorum-sensing genes (vfmC, vfmD, vfmE) in IFB0099 strain compared to IFB0223. Additionally, the genes related to the secretion system II (T2SS) (gspJ, nipE) displayed higher induction of expression in IFB0099. Furthermore, IFB0099 showed more elevated expression of genes involved in flagella formation, which coincides with enhanced motility and pathogenicity of this strain compared to IFB0223. To sum up, differential expression analysis of genes important for the virulence of D. solani indicated candidate genes, which might be crucial for the pathogenicity of this species.

Development and validation of genome-informed and multigene-based qPCR and LAMP assays for accurate detection of Dickeya solani: a critical quarantine pathogen threatening the potato industry

Dickeya solani one of the most aggressive pectinolytic phytopathogens, causes blackleg disease in potatoes, resulting in significant economic losses and adversely impacting one of the world's most important food crops. The diagnostics methods are critical in monitoring the latent infection for international trade of potato seed tubers and in implementation of control strategies. Our study employed a whole-genome comparative approach, identifying unique target gene loci (LysR and TetR family of transcriptional regulators gene regions) and designing loop-mediated isothermal amplification (LAMP) and a multi-gene-based multiplex TaqMan qPCR assays for specific detection and differentiation of D. solani. Both methods underwent meticulous validation with extensive inclusivity and exclusivity panels, exhibiting 100% accuracy and no false positives or negatives. The LAMP method demonstrated the detection limit of 100 fg and 1 CFU per reaction using pure genomic DNA and crude bacterial cell lysate, respectively. The qPCR detection limit was 1 pg, 100 fg and 10 fg with quadplex, triplex, and singleplex, respectively. None of the assays were impacted by any inhibitory or competitive effects after adding host DNA (in qPCR) or crude lysate (in LAMP). The assays proved robust and reproducible in detecting the target pathogen in infected samples, with the LAMP assay being field-deployable due to its simplicity and rapid results acquisition within approximately 9 min. The reproducibility was confirmed by performing the assay in two independent laboratories. These assays offer a robust, rapid, and reliable solution for routine testing, with applications in phytosanitary inspection, disease diagnosis, and epidemiological studies.IMPORTANCEDickeya solani, one of the most aggressive soft rot causing bacteria and a quarantine pathogen, poses a severe threat to food security by causing substantial economic losses to the potato industry. Accurate and timely detection of this bacterium is vital for monitoring latent infections, particularly for international trade of potato seed tubers, and for implementing effective control strategies. In this research, we have developed LAMP and multi-gene-based multiplex TaqMan qPCR assays for specific detection of D. solani. These assays, characterized by their precision, rapidity, and robustness, are crucial for distinguishing D. solani from related species. Offering unparalleled sensitivity and specificity, these assays are indispensable for phytosanitary inspection and epidemiological monitoring, providing a powerful tool enabling management of this threatening pathogen.

Investigating the genomic and metabolic abilities of PGPR Pseudomonas fluorescens in promoting plant growth and fire blight management

Pseudomonas fluorescens is commonly found in diverse environments and is well known for its metabolic and antagonistic properties. Despite its remarkable attributes, its potential role in promoting plant growth remains unexplored. This study examines these traits across 14 strains residing in diverse rhizosphere environments through pangenome and comparative genome analysis, alongside molecular docking studies against Erwinia amylovora to combat fire blight. Whole genome analysis revealed circular chromosome (6.01-7.07 Mb) with GC content averaging 59.95-63.39%. Predicted genes included 16S rRNA and protein-coding genes ranging from 4435 to 6393 bp and 1527 to 1541 bp, respectively. Pangenome analysis unveiled an open pangenome, shedding light on genetic factors influencing plant growth promotion and biocontrol, including nitrogen fixation, phosphorus solubilization, siderophore production, stress tolerance, flagella biosynthesis, and induced systemic resistance. Furthermore, pyrrolnitrin, phenazine-1-carboxylic acid, pyoluteorin, lokisin, 2,4-diacetylpholoroglucinol and pseudomonic acid were identified. Molecular docking against key proteins of E. amylovora highlighted the high binding affinities of 2,4-diacetylphloroglucinol, pseudomonic acid, and lokisin. These findings underscore the multifaceted role of P. fluorescens in plant growth promotion and biocontrol, with key biomolecules showing promising applications in plant growth and defense against pathogens.

Deciphering the genomic character of the multidrug-resistant Staphylococcus aureus from Dhaka, Bangladesh

Staphylococcus aureus is one of the leading agents of nosocomial and community-acquired infections. In this study, we explored the genomic characterization of eight methicillin-resistant clinical isolates of S. aureus from Dhaka, Bangladesh. Notably, all strains were resistant to penicillin, cephalosporins, and monobactams, with partial susceptibility to meropenem and complete susceptibility to amikacin, vancomycin, and tigecycline antibiotics. The strains were found to have an average genome size of 2.73 Mbp and an average of 32.64% GC content. Multi-locus sequence typing analysis characterized the most predominant sequence type as ST361, which belongs to the clonal complex CC361. All isolates harbored the mecA gene, often linked to SCCmec_type IV variants. Multidrug resistance was attributed to efflux pumps NorA, NorC, SdrM, and LmrS alongside genes encoding beta-lactamase BlaZ and factors like ErmC and MepA. Additionally, virulence factors including adsA, sdrC, cap8D, harA, esaA, essC, isdB, geh, and lip were commonly identified. Furthermore, genes associated with heme uptake and clumping were present, highlighting their roles in S. aureus colonization and pathogenesis. Nine secondary metabolite biosynthetic gene clusters were found, of which six were common in all the strains. Numerous toxin-antitoxin systems were predicted, with ParE and ParB-like nuclease domains found to be the most prevalent toxin and antitoxin, respectively. Pan-genome analysis revealed 2007 core genes and 229 unique genes in the studied strains. Finally, the phylogenomic analysis showed that most Bangladeshi strains were grouped into two unique clades. This study provides a genomic and comparative insight into the multidrug resistance and pathogenicity of S. aureus strains, which will play a crucial role in the future antibiotic stewardship of Bangladesh.

Pangenome analyses of Clostridium butyricum provide insights into its genetic characteristics and industrial application

Clostridium butyricum is a Gram-positive anaerobic bacterium known for its ability to produce butyate. In this study, we conducted whole-genome sequencing and assembly of 14C. butyricum industrial strains collected from various parts of China. We performed a pan-genome comparative analysis of the 14 assembled strains and 139 strains downloaded from NCBI. We found that the genes related to critical industrial production pathways were primarily present in the core and soft-core gene categories. The phylogenetic analysis revealed that strains from the same clade of the phylogenetic tree possessed similar antibiotic resistance and virulence factors, with most of these genes present in the shell and cloud gene categories. Finally, we predicted the genes producing bacteriocins and botulinum toxins as well as CRISPR systems responsible for host defense. In conclusion, our research provides a desirable pan-genome database for the industrial production, food application, and genetic research of C. butyricum.

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.

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 .

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.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation

Background

Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation.

Results

Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K⁺ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair.

Conclusion

Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

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.