Quantification of enterohemorrhagic Escherichia coli O157:H7 protein abundance by high-throughput proteome

Abstract

Enterohemorrhagic Escherichia coli (EHEC) O157:H7 is a human pathogen responsible for diarrhea, hemorrhagic colitis and hemolytic uremic syndrome (HUS). To promote a comprehensive insight into the molecular basis of EHEC O157:H7 physiology and pathogenesis, the combined proteome of EHEC O157:H7 strains, Clade 8 and Clade 6 isolated from cattle in Argentina, and the standard EDL933 (clade 3) strain has been analyzed. From shotgun proteomic analysis a total of 2,644 non-redundant proteins of EHEC O157:H7 were identified, which correspond approximately 47% of the predicted proteome of this pathogen. Normalized spectrum abundance factor analysis was performed to estimate the protein abundance. According this analysis, 50 proteins were detected as the most abundant of EHEC O157:H7 proteome. COG analysis showed that the majority of the most abundant proteins are associated with translation processes. A KEGG enrichment analysis revealed that Glycolysis / Gluconeogenesis was the most significant pathway. On the other hand, the less abundant detected proteins are those related to DNA processes, cell respiration and prophage. Among the proteins that composed the Type III Secretion System, the most abundant protein was EspA. Altogether, the results show a subset of important proteins that contribute to physiology and pathogenicity of EHEC O157:H7.

About the Authors:

Wanderson Marques Da Silva

Contributed equally to this work with: Wanderson Marques Da Silva, Jinlong Bei

Roles Methodology, Writing – original draft

Affiliation: Institute of Biotechnology, CICVyA, National Institute of Agricultural Technology, Hurlingham, Buenos Aires, Argentina

Jinlong Bei

Contributed equally to this work with: Wanderson Marques Da Silva, Jinlong Bei

Roles Conceptualization, Methodology, Writing – original draft

Affiliation: AGRO-Biological Gene Research Center, Guangdong Academy of Agricultural Sciences (GDAAS), Guangzhou, China

Natalia Amigo

Roles Formal analysis, Methodology

Affiliation: Institute of Biotechnology, CICVyA, National Institute of Agricultural Technology, Hurlingham, Buenos Aires, Argentina

María Pía Valacco

Roles Conceptualization, Methodology

Affiliation: CEQUIBIEM (Mass Spectrometry Facility), Faculty of Exact and Natural Sciences, University of Buenos Aires and CONICET (National Research Council), Buenos Aires, Argentina

Ariel Amadio

Roles Conceptualization, Writing – original draft

Affiliation: Rafaela Experimental Station, National Institute of Agricultural Technology, Rafaela, Santa Fe, Argentina

Qi Zhang

Roles Investigation, Methodology

Affiliation: AGRO-Biological Gene Research Center, Guangdong Academy of Agricultural Sciences (GDAAS), Guangzhou, China

Xiuju Wu

Roles Investigation, Methodology

Affiliation: AGRO-Biological Gene Research Center, Guangdong Academy of Agricultural Sciences (GDAAS), Guangzhou, China

Ting Yu

Roles Investigation, Methodology

Affiliation: AGRO-Biological Gene Research Center, Guangdong Academy of Agricultural Sciences (GDAAS), Guangzhou, China

Mariano Larzabal

Roles Investigation, Writing – original draft

Affiliation: Institute of Biotechnology, CICVyA, National Institute of Agricultural Technology, Hurlingham, Buenos Aires, Argentina

Zhuang Chen

Roles Conceptualization, Investigation, Methodology

Affiliation: AGRO-Biological Gene Research Center, Guangdong Academy of Agricultural Sciences (GDAAS), Guangzhou, China

Angel Cataldi

Roles Conceptualization, Investigation, Writing – review & editing

* E-mail: [email protected]

Affiliation: Institute of Biotechnology, CICVyA, National Institute of Agricultural Technology, Hurlingham, Buenos Aires, Argentina

ORCID logo http://orcid.org/0000-0002-6744-3477

Introduction

Enterohemorrhagic Escherichia coli (EHEC) O157:H7 is a zoonotic pathogen belonging to Shiga toxin-producing E. coli (STEC) and responsible for different diseases as diarrhea, hemorrhagic colitis and hemolytic uremic syndrome (HUS). HUS is distributed worldwide and considered to be a public health problem in several countries [1,2]. Unfortunately, Argentina is the country with the highest incidence of HUS in the world, with approximately 14 cases per 100,000 in children under 5 and a report of 500 cases per year [3,4]. Cattle are the main reservoir of EHEC. Several studies have shown that most cases related to infection in human may be attributed to the high consumption of foods of bovine origin and especially ground beef is the main source of contamination [5].

Great efforts had been made to characterize strains of E. coli O157:H7 isolated from Argentinian cattle [6]. Using the analysis of simple nucleotide polymorphisms, we have classified 16 strains of STEC O157:H7 in clade 6 and 8, which are the most virulent clades [6]. In vitro and in vivo experimental results showed that the strains Rafaela II (clade 8) and 7.1 Anguil (clade 6) have a high virulence potential when compared with other strains and the standard strain EHEC O157:H7 EDL933 [7]. These results enabled us to characterize the high prevalence of strains clade 6 and 8 in the Argentinian cattle. Importantly, these two clades might contribute to a high incidence of HUS in Argentina.

The availability of whole genome sequences of different EHEC strains has enabled genome-wide comparisons to identify factors that might be correlated to physiology and virulence of this pathogen [8]. In addition, the implementation of system biology approaches, such as prediction of protein-protein network, has contributed substantially in the understanding of the pathogen and interactions with its host [9].

