Skip to main content

Transcriptomic analysis reveals specific metabolic pathways of enterohemorrhagic Escherichia coli O157:H7 in bovine digestive contents

Abstract

Background

The cattle gastrointestinal tract (GIT) is the main enterohemorrhagic Escherichia coli (EHEC) reservoir. In order to identify nutrients required for the survival or multiplication of EHEC in the bovine GIT, we compared the transcriptomes of the EHEC O157:H7 reference strain EDL933 cultured in vitro in bovine digestive contents (DCs) (rumen, small intestine and rectum) using RNA-sequencing.

Results

Gene expression profiles showed that EHEC EDL933 activated common but also specific metabolic pathways to survive in the different bovine DCs. Mucus-derived carbohydrates seem important in EHEC nutrition in posterior DCs (small intestine and rectum) but not in rumen content. Additional carbohydrates (xylose, ribose, mannitol, galactitol) as well as gluconeogenic substrates (aspartate, serine, glycerol) would also be used by EHEC as carbon and/or nitrogen sources all along the bovine GIT including the rumen. However, xylose, GalNac, ribose and fucose transport and/or assimilation encoding genes were over-expressed during incubation in rectum content compared with rumen and intestine contents, and genes coding for maltose transport were only induced in rectum. This suggests a role for these carbohydrates in the colonization of the cattle rectum, considered as the major site for EHEC multiplication. In contrast, the transcription of the genes associated with the assimilation of ethanolamine, an important nitrogen source for EHEC, was poorly induced in EHEC growing in rectum content, suggesting that ethanolamine is mainly assimilated in the cattle rumen and small intestine. Respiratory flexibility would also be required for EHEC survival because of the redundancy of dehydrogenases and reductases simultaneously induced in the bovine DCs, probably in response to the availability of electron donors and acceptors.

Conclusion

EHEC EDL933 showed a high flexibility in the activation of genes involved in respiratory pathways and assimilation of carbon and nitrogen sources, most of them from animal origin. This may allow the bacterium to adapt and survive in the various bovine GIT compartments. Obtaining a better understanding of EHEC physiology in bovine GIT is a key step to ultimately propose strategies to limit EHEC carriage and shedding by cattle.

Background

Enterohemorrhagic Escherichia coli (EHEC) are Shiga toxin-producing E. coli (STEC) which represent a well-known group of foodborne zoonotic pathogens distributed worldwide. Their ability to produce Shiga toxins in human gut constitutes their main virulence attribute and leads to diseases ranging from uncomplicated diarrhea to hemorrhagic colitis (HC) and hemolytic and uremic syndrome (HUS), the principal cause of kidney failure in children [1]. Numerous outbreaks and sporadic cases of severe HC and HUS have been attributed to the O157:H7 serotype [2]. Epidemiological studies indicate that the healthy bovine gastrointestinal tract (GIT) is the main reservoir of EHEC O157:H7 [3]. Bovines are asymptomatic carriers because they lack globotriaosylceramide, the Shiga toxin receptor necessary for intestinal and renal endothelium pathophysiological effects [4, 5]. Bovine feces excreted in the environment are the principal source of herd contamination and human infection. Human infection is typically acquired through the consumption of contaminated food (undercooked meat, dairy products, fruits and vegetables) and water. On farm strategies to decrease EHEC carriage by cattle should efficiently reduce shedding and consequently human EHEC infection incidence [6]. Understanding of EHEC physiology during its transit through the bovine GIT and identifying nutrients preferentially used by EHEC to survive in the bovine gut are critical to propose nutritional strategies limiting EHEC shedding.

In the rumen, EHEC O157:H7 are able to survive but not grow, probably due to inhibition by the endogenous microbiota and/or physicochemical conditions [7,8,9]. Indeed, the rumen is not considered as an EHEC colonization site. However, the pathogen must be able to survive in the rumen, to cross the acidic barrier of the abomasum and to transit to the posterior digestive compartments, constituting more favorable environments for EHEC growth. E. coli O157 have been detected throughout the gut of experimentally inoculated or naturally infected animals even if the rectum is considered the main EHEC colonization site [10,11,12,13,14,15].

Little is known about the nutrients used and the metabolic pathways required for EHEC survival and/or multiplication in the bovine GIT. The GIT content is a complex ecosystem composed of a dense and diverse microbiota, and EHEC survival relies on competition for a variety of energy, carbon and nitrogen sources. According to the Freter niche theory, each bacterial species must use one or several nutrients more efficiently than other species to colonize the mammalian intestine [16]. In this complex ecosystem, available nutrients derive from the animal diet as well as eukaryotic cell secretion and debris released in the luminal content after eukaryotic and prokaryotic cell apoptosis. Such nutrients can be of different composition and are present at various concentrations depending on the intestinal site. We previously showed that carbohydrates released by the mucus layer covering the small intestine enterocytes (mannose, N-acetylglucosamine [GlcNAc], N-acetylneuraminic acid [Neu5Ac] and galactose) were used by EHEC to grow in bovine small intestine content [17]. Furthermore, each monosaccharide was used more rapidly by EHEC than by the endogenous microbiota, suggesting an effective assimilation of mucus-derived carbohydrates by the pathogen [17]. In addition, previous studies showed that i) the genes dctA and agaB, encoding the specific transport of C4-dicarboxylates and N-acetylgalactosamine (GalNac), respectively, are involved in EHEC colonization in the calf intestine and ii) fucose is an important carbon source for maintenance of EHEC in the cattle rectum [18, 19]. Also, ethanolamine (EA), released in the mammal gut by animal, plant and bacterial cell membrane phospholipids, is an important nitrogen source giving a competitive growth advantage to EHEC in bovine small intestine content [20]. Furthermore, the gluconeogenesis pathway is also involved in the survival and maintenance of EHEC in bovine small intestine content and several compounds such as glycerol, lactate and amino acids were identified as gluconeogenic substrates present in the bovine small intestine [21]. However, the presence of nutrients and their assimilation by EHEC in other digestive compartments are poorly documented.

The objective of this study was to identify nutrients and metabolic pathways potentially used by EHEC in different bovine digestive compartments, from rumen to rectum, in order to improve our knowledge of EHEC physiology in the animal GIT. We compared the transcriptomic profiles of an EHEC O157:H7 reference strain grown in vitro in rumen, small intestine and rectum contents and in minimal medium by the RNA-sequencing (RNA-seq) approach. Quantification of differentially expressed transcripts represents a method of choice to reveal bacterial transcriptomes in complex environments such as the mammalian intestine. Our results highlighted that EDL933 is well-adapted to the bovine gut environment and activates common but also specific metabolic pathways to survive in the contents of different bovine digestive compartments.

Results

Growth and survival of EHEC EDL933 in bovine digestive contents

The culture conditions of E. coli strains in bovine digestive contents (DCs) (rumen, small intestine, caecum, colon and rectum contents) were chosen to mimic the physiological environment of each compartment (see the Materials and Methods section). Bacterial growth of EDL933 RifR was analyzed in each DC and compared to the growth patterns of the spontaneous RifR mutants of the strains Sakai (O157:H7 EHEC), NV95 (bovine O157:H7 STEC) and BG1 (bovine commensal E. coli) (Fig. 1) (Additional file 1: Table S1). As shown in Fig. 1a, EDL933 RifR was able to grow in all DCs, except in rumen content. In the rumen, a ≈ 0.5 log CFU mL-1 decrease in the EDL933 RifR population was observed after 8h of incubation and the strain survived until 24h at a low level (< 104 CFU mL-1). In contrast, EDL933 RifR grew with a high yield (≈ 4 log CFU mL-1 increase) in small intestine content and remained at high concentration up to 24h of incubation. In caecum, colon and rectum contents, the EDL933 RifR population increased from ≈ 4 to ≈ 7 log CFU mL-1 and then remained stable up to 24h of incubation. Similar growth patterns were obtained with the other E. coli strains incubated in each DC (no statistically significant difference was observed between the strains tested) (Additional file 1: Figure S2).

Fig. 1
figure1

