Comparing Protein Expression in Erwinia amylovora Strain TS3128 Cultured under Three Sets of Environmental Conditions.

Erwinia amylovora, the causal agent of fire-blight disease in apple and pear trees, was first isolated in South Korea in 2015. Although numerous studies, including omics analyses, have been conducted on other strains of E. amylovora, studies on South Korean isolates remain limited. In this study, we conducted a comparative proteomic analysis of the strain TS3128, cultured in three media representing different growth conditions. Proteins related to virulence, type III secretion system, and amylovoran production, were more abundant under minimal conditions than in rich conditions. Additionally, various proteins associated with energy production, carbohydrate metabolism, cell wall/membrane/envelope biogenesis, and ion uptake were identified under minimal conditions. The strain TS3128 expresses these proteins to survive in harsh environments. These findings contribute to understanding the cellular mechanisms driving its adaptations to different environmental conditions and provide proteome profiles as reference for future studies on the virulence and adaptation mechanisms of South Korean strains.

Erwinia amylovora, which belongs to the phylum γ-Proteobacteria, is a gram-negative, motile, rod-shaped bacterium with peritrichous flagella (Mansfield et al., 2012). It causes fire blight, one of the most destructive diseases in the family Rosaceae, and poses a global threat to the production of apples and pears in particular (Malnoy et al., 2012). The virulence and biology of E. amylovora strains have been extensively studied worldwide using diverse omics technologies (Mansfield et al., 2012; Piqué et al., 2015). However, studies on E. amylovora strains isolated in South Korea remain limited because the fire blight disease was first reported in the country in 2015 (Myung et al., 2016). Thus, in this study, we investigated the cellular mechanisms driving the response of an E. amylovora strain isolated in South Korea to three different sets of conditions via a comparative proteomic analysis.

We used the virulent South Korean E. amylovora strain TS3128, which was isolated from pear trees in Anseong (Myung et al., 2016). The complete genome sequence of this strain was recently reported (Kang et al., 2021). For comparative proteomic analysis, we used a label-free shotgun proteomic technique, following a previously established protocol (Lee et al., 2021). The strain was cultured in three different liquid media: the Luria Bertani medium (LB; 10 g of tryptone, 5 g of yeast extract, and 10 g of NaCl in 1 l), an hrp inducing medium (HMM; 1 g of (NH4)2SO4, 0.246 g of MgCl2-6H2O, 0.1 g of NaCl, 8.708 g of K2HPO4, and 6.804 g of KH2PO4 in 1 l), and an amylovoran inducing medium (MBMA; 3 g of KH2PO4, 7 g of K2HPO4, 1 g of (NH4)2SO4, 2 ml of glycerol, 0.5 g of citric acid, and 0.03 g of MgSO4, and 10 g of sorbitol in 1 l). The LB medium was used to simulate rich conditions. The HMM and MBMA media were employed to stimulate minimal conditions as represented by two major factors of E. amylovora virulence: the type III secretion system (T3SS) and amylovoran, respectively (Akhlaghi et al., 2020). Nine bacterial cell samples (three replicates per set of conditions) were taken at an OD600 of 0.15 and disrupted via sonication. Total soluble proteins were collected, concentrated, and digested using trypsin to produce tryptic peptides.

Proteomic analysis was performed at the BT Research Facility Center, Chung-Ang University. One microgram of each sample was analyzed via liquid chromatography–tandem mass spectrometry (LC-MS/MS). The samples were injected into a split-free nano LC setup (EASY-nLC II, Thermo Fisher Scientific, Bremen, Germany), connected to an LTQ Velos Pro instrument (Thermo Fisher Scientific) using the nanospray ion mode. Peptides were identified and quantified using an established protocol (Kim et al., 2021). The mass spectra obtained from LC-MS/MS analysis were analyzed using the Proteome Discoverer software (ver. 1.3.0.399, Thermo Fisher Scientific) and SEQUEST algorithm. Whole proteome data on E. amylovora strain TS3128 (accession no. CP056034 for chromosome and CP056034 for plasmid), obtained from the National Center for Biotechnology Information, was used as the database for this analysis. The target-decoy method was employed to increase the confidence level (Elias and Gygi, 2007). The analyzed data were imported into Scaffold 4 (Proteome Software, Portland, OR, USA) for comparison.

