ORIGINAL RESEARCH article

Front. Microbiol., 13 February 2024

Sec. Food Microbiology

Volume 14 - 2023 | https://doi.org/10.3389/fmicb.2023.1249955

Retrospective survey of Dickeya fangzhongdai using a novel validated real-time PCR assay

  • National Institute of Biology, Ljubljana, Slovenia

Abstract

Dickeya fangzhongdai, an aggressive plant pathogen, causes symptoms on a variety of crops and ornamental plants including bleeding canker of Asian pear trees. Historical findings stress the need for a specific detection tool for D. fangzhongdai to prevent overlooking the pathogen or assigning it to general Dickeya spp. Therefore, a qualitative real-time PCR for specific detection of D. fangzhongdai has been developed and validated. The developed assay shows selectivity of 100%, diagnostic sensitivity of 76% and limit of detection with 95% confidence interval in plant matrices ranging from 311 to 2,275 cells/mL of plant extracts. The assay was successfully used in a retrospective survey of selected host plants of relevance to Europe and environmental niches relevant to D. fangzhongdai. Samples of potato tubers and plants, plants from the Malinae subtribe (apple, pear, quince, and Asian pear tree) and fresh surface water from Slovenia were analyzed. D. fangzhongdai was not detected in any plant samples, however, 12% of surface water samples were found to be positive.

1 Introduction

Dickeya fangzhongdai is one of the more recently described species of the pectinolytic genus Dickeya (). This genus groups diverse bacterial isolates that cause soft-rot disease in a variety of plant species, including economically important crops and ornamental plants, and cause wilting, black leg, and soft-rot symptoms (). Reports of soft-rot disease caused by the genus Dickeya have been limited to herbaceous plants, however, D. fangzhongdai was originally described as a causative agent of bleeding canker of Asian pear (Pyrus pyrifolia) in China (). The pathogen description was later expanded to soft rot strains isolated from various plant species, prevalently in monocots (,), and isolates from surface water (; ). Since the description of the species, multiple reports of the pathogen have been made in various host plants (see Table 1), including affecting fruits of additional tree species (). The species description has also been expanded to include strains isolated in the last century. Extension of the host range of Dickeya spp. to trees is a relatively new observation, however, it does not appear to be a unique characteristic of D. fangzhongdai, since it was also reported for D. dadantii (; ). Both species predominantly affect fruit trees such as pear, apple, peach, and jackfruit trees (; ; ; ; ). Based on the literature, the majority of the reported D. fangzhongdai strains and all tree infections have been reported in Asian countries (Table 1). There is little information regarding the economic damage and the extent of bacterial occurrence in host plants outside of Asia. The bacteria were reported as a causative agent of soft rot of orchids in commercial production in Europe () and Canada (), which originated with material from Asia and Europe, and as the causative agent of soft rot of onions in USA (). However, based on the outbreak reports, it has a more significant impact on agriculture in Asia. While data on losses is limited, several authors report outbreaks on various economically important plants such as orchids (; ; ; ), banana (), onion (; ), jack fruit (), Asian pear (; ), and even staple food crops such as taro (; ; ). Based on the reports, it can be surmised that D. fangzhongdai is well established in Asia.

TABLE 1