Information about the functions and activities of the individual proteins and pathways that control these systems is essential to understand complex processes occurring in living cells. Large scale quantitative proteomics is a powerful approach used to understand global proteomic dynamics in a cell, tissue or organism, and has been widely used to study protein profiles in the field of microbiology [10]. Furthermore, the study of the abundance of proteins in different conditions or during different stages of growth or disease can provide important information about the activities of individual protein components or protein networks and pathways. The rapid growth of proteomic and genomic methods and tools has managed to reveal the basic protein inventory of a few hundred different organisms. Quantitative proteomic approaches have been applied to determine the absolute or relative abundance of proteins. This information gives insights about the biological function and properties of the cell as well as how cells respond to environmental or metabolic changes or stresses [11, 12]. Quantitative proteomics analysis can contribute to the generation of datasets that are critical for our understanding of global proteins expression and modifications underlying the molecular mechanism of biological processes and disease states.

In a previous study, we reported the use of isobaric tags for comparative quantitation (TMT) method to identify the differentially expressed proteins among three EHEC O157:H7 isolates: Rafaela II (Clade 8), Anguil 7.1 (Clade 6) and EDL933 (Clade 3) [7]. The proteome differences observed among these strains are related mainly to proteins involved in both virulence and cellular metabolism; which might reflect the virulence potential of each strain [7]. The aim of the present study was to promote a more comprehensive insight into the molecular basis of EHEC O157:H7 physiology. For this purpose, we applied high-throughput proteomics to combine the proteome of three EHEC O157:H7 isolates: Rafaela II, Anguil 7.1 and EDL933 and normalized spectrum abundance factor (NSAF) approach [13] to quantify the EHEC O157:H7 proteome.

Material and methods

Bacterial strain and growth conditions

The EHEC O157:H7 strains Rafaela II (clade 8) and 7.1 Anguil (clade 6) isolated from cattle in Argentina and EDL933 (clade 3) strain recovered from a patient in USA were routinely maintained in Luria-Bertani broth (LB, Difco Laboratories, USA) or in LB 1.5% bacteriological agar plates, at 37°C. For the proteomic studies, bacterial strains were cultured as previously described by Amigo et al. [7]. Overnight cultures of the different EHEC O157:H7 strains growth in LB were inoculated (1:50) in Dulbecco’s modified Eagle’s medium (DMEM)-F12 nutrient until reach the mid-exponential growth phase (OD600 nm = 0.6) under a 5% CO2 atmosphere at 37°C.

Protein extraction and preparation of whole bacterial lysates for LC-MS/MS

After bacterial growth, protein extractions were performed according to Amigo et al. [7]. Three biological replicates of each culture were centrifuged at 5000 x g for 20 min at 4°C. The cell pellets were resuspended in ice-cold lysis buffer (50 mM Tris-HCl, pH 7.5, 25 mM NaCl, 5 mM DTT and 1 mM PMSF) and disrupted by three cycles in liquid N2 and subsequently placed in boiling water. The resulting lysates were centrifuged at 30,000 × g for 10 min and precipitated with 5 volumes of ice-cold acetone at -20°C overnight. Next, the protein pellets were resuspended in buffer containing 8 M urea, 2 M thiocarbamide and 200 mM tetraethylammonium bromide at pH 8.5. The protein concentration was determined by the Bradford assay using BSA curve as a standard. Subsequently, the samples were reduced with tris-(2-carboxyethyl)-phosphine (200 mM), alkylated with iodoacetamide (375 mM) and enzymatically digested with sequencing grade trypsin. Finally, the samples were labeled with TMT Reagents 6-plex Kit according to the manufacturer's instructions.

Liquid chromatography and mass spectrometry

The proteomic analyses were performed using High pH Reverse Phase Fractionation and Nano LC-MS/MS Analysis by Orbitrap Fusion. Firstly, the labeled peptides were pooled together and desalted using Sep-Pak SPE (Waters) to remove salt ions. The hpRP chromatography was performed with Dionex UltiMate 3000 model on an Xterra MS C18 column (3.5 um, 2.1 × 150 mm, Waters). The sample were dissolved in buffer A (20 mM ammonium formate, pH 9.5) and then eluted with a gradient of 10 to 45% buffer B (80% acetonitrile (ACN)/20% 20 mM NH4HCO2) for 30 min, followed by 45% to 90% buffer B for 10 min, and a 5-min hold at 90% buffer B. Forty-eight fractions collected at 1 min intervals were merged into 12 fractions. The nano LC MS/MS analysis was carried out using a Orbitrap Fusion tribrid (Thermo-Fisher Scientific, San Jose, CA) mass spectrometer with an UltiMate 3000 RSLC nano system (Thermo-Dionex, Sunnyvale, CA). The fraction was injected onto a PepMap C18 trapping column (5 μm, 200 μm × 1 cm, Dionex) and separated on a PepMap C18 RP nano column (3 μm, 75 μm × 15 cm, Dionex). For all the analysis, the mass spectrometer was operated in positive ion mode, MS spectra were acquired across 350–1550 m/z scan mass range, at a resolution of 12,0000 in the Orbitrap with the max injection time of 50 ms. Tandem mass spectra were recorded in high sensitivity mode (resolution >30000) and made by HCD at normalized collision energy of 40. Each cycle of data-dependent acquisition (DDA) mode selected the top10 most intense peaks for fragmentation. The data were acquired with Xcalibur 2.1 software (Thermo-Fisher Scientific).

Database searching, protein identification and abundance estimation