Growth curves of EDL933 RifR incubated in bovine digestive contents. (a) Unfiltered DCs; (b) Filtered DCs. CFU enumeration at each time point is the mean of at least three independent experiments. The presented values are the log10 mean number of CFU mL-1 ± the standard error of the mean (SEM). Statistical analyses were performed using the two-way ANOVA with the Bonferroni post-hoc test. In unfiltered DCs, growth yield in small intestine content was significantly different from that in the other DCs (p < 0.001); growth yield in caecum, colon and rectum contents was significantly different from that in rumen and small intestine contents (p < 0.001). In filtered DCs, growth yield in all DCs was significantly different from that in rumen content (p < 0.001); growth yield in small intestine content was significantly different from that in caecum and colon contents (p < 0.001); growth yield in rectum content was significantly different from that in caecum and colon contents (p < 0.05)

Short chain fatty acid (SCFA) concentration and pH were measured before and after incubation in the DC. In rumen content, the pH became more acidic than in the other DCs and decreased about 1.4 pH units after 24h of incubation, indicating a fermentative activity of the endogenous ruminal microbiota (Table 1). In the other DCs, pH values were all between 7 and 8 whatever the incubation time (Table 1) which corresponded to more favorable conditions for E. coli multiplication. Before incubation, the total SCFA concentration was higher in rumen content than in the other DC’s (Table 1) with a high acetate concentration (84.51 ± 6.93 mM in rumen content, result not shown). Otherwise, total SCFA concentration increased in all DCs after 24h of incubation, indicating active endogenous microbiota (Table 1).

Table 1 pH and SCFA concentrations in unfiltered bovine digestive contents before incubation (t=0) and after 6h and 24h of EHEC EDL933 incubation

EHEC EDL933 growth was also monitored in the absence of endogenous microbiota (filtered DC, Fig. 1b). The growth rate and yield of the bacteria grown in filtered or unfiltered DC (Fig. 1a and b) were similar except a slightly but not statistically significant lower rate in filtered colon content. This suggests that the lack of growth of EDL933 in rumen content is due to the physicochemical conditions (more acidic pH and higher SCFA concentration) and/or the presence of inhibitory compounds in rumen fluids rather than nutritional competition with the endogenous microbiota.

Transcriptome profiling of EHEC EDL933 incubated in filtered bovine digestive contents

To provide insight about nutrients and metabolic pathways preferentially used by EHEC to survive in the bovine gut, we performed RNA-seq analysis to compare global gene expression changes in EHEC EDL933 grown in filtered DCs (rumen, small intestine and rectum contents) and in M9 minimal medium supplemented with glucose as the sole carbon source (M9-Glc). Note that these DCs were chosen because (i) the rumen is the first digestive compartment encountered by EHEC during their transit, (ii) the small intestine is the central compartment in intestinal absorption and (iii) the rectum is described as the primary site of EHEC colonization and multiplication and is the last compartment before fecal excretion. M9-Glc was chosen as a reference medium because i) it was used as such in previous transcriptomic studies in our laboratory [17, 20, 21] and ii) M9-Glc is a completely defined medium, with glucose as sole carbon source; iii) glucose concentration is very low in the bovine GIT [17].

Transcriptomic profiles were determined from mRNA samples collected during mid- and late-exponential growth phases (3h and 6h of incubation, respectively), except for rumen content in which RNA samples were collected only after 6h of incubation (Additional file 1: Figure S1). The whole transcriptomes of EDL933 grown in the different DCs were first analyzed by a principal component analysis. The first principal component (PC1, accounting for 41% of the variation) and the second principal component (PC2, accounting for 22% of the variation) together explained 63% of the variance in the RNA-Seq data (Fig. 2). The transcriptomes of EDL933 grown in M9-Glc and harvested at the different growth phases formed two clusters which were distant from the other transcriptomes. In addition, the global gene expression profiles of EDL933 grown in rumen and rectum contents also showed different patterns.

Fig. 2
figure2

Principal component analysis results for the transcriptomes of EDL933 grown in digestive contents. Transcriptomes of EDL933 incubated in rumen content during 6h (purple circles), and in small intestine content (red), rectum content (blue) and M9-Glc (green) at the mid-exponential (triangle) and the late-exponential (circle) growth phases. Three biological replicates were analyzed for each condition

Statistical analysis revealed that from 540 to 701 genes carried by the EDL933 genome were differentially expressed (log2 fold-change [log2FC] > │2│; q-value < 0.05) in the three DCs compared to M9-Glc (Table 2). We also compared the number of genes up and down-regulated between the digestive contents (Table 2). The highest differences were found between rumen and rectum contents after 6h incubation (678 genes differentially expressed). To validate the data, the expression level of 32 genes found to be differentially expressed using RNA-Seq measures was quantified by RT-qPCR on samples collected under all experimental conditions. The results confirmed alteration of the targeted gene transcription in the different DC relative to M9-Glc (Additional file 1: Table S4). The sets of genes up and down-regulated in the three different DCs relative to M9-Glc are presented in Additional file 2: Table S5.

Table 2 Number of differentially expressed genes (log2 FC > │2│; q < 0.05) in EHEC EDL933 incubated in filtered digestive contents and M9-Glc

To facilitate subsequent analyses, the differentially expressed genes were classified into cluster of orthologous groups (COGs) in order to assign a physiological function to each gene. Note that some genes were classified in more than one category. The results showed that 36.4% of the differentially expressed genes in EDL933 grown in the three DCs were found within the functional categories “Carbohydrate transport and metabolism” (n=333), “Amino acid transport and metabolism” (n=467) and “Energy production and conversion” (n=321) (Additional file 3: Table S6). Genes of unknown and unclassified functions represented 22.9% of the genes differentially expressed in the three DCs (Additional file 3: Table S6). Most of the other functional categories contained less than 10% of differentially expressed genes (Additional file 3: Table S6).

In depth analysis of the genes commonly or specifically up-regulated in EHEC EDL933 incubated in the bovine DCs and included in the main functional categories described above was conducted in order to predict i) the metabolic pathways used and compounds assimilated by EHEC in bovine DCs and ii) the sequential nutrition of EHEC during incubation in the bovine DCs.

Genes up-regulated relative to carbon sources

Carbon sources common to the rumen, small intestine and rectum contents

We first analyzed the genes commonly up-regulated in EHEC EDL933 incubated in the three DCs compared to M9-Glc. A Venn diagram (Fig. 3) shows that 95 genes were commonly up-regulated in EDL933 incubated in the three DCs after 6h of incubation. As shown in Fig. 4 and Table S7 (Additional file 3), 52 of the 95 genes (54.7%) were associated with the three main functional categories mentioned above. Expression of the genes coding for proteins involved in the transport and/or metabolism of carbon sources (carbohydrates, amino acids, glycerol, and C4-dicarboxylates) (Additional file 3: Table S7) were subsequently analyzed.

Fig. 3
figure3

Venn diagram of up-regulated genes in EHEC EDL933 incubated in bovine DCs. The cells were collected after 6h of incubation in the three DCs compared to M9-Glc

Fig. 4
figure4

Functional classification of the 95 genes commonly up-regulated in EDL933 after 6h of incubation in filtered bovine DCs relative to M9-Glc

The genes xylF and xylA (transport and catabolism of xylose, respectively), fucA and fucO (fucose catabolism) and rbsACD (ribose transporter) were commonly up-regulated in EDL933 grown in the three DCs (Table 3). The genes agaEVW and agaYZ encoding the transport and catabolism of N-acetylgalactosamine (GalNAc), respectively, as well as the genes cmtB (mannitol transporter) and z4875 and z4876 (galactitol transporter) also exhibited increased expression (Table 3). Note that most of these genes (except agaY, rbsACD, z4875 and z4876) were expressed at higher levels in EDL933 incubated in small intestine and rectum contents than in rumen content (Table 3).

Table 3 RNA-Seq data

The transcription of genes involved in the transport and assimilation of amino acids also appeared to be commonly induced in EHEC EDL933 in the three DCs (Fig. 4, Additional file 3: Table S7). The genes encoding i) the aerobic transport of aspartate and deamination of aspartate to fumarate (dctA and aspA, respectively), ii) the conversion of serine (sdaA) and tryptophane (tnaA) to pyruvate, iii) the degradation of threonine and serine via the propanoate pathway (tdcE and tdcD, respectively) and iv) the hydrolysis of asparagine to aspartate (asnA) were found commonly over-expressed in EHEC EDL933 (Table 3). Note that DctA is also the aerobic transporter of common C4-dicarboxylates (succinate, fumarate and malate). All these genes were more transcribed in EDL933 incubated in rumen and small intestine contents (log2FC > 3) than in rectum content, except dctA and sdaA (Table 3). In addition, glpF and glpT (transport of glycerol and glycerol 3-phosphate [G3P], respectively), and glpK (ATP-dependent glycerol kinase) were commonly induced in EDL933 incubated in the three DCs (Table 3).