It has been reported that E. amylovora strain TS3128 contains 3,237 and 11 predicted protein-coding sequences (CDSs) in chromosomes and plasmids, respectively (Kang et al., 2021). The number of proteins identified in the nine samples during LC-MS/MS ranged between 1,008 and 1,097, corresponding to approximately 31.03–33.77% of total CDSs (Table 1). The number of proteins shared between the three biological replicates of LB, HMM, and MBMA was 969, 998, and 1,033, respectively. These proteins were used for comparative analysis. The normalized average of peptide spectrum matches (PSMs) in the three biological replicates was calculated (Choi et al., 2008) and quantitatively analyzed to identify proteins with counts significantly different from the mean (over two-fold). A student’s t-test was performed using these PSM values. Proteins with a P-value less than 0.05 were considered statistically different.

Table 1

Number of proteins and PSM values of Erwinia amylovora cell samples cultured in LB, HMM, and MBMA media, measured via liquid chromatography-tandem mass spectrometry

Media	1st		2nd		3rd		Proteins shared between the three biological replicates
	Proteins	PSM	Proteins	PSM	Proteins	PSM
LB	1,008	39,942	1,015	39,922	1,020	39,936	969
HMM	1,046	34,000	1,043	33,977	1,046	33,998	998
MBMA	1,096	39,829	1,094	39,810	1,097	39,812	1,033

PSM, peptide spectral match; LB, Luria Bertani; HMM, hrp inducing media; MBMA, amylovoran inducing media.

A total of 91 and 120 proteins were unique to LB and HMM, respectively (Fig. 1A). Additionally, 59 and 100 proteins significantly differed in abundance (over two-fold) in LB and HMM, respectively. These results indicate that the expression of approximately 11.4% of the proteins present in strain TS3128 was altered under those conditions. Comparing LB and MBMA, the expression of 447 proteins was altered. Among these, 96 and 153 proteins were unique, and 82 and 116 proteins were significantly abundant (over two-fold) in LB and MBMA, respectively (Fig. 1B). Lastly, 65 and 90 proteins were unique, and 31 and 37 proteins were significantly abundant (over two-fold) in HMM and MBMA, respectively (Fig. 1C). Notably, the number of proteins that significantly differed in abundance between HMM and MBMA (Fig. 1C) was considerably lower than that between LB and HMM and between LB and MBMA (Fig. 1A and B, respectively). This difference is potentially due to the different conditions of the three media. In contrast to LB, a rich medium, both HMM and MBMA represent minimal conditions, unfavorable for bacterial growth.

Fig. 1

Comparative proteomic analyses of Erwinia amylovora cultured in Luria Bertani (LB), hrp inducing (HMM), and amylovoran inducing (MBMA) media. The Venn diagrams display the numbers of proteins that exhibited significant differences in abundance (over two-fold) between the (A) LB and HMM, (B) LB and MBMA, and (C) HMM and MBMA samples.

The proteins with significantly different levels of expression were subsequently classified using the Clusters of Orthologous Groups (COGs) analysis (Tatusov et al., 2000) to identify the cellular mechanisms driving E. amylovora responses to different conditions. The number of proteins in most categories was higher in HMM than in LB, except for groups J (translation), L (replication/repair), N (transcription), T (signal transduction), and V (defense mechanisms) (Fig. 2A, Supplementary Tables 1 and 2). The HMM medium stimulates hrp gene expression in E. amylovora (Lee and Zhao, 2017). Previous studies have shown that the expression of the hairpin protein (HrpN), a T3SS protein secreted by E. amylovora, is stimulated in hrp-inducing minimal media (Gaudriault et al., 1998; Jin et al., 2001). In another study, rich and minimal media respectively suppressed and stimulated hrp gene expression in E. amylovora (Wei et al., 1992). Reflecting these findings, various T3SS proteins were identified using comparative proteomic analysis. The expression of HrpN (QKZ08390), the type III effector protein HrpK (QKZ08359), Hrp pili protein HrpA (QKZ10963), T3SS outer membrane ring subunit SctC (QKZ08387), and T3SS protein (QKZ08595) was also higher in HMM than LB (Supplementary Table 2). This indicates a successful comparative proteomic analysis of LB and HMM.