Tandem mass spectra were extracted and charge state deconvolution and deisotoping were not performed. All MS/MS samples were analyzed using Mascot (Matrix Science, London, UK; version 2.4.1). Mascot was set up to search the EDL933_NCBI_20141031.fasta; TW14539_exclusive_20150310 database (unknown version, 6341 entries) assuming the digestion enzyme trypsin. Mascot was searched with a fragment ion mass tolerance of 0.020 Da and a parent ion tolerance of 8.0 PPM. Carbamidomethyl of cysteine and TMT-6plex of lysine and the n-terminus were specified in Mascot as fixed modifications. Deamidated of asparagine and glutamine and oxidation of methionine were specified in Mascot as variable modifications. Scaffold (version Scaffold_4.8.8, Proteome Software Inc., Portland, OR) was used to validate MS/MS based peptide and protein identifications. Peptide identifications were accepted if they could achieve an FDR less than 1.0% by the Scaffold Local FDR algorithm and contained at least 1 identified peptide. Label free quantification value was calculated by Normalized spectrum abundance factor (NSAF) algorithm [14].

Bioinformatics analysis

Functional annotations were assigned by the COG database [15]. Metabolic pathways were determined by analyzing proteins with the Kyoto Encyclopedia of Genes pathways and Genomes (KEGG) [16].

Results and discussion

Global proteomic analysis and functional classification of Escherichia coli (EHEC) O157:H7 proteome

In this study we have promoted insights into EHEC O157:H7 proteome from a dataset generated with strains E. coli O157:H7 Rafaela II, Anguil 7.1 and EDL933. The strains were grown in D-MEM media and then, proteins from total bacterial lysates were extracted and digested in solution. The resulting peptides were analyzed by 2D-LC MS/MS. From this proteomic analysis, we detected 2,644 non-redundant EHEC O157:H7 proteins (S1 Table). When comparing this result with in silico data of EHEC O157:H7 genome, approximately 47% of the predicted proteome of this pathogen was identified (Fig 1A). To determine the abundance of the identified proteins, the NSAF approach [13] was used. This approach determines to significance of the expression changes based on individual protein intensities (Zybailov et al., 2006) (S1 Table).

[Figure omitted. See PDF.]

Fig 1. Characterization of EHEC O157:H7 proteome and correlation with in silico data.

(A) Correlation of the proteomic results with in silico data of EHEC O157:H7 genome. (B) Dynamic range based on the emPAI value of the proteins identified by LC-MS analysis; pink, most abundant proteins; green, less abundant proteins and red, proteins related to EHEC O157:H7 and that are present in the LEE pathogenicity island.

https://doi.org/10.1371/journal.pone.0208520.g001

According this analysis a dynamic range of protein abundance was generated (Fig 1B). Fifty proteins were identified as most abundant in EHEC O157:H7 proteome (Table 1). Of the total of proteins identified, 25 proteins are encoded by genes that are present in the pO157 plasmid; however, these proteins did not show a high abundance level (S1 Table)

[Figure omitted. See PDF.]

Table 1. List of the most abundant proteins of EHEC O157:H7 proteome.

https://doi.org/10.1371/journal.pone.0208520.t001

We subsequently performed functional annotation of the identified proteins using gene ontology [15]. Cluster of orthologous group (COG) analysis grouped the identified proteins into four important functional groups: (i) metabolism, (ii) information storage and processing, (iii) cellular processes and signaling, and (iv) poorly characterized (Fig 2A). Although most of the identified proteins are related to cellular metabolism, the most abundant proteins are involved in the translation process, followed by energy metabolism and posttranslational modification, protein turnover and chaperones, which shows an intense metabolic activity mainly in the protein synthesis (Fig 2B). On the other hand, most of the less abundant proteins are involved in replication, recombination and repair (Fig 2B).

[Figure omitted. See PDF.]

Fig 2. Functional analysis of the EHEC O157:H7 proteome.

(A) Proteins classified by COG functional categories (B) Categorization of the proteins identified into biological processes. [C] Energy production and conversion; [E] Amino acid transport and metabolism; [D] Cell cycle control, cell division, chromosome partitioning; [G] Carbohydrate transport and metabolism; [I] Lipid transport and metabolism; [J] Translation, ribosomal structure and biogenesis; [K] Transcription; [L] Replication, recombination and repair; [M] Cell wall/membrane/envelope biogenesis; [N] Cell motility; [O] Posttranslational modification, protein turnover, chaperones; [P] Inorganic ion transport and metabolism; [R] General function prediction only; [S] Function unknown; [T] Signal transduction mechanisms; [U] Intracellular trafficking, secretion, and vesicular transport; [V] Defense mechanisms; [X] Mobilome: prophages, transposons. (C) KEGG pathway enrichment analysis, the colors are based on the protein abundance; blue, most abundant and green, less abundant.

https://doi.org/10.1371/journal.pone.0208520.g002

Pieper et al. [17] and Ishihama et al. [18] also conducted proteomic studies on E. coli K-12 and EHEC O157:H7 strain 86–24, respectively, to determine the absolute abundance of proteins. Thirteen proteins of the most abundant proteins in our study were also found as the most abundant proteins in E. coli K-12 (Table 2). Those proteins are related to carbohydrate metabolism, transcription, translation, posttranslational modification and signal transduction mechanisms [18]. On the other hand, only 11 proteins of the most abundant group (Table 2) were the most abundant ones in the data obtained from quantitative proteome of EHEC O157:H7 strain 86–24 [17]. Some of those proteins (e. g. TerD, TerE, EspA and DNA-damage-inducible protein I) are absent from E. coli K-12. Interestingly, when comparing our results with those of Pieper et al. [17] and Ishihama et al. [18], the E. coli proteome was evaluated in different grown condition. Despite the different growth conditions, glyceraldehyde-3-phosphate dehydrogenase, translation elongation factor Tu, DNA-binding protein H-NS, alkyl hydroperoxidereductase protein C, GroEL chaperone and 50S ribosomal protein L7/L12 were detected as the most abundant proteins as well (Table 2). These results suggest a set of proteins that may play an important role in the biology of E. coli.