Carbon sources common to the posterior intestinal contents (small intestine and rectum contents)

The rumen is not considered as a site of EHEC colonization while the small intestine and the rectum constitute more favorable environments for EHEC multiplication. Therefore, we focused here on the genes commonly up-regulated in the two posterior DCs. Genes included in “Carbohydrate transport and metabolism”, “Amino acid transport and metabolism” and “Energy production and conversion” COG categories were found commonly up-regulated in EHEC at the mid- (81 genes) and late-exponential (79 genes) growth phases.

Up-regulation of the genes encoding proteins involved in the transport or metabolism of N-acetylneuraminic acid (Neu5Ac, nanAT), N-acetylglucosamine (GlcNAc, nagABE) and galactose (galKMPT and mglABC) was unique to EDL933 incubated in posterior DCs (Table 3). The genes fucAIOP and z0461 encoding proteins involved in the transport and assimilation of fucose were also over-expressed in both small intestine and rectum contents during the two growth phases, while only fucA and fucO appeared over-expressed in rumen content (Table 3). In addition to these genes associated with mucus-derived carbohydrate assimilation, genes encoding proteins involved in the transport and/or utilization of arabinose (araABDFG), fructose (fruABK), gluconate (gntKPTU) and hexuronate (galacturonate and glucuronate) (uxaABC and uxuAB) also exhibited increased expression in the posterior compartments. Noteworthy, i) uxaABC and gntKPTU were exclusively expressed during the mid-exponential growth phase and ii) the fru genes were more induced during incubation in small intestine than in rectum content. In addition to carbohydrates, EHEC EDL933 also induced pathways involved in the transport (sdaC) and catabolism (sdaB) of serine (Table 3).

Specific carbon sources in rectum contents

Because the rectum is the major site of EHEC colonization in the bovine GIT, particular attention was paid to the genes specifically up-regulated in EHEC during incubation in rectum content. The transcription of the genes malEFGK and manXYZ encoding maltose and mannose uptake, respectively, was exclusively induced in EDL933 during incubation in rectum content, mainly during the mid-exponential growth phase (Table 3). The genes encoding proteins required for the transport and assimilation of gluconate were all up-regulated in rectum content at mid-exponential growth phase (Table 3). In addition, the transcription of the genes encoding the transport (xylF) and catabolism (xylAB) of xylose was induced in EHEC EDL933 during the mid-exponential growth phase in rectum but not in small intestine contents (Table 3). Note also that i) all the aga genes were up-regulated in both mid- and late-exponential growth phases in rectum content but only during the late-exponential growth phase in small intestine content, ii) the fold change increase in aga genes expression during the late-exponential growth phase was higher in rectum than in small intestine or rumen contents, and iii) the genes xylF, xylA, fucA, fucO and rbsACD (transport or assimilation of xylose, fucose or ribose) had higher transcripts levels in EDL933 incubated in rectum content than in the two other DCs during the late-exponential growth phase (Table 3). The genes coding for the transport and assimilation of galactose were expressed at a higher level in rectum than in small intestine contents (except mglC during the late-exponential growth phase). Importantly, the genes encoding the fucose-sensing two component system fusK and fusR (z0462 and z0463) were over-expressed in EHEC grown in rectum content during both growth phases while only the regulator fusR was up-regulated in EHEC incubated in small intestine content during the mid-exponential growth phase.

To broaden this analysis, the expression of the genes mentioned above was also directly compared in bacteria grown in rectum vs small intestine content and rectum vs rumen content (Additional file 3: Table S8). The results confirmed the specific up-regulation of several genes involved in the transport or assimilation of GalNac, fucose, maltose, ribose and xylose (Additional file 3: Table S8). However, this analysis did not allow us to confirm the specific up-regulation of mannose and gluconate transport and/or assimilation encoding genes in the rectum content.

Taken together, these results suggest an important role for xylose, fucose, GalNac, ribose and maltose in the colonization of the cattle rectum by EHEC.

Genes up-regulated relative to transport and assimilation of ethanolamine

Ethanolamine (EA) is an important nitrogen source that confers a growth advantage to EHEC EDL933 in bovine small intestine content [20]. Here, we focused on the genes encoding EA transport and assimilation in the different DCs. As expected, the transcription of the eut genes was induced in small intestine content during the late-exponential growth phase (Table 3). Interestingly, most of the eut genes also exhibited an increased expression in EHEC incubated in rumen content at 6h of incubation whereas only 7 of the 14 eut genes were up-regulated in EHEC incubated in rectum content (Table 3). Note that eutH, encoding an active transport of EA across the bacterial membrane, was up-regulated in both small intestine and rumen contents but not in rectum content. None of the eut genes were found differentially expressed in EHEC during the mid-exponential growth phase in small intestine and rectum contents.

Genes up-regulated relative to C4-dicarboxylate transporters

The transport of common C4-dicarboxylates (fumarate, malate, succinate) as well as aspartate under aerobic and anaerobic conditions is well documented [22]. In this study, the gene dcuB (anaerobic antiporting of succinate against fumarate) was up-regulated in EHEC incubated in the three DC whereas the transcription of dcuA (anaerobic antiporting of succinate against aspartate) was only induced during incubation in small intestine and rectum contents (Table 3). Surprisingly, dctA (aerobic uptake of C4-dicarboxylates) was up-regulated in EHEC during incubation in all the DCs tested, including rumen and rectum contents under anaerobiosis, suggesting differences in dctA regulation depending on the bacterial growth medium. Interestingly, succinate was excreted by EHEC EDL933 in small intestine content (1.1 mM ± 0.13 and 5.3 mM ± 0.13 succinate at 3 and 6h of incubation, respectively), but not in rumen or rectum contents, suggesting that antiporting of aspartate against succinate only occurred in the bovine small intestine.

Genes relative to energy production and conversion

An in-depth analysis of the genes associated with respiratory pathways activated by EHEC EDL933 in the three DCs was undertaken. As shown in Table 3, common and distinct respiratory systems were induced by the bacteria according to the bovine DCs. The genes fdnG and fdoHI (formate dehydrogenase-N and -I respectively), frdABCD (fumarate reductase), glpABD (G3P dehydrogenase), hyaABCD and hybABC (hydrogenase-1 and -2 respectively) were commonly up-regulated in EHEC EDL933 in the three DCs. Otherwise, appBC (cytochrome bd oxydase), dmsABC (DMSO reductase), napAH and narGHIJ (nitrate reductases) as well as nrfEFG (formate-dependent nitrite reductase) were up-regulated in both rumen and small intestine contents (Table 3). Furthermore, the gene cydB (cytochrome bd oxidase) was up-regulated in both small intestine and rectum contents. The transcription of the genes sdhABC (succinate dehydrogenase) and nuoHIJKLM (NADH oxidoreductase) was only activated in rectum content except nuoI and nuoJ which were also up-regulated in small intestine content. Note i) the redundancy of formate dehydrogenases (fdn and fdo), hydrogenases (hya and hyb) and nitrate reductases (nap and nar) encoding genes and ii) that, surprisingly, genes encoding proteins involved in aerobic respiration were found up-regulated in EHEC EDL933 incubated in rumen (fdoHI, appBC and glpD) and rectum contents (cydB, fdoHI, sdhABC, nuoHIJKLM and glpD) under anaerobic culture conditions.

Discussion