D. fangzhongdai strainGenBank accession1Host2Origin2Year of isolationReferences3
DSM 101947 (JS5)TCP0250034Pyrus pyrifoliaChina2009–2010
LN1CP0315054Pyrus pyrifoliaChina2009–2010
QZH3CP0315074Pyrus pyrifoliaChina2009–2010
ECM-1MT820458Pyrus pyrifoliaKorea2019
ECM-2MT820459Pyrus pyrifoliaKorea2019
ECM-3MT820460Pyrus pyrifoliaKorea2019
B16CP0872264Phalaenopsis sp.Slovenia2010
S1JXBO000000004Phalaenopsis sp.Slovenia2012
FSPAD1MK394174Phalaenopsis aphroditeChina2017
Kot1MN400213Dendrobium nobileIndia2018
Kot2MN400214Dendrobium nobileIndia2018
Kot5MN400217Dendrobium nobileIndia2018
Ph1MZ081223Phalaenopsis sp.TaiwanNA
Ph2MZ081209Phalaenopsis sp.TaiwanNA
Ph3MZ081210Phalaenopsis sp.TaiwanNA
Ph20MZ081228Phalaenopsis sp.TaiwanNA
Ph21MZ081208Phalaenopsis sp.TaiwanNA
Ph22MZ081211Phalaenopsis sp.TaiwanNA
Ph29MZ081212Phalaenopsis sp.TaiwanNA
Ph16MZ081224Phalaenopsis sp.TaiwanNA
Ph17MZ081225Phalaenopsis sp.TaiwanNA
Ph18MZ081226Phalaenopsis sp.TaiwanNA
Ph19MZ081227Phalaenopsis sp.TaiwanNA
Ph4MZ081204Phalaenopsis sp.TaiwanNA
Ph5MZ081202Phalaenopsis sp.TaiwanNA
Ph6MZ081213Phalaenopsis sp.TaiwanNA
Ph7MZ081214Phalaenopsis sp.TaiwanNA
Ph8MZ081215Phalaenopsis sp.TaiwanNA
Ph9MZ081216Phalaenopsis sp.TaiwanNA
Ph10MZ081205Phalaenopsis sp.TaiwanNA
Ph11MZ081217Phalaenopsis sp.TaiwanNA
Ph12MZ081200Phalaenopsis sp.TaiwanNA
Ph13MZ081218Phalaenopsis sp.TaiwanNA
Ph14MZ081203Phalaenopsis sp.TaiwanNA
Ph15MZ081201Phalaenopsis sp.TaiwanNA
Ph23MZ081229Phalaenopsis sp.TaiwanNA
Ph24MZ081221Phalaenopsis sp.TaiwanNA
Ph25MZ081198Phalaenopsis sp.TaiwanNA
Ph26MZ081199Phalaenopsis sp.TaiwanNA
Ph27MZ081220Phalaenopsis sp.TaiwanNA
Ph28MZ081222Phalaenopsis sp.TaiwanNA
VNO1 (LC.H1)MW281723Paphiopedilum concolorVietnam2020
VNO2 (YB.H6)MW281724Paphiopedilum concolorVietnam2020
VNO3 (LS.HD9)MW281726Phalaenopsis amabilisVietnam2020
VNO4 (TN.PD11)MW281727Dendrobium anosmumVietnam2020
VNO5 (HN.PD16)MW281728Dendrobium anosmumVietnam2020
VNO6 (SL.PD20)MW281729Dendrobium anosmumVietnam2020
VNO7 (QN.PD24)MW281730Dendrobium anosmumVietnam2020
VNO8 (GL.PD26)MW281731Dendrobium anosmumVietnam2020
VNO9 (BP.HD29)MW281732Phalaenopsis amabilisVietnam2020
908CJADCNJ000000004OrchidCanada2020
ZXC1MN853405Colocasia esculentaChinaNA
MPC2MN853406Colocasia esculentaChinaNA
Orc3-1MT613404Colocasia esculentaTaiwanNA (GenBank submission)
Orc6-2MT613403Colocasia esculentaTaiwanNA (GenBank submission)
M1L2b-2MT613402Colocasia esculentaTaiwanNA (GenBank submission)
TH11MT613401Colocasia esculentaTaiwanNA (GenBank submission)
M1L1MT613400Colocasia esculentaTaiwanNA (GenBank submission)
C2MT613399Colocasia esculentaTaiwanNA (GenBank submission)
M1D3-2MT613398Colocasia esculentaTaiwanNA (GenBank submission)
20-1MT613397Colocasia esculentaTaiwanNA (GenBank submission)
M1O1-2MT613396Colocasia esculentaTaiwanNA (GenBank submission)
M1A1-2MT613395Colocasia esculentaTaiwanNA (GenBank submission)
918-9-1MT613394Colocasia esculentaTaiwanNA (GenBank submission)
918-8-2MT613393Colocasia esculentaTaiwanNA (GenBank submission)
918-9-2MT613392Colocasia esculentaTaiwanNA (GenBank submission)
B7-15 16SMT613391Colocasia esculentaTaiwanNA (GenBank submission)
PL145MN812278Colocasia esculentaHawaii, USANA
PL146MN812277Colocasia esculentaHawaii, USANA
NCPPB 2929MZ611617Colocasia esculentaSolomon IslandsNA
NCPPB 3274CM0019794AglaonemaSt. Lucia1983
643bCP0924584Aglaonema sp.USA2020 (GenBank submission)
CAS9MZ081207Allium fistulosumTaiwanNA
IAS4MZ081219Allium fistulosumTaiwanNA
TAS1MZ081206Allium fistulosumTaiwanNA
N1MK256333Allium fistulosumTaiwan2018
AP6CP0924604Allium cepaUSA2014
GZF2-2MZ76892Musa spp.ChinaNA
GZF1-8OK668082Musa spp.ChinaNA
ZG5MW332472Pinellia ternataChina2020
MK7CM0019844River waterScotland (UK)NA
ND14b5CP0094604WaterfallMalaysia2013 (GenBank submission)
M0055JSXD000000004WaterfallMalaysia2013 (GenBank submission)
M0745JRWY000000004WaterfallMalaysia2013 (GenBank submission)
GR29MH429934Estuarine waterIndia2017 (GenBank submission)
631dMH842153Artocarpus heterophyllusMalaysia2018
131MH842152Artocarpus heterophyllusMalaysia2018
241MH197139Artocarpus heterophyllusMalaysia2018
YZY-SG-17MW160421Belamcandae RhizomaChinaNA (GenBank submission)
Secpp 1600CP0234844RadishChina2016 (GenBank submission)
VNO2RMW281725NAVietnam2020
Onc5CP0804004NAChina2021 (GenBank submission)
IPO4215OM809171Solanum tuberosumNetherlands2020
IPO4216OM809172Solanum tuberosumNetherlands2020

Reports of D. fangzhongdai species in the literature to date. Most of the reported D. fangzhongdai strains were isolated in Asia.

1GenBank accession numbers are provided for genomic (whole genome data) sequences or, if those are not available, accession numbers for 16S ribosomal RNA gene, dnaX, gapA or other accessible gene are provided.

2Sample description for GenBank submissions were derived from the sample metadata or the submission title if metadata were not available.

3For the GenBank submissions the first author and the year of the submission is stated. Each GenBank submission is clearly noted by the term “GenBank submission” in the brackets.

4Genome sequence included in the ANI analysis.

5Strains were originally described as other species, but later identified as D. fangzhongdai ().

There are currently no data on whether the strains causing bleeding canker can infect and cause soft rot symptoms of herbaceous plants in nature, however, based on inoculation experiments, pear isolates can cause symptoms on potato, tomato, cabbage, and orchids (). The ability to persist in potato plants was also reported for the orchid pathogens D. fangzhongdai S1 and B16 (). In addition, D. fangzhongdai was isolated from field grown potato tubers in a 2020 study in the Netherlands. The isolates showed high virulence in a field bioassay, capable of causing blackleg to a similar extent as Pectobacterium brasiliense (). As yet, the species is not associated with diseases of food crops outside of Asia. However, the high aggressiveness of some D. fangzhongdai strains (; ; ), the general adaptability of the genus Dickeya () and the plasticity of the species phenotype () call for caution.