[Figure omitted. See PDF.]

Table 2. List of the most abundant proteins detected in E. coli K-12 and EHEC 86–24.

https://doi.org/10.1371/journal.pone.0208520.t002

We also detected shiga-toxin subunits such as StxA, StxB, Stx2a and Stx2cb; these proteins, however, were not among the most abundant proteins (Fig 1). Pieper et al. [17] also obtained similar results in EHEC 86–24 proteome. This low abundance can be associated with environmental or nutritional conditions that contribute to the bacterial lysis and consequently to the production of the toxin [17, 19, 20].

Metabolic network analysis

To identify the most relevant biological pathways of the identified proteins, we performed a KEGG enrichment analysis. This analysis provides a comprehensive understanding about pathways that might contribute to cellular physiology [16]. When we evaluated the most abundant proteins, we identified 10 pathways that were considered significant (p < 0.05), among them the Glycolysis / Gluconeogenesis was the most significant (Fig 2C). On the other hand, among the less abundant proteins were proteins related to ABC transport. Different studies have reported that glycolysis / gluconeogenesis pathway might influence in the colonization process of EHEC in the gastrointestinal tract of both mouse and bovine [21, 22]. Although glycolysis substrates inhibit the expression of genes that are localized in locus of enterocyte effacement (LEE), this pathway plays an important role in the initial colonization and maintenance of EHEC in the mouse intestine. In addition, gluconeogenesis not only induces LEE gene expression, but contributes also to the later stages of EHEC colonization in mouse [21, 23].

In our proteomic analysis, 23 proteins that composed the Glycolysis / Gluconeogenesis pathway of E. coli were identified (Fig 3). NAD-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the second most abundant protein of EHEC O157:H7 proteome (Table 1). This important cytoplasmic protein of the Glycolysis pathway is also described as a moonlight protein, owing to the distinct functions performed by this enzyme in different cellular localization [24]. Some studies showed that GAPDH secreted by EHEC and enteropathogenic E. coli (EPEC) strains can bind to fibrinogen and epithelial cell, which could contribute to the pathogenesis of this bacterium mainly through cell adhesion [25, 26]. Another protein that is also described as a moonlight protein and was detected among the most abundant proteins of the EHEC proteome is enolase (Table 1) [27]. This glycolytic enzyme that plays an important role in the carbon metabolism also acts in the RNA degradosome process, mainly in the RNA processing and gene regulation. In E. coli, enolase-RNase E/ degradosome complex regulates bacterial morphology under anaerobic condition by inducing a filamentous form, which is observed by some pathogenic E. coli strains under oxygen limiting conditions [27].

[Figure omitted. See PDF.]

Fig 3. Overview of the glycolysis / gluconeogenesis pathway of EHEC O157:H7.

Enzymes of the Glycolysis / Gluconeogenesis metabolism that were identified at the proteome level. Blue, proteins detected in our proteomic analysis; Green, proteins not identified in our study and Red, proteins detected as most abundant.

https://doi.org/10.1371/journal.pone.0208520.g003

Although, the lipid metabolism not has been detected among the KEEG pathway with significant p-value, the acyl carrier protein (ACP) was detected as the more abundant protein of EHEC O157:H7 proteome. Studies showed that E. coli contain in its genome only a copy of acpP gene, which codify to an ACP, this protein plays key role in the fatty acid biosynthesis and is also required for the growth of E. coli (Rawlings and Cronan, 1992; De Lay and Cronan, 1996). During the fatty acid biosynthesis, ACP is post translationally modified by 4’-phosphopantetheinyl (4’-PP). The acyl intermediates generated are bound to the 4’-PP thiol through a thioester linkage, which allows ACP to transport intermediaries among the fatty acid synthetic enzymes [28, 29].

Information storage and processing

Most proteins described as the most abundant are involved in translation processes. Similar results had been observed in E. coli K-12 [17]. In addition, according to the KEGG enrichment analysis, the ribosome was strongly enriched (Fig 2C). We identified proteins involved in structural elements of the ribosome as well as related to initiation, elongation and terminations steps, which are required to the translation process [30]. These results show an intense metabolic activity of EHEC mainly in protein synthesis. Among these proteins, the translation elongation factor Tu was identified (EF-Tu) (Table 1). EF-Tu could play a role in the resistance process of this bacterium in the gastrointestinal tract [31], as well as against cellular damage generated by the bile salt sodium deoxycholate [32]. Unlike E. coli K-12 [17], the proteins involved in transcription process in EHEC were identified as most abundant. CspA was identified to be among the most abundant proteins as well. This RNA chaperone is described as the major cold shock protein of E. coli. CspA binds to RNA molecules and destabilizes stem loop structures to prevent and resolve misfolding of RNA [33].

Cellular processes and signaling

Flagella are filamentous structures that contribute to pathogenesis of pathogenic E. coli, mainly in motility, adhesion and biofilm production [34]. Generally, this organelle is constituted by basal body, hook and a filament that is composed by flagelin or flagellar antigen FliC, which belongs to the H-antigens group [35, 36]. FliC was detected as highly abundant (Table 1). In addition, a study performed in EPEC showed that FliC might be involved in the inflammatory response during the EPEC infection, due to the capacity of flagelin to induce interleukin-8 (IL-8) release in T84 cells [36].