Mucin, the main constituent of the mammalian intestinal mucus, is a filamentous glycoprotein rich in oligosaccharides (up to 80% of the mucin biomass) and represents a niche for bacterial colonization of the intestine [23, 24]. Indeed, mucus is constantly released into the mammal GIT and degraded by a subset of commensal bacteria able to breakdown mucin into sugars. Monomers of the six main carbohydrates present in the mucus covering the bovine small intestine epithelium (galactose, GlcNAc, GalNAc, fucose, mannose and Neu5Ac) [25] have been previously detected in bovine small intestinal contents [17]. Importantly, galactose, GlcNAc, mannose and Neu5Ac conferred a competitive growth advantage to EHEC in bovine small intestine content [17]. However, the carbohydrate composition of the mucus layer covering other bovine digestive compartments has not yet been documented. Transcriptome analyses suggested that EHEC EDL933 could assimilate galactose, GlcNAc and Neu5Ac in posterior bovine DCs (small intestine and rectum contents) whereas fucose and GalNac could be used by the bacteria all along the bovine GIT. In addition, our results suggest that fucose and GalNac would be more important substrates for EHEC in the rectum. Altogether, these data suggest differences in mucus-derived carbohydrate composition and availability depending on bovine digestive compartments. Accordingly, renewal rate, carbohydrate composition and thickness of the mucus layer vary for different mammal gut regions. For example, nitrogen and organic matter content of bovine mucus is known to be different between the ileum and colon [26].

Survival of E. coli in the stationary growth phase depends, at least partially, on nutrients that are still available after the bacterial exponential growth phase. In this study, comparison of gene expression between mid- and late-exponential growth phases allows us to propose a sequential utilization of mucus-derived carbohydrates by EHEC in the small intestine content: galactose and GlcNAc would be rapidly assimilated by EHEC and exhausted when the bacteria reached the late-exponential growth phase; fucose and Neu5Ac would be used by EHEC up to the late-exponential growth phase whereas GalNAc would be assimilated in a second step when galactose and GlcNAc are exhausted. Accordingly, it has been shown that E. coli MG1655 also displays a preference for specific mucus derived monosaccharides at different stages of mouse intestine colonization: genes involved in GlcNAc catabolism were essential for colonization initiation whereas genes involved in fucose catabolism were required for maintenance [27]. Taken together, our results strongly suggested that mucus-derived carbohydrates i) are important carbon sources for successful survival of EHEC in the bovine gut and ii) would be assimilated by EHEC in order of preference for sugars.

Additional carbohydrates can be assimilated by EHEC in the mammal GIT. For example, arabinose, hexuronate and ribose are known to be used by EHEC EDL933 to colonize streptomycin-treated mice [28]. The concentration of free mono- and di-saccharides released from the degradation of the animal diet is generally low in the rumen of cattle because the endogenous microbiota is very efficient in assimilating these substrates [29, 30]. However, in addition to mucus-derived carbohydrates, these sugars may originate from other host derived glycans such as epithelial glycocalyx present underneath the mucus layer or from prokaryotic and eukaryotic lysed cell components. They may also be released from the degradation of dietary polysaccharides that have escaped ruminal hydrolysis. In this study, we hypothesized that EHEC EDL933 assimilated xylose, ribose, mannitol and galactitol all along the bovine GIT whereas arabinose, fructose, gluconate and hexuronates would be exclusively used by the bacteria in small intestine and rectum contents. Ribose is one of the preferred sugars used by E. coli and the catabolism of xylose and ribose allows the synthesis of precursors constituting important substrates for the pentose phosphate pathway [27, 28]. Interestingly, mutation in the ribose assimilation pathway causes a defect in the ability to colonize for EHEC, but not for commensal E. coli in the mouse intestine [28]. Gluconate is also metabolized by E. coli in the mice intestine and its distribution in drinking water impacts mice colonization by EHEC EDL933 [31, 32]. Noteworthy, gluconate transport and utilization were also induced by the probiotic bacterium Propionibacterium freudenreichii in the colonic environment in vivo compared with artificial medium in vitro [33]. Gluconic acid is present in plant products, but also in mammal intestinal mucus [27, 33]. Therefore, it is reasonable to speculate that gluconate could be present in the bovine small intestine and rectum and potentially assimilated by EHEC at these sites.

Gene transcription analysis during incubation in the bovine rectum provided valuable information since this intestinal section is referred to as the major site of O157:H7 EHEC colonization in adult cattle [19]. As mentioned above, EHEC activated the transcription of genes required for the assimilation of galactose, GlcNAc, Neu5Ac, fucose and GalNac in rectum contents. Importantly, we showed that the genes encoding PTS for uptake or assimilation of maltose were exclusively over-expressed in EHEC EDL933 during incubation in rectum content, suggesting a specific role for this sugar at this site. Maltose also plays an important role in the colonization of mice by EHEC and the utilization of this sugar confers a competitive advantage to the pathogen during the initiation of colonization when competing with commensal E. coli [31]. Maltose can have a dietary origin and be generated through the degradation of starch present in the bovine ration by the endogenous microbiota. However, in the bovine rectum, maltose would be provided by lysed cells releasing their endogenous glycogen. Indeed, glycogen storage is widespread in gut microbiota because it represents a useful strategy for rapid growth in a very competitive ecosystem [31, 34]. Also, our study strongly suggested an important role for xylose, GalNac and fucose in EHEC colonization of the rectum. Fucose was previously described as a critical carbon source for the maintenance of EHEC in the rectum of adult cattle [19]. Accordingly, our results showed that the transcription of most fuc genes, including fucP (fucose/H+ symporter) and z0461 (optimal fucose import), were activated in EHEC EDL933 incubated in rectum content. Note that i) a z0461 deletion reduces the expression of fuc genes and decreases the capacity of EHEC to use fucose to grow [35], and ii) fusK and fusR, encoding a fucose sensing two-component system, were up-regulated in EHEC during incubation in rectum content. The FusKR system is known to repress the expression of z0461 and the fuc genes during incubation of EHEC in minimal media supplemented with fucose as a sole carbon source [35]. The up-regulation of fuc, z0461 and fus genes in rectum content, a much more complex medium, strongly suggested differences in FusKR regulation depending on compounds or chemical conditions in bacterial environments. Further studies are needed to analyze the regulation of fuc gene expression under our experimental conditions.

The gluconeogenesis pathway contributes to bacterial fitness in vivo and utilization of gluconeogenic substrates is required for a maximal colonization of the mouse intestine by EHEC and for maintaining EHEC in bovine small intestine content [21, 36]. Previous studies from our laboratory showed that i) EHEC EDL933 uses preferentially mucus-derived sugars during the exponential growth phase in bovine small intestine content and then assimilates gluconeogenic substrates when sugars were exhausted, and ii) gluconeogenic substrates such as aspartate, serine and glycerol are present in small intestine content [17, 21]. The present work showed that amino acids (aspartate and serine), glycerol and G3P would also probably be used all along the bovine GIT as genes encoding their transport and/or metabolism were activated in all the DCs tested. It is well documented that bovine mucin contains a high concentration of glutamic and aspartic acids, threonine and serine [25]. Furthermore, serine and glycerol are part of phospholipids constituting eukaryotic and prokaryotic cell membranes and thus could be released into bovine gut contents during cell renewal.

Ethanolamine (EA) is included in phosphatidylethanolamine, the most abundant phospholipid in eukaryotic and prokaryotic cell membranes. Free EA, used as a nitrogen source, confers a nutritional advantage to EHEC and contributes to the survival of the pathogen in bovine small intestine content [20]. In a previous study, we showed that the concentration of free EA is much more important in bovine small intestine (2.2 mM) than in rumen (7.7 nM) or colon (17.3 nM) contents [20]. In the present report, most of the eut genes were found up-regulated in EHEC incubated in small intestine and rumen contents but only half of them in rectum content. However, the weak concentration of EA in rumen content strongly suggested that EA probably cannot promote EHEC growth in the rumen of cattle. Kendall et al. [37] showed that EA could be a signal for gene regulation, and could then play complex roles in metabolism, cell-to-cell signaling and bacterial virulence. However, the specific mechanisms of gene regulation by EA remain to be completely elucidated. The results of our study suggest that EA could positively regulate the expression of genes whose products are involved in EHEC survival in the rumen, but further studies will be necessary to elucidate its role in the cattle rumen.