Multiple D. fangzhongdai outbreaks (; ; ; ; ; ; ; ; ; ) in recent years suggest that the pathogen is present and spreading in Asia and also other continents. It was indicated before that trade plays can promote the spread of the Dickeya disease (). Moreover, there is evidence that ornamental plants may represent potential routes for the introduction of additional Dickeya species and strains with a broad host range (; ), as observed in D. solani spread to potato. Therefore, close consideration should be given to potential infection routes and the adaptability of these pathogens to other plant hosts and environments. Especially so for pathogens like D. fangzhongdai with a broad host range that predominantly consists of economically important plants with high trade rates (). However, only adequately validated diagnostic tools enable timely identification of D. fangzhongdai infected plants to support prevention of its introduction and aid epidemiological investigations. Therefore, there is a great need for the development of specific and reliable diagnostic tools that enable monitoring the presence and spread of D. fangzhongdai. Two detection tests specific to D. fangzhongdai species were previously developed: a real-time PCR test described by specifically developed and tested on D. fangzhongdai strains from Asian pear trees, and a loop-mediated isothermal amplification (LAMP) test described by that was comprehensively validated on taro, onion, and orchid matrices.

The aims of this study were to (i) develop a qualitative real-time PCR for specific detection of D. fangzhongdai in various plant matrices (and extend validation to novel matrices including potato and orchids) and ecological niches, (ii) to validate the developed test according to the guidelines of the European and Mediterranean Plant Protection Organization (EPPO) (), and compare its performance with the previously developed test described by , and (iii) to use the newly developed test in a retrospective survey to assess the presence of these bacteria in potato plants with soft rot symptoms, and surface water in Slovenia.

2 Materials and methods

2.1 Selection of host plants and niches relevant for D. fangzhongdai survey

Based on previous experience with the occurrence of soft-rot and a literature search, an informed selection of potential target plants and environmental niches was made. The literature search was performed using Google Scholar, to find reports of novel D. fanzhongdai outbreaks and isolates. Further, the GenBank database () was searched for D. fanzhongdai nucleotide sequences that have not yet been reported in any publication. Based on the search results, a table (Table 1) was created, summarizing the currently reported D. fangzhongdai strains, hosts, geographical origin, and year of isolation.

Average nucleotide identity (ANI) was calculated for all reported D. fangzhongdai isolates with publicly available genome sequences to confirm isolate identification using the Genome-based distance matrix calculator (; ).

2.2 Bacterial strains

Bacterial strains used in the study are shown in Table 2. Dickeya spp. were grown overnight on Casamino acid-Peptone-Glucose (CPG; ) medium at 28°C. Bacteria of other genera were grown overnight on yeast-extract peptone glucose agar (YPGA; ), CPG or nutrient agar (NA; ) medium at 25°C.

TABLE 2

IsolateHost plantOrigin and year of isolation
Target species
Dickeya fangzhongdai
DSMS 101947 (JS5)TPyrus pyrifoliaChina, 2009
B16Phalaenopsis sp.Slovenia, 2010
MK7River waterScotland, NA
NCPPB 3274Aglaonema sp.St. Lucia, 1983
S1Phalaenopsis sp.Slovenia, 2012
Non-target species
Dickeya solani
IPO 2222TSolanum tuberosumThe Netherlands, 2007
RNS 08.23.3.1ASolanum tuberosumFrance, 2008
7044NANA
GBBC 2040Solanum tuberosumBelgium, 2007
GBBC 500Solanum tuberosumBelgium, NA
GBBC 1021Solanum tuberosumBelgium, NA
Dickeya dadantii NCPPB 898Pelargonium capitatumComoro Islands, 1961
Dickeya dadantii subsp. dieffenbachiae LMG 25992Dieffenbachia sp.USA, 1957
Dickeya dianthicola
LMG 2485TDianthus caryophyllusUK, 1956
8823NANA
RNS 04.9Solanum tuberosumFrance, 2004
Dickeya chrysanthemi
LMG 2804TChrysanthemum morifoliumUSA, 1956
NCPPB 402Chrysanthemum morifoliumUSA, 1956
Dickeya zeae
LMG 2497Zea mays var. rugosaUSA, NA
LMG 2505TZea maysUSA, 1970
Dickeya paradisiaca LMG 2542Musa paradisiacaColombia, 1973
Dickeya aquatica NCPPB 4589River waterUK, 2008
Pectobacterium wasabiae LMG 25890Solanum tuberosumNew Zealand, 2005
Pectobacterium carotovorum pv. brasiliense PRI 3710NANA
Pectobacterium atrosepticum LMG 2386TSolanum tuberosumUK, 1957
Pectobacterium carotovorum subsp. carotovorum NCPPB 1848Cattleya sp.Brazil, 1966
Clavibacter michiganensis subsp. sepedonicus NCPPB 4053Solanum tuberosumSweden, 1994
Pseudomonas sp. 183/03-2Pyrus communisSlovenia, 2003
Escherichia coli GSPB 48NANA
Brenneria alni CFBP 3923Alnus cordataItaly, 1990
Brenneria quercina NCPPB 1852Quercus sp.USA, 1966
Proteus vulgarisenvironmental bacteriaNA
Pantoea ananatis 940/18-11Zea maysSlovenia, 2011
Pantoea agglomerans 363/18-2Zea maysSlovenia, 2018
Erwinia amylovora 106/13-1Pyrus sp.Slovenia, 41395
Serratia liquefaciens 103/14-8Lycopersicon esculentumSlovenia, 41760
Enterobacter sp. NCCPB 4168NANA, 2001
Burkholderia gladioli pv. gladioli NCPPB 1891Gladiolus sp.NA, 1966
Acidovorax avenae subsp. cattleyae NCPPB 4196Phalaenopsis sp.Brazil, 2000
Paraburkholderia caryophylli NCCPB 353NAUSA, 1954
Bacillus polymyxa NCCPB 4162Solanum tuberosumFrance, 2001
Ralstonia solanacearum 12/19-3Solanum tuberosumSlovenia, 2019