Fig. 2

Classification of differentially abundant proteins from the comparative proteomic analysis using clusters of orthologous groups. The bar graphs show numbers of classified proteins from different media (A) Luria Bertani (LB) vs. hrp inducing (HMM), (B) LB vs. amylovora inducing (MBMA), and (C) HMM vs. MBMA. Functional groups: C, energy production and conversion; D, cell cycle control and mitosis; E, amino acid metabolism and transport; F, nucleotide metabolism and transport; G, carbohydrate metabolism and transport; H, coenzyme metabolism; I, lipid metabolism; J, translation; K, transcription; L, replication and repair; M, cell wall/membrane/envelop biogenesis; N, cell motility; O, post-translational modification, protein turnover, chaperone functions; P, inorganic ion transport and metabolism; Q, secondary structure; S, function unknown; T, signal transduction; U, intracellular trafficking and secretion; V, defense mechanisms.

The HMM medium contained 10 mM galactose. Therefore, we hypothesized that galactose metabolism-related proteins would be abundant in E. amylovora cells cultured in HMM. As expected, the expression of several galactose-related proteins was higher in HMM than that in LB (Supplementary Table 2). This held true for MglB (QKZ09819) and MglA (QKZ09818), which transport galactose, galactokinase (QKZ08930), which produces phosphorylated galactose, and multicopper oxidase CueO (QKZ10250) and peroxiredoxin (QKZ08425), which are related to virulence in gram-negative bacteria (Achard et al., 2010; Kaihami et al., 2014). Notably, the number of proteins in group N (motility) was much higher in LB than that in HMM. The expression of T3SS- and flagellar-related proteins is negatively cross-controlled in Pseudomonas aeruginosa (Soscia et al., 2007). Therefore, it is likely that motility-related protein levels were lower in HMM than that in LB because of the high expression of T3SS proteins in HMM. The expression of glycosyl transferase (QKZ09182) and the omptin family outer membrane protease (QKZ11062) was also higher in HMM than that in LB. Glycosyl transferase and bacterial protease are related to virulence and other cellular mechanisms in pathogenic bacteria (Brannon et al., 2015; Li and Wang, 2012). The comparative proteomic analysis of this study reveals that E. amylovora likely secretes higher levels of various proteins associated with virulence in HMM than those in LB.

In addition, 28 proteins related to energy production and conversion were identified in LB and HMM. For example, a formate dehydrogenase accessory sulfurtransferase, FdhD (QKZ07880), was activated in HMM. In a previous study, FdhD was determined to be essential for formate dehydrogenase in Escherichia coli and contribute to oxidative stress tolerance (Iwadate et al., 2017). Moreover, nicotinamide adenine dinucleotide and hydrogen (NADH)-quinone oxidoreductase subunit L (QKZ09880), subunit NuoH (QKZ09884), and quinone oxidoreductase (QKZ10653) were more abundant in HMM than LB. The NADH-quinone oxidoreductase protein complex is the first in the electron transport chain (Spero et al., 2015). In addition, numerous organic and inorganic ATP-binding cassette transporters were detected in LB and HMM (Supplementary Table 2). Such ATP-binding cassette transporters are important factors of virulence, and their expression is likely stimulated by exposure to harsh environments and host conditions (Zeng and Charkowski, 2021). Therefore, E. amylovora likely expressed proteins related to energy production, stress tolerance, and uptake to survive in the harsh, growth-inhibiting conditions of HMM.