The mammalian intestine is not strictly anaerobic but contains distinct microaerobic and anaerobic regions due to i) dynamic cycles of oxygen diffusion from intestinal epithelium and ii) oxygen consumption by facultative anaerobes [38,39,40,41]. It is well documented that nitrate is the preferred terminal electron acceptor during nitrate respiration in E. coli when oxygen became limited in the mammal intestine [22, 40, 42]. In contrast, fumarate constitutes the terminal electron acceptor during fumarate respiration, when oxygen is absent and nitrate is present at low levels [22, 40, 42]. In this study, EHEC EDL933 up-regulated simultaneously several aerobic and anaerobic respiratory systems in all the DCs tested, suggesting that EHEC colonization of the bovine GIT requires respiratory flexibility. Accordingly, microaerobic (oxygen) and anaerobic (nitrate and fumarate) niches are equally crucial for successful EHEC colonization of the mouse intestine [40, 41]. Taken together, our transcriptome analyses strongly suggested that EHEC could potentially use distinct dehydrogenases and reductases during its transit in the bovine GIT. This could be due to availability and fluctuation of electron donors (H2, formate, succinate, G3P) and acceptors (oxygen, nitrate, fumarate) in the animal gut. The redundancy and complexity in the respiratory systems activated by EHEC could facilitate a rapid adaptation of the bacteria in changing intestinal environments.

Fumarate, detected in the mammalian intestine at low levels, is endogenously generated by bacteria via central metabolism. Interestingly, the transcription of the genes encoding anaerobic antiporting of succinate against aspartate or fumarate described by Unden et al. [22] (aspA, dcuA, dcuB) was activated in EHEC incubated in the three DCs. We have previously shown that EHEC assimilate aspartate present in bovine small intestinal content [21] and we demonstrated in this study that EHEC EDL933 excreted succinate in small intestine content (but not in rumen or rectum contents) whereas fumarate was undetectable before and after EHEC incubation (results not shown). Taken together, these results strongly suggested i) the anaerobic antiporting of external aspartate against succinate due to the DcuA antiporter and ii) a supply of intracellular fumarate due to aspartate deamination under anaerobiosis. We speculated that the supply of intracellular fumarate via aspartate deamination might represent an additional source of electron acceptors used by EHEC during anaerobic fumarate respiration in the bovine small intestine. This hypothesis was recently confirmed during incubation of EHEC in small intestine content [43].

Conclusions

In conclusion, EHEC EDL933 showed a high flexibility in activation of genes involved in respiratory pathways and assimilation of carbon and nitrogen sources, most of them from animal origin. Such a flexibility may allow the bacterium to adapt and survive in the various segments of the bovine GIT with different nutrient composition. Our results open new avenues to find strategies designed to specifically reduce intestinal nutrients required for E. coli O157:H7 colonization. To this end, administration of probiotics with a high efficiency in the assimilation of the nutrients preferred by EHEC to ruminants may be an effective pre-harvest strategy to limit EHEC carriage and shedding by cattle.

Methods

Bacterial strains characteristics

The wild type E. coli strains and the rifampicin-resistant (RifR) mutants used in this study were described in Additional file 1: Table S1. Bacteria were routinely cultured on Luria Bertani (LB) agar plates. Spontaneous RifR mutants of E. coli NV95 and Sakai strains were isolated by culturing the wild-type strain on LB agar plates containing 100 μg mL-1 of rifampicin. The wild-type strains and their corresponding spontaneous RifR mutants showed identical growth curves when incubated in LB broth at 37°C with aeration (results not shown).

Origin and sampling of digestive contents

Eight healthy young bulls (Bos taurus) from the “Herbipole experimental INRA farm”, approximately 2 years of age and 550 kg mean weight, were slaughtered in the experimental slaughterhouse of the “Herbipole” (National Institute for Agronomic Research [INRA], Saint-Genès-Champanelle, France) (Slaughterhouse Permit number: 63345001). The bulls were raised according to current INRA ethical guidelines for animal welfare and the experiments were approved by the local ethics committee (Permit Number: C6334517). Digestive contents (DCs) were collected from each bull fed a mixed diet containing hay (80%) and concentrate (20%) (major compounds: rapeseed [26%], maize [20%], barley [14%], sunflower [8.5%], triticale colza [8.5%]). All DCs from rumen, small intestine, caecum and colon compartments were collected from each animal at slaughter while rectum contents were collected two days before slaughter by rectal palpation. Small intestine contents were collected from the jejunum and ileum removed in a single piece as previously described [20]. All DCs were rapidly collected and immediately brought to the laboratory. Small intestine contents were directly distributed in sterile tubes without any particular attention paid regarding anaerobiosis. The other contents were processed under strictly anaerobic conditions as previously described [44]. Briefly, rumen contents were filtered through four layers of cheesecloth to remove large feed particles. Caecum, colon and rectum contents were diluted 1:1 in reduced potassium phosphate buffer (50 mM potassium phosphate, resazurin 0.1%, 40 mM Na2CO3, 3 mM cysteine, pH 7.6) which was prepared in order to maintain a low redox potential [44]. Rumen, caecum, colon and rectum samples were distributed in sterile O2-free CO2-saturated Hungate tubes (Bellco). All these DC samples, containing the live endogenous microbiota and noted as unfiltered samples, were frozen at -80°C until use. These conditions were previously described as an appropriate way to maintain the viability of the bovine intestinal microbiota [7]. DC samples from the eight animals were pooled before use. Endogenous microbiota rifampicin sensibility was confirmed by spotting each pooled DC (100 μL) on Sorbitol Mac Conkey (SMAC) agar plates supplemented with 100 μg mL-1 rifampicin before incubation at 37°C for 24h.

To remove endogenous microbiota, the frozen unfiltered DC samples were thawed, centrifuged twice for 15 min at 10,000 × g, and supernatants were filtered through a Steritop and Stericup system with a membrane pore size 0.22 μm (Millipore). Filtrates, except small intestine filtrates, were placed in Hungate tubes which were left without stoppers in an anaerobic chamber (JACOMEX, Lyon, France) under 80-90 ppm of oxygen during three days at room temperature. This allowed the filtrates to be under anaerobic conditions. Tubes were then recapped before being removed from the anaerobic chamber, refiltered through 0.22 μm pore-size filters (Millipore) and placed into new O2-free CO2-saturated sterile Hungate tubes. The filtration efficacy was checked by inoculating an aliquot of the filtrates into LB agar plates before overnight incubation at 37°C. Filtered DCs were stored at room temperature until use.

Inoculation of E. coli strains in unfiltered and filtered digestive contents

The wild type strains EDL933, Sakai, NV95 and BG1, and their respective spontaneous RifR mutants were inoculated from a single colony and incubated in LB medium without antibiotic or supplemented with rifampicin (100 μg mL-1), respectively, for 7h at 37°C with aeration. The cultures were then 50-fold diluted in filtered DCs and grown overnight at 39°C without aeration. The day after, the bacterial concentration was spectrophotometrically adjusted at 600 nm to ≈ 104 bacteria mL-1 before inoculating unfiltered and filtered DCs. DCs were finally incubated at 39°C (internal bovine temperature) and i) under strict anaerobiosis with gentle shaking (rumen content) or without shaking (caecum, colon and rectum contents) or ii) under oxygen-limited conditions without shaking in small intestine content as previously described [21]. These conditions (temperature, shaking and oxygenation) were chosen to reflect the in vivo conditions for each bovine digestive compartment [7, 20, 21, 44]. At each time point (every two hours), an aliquot of unfiltered and filtered DC was 10-fold serially diluted in sterile phosphate buffered saline buffer (pH 7.4) and spotted on SMAC agar plates supplemented with rifampicin (100 μg mL-1) or without antibiotic, respectively. The plates were incubated overnight at 37°C before counting colony forming unit (CFU). Each experiment was replicated at least three times.

pH and metabolite quantification

Metabolite concentration was quantified in unfiltered and filtered DCs. The samples were first centrifuged (10,000 × g for 10 min). For short chain fatty acid (SCFA) quantification in unfiltered DCs, supernatants were filtered (0.22 μm) and 30 μL of orthophosphoric acid (75%) were added to 1 mL of supernatant. Total SCFA concentrations were determined by gas chromatography. For succinate quantification in filtered DCs, 200 μL of ZnSO4 5% and Ba (OH)2 3 M were added to 800 μL of supernatant before quantification determined by high-performance liquid chromatography. All the analysis were performed by AFYREN INVESTMENT (Biopole Clermont Limagne, Saint Beauzire, France). pH measurements were performed for unfiltered DCs using a HI-8424N pH meter (HANNA instruments).

RNA extraction and mRNA enrichment