Table of bacteria isolates used to determine analytical specificity of the developed real-time PCR test.

2.3 Samples and sample preparation

2.3.1 Analytical specificity

Suspensions of target and non-target bacteria (Table 2) were prepared from overnight cultures in 10 mM phosphate buffer (PB; 1.07 g Na2HPO4, 0.4 g NaH2PO4 × 2H2O per liter of water, pH 7.2) to an approximate concentration of 106 cells/mL (Densitometer DEN-1, Biosan). Inclusivity was determined on five D. fangzhongdai isolates from three different geographical regions (Asia, America, Europe) and 3 different host niches (herbaceous plants, trees, water). Exclusivity was determined on other Dickeya spp. (17 isolates), selected bacteria from Enterobacteriaceae family (15 strains), and bacteria colonizing the same host plant niches (5 isolates) as the target bacteria (Table 2).

2.3.2 Analytical sensitivity

A D. fangzhongdai B16 and JS5T bacterial suspension with concentration of 107 cells/mL was prepared in a 10 mM PB with 30% (V/V) glycerol. DNA was extracted and standard curves were prepared by 10-fold dilutions in TE buffer (Sigma-Aldrich, Merck, Germany) with the addition of salmon sperm DNA (25 μg/mL).

2.3.3 Diagnostic sensitivity

Plant extracts were prepared from relevant asymptomatic plants, namely orchids (genus Phalaenopsis), potato (Solanum tuberosum, cultivar Carrera), and apple (Malus domestica). For the preparation of Phalaenopsis extracts, leaf material was collected in July 2019 and surface sterilized with 70% ethanol. One gram of leaf tissue was macerated in 3.5 mL of sterile 10 mM phosphate buffered saline (10 mM PBS; 1.08 g Na2HPO4, 0.4 g NaH2PO4 × 2H2O, 8 g NaCl, 1 L distilled water, pH 7.2). The supernatant was separated from the plant tissue by pipetting. A field sample of an asymptomatic potato plant, cultivar Carrera, was collected in July 2016. The surface of the sampled plant stems was cleaned and surface sterilized with 70% ethanol. The asymptomatic potato stem was cut into smaller pieces (approximately 2 cm size) and covered with sterile 10 mM PBS buffer, vortexed, and incubated for several minutes (up to 20 min) at room temperature. The supernatant was separated from the plant tissue by pipetting. An apple extract was prepared from Malus domestica asymptomatic plant material collected in July 2018. The sample material consisted of twigs that were surface sterilized with 70% ethanol. Vascular tissue was scraped from sampled twigs and covered with sterile 10 mM PBS containing 0.1% Tween 20, vortexed, and incubated with shaking for 90 min at room temperature. The supernatant was then separated from the plant tissue and centrifuged at 1,500 g for 10 min, transferred to a new tube and centrifuged at 7,000 g for 20 min. The pellets were suspended in 2 mL of 10 mM PBS.

The health status of plant extracts was confirmed with real-time PCR analysis using generic Dickeya spp. assay ().

Surface water was sampled from the Pivka River in a western part of Slovenia in August 2017. Temperature and pH of the water at the time of sampling were 22°C and 7, respectively. One liter of water was aliquoted to 250 mL and centrifuged for 20 min at 10,000 g at 4–10°C. Pellets were resuspended in 1 mL of 10 mM PB buffer. The absence of Dickeya spp. in the surface water extract was confirmed by real-time PCR analysis using generic Dickeya spp. assay ().

Standard curves of Dickeya fangzhongdai B16 in plant extracts and surface water extract were prepared by mixing bacterial suspensions with aliquots of extracts to final concentrations ranging from 107 to 101 cells/mL of plant or surface water extract.

2.3.4 Retrospective survey

A retrospective survey was performed on the collection of DNA extracts from sample material selected as described in Section “2.1 Selection of host plants and niches relevant for D. fangzhongdai survey.” Potato samples, Malinae samples, and surface water samples previously obtained in diagnostic activity in the years 2017–2021 were included in the survey. Altogether, 278 plant samples were analyzed, consisting of 130 potato samples, 148 Malinae samples and 53 surface water samples.

2.3.4.1 Samples of potato plants and tubers

Samples of potato plants and tubers with soft rot symptoms were analyzed. The surface of the sampled plants was cleaned, and surface sterilized with 70% ethanol. Symptomatic material was covered with sterile 10 mM PBS buffer, vortexed, and incubated for several minutes (up to 20 min) at room temperature. The supernatant was separated from plant tissue. DNA was extracted as described in Section “2.3 Samples and sample preparation.” Extracted DNA was stored below −15°C until analysis. Potato samples comprised of potato plants with soft rot symptoms (119 samples) and potato tubers with soft rot (11 samples).

2.3.4.2 Malinae samples

Tree samples (from the Malinae subtribe) were surface sterilized with 70% ethanol. Vascular tissue was scraped from sampled twigs and covered with sterile 10 mM PBS containing 0.1% Tween 20, vortexed, and incubated for 90 min at room temperature. The supernatant was then separated from the plant tissue and centrifuged at 1,500 g for 10 min, transferred to a new tube and centrifuged at 7,000 g for 20 min. The pellets were suspended in 2 mL of 10 mM PBS. DNA was extracted as described in Section “2.4 DNA extraction and purification.” Extracted DNA was stored below −15°C until analysis. The Malinae samples were included apple samples (84 samples), pear samples (8 samples) and Asian pear sample (1 sample) with fire blight symptoms, and mixed Malinae samples (55 samples) sampled for latent testing to fire blight.