During infection, E. coli is subject to different environmental conditions, for example, temperature changes that occur both in external ambient and within host. In our proteomic analysis, DnaK, GroEL and GroES were detected among the most abundant proteins (Table 1). Studies have shown that these proteins contribute to the resistance process of EHEC under elevated temperature [37, 38]. In addition, Kudva et al. [39] demonstrated that DnaK and GroEL were induced when EHEC was grown in bovine rumen fluid, thus showing the contribution of these proteins in the adaptation of EHEC to the bovine rumen.

Other type of stress commonly found by EHEC during the infection process is oxidative stress, which is generated by reactive oxygen species (ROS) such as superoxide anion (O2-) hydrogen peroxide (H2O2) and the hydroxyl radical (OH·) produced mainly by host immune response [40]. Thus, to adapt and survive under this stress condition, this bacterium presents different anti-oxidant systems. We detected two members of peroxiredoxins (Prxs) family: periplasmic thiol peroxidase (Tpx) and alkyl hydroperoxide reductase C (AhpC) system (Table 1). These two antioxidant systems play an important role in the scavengers of H2O2 and organic hydroperoxides [41, 42]. Glutaredoxin 3 (Grx3) was also among the most abundant proteins (Table 1). Grx3 is associated with Glutaredoxin (Grx) system, whose function is to reduce disulfide bond in target proteins to control the intracellular redox environment [43]. In addition, Smirnova et al. [44] showed that glutaredoxin proteins might be involved in the resistance of E. coli to antibiotics as ampicillin. Altogether, these different systems promote an efficient pathway of antioxidant defense in EHEC that contributes to the pathogenesis of this bacterium.

The ter operon related to tellurite resistance is widely spread in several Gram positive and Gram negative pathogenic species [45, 46]. In EDL933, this operon is composed by six genes (terZABCDE). Among the proteins expressed by that operon, only TerC was absent from our proteomic analysis. Interestingly, TerD and TerE proteins were among the most abundant proteins of EHEC O157:H7 proteome (Table 1). A study performed with Uropathogenic E. coli (UPEC) isolates showed that the introduction of the ter gene cluster contributes to improve bacterial fitness inside macrophages [47]. On the other hand, Yin et al. [48] demonstrated that ter genes contribute to adherence of EHEC O157:H7 to epithelial cells. However, the true role of these genes in the EHEC pathogenesis remains unclear. Although tellurium is absent from the EHEC niche, interestingly, proteomic studies have detected tellurium resistance proteins in EHEC O157 proteome in different media and growth conditions such as D-MEM [49], minimal medium [50], CHROMagar STEC [51], bovine fluid rumen [39] and under conditions that stimulate the quorum sensing pathway [52]. Despite the several studies in this area, more efforts are necessary to unveil the true role of the tellurium resistance proteins in EHEC pathogenesis.

Locus of Enterocyte Effacement (LEE)

The LEE is a pathogenicity island of 35.6 kb that is organized into five polycistronic operons (LEE1 to LEE5) and is an additional bicistronic operon of glr regulatory proteins [53]. LEE is related to intimate adherence of EHEC to cell host and is required for attaching and effacing (A/E) lesions, followed by the translocation of effector proteins that contribute mainly to host modulation of the immune system [54]. In addition, LEE contains the genes that encode the Type III secretion system (T3SS) as well as some effectors molecules that are exported by this system. The T3SS is responsible for the translocation of effectors from within the host cell, whose are directly involved in the EHEC pathogenesis, mainly in the host modulation of the immune system [54]. In this study the EHEC strains were grown in D-MEM, a medium known to induce expression of genes encoding T3SS [55]. We identified 24 LEE-encoded proteins (Fig 1C, S1 Table). Among these proteins, the most abundant were EspA (filamentous structure of the T3SS), Tir (translocated intimin receptor), EspB (pore formation and effector activity) and EspD (outer membrane adhesin) (Fig 1C).

Interestingly, these proteins play an important role in the E. coli O157 adhesion [56, 57]. On the other hand, EspA, EspB, Tir and Intimin are potential vaccine candidates against EHEC infection [58, 59]. EspA, which was detected as the most abundant protein of LEE, forms a channel that connect the bacterial cytoplasm with the host cell; this exportation conduct allows the translocation of effectors from within the host cell [60]. EspB together with EspD are responsible for the formation of the translocation pore and for the effector translocation of Tir. In addition, EspB can inhibit the interaction between myosin and actin, which promotes loss of microvilli and consequently contributes to the induction of diarrhea [61]. The interaction between Tir and Intimin contributes directly to EHEC O157:H7 persistence during the infection process [62, 63]. Furthermore, Tir and Intimin are involved in the modulation of host immunity. Tir might inhibit tumor necrosis factor receptor-associated factor 6 (TRAF-6)-mediated by NF-κB activation [64]. Instead, intimin can induce a T-helper cell type 1 response as well as to stimulate the proliferation of spleen CD4+ T lymphocytes and cells from lymphoid tissues [65, 66].

Conclusion

In this work, we applied the quantitative proteomic (TMT)-based and emPAI analyses to estimate the quantification of EHEC O157:H7 proteome of combined proteomes of two EHEC O157:H7 isolates from Argentinian cattle and of the standard strainEDL933. These comprehensive proteomic analyses generated a quantitative dataset of EHEC proteome composed of a subset of proteins involved in different biological processes. All these proteins together might form a network of factors that play an important role in the pathogenesis and physiology of this pathogen. Altogether, the results presented in this study provide insights into the functional genome of EHEC O157:H7 at the protein level and could contribute to the understating of the factors associated with the biology of this pathogen.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Total list of proteins identified and quantified by NSAF approach.

