ORIGINAL RESEARCH article

Front. Plant Sci., 20 June 2023

Sec. Plant Pathogen Interactions

Volume 14 - 2023 | https://doi.org/10.3389/fpls.2023.1168480

Insight into biodiversity of the recently rearranged genus Dickeya

  • 1. Microbiologie Adaptation et Pathogénie, UMR 5240 CNRS, University Lyon, INSA Lyon, Villeurbanne, France

  • 2. Institute of Ecology and Environmental Sciences, Sorbonne University, CNRS, INRAE, Paris, France

Abstract

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.

Introduction

Members of the genus Dickeya belong to the family Pectobacteriaceae in the order Enterobacterales (Adeolu et al., 2016). Most of them are phytopathogenic bacteria and several species have a broad host range; they infect numerous vegetable crops and ornamental plants, including both monocot and dicot plants (Charkowski et al., 2012; Toth et al., 2021; van der Wolf et al., 2021). These pathogens cause either soft rots or vascular wilts in their plant hosts in temperate, tropical and subtropical climates (Charkowski et al., 2012). The soft rot symptoms are due to the action of bacterial pectinases associated with other plant cell wall-degrading enzymes (PCWDEs), which degrade the main structural components of the middle lamella and primary plant cell wall (Van Gijsegem et al., 2021). These devastating plant pathogens have a significant impact on agriculture, causing crop losses both in fields and during storage. Dickeya was classified among the top ten most important bacterial plant pathogens based on its economic and scientific impact (Mansfield et al., 2012). While most Dickeya strains have been isolated from diseased plants, some strains have also been found in surface water (Parkinson et al., 2014; Potrykus et al., 2016; Hugouvieux-Cotte-Pattat et al., 2019; Oulghazi et al., 2019; Pedron and Van Gijsegem, 2019; Ben Moussa et al., 2022).

The history of Dickeya classification began in 2005 with the proposal to place strains formerly designated Pectobacterium chrysanthemi (formerly Erwinia chrysanthemi) or Brenneria paradisiaca into the new genus Dickeya (Samson et al., 2005). At this time, the genus Dickeya comprised six recognized species, namely Dickeya chrysanthemi, Dickeya dadantii, Dickeya dianthicola, Dickeya dieffenbachiae, Dickeya paradisiaca, and Dickeya zeae (Samson et al., 2005). Thereafter, new changes were proposed in the genus Dickeya. Members of the species D. dieffenbachiae were reclassified as a subspecies of D. dadantii (i.e. D. dadantii subsp. dieffenbachiae) (Brady et al., 2012). Then, classification largely evolved with the identification of new Dickeya species and the contribution of genomics. Dickeya solani was identified as the causal agent of severe disease outbreaks of potatoes in Europe (van der Wolf et al., 2014). The novel species Dickeya fangzhongdai was isolated from pear trees in China (Tian et al., 2016) and from orchids in different countries (Alic et al., 2018). Recently, some members of the heterogeneous species D. zeae were reclassified into Dickeya oryzae including rice strains (Wang et al., 2020) and Dickeya parazeae (Hugouvieux-Cotte-Pattat and Van Gijsegem, 2021). Three new Dickeya species were also isolated from water sources. Dickeya aquatica was isolated from freshwaters in Scotland and Finland (Parkinson et al., 2014). Dickeya lacustris was found in lake water and in the rhizosphere of waterside plants in France (Hugouvieux-Cotte-Pattat et al., 2019). Dickeya undicola was found in water samples collected in Malaysia and France (Oulghazi et al., 2019). In addition, old strains from collection isolated in Australia from sugarcane or other Poaceae were classified as D. poaceiphila (Hugouvieux-Cotte-Pattat et al., 2020). More recently, the species D. paradisiaca has been reassigned to the new genus Musicola, including two species Musicola paradisiaca and Musicola keenii (Hugouvieux-Cotte-Pattat et al., 2021). Thus, the genus Dickeya currently comprises twelve validly accepted species: D. aquatica, D. chrysanthemi, D. dadantii, D. dianthicola, D. fangzhongdai, D. lacustris, D. oryzae, D. paradisiaca, D. poaceiphila, D. solani, D. undicola, and D. zeae.

Assigning the correct bacterial species to new and old Dickeya strains is necessary. The large host-range and non-specific symptoms of the soft-rot diseases increase the importance of taxonomy for accurate identification of the causal pathogens (Toth et al., 2021). Phenotypic analyses, mainly based on biochemical and nutritional traits, were previously used to differentiate species. This method remains interesting for a first classification of a large number of strains. However, the taxonomy has rapidly evolved with the use of genetic and genomic tools. Sequencing of the housekeeping genes recA or gapA is routinely used to refine interspecific phylogenetic positions of strains from the genus Dickeya (Parkinson et al., 2009; Suharjo et al., 2014; Cigna et al., 2017). Whole-genome sequence data are now widely used to understand the evolutionary and taxonomic relationships in bacteria. Thus, both phenotypic analysis and comparative genomics are key tools to improve the taxonomy of cultivable strains.

The virulence equipment of the phytopathogenic Dickeya species is based on the production and secretion of a battery of PCWDEs (Van Gijsegem et al., 2021). These bacteria have a high capacity of production and secretion of enzymes possessing pectinase, cellulase or protease activity. The maceration symptom is mainly due to the activity of pectate lyases of the PL1 family (Hugouvieux-Cotte-Pattat et al., 2014), acting together with several accessory pectinases, the cellulase CelZ and a few metalloproteases (Van Gijsegem et al., 2021). The multiplication of PCWDEs in D. dadantii 3937 is partly due to gene duplications leading to clusters of homologous genes, namely pelADE, pelBC, pehXVW, and prtABC (Hugouvieux-Cotte-Pattat et al, 1996; Hugouvieux-Cotte-Pattat et al, 2014). However, virulence is multifactorial and additional virulence determinants have been identified (Van Gijsegem et al., 2021). Different protein secretion systems (T1SS to T6SS) are present in Dickeya. The T2SS Out drives the secretion of extracellular pectinases and of the cellulase CelZ; it is essential for virulence and present in all Dickeya members. Most Dickeya members possess the T1SS Prt devoted to the secretion of the proteases and the T3SS Hrp secreting a few effectors. In contrast, the equipment in T4SS, T5SS or T6SS largely varies among Dickeya species (Van Gijsegem et al., 2021). Dickeya members produce an array of secondary metabolites involved in plant-bacteria interactions, competition to other microorganisms, or adaptation to diverse environments. In particular, they carry large gene clusters that encode complex non-ribosomal peptide synthetases (NRPS) and polyketide synthases (PKS) (Van Gijsegem et al., 2021).