2.3.4.3 Surface water samples

Fifty diagnostic samples of surface water were collected from different freshwater. All samples were collected in summer, and the temperature and pH of the water at the sampling site were recorded (Supplementary Table 1). Samples were analyzed as follows: one liter of water was aliquoted to 250 mL and centrifuged for 20 min at 10,000 g at 4–10°C. Pellets were resuspended in collectively 1 mL of 10 mM PB buffer. DNA was extracted as described in Section “2.3 Samples and sample preparation.” Extracted DNA was stored below −15°C until analysis.

2.4 DNA extraction and purification

DNA was extracted from 100 μL aliquots of pure bacterial suspensions, spiked plant extracts, spiked surface water, field plant samples, and surface water samples using magnetic beads-based DNA extraction on QuickPick SML Plant DNA kits (BioNobile, Finland), according to , with the minor modification of using 440 μL lysate in the downstream purification.

DNA used for analytical specificity was extracted from 500 μL of pure bacterial suspension in PB buffer using heat lysis. Samples were incubated at 95°C for 10 min in a thermoblock, and then immediately put on ice for 3 min. After centrifugation for 1 min at 6,000 rpm supernatant was collected.

2.5 Real-time PCR assay design

A D. fangzhongdai specific real-time PCR assay was designed according to . Unique diagnostic markers of D. fangzhongdai strains were identified by RUCS (). A positive dataset comprised of 10 D. fangzhongdai genomic sequences was compared to a negative dataset of 39 Dickeya spp. genomic sequences, including D. solani, D. dadantii, D. dianthicola, D. chrysanthemi, D. undicola, D. aquatica, D. zeae, and D. paradisiaca (Supplementary Table 2). The complete genome sequence of D. fangzhongdai ND14b was selected as the positive reference genome.

Specificity of the identified unique sequences was confirmed by Blastn () analysis against the whole GenBank database. Altogether, nine suitable unique sequences of sufficient length (above 100 bp) were identified. Primers and hydrolysis probes for real-time PCR were designed using Primer Express version 2.0 (Applied Biosystems). The quality of the designed assays was evaluated in silico by OligoAnalyzer Tool (IDT) and Blastn (), and experimentally.

The optimal assay, assay Df_tr (Table 3; Supplementary Table 3), designed against a transcriptional regulator gene (Dickeya_fangzhongdai_ND14b.0976; GenBank locus tag LH89_04605), was selected for validation (). Assays with poor performance or those targeting hypothetical genes or genes of extrachromosomal origin were omitted from further analysis.

TABLE 3

AssayNameSequence (5′-3′)Amplicon length
Df_trDf_tr_FGGCCGCGTCTAT GTTCTCA76 bp
Df_tr_PFAM-ACTGCATGGCGTCAATAT TTCCCCC-BHQ1
Df_tr_RACATACATTTGACACCGT CATATTTGT

Primers and probes used in real-time PCR assays designed and evaluated in this study.

2.6 Setup of the real-time PCR experiment

Real-time PCR reactions were performed on a QuantStudio 7 (Applied Biosystems, Thermo Fisher) using universal cycling conditions (2 min at 50°C, 10 min at 95°C, followed by 45 cycles of 15 s at 95°C and 1 min at 60°C, with 1.6°C/s ramping speed) according to the PCR Master Mix manufacturer’s recommendations. The reaction volume of 10 μL contained, in final concentrations: 1x TaqMan™ Universal PCR Master Mix (Applied Biosystems, Thermo Fisher), 900 nM primers (Eurofins), 200 nM probe (Eurofins), and 2 μL DNA. The QuantStudio™ real-time PCR Software 1.3 and 1.6 (Applied Biosystems, Thermo Fisher Scientific) were used for fluorescence acquisition and calculation of the threshold cycles (Cq). The baseline was set automatically, and the fluorescence threshold was set manually to intersect with the linear part of the amplification curves of all real-time PCR assays.

Analysis parameters in Df_tr validation procedure included the automatic baseline setting, and the fluorescence threshold set manually to 0.05.

Amplification of the plant endogenous sequence COX was used as an extraction and amplification control [COX; , forward primer and probe and , reverse primer)]. real-time PCR assays for non-specific detection of Dickeya spp., assay ECH () was used as a control for presence of Dickeya spp. The standard curves prepared by mixing target bacteria and plant extracts were used to determine analytical sensitivity of the novel assay and the real-time assay described by . Fluorescence thresholds for those assays were manually set to 0.1, 0.1, and 0.06 for COX, ECH, and Df_tr assay, respectively. A reaction was interpreted as positive if it produced an amplification curve and a fluorescence signal that exceeded the threshold.

Positive amplification controls and negative amplification controls were included in every real-time PCR experiment for each assay.

2.7 Validation of D. fangzhongdai specific real-time PCR assays

2.7.1 Analytical specificity and selectivity

The analytical specificity of the real-time PCR assay was tested in silico by Blastn () and experimentally by amplification of five target D. fangzhongdai strains and 37 non-target strains, including strains from eight different Dickeya genera (Table 2). Selectivity of the assay was tested on relevant plant matrixes, namely plant extracts from orchid plants, potato plants, and apple tree bark scrapings, free of disease symptoms.

2.7.2 Analytical and diagnostic sensitivity

Analytical sensitivity was determined in dilutions of DNA from pure cultures of D. fangzhongdai B16 and D. fangzhongdai JS5T, and diagnostic sensitivity was determined on standard curves of D. fangzhongdai B16 in plant extracts of Phalaenopsis, potato and apple, and surface water. Each standard curve was analyzed in triplicate. The following control systems were used to assure reliability of results and provide further information on method performance: (i) use of undiluted and diluted (1:10 in molecular grade water) DNA extracts from spiked plant extracts and surface water, and (ii) amplification of plant endogenous sequence as an extraction and amplification control (COX) (; ).

