Genomic analysis of the Phalaenopsis pathogen Dickeya sp. PA1, representing the emerging species Dickeya fangzhongdai

Dickeya fangzhongdai was prevalent and caused taro soft rot when coexisting with the Pectobacterium complex, with a preference for Araceae plants

Bacterial soft rot caused by coinfection with Dickeya spp. and Pectobacterium spp. in hosts can cause successive changes in fields, and it is difficult to prevent the spread of and control the infection. Pectobacterium spp. are prevalent in the growing areas of tuberous crops, including taro and potato. Recently, Dickeya fangzhongdai has emerged as a virulent pathogen in taro. To determine the prevalence status of the causal agents and evaluate the potential spreading risks of D. fangzhongdai, screening and taxonomic classification were performed on phytopathogenic bacteria collected from different taro-growing areas in Guangdong Province, China, and biological and genomic characteristics were further compared among typical strains from all defined species. The causative agents were verified to be phytobacterial strains of D. fangzhongdai, Pectobacterium aroidearum and Pectobacterium colocasium. P. aroidearum and P. colocasium were found to form a complex preferring Araceae plants and show intensive genomic differentiation, indicating their ancestor had adapted to taro a long time prior. Compared with Pectobacterium spp., D. fangzhongdai was more virulent to taro corms under conditions of exogenous infection and more adaptable at elevated temperatures. D. fangzhongdai strains isolated from taro possessed genomic components of additional T4SSs, which were accompanied by additional copies of the hcp-vgrG genes of the T6SS, and these contributed to the expansion of their genomes. More gene clusters encoding secondary metabolites were found within the D. fangzhongdai strains than within the Pectobacterium complex; interestingly, distinct gene clusters encoding zeamine and arylpolyene were both most similar to those in D. solani that caused potato soft rot. These comparisons provided genomic evidences for that the newly emerging pathogen was potentially equipped to compete with other pathogens. Diagnostic qPCR verified that D. fangzhongdai was prevalent in most of the taro-growing areas and coexisted with the Pectobacterium complex, while the plants enriching D. fangzhongdai were frequently symptomatic at developing corms and adjacent pseudostems and caused severe symptoms. Thus, the emerging need for intensive monitoring on D. fangzhongdai to prevent it from spreading to other taro-growing areas and to other tuberous crops like potato; the adjustment of control strategies based on different pathopoiesis characteristics is recommended.

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

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.

Comparative Pathogenomic Analysis of Two Banana Pathogenic Dickeya Strains Isolated from China and the Philippines

Dickeya is a major and typical member of soft rot Pectobacteriaceae (SRP) with a wide range of plant hosts worldwide. Previous studies have identified D. zeae as the causal agent of banana soft rot disease in China. In 2017, we obtained banana soft rot pathogen strain FZ06 from the Philippines. Genome sequencing and analysis indicated that FZ06 can be classified as D. dadantii and represents a novel subspecies of D. dadantii, which we propose to name as subsp. paradisiaca. Compared with Chinese banana soft rot pathogenic strain D. zeae MS2, strain FZ06 has a similar host range but different virulence; FZ06 is significantly less virulent to banana and potato but more virulent to Chinese cabbage and onion. Characterization of virulence factors revealed obviously less production of pectate lyases (Pels), polygalacturonases (Pehs), proteases (Prts), and extrapolysaccharides (EPSs), as well as lower swimming and swarming motility and biofilm formation in strain FZ06. Genomic comparison of the two strains revealed five extra gene clusters in FZ06, including one Stt-type T2SS, three T4SSs, and one T4P. Expression of cell wall degrading enzyme (CWDE)-encoding genes is significantly lower in FZ06 than in MS2.

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

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

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