https://doi.org/10.1371/journal.pone.0208520.s001

(XLSX)

Acknowledgments

The present study was supported by grants PICT #0211 from FONCYT, Argentina, PNBIO1131034 INTA, Argentina; Presidential Foundation of Guangdong Academy of Agricultural Sciences (No.: 201320), and Science and Technology Program of Guangdong Province, (2016B070701013) China.

WMdS, AA, MPV, ML and AC are CONICET fellows. NA holds a CONICET fellowship. We thank Julia Sabio y García for her critical reading.

Citation: Da Silva WM, Bei J, Amigo N, Valacco MP, Amadio A, Zhang Q, et al. (2018) Quantification of enterohemorrhagic Escherichia coli O157:H7 protein abundance by high-throughput proteome. PLoS ONE 13(12): e0208520. https://doi.org/10.1371/journal.pone.0208520

References

1. Noris M, Remuzzi G. Hemolytic uremic syndrome. J Am Soc Nephrol. 2005;16:1035–1050. pmid:15728781

2. Scallan E, Hoekstra RM, Angulo FJ, Tauxe RV, Widdowson MA, Roy SL et al. Foodborne illness acquired in the United States—major pathogens. Emerg Infect Dis. 201;117:7–15.

3. Ferraris JR, Ramirez JA, Ruiz S, Caletti MG, Vallejo G, Piantanida JJ et al. Shiga toxin-associated hemolytic uremic syndrome: absence of recurrence after renal transplantation. Ped Nephrol. 2002;17:809–814.

4. Rivas M, Padola NL, Lucchesi PMA, Massana M. Diarrheagenic Escherichia coli in Argentina In: Torres AG, editor. Pathogenic Escherichia coli in Latin America. Oak Park (IL), EstadosUnidos: Bentham Science Publishers Ltd.; 2010. pp. 142–61.

5. Etcheverría AI, Padola NL. Shiga toxin-producing Escherichia coli: factors involved in virulence and cattle colonization. Virulence. 2013;1:366–372.

6. Amigo N, Mercado E, Bentancor A, Singh P, Vilte D, Gerhardt E et al. Clade 8 and Clade 6 Strains of Escherichia coli O157:H7 from Cattle in Argentina have Hypervirulent-Like Phenotypes. PLoS One. 2015;10:e0127710. pmid:26030198

7. Amigo N, Zhang Q, Amadio A, Zhang Q, Silva WM, Cui B et al. Overexpressed Proteins in Hypervirulent Clade 8 and Clade 6 Strains of Escherichia coli O157:H7 Compared to E. coli O157:H7 EDL933 Clade 3 Strain. PLoS One. 2016;23:e0166883.

8. Sadiq SM, Hazen TH, Rasko DA, Eppinger M. EHEC Genomics: Past, Present, and Future. Microbiol Spectr. 2014;2:EHEC-0020-2013.

9. Bose T, Venkatesh KV, Mande SS. Computational Analysis of Host-Pathogen Protein Interactions between Humans and Different Strains of Enterohemorrhagic Escherichia coli. Front Cell Infect Microbiol. 2017;19:128.

10. Otto A, Becher D, Schmidt F. Quantitative proteomics in the field of microbiology. Proteomics. 2014;14:547–565. pmid:24376008

11. Otto A, Bernhardt J, Hecker M, Becher D. Global relative and absolute quantitation in microbial proteomics. Curr Opin Microbiol. 2012;15:364–372. pmid:22445110

12. Li H, Han J, Pan J, Liu T, Parker CE, Borchers CH. Current trends in quantitative proteomics—an update. J Mass Spectrom. 2017;52:319–341. pmid:28418607

13. Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP. Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae. J Proteome Res. 2006: 5:2339–2347. pmid:16944946

14. Zhang B, VerBerkmoes NC, Langston MA, Uberbacher E, Hettich RL, Samatova NF. Detecting differential and correlated protein expression in label-free shotgun proteomics J. Proteome Res., 2006; 5:2909–2918.

15. Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43:D261–9. pmid:25428365

16. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res 2012;40: D109–D114. pmid:22080510

17. Pieper R, Zhang Q, Clark DJ, Huang ST, Suh MJ, Braisted JC et al. Characterizing the Escherichia coli O157:H7 Proteome Including Protein Associations with Higher Order Assemblies. PLoS One. 2011;6:e26554. pmid:22087229

18. Ishihama Y, Schmidt T, Rappsilber J, Mann M, Hartl FU, Kerner MJ et al. Protein abundance profiling of the Escherichia coli cytosol. BMC Genomics. 2008;27:102.

19. Harris SM, Yue WF, Olsen SA, Hu J, Means WJ, McCormick RJ et al. Salt at concentrations relevant to meat processing enhances Shiga toxin 2 production in Escherichia coli O157:H7. Int J Food Microbiol. 2012;15:186–192.

20. Mei GY, Tang J, Carey C, Bach S, Kostrzynska M. The effect of oxidative stress on gene expression of Shiga toxin-producing Escherichia coli (STEC) O157:H7 and non-O157 serotypes. Int J Food Microbiol 2015;23:7–15.

21. Miranda RL, Conway T, Leatham MP, Chang DE, Norris WE, Allen JH et al. Glycolytic and gluconeogenic growth of Escherichia coli O157:H7 (EDL933) and E. coli K-12 (MG1655) in the mouse intestine. Infect Immun. 2004;72:1666–1676. pmid:14977974