The limit of detection LOD95 was defined as the target amount giving positive results with 95% confidence and was calculated using drc package in R (; ). The slope (k) of the linear regression line between logarithmic values of cell numbers (independent variable) and Cq values (dependent variable) was used to calculate the amplification efficiency, E = (10[-1/k])–1, where a value of one corresponds to 100% amplification efficiency (). The dynamic range, i.e., the range of concentrations for which Cq values were in linear relationship with logarithms of concentrations, was determined by visually exploring the slope across sections of the Cq values × log concentration plot.

Performance of the developed real-time PCR Df_tr assay was compared to real-time PCR assay described by .

3 Results

3.1 Selection of host plants and niches relevant for D. fangzhongdai survey

Dickeya fangzhongdai isolates collected from reports in publications and the GenBank database are shown in Table 1. The majority of the species (87%) were reported from Asia. The reported isolates were predominantly isolated from soft rot symptoms on orchids (48%) and taro plants (20%). Six isolates (6%) were reported to be isolated from Asian pear trees, causing bleeding canker disease and 5 isolates (5%) were isolated from water sources. Overall, 80% of isolates were isolated from monocot plants, suggesting that D. fangzhongdai might have preference based on cotyledon types. Since D. fangzhongdai species description is relatively new, it is very likely that isolates found before the species description were assigned only to Dickeya spp. level, as was the case for NCPPB 3274. Therefore, the true list of D. fangzhongdai isolates is likely to be far more substantial. For example, it was indicated that Dickeya spp. isolates from several host plants described by correspond to D. fangzhongdai (,).

Dickeya fangzhongdai isolates with known whole genome sequence (Table 1) share above 96% average nucleotide identity (ANI) and above 86% coverage, regardless of the geographical origin or host.

Based on the literature search and previous experiences, the survey was focused on plants of agricultural importance (e.g., potato and members of the Malinae subtribe) and water samples. The latter give broader environment representation compared to individual plant samples.

3.2 Validation of the real-time PCR Df_tr assay

3.2.1 Analytical specificity and selectivity

The real-time PCR assay Df_tr, targeting a transcriptional regulator containing an amidase domain and an AraC-type DNA-binding HTH domain, was found to be specific for detection of D. fangzhongdai species. The assay exhibited 100% inclusivity (5/5 isolates) since all D. fangzhongdai isolates were reliably detected regardless of their host or geographical origin. Moreover, no cross reactivity with any of the tested non-target bacteria (37 isolates) was observed, meaning that the test is 100% exclusive to D. fangzhongdai.

Lastly, no non-specific amplification was obtained from healthy plant matrices, therefore selectivity of the assay was determined to be 100% on tested matrices, namely potato, orchid, and apple tree.

3.2.2 Analytical sensitivity

Analytical sensitivity was determined on DNA standard curves of two different D. fangzhongdai isolates, B16 and JS5T, from different environmental niches. The newly developed assay showed high analytical sensitivity (Table 4) with LOD95 below 104 cells/mL of plant extracts in all three matrices. The performance characteristics of the real-time PCR Df_tr assay were very similar to the real-time PCR described by (Supplementary Table 4). Both assays gave almost identical results tested on DNA standard curves.

TABLE 4

Dynamic range (cells/mL)1Linear regression2LOD953
FromToSlope (k)R2ELog. conc [log(cells/ml)]Cells/
mL
Residual error
DNA standard curve
D. fangzhongdai B16104107−3.61.000.913.751646.22 × 10–02
D. fangzhongdai JS5T103107−3.41.000.982.53117.21 × 10–10
Spiked plant matrix
Potato plant104106−3.30.991.012.53117.21 × 10–10
Orchids103107−3.51.000.932.53253.03 × 10–02
Apple tree104107−3.61.000.913.422753.03 × 10–02
Spiked water
Surface water104107−3.00.991.143.637763.03 × 10–02

Performance characteristics of real-time PCR Df_tr assays evaluated on bacterial suspension, spiked plant matrices and spiked surface water.

1The range of concentrations for which Cq values were in linear relationship with logarithms of concentrations.

2Linear regression of all positive samples of Cq values against logarithmic number of D. fangzhongdai cells; k: slope of the determined linear regression line; R2: average square regression coefficient; E: efficiency of amplification calculated from k.

3LOD95: limit of detection was defined as the target amount giving positive results with 95% confidence.

3.2.3 Diagnostic sensitivity

Diagnostic sensitivity was determined on spiked plant matrices and surface water. No inhibition of amplification could be observed in any of the matrices tested; however, the sensitivity of the assay varied from matrix to matrix, suggesting a matrix effect on DNA extraction procedure (Figure 1). Of all the plant matrices tested, sensitivity was most affected by the apple tree matrix, with an LOD95 of 2,275 cells/mL (Supplementary Figure 1) compared to potato and orchid matrices with LOD95 in the range of 311 – 325 cell/mL of sample. The diagnostic sensitivity of the assay was slightly lower in surface water than in plant matrices. The LOD95 was of 3,776 cells/mL of sample (Supplementary Figure 1). However, the performance characteristics of the matrices tested were not significantly different (Table 4). The Cq values were consistently below 37 at the detection limit where all parallel reactions were positive for all samples tested. Inter-run repeatability was high for all the samples and matrices tested, with coefficients of variation of Cq values below 7% within the dynamic range. The performance characteristics of the real-time PCR Df_tr assay were better compared to performance characteristics of the real-time PCR described by in all spiked matrices (Supplementary Table 4).

FIGURE 1