The M9 minimal medium [45] was supplemented with glucose (40 mM), MgSO4 (1 mM), CaCl2 (0.1 mM) and trace metals (M9-Glc), and adjusted to ≈ pH 7.4. EHEC EDL933 inoculated from a single colony was incubated in filtered DCs under growth conditions described above, except that the initial bacterial concentrations were adjusted to ≈ 108 cells mL-1 before inoculating rumen content, ≈ 106 cells mL-1 before inoculating small intestine and rectum contents and ≈ 107 cells mL-1 before inoculating M9-Glc (Additional file 1: Figure S1). The M9-Glc bacterial cultures were grown under oxygen-limited conditions. Three biological replicates were performed for each culture condition. After incubation, bacterial suspensions were centrifuged at 10,000 × g for 15 min. The supernatants were stored at -20°C for further investigation and the bacterial pellets were flash frozen in liquid nitrogen and rapidly stored at -80°C. The next day, total RNA was extracted as previously described [44]. Briefly, the bacterial pellet was resuspended in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8) in 2 mL Graduated Skirted tube with tethered screw cap (BioSpec, USA) containing 600 mg 0.1 mm diameter zirconia/silica beads (Biospec, USA), 1 volume of AquaPhenol™ (pH 4.5) (MP Biomedicals), 1/10 volume of 10% sodium dodecyl sulfate, and 3.5 μL β-mercaptoethanol (PROLABO, France). The lysis of bacterial cells was performed using a FastPrep®-24 instruments (MP Biomedicals) twice for 30 seconds with a speed of 6 m/s. Total RNA was then purified from the bacterial pellet using the Nucleospin® RNA (Macherey Nagel) according to the manufacturer recommendations. To assess DNA contamination, PCR control on each sample was performed using the tufA primers (TGGTTGATGACGAAGAGCTG and GCTCTGGTTCCGGAATGTAG). The RNA concentration was measured using a Nanodrop ND-1000 spectrophotometer (Nanodrop technologies, France) and the RNA quality was analyzed using an Agilent 2100 Bioanalyser (Agilent technologies, France). The high quality of each RNA samples was confirmed with 23S/16S rRNA ratio ≈ 2 and RNA Integrity Number ≥ 8.

Enriched fractions of mRNA were prepared using the MicrobExpress™ Bacterial mRNA Purification kit (Ambion) according to the manufacturer instructions. The reduction in 16S and 23S rRNA in mRNA enriched fractions was confirmed using an Agilent 2100 Bioanalyser (Agilent technologies, France).

RNA Sequencing

RNA-Seq was performed at the GeT-PlaGe core facility, INRA Toulouse, France. RNA-seq libraries were prepared according to Illumina’s protocols using the Illumina TruSeq Stranded mRNA sample prep kit. Briefly, after mRNA enrichment, 200 ng of mRNA were fragmented to generate double stranded cDNA and adapters were ligated. A total of 10 cycles of PCR were applied to amplify libraries. Library quality was assessed using an Advanced Analytical Fragment Analyzer and libraries were quantified by qPCR using the Kapa Library Quantification kit. RNA-seq experiments were performed on an Illumina HiSeq3000 using a paired-end read length of 2x150 bp with the Illumina HiSeq3000 chemistry. A total of 11 to 42 million paired-reads per sample was obtained, except for the sample EDL933_21 (122 million paired-reads) (Additional file 1: Table S2).

RNA-Seq data analysis and bioinformatics

Raw read files have been stored in ng6 [46] and were checked using fastQC [47]. The sequencing adapters were removed using cutadapt (version 1.8.3, standard parameters) [48]. The reads were then aligned to EDL933 genome (Genbank accession numbers NZ_CP008957.1 and NZ_CP008958.1) [49] using bwa mem (version 0.7.12-r1039, standard options) [50]. Reads were counted using featureCount (version v1.4.5-p1) [51]. Read counts corresponding to the 131 rRNA, tRNA and ncRNA genes were excluded [49]. Differential gene expression was then performed using DESeq2 version 1.12.4 [52] with R version 3.3.2 [53] following the standard workflow. The DESeq2 method internally corrects for library size and uses negative binomial generalized linear models to test for differential expression. All genes with a log2 fold-change (log2FC) in expression greater than 2 and a Benjamini-Hochberg adjusted p-value (q-value) smaller than 0.05 were reported as differentially expressed. For further analysis, additional gene annotations from older published EDL933 chromosome and plasmid sequences were collected (Genbank accession numbers CP008957.1 and CP008958.1 [49]; AE005174.2 [54] and AF074613.1 [55]). Additional gene annotations were also performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG). Differentially expressed genes were assigned to functional categories of Clusters of Orthologous Groups (COGs) of proteins using blastp against the NCBI COG 2014 database [56]. RNA-seq data have been deposited under SRA accession SRP136076.

Note that the fold change was first calculated by comparing the expression ratio of each gene from a specific DC relative to M9-Glc. This allowed us to identify the genes commonly up-regulated in several DCs. The fold change was also calculated for each gene by directly comparing its expression level between two digestive compartments, in order to identify genes specifically up-regulated in each DC.

Reverse transcription followed by quantitative PCR (RT-qPCR)

One microgram of each RNA sample (in triplicates) was reverse transcribed using the SuperScript II Reverse Transcriptase kit (Invitrogen) with 3 μg of random primer and 100 units of SuperScript II Rnase H. Quantitative PCR runs were carried out using the Mastercycler ep realplex apparatus (Eppendorf) with 20 ng of cDNA, 0.5 mM of each primer, 10 μL of SYBR Premix Ex Taq mix (Takara Bio Inc.) in a final volume of 20 μL. Amplification conditions were as follows: 95°C for 15 s, 55°C for 15 s, and 72°C for 20 s. The house keeping gene mdh was used for normalization of mRNA quantification (Additional file 1: Table S3). The relative RNA quantification was performed using primers designed to specifically amplify fragments of 90 to 200 bp (Additional file 1: Table S3). Triplicate samples were amplified in each case. Results were calculated using the comparative cycle threshold method. The results presented are average from at least duplicate experiments.

Statistical Analysis

Statistical analyses of the growth and survival of E. coli strains in DCs were analysed by the two-way ANOVA with the Bonferroni post-hoc test with p < 0.05 considered as significant. For RT-qPCR, Student t-test and the one-way ANOVA with the Tukey post-hoc test were used to identify significant differences in gene expression (p < 0.05 was considered as significantly different). The statistical analysis of RNA-seq data is presented above. PCA plot was generated using plotPCA function from DESeq2 package version 1.12.4 and ggplot2 version 2.2.0.

Abbreviations

CFU:

Colony forming unit

COG:

Cluster of orthologous group

DC:

Digestive content

EA:

Ethanolamine

EHEC:

Enterohemorrhagic Escherichia coli

GalNAc:

N-acetylgalactosamine

GIT:

Gastrointestinal tract

GlcNAc:

N-acetylglucosamine

HC:

Hemorrhagic colitis

HUS:

Hemolytic and uremic syndrome

Log2FC:

Log2 fold-change

M9-Glc:

M9 minimal medium supplemented with glucose

Neu5Ac:

N-acetylneuraminic acid

RifR :

Rifampicin-resistant

RNA-seq:

RNA sequencing

SCFA:

Short chain fatty acid

STEC:

Shiga toxin producing Escherichia coli

References

  1. 1.

    Gomes TA, Elias WP, Scaletsky IC, Guth BE, Rodrigues JF, Piazza RM, et al. Diarrheagenic Escherichia coli. Braz J Microbiol. 2016;47(Suppl 1):3–30.

    CAS  Article  Google Scholar 

  2. 2.

    Kaper JB, O'Brien AD. Overview and historical perspectives. Microbiol Spectr. 2014;2(6):EHEC-0028-2014.

  3. 3.

    Persad AK, LeJeune JT. Animal reservoirs of shiga toxin producing Escherichia coli. Microbiol Spectr. 2014;2(4):EHEC-0027-2014.

    PubMed  Google Scholar 

  4. 4.

    Hoey DE, Currie C, Else RW, Nutikka A, Lingwood CA, Gally DL, et al. Expression of receptors for verotoxin 1 from Escherichia coli O157 on bovine intestinal epithelium. J Med Microbiol. 2002;51(2):143–9.

    CAS  Article  Google Scholar 

  5. 5.

    Pruimboom-Brees IM, Morgan TW, Ackermann MR, Nystrom ED, Samuel JE, Cornick NA, et al. Cattle lack vascular receptors for Escherichia coli O157:H7 shiga toxins. Proc Natl Acad Sci U S A. 2000;97(19):10325–9.

    CAS  Article  Google Scholar 

  6. 6.

    Saeedi P, Yazdanparast M, Behzadi E, Salmanian AH, Mousavi SL, Nazarian S, et al. A review on strategies for decreasing E. coli O157:H7 risk in animals. Microb Pathog. 2017;103:186–95.

    CAS  Article  Google Scholar 

  7. 7.

    Chaucheyras-Durand F, Faqir F, Ameilbonne A, Rozand C, Martin C. Fates of acid-resistant and non-acid-resistant shiga toxin-producing Escherichia coli strains in ruminant digestive contents in the absence and presence of probiotics. Appl Environ Microbiol. 2010;76(3):640–7.

    CAS  Article  Google Scholar 

  8. 8.

    Chaucheyras-Durand F, Madic J, Doudin F, Martin C. Biotic and abiotic factors influencing in vitro growth of Escherichia coli O157:H7 in ruminant digestive contents. Appl Environ Microbiol. 2006;72(6):4136–42.

    CAS  Article  Google Scholar 

  9. 9.

    de Vaux A, Morrison M, Hutkins RW. Displacement of Escherichia coli O157:H7 from rumen medium containing prebiotic sugars. Appl Environ Microbiol. 2002;68(2):519–24.

    CAS  Article  Google Scholar 

  10. 10.

    Grauke LJ, Kudva IT, Yoon JW, Hunt CW, Williams CJ, Hovde CJ. Gastrointestinal tract location of Escherichia coli O157:H7 in ruminants. Appl Environ Microbiol. 2002;68(5):2269–77.

    CAS  Article  Google Scholar 

  11. 11.

    Keen JE, Laegreid WW, CG C-MK, Durso LM, Bono JL. Distribution of shiga-toxigenic Escherichia coli O157 in the gastrointestinal tract of naturally O157-shedding cattle at necropsy. Appl Environ Microbiol. 2010;76(15):5278–81.

    CAS  Article  Google Scholar 

  12. 12.

    Kudva IT, Dean-Nystrom EA. Bovine recto-anal junction squamous epithelial (RSE) cell adhesion assay for studying Escherichia coli O157 adherence. J Appl Microbiol. 2011;111(5):1283–94.

    CAS  Article  Google Scholar 

  13. 13.

    Lim JY, Li J, Sheng H, Besser TE, Potter K, Hovde CJ. Escherichia coli O157:H7 colonization at the rectoanal junction of long-duration culture-positive cattle. Appl Environ Microbiol. 2007;73(4):1380–2.

    CAS  Article  Google Scholar 

  14. 14.

    Low JC, McKendrick IJ, McKechnie C, Fenlon D, Naylor SW, Currie C, et al. Rectal carriage of enterohemorrhagic Escherichia coli O157 in slaughtered cattle. Appl Environ Microbiol. 2005;71(1):93–7.

    CAS  Article  Google Scholar 

  15. 15.

    Naylor SW, Low JC, Besser TE, Mahajan A, Gunn GJ, Pearce MC, et al. Lymphoid follicle-dense mucosa at the terminal rectum is the principal site of colonization of enterohemorrhagic Escherichia coli O157:H7 in the bovine host. Infect Immun. 2003;71(3):1505–12.

    CAS  Article  Google Scholar 

  16. 16.

    Freter R, Brickner H, Botney M, Cleven D, Aranki A. Mechanisms that control bacterial populations in continuous-flow culture models of mouse large intestinal flora. Infect Immun. 1983;39(2):676–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Bertin Y, Chaucheyras-Durand F, Robbe-Masselot C, Durand A, de la Foye A, Harel J, et al. Carbohydrate utilization by enterohaemorrhagic Escherichia coli O157:H7 in bovine intestinal content. Environ Microbiol. 2013;15(2):610–22.

    CAS  Article  Google Scholar 

  18. 18.

    Dziva F, van Diemen PM, Stevens MP, Smith AJ, Wallis TS. Identification of Escherichia coli O157:H7 genes influencing colonization of the bovine gastrointestinal tract using signature-tagged mutagenesis. Microbiology. 2004;150(11):3631–45.

    CAS  Article  Google Scholar 

  19. 19.

    Snider TA, Fabich AJ, Conway T, Clinkenbeard KD. E. coli O157:H7 catabolism of intestinal mucin-derived carbohydrates and colonization. Vet Microbiol. 2009;136(1-2):150–4.

    CAS  Article  Google Scholar 

  20. 20.

    Bertin Y, Girardeau JP, Chaucheyras-Durand F, Lyan B, Pujos-Guillot E, Harel J, et al. Enterohaemorrhagic Escherichia coli gains a competitive advantage by using ethanolamine as a nitrogen source in the bovine intestinal content. Environ Microbiol. 2011;13(2):365–77.

    CAS  Article  Google Scholar 

  21. 21.

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

    Article  Google Scholar 

  22. 22.

    Unden G, Strecker A, Kleefeld A, Kim OB. C4-dicarboxylate utilization in aerobic and anaerobic growth. EcoSal Plus. 2016;7(1). ESP-0021-2015.

  23. 23.

    Deplancke B, Gaskins HR. Microbial modulation of innate defense: goblet cells and the intestinal mucus layer. Am J Clin Nutr. 2001;73(6):1131S–41S.

    CAS  Article  Google Scholar 

  24. 24.

    Sicard JF, Le Bihan G, Vogeleer P, Jacques M, Harel J. Interactions of intestinal bacteria with components of the intestinal mucus. Front Cell Infect Microbiol. 2017;7:387.

    Article  Google Scholar 

  25. 25.

    Montagne L, Toullec R, Lalles JP. Calf intestinal mucin: isolation, partial characterization, and measurement in ileal digesta with an enzyme-linked immunosorbent assay. J Dairy Sci. 2000;83(3):507–17.

    CAS  Article  Google Scholar 

  26. 26.

    Aperce CC, Heidenreich JM, Drouillard JS. Capacity of the bovine intestinal mucus and its components to support growth of Escherichia coli O157:H7. Animal. 2014;8(5):731–7.

    CAS  Article  Google Scholar 

  27. 27.

    Chang DE, Smalley DJ, Tucker DL, Leatham MP, Norris WE, Stevenson SJ, et al. Carbon nutrition of Escherichia coli in the mouse intestine. Proc Natl Acad Sci U S A. 2004;101(19):7427–32.

    CAS  Article  Google Scholar 

  28. 28.

    Fabich AJ, Jones SA, Chowdhury FZ, Cernosek A, Anderson A, Smalley D, et al. Comparison of carbon nutrition for pathogenic and commensal Escherichia coli strains in the mouse intestine. Infect Immun. 2008;76(3):1143–52.

    CAS  Article  Google Scholar 

  29. 29.

    Bertin Y, Habouzit C, Dunière L, Laurier M, Durand A, Duchez D, et al. Lactobacillus reuteri suppresses E. coli O157:H7 in bovine ruminal fluid: Toward a pre-slaughter strategy to improve food safety? PLoS One. 2017;12(11):e0187229.

    Article  Google Scholar 

  30. 30.

    McSweeney CS, Mackie R. Micro-organisms and ruminant digestion: state of knowledge, trends and future prospects. In: Food and Agriculture Organization of the United Nation Background study paper N°61; 2012.

    Google Scholar 

  31. 31.

    Jones SA, Jorgensen M, Chowdhury FZ, Rodgers R, Hartline J, Leatham MP, et al. Glycogen and maltose utilization by Escherichia coli O157:H7 in the mouse intestine. Infect Immun. 2008;76(6):2531–40.

    CAS  Article  Google Scholar 

  32. 32.

    Leatham MP, Banerjee S, Autieri SM, Mercado-Lubo R, Conway T, Cohen PS. Precolonized human commensal Escherichia coli strains serve as a barrier to E. coli O157:H7 growth in the streptomycin-treated mouse intestine. Infect Immun. 2009;77(7):2876–86.

    CAS  Article  Google Scholar 

  33. 33.

    Saraoui T, Parayre S, Guernec G, Loux V, Montfort J, Le Cam A, et al. A unique in vivo experimental approach reveals metabolic adaptation of the probiotic Propionibacterium freudenreichii to the colon environment. BMC Genomics. 2013;14:911.

    Article  Google Scholar 

  34. 34.

    Matheron C, Delort AM, Gaudet G, Forano E, Liptaj T. 13C and 1H nuclear magnetic resonance study of glycogen futile cycling in strains of the genus Fibrobacter. Appl Environ Microbiol. 1998;64(1):74–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Pacheco AR, Curtis MM, Ritchie JM, Munera D, Waldor MK, Moreira CG, et al. Fucose sensing regulates bacterial intestinal colonization. Nature. 2012;492(7427):113–7.

    CAS  Article  Google Scholar 

  36. 36.

    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(3):1666–76.

    CAS  Article  Google Scholar 

  37. 37.

    Kendall MM, Gruber CC, Parker CT, Sperandio V. Ethanolamine controls expression of genes encoding components involved in interkingdom signaling and virulence in enterohemorrhagic Escherichia coli O157:H7. MBio. 2012;3(3).

  38. 38.

    Espey MG. Role of oxygen gradients in shaping redox relationships between the human intestine and its microbiota. Free Radic Biol Med. 2013;55:130–40.

    CAS  Article  Google Scholar 

  39. 39.

    He G, Shankar RA, Chzhan M, Samouilov A, Kuppusamy P, Zweier JL. Noninvasive measurement of anatomic structure and intraluminal oxygenation in the gastrointestinal tract of living mice with spatial and spectral EPR imaging. Proc Natl Acad Sci U S A. 1999;96(8):4586–91.

    CAS  Article  Google Scholar 

  40. 40.

    Jones SA, Chowdhury FZ, Fabich AJ, Anderson A, Schreiner DM, House AL, et al. Respiration of Escherichia coli in the mouse intestine. Infect Immun. 2007;75(10):4891–9.

    CAS  Article  Google Scholar 

  41. 41.

    Jones SA, Gibson T, Maltby RC, Chowdhury FZ, Stewart V, Cohen PS, et al. Anaerobic respiration of Escherichia coli in the mouse intestine. Infect Immun. 2011;79(10):4218–26.

    CAS  Article  Google Scholar 

  42. 42.

    Unden G, Bongaerts J. Alternative respiratory pathways of Escherichia coli: energetics and transcriptional regulation in response to electron acceptors. Biochim Biophys Acta. 1997;1320(3):217–34.

    CAS  Article  Google Scholar 

  43. 43.

    Bertin Y, Segura A, Jubelin G, Dunière L, Durand A, Forano E. Aspartate metabolism is involved in the maintenance of enterohaemorrhagic Escherichia coli O157:H7 in bovine intestinal content. Environ Microbiol. 2018. https://doi.org/10.1111/1462-2920.14380.

  44. 44.

    Segura A, Auffret P, Bibbal D, Bertoni M, Durand A, Jubelin G, et al. Factors involved in the persistence of a shiga toxin-producing Escherichia coli O157:H7 strain in bovine feces and gastro-intestinal content. Front Microbiol. 2018;9:375.

    Article  Google Scholar 

  45. 45.

    Sambrook J, Fritsch EF, Maniatis T. Current protocols in molecular biology. In: Ausubel FM, Brent M, Kingston RE, Moore JG, editors. Molecular cloning: A laboratory manual. 1 and 2. Seconde ed. New York: Greene Publishing Associates and John Wiley & Sons; 1989. p. 1626.

    Google Scholar 

  46. 46.

    Mariette J, Escudie F, Allias N, Salin G, Noirot C, Thomas S, et al. NG6: Integrated next generation sequencing storage and processing environment. BMC Genomics. 2012;13:462.

    CAS  Article  Google Scholar 

  47. 47.

    Andrews S. FastQC: a quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (2010). Accessed 10 Nov 2016.

  48. 48.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.

    Article  Google Scholar 

  49. 49.

    Latif H, Li HJ, Charusanti P, Palsson BO, Aziz RK. A gapless, unambiguous genome sequence of the enterohemorrhagic Escherichia coli O157:H7 strain EDL933. Genome Announc. 2014;2(4):e00821–14.

    Article  Google Scholar 

  50. 50.

    Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    CAS  Article  Google Scholar 

  51. 51.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

    CAS  Article  Google Scholar 

  52. 52.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  Google Scholar 

  53. 53.

    R Core Team, R: A language and environment for statistical computing The R Project for Statistical Computing, Vienna, Austria. https//www.r-project.org (2016). Accessed 20 Dec 2016.

  54. 54.

    Perna NT, Plunkett G 3rd, Burland V, Mau B, Glasner JD, Rose DJ, et al. Genome sequence of enterohaemorrhagic Escherichia coli O157:H7. Nature. 2001;409(6819):529–33.

    CAS  Article  Google Scholar 

  55. 55.

    Burland V, Shao Y, Perna NT, Plunkett G, Sofia HJ, Blattner FR. The complete DNA sequence and analysis of the large virulence plasmid of Escherichia coli O157:H7. Nucleic Acids Res. 1998;26(18):4196–204.

    CAS  Article  Google Scholar 

  56. 56.

    Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, et al. The NCBI BioSystems database. Nucleic Acids Res. 2010;38(Database issue):D492–6.

    CAS  Article  Google Scholar 