22. Bertin Y, Deval C, de la Foye A, Masson L, Gannon V, Harel J et et al. The gluconeogenesis pathway is involved in maintenance of enterohaemorrhagic Escherichia coli O157:H7 in bovine intestinal content. PLoS One 2014;2:e98367.

23. Njoroge JW, Nguyen Y, Curtis MM, Moreira CG, Sperandio V. Virulence meets metabolism: Cra and KdpE gene regulation in enterohemorrhagic Escherichia coli. MBio. 2012;16:e00280–12.

24. Henderson B, Martin AC. Protein moonlighting: a new factor in biology and medicine. Biochem Soc Trans. 2014;42:1671–1678. pmid:25399588

25. Egea L, Aguilera L, Giménez R, Sorolla MA, Aguilar J, Badía J et al. Role of secreted glyceraldehyde-3-phosphate dehydrogenase in the infection mechanism of enterohemorrhagic and enteropathogenic Escherichia coli: interaction of the extracellular enzyme with human plasminogen and fibrinogen. Int J Biochem Cell Biol. 2007; 39:1190–1203. pmid:17449317

26. Aguilera L, Giménez R, Badia J, Aguilar J, Baldoma L. NAD+-dependent post-translational modification of Escherichia coli glyceraldehyde-3-phosphate dehydrogenase. Int Microbiol. 2009;12:187–192. pmid:19784925

27. Murashko ON, Lin-Chao S. Escherichia coli responds to environmental changes using enolasicdegradosomes and stabilized DicF sRNA to alter cellular morphology. Proc Natl Acad Sci U S A. 2017;19:E8025–E8034.

28. Rawlings M, Cronan JE Jr. The gene encoding Escherichia coli acyl carrier protein lies within a cluster of fatty acid biosynthetic genes. J Biol Chem. 1992; 25:5751–5754.

29. De Lay NR1, Cronan JE. Gene-specific random mutagenesis of Escherichia coli in vivo: isolation of temperature-sensitive mutations in the acyl carrier protein of fatty acid synthesis. J Bacteriol. 2006; 188:287–96. pmid:16352845

30. Gingold H, Pilpel Y. Determinants of translation efficiency and accuracy. Mol Syst Biol. 2011;12:481.

31. Alpert C, Scheel J, Engst W, Loh G, Blaut M. Adaptation of protein expression by Escherichia coli in the gastrointestinal tract of gnotobiotic mice. Environ Microbiol. 2009;11:751–761. pmid:19175791

32. Ribeiro CB, Sobral MG, Tanaka CL, Dallagassa CB, Picheth G, Rego FG et al. Proteins differentially expressed by Shiga toxin-producing Escherichia coli strain M03 due to the biliar salt sodium deoxycholate. Genet Mol Res. 2013;24:4909–4917.

33. Jiang W, Hou Y, Inouye M. CspA, the major cold-shock protein of Escherichia coli, is an RNA chaperone. J Biol Chem. 1997;3:196–202.

34. Kakkanat A, Phan MD, Lo AW, Beatson SA, Schembri MA. Novel genes associated with enhanced motility of Escherichia coli ST131. PLoS One. 2017;10:e0176290.

35. Fratamico PM, DebRoy C, Liu Y, Needleman DS, Baranzoni GM, Feng P. Advances in Molecular Serotyping and Subtyping of Escherichia coli. Front Microbiol. 2016;3:644.

36. Zhou X, Girón JA, Torres AG, Crawford JA, Negrete E, Vogel SN et al. Flagellin of enteropathogenic Escherichia coli stimulates interleukin-8 production in T84 cells. Infect Immun 2003;71:2120–2129. pmid:12654834

37. Carruthers MD, Minion C. Transcriptome analysis of Escherichia coli O157:H7 EDL933 during heat shock. FEMS Microbiol Lett 2009; 295:96–102. pmid:19473256

38. Singh R, Jiang X. Expression of stress and virulence genes in Escherichia coli O157:H7 heat shocked in fresh dairy compost. J Food Prot. 2015; 78:31–41. pmid:25581175

39. Kudva IT, Stanton TB, Lippolis JD. The Escherichia coli O157:H7 bovine rumen fluid proteome reflects adaptive bacterial responses. BMC Microbiol. 2014;21;14:48.

40. Vidovic S, Korber DR. Escherichia coli O157: Insights into the adaptive stress physiology and the influence of stressors on epidemiology and ecology of this human pathogen. Crit Rev Microbiol. 2016; 42:83–93. pmid:24601836

41. Cha MK, Kim HK, Kim IH. Mutation and Mutagenesis of thiol peroxidase of Escherichia coli and a new type of thiol peroxidase family. J Bacteriol 1996; 178:5610–5614. pmid:8824604

42. Seaver LC, Imlay JA. Alkyl hydroperoxide reductase is the primary scavenger of endogenous hydrogen peroxide in Escherichia coli. J Bacteriol. 2001; 183:7173–7181. pmid:11717276

43. Meyer Y, Buchanan BB, Vignols F, Reichheld JP. Thioredoxins and glutaredoxins: unifying elements in redox biology. Annu Rev Gene. 2009; 43:335–367.

44. Smirnova G, Muzyka N, Lepekhina E, Oktyabrsky O. Roles of the glutathione- and thioredoxin-dependent systems in the Escherichia coli responses to ciprofloxacin and ampicillin. Arch Microbiol. 2016; 198:913–921. pmid:27277520

45. Taylor DE. Bacterial tellurite resistance. Trends Microbiol 1999; 7:111–115. pmid:10203839

46. Kormutakova R, Klucar L, Turna J. DNA sequence analysis of the tellurite-resistance determinant from clinical strain of Escherichia coli and identification of essential genes. Biometals 2000; 13:135–159. pmid:11016400