The distribution of proteins observed in the COGs comparison of LB and MBMA was similar to that of LB and HMM (Fig. 2B). The synthetic MBMA medium creates conditions conducive to amylovoran biosynthesis (Lee and Zhao, 2017). As expected, five proteins related to amylovoran biosynthesis (QKZ09777, QKZ09779, QKZ09783, QKZ09861, and QKZ08087) were observed exclusively in MBMA (Supplementary Tables 3 and 4). In addition, the transcriptional regulator RcsB (QKZ09861) was more abundant in MBMA than in LB (Supplementary Table 4). In previous studies, RcsB positively regulated amylovoran production and contributed to the virulence of E. amylovora (Ancona et al., 2015; Wang et al., 2009). Proteins in group M (cell wall/membrane/envelope biogenesis) were less abundant in LB than in MBMA. Among these proteins, two TonB-dependent receptors (QKZ08834 and QKZ09398) were detected exclusively in MBMA. Other studies have shown that the TonB-dependent receptor is important for virulence and utilizing iron substrates (Kachadourian et al., 1996). Some outer membrane-related proteins-OmpA (QKZ07876), OmpC (QKZ10993), and OmpW (QKZ11016)-were expressed more or detected exclusively in MBMA. As shown in a previous study, OmpA and OmpC influence the virulence of Gram-negative bacteria (Vila-Farrés et al., 2017) and E. coli (Hejair et al., 2017), respectively. Meanwhile, OmpW has been widely found in Gram-negative bacteria responding to antibiotic stress and environmental stimuli (Xiao et al., 2016). In addition, outer membrane assembly proteins of the AsmA family (QKZ10885, QKZ07930, and QKZ09789) were more abundant in MBMA than in LB. The AmsA protein family is related to membrane biosynthesis and lipopolysaccharide production, a major virulence factor (Deng and Misra, 1996). This likely explains why E. amylovora produced more proteins related to virulence in MBMA than in LB.

In a previous study, sorbitol, a component of MBMA, triggered a high level of amylovoran production in E. amylovora (Bellemann et al., 1994). Sorbitol-6-phosphate dehydrogenase (QKZ10463) can be used by E. amylovora to utilize sorbitol for carbohydrate transport in the host plant and the interconversion of sorbitol to fructose (Salomone-Stagni et al., 2018). Therefore, because the expression of sorbitol-6-phosphate dehydrogenase was higher in MBMA than that in LB, it is likely that sorbitol was used as a carbon source by E. amylovora cells cultured in the former. Additionally, the phosphotransferase transport system (PTS) glucitol/sorbitol transporter subunits IIA (QKZ10464) and IIB (QKZ10465) and glucitol/sorbitol permease IIC (QKZ10466) were expressed more in MBMA than in LB. Many carbohydrates, including sorbitol, are transported into bacterial cells via the PTS (Kotrba et al., 2001). The various proteins related to the transport of sorbitol as a sugar source expressed in MBMA may contribute to the production of amylovoran. Moreover, the universal stress protein UspA (QKZ07979) was more abundant in MBMA than in LB. A previous study reported that UspA is crucial for the virulence of Salmonella (Liu et al., 2007). These results indicate that E. amylovora cultured in MBMA expresses various proteins associated with virulence.

The distribution of categorized proteins in HMM and MBMA differed from that of the other two comparisons (Fig. 2C). The number of proteins expressed in most groups in HMM and MBMA was similar. Notably, more proteins in group I (lipid metabolism) were expressed in HMM than in MBMA. In contrast, the number of proteins in groups J (translation), M (cell wall/membrane/envelope biogenesis), and N (cell motility) in MBMA was considerably higher than those in HMM (Fig. 2C). Overall, the expression of proteins associated with sorbitol metabolism was higher in MBMA, whereas galactose-related proteins were more abundant in HMM (Supplementary Tables 5 and 6). This difference was possibly due to the composition of the media.

In this study, we used an E. amylovora strain isolated from South Korea to investigate the mechanisms driving the expression of different proteins under three different sets of conditions. Many virulence-related proteins, including those involved in T3SS and amylovoran production, were expressed more under minimal conditions than in rich condition (Fig. 3). In addition, various proteins expressed in the three groups of samples were related to carbohydrate metabolism, energy production, ion uptake, transcription, and cell wall/membrane/envelope biogenesis. These results serve as a valuable point of reference for studying the virulence and adaptation mechanisms of South Korean strains of E. amylovora.

Fig. 3

Diagram of putative mechanisms driving Erwinia amylovora responses to minimal and rich conditions. Red and green indicate more abundant proteins in hrp inducing (HMM) and amylovoran inducing (MBMA; minimal) and Luria Bertani (LB; rich) media, respectively. The diagram shows selected proteins from the list.