(B). The curves show bacterial standard curves prepared in potato matrix (green curves), in apple matrix (orange curves), in orchid matrix (blue curves) and in surface water (gray curves). The threshold line for the real-time PCR assay Df_tr is shown in orange and for in red.

The greatest difference in sensitivity between assays was observed in the plant matrices. The LOD95 of the Df_tr assay was 311 and 325 cells/mL, compared to the real-time PCR described by , which had LOD95 of 2,275 and 2,438 cells/mL, for potato and orchid matrices, respectively. Nonetheless, the difference in sensitivity was less pronounced for the apple tree matrix (LOD95 of 2,275 for the Df_tr assay and 3,776 for the real-time PCR described by . A similar difference in sensitivity was observed in surface water. The Df_tr assay showed higher sensitivity, with LOD95 of 3,776 cells/mL than the real-time PCR described by with LOD95 of 15,241 cells/mL (Figure 2). Overall, the sensitivity of Df_tr assay was better than the real-time PCR described by in plant matrices and surface water (Figure 2). In addition, fluorescence (ΔRn; Figure 1) was consistently higher in the Df_tr assay compared to real-time PCR described by .

FIGURE 2

(shown in green). The concentrations shown are expressed as log(cells/mL of plant extract), and in the brackets as cells/mL of plant extracts. The model used for both assays on spiked plant matrixes (A) is two-parameter log-logistic function (LL.2), and models used on spiked surface water (B) is log-logistic function (LL.2) for Df_tr real-time PCR assay and two-parameter Weibull function (W2.2) for real-time PCR assay described by . The dotted line denotes 95% probability of detection.

No false positives were observed for Df_tr. As expected, only samples with a bacterial concentration below LOD95 gave false negative results when compared with their known health status (Table 5). Accuracy of the test based on spiked samples was 79% and diagnostic sensitivity 76%.

TABLE 5

Known status
InfectedNon-infectedTotal
Result of testposTPaFPbTP + FP
16.00.016.0
negFNcTNdFN + TN
5.03.08.0
totalTP + FNFP + TNNe
21.03.024.0

Contingency table for real-time PCR Df_tr calculated on spiked plant samples.

If a sample was spiked with D. fangzhongdai B16 suspension, its health status was considered “infected” even if the concentration was below the expected LOD. The table combines the results for all 3 tested plant matrices.

aTrue positive;

bFalse positive;

cFalse negative;

dTrue negative;

eTotal sample count.

Samples of potato plants and tubers, and samples of Malinae trees were tested for a general presence of Dickeya spp. and D. fangzhongdai strains. The retrospective assay did not confirm the presence of D. fangzhongdai in any of the tested plant samples. However, 12% (16/130 samples) of tested samples with soft rot symptoms were positive using a non-specific Dickeya spp. real-time assay (), indicating the presence of other Dickeya species in potato. Similarly, no D. fangzhongdai nor other Dickeya spp. were detected in any of the samples of Malinae members. The Malinae samples were collected from trees that are of economic importance in the Slovenian environment, therefore the majority of the samples represent the genera Pyrus and Malus. The general presence of Dickeya spp. was confirmed in samples of surface water using real-time PCR. Dickeya spp. were detected in 70% (35 out of 50) surface water samples in low concentrations (32 ≤ Cq ≤ 39). In 6 of these samples, we also detected D. fangzhongdai, in 4 samples from August 2018 and 2 samples from August 2021. In all samples, concentrations were relatively low (34 ≤ Cq ≤ 38; Figure 3), close to the limit of detection of the assay (LOD95 of 3,776 cells/mL of sample water extract; Figure 2). The sample was considered positive if at least one reaction produced a signal above threshold and a characteristic amplification curve was present.

FIGURE 3

Positive water samples were from different freshwater sources and of different types (Table 6). The samples were collected from different parts of Slovenia, but most samples had slightly acidic pH of 6 (4 of 6 samples). The temperature of the water ranged from 17°C to 28.5°C, measured at the sampling site.

TABLE 6

Sampling location, YearType of surface waterpHWater temperature [°C]
Manče, 2021Spring7.019.0
Dobruška vas, 2021River6.018.0
Gradišče, 2018Lake6.027.0
Radehova, 2018Lake6.528.5
Vanganel, 2018Stream6.017.0
Ajdovščina, 2018Stream6.024.0

Metadata of the water samples positive presence of Dickeya spp. (real-time PCR ECH) and D. fangzhongdai (real-time PCR Df_tr).

4 Discussion

In this study, a real-time PCR assay was developed for specific detection of D. fangzhongdai, along with its validation in matrices of orchids, potatoes, and Malinae members. The assay was used in a retrospective survey of relevant ecological niches in Slovene environments.

Designed assay exhibited very good performance characteristics in the validation, which proves its suitability for the detection of D. fangzhongdai with 100% inclusivity and exclusivity and good analytical and diagnostic sensitivity. Diagnostic sensitivity of the test ranged from 1 to 10 cells per reaction (LOD95), showing that the sensitivity of the assays is close or equal to the theoretical sensitivity of the method (). No inhibition of the real-time PCR reaction could be observed in any of the tested plant matrices and the reaction efficiency was close to optimal. Compared to the real-time assay described by , the new test exhibited better diagnostic sensitivity in samples that contained plant matrices, showing better suitability for diagnostic purposes.

Dickeya fangzhongdai is the first known member of Dickeya spp. that causes disease not only on herbaceous plants, but also on trees. The majority of Dickeya fangzhongdai isolates originate from Asia, and few occurrences of D. fangzhongdai have been reported in Europe or America. This species has not yet been found to be associated with any significant plant disease in Europe’s open environment. It has been isolated from asymptomatic potato tubers in the Netherlands, but never from symptomatic plants in farmers’ fields in Europe (). Water sources seem to be an alternative habitat for Dickeya spp., as three species, namely D. aquatica, D. undicola and D. lacustris, are limited to water habitats and many others, including D. fangzhongdai, were also isolated from various water sources (; ; ; ; ). Furthermore, presence of D. fangzhongdai was indirectly confirmed in wastewaters in Slovenia by isolation of D. fangzhongdai specific bacteriophages (). The results of the screening test performed in this study correspond with the described findings. However, the importance of water as an ecological niche is not yet understood. It is not known whether water presents a transmission source or only transient ecological niche that the bacteria is able to persist in. Nevertheless, virulence genes and genes involved in virulence regulation are also conserved in isolates from water ().

In the screening test of potato plants and tubers with and without soft rot symptoms, and trees from the Malinae subtribe, D. fangzhongdai could not be detected in the tested samples. In Asian pear trees, bark tissue has been shown to be affected by D. fangzhongdai, therefore vascular tissue from twigs was selected as sample material for Malinae samples (; ). However, it is not known which tissue would be most suitable for testing asymptomatic trees for the presence of D. fangzhongdai. Sixty-five of the surface water samples tested in this study were positive for presence of Dickeya spp., and 11% of those samples also contained D. fangzhongdai. D. fangzhongdai was present at low concentration, and the limit of detection (LOD95) of the assay in surface water is 3,776 cells/mL sample water extract. Other Dickeya spp. was detected in some potato samples with soft rot symptoms, however its prevalence is approximately 5 times lower compared to water samples. Based on the results, Dickeya spp. including D. fangzhongdai have not yet entered the agricultural environment but is present at low concentrations in some water sources in Slovenia.

Repeated reports of D. fangzhongdai outbreaks in Asia indicate that the pathogen is posing a threat to cultivation of various crops, ornamental plants and trees (; ; ; ; ; ; ; ; ). There have been no reports of D. fangzhongdai outbreaks or infections of plants in the open environment in Europe, however due to lack of specific testing for D. fangzhongdai species, isolates can be overlooked or assigned to Dickeya spp. The newly developed real-time PCR is reliable, sensitive and adequately validated, and therefore a suitable detection test for D. fangzhongdai detection, identification, and monitoring. Based on the results of the retrospective survey, D. fangzhongdai seems to be present in some water sources in Slovene environment. Presence of D. fangzhongdai was not confirmed in any tested plant species, however, its persistence cannot be excluded from hosts that were not included in this study. Specific testing for D. fangzhongdai presence and accordingly implementing preventive measures, is currently the only mechanism to prevent establishment of the species in new environments and environments in which the species had been sporadically detected.

Since the beginning of the 21st century the most detrimental Dickeya spp. for in European agriculture was D. solani. In 2012 it was listed among the 10 most important bacterial pathogens because of its sudden clonal spread and impact on the potato industry under higher temperatures (). The pathogen was first isolated in 2005 and then in 2009 () but was recognized as a species only in 2014 (). However, in more recent studies it was shown that D. solani was present in potato more than a decade before the first reported outbreak. The early strains are genetically very close to the epidemic clones isolated during the 2000s outbreaks. Potentially aggressive D. solani strains in potato seeds were already present in the last century (), therefore it does not seem that genetics played an exclusive role in promoting pathogenicity, but rather an additional factor to the environmental conditions. Based on the reports, D. fangzhongdai seems to be as or even more aggressive than D. solani (; ). Even if there is currently no association of D. fangzhongdai with any host in the open European environment, it is highly likely that the pathogen will not have a problem finding a host in favorable environmental conditions. Previous experience with D. solani showed that sporadic detection of such an aggressive pathogen in the environment might be a warning sign supporting the need for specific surveying of D. fangzhongdai.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

TD and ŠA: conceptualization, methodology, validation, and funding acquisition. ŠA and KB: formal analysis, investigation, and data curation. ŠA, KB, and TD: resources. ŠA: writing—original draft preparation and visualization. TD: writing—review and editing, supervision, and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financed by the Slovenian Research and Innovation Agency (contract numbers: P4-0165 and Z1-3190) and by the Euphresco Phytosanitary ERA-NET project on taxonomy and epidemiology of Pectobacterium and Dickeya spp. in Europe, North America, and South Africa (2020-F-336). The plant material used in this study was collected by the Administration of the Republic of Slovenia for Food Safety, Veterinary and Plant Protection, of the Ministry of Agriculture and Environment and Phytosanitary Inspectorate.

Acknowledgments

We thank Prof. Ian Toth and Dr. Sonia Humphries of the James Hutton Institute for the Dickeya sp. MK7 strain; Dr. Manca Pirc, Neža Turnšek, Lidija Matičič, Špela Prijatelj Novak, and Aleš Blatnik of the NIB for diagnostic samples preparation and technical assistance, and Dr. Carissa Robyn Bleker for critical reading and English language editing.

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/fmicb.2023.1249955/full#supplementary-material

References

Summary

Keywords

molecular testing, diagnostics, plant pathogen, real-time PCR, Dickeya, survey, water

Citation

Alič Š, Bačnik K and Dreo T (2024) Retrospective survey of Dickeya fangzhongdai using a novel validated real-time PCR assay. Front. Microbiol. 14:1249955. doi: 10.3389/fmicb.2023.1249955

Received

29 June 2023

Accepted

27 December 2023

Published

13 February 2024

Volume

14 - 2023

Edited by

Mohammad Arif, University of Hawai’i at Manoa, United States

Reviewed by

Emilio Marco Stefani, University of Modena and Reggio Emilia, Italy

Gamze Boluk-Sari, Ministry of Agriculture and Forestry, Türkiye

Updates

Copyright

*Correspondence: Tanja Dreo,

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