47. Valkova D, Valkovicova L, Vavrova S, Kovacova E, Mravec J, Turna J. The contribution of tellurite resistance genes to the fitness of Escherichia coli uropathogenic strains. Cent Eur J Biol. 2007; 2:182–191.

48. Yin X, Wheatcroft R, Chambers JR, Liu B, Zhu J, Gyles CL. Contributions of O island 48 to adherence of enterohemorrhagic Escherichia coli O157:H7 to epithelial cells in vitro and in ligated pig ileal loops. Appl Environ Microbiol. 2009;75:5779–5786 pmid:19633120

49. Kudva IT, Griffin RW, Krastins B, Sarracino DA, Calderwood SB, John M. Proteins other than the locus of enterocyte effacement-encoded proteins contribute to Escherichia coli O157:H7 adherence to bovine rectoanal junction stratified squamous epithelial cells. BMC Microbiol. 2012; 12:103. pmid:22691138

50. Islam N, Nagy A, Garrett WM, Shelton D, Cooper B, Nou X. Different Cellular Origins and Functions of Extracellular Proteins from Escherichia coli O157:H7 and O104:H4 as Determined by Comparative Proteomic Analysis. Appl Environ Microbiol. 2016; 30:4371–4378.

51. Kalule JB, Fortuin S, Calder B, Robberts L, Keddy KH, Nel AJM et al. Proteomic comparison of three clinical diarrhoeagenic drug-resistant Escherichia coli isolates grown on CHROMagar STEC media. J Proteomics. 2017;5:S1874–3919.

52. Kim Y, Oh S, Ahn EY, Imm JY, Oh S, Park S et al. Proteome Analysis of Virulence Factor Regulated by Autoinducer-2–like Activity in Escherichia coli O157:H7. J Food Prot. 2007; 70:300–307. pmid:17340862

53. Elliott SJ, Wainwright LA, McDaniel TK, Jarvis KG, Deng YK, Lai LC et al. The complete sequence of the locus of enterocyte effacement (LEE) from enteropathogenic Escherichia coli E2348/69. Mol Microbiol. 1998; 28:1–4. pmid:9593291

54. Stevens MP, Frankel GM. The Locus of Enterocyte Effacement and Associated Virulence Factors of Enterohemorrhagic Escherichia coli. Microbiol Spectr 2014; 2:EHEC-0007-2013.

55. Tobe T, Beatson SA, Taniguchi H, Abe H, Bailey CM, Fivian A et al. An extensive repertoire of type III secretion effectors in Escherichia coli O157 and the role of lambdoid phages in their dissemination. Proc Natl Acad Sci U S A. 2006; 3:14941–14946.

56. Naylor SW, Roe AJ, Nart P, Spears K, Smith DG, Low JC et al. Escherichia coli O157: H7 forms attaching and effacing lesions at the terminal rectum of cattle and colonization requires the LEE4 operon. Microbiology 2005; 151:2773–2781. pmid:16079353

57. Sharma VK, Kudva IT, Bearson BL, Stasko JA. Contributions of EspA Filaments and Curli Fimbriae in Cellular Adherence and Biofilm Formation of Enterohemorrhagic Escherichia coli O157:H7. PLoS One. 2016; 22:e0149745.

58. Li Y, Frey E, Mackenzie AM, Finlay BB. Human response to Escherichia coli O157:H7 infection: antibodies to secreted virulence factors. Infect Immun 2000; 68:5090–5095. pmid:10948130

59. Vilte DA, Larzábal M, Garbaccio S, Gammella M, Rabinovitz BC, Elizondo AM et al. Reduced faecal shedding of Escherichia coli O157:H7 in cattle following systemic vaccination with γ-intimin C₂₈₀ and EspB proteins. Vaccine. 2011; 23:3962–3968.

60. Ebel F, Podzadel T, Rohde M, Kresse AU, Krämer S, Deibel C et al. Initial binding of Shiga toxin-producing Escherichia coli to host cells and subsequent induction of actin rearrangements depend on filamentous EspA-containing surface appendages. Mol Microbiol. 1998; 30:147–161. pmid:9786192

61. Iizumi Y, Sagara H, Kabe Y, Azuma M, Kume K, Ogawa M et al. The enteropathogenic E. coli effector EspB facilitates microvillus effacing and antiphagocytosis by inhibiting myosin function. Cell Host Microbe. 2007; 13:383–392.

62. Cornick NA, Booher SL, Moon HW. Intimin facilitates colonization by Escherichia coli O157:H7 in adult ruminants. Infect Immun. 2002; 70:2704–2707. pmid:11953416

63. Dean-Nystrom EA, Bosworth BT, Moon HW, O'Brien AD. Escherichia coli O157:H7 requires intimin for enteropathogenicity in calves. Infect Immun. 1998; 66:4560–4563. pmid:9712821

64. Yan D, Quan H, Wang L, Liu F, Liu H, Chen J et al. Enteropathogenic Escherichia coliTir recruits cellular SHP-2 through ITIM motifs to suppress host immune response. Cell Signal. 2013; 25:1887–1894. pmid:23707390

65. Higgins LM, Frankel G, Connerton I, Gonçalves NS, Dougan G, MacDonald TT. Role of bacterial intimin in colonic hyperplasia and inflammation. Science. 1999;285:588–591. pmid:10417389

66. Gonçalves NS, Hale C, Dougan G, Frankel G, MacDonald TT. Binding of intimin from enteropathogenic Escherichia coli to lymphocytes and its functional consequences. Infect Immun. 2003; 71:2960–2965. pmid:12704179