Many Dickeya strains have been stored in international bacterial collections but this material deserves to be better valued (Broders et al., 2022). For example, the analysis of old Dickeya strains from the French Collection of Phytopathogenic Bacteria, CFBP (https://cirm-cfbp.fr/) contributed to clarify the taxonomy and epidemiology of the Pectobacterium species (Portier et al., 2020) and allowed the description of the two species D. poaceiphila and M. keenii (Hugouvieux-Cotte-Pattat et al., 2020; Hugouvieux-Cotte-Pattat et al., 2021). Most Dickeya strains isolated and analyzed before 2010s came from diseased crop plants, and there was a need to expand exploration to other environments. Dickeya strains are more difficult to isolate from environmental samples which contain many other bacteria, whereas rotten tissues contain almost only the causal agent. This approach was however possible thanks to a powerful semi-selective medium (Hélias et al., 2012) which allowed the isolation of the new aquatic species D. aquatica, D. lacustris, and D. undicola (Parkinson et al., 2014; Hugouvieux-Cotte-Pattat et al., 2019; Oulghazi et al., 2019).

The objective of this study was to present an updated view of the taxonomy of the genus Dickeya and to reassign the correct species to Dickeya strains isolated before the current classification. We resumed knowledge on the genus Dickeya and added new data from phenotypic, genomic and phylogenetic analyses. We tried to favour inter-species comparison by applying the analyses uniformly to all the Dickeya species. As the taxonomy within the Dickeya genus has been subject to several changes over the past 10 years, this update on the genomic and phenotypic biodiversity within this taxon is supported by information on the strains belonging to each Dickeya species.

Materials and methods

Strains

Strains of each species are listed in Tables 1–7; S1, with the plant hosts, country of origin and isolation date when these data were available. Musicola strains are listed in Table 1, D. aquatica, D. lacustris, D. poaceiphila and D. undicola strains in Table 2, D. dadantii strains in Tables 3; S1, D. dianthicola and D. solani strains in Tables 4; S1, D. fangzhongdai strains in Tables 5; S1, D. chrysanthemi in Table 6, andD. oryzae, D. parazeae and D. zeae strains in Tables 7; S1.

Table 1

Strain designationsOrigin: Country, year, plantIdentification criteria
Musicola paradisiaca
Species authority: Hugouvieux-Cotte-Pattat et al., 2021
CFBP4178T (NCPPB2511T)Colombia 1970, Musa paradisiacagenome, phenotype
Ech703Australia, Solanum tuberosumgenome
CFBP3696 (NCPPB4430)Cuba 1987, Musa paradisiacagenome, phenotype
CFBP3699Cuba 1987, Zea maysgenome, phenotype
CFBP1445Colombia 1972, Musa paradisiacagenes gapA, recA, dnaX, leuS
CFBP1446Colombia 1972, Musa paradisiacagenes gapA, recA, dnaX, leuS
CFBP1451Colombia 1972, Musa paradisiacagenes gapA, recA, dnaX, leuS
CFBP3477Colombia 1968, Musa paradisiacagenes gapA, recA, dnaX, leuS
CFBP2811 (NCPPB2512)Colombia 1973, Musa paradisiacagenes recA, dnaX, leuS
CFBP 4179 (NCPPB2513)Colombia 1973, Musa paradisiacagenes recA, dnaX, leuS
NCPPB2477Jamaica 1972, Musa paradisiacaphenotype
NCPPB2915Panama (<1977), Musa sp.phenotype
NCPPB2924Panama (<1977), Musa sp.phenotype
E353PRC (China)gene RNA 16S
572 (= Dickey 141)Musa paradisiacagene RNA 16S
Musicola keenii
Species authority: Hugouvieux-Cotte-Pattat et al., 2021
A3967T (CFBP722,
CFBP8732T, LMG31880T)France 1965, Solanum lycopersicongenome, phenotype

The characterized Musicola strains.

Table 2

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. aquatica
Species authority: Parkinson et al., 2014
174/2T (LMG 27354T, NCPPB 4580T)UK 2012, river watergenome, phenotype
DW0440Finland 2005, river watergenome
CSL RW240UK, river watergenome
181/2UK 2012, river watergenes gyrB, infB, rpoB
Dw054Finland 2005, river water16S rRNA gene
Dw0431Finland, river water16S rRNA gene
Dw0512Finland 2005, river water16S rRNA gene
JDA74 (CFBP8722)France 2018, Solanum dulcamara rhizospheregene gapA, phenotype
Not clearly identified strains
Ca3A, Ca3BIreland 2016, Daucus carotagene recA
D. lacustris
Species authority: Hugouvieux-Cotte-Pattat et al., 2019
S29T (CFBP8647T, LMG308990T)France 2017, lake watergenome, phenotype
isb1Pakistan 2021, Homo sapiens stoolgenome
S12 (CFBP8648)France 2017, lake watergene gapA, phenotype
S39 (CFBP8649)France 2017, Solanum dulcamara rhizospheregene gapA, phenotype
J114 (CFBP8721)France 2018, lake watergene gapA, phenotype
S15France 2017, lake watergene gapA, phenotype
S24France 2017, lake watergene gapA, phenotype
D. poaceiphila
Species authority: Hugouvieux-Cotte-Pattat et al., 2020
NCPPB569T (CFBP8731T)Australia 1958, Saccharum officinarumgenome, phenotype
CFBP2040Australia 1980, Megathyrsus maximusgenome, phenotype
CFBP1537 (SR149)Australia 1958, Saccharum officinarumgene gapA, phenotype
D. undicola
Species authority: Oulghazi et al., 2019
2B12T (CFBP8650T) LMG30903TMalaysia 2014, lake watergenome, phenotype
FVG1-MFV-O17France, 2017, surface watergenome, phenotype
FVG10-MFV-A16France, 2016, surface watergenome, phenotype
Ca8RC (Taiwan), Daucus carotagenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP7083 (CITA C-29)Spain 2005, Allium cepagene gapA

The characterized strains of D. aquatica, D. lacustris, D. poaceiphila and D. undicola.

Table 3

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. dadantii
Species authority: Brady et al., 2012
D. dadantii subsp. dadantii
DSM18020T (NCPPB898T, CFBP1269T)Comores 1960, Pelargoniumgenome
3937 (CFBP3855)France 1977, Saintpauliagenome
NCPPB3537 (IPO598)Peru 1985, potatogenome
BI3-1, Housui2-1, BI1-1Japan 2016, apple treegenome
Kousui1-1Japan 2016, japanese pear treegenome
Yana2-2, Aka1-1Japan 2016, peach treegenome
Kunimi-3Japan 2018, peach treegenome
CZ1501PRC (China) 2016, sweet potatogenome
M2-3PRC 2019, potatogenome
ICMP 9290Papua New Guinea 2005, sweet potatogenome
CFBP1245USA 1945, Philodendrongenes dnaX, leuS, recA
CFBP3695Cuba 1987, corngenes dnaX, leuS, recA
CFBP1443 (NCPPB2351)USA 1969, Syngoniumgenes dnaX, leuS, recA
CFBP1444Honduras 1973, Syngoniumgenes dnaX, leuS, recA
CFBP1449USA 1971, Aglaonemagenes dnaX, leuS, recA
CFBP1613France 1974, Euphorbiagenes dnaX, leuS, recA
CFBP1891USA, tabacogenes dnaX, leuS, recA
CFBP2014France 1974, Dieffenbachiaegenes dnaX, leuS, recA
CFBP2593Peru, potatogenes dnaX, leuS, recA
CFBP3694 (CFBP5649)Cuba 1987, tomatogenes dnaX, leuS, recA
CFBP3697Cuba 1987, sweet potatogenes dnaX, leuS, recA
CFBP3698Cuba 1987, bananagenes dnaX, leuS, recA
CFBP3780Italy, carnationgenes dnaX, leuS, recA
CFBP4152Greece 1985, Philodendrongenes dnaX, leuS, recA
CFBP4177Jamaica 1970, bananagenes dnaX, leuS, recA
MAFF301767eggplantgenes recA, dnaX, gyrB, rpoD
SUPP877carrotgenes recA, dnaX, gyrB, rpoD
SUPP2162strawberrygenes recA, dnaX, gyrB, rpoD
Dri1RC (Taiwan), Drimiopsisgenes dnaA, dnaJ, dnaX, gyrB, recN
PD598Netherlands 1885, Kalanchoegene recA
PD1132Netherlands 1985, Gymnocaliciumgene recA
PD552Netherlands 1988, Scindapsus pictusgene recA
NCPPB2477Jamaica 1972, bananagene recA
NCPPB3065Brazil 1978, potatogene recA
NCPPB3458Hungary 1986, Dieffenbachiagene recA
PD753Netherlands 1986, Eryngium alpinumgene recA
PD1713Netherlands 1990, Euphorbia miliigene recA
PD916RC (Taiwan) 1987, Packeragene recA
PD168Netherlands 1979, Gymnocaliciumgene recA
NCPPB2958USA 1977, sweet potatogene recA
NCPPB3476Papua New Guinea1986, sweet potatogene recA
NCPPB4097Denmark 2000, Euphorbiagene recA
IPO1260Germany, 2007, potatogene dnaX
IPO2017Netherlands, hyacinthgene dnaX
D. dadantii subsp. dieffenbachiae
NCPPB2976T (CFBP2051T)USA 1957, Dieffenbachiagenome
CFBP1152Italy 1962, Dieffenbachiagenes dnaX, leuS, recA
CFBP1237Germany 1959, Dieffenbachiagenes dnaX, leuS, recA
CFBP1360France 1970, Dieffenbachiagenes dnaX, leuS, recA
CFBP1870Ivory Coast 1976, Dieffenbachiagenes dnaX, leuS, recA
CFBP2597Switzerland 1963, Dieffenbachiagenes dnaX, leuS, recA
NCPPB2454UK 1972, Dieffenbachiagene recA
IPO1259Germany, potatogene dnaX
D. dadantii (unspecified subspecies)
S3-1RC (Taiwan) 2002, calla lily (arum)genome
FZ06Philippines 2018, bananagenome
A622-S1-A17France 2017, watergenome

Characterized strains of D. dadantii.

As they are too numerous, only a selection of strains is given, the complete list of D. dadantii strains and the latin name of the host plants are given en Table S1.

Table 4

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. dianthicola
Species authority: Samson et al., 2005
NCPPB453T(CFBP1200T)UK 1956, carnationgenome
NCPPB3534Netherlands 1987, potatogenome
RNS04.9France 2004, potatogenome
WV516USA 2016, potatogenome
SS70Pakistan 2017, potatogenome
ME23USA, potatogenome
GBBC2039Belgium, potatogenome
MIE34Switzerland 2013, potatogenome
S4.16.03.P2.4Morocco 2016, potatogenome
67.19USA 2019, New Guinea Impatiensgenome
CFBP1805Denmark 1977, Kalanchoegenome
CFBP1984France 1972, Dianthusgenome
CFBP2598Switzerland, 1982, Kalanchoegenome
CFBP2982France 1978, Kalanchoegenome
CFBP3706Switzerland 1986, chicorygenome
CFBP6548France 1994, chicorygenome
59WUSA 2016, watergenome
NY1528BUSA 2016, potatogenome
A260-S21-A16France 2016, watergenome
Dd31Serbia 2018, potatogenes acnA, gapA, icdA, mdh
MAFF302984Japan, 1992, yacongenes recA, dnaX, gyrB, rpoD
MAFF311149Japan 1996, Kalanchoegenes recA, dnaX, gyrB, rpoD
SUPP2525Japan 2006, carnationgenes recA, dnaX, gyrB, rpoD
SUPP2565Japan 2006, potatogenes recA, dnaX, gyrB, rpoD
CFBP1150Italy 1967, carnationgenes dnaX, leuS, recA
CFBP1244 (NCPPB1956)Netherlands 1966, Dahliagenes dnaX, leuS, recA
CFBP1274 (NCPPB429)UK 1956, carnationgenes dnaX, leuS, recA
CFBP1276 (NCPPB1385)Romania 1962, Dahliagenes dnaX, leuS, recA
CFBP1353Netherlands 1969, Begoniagenes dnaX, leuS, recA
CFBP1358France, 1970, Dahliagenes dnaX, leuS, recA
CFBP3702France, 1984, artichokegenes dnaX, leuS, recA
D18-A1Japan 2018, fleabanegenes recA, dnaX
D18-B1Japan 2018, butterburgenes recA, dnaX
D19-W1Japan 2019, watergenes recA, dnaX
S4.16.03.P2.18Morocco 2016, potatogene gapA
CH85/54Switzerland 1985, potatogene gapA
MG717687Australia 2017, potatogene recA
NCPPB394USA, 1957, Chrysanthemumgene recA
NCPPB2421 (PD823)Netherlands, 1969, Begoniagene recA
PD554Netherlands 1985, Kalanchoegene recA
PD788Netherlands 1987, Chicorygene recA
PD1325Netherlands, 1989, Kalanchoegene recA
PD1077Bangladeshs, 1988, potatogene recA
IPO1302Spain, potatogene dnaX
IPO2096Finland 2005, potatogene dnaX
D. solani
Species authority: Van der Wolf et al., 2014
IPO2222T(NCPPB4479T)Netherlands 2007, potatogenome
Ds0432-1Finland 2004, potatogenome
RNS 05.1.2AFrance 2005, potatogenome
RNS 08.23.3.1.A (3337)France 2008, potatogenome
RNS 07.7.3BFrance 2007, potatogenome
GBBC 2040Belgium 2007, potatogenome
MK16UK (Scotland), watergenome
IPO2019Netherlands 2009, hyacinthgenome
PPO 9019Netherlands 2006, muscarigenome
PPO 9134Netherlands 2008, hyacinthgenome
IFB0099Poland 2005, potatogenome
IFB0158, IFB0167Poland 2009, potatogenome
IFB0212Poland 2010, potatogenome
IFB0221, IFB0223Germany 2005, potato rhizospheregenome
IFB0231Finland 2008, potatogenome
IFB0311Poland 2011, potatogenome
IFB0487Poland 2013, potatogenome
IFB0695Poland 2014, potatogenome
IFB0417, IFB0421Portugal 2012, potatogenome
CH05026-1Switzerland 2005, potatogenome
CH07044Switzerland 2007, potatogenome
CH9635-1Switzerland 1996, potatogenome
CH9918-774Switzerland 1999, potatogenome
MIE35Switzerland 2005, potatogenome
D12, F012Russia 2010, potatogenome
M21aFrance 2014, potatogenome
Am3aFrance 2015, potatogenome
MK10Israel, potatogenome
A623-S20-A17France 2017, watergenome
CFBP7085Spain 2002, potatogenes dnaX, leuS, recA
CFBP7373Syria 2004, potatogenes dnaX, leuS, recA
CFBP5647France (Guadeloupe), tomatogene gapA
20711883UK 2007, potatogene recA
G-115Israel 2007, potatogene recA
IPO2093Finland 2005, potatogene dnaX

Characterized strains of D. dianthicola and D. solani.

Only a selection of strains is given, the complete list of D.dianthicola and D. solani strains and the latin name of the host plants are given en Table S1.

Table 5

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. fangzhongdai
Species authority: Tian et al., 2016
DSM101947T(CFBP8607T, JS5T)PRC (China) 2009, pear treegenome, phenotype
B16 (CFBP8496)Slovenia 2010, Phalaenopsisgenome, phenotype
S1Slovenia 2012, Phalaenopsisgenome, phenotype
MK7UK (Scotland), watergenome, phenotype
NCPPB3274St. Lucia 1983, Aglaonemagenome, phenotype
PA1PRC 2011, Phalaenopsisgenome
Onc5PRC 2021genome
M005Malaysia 2013, waterfallgenome
M074Malaysia 2013, waterfallgenome
ND14bMalaysia 2013, waterfallgenome
908CCanada 2020, marketgenome
AP6USA 2014, oniongenome
643bUSA 2020, Aglaonemagenome
LN1PRC 2014, pear treegenome
QZH3PRC 2014, pear treegenome
Secpp1600PRC 2016, radishgenome
Ph1 to Ph29RC (Taiwan), Phalaenopsisgenes dnaA, dnaJ, dnaX, gyrB, recN
CAS9, IAS4 and TAS1RC, Welsh oniongenes dnaA, dnaJ, dnaX, gyrB, recN
SUPP40Japan 1982, Welsh oniongenes recA, dnaX, gyrB, rpoD
SUPP420Japan 1985, Cliviagenes recA, dnaX, gyrB, rpoD
SUPP1034Japan 1988, Phalaenopsisgenes recA, dnaX, gyrB, rpoD
SUPP1152Japan 1989, Oncidiumgenes recA, dnaX, gyrB, rpoD
SUPP1164Japan 1989, Vandagenes recA, dnaX, gyrB, rpoD
SUPP1352Japan 1990, Dracaenagenes recA, dnaX, gyrB, rpoD
SUPP1399Japan 1990, Cattleyagenes recA, dnaX, gyrB, rpoD
SUPP1539Japan 1992, Irisgenes recA, dnaX, gyrB, rpoD
SUPP2451Japan 2004, Welsh oniongenes recA, dnaX, gyrB, rpoD
SUPP2586Japan 2004, tarogenes recA, dnaX, gyrB, rpoD
SUPP2737Japan 1989, Phalaenopsisgenes recA, dnaX, gyrB, rpoD
SUPP2738Japan 1990, Oncidiumgenes recA, dnaX, gyrB, rpoD
MAFF311172Japan 1998, tarogenes recA, dnaX, gyrB, rpoD
PD813Netherlands 1987, Phalaenopsisgene recA
PD1750Netherlands 1990, Yuccagene recA
NCPPB2915Panama 1977, bananagene recA
NCPPB2929Solomon Islands 1977, tarogene recA
NCPPB3211Sri Lanka 1982, orchidgene recA
NCPPB3306UK 1984, Polysciasgene recA

Characterized strains of D. fangzhongdai.

Only a selection of strains is given, the complete list of D. fangzhongdai strains and the latin name of the host plants are given en Table S1.

Table 6

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. chrysanthemi
Species authority: Samson et al., 2005
NCPPB402T (CFBP2048T)USA 1956, Chrysanthemum morifoliumgenome
NCPPB516 (CFBP1270)Danemark 1957, Parthenium argentatumgenome
Ech1591genome
NCPPB 3533USA 1985, Solanum tuberosumgenome
EC16, ATCC11662USA 1956, Chrysanthemum morifoliumgenome
L11Malaysia 2014, watergenome
ws52PRC 2017, Nicotiana tabacumgenome
A604-S21-A17France 2017, watergenome
CH91/70-1Switzerland 1991, Solanum tuberosumgene gapA
CH93/38-317-4Switzerland 1993, Solanum tuberosumgene gapA
CH93/40-24-1Switzerland 1993, Solanum tuberosumgene gapA
CH98/10Netherlands 1998, Solanum tuberosumgene gapA
CFBP1236 (NCPPB1861)USA 1945, Parthenium argentatumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP1242 (NCPPB427)USA 1957, Chrysanthemum maximumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP1275 (NCPPB1111)UK 1961, Dianthus caryophylusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP 1346Italy 1969, Chrysanthemum maximumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP1347Italy 1969, Chrysanthemum maximumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP1348Italy 1969, Chrysanthemum maximumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP1441USA, Dianthus caryophylusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP3262France 1981, Cichorum intybusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP3263France 1982, Cichorum intybusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP3701France 1981, Lycopersicon esculantumgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP3703France 1986, Helianthus annuusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP3704France (La Réunion) 1986, Cynara scolymusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP5847Brazil 1994, Daucus carotagenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP6689France 2002, Cichorum intybusgenes dnaA, dnaJ, dnaX, gyrB, recA
CFBP7086Spain 2003, Solanum tuberosumgenes dnaA, dnaJ, dnaX, gyrB, recA
NCPPB1849USA 1966, Parthenium argentatumgene recA
NCPPB2148UK 1968, Euphorbia pulcherrimagene recA
NCPPB2149UK 1968, Euphorbia pulcherrimagene recA
NCPPB2227UK 1969, Chrysanthemum morifoliumgene recA
NCPPB2309Italy, Chrysanthemum morifoliumgene recA
NCPPB2899USA 1976, Daucus carotagene recA
NCPPB3930Brazil, Lycopersicon esculentumgene recA
SUPP20Japan 1983, Chrysanthemum sp.gene recA
SUPP1844Japan 1998, Chrysanthemum sp.gene recA
MAFF302132Japan 1989, Solanum melongenagene recA
MAFF311043Japan 1992, Chrysanthemum sp.gene recA
MAFF311151Japan 1990, Cichorum intibusgene recA
IFB0284Poland 2011, watergene recA
IFB0320Poland 2011, watergene recA
IFB0336Poland 2011, watergene recA
PD720Kalanchoegene recA
PD806gene recA
IPO2117gene recA
SD17-1 to SD17-11gene recA

Characterized strains of D. chrysanthemi.

Table 7

Strain designationsOrigin: Country, year, plant/habitatIdentification criteria
D. oryzae
Species authority: Wang et al., 2020
ZYY5TPRC (China), ricegenome, phenotype
DZ2Q (CFBP8738)Italy, ricegenome, phenotype
S20 (CFBP8715)France 2017, watergenome, phenotype
FVG03France 2017, watergenome, phenotype
NCPPB3531 (CFBP8729)Australia, potatogenome, phenotype
EC1, EC2, ZJU1202PRC, ricegenome
CSL RW192UK, watergenome
BRIP64262Australiagenome
A003-S1-M15France 2015, watergenome
A642-S2-A17France 2017, watergenome
NCPPB2547India 1969, corngene gapA, phenotype
CFBP1271Egypt 1961, corngene gapA, phenotype
CFBP3707Israel 1986, watergene gapA, phenotype
CFBP4148Japan 1978, ricegene gapA, phenotype
CH91/71-2Switzerland 1991, potatogene gapA
A10-S1-M15 + 6 stFrance 2015, watergene gapA
A223-S2-A16 + 3 stFrance 2016, watergene gapA
A443-S1-J17 + 9 stFrance 2017, watergene gapA
IFB0324, 0330, 0334Poland 2011, watergene recA
IPO648 + 3 stNetherlands, potatogene recA
SUPP410Japan 1985, milletgene recA
SUPP739Japan 1977, ricegene recA
MAFF106502Japan 1984, ricegene recA
SUPP3076Japan 2014, ricegene recA
BC2880Korea, corngene recA
KFB414, 415, 417Serbia 2019gene recA
B1B2Australia 2017gene recA
DZ15SB01 + 3 stThailand 2015, corngene recA
IMI389157India, aloegene recA
SR120corngene recA
D. parazeae
Species authority: Hugouvieux-Cotte-Pattat et al., 2021
S31T (CFBP8716T, LMG8719T)France 2017, watergenome, phenotype
Ech586USA, Philodendrongenome
A586-S18-A17France 2017, watergenome
CFBP1531USA 1966, corngene gapA
CFBP1596 (NCPPB3731)France 1974, corngene gapA
NCPPB2540USA, corngene recA
SUPP27Japan 1980, corngene recA
MAFF311098Japan 1990, corngene recA
PD1619gene recA
D. zeae
Species authority: Samson et al., 2005
NCPPB2538T (CFBP2052T)USA 1970, corngenome, phenotype
NCPPB3532Australia, potatogenome, phenotype
MK19UK, watergenome
MS1PRC 2009, bananagenome
MS2PRC 2012, bananagenome
MS2014PRC 2014, bananagenome
MS2018PRC 2018, bananagenome
WH1PRC 2021, rice rootsgenome
A661-S21-A17France 2017, watergenome
A5272USA (Hawai)genome
BRIP64263Australiagenome
CFBP1268 (NCPPB1851)USA 1966, corngene gapA
CFBP6466France (Martinique) 1991, pineapplegene gapA
CFBP7084Spain 2005, watergene gapA
MS3PRC 2012, bananagene recA
NCPPB2339 (CFBP4176)USA, 1970, Chrysanthemumgene recA
NCPPB2340USA, 1970, Chrysanthemumgene recA
NCPPB2347Italy 1971, corngene recA
SUPP1158Japan 1989, Calanthegene recA
IPO651potatogene recA
SR171corngene recA
Potential new classification
FVG08France 2017, watergenome, phenotype
CE1PRC, Cannagenome
JZL7PRC 2017, Cliviagenome
PL65USA (Hawai) 2018, tarogenome
A5410USA (Hawai) 2007, pineapplegenome

Characterized strains of the Dickeya zeae complex (D. oryzae, D. parazeae, D. zeae).

Only a selection of strains is given, the complete list of strains and the latin name of the host plants are given en Table S1.

Growth conditions and phenotypic analysis

To test their growth with different carbon sources, strains were inoculated onto M63 minimal medium plates supplemented with a sole carbon source (2 g l-1). The enzyme secretion was assessed on plates containing an enzyme substrate (Hugouvieux-Cotte-Pattat et al., 2019). For the detection of pectinase activity, bacteria were grown on M63 agar plates supplemented with 2 g l-1 glycerol and 4 g l-1 polygalacturonic acid. After incubation for 24 h at 30°C, plates were flooded with a saturated solution of copper acetate. Clear haloes appear around colonies secreting pectate-lyases. Cellulase activity was detected on M63 agar plates supplemented with 2 g l-1 glycerol and 10 g l-1 carboxymethylcellulose. After incubation for 24 h at 30°C, the plates were overflowed with 10 mg ml-1 Congo red solution for 10 min and washed for 5 min with 1 M NaCl. Clear haloes surround colonies secreting cellulases. Protease production was tested on LB plates supplemented with skim milk (12.5 g l-1). After incubation at 30°C for 24 to 48 h, clear haloes are visible around colonies secreting proteases.

PCR, amplicon sequencing and recA phylogenetic analysis

The strains were subjected to PCR with primers specific for the gene gapA (gapAF and gapAR, 0.8 kb amplicon) (Cigna et al., 2017). The PCR products were sequenced in both directions with the same set of primers, using a commercial service (Microsynth France, Vaulx-en-Velin France). These gapA sequences were used for strain identification using BlastN comparison and construction of phylogenetic trees including the type strains. In addition, several recA and 16S rRNA sequences were retrieved from NCBI databases and used to define or verify the species affectation of the strain. Several strains firstly identified as E. chrysanthemi or Dickeya sp. were assigned to species that were not recognized at the time of their isolation.

The phylogenetic tree based on the recA gene was constructed using the pipeline Phylogeny.fr (http://www.phylogeny.fr/phylogeny.cgi); the nucleotide sequences are aligned with MUSCLE (Edgar, 2004), the phylogenetic tree is reconstructed using the maximum likelihood method implemented in the PhyML program (Guindon and Gascuel, 2003) and graphical representation is performed with TreeDyn (Dereeper et al., 2008). The recA sequences were retrieved from the NCBI database; accession numbers are given in the tree.

Comparative genomic analysis

The dDDH values were calculated using the webserver Genome-to-Genome Distance Calculator (GGDC) version 2.1 with the formula 2 (Meier-Kolthoff et al., 2013). The Average Nucleotide Identity (ANI) based on the Nucleotide MUMmer algorithm (Delcher et al., 2002) was calculated using the JSpecies Web Server with default parameters (Richter and Rosselló-Móra, 2009). The dDDH value of 70% and the cut-off ANI values of 95–96% between two strains were considered for species delineation (Goris et al., 2007).

To determine the phylogenetic position of Dickeya species, the phylogenomic tree was constructed from concatenated protein sequences of 963 unique homologous proteins (293566 sites). It was computed using the BioNJ distance method (Criscuolo and Gascuel, 2008). Two hundred bootstrap replicates were performed to assess the statistical support of each node.

Results and discussion

The new genus Musicola, analysis of old strains from collections

After analysis of a large panel of strains, Samson et al. (2005) proposed the species D. paradisiaca which comprised the former Brenneria paradisiaca type strain CFBP 4178T and six E. chrysanthemi strains (CFBP 1445, CFBP 1446, CFBP 1451, CFBP 3477, CFBP 3696, CFBP 3699). This classification was mainly based on phylogenetic analyses of 16S rRNA gene sequences which suggested that these strains were phylogenetically distant to all Brenneria species. Recently, genomic analyses have suggested that the differences between D. paradisiaca and other Dickeya species justify the creation of a separate genus (Pritchard et al, 2016; Pedron and Van Gijsegem, 2019). The genus Musicola was recognized in 2021, on the basis of phenotypic, phylogenetic and genomic arguments (Hugouvieux-Cotte-Pattat et al., 2021). Most D. paradisiaca strains were isolated from banana trees (Musa paradisiaca) in tropical or subtropical countries (Dickey and Victoria, 1980) (Table 1).

The species previously named D. paradisiaca was considered as the earliest branching lineage in the Dickeya genus in previous evolution studies (Duprey et al., 2019). However, data on the genetic diversity of this species was scarce, and it was poorly characterized at the genomic level. Only two genome sequences were available, including those of the type strain CFBP 4178T (Pritchard et al., 2013). Calculation of aligned fraction values between Dickeya and Musicola genomes showed that less than 33% of the genomes of these two genera could be aligned together (Hugouvieux-Cotte-Pattat et al., 2021).

A new Musicola strain was identified among poorly characterized strains of the CFBP collection. This strain, A3967 (CFBP 722), is related to M. paradisiaca CFBP 4178T but showed atypical phenotypes for sugar assimilation. Genomic comparison (ANI and dDDH values of 96.21% and 68.3%, respectively) confirmed the divergence and justified the description of a novel species Musicola keenii (Hugouvieux-Cotte-Pattat et al., 2021). Since the main phenotypic differences between M. keenii and M. paradisiaca members concern sugar assimilation, a simple distinction between the two species can be obtained by testing the bacterial growth in the presence of myo-inositol, melibiose or raffinose as the sole carbon source (Table 8). Sequencing of the genomes of two other Musicola strains, CFBP 3477 and CFBP 3699, confirms their assignation to the species M. paradisiaca, as suggested by phenotypic analysis (Table 1).

Musicola strains are less equipped than Dickeya members in PCWDEs and other virulence factors. They lack genes encoding PelA, PelI, PnlH, PehN, Rhi, PemB, GanA and Prt (Table 9). The reduced number of PCWDE genes could explain the low pectinase activity and the absence of protease secretion observed for Musicola members (Table 8). While the number of T1SS, T3SS, T4SS and T6SS varies among Dickeya species, Musicola members have none of them. In addition, they are deprived of factors important to deal with the plant stress responses occuring during infection, such as the nitric oxide dioxygenase HmpX that detoxify nitric oxide produced by plants, the suf cluster involved in the repair of damaged Fe/S clusters, and the ind cluster encoding the ROS scavenging pigment indigoïdine. These dissimilarities suggest notable differences in the virulence strategies of Dickeya and Musicola members. Indeed, virulence tests showed that the Musicola strains have a weak maceration activity on potato tubers and chicory leaves, two dicot models classically used to evaluate the virulence of Dickeya strains (Hugouvieux-Cotte-Pattat et al., 2021). With a low pectate lyase activity, Musicola strains are less equipped than Dickeya strains to degrade the high pectin content of dicot cell wall. Musicola strains may be better adapted to monocots such as Poaceae whose primary cell wall has a low pectin content (Jarvis et al., 1988).

Table 8

Carbon sourcesEnzyme secretionTemp.
SpeciesStrainDaraMelMtlXylClbInoG6PGluAGalADpsiLtartPGAPelCelPrt39°C
D. dadantii ssp. dad.3937w+++++++___++ww+
D. dadantii ssp. dief.CFBP 2051Tw_++++++___+++++
D. solaniIPO2222Tw+++w++++__+++++
D. dianthicolaCFBP 1200T_++++++++__+++w_
D. fangzhongdaiDSM101947Tw+++++++_ww+++++
D. undicola2B12Tw+++ND++++_++++++
D. poaceiphilaNCPPB 569Tw+++_____+__w_w+
D. oryzaeDZ2Qw++++++w_w_+++++
D. parazeaeS31Tw+++_++w__w+++w+
D. zeaeNCPPB 2538Tw+++w++w__w+++++
D. chrysanthemiNCPPB 402T_+++w+++__+++++
D. aquatica174/2T_+___+ND+__++++w+
D. lacustrisS29T_+_+_+ND+__++++++
M. paradisiacaNCPPB 2511Tw+_____+___+w+_+
M. keeniiCFBP 8732Tw__+_+++_+_+w+_+

Phenotypic comparison of the different Dickeya species.

For carbon source assimilation, the sign - indicates no growth at 72h; +, indicates growth at 24 h; w, indicates weak growth (visible after 48 or 72 h). Dara, D-arabinose; Mel, D-melibiose; Mtl, mannitol; Xyl, D-xylose; Clb, cellobiose; Ino, myo-inositol; G6P, glucose-6-phosphate; GalA, galactonic acid; GluA, gluconic acid; Dpsi, D-psicose; Ltart, L-tartaric acid; PGA, polygalacturonate (pectin backbone).

For PCWDE secretion: +, positive; -, negative (at 24 h). Pel, pectinase; Cel, cellulase; Prt, protease. +, positive; w, weak; -, negative. Growth at the temperature of 39°C was observed after 48 h in LB medium. ND, nt determined.

Only 16 Musicola strains have been characterized up to now (Table 1). These strains were mainly isolated from banana and in a few cases from another monocot, corn, and from dicots of the Solanaceae family, potato and tomato. They were isolated between 1968 and 1987 (Table 1). Unfortunately, no Musicola strain has been isolated more recently.

The water species D. aquatica and D. lacustris, true Dickeya or members of another genus?

Some water isolates belong to phylogenetic groups differing from plant-related Dickeya species. The two species D. aquatica and D. lacustris form a clade distinct from the other Dickeya species in phylogenomic studies (Duprey et al., 2019; Hugouvieux-Cotte-Pattat et al., 2019) (Figure 1). Some D. lacustris strains were found associated to the rhizosphere of bittersweet nightshade (Solanum dulcamara), a weed whose roots are in contact with water (Hugouvieux-Cotte-Pattat et al., 2019). Since water strains deserve less interest than pathogenic isolates, they were less studied. Only a few genomes have been sequenced, three for D. aquatica and one for D. lacustris (Table 2). Very recently, a new D. lacustris genome sequence (GCA_027571425.1) became available, it corresponds to strain isb1 isolated in Pakistan from Human stool. The surprising origin of this strain needs further investigation.

Figure 1

Although there is no formal definition of genus delineation on the basis of genome similarity, recent approaches used the genome aligned fraction to discriminate at genus level, showing that genus boundary corresponds to aligned fraction values around 60% (Barco et al., 2020). Indeed, for strains of the ten species D. chrysanthemi, D. dadantii, D. dianthicola, D. fangzhongdai, D. oryzae, D. parazeae, D. poaceiphila, D. solani, D. undicola and D. zeae, paired aligned fractions correspond to at least 58% of the genome lengths (Hugouvieux-Cotte-Pattat et al., 2021). These ten species constitute a coherent group corresponding to the “core Dickeya genus” (Figure 1). With aligned fractions of 73% between D. lacustris and D. aquatica genomes, these two species correspond to a homogeneous group. In contrast, the aligned fraction of D. lacustris and D. aquatica genomes is less than 41% with other Dickeya species (Hugouvieux-Cotte-Pattat et al., 2021). At the phenotypic level, D. lacustris and D. aquatica showed some metabolic particularities in comparison to the ten species of the “core Dickeya genus”, such as the non-assimilation of xylose or mannitol (Table 8). However, their PCWDE equipment is quite similar to that of other Dickeya species. They lack a few accessory pectinases PnlG, PehN and PemB (Table 9), but they have the capacity to macerate different plants in laboratory conditions (Duprey et al., 2019). A notable difference is their ecology as D. lacustris and D. aquatica seem to inhabit only water. An exception could be an isolate from carrot in Northern Ireland which was putatively identified as D. aquatica on the basis of a partial recA sequence (Zaczek-Moczydłowska et al., 2019). However, phylogenetic analysis of the corresponding sequences (MH688057, MH688058) suggests that strain Ca3 is quite distant from D. aquatica (Figure 2). Further analysis will be required to clarify the classification of this strain. The isolation of new strains and sequencing of new genomes is essential in order to obtain more information on the phenotypic and genetic diversity of these water species.

Table 9

Dickeya sp.Musicola sp.
Enzymatic activityCAZY familyGene nameDdadDsolDdianDfangDundDpoaDorDparaDzeaeDchrysDaqDlacMparadMkeen
Pectin lyasePL1pnlG11111010100000
pnlH00100000010111
Pectate lyasePL1pelA/D/E33233133322322
pelB/C22222122222222
pelZ11111011111111
PL2pelW11111111111111
PL3pelI11111111111100
PL9pelL11111011111111
pelN10111111111111
pelX11111111111111
PL10pel1000010000001000
RhamnogalacturonatePL4rhiE11111000110100
lyasePL26rhiF11110000010100
PolygalacturonaseGH28pehV/W/X33211111121111
pehK10111011101111
pehN11111000110000
Pectin methylesteraseCE8pemA11111111111111
pemB11111000010000
Pectin acetylesteraseCE10paeY11111011111111
CE12paeX11111111111111
Ferruloyl esteraseCE10faeD11111111111111
faeT11111011111111
CellulaseGH5celZ11111011111111
GalactanaseGH 53ganA11111011111100
ProteaseprtA/B/C33233133333300
prtG11010011111100

Repartition of the main PCWDEs in the different Dickeya species.

Data are given for the type strain of each species. The families of enzymes that degrade or modify polysaccharides are described in the CAZy database (http://www.cazy.org/) (Drula et al., 2022). GH, glycoside hydrolase; PL, polysaccharide lyase; CE, carbohydrate esterase.

Figure 2

Analysis of old strains from collections, the species Dickeya poaceiphila

The phenotypic and genetic analysis of poorly characterized strains from collections allowed the characterization of the species Dickeya poaceiphila, including a group of strains isolated in Australia from sugarcane and other Poaceae (commonly known as grasses) (Hugouvieux-Cotte-Pattat et al., 2020). Between 1955 and 1957, a plant disease affecting sugarcane in Australia was designated by bacterial mottle (Steindl, 1964). The symptoms consisted of chlorotic striping of the leaves which were heavily invaded by bacteria as were also the stems, causing severe stunting and wilting (Dowson and Hayward, 1960). Symptoms of chlorosis and wilting were also observed on various grasses occurring near sugarcane fields, such as Pennisetum purpureum (elephant grass), Megathyrsus maximus (Guinea grass), and Brachiaria mutica (para grass). The strain NCPPB 569T was sampled during this outbreak in Australia; its genome sequence was reported in 2013 (Pritchard et al., 2013). Genomic comparisons indicated that this strain was a candidate for assignment to a novel Dickeya species (Pritchard et al., 2016; Duprey et al., 2019) and the species D. poaceiphila was recognized in 2020 (Hugouvieux-Cotte-Pattat et al., 2020). In addition to the type strain, it includes strain CFBP 2040 whose genome was also sequenced (dDDH values of 90.2%) and strain CFBP 1537 also isolated in Australia (Table 2).

The three D. poaceiphila strains showed substantial phenotypic differences compared to other Dickeya species. They produce a low level of pectinases and proteases and no cellulase activity (Table 8). This low secretion of PCWDEs is associated with a low capacity to macerate potato tubers and chicory leaves (Hugouvieux-Cotte-Pattat et al., 2020). The absence of extracellular cellulase activity in D. poaceiphila results from the absence of the cellulase gene celZ present in all other Dickeya species (Table 9). The weak protease production could be due to the presence of a single prt gene in D. poaceiphila, while several prt genes are present in other Dickeya species (Table 9). Similarly, the low level of pectinase activity could be due to a low number of pectinase genes. D. poaceiphila has only two genes encoding secreted pectate lyases of the polysaccharide lyase family 1 (PL1), one of the pelADE cluster and one of the pelBC cluster. PL1 are responsible for the major pectate lyase activity of Dickeya (Hugouvieux-Cotte-Pattat et al., 2014) and the other Dickeya species contain 5 to 6 genes encoding PL1 pectate lyases (Table 9). D. poaceiphila is the only Dickeya species lacking the pectate lyase genes pelL and pelZ. It also lacks the accessory pectinase genes pnlG, pnlH, pehK, pehN, pemB and the galactanase gene ganA. It has only one polygalacturonase gene of the cluster pehXVW (Table 9). Several gene duplications leading to clusters of pel, peh or prt genes in other Dickeya species are not found in D. poaceiphila. This restricted number of PCWDE activities and genes could explain why D. poaceiphila strains are poorly effective in causing soft rot symptoms (Tables 8, 9) and instead produce other types of symptoms, such as stunting or wilting (Dowson and Hayward, 1960).

Additional phenotypic differences between D. poaceiphila and other Dickeya species were observed for sugar assimilation, including an absence of growth with D-gluconic acid, D-fructose-6-phosphate, D-glucose-1-phosphate, D-glucose-6-phosphate, myo-inositol, or polygacturonate (PGA), the pectin backbone (Table 8). The absence of growth of D. poaceiphila on pectin or polygalacturonate may result from the inactivation of kduI, a gene involved in the intracellular pectin catabolic pathway (Hugouvieux-Cotte-Pattat et al., 1996). This gene is annotated as a pseudogene in D. poaceiphila because it has an internal deletion of 555 bp. A defect in pectin assimilation is clearly a factor that could affect the bacterial growth in macerated plant tissues. It is less troublesome for a species that does not cause maceration of its host plants but other types of symptoms.

While D. poaceiphila strains show clear differences in phenotypic and genomic features from strains of the other characterized Dickeya species, phylogenomic studies confirmed that these strains clearly belong to the genus Dickeya (Hugouvieux-Cotte-Pattat et al., 2020). However, the D. poaceiphila genomes appear to be the smallest ones (4.22-4.32 Mb) among the Dickeya species and they have the lowest GC content (52.60-52.70%) (Table 10). A reduced genome size is consistent with the numerous gene losses observed in this species and suggests an evolution of these genomes towards an adaptation to a more restricted habitat. This could correspond to an adaptation to a specialized host range and type of symptom as D. poaceiphila has only been shown to cause wilt disease on monocot herbaceous plants of the Poaceae family.

Table 10

SpeciesGenome numberGenome size (Mb)GC%Intraspecies ANIIntraspecies dDDH
D. dadantii174.66-5.3555.90-56.50>96.3>68.9
D. solani444.81-5.0756.10-56.38>98.7>89.2
D. dianthicola764.68-4.9155.60-56.00>97.3>76.4
D. fangzhongdai164.93-5.1856.44-56.90>96.1>68.5
D. undicola34.35-4.6154.50>98.9>92.9
D. poaceiphila24.02-4.3252.60-52.8098.890.3
D. oryzae104.53-4.7553.30-53.70>96.0>67.9
D. parazeae34.71-4.8253.60-53.70>98.6>88.1
D. zeae214.56-4.9353.30-53.70>96.1>68.3
D. chrysanthemi84.62-4.8154.19-54.51>96.2>68.9
D. aquatica34.34-4.5053.35-53.60>99.9>99.1
D. lacustris24.30-4.3153.10100100
M. paradisiaca24.63-4.6855.0010099.9
M. keenii14.4054.40__

Range of genome size and GC% for the Dickeya and Musicola species.

For each species, are given the number of available genomes, their size range and their GC% range (data from the NCBI genome database, February 2023) and the intraspecies ANI and dDDH values.

Only three D. poaceiphila strains have been characterized and only two genome sequences have been reported (Table 2). Since the three strains were isolated from the same country, isolation and characterization of new strains from different origins will be necessary to access to the diversity of the species D. poaceiphila and to confirm its genomic divergence and potential specialization within the genus Dickeya.

The main Dickeya clade includes D. dadantii, D. dianthicola, D. fangzhongdai, D. solani and D. undicola

Using phylogenetic studies, two clades are observed in the core Dickeya genus (Figure 1). One includes D. chrysanthemi, D. poaceiphila, D. oryzae, D. parazeae and D. zeae. The other, called the “main Dickeya clade”, regroups the five species D. dadantii, D. dianthicola, D. fangzhongdai, D. solani and D. undicola (Figure 1). These five species form a clade sharing high ANI values (89 to 94%) and share more than 50% of their protein families (around 2600) including most virulence genes (Pedron and Van Gijsegem, 2019). They have a large equipment in PCWDEs (Table 9). Their genomes encode at least ten pectate lyases (PelB, PelC, PelD, PelE, PelI, PelL, PelN, PelW, PelX and PelZ), the rhamnogalacturonate lyase RhiE, the predicted pectin lyase PnlG, the pectin acetyl esterases PaeX and PaeY, the feruloyl esterases FaeD and FaeT, and the two pectin methyl esterases PemA and PemB (Table 9). The arsenal of polygalacturonases is more variable, ranging from three to five, including at least PehN and PehX. All species of the main Dickeya clade produce the cellulase CelZ and the galactanase GanA. Their genome harbor at least two protease genes (Table 9). While the siderophore clusters cbs and acs belongs to the Dickeya core genome, most clusters encoding secondary metabolite biosynthesis are members of the Dickeya accessory genomes, as they variably distributed among species and sometimes among strains of the same species (Table 11).

Despite their high genetic proximity and sharing of virulence related genes, the closeness and host ranges of these five Dickeya species differ, ranging from narrow to wide host ranges and genomic diversities. The two species D. dadantii and D. fangzongdai are even able to cause tree diseases in orchards, such as the bleeding canker necrosis of pear trees (Tian et al., 2016) and the bacterial quick decline of apple trees (Fujikawa et al., 2019).

D. dadantii, a wide host range, the model strain 3937

Noticeably, the D. dadantii species includes the strain 3937 that has been used for a long time as a model of the genus Dickeya for genetic studies leading to the characterization of virulence factors, including the different enzymes constituting the PCWDE equipment, and the regulatory network controlling expression of the virulence genes (Hugouvieux-Cotte-Pattat et al., 1996; Hugouvieux-Cotte-Pattat et al., 2014; Reverchon et al., 2016). The 3937 genome was the first Dickeya genome to be sequenced (Glasner et al., 2011).

Among all Dickeya species, D. dadantii has the widest host range and geographical distribution. It was isolated in all continents, from several dicots and monocot plants. More than 100 strains have been isolated from at least 25 countries and 26 different plants, including apple, peach and pear trees, vegetables (potato, tomato, eggplant, carrot, sweet potato), and numerous ornamentals (Aglaonema, Dieffenbachia, Drimiopsis, Euphorbia, Kalenchoe, Pelargonium, Phylodendron, Saintpaulia, Syngonium, etc). It was occasionally found in corn (Zea mays), cactus (Gymnocalicium), alpine plants (Erygium alpinum) and water (Table 3; Table S1).

First described as a separate species, D. dieffenbachiae was included as a subspecies of D. dadantii (D. dadantii subsp. dieffenbachiae) since both species type strains shared 96.6 ANI and 71.4% dDDH values. The main D. dadantii clade was renamed D. dadantii subsp. dadantii (Brady et al., 2012). The main phenotypic difference between the two subspecies is the assimilation of melibiose (Table 8). D. dadantii subsp. dieffenbachiae strains were nearly all isolated from ornamentals of the genus Dieffenbachia (Araceae family). With the isolation of new Asiatic strains of D. dadantii (Fujikawa et al., 2019; Wei et al., 2021; Tan et al., 2022), the heterogeneity within this species may increase. The strain S3-1 isolated from an Araceae (Zantedeschia aethiopica) in Taiwan (Wei et al., 2021) is distinct from members of both D. dadantii subspecies, with 96.3-96.7 ANI values and 68-72% dDDH values. This raises the possibility of an additional D. dadantii subspecies including strains S3-1, FZ06 and A622-S1-A17 (Table 3, Figure 1). A classification at the subspecies level can be justified when each subspecies has a specific feature or is linked to a particular host. However, newly isolated strains are rarely classified at this level, except when there is a clear phenotypic difference between the subspecies. So far, no specific traits have been described for the three D. dadantii strains not assigned to a recognized subspecies. Nevertheless, with 17 genome sequences available (Table 3), the genomic diversity of the D. dadantii species is exemplified by the large range of genome sizes, from 4.66 to 5.35 Mb (Table 10).

Another intriguing type of diversity among D. dadantii members concerns the vfm cluster which is present in all Dickeya strains (Table 11). The VFM molecule plays a major role in intercellular communication by acting as a quorum sensing molecule (Nasser et al., 2013). However, differences in the enzymatic specificity of the proteins VfmO and VfmP generate the production of different analogs of the VFM molecule, which vary according to the strain (Hugouvieux-Cotte-Pattat et al., 2022). As these analogs are differentially recognized, signaling through the Vfm quorum sensing system is limited to strains belonging to compatible groups. For instance, the D. dadantii strains are divided in two groups and the production of two different VFM analogs affects the intra-species communication (Table 11).

Table 11

Dickeya sp.Musicolasp.
ClusterGene nbrProductActivityDdadDsolDdianDfangDundDpoaDorDparaDzeaeDchrysDaqDlacMparadMkeen
vfm26VFM moleculequorum-sensing++++++++++++++
VFM group (s)VFM analogsI, IVIII, IVII, III, IVIII, IVVII, IIIIII, IVIII, IVII, III, IVIVIVVIVI
cbs16chrysobactinsiderophore++++++++++++++
acs6achromobactinsiderophore++++++++++++
ind3indigoidinantioxidant+++++++++++___
sol/ssm11solanimycinantifungalv++_+v++v+++_+___
ooc20oocydin Aantifungal_++v-__v__vv-_+_
zms17zeaminephytotoxin, antibiotic_+_+__v-_______
car8carbapenemantibioticv-___+_v-_vv+++__
cfa10coronafacic acidphytotoxinv-____+________
cluster A16unknownunknown__________+___
cluster B9unknownunknown____________+_
cluster C37unknownunknown____________+_

Clusters encoding secondary metabolites in the Dickeya species.

Data are given for the sequenced genomes of each species. +, present; - absent; v, variable; v+, variable but mostly present; v- variable but mostly absent.

D. dianthicola, a potato pathogen known since 1950s

The two species D. dianthicola and D. solani deserved special attention because they caused severe damages in the important crop potato. Most analyzed genomes are from strains isolated from diseased potatoes introducing a severe bias for diversity analyses. Despite its reputation of potato pathogen, D. dianthicola has a broad host range (van der Wolf et al., 2021). It was first described in an outbreak on carnation in Europe in the 1950s. Nearly 300 strains have now been characterized, including about 80% potato isolates (Table 4; Table S1). Other isolates were reported from vegetables, such as chicory or artichoke, and ornamentals, frequently carnation, kalanchoe, begonia and dahlia (Table 4). Potato infections in Europe have been detected since the 1970s but losses caused by D. dianthicola remained generally low and sporadic (Toth et al., 2011) except in Switzeland and the Netherlands where D. dianthicola dominated till the advent of D. solani in the 2000s (Pedron et al., 2021; Pedron et al., 2022). Since 2015, D. dianthicola caused a severe blackleg outbreak in the US that originated from Maine and further spread in at least eighteen US states (Charkowski, 2018). D. dianthicola also caused a blackleg outbreak in Western Australia in 2017 (Wright et al., 2018).

All D. dianthicola genomes contain the cluster ooc involved in the biosynthesis of an oocydin A-like molecule that may contribute to the maintenance of these species in the potato environment, by favouring the competition with other microorganisms (Brual et al., 2023). At the genomic level, D. dianthicola has a truncated pectate lyase gene pelA and lacks two protease genes (Table 9). At the phenotypic level, it is the sole member of the main Dickeya clade unable to grow with D-arabinose as the sole carbon source and it is the sole Dickeya species unable to grow at 39°C (Table 8). Such phenotypic differences could reflect a divergent adaptation of D. dianthicola to environmental conditions.

Due to the large interest for this potato pathogen, more than 70 D. dianthicola genome sequences are now available (Tables 4, 10). Most show strong relatedness illustrated by ANI values greater than 99% and clonal origin of the strains isolated during the US potato outbreak (Ge et al., 2021). Analysis of additional genomes highlights the diversity in D. dianthicola, with two strains isolated from potato in the Netherlands in 1975 and from impatiens in the US in 2019 (Liu et al, 2021) (Table 4), sharing 97 to 98% ANI values with other D. dianthicola genomes (Pedron et al., 2022). A recent report of D. dianthicola presence in Asteraceae weeds (fleabane and butterbur) close to potato fields points to another possible route of transmission from weeds to potato through surface water flow (Aono et al., 2022).

D. solani, a potato pathogen emerging in the 2000s

In contrast to D. dianthicola, D. solani has a narrow host range, it was recognized as the agent of a potato outbreak in Europe in the 2000s (Slawiak et al., 2009; van der Wolf et al., 2014), then it spread in Asia (Israel, Syria, Turkey, Georgia, oriental Russia), North Africa (Morocco) and Latin America (Brazil) (van der Wolf et al., 2021). On 65 D. solani characterized strains so far (Table 4; Table S1), only two other hosts have been identified apart to potato, the ornamentals muscari (one strain) and hyacinth (2 strains), and two strains were isolated from water (Table 4). A large majority of D. solani isolates collected in Europe and the Mediterranean Basin are clonal (Khayi et al., 2015; Golanowska et al., 2018; Blin et al., 2021). Even strains isolated more than 20 years apart are highly related, indicating a high genomic stability (Pedron et al., 2021). Genomic variations in D. solani genomes consist mainly in the presence of genes related to phages and even complete prophages (Golanowska et al., 2018; Khayi et al., 2022). The first exception was observed with strain RNS05.1.2A isolated from potato in France, which showed an ANI value of only 98% with other D. solani genomes (Khayi et al., 2015). Recently, two other French strains, RNS10-105-1A and A623S-20A-17, isolated from potato and surface water respectively, were shown to belong to the same rare subclade (Khayi et al., 2022). The other genetic variations observed in D. solani results from horizontal gene transfers leading to gene replacement. A few potato strains collected in France (RNS07.7.3B, RNS13-30-1A, RNS13-31-1A and RNS13-48-1A) acquired genomic regions related to the RNS05.1.2A subclade (Khayi et al., 2022). Two D. solani strains isolated from ornamentals (PPO9019 and PPO9134) acquired gene clusters from D. dianthicola. The strain PPO9019 also possesses a 45 kb plasmid identical to a plasmid of Burholderia ambifaria (Khayi et al., 2015). Interestingly, D. solani and B. ambifaria beside sharing an identical plasmid, exhibit similarities in the O-polysaccharide compositions that both contain 6-deoxyaltropyranose (Ossowska et al., 2017). Even if divergences remain rather limited, the D. solani strains from ornamentals or water seem to show more diversity than the potato isolates. Strains showing genomic variations were collected in The Netherlands and France and it would be interesting to further analyze the potential D. solani diversification in other countries.

A genomic comparison revealed the presence of three genomic regions encoding NRPS, PKS and associated proteins which are present in all D. solani genomes but absent in D. dadantii 3937 (Garlant et al., 2013; Pedron et al., 2014). First described for some D. oryzae strains (Zhou et al., 2011), the cluster zms is found in all D. solani and D. fangzhongdai genomes (Table 11). It encodes proteins involved in the biosynthesis of toxins of the zeamine family. The zeamine produced by D. solani is involved in antibacterial activity (Effantin et al., 2021; Brual et al., 2023). The D. solani cluster ooc comprises genes involved in the biosynthesis of an oocydin A-like molecule, a compound inhibiting Ascomycetes growth (Effantin et al., 2021; Brual et al., 2023). This cluster is also found in all the D. dianthicola genomes and some strains of D. fangzongdai, D. oryzae, D. chryanthemi and D. aquatica (Table 11). The cluster sol (or ssm) produces an antifungal compound recently named solanimycin (Matilla et al., 2022). Solanimycin is active against a broad range of plant-pathogenic fungi and the human pathogen Candida albicans (Effantin et al., 2021; Matilla et al., 2022; Brual et al., 2023). The sol cluster is conserved in most Dickeya species, including all the genomes of D. fangzongdai, D. poaceiphila, D. parazeae and D. zeae and in some strains of D. dadantii, D. undicola and D. oryzae (Table 11). The presence of the sol cluster appears to be correlated with signs of horizontal genetic transfer (Matilla et al., 2022).

D. fangzhongdai, a tree and orchid pathogen

D. fangzhongdai was the first Dickeya species found to infect trees, where it causes bleeding canker necrosis (Tian et al., 2016). Most D. fangzhongdai isolates originated from Asia (Table 5, Table S1), but the species was also identified in Europe, mainly in orchids (Alic et al., 2017). Analysis of old bacterial collections showed that D. fangzhongdai was introduced to Europe as early as 1985, through the propagation of ornamental plants (Parkinson et al., 2009; Van Vaerenbergh et al., 2012; Suharjo et al., 2014; Alic et al., 2018). Beside pear trees, D. fangzhongdai is able to cause soft rot on orchids (Phalenopsis, Oncidium, Vanda, Cattleya), and it is also found in water (Alic et al., 2018) (Table 5). Unlike other members of the main Dickeya clade, D. fangzhongdai strains are able to assimilate the rare sugar D-psicose (Table 8). The D-psicose catabolic pathway was not yet characterized but this capacity was also observed for D. oryzae and D. poaceiphila (Table 8). Genome mining in the three Dickeya species able to utilize D-psicose demonstrates the correlation between the presence of the gene cluster alsRBACDEKS and the bacterial utilization of D-psicose (Figure 3A). The importance of this pathways during pathogenesis is unknown but D-psicose was shown to induce some plant defense genes, conferring plant resistance to diseases (Kano et al., 2011).

Figure 3

When analyzed for their maceration ability in laboratory conditions, D. fangzhongdai members appear to be more efficient than any other Dickeya species (Alic et al., 2018). In addition to classical PCWDE equipment, the D. fangzhongdai genome contains a gene encoding a pectate lyase of the rare family PL10 (Table 9). PL10 enzymes exhibit catalytic properties similar to those of pectate lyases of the PL1 family such as an alkaline optimal pH, a Ca2+ requirement and a preference for low methylated pectin as substrate (Hugouvieux-Cotte-Pattat et al., 2014). They adopt an (α/α)3 barrel topology in place of the parallel β-helix topology found for PL1 pectate lyases, but the catalytic machinery used by PL10 and PL1 displays a striking resemblance, presumably due to convergent evolution (Charnock et al., 2002).

D. fangzhongdai clearly shows intra-species diversity with ANI values ranging from 96 to 99% between its members. The D. fangzhongdai strains are distributed in three clades in phylogenetic trees (Figure 1). One strain, NCPPB3274, is at the limit of being part of the species, with ANI values of ~96% and dDDH values of 68.6-70.5% with other D. fangzhongdai isolates. The D. fangzhongdai strain S1 carries a 23 kb plasmid similar to an Acidovorax plasmid. Besides genes involved in plasmid replication, stabilization and conjugative transfer, this plasmid contains two streptomycin kinase genes responsible for the resistance to the antibiotic streptomycin (Alic et al., 2019).

D. fangzhongdai is particularly well armed in secondary metabolites. In addition to the clusters zms and sol involved in the biosynthesis of the antibiotic zeamine and the antifungal compound solanimycin, respectively (Table 11), all the D. fangzhongdai genomes possess additional NRPS/PKS clusters (Alic et al., 2019; Van Gijsegem et al., 2021). Two of them are predicted to be involved in the biosynthesis of cyanobactin-related and thiopeptide-related metabolites, respectively. The atypical strain NCPPB3274 harbours the cluster ooc involved in the biosynthesis of an oocydin A-like molecule.

D. undicola, only a water species?

The species D. undicola was only recently described after isolation of three strains from surface waters, one from a lake in Malaysia and two from an irrigation canal in the South of France (Oulghazi et al., 2019) (Table 2). Phylogenetic analyses identified D. fangzhongdai as the closest Dickeya species (Figure 1). Carbon sources utilization analysis identified galactonic acid as a potential discriminative character for distinguishing between these two species (Table 8). The PCWDE repertoire of D. undicola is identical to that of D. dadantii, except that it produces only one polygalacturonase of the pelXVW cluster (Table 9). The genes involved in resistance to stresses encountered during plant infection are shared by all species of the main Dickeya clade, except that D. undicola harbors a truncated flavohemoglobin gene hmpX (Van Gijsegem et al., 2021). D. undicola genomes contain the car locus involved in the synthesis of the β-lactam antibiotic, 1-carbapen-2-em-3-carboxylic acid (carbapenem) as well as genes encoding the carbapenem immunity proteins. This cluster is also found in the water species D. aquatica and D. lacustris, some strains of D. chrysanthemi, D. oryzae and D. zeae (Table 11) and some Pectobacterium species (Van Gijsegem et al., 2021).

Despite they were isolated from far regions, the three D. undicola strains are closely related as they shared ANI and dDDH values higher than 99 and 92%, respectively, but each genome carries hundreds of specific genes (Pedron and Van Gijsegem, 2019). The low number of genomes sequenced, only three, remains insufficient to relate the closeness of the accessory genomes to the geographical origin. Recently, a new D. undicola strain was isolated from carrot in Taiwan (Wei et al., 2021) and we identified as D. undicola a strain from the CFBP collection isolated from onion in 2005 in Spain (Table 2). These two isolates suggest that D. undicola is not only a water species but can also infect plants, perhaps after irrigation with contaminated water. To better understand the D. undicola diversity and habitat, it will clearly be interesting to analyze new strains and to sequence the genomes of isolates from diverse origins.

D. chrysanthemi, the type species of the genus Dickeya

The species D. chrysanthemi is the type species of the genus Dickeya (Samson et al., 2005). The type strain was isolated in 1956 from Chrysanthemum morifolium displaying signs of soft rot and wilt. Other members of the species have been isolated from various plants, mostly Chrysanthemum sp., Parthenium argentatum, carnation, sunflower, tobacco, and from several vegetables including potato, tomato, carrot, eggplant, chicory, and artichoke. Several members were also isolated from surface water (Table 6).

At the phenotypic level, D. chrysanthemi is unable to grow with D-arabinose as the sole carbon source, a feature shared by a few other Dickeya species such as D. dianthicola, D. aquatica, and D. lacustris (Table 8). The D. chrysanthemi PCWDE equipment is quite typical of the genus Dickeya; the genome just lacks the secondary pectinase genes pnlG and pehK, one pel gene of the cluster pelADE and the protease gene prtG (Table 9).

Phylogenomic data suggest that the species D. chrysanthemi contains at least two subclades. Since dDDH values of 69%−70% were observed between most distant members, at the border of the species delimitation; this species clearly includes a high genetic diversity (Table 10). However, although this species has been known for a long time, only eight D. chrysanthemi genome sequences are currently available (Table 6); this is insufficient to propose a new classification in this heterogeneous species.

The Dickeya zeae complex: D. oryzae, D. parazeae, D. zeae and maybe others

Since classification of the strains previously classified as D. zeae has recently evolved, the term “D. zeae complex” was used to describe all these isolates (Hugouvieux-Cotte-Pattat and Van Gijsegem, 2021). In addition to the species D. zeae, the D. zeae complex comprises the species D. oryzae, initially described for a rice isolate (Wang et al., 2020), and D. parazeae whose type strain is a water isolate (Hugouvieux-Cotte-Pattat and Van Gijsegem, 2021). Several genome sequences of strains belonging to the D. zeae complex are avaible in public databases but most of the corresponding strains have not been deposited in collections and are not publicly available for phenotypic analyses, restraining their microbiological characterization.

A high diversity was noticed within genomes of the D. zeae complex which form distinct clades in phylogenetic analyses (Pedron and Van Gijsegem, 2019). The two main clades includes a group of genomes clustering with the D. zeae type strain, and another group clustering with the D. oryzae type strain (Figure 1). It should be noted that the majority of strains previously classified as D. zeae are reclassified as D. oryzae (Table 7; Table S1). In addition, a separate branch is formed by strains of the species D. parazeae. The strain CE1 clusters with the strains PL65 and JZL7 isolated from taro and clivia, respectively (Hu et al., 2018; Boluk et al., 2021) (subclade II, Figure 1). These three strains share dDDH values of about 90%, and values of ~68%, ~65% and ~56% with D. zeae subclade I (Figure 1), D. parazeae and D. oryzae, respectively. Thus, they may correspond to a novel species of the D. zeae complex and, indeed, strain PL65 was recently proposed as the type strain of a new species, Dickeya colocasiae (Boluk et al., 2022), but this classification is not officially recognized. The D. zeae strain A5410, which was isolated from pineapple (Boluk et al., 2021), also forms a separate branch (Figure 1); it shares dDDH similarities of 66−70% to D. oryzae, and 57-58% to D. zeae and D. parazeae. Strain FVG08, isolated from water is another example of strains not clearly positionned in the phylogenetic tree (Figure 1); it is at the limit of belonging to the species D. oryzae with dDDH values of 67.1-69.2% with other D. oryzae strains. These examples highlight the residual heterogeneity of the D. zeae complex. Interesting, the different D. oryzae subclades show some host preference, suggesting an evolution towards host adaptation. The rice strains form a homogeneous subclade also including a millet strain. A second D. oryzae branch brings together mostly maize strains. A third branch includes isolates from water and potato; these later may result from contaminations of potato crops by irrigation water (Hugouvieux-Cotte-Pattat and Van Gijsegem, 2021).

The D. oryzae, D. parazeae and D. zeae strains exhibit high secreted pectinase, protease and cellulase activities (Table 8). Their repertoire of PCWDEs is similar to those of the main Dickeya clade (Table 9). Differences were observed only for accessory pectinases; the genes encoding PnlH, PemB and RhiF are absent in all genomes of the D. zeae complex. The genes encoding PehN and RhiE are absent in D. parazeae and D. oryzae. Some differences in their virulence equipment disclose particularities of rice strains among D. oryzae members. The D. oryzae rice strains possess the zeamine cluster zms, the gene pnlG, and the genes cyt encoding entomotoxins but not the glycopeptidase gene avrL. In contrast, the other D. oryzae strains possess avrL but not zms, pnlG, or cyt. The toxin zeamine was shown to be involved in the rice foot rot disease caused by D. oryzae (Zhou et al., 2011). Despite its name, zeamine is not produced by D. zeae but by D. solani, D. fangzhondai and some D. oryzae strains isolated from rice (Table 11).

A phenotypic difference was observed between D. oryzae and D. zeae for the assimilation of L-tartarate, a compound abundant in many fruits (Table 8) (Hugouvieux-Cotte-Pattat and Van Gijsegem, 2021). Genomic comparison allowed the identification of the gene cluster ttdRMAB-odxA encoding the L-tartarate catabolic pathway in Dickeya (Figure 3B). This pathway was previously described in Salmonella typhimurium and Escherichia coli (Hurlbert and Jakoby, 1965), but the Dickeya cluster encodes proteins of different families sharing the same function, indicating a convergent evolution among Enterobacterales. The six species D. zeae, D. parazeae D. fangzhongdai, D. aquatica, D. lacustris and D. undicola are able to assimilate L-tartarate (Table 8). Their genome contains the cluster ttdRMAB-odxA, confirming the correlation between phenotypic and genomic data.

To reassign strains previously classified as D. zeae, we took advantage of the many available recA sequences and we sequenced the gapA gene from available strains (Table 7). This study also gives information into the host range or habitat of each species. Most D. oryzae isolates come from corn and they are frequently found in potato, rice and water. Most of the D. zeae isolates originate from corn, banana, potato, and water. Similarly, the D. parazeae strains were mostly isolated from corn and water (Table 7). Thus, members of the D. zeae complex are noticeably most common on monocots (corn, rice, banana, pineapple), but also frequent in surface water and in a few dicots, such as potato.

Conclusion

Since the definition of the genus Dickeya and the six first historic species D. chrysanthemi, D. dadantii, D. dianthicola, D. dieffenbachiae, D. paradisiaca, and D. zeae (Samson et al., 2005), mainly based on phenotypes and 16S rRNA comparison, we recently faced numerous assignation modifications and new species occurrence due to sharp genomic comparisons of known Dickeya members, new disease outbreaks, data mining in old collections and enlargement of habitat surveys. Genomic analyses resulted in the re-assignation of D. paradisiaca to the new genus Musicola. They also led to the split of D. zeae into three species and putatively more, given the high diversity among the D. zeae members. The severe outbreak in potato crops in Europe in the 2000s led to the identification of the new species D. solani. Similarly, D. fangzhongdai was described in 2016 after outbreaks in apple trees and in orchids. The atypical species D. poaceiphila was described after analysis of the diversity of old strains present in collections. Surveys of surface waters allowed the characterization of three new species, D. aquatica, D. lacustris and D. undicola, absent or rarely present on plants. New investigations in other environments such as soils but also alternative hosts such as insects or nematodes could also expand the range of Dickeya diversity. It should also be noted that, apart from the Musicola strains isolated almost exclusively in Latin America, international bacterial collections contain very few isolates from Central and South America or Central Africa. Given their climatic particularities, investigations in these underscored geographical locations should also broaden our vision of the Dickeya diversity.

Even if they are usually described as pathogens presenting a large host range, some Dickeya species seem to have a preference for herbaceous monocots. For instance, D. poaceiphila was only isolated from Poaceae and members of D. oryzae, D. parazeae and D. zeae were more often isolated from monocots than from dicots (Tables 2, 7). Similarly, Musicola members were mostly identified on monocots (banana trees) (Table 1). Pectate lyases are critical virulence factors for causing the soft rot symptoms; as these enzymes have no strict substrate specificity (Hugouvieux-Cotte-Pattat et al, 2014) they allow Dickeya to attack a large range of plants. Adaptation to a preferential host may depend on the variety and abundance of PCWDEs. For instance, D. poaceiphila has the poorest PCWDE repertoire among Dickeya members and shows a low production of pectate lyases that would disfavor the breakdown of the cell walls rich in pectin. D. poaceiphila is less well equipped to degrade the high pectin content of dicot cell wall and may be better adapted to monocots whose primary cell wall has low pectin content, such as Poaceae (Jarvis et al., 1988). The inactivation of pectin catabolism in D. poaceiphila is also a strong argument indicating that pectin degradation has become a secondary virulence factor for this species that produces mostly symptoms of chlorosis and wilting.

Genomic analyses are powerful but not suitable to the first screening of a large number of isolates as they need sequencing facilities and remains costly in time and resources. Thus, the phenotypic analysis of simple growth traits should not be forgotten for a preliminary strain classification (Table 8). In addition, the sequencing of a quality marker gene, such as recA or gapA is recommended in order to obtain a reliable taxonomic assignation of newly isolated strains. In any case, a correct classification of phytopathogenic bacteria in genera and species requires many strains of various origins (in terms of host plants, habitats and countries), a phenotypic analysis of these strains and a sufficient number of sequenced genomes to propose a solid phylogeny.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

NH-C-P – conceptualization of the study, phenotypic analysis, writing, reviewing, and editing of the manuscript. JP – comparative genomics, phylogenetic analysis. FV – genome analysis, writing, reviewing, and editing of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

The work was financially supported by funding of CNRS, University Lyon 1, and INSA Lyon to UMR 5240.

Acknowledgments

We thank Veronique Utzinger for medium preparation and contribution to phenotypic analysis, and the members of the team MTSB (MAP, UMR 5240) for sharing of unpublished information. We also thank Sebastien Santini (CNRS UMR7256) and the PACA Bioinfo platform for the availability and management of the phylogeny.fr website used to construct the recA tree.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1168480/full#supplementary-material

Supplementary Table 1

Complete list of D. dadantii, D. dianthicola, D. solani, D. fangzhongdai D. oryzae, and D. zeae strains.

Supplementary Table 2

List of the genome sequences used in this study with their accession numbers and references.

Abbreviations

NCPPB, National Collection of Plant Pathogenic Bacteria, UK; CFBP, Collection Française de Bacteíries Phytopathogènes, France; IPO, Institute for Plant Protection, Netherlands; PD, Plant Protection Service, Netherlands; IFB, Intercollegiate Faculty of Biotechnology, Poland; SUPP, Shizuoka University of Plant Pathology, Japan; MAFF, Ministry of Agriculture, Forestry and Fisheries, Japan; LMG, Laboratory of Microbiology of Ghent University, Belgium; DSM, DSMZ-German Collection of Microorganisms and Cell Cultures, Germany; ICMP, International Collection of Microorganisms from Plants, New Zealand; PCWDE, plant cell wall degrading enzyme; NRPS, non-ribosomal peptide synthetase; PKS, polyketide synthase.

References

  • 1

    AdeoluM.AlnajarS.NaushadS.GuptaR. S. (2016). Genome based phylogeny and taxonomy of the ‘Enterobacteriales’: proposal for enterobacterales ord. nov. divided into the families Enterobacteriaceae, erwiniaceae fam. nov., Pectobacteriaceae fam. nov., Yersiniaceae fam. nov., Hafniaceae fam. nov., Morganellaceae fam. nov., and Budviciaceae fam. nov. Int. J. Syst. Evol. Microbiol.66, 5575–5599. doi: 10.1099/ijsem.0.001485

  • 2

    AlicS.PedronJ.DreoT.Van GijsegemF. (2019). Genomic characterisation of the new dickeya fangzhongdai species regrouping plant pathogens and environmental isolates. BMC Genomics20, 34. doi: 10.1186/s12864-018-5332-3

  • 3

    AlicS.Tusek-ZnidaricM.PeterkaM.RavnikarM.DreoT. (2017). Putative new species of the genus dickeya as major soft rot pathogens in Phalaenopsis orchid production. Plant Pathol.66, 1357–1368. doi: 10.1111/ppa.12677

  • 4

    AlicS.Van GijsegemF.PedronJ.RavnikarM.DreoT. (2018). Diversity within the novel dickeya fangzhongdai sp., isolated from infected orchids, water and pears. Plant Pathol.67, 1612–1620. doi: 10.1111/ppa.12866

  • 5

    AonoY.NakayamaT.OgawaS.FujimotoT.OhkiT.OkaN.et al. (2022). Asteraceae weeds may be an alternative host of dickeya dianthicola, a causal agent of potato blackleg in Japan. Eur. J. Plant Pathol.63, 257–268. doi: 10.1007/s10658-022-02474-1

  • 6

    BarcoR. A.GarrityG. M.ScottJ. J.AmendJ. P.NealsonK. H.EmersonD. (2020). A genus definition for bacteria and archaea based on a standard genome relatedness index. mBio11 (1), e02475–e02419. doi: 10.1128/mBio.02475-19

  • 7

    Ben MoussaH.BertrandC.Rochelle-NewallE.FioriniS.PedronJ.BarnyM. A. (2022). The diversity of soft rot Pectobacteriaceae along the durance river stream in the south-east of France revealed by multiple seasonal surveys. Phytopathology112 (8), 1676–1685. doi: 10.1094/phyto-12-21-0515-r

  • 8

    BlinP.RobicK.KhayiS.CignaJ.MunierE.DewaegeneireP.et al. (2021). Pattern and causes of the establishment of the invasive bacterial potato pathogen Dickeya solani and of the maintenance of the resident pathogen d. dianthicola. Mol. Ecol.30, 608–624. doi: 10.1111/mec.15751

  • 9

    BolukG.ArizalaD.DobhalS.ZhangJ.HuJ.AlvarezA. M.et al. (2021). Genomic and phenotypic biology of novel strains of dickeya zeae isolated from pineapple and taro in Hawaii: insights into genome plasticity, pathogenicity, and virulence determinants. Front. Plant Sci.12, 663851. doi: 10.3389/fpls.2021.663851

  • 10

    BolukG.DobhalS.ArizalaD.AlvarezA. M.ArifM. (2022). Dickeya colocasiae sp. nov. isolated from wetland taro, Colocasia esculentum. bioRxiv. doi: 10.1101/2022.01.14.476417

  • 11

    BradyC. L.CleenwerckI.DenmanS.VenterS. N.Rodríguez-PalenzuelaP.CoutinhoT. A.et al. (2012). Proposal to reclassify Brenneria quercina (Hildebrand & schroth 1967) hauben et al. 1999 into a novel genus, Lonsdalea gen. nov., as Lonsdalea quercina comb. nov., descriptions of Lonsdalea quercina subsp. quercina comb. nov., Lonsdalea quercina subsp. iberica subsp. nov., and Lonsdalea quercina subsp. britannica subsp. nov., emendation of the description of the genus Brenneria, reclassification of Dickeya dieffenbachiae as Dickeya dadantii subsp. dieffenbachiae comb. nov., and emendation of the description of Dickeya dadantii. Int. J. Syst. Evol. Microbiol.62, 1592–1602. doi: 10.1099/ijs.0.035055-0

  • 12

    BrodersK.AspinA.BaileyJ.ChapmanT.PortierP.WeirB. S. (2022). Building more resilient culture collections: a call for increased deposits of plant-associated bacteria. Microorganisms10, 741. doi: 10.3390/microorganisms10040741

  • 13

    BrualT.EffantinG.BaltenneckJ.AttaiechL.GrosboisC.RoyerM.et al. (2023). A natural single nucleotide mutation in the small regulatory RNA ArcZ of dickeya solani switches off the antimicrobial activities against yeast and bacteria. PloS Genet.19 (4), e1010725. doi: 10.1371/journal.pgen.1010725

  • 14

    CharkowskiA. O. (2018). The changing face of bacterial soft-rot diseases. Annu. Rev. Phytopathol.56, 269–288. doi: 10.1146/annurev-phyto-080417-045906

  • 15

    CharkowskiA.BlancoC.CondemineG.ExpertD.FranzaT.HayesC.et al. (2012). The role of secretion systems and small molecules in soft-rot Enterobacteriaceae pathogenicity. Annu. Rev. Phytopathol.50, 425–449. doi: 10.1146/annurev-phyto-081211-173013

  • 16

    CharnockS. J.BrownI. E.TurkenburgJ. P.BlackG. W.DaviesG. J. (2002). Convergent evolution sheds light on the anti-beta -elimination mechanism common to family 1 and 10 polysaccharide lyases. Proc. Natl. Acad. Sci. U. S. A.99, 12067–12072. doi: 10.1073/pnas.182431199

  • 17

    CignaJ.DewaegeneireP.BeuryA.GobertV.FaureD. (2017). A gapA PCR-sequencing assay for identifying the dickeya and Pectobacterium potato pathogens. Plant Dis.101, 1278–1282. doi: 10.1094/PDIS-12-16-1810-RE

  • 18

    CriscuoloA.GascuelO. (2008). Fast NJ-like algorithms to deal with incomplete distance matrices. BMC Bioinf.9, 166. doi: 10.1186/1471-2105-9-166

  • 19

    DelcherA. L.PhillippyA.CarltonJ.SalzbergS. L. (2002). Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res.30, 2478–2483. doi: 10.1093/nar/30.11.2478

  • 20

    DereeperA.GuignonV.BlancG.AudicS.BuffetS.ChevenetF.et al. (2008). Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res.36, W465–W469. doi: 10.1093/nar/gkn180

  • 21

    DickeyR. S.VictoriaJ. I. (1980). Taxonomy and emended description of strains of erwinia isolated from Musa paradisiaca lineaus. Int. J. Syst. Bacteriol.30, 129–134. doi: 10.1099/00207713-30-1-129

  • 22

    DowsonW. J.HaywardA. C. (1960). The bacterial mottle pathogen of Queensland sugar cane. Intern. Sugar J.62, 275.

  • 23

    DrulaE.GarronM. L.DoganS.LombardV.HenrissatB.TerraponN. (2022). The carbohydrate-active enzyme database: functions and literature. Nucleic Acids Res.50 (D1), D571–D577. doi: 10.1093/nar/gkab1045

  • 24

    DupreyA.TaibN.LeonardS.GarinT.FlandroisJ. P.NasserW.et al. (2019). The phytopathogenic nature of Dickeya aquatica 174/2 and the dynamic early evolution of dickeya pathogenicity. Environ. Microbiol.21, 2809–2835. doi: 10.1111/1462-2920.14627

  • 25

    EdgarR. C. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinf.32, 1792–1797. doi: 10.1186/1471-2105-5-113

  • 26

    EffantinG.BrualT.RahbéY.Hugouvieux-Cotte-PattatN.GueguenE. (2021). Dickeya solani d s0432-1 produces an arsenal of secondary metabolites with anti-prokaryotic and anti-eukaryotic activities against bacteria, yeasts, fungi, and aphids. bioRxiv. doi: 10.1101/2021.07.19.452942

  • 27

    FujikawaT.OtaN.SasakiM.NakamuraT.IwanamiT. (2019). Emergence of apple bacterial quick decline caused by Dickeya dadantii in Japan. J. Gen. Plant Pathol.85, 314–319. doi: 10.1007/s10327-019-00852-y

  • 28

    GarlantL.KoskinenP.RouhiainenL.LaineP.PaulinL.AuvinenP.et al. (2013). Genome sequence of dickeya solani, a new soft rot pathogen of potato, suggests its emergence may be related to a novel combination of non-ribosomal peptide/polyketide synthetase clusters. Diversity5, 824–842. doi: 10.3390/d5040824

  • 29

    GeT.JiangH.TanE. H.JohnsonS. B.LarkinR. P.CharkowskiA. O.et al. (2021). Pangenomic analysis of dickeya dianthicola strains related to the outbreak of blackleg and soft rot of potato in USA. Plant Dis.105, 3946–3955. doi: 10.1094/PDIS-03-21-0587-RE

  • 30

    GlasnerJ. D.YangC. H.ReverchonS.Hugouvieux-Cotte-PattatN.CondemineG.BohinJ. P.et al. (2011). Genome sequence of the plant-pathogenic bacterium dickeya dadantii 3937. J. Bacteriol.193, 2076–2077. doi: 10.1128/JB.01513-10

  • 31

    GolanowskaM.PotrykusM.Motyka-PomagrukA.KabzaM.BacciG.GalardiniM.et al. (2018). Comparison of highly and weakly virulent Dickeya solani strains, with a view on the pangenome and panregulon of this species. Front. Microbiol.9, 1940. doi: 10.3389/fmicb.2018.01940

  • 32

    GorisJ.KonstantinidisK. T.KlappenbachJ. A.CoenyeT.VandammeP.TiedjeJ. M. (2007). DNA–DNA Hybridization values and their relationship to whole genome sequence similarities. Int. J. Syst. Evol. Microbiol.57, 81–91. doi: 10.1099/ijs.0.64483-0

  • 33

    GuindonS.GascuelO. (2003). A simple, fast and accurate algorithm to estimate larges phylogenies by maximum likelihood. Syst. Biol.52, 696–704. doi: 10.1080/10635150390235520

  • 34

    HéliasV.HamonP.HuchetE.van der WolfJ.AndrivonD. (2012). Two new effective semiselective crystal violet pectate media for isolation of pectobacterium and Dickeya: isolating pectolytic bacteria on CVP. Plant Pathol.61, 339–345. doi: 10.1111/j.1365-3059.2011.02508.x

  • 35

    HuM.LiJ.ChenR.LiW.FengL.ShiL.et al. (2018). Dickeya zeae strains isolated from rice, banana and clivia rot plants show great virulence differentials. BMC Microbiol.18, 136. doi: 10.1186/s12866-018-1300-y

  • 36

    Hugouvieux-Cotte-PattatN.Brochier-ArmanetC.FlandroisJ. P.ReverchonS. (2020). Dickeya poaceaphila sp. nov., a plant-pathogenic bacterium isolated from sugar cane (Saccharum officinarum). Int. J. Syst. Evol. Microbiol.70, 4508–4514. doi: 10.1099/ijsem.0.004306

  • 37

    Hugouvieux-Cotte-PattatN.CondemineG.NasserW.ReverchonS. (1996). Regulation of pectinolysis in erwinia chrysanthemi. Annu. Rev. Microbiol.50, 213–257. doi: 10.1146/annurev.micro.50.1.213

  • 38

    Hugouvieux-Cotte-PattatN.CondemineG.ShevchikV. E. (2014). Bacterial pectate lyases, structural and functional diversity. Environ. Microbiol. Rep.6, 427–440. doi: 10.1111/1758-2229.12166

  • 39

    Hugouvieux-Cotte-PattatN.des-CombesC. J.BriolayJ.PritchardL. (2021). Proposal for the creation of a new genus musicola gen. nov., reclassification of Dickeya paradisiaca (Samson et al. 2005) as Musicola paradisiaca comb. nov. and description of a new species Musicola keenii sp. nov. Int. J. Syst. Evol. Microbiol.71, 005037. doi: 10.1099/ijsem.0.005037

  • 40

    Hugouvieux-Cotte-PattatN.Jacot-des-CombesC.BriolayJ. (2019). Dickeya lacustris sp. nov., a water-living pectinolytic bacterium isolated from lakes in France. Int. J. Syst. Evol. Microbiol.69, 721–726. doi: 10.1099/ijsem.0.003208

  • 41

    Hugouvieux-Cotte-PattatN.RoyerM.GueguenE.Le GuenP.SüssmuthR. D.ReverchonS.et al. (2022). Allelic variations within the vfm locus define quorum sensing specificity groups in plant pathogenic bacteria of the genus Dickeya. Environ. Microbiol.24, 1467–1483. doi: 10.1111/1462-2920.15889

  • 42

    Hugouvieux-Cotte-PattatN.Van GijsegemF. (2021). Diversity within the Dickeya zeae complex, identification of dickeya zeae and Dickeya oryzae members, proposal of the novel species Dickeya parazeae sp. nov. Int. J. Syst. Evol. Microbiol.71, 005059. doi: 10.1099/ijsem.0.005059

  • 43

    HurlbertR. E.JakkobyW. B. (2022). Tartaric acid metabolism. I. Subunits of L(+)-tartaric acid dehydrase. J. Biol. Chem.240, 2772–2777. doi: 10.1016/0926-6569(64)90298-6

  • 44

    JarvisM. C.ForsythW.DuncanH. J. (1988). A survey of the pectic content of nonlignified monocot cell walls. Plant Physiol.88, 309–314. doi: 10.1104/pp.88.2.309

  • 45

    KanoA.HosotaniK.GomiK.Yamasaki-KodukoY.ShirakawaC.FukumotoT.et al. (2011). D-psicose induces upregulation of defense-related genes and resistance in rice against bacterial blight. J. Plant Physiol.168, 1852–1857. doi: 10.1016/j.jplph.2011.04.003

  • 46

    KhayiS.BlinP.PedronJ.ChongT. M.ChanK. G.MoumniM.et al. (2015). Population genomics reveals additive and replacing horizontal gene transfers in the emerging pathogen Dickeya solani. BMC Genom.16, 788. doi: 10.1186/s12864-015-1997-z

  • 47

    KhayiS.ChanK. G.FaureD. (2022). Patterns of genomic variations in the plant pathogen dickeya solani. Microorganisms10, 2254. doi: 10.3390/microorganisms10112254

  • 48

    LiuY.TylerC.HelmannT. C.StodghillP.FiliatraultM. J. (2021). Complete genome sequence resource for the necrotrophic plant pathogenic bacterium Dickeya dianthicola 67-19 isolated from new Guinea impatiens. Plant Dis.105, 1174–1176. doi: 10.1094/PDIS-09-20-1968-A

  • 49

    MansfieldJ.GeninS.MagoriS.CitovkyV.SriariyanumM.RonaldP.et al. (2012). Top 10 plant pathogenic bacteria in molecular plant pathology. Mol. Plant Pathol.13, 614–629. doi: 10.1111/j.1364-3703.2012.00804.x

  • 50

    MatillaM. A.MonsonR. E.MurphyA.SchicketanzM.RawlinsonA.DuncanC.et al. (2022). Solanimycin: biosynthesis and distribution of a new antifungal antibiotic regulated by two quorum-sensing systems. mBio13 (6), e0247222. doi: 10.1128/mbio.02472-22

  • 51

    Meier-KolthoffJ. P. M.AuchA. F.KlenkH. P.GökerM. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinf.14, 60. doi: 10.1186/1471-2105-14-60

  • 52

    NasserW.DorelC.WawrzyniakJ.Van GijsegemF.GroleauM. C.DézielE.et al. (2013). Vfm a new quorum sensing system controls the virulence of Dickeya dadantii. environ. Microbiol15, 865–880. doi: 10.1111/1462-2920.12049

  • 53

    OssowskaK.CzerwickaM.SledzW.ZoledowskaS.MotykaA.GolanowskaM.et al. (2017). The uniform structure of O-polysaccharides isolated from dickeya solani strains of different origin. Carbohydr. Res.445, 40–43. doi: 10.1016/j.carres.2017.04.001

  • 54

    OulghaziS.PedronJ.CignaJ.LauY. Y.MoumniM.Van GijsegemF.et al. (2019). Dickeya undicola sp. nov., a novel species for pectinolytic isolates from surface waters in Europe and Asia. Int. J. Syst. Evol. Microbiol.69, 2440–2444. doi: 10.1099/ijsem.0.003497

  • 55

    ParkinsonN.DeVosP.PirhonenM.ElphinstoneJ. (2014). Dickeya aquatica sp. nov., isolated from waterways. Int. J. Syst. Evol. Microbiol.64, 2264–2266. doi: 10.1099/ijs.0.058693-0

  • 56

    ParkinsonN.SteadD.BewJ.HeeneJ.Tsror LahkimL.ElphinstoneJ. (2009). Dickeya species relatedness and clade structure determined by comparison of recA sequences. Int. J. Syst. Evol. Microbiol.59, 2388–2393. doi: 10.1099/ijs.0.009258-0

  • 57

    PedronJ.MondyS.Raoul des EssartsY.Van GijsegemF.FaureD. (2014). Genomic and metabolic comparison with Dickeya dadantii 3937 reveals the emerging dickeya solani potato pathogen to display distinctive metabolic activities and T5SS/T6SS-related toxin repertoire. BMC Genom.15, 283. doi: 10.1186/1471-2164-15-283

  • 58

    PedronJ.SchaererS.KellenbergerI.Van GijsegemF. (2021). Early emergence of dickeya solani revealed by analysis of Dickeya diversity of potato blackleg and soft rot causing pathogens in Switzerland. Microorganisms9, 1187. doi: 10.3390/microorganisms9061187

  • 59

    PedronJ.van der WolfJ. M.PortierP.CaullireauE.Van GijsegemF. (2022). The broad host range plant pathogen dickeya dianthicola shows a high genetic diversity. Microorganisms10, 1024. doi: 10.3390/microorganisms10051024

  • 60

    PedronJ.Van GijsegemF. (2019). Diversity in the bacterial genus dickeya grouping plant pathogens and waterways isolates. OBM Genet.3 (4), 098. doi: 10.21926/obm.genet.1904098

  • 61

    PortierP.PedronJ.TaghoutiG.DutrieuxC.BarnyM.-A. (2020). Updated taxonomy of pectobacterium genus in the CIRM-CFBP bacterial collection: when newly described species reveal “old” endemic population. Microorganisms8, 1441. doi: 10.3390/microorganisms8091441

  • 62

    PotrykusM.GolanowskaM.SledzW.ZoledowskaS.MotykaA.KolodziejskaA.et al. (2016). Biodiversity of dickeya spp. isolated from potato plants and water sources in temperate climate. Plant Dis.100, 408–417. doi: 10.1094/PDIS-04-15-0439-RE

  • 63

    PritchardL.GloverR. H.HumphrisS.ElphinstoneJ. G.TothI. K. (2016). Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. Anal. Methods8, 12–24. doi: 10.1039/C5AY02550H

  • 64

    PritchardL.HumphrisS.SaddlerG. S.ElphinstoneJ. G.PirhonenM.TothI. K. (2013). Draft genome sequences of 17 isolates of the plant pathogenic bacterium Dickeya. Genome Announc1 (6), e00978–e00913. doi: 10.1128/genomeA.00978-13

  • 65

    ReverchonS.MuskhelisviliG.NasserW. (2016). Virulence program of a bacterial plant pathogen: the dickeya model. Prog. Mol. Biol. Transl. Sci.142, 51–92. doi: 10.1016/bs.pmbts.2016.05.005

  • 66

    RichterM.Rosselló-MóraR. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. U.S.A.106, 19126–19131. doi: 10.1073/pnas.0906412106

  • 67

    SamsonR.LegendreJ. B.ChristenR.Fischer-Le SauxM.AchouakW.GardanL. (2005). Transfer of Pectobacterium chrysanthemi (Burkholder et al 1953) Brenner et al. 1973 and Brenneria paradisiaca to the genus Dickeya gen. nov. as Dickeya chrysanthemi comb. nov. and Dickeya paradisiaca comb. nov. and delineation of four novel species, Dickeya dadantii sp. nov., Dickeya dianthicola sp. nov., Dickeya dieffenbachiae sp. nov. and Dickeya zeae sp. nov. Int. J. Syst. Evol. Microbiol.55, 1415–1427. doi: 10.1099/ijs.0.02791-0

  • 68

    SlawiakM.van BeckhovenJ. R. C. M.SpeksnijdeA. G. C. L.CzajkowskiR. L.GrabeG.van der WolfJ. M. (2009). Biochemical and genetical analysis reveal a new clade of biovar 3 dickeya spp. strains isolated from potato in Europe. Eur. J. Plant Pathol.125, 245–261. doi: 10.1007/s10658-009-9479-2

  • 69

    SteindlD. R. L. (1964). “Bacterial mottle,” in Sugar cane diseases of the world. Eds. HughesC. G.AbbottE. V.WismerC. A. (New York: Elsevier Publishing Company) 2, 2–11.

  • 70

    SuharjoR.SawadaH.TakikawaY. (2014). Phylogenetic study of Japanese dickeya spp. and development of new rapid identification methods using PCR–RFLP. J. Gen. Plant Pathol.80, 237–254. doi: 10.1007/s10327-014-0511-9

  • 71

    TanC.LiC.HuM.HuA.XueY.ZhouX.et al. (2022). Comparative pathogenomic analysis of two banana pathogenic Dickeya strains isolated from China and the Philippines. Int. J. Mol. Sci.23, 12758. doi: 10.3390/ijms232112758

  • 72

    TianY.ZhaoY.YuanX.YiJ.FanJ.XuZ.et al. (2016). Dickeya fangzhongdai sp. nov., a plant-pathogenic bacterium isolated from pear trees (Pyrus pyrifolia). Int. J. Syst. Evol. Microbiol.66, 2831–2835. doi: 10.1099/ijsem.0.001060

  • 73

    TothI. K.BarnyM.BrurbergM. B.CondemineG.CzajkowskiR.ElphinstoneJ. G.et al. (2021). “Pectobacterium and dickeya: environment to disease development,” in Plant diseases caused by dickeya and pectobacterium species. Eds. Van GijsegemF.van der WolfJ. M.TothI. K. (Cham, Switzerland: Springer International Publishing), 39–84.

  • 74

    TothI. K.van der WolfJ. M.SaddlerG.LojkowskaE.HéliasV.PirhonenM.et al. (2011). Dickeya species: an emerging problem for potato production in Europe. Plant Pathol.60, 385–399. doi: 10.1111/j.1365-3059.2011.02427.x

  • 75

    van der WolfJ.AcuñaI.De BoerS. H.BrurbergM. B.CahillG.CharkowskiA. O.et al. (2021). “Diseases caused by pectobacterium and dickeya species around the world,” in Plant diseases caused by dickeya and pectobacterium species. Eds. Van GijsegemF.van der WolfJ. M.TothI. K. (Cham, Switzerland: Springer International Publishing), 215–270.

  • 76

    van der WolfJ. M.NijhuisE. H.KowalewskaM. J.SaddlerG. S.ParkinsonN.ElphinstoneJ. G.et al. (2014). Dickeya solani sp. nov., a pectinolytic plant-pathogenic bacterium from potato (Solanum tuberosum). Int. J. Syst. Evol. Microbiol.64, 768–774. doi: 10.1099/ijs.0.052944-0

  • 77

    Van GijsegemF.Hugouvieux-Cotte-PattatN.KraepielY.LojkowskaE.MolelekiL. N.GorshkovV.et al. (2021). “Molecular interactions of pectobacterium and dickeya with plants,” in Plant diseases caused by dickeya and pectobacterium species. Eds. Van GijsegemF.van der WolfJ. M.TothI. K. (Cham, Switzerland: Springer International Publishing), 85–147.

  • 78

    Van VaerenberghJ.BaeyenS.De VosP.MaesM. (2012). Sequence diversity in the dickeya fliC gene: phylogeny of the Dickeya genus and TaqMan PCR for ‘D. solani’, new biovar 3 variant on potato in Europe. PloS One7, e35738. doi: 10.1371/journal.pone.0035738

  • 79

    WangX.HeS. W.GuoH. B.HanJ. G.ThinK. K.GaoJ. S.et al. (2020). Dickeya oryzae sp. nov., isolated from the roots of rice. Int. J. Syst. Evol. Microbiol.70, 4171–4178. doi: 10.1099/ijsem.0.004265

  • 80

    WeiX.-Y.DengW.-L.ChuC.-C. (2021). Phylogenetic and phenotypic analyses on dickeya spp. isolated from different host plants in Taiwan. J. Phytopathol.169, 678–691. doi: 10.1111/jph.13038

  • 81

    WrightD.BwyeA.BanovicM.BaulchJ.WangC.HairS.et al. (2018). First report of dickeya dianthicola in potatoes in Australia. Plant Dis.102, 2029. doi: 10.1094/PDIS-01-18-0094-PDN

  • 82

    Zaczek-MoczydłowskaM. A.FlemingC. C.YoungG. K.CampbellK.HanlonR. O. (2019). Pectobacterium and Dickeya species detected in vegetables in northern Ireland. Eur. J. Plant Pathol.154, 635–647. doi: 10.1007/s10658-019-01687-1

  • 83

    ZhouJ.ZhangH.WuJ.LiuQ.XiP.LeeJ.et al. (2011). A novel multidomain polyketide synthase is essential for zeamine production and the virulence of dickeya zeae. Mol. Plant Microbe Interact.24, 1156–1164. doi: 10.1094/MPMI-04-11-0087

Summary

Keywords

Pectobacteriaceae, soft-rot, Musicola, Dickeya, phytopathogens, pectinases

Citation

Hugouvieux-Cotte-Pattat N, Pédron J and Van Gijsegem F (2023) Insight into biodiversity of the recently rearranged genus Dickeya. Front. Plant Sci. 14:1168480. doi: 10.3389/fpls.2023.1168480

Received

17 February 2023

Accepted

22 May 2023

Published

20 June 2023

Volume

14 - 2023

Edited by

Vittoria Catara, University of Catania, Italy

Reviewed by

David J. Studholme, University of Exeter, United Kingdom; Agata Motyka-Pomagruk, Intercollegiate Faculty of Biotechnology of University of Gdańsk and Medical University of Gdańsk, Poland

Updates

Copyright

*Correspondence: Nicole Hugouvieux-Cotte-Pattat,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics