Transcriptomic Analysis of Shiga Toxin-Producing Escherichia coli FORC_035 Reveals the Essential Role of Iron Acquisition for Survival in Canola Sprouts and Water Dropwort

Enterohemorrhagic Escherichia coli (EHEC) is a foodborne pathogen that poses a serious threat to humans. Although EHEC is problematic mainly in food products containing meat, recent studies have revealed that many EHEC-associated foodborne outbreaks were attributable to spoiled produce such as sprouts and green leafy vegetables. To understand how EHEC adapts to the environment in fresh produce, we exposed the EHEC isolate FORC_035 to canola spouts (Brassica napus) and water dropwort (Oenanthe javanica) and profiled the transcriptome of this pathogen at 1 and 3 h after incubation with the plant materials. Transcriptome analysis revealed that the expression of genes associated with iron uptake were down-regulated during adaptation to plant tissues. A mutant strain lacking entB, presumably defective in enterobactin biosynthesis, had growth defects in co-culture with water dropwort, and the defective phenotype was complemented by the addition of ferric ion. Furthermore, gallium treatment to block iron uptake inhibited bacterial growth on water dropwort and also hampered biofilm formation. Taken together, these results indicate that iron uptake is essential for the fitness of EHEC in plants and that gallium can be used to prevent the growth of this pathogen in fresh produce.


INTRODUCTION
Enterohemorrhagic Escherichia coli (EHEC) is a virulent foodborne pathogen that causes outbreaks associated with serious health problems worldwide (Sperandio and Nguyen, 2012). It causes bloody diarrhea, abdominal cramps, and hemolytic uremic syndrome in human, which may be lifethreatening, especially in young children and the elderly (Karmali et al., 1983;Karmali, 1989;Griffin and Tauxe, 1991). EHEC is commonly found in animal-based foods made with beef, pork, or poultry (Meng and Doyle, 1997). However, fresh produce (e.g., lettuce, spinach, alfalfa sprouts) has been identified as another notable source of EHEC outbreaks (Taormina et al., 1999;Sivapalasingam et al., 2004). The increase in number of outbreaks associated with fresh produce has been attributed in part to increased consumption of fruits and vegetables in recent decades (Lynch et al., 2009). Furthermore, fresh produce is generally consumed without thermal treatment, which is likely to aggravate the occurrence of produce-associated foodborne diseases (Herman et al., 2015).
A number of studies have shown that many foodborne pathogens can survive and replicate in plants, taking advantage of plant-derived amino acids and carbon compounds as nutrients (Cooley et al., 2003;Jablasone et al., 2005;Schikora et al., 2008;Deering et al., 2012). The interactions between bacterial pathogens and their host plants in the context of unfavorable plant conditions including desiccation, UV radiation, and plant immune responses have also been investigated extensively (Shi et al., 2007;Erickson, 2012;Fink et al., 2012). Interestingly, bacterial determinants required for persistence in plants are variable depending on bacterial genomic features and reservoir host plants (O'brien and Lindow, 1989;Heaton and Jones, 2008;Junker et al., 2011). Therefore, to decipher the mechanisms of persistence and replication of foodborne pathogens in plants, the behavior of diverse bacterial pathogens in contact with a variety of plants should be documented. From a food safety and hygiene perspective, although any produce has the potential to carry foodborne pathogens in the field or during postharvest management, it is critical to understand bacterial survival strategies in fresh produce, which is consumed without thermal treatments.
Recent studies of plant-microbe interactions have been biased toward lettuce, spinach, and tomato, which are associated with foodborne disease outbreaks (Carey et al., 2009;Kyle et al., 2010). To extend our understanding of the interactions between foodborne bacteria and plants, we investigated canola sprouts (Brassica napus) and water dropwort (Oenanthe javanica) in this study. Sprouted seeds of plants like canola are cultivated using warm water and therefore provide favorable niches for growth of pathogens such as Salmonella enterica and E. coli O157:H7 (Andrews et al., 1982;Taormina et al., 1999;Erickson, 2012). A multistate outbreak of EHEC occurred in the United States in 1997 because of alfalfa sprout contamination (Breuer et al., 2001). Water dropworts are members of the genus Oenanthe in the plant family Apiaceae and the species O. javanica is mostly consumed raw in eastern Asia. Water dropwort is usually cultivated in flooded soil and is prone to contamination with sewage containing bacterial pathogens.
Transcriptome analysis can lead to a systematic understanding of bacterial responses to environmental changes such as contact with plant tissues (Schenk et al., 2012). In this study, to gain insight into the interaction between fresh produce and EHEC, we analyzed the transcriptome of EHEC FORC_035 when exposed to canola sprouts and water dropwort using an RNA-Sequencing (RNA-Seq) approach, as this provides higher resolution data than microarrays or related techniques (Matkovich et al., 2010). FORC_035, a strain isolated from kimchi in South Korea, possesses multiple genes encoding virulence factors including Stx2 and hemolysin. It causes attaching and effacing (A/E) lesions, and has the potential to cause severe diseases in humans when consumed in the form of contaminated fresh produce. Comprehensive transcriptome analysis of FORC_035 in contact with plants provided informative data for unraveling the interaction between pathogens and plants and revealed that iron acquisition was critical for EHEC survival in plants. Iron is an essential element for all living organisms, including bacteria. Therefore, bacteria possess a variety of complex systems for iron acquisition and utilization (Doherty, 2007). Among multiple iron uptake systems, siderophore-mediated iron acquisition is the most common form in bacteria (Visca et al., 2002). E. coli has been reported to utilize diverse forms of siderophores such as ferrichrome, enterobactin, yersiniabactin, and catecholate that are transported through highly specific proteins encoded by fhu operon, fep operon, ybt operon, and fiu/cir/iroN, respectively (Braun, 2003). The role of iron in bacterial survival and growth has been extensively investigated in bacterial pathogens that encounter iron starvation stress upon entering their hosts (Skaar, 2010). As a critical growth determinant in plants, we further investigated the role of iron in E. coli adaptation to plant environments.

Bacterial Strains and Culture Conditions
The EHEC isolate strain FORC_035 and isogenic mutant derivatives used in this study are described in Table 1. FORC_035 is EHEC (O8:H30) producing Stx2 and hemolysin E. FORC_035 strain was isolated from kimchi by the Incheon Institute of Health and Environment (Incheon, South Korea). The genome sequence of FORC_035 was deposited in the database of the Foodborne pathogen Omics Research Center 1 . Unless stated, all bacterial strains were grown in Luria-Bertani (LB) medium at 37 • C, and the following antibiotics were added if needed: kanamycin (25 µg/mL) and ampicillin (100 µg/mL). The selective medium of tryptone bile X-glucuronide (TBX, Oxoid, Cambridge, United Kingdom) was used to enumerate FORC_035 populations. As a nutrient-limitation condition, M9 minimal medium supplemented with 0.4% glucose was used. Fe 2 (SO 4 ) 3 and Ga(NO 3 ) 3 (Sigma-Aldrich, St. Louis, MO, United States) at pH 7.0 in sterile water were added to bacterial cultures to examine the effect of iron on bacterial growth.
To study the effect of two siderophores, enterobactin and yersiniabactin, entB (deletion mutant of entB encoding enterobactin synthase component B) and ybtS (deletion mutant of ybtS encoding yersiniabactin biosynthesis salicylate synthase) were constructed. An isogenic mutant of FORC_035, entB, was constructed using the λ Red recombination system. PCR products containing a kanamycin resistance gene were generated using pKD13 as a template (Datsenko and Wanner, 2000;Murphy and Campellone, 2003). For homologous recombination, PCR primers were designed to introduce nucleotides identical to the flanking regions of entB at each end of the PCR products (Supplementary Table S1). FORC_035 strain harboring the Red helper plasmid pKD46 was grown in the presence of L-arabinose (50 mM) at 30 • C with shaking (220 rpm) until the OD 600 reached 0.6. Bacterial cells were then centrifuged (10,000 × g, 5 min, 4 • C), resuspended in ice-cold water, and transformed by electroporation (Bio-Rad Laboratories, Hercules, CA, United States). Recombinant cells were selected by kanamycin resistance (Jandu et al., 2009), and subjected to diagnostic PCR using the primers listed in Supplementary Table S1. The kanamycin resistance cassette was subsequently removed using the FLP recombinase of pCP20 (Cherepanov and Wackernagel, 1995). A ybtS mutant strain was also constructed in the same way using ybtS-Lambda-F/ybtS-Lambda-R primer set (Supplementary Table S1).

Preparation of Canola Sprouts and Water Dropwort
Canola sprouts (Brassica napus) and water dropwort (Oenanthe javanica) were purchased from commercial markets in Seoul, South Korea. To prepare water dropwort, only stems around 15 cm in length were cut into pieces, because these are more likely to be contaminated by bacteria than other parts of the plant. Five grams of each plant were washed using 50 mL of sterile water three times with agitation (220 rmp) for 10 min each time prior to contact with bacteria.

In planta Assay
To prepare bacterial total RNA after contact with plants, the FORC_035 strain was cultured in LB until OD 600 reached 1.0, then bacterial cells were washed twice using M9 minimal medium and resuspended in fresh M9 minimal medium (15 mL) at 10 9 CFU/mL. Five-gram of canola sprouts or water dropwort was added to the bacterial culture and cells were incubated at 37 • C with shaking at 220 rpm. The bacterial culture grown in the same condition without adding plants in 20 mL M9 minimal medium was used as a control. At 1-and 3-h post-inoculation, total RNAs were prepared in biological duplicate.

Total RNA Extraction, Sequencing, and Analysis
One hundred microliter of bacterial cultures was treated with RNAprotect bacterial reagent (Qiagen, Hilden, Germany) to quench RNA degradation, and bacterial total RNAs were isolated using an RNeasy mini kit (Qiagen) according to the manufacturer's instructions. RNA samples were treated with Ambion Turbo DNA-free TM (Ambion, Austin, TX, United States) to remove all genomic DNA. The quantity and quality of total RNA samples were examined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States), and RNA samples with RNA integrity numbers (RIN) larger than 9 were used for further analysis (Schroeder et al., 2006). Total RNA samples were stored at −81 • C until use. Five micrograms of total RNA from each sample was used as a starting material and subjected to an rRNA-removal process based on the subtractive hybridization/bead capture system of the Ribo-Zero kit (Epicentre Biotechnologies, Madison, WI, United States). Purified RNA samples were used for mRNA-Seq library construction using the Illumina TruSeq RNA Sample Preparation kit v.2 (Illumina, San Diego, CA, United States). RNA-Seq was performed by two runs of Illumina Hiseq to generate single-end-reads around 100 bp in length. All RNA-Seq data analyzed in this study, including whole transcriptome profiles, are provided as supplementary information files (Supplementary Tables S2, S3) Table S4). Because the coefficient of variation (CV) value for the TMM method was lower than those of RPKM and RLE, TMM was used for the normalization of expression level of genes. The p-value was calculated using edgeR and the [fold change] value was calculated as [TMM Freshproduce /TMM Control ]. For further experiments, differentially expressed genes (DEGs) with an absolute log 2 [fold change] larger than 2 were filtered and visualized using the CLRNAseq program (Chunlab, Seoul, South Korea). Clusters of orthologous groups (COG) analysis (Tatusov et al., 1997) was used for functional grouping of all genes of EHEC FORC_035. The proportion of DEGs in each functional group was calculated.

Heat Map Generation
Heat maps were drawn to graphically visualize global response patterns. Log 2 [fold change] values were used to determine relative gene expression between FORC_035 only and FORC_035 co-cultured with canola sprouts and water dropwort. Heat maps and hierarchical clusters were then generated using Gitools v2.2.2 (Perez-Llamas and Lopez-Bigas, 2011).

Quantitative Real-Time RT-PCR (qRT-PCR)
The expression level of specific genes was validated using qRT-PCR as described in Yoo et al. (2017). For all results of qRT-PCR, GAPDH was used as normalization control. The primers used for the detection of genes are listed in Supplementary Table S5.

Construction of pEntB Plasmid
To construct pEntB plasmid expressing entB gene, pUHE21-2lacI q plasmid was used (Soncini et al., 1995). The entB gene was amplified by PCR using the primer set of entB-F-pUHE and entB-R-pUHE (Supplementary Table S1). The purified PCR product was inserted between the BamHI and HindIII sites of the vector.

Bacterial Cell Count Assay for Growth Comparison
Bacterial cells cultured overnight were washed with M9 minimal media and added in 15 mL of fresh M9 minimal media at a 1:100 ratio with or without 5 g of fresh produce. Bacterial viability was examined by plating on TBX selective agar and represented as log(N/N 0 ) (CFU/mL) values, where N 0 is the initial number of cells at 0 h and N is the number of cells after 5 h incubation. Ferric sulfate or gallium nitrate at indicated concentrations was added to M9 minimal media when required. For strains harboring empty plasmid or pEntB plasmid, 50 µM of IPTG was used for induction. All experiments were performed in biological triplicate. P-value was calculated by Student t-test.

Biofilm Reduction Assay
Bacterial cells cultivated in LB medium overnight were diluted into fresh LB medium (200 µL) at a 1:100 ratio in a polystyrene 96-well plate (SPL, Seoul, South Korea), and grown at 37 • C for 48 h without shaking. Ga(NO 3 ) 3 dissolved in sterile water was added at different concentrations when needed. Considering noobservable-adverse-effect level (NOAEL) and acceptable-dailyintake (ADI) values of gallium, Ga(NO 3 ) 3 was used up to 1 mM (Gómez et al., 1992). Biofilm formation was quantified by staining with 0.1% crystal violet for 5 min followed by washing with phosphate buffer (0.1 M, pH 7.4). All dye attached to the biofilm was dissolved with 200 µL of 33% glacial acetic acid, and OD 570 was measured to quantify the total biofilm mass (Vikram et al., 2010). The assay was performed in biological triplicate, and the p-value was calculated by Student t-test.

Comprehensive Transcriptome Analysis of FORC_035 Cultivated With Canola Sprouts or Water Dropwort
Bacterial transcriptome analysis is an effective approach to understand the interactions between foodborne pathogens and plants (Kyle et al., 2010;Fink et al., 2012;Landstorfer et al., 2014). In this study, the transcriptome of FORC_035 grown with canola sprouts or water dropwort was analyzed using RNA-Seq technology. Bacterial growth was increased in the presence of plant tissues (Figure 1), which was consistent with a previous study showing that bacteria could take advantage of plantsderived nutrients for their growth (Cooley et al., 2003;Jablasone et al., 2005;Schikora et al., 2008;Deering et al., 2012). To understand the bacterial response and adaptations in contact with plants, FORC_035 cells were harvested at 1 and 3 h post-inoculation and their RNAs were isolated for RNA-Seq analysis. RNA-Seq data were acquired, mapped, and normalized as described in the methods (Supplementary Tables S2-S4). After mRNA expression levels were compared between plantexposed samples and not-exposed samples, genes with a log 2 (fold change) greater than 2 or less than −2 were selected and considered to be DEGs (Supplementary Figure S1).
To delineate the physiological changes that occurred in the bacterial cells during contact with plants, 4,977 genes of FORC_035 were categorized according to clusters of orthologous FIGURE 1 | Growth of FORC_035 in M9 minimal media with or without fresh produce. Bacterial cells were inoculated at 10 9 CFU/mL to fresh M9 minimal media containing canola sprouts or water dropwort, respectively. As a control, FORC_035 was cultivated in M9 minimal medium without plant material. In order to discriminate E. coli among indigenous bacteria of plant tissues, culture suspension was diluted and spread onto E. coli selective TBX agar (dotted line) and LB agar (solid line), simultaneously. When only plants were incubated, growth of E. coli was not detected (Supplementary Figure S2). Each symbol indicates the mean value from triplicate measurements.
Frontiers in Microbiology | www.frontiersin.org Frontiers in Microbiology | www.frontiersin.org groups (COG) (Tatusov et al., 1997) and the proportion of DEGs in each functional group was calculated (Figure 2). In FORC_035 co-cultured with canola sprouts, the expression of 8.14% (405/4,977) and 12.54% (624/4,977) genes changed significantly at 1 and 3 h, respectively. In cells exposed to water dropwort, 5.65% (281/4,977) and 5.67% (282/4,977) of genes exhibited significant differences in their expression at 1 and 3 h, respectively. In both of canola sprouts and water dropwort, co-cultivation for 1 h induced significant expression changes in genes associated with translation, ribosomal structure, and biogenesis (30.56% for canola sprouts and 19.44% for water dropwort), cell motility (20.79% for canola sprouts and 10.89% for water dropwort), energy production and conversion (14.98% for canola sprouts and 10.10% for water dropwort), and carbohydrate transport and metabolism (12.70% for canola sprouts and 9.26% for water dropwort). Meanwhile at 3 h postinoculation, groups of genes related to secondary metabolites biosynthesis, transport and catabolism (42.55% for canola sprouts and 10.64% for water dropwort), nucleotide transport and metabolism (31.00% for canola sprouts and 11.00% for water dropwort), and energy production and conversion (27.36% for canola sprouts and 19.54% for water dropwort) showed significant differences in their expression in common. Functional groups of genes that showed significant expression changes upon contact with both plants were graphically compared at 1 and 3 h based on COG (Figures 2C,D).

RNA-Seq Analysis Identifies Genes Differentially Expressed in Response to Canola Sprouts or Water Dropwort
Based on their predicted functions, genes up-or down-regulated as a cluster or an operon in response to plant tissues were further defined, because these genes might be involved in bacterial adaptation to plant environments. Genes specifically regulated by exposure to fresh produce and their functions are presented in Table 2. Representative genes from each group were validated by qRT-PCR (Supplementary Figure S3).

(i) Adhesion
In this study, genes encoding various types of fimbriae, including curli and type 1 fimbriae, were down-regulated upon contact with canola sprouts and water dropwort ( Table 2). Bacteria switch between planktonic and sessile states depending on environmental conditions. When they encounter stresses such as nutrient limitation, they adhere to surfaces and develop extracellular polymeric substances (EPSs) to form a bacterial aggregation called a biofilm (Oh et al., 2007), which enables them to resist various stresses and acquire scarce nutrients (Enos and Taylor, 1996). Bacteria utilize diverse cell surface appendages, including fimbriae, lipopolysaccharides, O-antigen capsule, flagellae, and non-fimbrial adhesions for initial attachment to surfaces (Kierek-Pearson and Karatan, 2005). For example, enteric pathogens such as E. coli and Salmonella enterica possess several different kinds of fimbrial structures and exploit them to colonize salad and alfalfa sprouts (Barak et al., 2005;Fink et al., 2012). In contrast, bacteria favor a planktonic lifestyle under nutrient-rich conditions (Marshall et al., 1971;Brown et al., 1977). Therefore, consistent with the increased growth of FORC_035 in the presence of plants (Figure 1), the downregulation of fimbriae genes in response to plants implies that free-living cells are preferred due to the abundance of nutrients exuded from plant tissues.

(ii) Amino Acid and Nitrogen Metabolism
In accordance with the profound expression changes of genes associated with energy production and conversion (Figure 2), numerous genes important for amino acid and nitrogen metabolism showed significantly altered expression upon contact with plants, intimating the efflux of nutrients from co-inoculated plant tissues.
For example, genes in the ast and puu operons, which are required for arginine and proline metabolism, respectively, tended to be induced during adaptation in plant-supplemented environments ( Table 2): the log 2 [fold change] was higher at 3 h than at 1 h in both plants. Arginine succinyltransferase (AST) pathway comprises an ast operon that allows utilization of arginine as a nitrogen source in E. coli (Schneider et al., 1998). Genes in the puu operon encode enzymes required for degradation of putrescine, a major polyamine molecule involved in a variety of biological processes, including growth rate control (Tweeddale et al., 1998). The induction of genes involved in arginine and putrescine degradation during bacterial adaptation to plant tissues might be indicative of abundant nitrogen during early exposure, but the depletion of preferable nitrogen sources due to bacterial overgrowth in the medium at later time points. In a similar context, considering that canola sprouts are enriched in arginine and proline (Chung et al., 1989), the overproduction of AST enzymes in canola sprout-supplemented medium at 3 h might be a bacterial strategy to exploit plant arginine as an alternative nitrogen source. Bacteria utilize putrescine not only as a sole nitrogen source, but also as a signaling molecule in response to unfavorable conditions such as oxidative stress, which is a plausible stress in the context of contact with plant tissues (Schneider and Reitzer, 2012). Interestingly, arginine can be converted to putrescine, suggesting that the two pathways are closely related (Schneider et al., 2013).
Bacteria colonizing plants are able to scavenge nitrate, a source of nitrogen in plants ( Mantelin and Touraine, 2004), and utilize it as an alternative electron acceptor for energy generation (Bonnefoy and Demoss, 1994). E. coli produces two different membrane-bound respiratory nitrate reductases, nitrate reductase A (NRA) and nitrate reductase Z (NRZ), for dissimilative nitrate reduction (Bonnefoy and Demoss, 1994). Transcriptome analysis revealed that the expression of genes encoding NRA (narG and narK) and NRZ (narUZYWV) were both significantly altered during adaptation to a plant environment ( Table 2); there was an increase in narGK expression at both time points, while narUZYWV expression decreased at 1 h but increased at 3 h. These results suggest that FORC_035 might operate an aerobic respiration system for shorter exposure times to plant material, but convert to anaerobic respiration for longer exposure times, using nitrate leaked from plant tissues as an electron acceptor.

(iii) Carbohydrate Transport
Given the role of sugars and carbohydrates as carbon sources, E. coli employs a range of different transporter systems for efficient sugar acquisition, including the phosphotransferase system (PTS) and ATP-binding cassette (ABC) transporter system (Braun, 2003). The genes sgcCBX, mglB, ytfQRT, and ytfF were up-regulated during FORC_035 adaptation to both plant types ( Table 2). Genes of the sgc operon encode a PTS and are involved in uptake and breakdown of pentose and pentitol (Reizer et al., 1996;Tchieu et al., 2001). Products of mglB, ytfQRT, and ytfF comprise ABC transporters for galactopyranose (mglB) or galactofuranose (ytfQRT and ytfF) (Horler et al., 2009). When preferred carbon sources such as glucose are in limited supply, bacteria undergo catabolite repression to consume secondary carbon sources and improve their fitness in response to changes in nutrients (Stülke and Hillen, 1999). E. coli deprived of glucose is known to up-regulate the expression of genes required for uptake of alternative sugars such as galactose, ribose, and mannose (Busby and Kolb, 1996;Raman et al., 2005). Upregulation of genes involved in diverse carbohydrate transport at 3 h post-inoculation in this study suggests that FORC_035 might experience nutrient starvation at 3 h after contact with plants and switch metabolic pathways to utilize alternative carbon sources abundant in plant tissues.

(iv) Iron Acquisition
Iron is a critical nutrient for bacterial growth and also important for virulence in pathogenic bacteria (Miethke and Marahiel, 2007). Therefore, bacteria possess many different iron transport systems (Porcheron et al., 2013). E. coli also exploits a variety of chelating compounds for iron uptake and siderophores are the primary chelators responsible for ferric iron acquisition (Miethke and Marahiel, 2007). E. coli produces redundant siderophores including enterobactin, salmochelin, and aerobactin. In addition, yersiniabactin, a bacterial siderophore found in various pathogenic Yersinia species, is also present in several pathogenic E. coli strains (Garcia et al., 2011;Johnstone and Nolan, 2015). The complex of ferric ion and siderophore is actively transported across the outer membrane using ATP hydrolysis mediated by the TonB-ExbB-ExbD system (Braun, 2003). Once in the periplasm, periplasmic binding proteins (PBPs) capture ferric-siderophores and deliver them to cytoplasmic membrane ABC transporters to release ferric iron into the cytosol for subsequent utilization as a nutrient (Braun, 2003). Transcriptome analysis revealed a significant change in the expression of genes associated with iron uptake. In particular, the expression of genes required for the production of yersiniabactin and enterobactin siderophores declined during bacterial adaptation to plant tissues ( Table 2). Taking into account the high binding affinities of siderophores for ferric iron under iron-deficient conditions, decreased expression of siderophore genes at 3 h post-inoculation might be indicative of an abundant supply of iron from plant exudates.

Enterobactin Is Required for Survival of FORC_035 in Water Dropwort
Iron is an essential nutrient involved in a variety of biological processes in bacteria, including DNA synthesis, electron transport system, ATP synthesis, and oxygen transport (Andrews et al., 2003;Krewulak and Vogel, 2008). Multiple genes relevant to iron uptake showed an opposing expression tendency between 1 and 3 h post-inoculation; up-regulation at 1 h versus downregulation at 3 h ( Table 2 and Supplementary Figure S4). The expression levels were quantified by qRT-PCR, supporting consistency with the RNA-Seq results (Supplementary Figure S5). A similar gene expression pattern was observed when bacterial cells were cultured at room temperature (25 • C) in the presence of fresh produce (Supplementary Figures S6, S7). These results suggest that FORC_035 might encounter a dynamic change in iron abundance during its adaptation to the plant environment and control its iron uptake ability to improve its fitness. Especially, genes associated with two siderophores, enterobactin and yersiniabactin, were down-regulated by 3 log-fold or more when bacterial cells were exposed to fresh produce for 3 h ( Table 2). Considering the indispensable roles of iron in numerous biological processes, we evaluated the importance of iron uptake via enterobactin and yersiniabactin during bacterial adaptation to the plants by comparing the growth of entB and ybtS strains in the presence and absence of plant tissues. Deletion of entB did not compromise bacterial growth significantly in M9 minimal medium supplemented with canola sprouts. Taking into account the observation that co-culturing with plant tissues stimulated FORC_035 to alter transcription of diverse genes required for iron transport, including those involved in ferrous ion and ferric citrate transport (Supplementary Figure S4), we reasoned that the entB strain might be able to circumvent the lack of enterobactin by activating iron transport systems other than enterobactin to acquire iron released from canola sprouts. However, lack of entB dampened bacterial growth significantly in the presence of water dropwort, which has a lower iron content than sprouts (Park and Kim, 1996) (Figure 3A and Supplementary Figure S8). This result suggests that the entB strain suffered from a shortage of iron when interacting with water dropwort due to its lack of enterobactin. To test this possibility, ferric sulfate was added to the culture medium at 10 µM as an iron supplement and bacterial growth was compared. As expected, addition of ferric ion complemented the growth defect of the entB strain in water dropwort-containing medium ( Figure 3A). Moreover, the growth defect was complemented with trans-encoded entB by the pEntB plasmid (Figure 3B), indicating that iron uptake via enterobactin is required for the growth of bacterial cells exposed to water dropwort. On the other hand, the mutation on yersiniabactin showed no defects in bacterial growth. These results suggested that enterobactin plays a greater role than yersiniabactin does for iron acquisition. Differential contribution to bacterial growth by enterobactin and yersiniabactin has been reported previously in Klebsiella pneumonia (Lawlor et al., 2007).

Gallium as an Iron Mimetic Inhibits FORC_035 Growth and Biofilm Formation
The transition metal gallium (Ga 3+ ), which has an ionic radius nearly identical to iron (Fe 3+ ), can act as a "Trojan horse" to disturb bacterial iron uptake, as many biologic systems, including bacteria, cannot distinguish between Ga 3+ and Fe 3+ (Chitambar and Narasimhan, 1991;Kaneko et al., 2007). Taking into account that the entB strain suffering from iron shortage due to its lack of enterobactin slowed down its growth when in contact with water dropwort, disturbance of cellular iron uptake by the addition of gallium could be a promising strategy to control the growth and survival of EHEC FORC_035 on water dropwort. To examine this possibility, FORC_035 in M9 minimal medium containing water dropwort was supplemented with Ga(NO 3 ) 3 at 1 mM and its growth was evaluated (Figure 4). FORC_035 cocultured with water dropwort multiplied faster than FORC_035 grown in the absence of plant tissues, but the addition of gallium led to a 1.7 log reduction in bacterial numbers, indicating that FIGURE 3 | Comparison of bacterial growth between FORC_035 wild-type (WT) and entB mutant strains. (A) Ferric sulfate (10 µM) was added to M9 minimal media with or without water dropwort. (B) Four strains of WT, entB, and entB with either empty plasmid or pEntB were cultivated in the presence of water dropwort. pUHE21-2lacI q plasmid was used and 50 µM of IPTG was used for induction. In both graphs, live bacterial numbers were measured in M9 minimal medium broth supplemented with water dropwort after a 5-h incubation. Viability of bacterial cells was measured by comparing log(N/N 0 ) values, where N 0 is the initial number of cells at 0 h and N is the number of cells after a 5-h incubation. The t-test was used to evaluate the significance of differences in viability, and significance is indicated as follows: * P < 0.05; * * P < 0.01; ns, not significant.
gallium had antimicrobial activity against EHEC FORC_035. Gallium had a growth inhibitory effect regardless of the presence of water dropwort, but gallium exerted greater growth inhibition when bacteria were in contact with water dropwort tissues than when they were not (Figure 4). The greater growth inhibition of gallium in the presence of water dropwort was probably because the growth-promoting effect of the plant (Figure 1) was abolished by gallium, which functioned as a competitive antagonist of iron in the plant-supplemented environment. Antimicrobial effect of gallium has been observed in other bacteria as FIGURE 4 | Effect of gallium nitrate on the growth of FORC_035. Wild-type FORC_035 was cultivated in M9 minimal medium broth with or without water dropwort and Ga(NO 3 ) 3 was added to the medium at 1 mM. Bacterial numbers were estimated at 9 h post-inoculation. Note that log(N/N 0 ) is plotted on y-axis, where N 0 is the initial number of cells at 0 h and N is the number of cells after a 9-h incubation. Asterisks indicate significant differences ( * * P < 0.01).
well, including Pseudomonas aeruginosa and Mycobacterium tuberculosis wherein gallium inhibits cellular Fe-dependent metabolic pathways by substituting Fe 3+ (Olakanmi et al., 2000;Kaneko et al., 2007). In mammalian cells, Ga 3+ makes up Gatransferrin complex instead of Fe 3+ , suggesting the possibility that gallium exerts its role by combining with iron-chelating proteins. Biofilm formation is achieved through a complex process influenced by multiple environmental signals. Among diverse signals, iron can also regulate biofilm formation and this regulation depends on bacterial species. The availability of iron influences the ability of bacteria to form biofilms effectively; high concentrations of iron promote bacterial biofilm development in E. coli and P. aeruginosa (Banin et al., 2005;Wu and Outten, 2009). Therefore, the influence of gallium on bacterial biofilm formation was investigated. The presence of Ga(NO 3 ) 3 at 10 µM inhibited biofilm formation by threefold without changing the growth of FORC_035 (Figure 5). The inhibitory effect of gallium on biofilm formation was not influenced by temperature changes (Supplementary Figure S9). Biofilms enable bacteria to resist stressful conditions such as disinfectants and the defense responses of plants and animals, which makes biofilms an important issue to deal with in food hygiene (Srey et al., 2013). Gallium was approved for hypercalcemia treatment by the FDA, and has low toxicity at appropriate doses (Chitambar and Narasimhan, 1991;Gómez et al., 1992). The inhibitory effects of gallium on bacterial growth and biofilm formation suggest that gallium can potentially be used as an antimicrobial agent.

CONCLUSION
In conclusion, the transcriptomes of FORC_035 in contact with canola sprouts and water dropwort provided insights into the overall bacterial transcriptional response to plants. There were significant changes in the expression of genes associated with fitness and survival upon co-culture with plant materials. In particular, genes required for iron uptake were found to be critical determinants of the adaptation of EHEC FORC_035 to plants.

DATA AVAILABILITY
The database of the Foodborne pathogen Omics Research Center (http://forcdb.snu.ac.kr) in which the datasets for this manuscript were deposited are not publicly available due to policy of the institution. Requests to access the datasets should be directed to SR.

AUTHOR CONTRIBUTIONS
HN, HY, and SR conceived and designed the experiments. HN, YK, and DK performed the experiments. HN, YK, DK, and HY analyzed the data and wrote the paper. HY and SR revised the manuscript.

FUNDING
This study was supported by a grant (14162MFDS972) from the Ministry of Food and Drug Safety, South Korea, in 2018. This work was also supported by the BK21 Plus Program of the Department of Agricultural Biotechnology, Seoul National University, Seoul, South Korea.