Dickeya zeae, a bacterial plant pathogen of the family Pectobacteriaceae, is responsible for a wide range of diseases on potato, maize, rice, banana, pineapple, taro, and ornamentals and significantly reduces crop production. D. zeae causes the soft rot of taro (Colocasia esculenta) and the heart rot of pineapple (Ananas comosus). In this study, we used Pacific Biosciences single-molecule real-time (SMRT) sequencing to sequence two high-quality complete genomes of novel strains of D. zeae: PL65 (size: 4.74997 MB; depth: 701x; GC: 53.6%) and A5410 (size: 4.7792 MB; depth: 558x; GC: 53.5%) isolated from economically important Hawaiian crops, taro, and pineapple, respectively. Additional complete genomes of D. zeae representing three additional hosts (philodendron, rice, and banana) and other species used for a taxonomic comparison were retrieved from the NCBI GenBank genome database. Genomic analyses indicated the truncated type III and IV secretion systems (T3SS and T4SS) in the taro strain, which only harbored one and two genes of T3SS and T4SS, respectively, and showed high heterogeneity in the type VI secretion system (T6SS). Unlike strain EC1, which was isolated from rice and recently reclassified as D. oryzae, neither the genome PL65 nor A5410 harbors the zeamine biosynthesis gene cluster, which plays a key role in virulence of other Dickeya species. The percentages of average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) between the two genomes were 94.47 and 57.00, respectively. In this study, we compared the major virulence factors [plant cell wall-degrading extracellular enzymes and protease (Prt)] produced by D. zeae strains and evaluated the virulence on taro corms and pineapple leaves. Both strains produced Prts, pectate lyases (Pels), and cellulases but no significant quantitative differences were observed (p > 0.05) between the strains. All the strains produced symptoms on taro corms and pineapple leaves, but the strain PL65 produced symptoms more rapidly than others. Our study highlights the genetic constituents of pathogenicity determinants and genomic heterogeneity that will help to understand the virulence mechanisms and aggressiveness of this plant pathogen.

First Report of Bacterial Soft Rot Disease on Taro Caused by Dickeya fangzhongdai in China

Taro [Colocasia esculenta (L.) Schott.] is an important root crop in the world with great economic value. In recent years, outbreaks of soft rot were observed on taro plants in several plantation areas located in Shaoguan, Guangdong Province, China (25°7'57" N, 113°19'5" E). Root tubers of taro (Paodan variety) infected by soft rot had water-soaked lesions with a dark brown-black margin including a rotten smell, they also had internal rot that was also found in root tubers with no external symptoms. In some areas, the incidence of soft rot can reach up to 30%. To isolate the causal agent, ten pieces of taro root tubers with typical symptoms were surface-sterilized with 75% ethanol and 0.1% HgCl2 solution and then washed thrice with sterile water. The tuber slices were soaked in 50 ml sterile water and shaken at 28°C, 200 rpm for 2 h, and 100 µl was streaked onto the modified Yeast Extract Beef (YEB) agar medium (1% peptone, 0.5% yeast extract, 0.5% sucrose, 0.5% NaCl, 1 Mmol/L MgSO4•7H2O, 1.5% agar, pH 7.0) plates (Zhou et al. 2011) and incubated at 28°C for 24 h. Single colonies grown on YEB were selected for preliminary inoculation onto healthy taro (Paodan variety) slices. Two of the Gram-negative bacteria, named as ZXC1 and MPC2, developed symptoms consistent in rotted decay inside the root tubers after incubation for 24h at 30°C. ZXC1 and MPC2 were biochemically profiled using a Biolog Gen III MicroPlate (Microlog 3, 5.2) (Shen et al. 2019) and resulted Dickeya sp. (SIM 0.856 and 0.704). To determine the species of the Dickeya isolates, 16S rRNA sequences were amplified by primers 27F and 1492R (Hauben et al. 1998). Housekeeping genes including gyrB, atpD, rpoB, and infB were also amplified using degenerate primers (Brady et al. 2008). Results from the BLASTn analysis of the 16S rRNA (GenBank accession numbers MN853405, MN853406), gyrB (GenBank accession numbers MN866299, MN866303), atpD (GenBank accession numbers MN866298, MN866302), rpoB (GenBank accession numbers MN866301, MN866305), and infB (GenBank accession numbers MN866300, MN866304) genes in the isolates ZXC1 and MPC2 showed 99% identities to those of the previously reported D. fangzhongdai isolates from Phalaenopsis (Zhang et al. 2018). Multilocus sequence analysis (MLSA) by MEGA 7.0 performed with four housekeeping genes (gyrB, atpD, rpoB, infB) showed that they clustered with D. fangzhongdai isolates. Analyses using scanning and transmission electron microscopy showed that ZXC1 and MPC2 bacteria were rod-shaped, 0.5-1.0 μm × 1.0-3.0 µm, with peritrichous flagella. Pathogenicity tests were performed thrice using surface-sterilized 2-month-old taro seedlings (Paodan variety). Six individual seedlings were inoculated using a sterile syringe with ten microliters of bacterial suspension (108 CFU/ml) in Tris buffer (0.1 mol/L Tris and 0.1 mol/L HCl, pH 7.4). Taro seedlings injected with sterile Tris buffer were used as the negative control. These taro seedlings were grown in the greenhouse (30 ± 2°C, 90 ± 5% relative humidity). At the 25th day post inoculation, soft rot symptoms were observed in inoculated taro, while all control taro plants remained symptom-free. Small and pale yellow with irregular margins colonies consistent with morphological characteristics of those of D. fangzhongdai were re-isolated from symptomatic taro tubers and the housekeeping genes presence was verified by sequencing as described above, fulfilling Koch's postulates. D. fangzhongdai is a newly emerging bacterial pathogen, which causes bleeding cankers in pear trees (Tian et al. 2016), and soft rot of Phalaenopsis (Zhang et al. 2018). This is the first report of D. fangzhongdai causing soft rot disease in taro. Considering the high incidence of soft rot, this pathogen might pose a significant threat to taro and other economically important crops. Therefore, further researches are needed to investigate host range of the pathogen and develop appropriate integrated management to contain this disease spreading.

Genomic divergence between Dickeya zeae strain EC2 isolated from rice and previously identified strains, suggests a different rice foot rot strain

Rice foot rot caused by Dickeya zeae is an important bacterial disease of rice worldwide. In this study, we identified a new strain EC2 from rice in Guangdong province, China. This strain differed from the previously identified strain from rice in its biochemical characteristics, pathogenicity, and genomic constituents. To explore genomic discrepancies between EC2 and previously identified strains from rice, a complete genome sequence of EC2 was obtained and used for comparative genomic analyses. The complete genome sequence of EC2 is 4,575,125 bp in length. EC2 was phylogenetically closest to previously identified Dickeya strains from rice, but not within their subgroup. In terms of secretion systems, genomic comparisons revealed that EC2 harbored only type I (T1SS), typeⅡ (T2SS), and type VI (T6SS) secretion systems. The flagella cluster of this strain possessed specific genomic characteristics like other D. zeae strains from Guangdong and from rice; within this locus, the genetic diversity among strains from rice was much lower than that of within strains from non-rice hosts. Unlike other strains from rice, EC2 lost the zeamine cluster, but retained the clustered regularly interspaced short palindromic repeats-1 (CRISPR-1) array. Compared to the other D. zeae strains containing both exopolysaccharide (EPS) and capsular polysaccharide (CPS) clusters, EC2 harbored only the CPS cluster, while the other strains from rice carried only the EPS cluster. Furthermore, we found strain MS1 from banana, carrying both EPS and CPS clusters, produced significantly more EPS than the strains from rice, and exhibited different biofilm-associated phenotypes. Comparative genomics analyses suggest EC2 likely evolved through a pathway different from the other D. zeae strains from rice, producing a new type of rice foot rot pathogen. These findings emphasize the emergence of a new type of D. zeae strain causing rice foot rot, an essential step in the early prevention of this rice bacterial disease.