Download references

Acknowledgments

The authors thank Laura Quentin for excellent technical assistance. We thank the Herbipole (Abattoir) personnel for their help in collecting the bovine gastrointestinal contents. The sequencing was performed in collaboration with the GeT core facility, Toulouse, France (http://get.genotoul.fr). We are also grateful to the Genotoul bioinformatics platform Toulouse Midi-Pyrenées (Genotoul Bioinfo) for providing computing and storage resources. The authors thank Dr RJ Forster for careful reading and editing of the manuscript.

Funding

This work was supported by INRA. AS was the recipient of an INRA-Région Auvergne PhD fellowship (contract 23000731). The GeT sequencing platform was supported by France Genomique National infrastructure, funded as part of “Investissement d’Avenir” program managed by Agence Nationale pour la Recherche (contract ANR-10-INBS-09). The funding bodies had no role either in the design, analysis and interpretation of data, or in the writing of the manuscript.

Availability of data and materials

The datasets generated during the current study are available in the NCBI SRA repository under project accession number PRJNA439245 (https://www.ncbi.nlm.nih.gov/Traces/study/?acc=PRJNA439245). The R script used for the RNAseq analysis is publicly available at https://github.com/paulineauffret/TransEHEC.

Author information

Affiliations

Authors

Contributions

EF, AS and YB conceived and designed the study. AS, MB and AD performed the biological experiments. OB and CG performed the sequencing experiment. PA and CK performed the bioinformatics analysis. AS, YB and EF analyzed the data and wrote the manuscript. All authors read, improved and approved the final manuscript.

Corresponding author

Correspondence to Evelyne Forano.

Ethics declarations

Ethics approval consent to participate

Healthy young bulls were slaughtered in the experimental slaughterhouse of the “Herbipole” (National Institute for Agronomic Research [INRA], Saint-Genès-Champanelle, France) (Slaughterhouse Permit number: 63345001). The bulls were raised according to current INRA and national ethical guidelines for animal welfare and the experiments were approved by the local ethics committee C2E2A (Comité d'Ethique pour l'Expérimentation Animale Auvergne) (Permit Number: C6334517).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Supplemental Figures showing the growth of E. coli strains in digestive contents and minimal media and supplemental tables describing the strains used in this study, the RNA-seq mapping assessment, the PCR primers and RT-qPCR results. (DOCX 193 kb)

Additional file 2:

Supplemental Tables showing the results of RNAseq statistical analysis of EDL933 incubated in the digestive contents versus minimal medium. (XLSX 557 kb)

Additional file 3:

 Supplemental Tables showing the COG classification of differentially expressed genes in EDL933 incubated in digestive contents compared to M9-Glc and genes up-regulated between digestive contents. (DOCX 31 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Segura, A., Bertoni, M., Auffret, P. et al. Transcriptomic analysis reveals specific metabolic pathways of enterohemorrhagic Escherichia coli O157:H7 in bovine digestive contents. BMC Genomics 19, 766 (2018). https://doi.org/10.1186/s12864-018-5167-y

Download citation

Keywords

  • EHEC O157:H7
  • Bovine gastrointestinal tract
  • RNA-seq
  • Metabolism
  • Carbon sources
  • Mucus-derived carbohydrates
  • Ethanolamine
  • Respiration