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

Audrey Segura, Marine Bertoni, Pauline Auffret, Christophe C. Kropp, Olivier Bouchez, Clémence Genthon, Alexandra Durand, Yolande Bertin, Evelyne Forano

► To cite this version:

Audrey Segura, Marine Bertoni, Pauline Auffret, Christophe C. Kropp, Olivier Bouchez, et al. Transcriptomic analysis reveals specific metabolic pathways of enterohemorrhagic Escherichia coli O157:H7 in bovine digestive contents. BMC Genomics, 2018, 19 (1), 10.1186/s12864-018-5167-y. hal-02628554

HAL Id: hal-02628554

https://hal.inrae.fr/hal-02628554

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

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

Audrey Segura¹, Marine Bertoni¹,⁴, Pauline Auffret¹,⁵, Christophe Klopp², Olivier Bouchez³, Clémence Genthon³, Alexandra Durand¹, Yolande Bertin¹ and Evelyne Forano¹*

<figure><img src="image_3.png" /></figure>

**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.

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

<figure><img src="image_4.png" /></figure>

BMC

© The Author(s). 2018 **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.

<footnote>* Correspondence: evelyne.forano@inra.fr
¹Université Clermont Auvergne, INRA, UMR 0454 MEDIS, F-63000 Clermont-Ferrand, France
Full list of author information is available at the end of the article</footnote>
## 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–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–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 C₄-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 Rif® 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⁻¹ decrease in the EDL933 RifR population was observed after 8h of incubation and the strain survived until 24h at a low level (≈ 10⁴ CFU mL⁻¹). In contrast, EDL933 RifR grew with a high yield (≈ 4 log CFU mL⁻¹ 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⁻¹ 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).

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).

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].

<figure><img src="image_5.png" /><figcaption>**Fig. 1** 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 log₁₀ mean number of CFU mL⁻¹ ± 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)</figcaption></figure>
<table><caption>**Table 1** pH and SCFA concentrations in unfiltered bovine digestive contents before incubation (t=0) and after 6h and 24h of EHEC EDL933 incubation</caption><thead><tr><th></th><th>Time</th><th>Rumen</th><th>Small intestine</th><th>Caecum<sup>a</sup></th><th>Colon<sup>a</sup></th><th>Rectum<sup>a</sup></th></tr></thead><tbody><tr><td rowspan="3">pH</td><td>0</td><td>7.25 (± 0.05)</td><td>7.22 (± 0.01)</td><td>7.29 (± 0.23)</td><td>7.42 (± 0.39)</td><td>7.38 (± 0.42)</td></tr><tr><td>6h</td><td>6.30 (± 0.04)</td><td>7.14 (± 0.02)</td><td>7.90 (± 0.20)</td><td>7.53 (± 0.03)</td><td>7.46 (± 0.00)</td></tr><tr><td>24h</td><td>5.85 (± 0.02)</td><td>7.66 (± 0.02)</td><td>7.88 (± 0.33)</td><td>7.38 (± 0.17)</td><td>7.59 (± 0.06)</td></tr><tr><td rowspan="3">Total SCFAs (mM)</td><td>0</td><td>115.78 (± 9.56)</td><td>15.75 (± 0.58)</td><td>36.27 (± 0.69)</td><td>38.49 (± 2.28)</td><td>19.35 (± 1.82)</td></tr><tr><td>6h</td><td>90.80<sup>b</sup></td><td>31.94<sup>b</sup></td><td>27.66<sup>b</sup></td><td>37.36<sup>b</sup></td><td>25.01 (± 0.33)</td></tr><tr><td>24h</td><td>129.53<sup>b</sup></td><td>50.31<sup>b</sup></td><td>37.01<sup>b</sup></td><td>70.61<sup>b</sup></td><td>50.21 (± 2.35)</td></tr></tbody></table>

The presented values are the mean of at least three measures (± SEM).

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.

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.

<figure><img src="image_6.png" /><figcaption>**Fig. 2** 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</figcaption></figure>

<footnote><sup>a</sup>pH and SCFA concentrations were measured in the incubation media (i.e. after dilution of DC for those which were diluted, see materials and methods)</footnote>

<footnote><sup>b</sup>SCFA concentrations were measured only once for these samples</footnote>
<table><caption>**Table 2** Number of differentially expressed genes (log2 FC > 2) ; q < 0.05) in EHEC EDL933 incubated in filtered digestive contents and M9-Glc</caption><thead><tr><th>Digestive contents</th><th>up<sup>a</sup></th><th>down<sup>b</sup></th><th>Total</th></tr></thead><tbody><tr><td>After 3 hours of incubation</td><td></td><td></td><td></td></tr><tr><td>Small intestine/M9</td><td>347</td><td>354</td><td>701</td></tr><tr><td>Rectum/M9</td><td>290</td><td>250</td><td>540</td></tr><tr><td>Small intestine/rectum</td><td>109</td><td>171</td><td>280</td></tr><tr><td>After 6 hours of incubation</td><td></td><td></td><td></td></tr><tr><td>Rumen/M9</td><td>337</td><td>280</td><td>617</td></tr><tr><td>Small intestine/M9</td><td>346</td><td>254</td><td>600</td></tr><tr><td>Rectum/M9</td><td>305</td><td>320</td><td>625</td></tr><tr><td>Rectum/rumen</td><td>277</td><td>401</td><td>678</td></tr><tr><td>Small intestine/rectum</td><td>283</td><td>173</td><td>456</td></tr></tbody></table>

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 $C_4$-dicarboxylates) (Additional file 3: Table S7) were subsequently analyzed.

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).

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 $C_4$-dicarboxylates (succinate, fumarate and malate). All these genes were more transcribed in EDL933 incubated in rumen and small intestine contents (log2FC > 3) than

<figure><img src="image_7.png" /><figcaption>**Fig. 3** 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</figcaption></figure>

<footnote><sup>a</sup>number of up-regulated genes ; <sup>b</sup>number of down-regulated genes
The comparison between small intestine and rumen after 6h of incubation
does not appear in the table because of the low significance of the
statistical analysis</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 4** Functional classification of the 95 genes commonly up-regulated in EDL933 after 6h of incubation in filtered bovine DCs relative to M9-Glc</figcaption></figure>

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. Note-worthy, 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
<table><caption>**Table 3 RNA-Seq data**</caption><thead><tr><th rowspan="3">Gene tag</th><th rowspan="3">Gene</th><th rowspan="3">Gene product</th><th colspan="4">3 hours incubation</th><th colspan="4">6 hours incubation</th></tr><tr><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Rumen</th></tr><tr><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th></tr></thead><tbody><tr><td>Z2315</td><td>acpD</td><td>FMN-dependent NADH-azoreductase</td><td>2.53</td><td>7.91E-05</td><td>5.12</td><td>7.34E-17</td><td>3.16</td><td>6.19E-07</td><td>4.77</td><td>7.25E-15</td><td>2.38</td><td>2.65E-04</td></tr><tr><td>Z4492</td><td>agaB</td><td>PTS system galactosamine-specific IIB component</td><td>NDE</td><td></td><td>3.01</td><td>2.87E-05</td><td>NDE</td><td></td><td>3.27</td><td>2.44E-06</td><td>NDE</td><td></td></tr><tr><td>Z4493</td><td>agaC</td><td>PTS system galactosamine-specific IIC component</td><td>NDE</td><td></td><td>2.11</td><td>4.41E-03</td><td>NDE</td><td></td><td>2.15</td><td>1.61E-03</td><td>NDE</td><td></td></tr><tr><td>Z4487</td><td>agaE</td><td>PTS system N-acetylgalactosamine-specific IID component</td><td>NDE</td><td></td><td>3.97</td><td>1.62E-10</td><td>2.60</td><td>4.74E-05</td><td>3.99</td><td>6.14E-11</td><td>2.44</td><td>1.35E-04</td></tr><tr><td>Z4488</td><td>agaF</td><td>PTS system N-acetylgalactosamine-specific IIA component</td><td>NDE</td><td></td><td>3.23</td><td>1.41E-08</td><td>NDE</td><td></td><td>3.20</td><td>7.02E-09</td><td>NDE</td><td></td></tr><tr><td>Z4485</td><td>agaV</td><td>PTS system N-acetylgalactosamine-specific IIB component</td><td>NDE</td><td></td><td>6.10</td><td>6.41E-23</td><td>3.22</td><td>5.89E-07</td><td>4.79</td><td>1.45E-14</td><td>2.52</td><td>1.40E-04</td></tr><tr><td>Z4486</td><td>agaW</td><td>PTS system N-acetylgalactosamine-specific IIC component</td><td>NDE</td><td></td><td>5.36</td><td>7.49E-22</td><td>3.37</td><td>5.24E-09</td><td>5.06</td><td>1.66E-19</td><td>3.07</td><td>1.25E-07</td></tr><tr><td>Z4491</td><td>agaY</td><td>Tagatose 1,6-diphosphate aldolase (GatY/KbaY subunit)</td><td>NDE</td><td></td><td>4.57</td><td>1.82E-13</td><td>2.54</td><td>7.96E-05</td><td>3.85</td><td>4.58E-10</td><td>2.69</td><td>2.25E-05</td></tr><tr><td>Z4484</td><td>agaZ</td><td>Tagatose 1,6-diphosphate aldolase (GatZ/KbaZ subunit)</td><td>NDE</td><td></td><td>5.81</td><td>2.92E-18</td><td>2.54</td><td>4.57E-04</td><td>4.24</td><td>5.68E-10</td><td>2.30</td><td>1.52E-03</td></tr><tr><td>Z1396</td><td>appB</td><td>Cytochrome bd oxidase subunit II</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.90</td><td>6.57E-08</td><td>NDE</td><td></td><td>3.18</td><td>1.60E-05</td></tr><tr><td>Z1395</td><td>appC</td><td>Cytochrome bd ubiquinol oxidase subunit I</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.33</td><td>7.02E-09</td><td>NDE</td><td></td><td>3.60</td><td>2.31E-06</td></tr><tr><td>Z0070</td><td>araA</td><td>L-arabinose isomerase</td><td>2.79</td><td>2.22E-03</td><td>NDE</td><td></td><td>3.14</td><td>6.37E-04</td><td>2.22</td><td>1.88E-02</td><td>NDE</td><td></td></tr><tr><td>Z0072</td><td>araB</td><td>L-ribulokinase</td><td>3.80</td><td>7.49E-05</td><td>3.49</td><td>4.33E-04</td><td>3.65</td><td>1.74E-04</td><td>3.82</td><td>7.55E-05</td><td>NDE</td><td></td></tr><tr><td>Z0069</td><td>araD</td><td>L-ribulose 5-phosphate 4-epimerase</td><td>2.81</td><td>1.03E-03</td><td>NDE</td><td></td><td>2.26</td><td>1.19E-02</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z2954</td><td>araF</td><td>L-arabinose transport system (substrate-binding protein)</td><td>4.06</td><td>1.61E-05</td><td>2.69</td><td>8.44E-03</td><td>4.56</td><td>1.30E-06</td><td>4.13</td><td>1.41E-05</td><td>NDE</td><td></td></tr><tr><td>Z2953</td><td>araG</td><td>L-arabinose transport system (ATP-binding protein)</td><td>2.57</td><td>6.40E-03</td><td>NDE</td><td></td><td>3.08</td><td>1.08E-03</td><td>2.09</td><td>3.18E-02</td><td>NDE</td><td></td></tr><tr><td>Z5245</td><td>asnA</td><td>Aspartate ammonia-ligase</td><td>NDE</td><td></td><td>2.68</td><td>9.66E-04</td><td>3.11</td><td>8.23E-05</td><td>2.84</td><td>3.50E-04</td><td>3.05</td><td>1.14E-04</td></tr><tr><td>Z0821</td><td>asnB</td><td>Asparagine synthase</td><td>NDE</td><td></td><td>2.84</td><td>2.54E-04</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5744</td><td>aspA</td><td>Aspartate ammonia-lyase</td><td>6.08</td><td>3.44E-23</td><td>5.03</td><td>8.65E-16</td><td>5.34</td><td>5.67E-18</td><td>2.84</td><td>1.33E-05</td><td>4.46</td><td>1.04E-12</td></tr><tr><td>Z5734</td><td>cadA</td><td>Lysine decarboxylase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.81</td><td>1.61E-03</td><td>NDE</td><td></td><td>3.33</td><td>1.21E-04</td></tr><tr><td>Z5735</td><td>cadB</td><td>Cadaverine:lysine antiporter</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.96</td><td>2.74E-05</td><td>NDE</td><td></td><td>4.15</td><td>9.05E-06</td></tr><tr><td>Z3712</td><td>cchB</td><td>Ethanolamine utilization protein EutN</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.3</td><td>4.58E-04</td><td>NDE</td><td></td><td>4.47</td><td>2.33E-04</td></tr><tr><td>Z0762</td><td>citC</td><td>[citrate (pro-3S)-lyase] ligase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.02</td><td>8.56E-03</td></tr><tr><td>Z4278</td><td>cmtB</td><td>Mannitol-specific IIA component</td><td>2.35</td><td>4.79E-03</td><td>3.04</td><td>2.28E-04</td><td>2.10</td><td>1.38E-02</td><td>3.02</td><td>1.74E-04</td><td>2.04</td><td>1.38E-02</td></tr><tr><td>Z0901</td><td>cydB</td><td>Cytochrome bd ubiquinol oxydase (subunit II)</td><td>2.12</td><td>1.29E-04</td><td>2.37</td><td>2.32E-05</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4942</td><td>dctA</td><td>Aerobic C<sub>4</sub>-dicarboxylates transporter</td><td>4.13</td><td>1.21E-09</td><td>3.91</td><td>1.73E-08</td><td>3.35</td><td>1.60E-06</td><td>3.65</td><td>1.29E-07</td><td>2.11</td><td>3.99E-03</td></tr><tr><td>Z5743</td><td>dcuA</td><td>Anaerobic C<sub>4</sub>-dicarboxylates transporter</td><td>3.21</td><td>2.46E-08</td><td>2.51</td><td>2.73E-05</td><td>2.12</td><td>5.02E-04</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5725</td><td>dcuB</td><td>Anaerobic C<sub>4</sub>-dicarboxylates transporter</td><td>2.98</td><td>2.64E-05</td><td>3.49</td><td>8.54E-07</td><td>3.26</td><td>4.39E-06</td><td>NDE</td><td></td><td>2.11</td><td>4.48E-03</td></tr><tr><td>Z1240</td><td>dmsA</td><td>Dimethyl sulfoxide reductase subunit A</td><td>3.01</td><td>3.54E-07</td><td>NDE</td><td></td><td>3.47</td><td>3.43E-09</td><td>NDE</td><td></td><td>3.07</td><td>2.25E-07</td></tr><tr><td>Z1241</td><td>dmsB</td><td>Dimethyl sulfoxide reductase subunit B</td><td>3.48</td><td>6.54E-08</td><td>NDE</td><td></td><td>3.56</td><td>3.63E-08</td><td>NDE</td><td></td><td>2.97</td><td>6.04E-06</td></tr><tr><td>Z1242</td><td>dmsC</td><td>Dimethyl sulfoxide reductase subunit C</td><td>2.56</td><td>3.33E-05</td><td>NDE</td><td></td><td>2.59</td><td>3.36E-05</td><td>NDE</td><td></td><td>2.29</td><td>2.76E-04</td></tr><tr><td>Z3707</td><td>eutA</td><td>Ethanolamine utilization protein EutA</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.52</td><td>2.76E-03</td><td>2.75</td><td>2.10E-02</td><td>4.04</td><td>4.15E-04</td></tr><tr><td>Z3706</td><td>eutB</td><td>Ethanolamine ammonia-lyase (large subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.68</td><td>1.68E-03</td><td>3.12</td><td>8.17E-03</td><td>4.02</td><td>4.60E-04</td></tr><tr><td>Z3705</td><td>eutC</td><td>Ethanolamine ammonia-lyase (small subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.37</td><td>4.18E-03</td><td>2.79</td><td>1.85E-02</td><td>3.62</td><td>1.71E-03</td></tr><tr><td>Z3711</td><td>eutE</td><td>Aldehyde dehydrogenase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.14</td><td>7.46E-04</td><td>NDE</td><td></td><td>4.50</td><td>2.01E-04</td></tr></tbody></table>
<table><caption>**Table 3 RNA-Seq data (Continued)**</caption><thead><tr><th rowspan="3">Gene tag</th><th rowspan="3">Gene</th><th rowspan="3">Gene product</th><th colspan="4">3 hours incubation</th><th colspan="4">6 hours incubation</th><th colspan="2">Rumen</th></tr><tr><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th rowspan="2">Log2FC</th><th rowspan="2">q-value</th></tr><tr><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th></tr></thead><tbody><tr><td>Z3709</td><td>eutG</td><td>Alcohol dehydrogenase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.87</td><td>1.59E-03</td><td>NDE</td><td></td><td>4.42</td><td>2.15E-04</td></tr><tr><td>Z3708</td><td>eutH</td><td>Ethanolamine transporter</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.42</td><td>3.79E-03</td><td>NDE</td><td></td><td>4.10</td><td>3.45E-04</td></tr><tr><td>Z3714</td><td>eutI</td><td>Phosphotransacetylase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.19</td><td>4.95E-04</td><td>NDE</td><td></td><td>4.51</td><td>1.43E-04</td></tr><tr><td>Z3710</td><td>eutJ</td><td>Ethanolamine utilization protein EutJ</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.64</td><td>3.76E-03</td><td>NDE</td><td></td><td>4.25</td><td>4.99E-04</td></tr><tr><td>Z3703</td><td>eutK</td><td>Ethanolamine utilization protein EutK</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.22</td><td>6.41E-03</td><td>3.01</td><td>1.02E-02</td><td>3.52</td><td>2.29E-03</td></tr><tr><td>Z3704</td><td>eutL</td><td>Ethanolamine utilization protein EutL</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.28</td><td>5.36E-03</td><td>2.82</td><td>1.70E-02</td><td>3.51</td><td>2.40E-03</td></tr><tr><td>Z3713</td><td>eutM</td><td>Ethanolamine utilization protein EutM</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.37</td><td>2.51E-04</td><td>NDE</td><td></td><td>4.64</td><td>8.61E-05</td></tr><tr><td>Z3717</td><td>eutP</td><td>Ethanolamine utilization protein EutP</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.80</td><td>1.76E-03</td><td>2.69</td><td>3.04E-02</td><td>4.33</td><td>2.60E-04</td></tr><tr><td>Z3716</td><td>eutQ</td><td>Ethanolamine utilization protein EutQ</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.87</td><td>1.56E-03</td><td>NDE</td><td></td><td>4.25</td><td>3.92E-04</td></tr><tr><td>Z3718</td><td>eutS</td><td>Ethanolamine utilization protein EutS</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.63</td><td>2.93E-03</td><td>3.27</td><td>7.36E-03</td><td>4.54</td><td>1.17E-04</td></tr><tr><td>Z2236</td><td>fdnG</td><td>Formate dehydrogenase-N (α subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.73</td><td>6.18E-06</td><td>3.48</td><td>3.63E-09</td><td>2.06</td><td>9.55E-04</td></tr><tr><td>Z2235</td><td>fdnH</td><td>Formate dehydrogenase-N (β subunit)</td><td>3.63</td><td>5.60E-08</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z2234</td><td>fdnI</td><td>Formate dehydrogenase-N (γ subunit)</td><td>3.02</td><td>2.25E-09</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5435</td><td>fdoH</td><td>Formate dehydrogenase-I (iron-sulfur subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.2</td><td>2.40E-06</td><td>3.92</td><td>3.75E-09</td><td>2.51</td><td>3.15E-04</td></tr><tr><td>Z5434</td><td>fdol</td><td>Formate dehydrogenase-I (γ subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.84</td><td>2.29E-05</td><td>3.53</td><td>6.94E-08</td><td>2.11</td><td>2.25E-03</td></tr><tr><td>Z5762</td><td>frdA</td><td>Fumarate reductase (flavoprotein subunit)</td><td>4.43</td><td>1.26E-12</td><td>2.24</td><td>9.63E-04</td><td>4.70</td><td>4.78E-14</td><td>2.02</td><td>3.04E-03</td><td>4.03</td><td>1.64E-10</td></tr><tr><td>Z5760</td><td>frdB</td><td>Fumarate reductase (iron-sulfur subunit)</td><td>4.90</td><td>4.39E-12</td><td>2.28</td><td>3.74E-03</td><td>4.88</td><td>6.63E-12</td><td>NDE</td><td></td><td>4.17</td><td>7.12E-09</td></tr><tr><td>Z5759</td><td>frdC</td><td>Fumarate reductase (C subunit)</td><td>4.92</td><td>3.51E-11</td><td>2.27</td><td>6.14E-03</td><td>4.57</td><td>1.13E-09</td><td>NDE</td><td></td><td>3.83</td><td>5.31E-07</td></tr><tr><td>Z5758</td><td>frdD</td><td>Fumarate reductase (D subunit)</td><td>4.85</td><td>1.93E-10</td><td>2.51</td><td>2.61E-03</td><td>4.04</td><td>1.93E-07</td><td>NDE</td><td></td><td>3.20</td><td>5.72E-05</td></tr><tr><td>Z3425</td><td>fruA</td><td>PTS system fructose-specific IIB component</td><td>2.95</td><td>3.71E-04</td><td>2.44</td><td>5.39E-03</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z3427</td><td>fruB</td><td>PTS system fructose-specific IIA component</td><td>3.24</td><td>4.31E-04</td><td>2.67</td><td>6.21E-03</td><td>2.59</td><td>7.71E-03</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z3426</td><td>fruK</td><td>1-phosphofructokinase</td><td>3.38</td><td>1.07E-04</td><td>2.76</td><td>2.79E-03</td><td>2.12</td><td>2.73E-02</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4117</td><td>fucA</td><td>L-fucose-phosphate aldolase</td><td>3.99</td><td>4.76E-07</td><td>4.06</td><td>4.41E-07</td><td>3.45</td><td>1.86E-05</td><td>4.24</td><td>7.00E-08</td><td>2.69</td><td>1.14E-03</td></tr><tr><td>Z4119</td><td>fucl</td><td>L-fucose/D-arabinose isomerase</td><td>4.12</td><td>5.83E-08</td><td>4.47</td><td>5.45E-09</td><td>2.87</td><td>3.07E-04</td><td>5.40</td><td>4.28E-13</td><td>NDE</td><td></td></tr><tr><td>Z4120</td><td>fucl</td><td>L-fuculokinase</td><td>2.50</td><td>3.90E-04</td><td>2.37</td><td>1.16E-03</td><td>NDE</td><td></td><td>3.34</td><td>1.32E-06</td><td>NDE</td><td></td></tr><tr><td>Z4116</td><td>fucO</td><td>Lactaldehyde reductase</td><td>3.22</td><td>1.17E-08</td><td>3.36</td><td>6.44E-06</td><td>3.39</td><td>4.25E-06</td><td>3.92</td><td>6.72E-08</td><td>2.91</td><td>9.95E-05</td></tr><tr><td>Z4118</td><td>fucP</td><td>L-fucose permease</td><td>3.62</td><td>5.13E-07</td><td>4.44</td><td>7.68E-10</td><td>2.62</td><td>5.68E-04</td><td>5.13</td><td>3.02E-13</td><td>NDE</td><td></td></tr><tr><td>Z4121</td><td>fucU</td><td>L-fucose mutarotase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.05</td><td>1.07E-03</td><td>NDE</td><td></td></tr><tr><td>Z2615</td><td>fumA</td><td>Fumarase</td><td>3.90</td><td>5.07E-19</td><td>3.35</td><td>7.45E-14</td><td>3.50</td><td>2.06E-15</td><td>3.53</td><td>1.57E-15</td><td>2.50</td><td>3.54E-08</td></tr><tr><td>Z5724</td><td>fumB</td><td>Fumarate hydratase class I</td><td>3.57</td><td>1.16E-05</td><td>3.72</td><td>6.56E-06</td><td>3.76</td><td>4.38E-06</td><td>NDE</td><td></td><td>2.68</td><td>1.56E-03</td></tr><tr><td>Z2614</td><td>fumC</td><td>Fumarate hydratase class II</td><td>2.99</td><td>4.63E-10</td><td>2.19</td><td>1.31E-05</td><td>2.46</td><td>4.60E-07</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z0462</td><td>fusK</td><td>Fucose sensing two-component system (sensor)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.45</td><td>5.35E-11</td><td>NDE</td><td></td></tr><tr><td>Z0463</td><td>fusR</td><td>Fucose sensing two-component system (regulator)</td><td>3.37</td><td>2.15E-09</td><td>3.81</td><td>1.78E-11</td><td>NDE</td><td></td><td>3.58</td><td>2.40E-10</td><td>NDE</td><td></td></tr><tr><td>Z0927</td><td>galK</td><td>Galactokinase</td><td>2.55</td><td>1.80E-05</td><td>3.50</td><td>2.71E-09</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z0926</td><td>galM</td><td>Aldolase-1-epimerase</td><td>2.04</td><td>6.59E-04</td><td>2.62</td><td>1.02E-05</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4288</td><td>galP</td><td>MFS transporter galactose:H<sup>+</sup> symporter</td><td>3.23</td><td>1.20E-04</td><td>4.66</td><td>1.82E-08</td><td>NDE</td><td></td><td>2.75</td><td>1.54E-03</td><td>NDE</td><td></td></tr><tr><td>Z0928</td><td>galT</td><td>Galactose-1-phosphate uridylyltransferase</td><td>3.24</td><td>7.98E-07</td><td>4.58</td><td>1.59E-12</td><td>NDE</td><td></td><td>2.92</td><td>1.32E-05</td><td>NDE</td><td></td></tr><tr><td>Z3499</td><td>glpA</td><td>Glycerol 3-phosphate dehydrogenase (a subunit)</td><td>4.13</td><td>3.90E-06</td><td>2.50</td><td>1.07E-02</td><td>5.21</td><td>4.03E-09</td><td>3.42</td><td>2.20E-04</td><td>4.15</td><td>4.64E-06</td></tr><tr><td>Z3500</td><td>glpB</td><td>Glycerol 3-phosphate dehydrogenase (b subunit)</td><td>4.33</td><td>4.60E-06</td><td>NDE</td><td></td><td>5.01</td><td>1.05E-07</td><td>2.20</td><td>3.19E-02</td><td>3.77</td><td>1.06E-04</td></tr></tbody></table>
<table><caption>**Table 3 RNA-Seq data (Continued)**</caption><thead><tr><th rowspan="3">Gene tag</th><th rowspan="3">Gene</th><th rowspan="3">Gene product</th><th colspan="6">3 hours incubation</th><th colspan="6">6 hours incubation</th></tr><tr><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Rumen</th></tr><tr><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th></tr></thead><tbody><tr><td>Z3501</td><td>glpC</td><td>Glycerol 3-phosphate dehydrogenase (c subunit)</td><td>3.91</td><td>2.15E-05</td><td>NDE</td><td></td><td>3.84</td><td>4.04E-05</td><td>NDE</td><td></td><td>2.51</td><td>1.03E-02</td></tr><tr><td>Z4786</td><td>glpD</td><td>Glycerol 3-phosphate dehydrogenase</td><td>3.62</td><td>8.99E-05</td><td>2.65</td><td>7.55E-03</td><td>4.98</td><td>4.13E-08</td><td>4.75</td><td>1.79E-07</td><td>4.23</td><td>4.61E-06</td></tr><tr><td>Z5472</td><td>glpF</td><td>Glycerol uptake facilitator protein</td><td>4.02</td><td>3.20E-09</td><td>3.66</td><td>1.38E-07</td><td>4.02</td><td>4.65E-09</td><td>5.33</td><td>2.65E-15</td><td>2.63</td><td>2.44E-04</td></tr><tr><td>Z5471</td><td>glpK</td><td>Glycerol kinase</td><td>3.80</td><td>3.61E-06</td><td>2.89</td><td>8.21E-04</td><td>3.87</td><td>3.11E-06</td><td>4.33</td><td>1.20E-07</td><td>2.38</td><td>6.39E-03</td></tr><tr><td>Z3498</td><td>glpT</td><td>Glycerol 3-phosphate transporter</td><td>5.65</td><td>1.72E-09</td><td>3.47</td><td>5.41E-04</td><td>5.44</td><td>8.86E-09</td><td>3.71</td><td>1.57E-04</td><td>4.05</td><td>3.23E-05</td></tr><tr><td>Z4805</td><td>gntK</td><td>Gluconokinase</td><td>NDE</td><td></td><td>2.44</td><td>3.04E-04</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5919</td><td>gntP</td><td>High-affinity gluconate transporter</td><td>5.85</td><td>5.74E-16</td><td>3.17</td><td>4.16E-05</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4770</td><td>gntT</td><td>Gnt-I system high-affinity gluconate transporter</td><td>NDE</td><td></td><td>2.52</td><td>1.95E-07</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4804</td><td>gntU</td><td>Gnt-I system low-affinity gluconate transporter</td><td>NDE</td><td></td><td>2.35</td><td>1.99E-03</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z1076</td><td>grxA</td><td>Glutaredoxin 1</td><td>2.97</td><td>1.14E-06</td><td>5.31</td><td>2.61E-19</td><td>3.99</td><td>2.46E-11</td><td>4.51</td><td>3.83E-14</td><td>4.00</td><td>2.26E-11</td></tr><tr><td>Z4879</td><td>hpR</td><td>Phosphocarrier protein</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.56</td><td>5.73E-04</td><td>2.14</td><td>7.02E-03</td><td>2.32</td><td>1.63E-03</td></tr><tr><td>Z1389</td><td>hyaA</td><td>Hydrogenase I (small subunit)</td><td>2.62</td><td>3.08E-03</td><td>2.51</td><td>6.62E-03</td><td>6.24</td><td>3.71E-14</td><td>2.38</td><td>8.91E-03</td><td>5.66</td><td>9.23E-12</td></tr><tr><td>Z1390</td><td>hyaB</td><td>Hydrogenase I (large subunit)</td><td>2.58</td><td>3.41E-03</td><td>2.48</td><td>6.98E-03</td><td>6.32</td><td>1.14E-14</td><td>2.72</td><td>2.20E-03</td><td>5.69</td><td>5.17E-12</td></tr><tr><td>Z1391</td><td>hyaC</td><td>Hydrogenase I (B-type cytochrome subunit)</td><td>2.31</td><td>9.30E-03</td><td>2.39</td><td>9.57E-03</td><td>5.77</td><td>1.92E-12</td><td>2.45</td><td>6.16E-03</td><td>5.05</td><td>1.15E-09</td></tr><tr><td>Z1392</td><td>hyaD</td><td>Hydrogenase I (maturation protease)</td><td>2.22</td><td>1.16E-03</td><td>2.42</td><td>7.60E-03</td><td>5.47</td><td>1.45E-11</td><td>2.29</td><td>9.99E-03</td><td>4.81</td><td>4.38E-09</td></tr><tr><td>Z4350</td><td>hybA</td><td>Hydrogenase-2 (small subunit)</td><td>4.01</td><td>6.41E-09</td><td>2.92</td><td>5.37E-05</td><td>2.77</td><td>1.24E-04</td><td>NDE</td><td></td><td>2.27</td><td>2.01E-03</td></tr><tr><td>Z4349</td><td>hybB</td><td>Ni/Fe component of hydrogenase-2</td><td>4.24</td><td>1.21E-09</td><td>3.13</td><td>1.59E-05</td><td>2.86</td><td>8.54E-05</td><td>NDE</td><td></td><td>2.27</td><td>2.36E-03</td></tr><tr><td>Z4348</td><td>hybC</td><td>Hydrogenase-2 (large subunit)</td><td>3.83</td><td>3.42E-08</td><td>2.8</td><td>1.09E-04</td><td>2.65</td><td>2.71E-04</td><td>NDE</td><td></td><td>2.00</td><td>7.53E-03</td></tr><tr><td>Z4347</td><td>hybD</td><td>Hydrogenase maturation protease</td><td>3.39</td><td>1.11E-06</td><td>2.33</td><td>1.68E-03</td><td>2.22</td><td>2.74E-03</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5632</td><td>malE</td><td>Maltose/maltodextrin transport system (substrate-binding protein)</td><td>NDE</td><td></td><td>3.22</td><td>4.60E-05</td><td>NDE</td><td></td><td>2.41</td><td>3.70E-03</td><td>NDE</td><td></td></tr><tr><td>Z5631</td><td>malF</td><td>Maltose/maltodextrin transport system (permease protein)</td><td>NDE</td><td></td><td>2.71</td><td>9.82E-06</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5630</td><td>malG</td><td>Maltose/maltodextrin transport system (permease protein)</td><td>NDE</td><td></td><td>2.79</td><td>3.81E-05</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5633</td><td>malK</td><td>Multiple sugar transport system (ATP-binding protein)</td><td>NDE</td><td></td><td>3.99</td><td>7.37E-08</td><td>NDE</td><td></td><td>2.36</td><td>3.38E-03</td><td>NDE</td><td></td></tr><tr><td>Z2860</td><td>manX</td><td>PTS system mannose-specific IIAB component</td><td>NDE</td><td></td><td>2.69</td><td>4.64E-08</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z2861</td><td>manY</td><td>PTS system mannose-specific IIC component</td><td>NDE</td><td></td><td>2.81</td><td>5.05E-09</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z2862</td><td>manZ</td><td>PTS system mannose-specific IID component</td><td>NDE</td><td></td><td>2.66</td><td>1.49E-07</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z3404</td><td>mglA</td><td>Methyl-galactoside transport system (ATP-binding protein)</td><td>2.98</td><td>5.30E-04</td><td>4.04</td><td>1.89E-06</td><td>2.84</td><td>1.28E-03</td><td>3.38</td><td>8.10E-05</td><td>NDE</td><td></td></tr><tr><td>Z3405</td><td>mglB</td><td>Methyl-galactoside transport system (substrate-binding protein)</td><td>3.75</td><td>1.06E-06</td><td>4.77</td><td>4.85E-10</td><td>2.76</td><td>6.21E-04</td><td>4.09</td><td>1.00E-07</td><td>NDE</td><td></td></tr><tr><td>Z3403</td><td>mglC</td><td>Methyl-galactoside transport system (permease protein)</td><td>3.02</td><td>4.18E-05</td><td>3.33</td><td>7.49E-06</td><td>3.06</td><td>4.16E-05</td><td>2.85</td><td>1.52E-04</td><td>NDE</td><td></td></tr><tr><td>Z0824</td><td>nagA</td><td>N-acetylglucosamine 6-phosphate deacetylase</td><td>2.37</td><td>7.39E-06</td><td>2.95</td><td>2.19E-08</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z0825</td><td>nagB</td><td>Glucosamine 6-phosphate deaminase</td><td>3.44</td><td>1.74E-07</td><td>4.49</td><td>7.20E-12</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z0826</td><td>nagE</td><td>PTS system N-acetylglucosamine-specific IIABC component</td><td>3.31</td><td>4.52E-08</td><td>4.11</td><td>9.49E-12</td><td>NDE</td><td></td><td>2.25</td><td>3.91E-04</td><td>NDE</td><td></td></tr><tr><td>Z4583</td><td>nanA</td><td>N-acetylneuraminate lyase</td><td>7.79</td><td>1.62E-18</td><td>5.14</td><td>2.72E-08</td><td>3.55</td><td>2.07E-04</td><td>5.90</td><td>8.61E-11</td><td>NDE</td><td></td></tr><tr><td>Z4582</td><td>nanT</td><td>Sialic acid transporter</td><td>6.32</td><td>1.03E-13</td><td>3.85</td><td>1.80E-05</td><td>2.30</td><td>1.71E-02</td><td>5.18</td><td>2.28E-09</td><td>NDE</td><td></td></tr></tbody></table>
<table><caption>**Table 3 RNA-Seq data (Continued)**</caption><thead><tr><th rowspan="3">Gene tag</th><th rowspan="3">Gene</th><th rowspan="3">Gene product</th><th colspan="6">3 hours incubation</th><th colspan="6">6 hours incubation</th></tr><tr><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Rumen</th></tr><tr><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th></tr></thead><tbody><tr><td>Z3463</td><td>napA</td><td>Nitrate reductase (catalytic subunit)</td><td>2.78</td><td>2.25E-07</td><td>NDE</td><td></td><td>2.66</td><td>1.03E-06</td><td>NDE</td><td></td><td>2.77</td><td>3.14E-07</td></tr><tr><td>Z3460</td><td>napB</td><td>Cytochrome v-type protein NapB</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.35</td><td>4.67E-03</td><td>2.24</td><td>7.15E-03</td><td>NDE</td><td></td></tr><tr><td>Z3461</td><td>napH</td><td>Ferredoxin-type protein NapH</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.33</td><td>3.08E-03</td><td>NDE</td><td></td><td>2.15</td><td>6.29E-03</td></tr><tr><td>Z2001</td><td>narG</td><td>Nitrate reductase (α subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>4.50</td><td>1.93E-23</td><td>NDE</td><td></td><td>4.97</td><td>1.46E-28</td></tr><tr><td>Z2002</td><td>narH</td><td>Nitrate reductase (β subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.97</td><td>2.34E-15</td><td>NDE</td><td></td><td>3.36</td><td>1.58E-19</td></tr><tr><td>Z2004</td><td>narI</td><td>Nitrate reductase (γ subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.95</td><td>8.43E-11</td><td>NDE</td><td></td><td>3.30</td><td>2.02E-13</td></tr><tr><td>Z2003</td><td>narJ</td><td>Nitrate reductase molybdenum cofactor</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.23</td><td>3.23E-11</td><td>NDE</td><td></td><td>3.67</td><td>2.51E-14</td></tr><tr><td>Z4726</td><td>nirB</td><td>Nitrite reductase (large subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.06</td><td>9.01E-07</td></tr><tr><td>Z5671</td><td>nrfC</td><td>Nitrite reductase (formate-dependent)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>5.35</td><td>1.49E-13</td><td>2.63</td><td>8.52E-04</td><td>5.20</td><td>7.59E-13</td></tr><tr><td>Z5673</td><td>nrfE</td><td>Cytochrome c-type biogenesis protein NrfE</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.43</td><td>2.08E-06</td><td>NDE</td><td></td><td>3.35</td><td>3.59E-06</td></tr><tr><td>Z5674</td><td>nrfF</td><td>Formate-dependent nitrite reductase complex subunit NrfF</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.02</td><td>1.58E-04</td><td>NDE</td><td></td><td>3.06</td><td>1.23E-04</td></tr><tr><td>Z5675</td><td>nrfG</td><td>Formate-dependent nitrite reductase complex subunit NrfG</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.68</td><td>4.04E-05</td><td>NDE</td><td></td><td>2.78</td><td>1.73E-05</td></tr><tr><td>Z3541</td><td>nuoH</td><td>NADH-quinone oxidoreductase subunit H</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.45</td><td>5.24E-04</td><td>NDE</td><td></td></tr><tr><td>Z3540</td><td>nuoI</td><td>NADH-quinone oxidoreductase subunit I</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.04</td><td>5.31E-03</td><td>2.72</td><td>9.55E-05</td><td>NDE</td><td></td></tr><tr><td>Z3539</td><td>nuoJ</td><td>NADH-quinone oxidoreductase subunit J</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.02</td><td>5.51E-03</td><td>2.82</td><td>1.86E-05</td><td>NDE</td><td></td></tr><tr><td>Z3538</td><td>nuoK</td><td>NADH-quinone oxidoreductase subunit K</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.89</td><td>4.64E-06</td><td>NDE</td><td></td></tr><tr><td>Z3537</td><td>nuoL</td><td>NADH-quinone oxidoreductase subunit L</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.76</td><td>1.98E-05</td><td>NDE</td><td></td></tr><tr><td>Z3536</td><td>nuoM</td><td>NADH-quinone oxidoreductase subunit M</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.04</td><td>1.26E-06</td><td>NDE</td><td></td></tr><tr><td>Z4758</td><td>pckA</td><td>Phosphoenolpyruvate carboxykinase</td><td>3.54</td><td>2.78E-10</td><td>2.60</td><td>8.68E-06</td><td>2.86</td><td>6.67E-07</td><td>3.44</td><td>1.30E-09</td><td>NDE</td><td></td></tr><tr><td>Z3401</td><td>preT</td><td>Dihydropyrimidine dehydrogenase</td><td>2.94</td><td>2.12E-06</td><td>2.48</td><td>1.19E-04</td><td>6.46</td><td>4.01E-28</td><td>2.32</td><td>3.81E-04</td><td>5.40</td><td>9.99E-20</td></tr><tr><td>Z3980</td><td>proW</td><td>Glycine betaine/proline transport system</td><td>2.36</td><td>8.30E-03</td><td>NDE</td><td></td><td>NDE</td><td></td><td>-2.58</td><td>4.47E-03</td><td>NDE</td><td></td></tr><tr><td>Z3981</td><td>proX</td><td>Glycine betaine/proline transport system</td><td>2.74</td><td>6.12E-04</td><td>NDE</td><td></td><td>NDE</td><td></td><td>-2.1</td><td>1.24E-02</td><td>NDE</td><td></td></tr><tr><td>Z3653</td><td>ptsN</td><td>PTS system fructose-specific IIB-like component</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.55</td><td>4.69E-03</td></tr><tr><td>Z5250</td><td>rbsA</td><td>Ribose transport system ATP-binding protein</td><td>6.24</td><td>4.46E-23</td><td>5.83</td><td>7.52E-20</td><td>2.84</td><td>2.63E-05</td><td>5.80</td><td>9.53E-20</td><td>4.60</td><td>9.20E-13</td></tr><tr><td>Z5252</td><td>rbsB</td><td>Ribose transport system substrate-binding protein</td><td>4.23</td><td>1.81E-08</td><td>2.16</td><td>9.34E-03</td><td>NDE</td><td></td><td>2.59</td><td>1.10E-03</td><td>NDE</td><td></td></tr><tr><td>Z5251</td><td>rbsC</td><td>Ribose transport system permease protein</td><td>6.49</td><td>5.50E-29</td><td>5.18</td><td>2.95E-18</td><td>4.04</td><td>1.99E-11</td><td>6.28</td><td>1.33E-26</td><td>5.00</td><td>3.27E-17</td></tr><tr><td>Z5249</td><td>rbsD</td><td>D-ribose pyranase</td><td>6.17</td><td>2.40E-23</td><td>5.28</td><td>5.82E-17</td><td>3.40</td><td>1.91E-07</td><td>6.10</td><td>2.29E-22</td><td>4.81</td><td>3.25E-14</td></tr><tr><td>Z5253</td><td>rbsK</td><td>Ribokinase</td><td>3.82</td><td>6.03E-12</td><td>2.13</td><td>3.67E-04</td><td>NDE</td><td></td><td>2.81</td><td>9.77E-07</td><td>NDE</td><td></td></tr><tr><td>Z2857</td><td>sdaA</td><td>L-serine dehydratase</td><td>2.18</td><td>2.59E-15</td><td>NDE</td><td></td><td>2.46</td><td>3.92E-19</td><td>2.47</td><td>4.52E-19</td><td>2.25</td><td>3.44E-16</td></tr><tr><td>Z4114</td><td>sdaB</td><td>L-serine dehydratase (deaminase)</td><td>4.36</td><td>2.42E-18</td><td>2.93</td><td>1.87E-08</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4113</td><td>sdaC</td><td>Serine transporter</td><td>4.59</td><td>1.85E-30</td><td>3.72</td><td>8.52E-20</td><td>NDE</td><td></td><td>2.39</td><td>1.43E-08</td><td>NDE</td><td></td></tr><tr><td>Z0877</td><td>sdhA</td><td>Succinate dehydrogenase (flavoprotein subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.12</td><td>5.16E-03</td><td>NDE</td><td></td></tr><tr><td>Z0878</td><td>sdhB</td><td>Succinate dehydrogenase (iron-sulfur subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.10</td><td>1.39E-02</td><td>NDE</td><td></td></tr><tr><td>Z0876</td><td>sdhD</td><td>Succinate dehydrogenase (membrane anchor subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.04</td><td>2.05E-03</td><td>NDE</td><td></td></tr><tr><td>Z0880</td><td>sucA</td><td>2-oxoglutarate dehydrogenase E1 component</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.60</td><td>6.80E-04</td><td>NDE</td><td></td></tr><tr><td>Z0881</td><td>sucB</td><td>2-oxoglutarate dehydrogenase E2 component</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.91</td><td>3.11E-04</td><td>NDE</td><td></td></tr><tr><td>Z0882</td><td>sucC</td><td>Succinyl-CoA synthetase (β subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.93</td><td>2.18E-04</td><td>NDE</td><td></td></tr></tbody></table>
<table><caption>**Table 3 RNA-Seq data (Continued)**</caption><thead><tr><th rowspan="3">Gene tag</th><th rowspan="3">Gene</th><th rowspan="3">Gene product</th><th colspan="4">3 hours incubation</th><th colspan="4">6 hours incubation</th><th colspan="2"></th></tr><tr><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Small intestine</th><th colspan="2">Rectum</th><th colspan="2">Rumen</th></tr><tr><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th><th>Log2FC</th><th>q-value</th></tr></thead><tbody><tr><td>Z0883</td><td>sucD</td><td>Succinyl-CoA synthetase (α subunit)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.89</td><td>1.78E-04</td><td>NDE</td><td></td></tr><tr><td>Z5501</td><td>talC</td><td>Fructose 6-phosphate aldolase 2</td><td>2.45</td><td>3.24E-05</td><td>NDE</td><td></td><td>2.61</td><td>1.21E-05</td><td>NDE</td><td></td><td>2.16</td><td>3.41E-04</td></tr><tr><td>Z4469</td><td>tdcb</td><td>Threonine dehydratase</td><td>6.04</td><td>1.04E-09</td><td>NDE</td><td></td><td>6.43</td><td>7.98E-11</td><td>NDE</td><td></td><td>4.56</td><td>8.29E-06</td></tr><tr><td>Z4468</td><td>tdcc</td><td>Threonine/serine transporter</td><td>4.99</td><td>3.98E-07</td><td>NDE</td><td></td><td>5.02</td><td>3.93E-07</td><td>NDE</td><td></td><td>3.14</td><td>2.54E-03</td></tr><tr><td>Z4467</td><td>tdcd</td><td>Propionate kinase</td><td>5.22</td><td>2.18E-07</td><td>NDE</td><td></td><td>5.89</td><td>4.20E-09</td><td>2.29</td><td>3.87E-02</td><td>4.22</td><td>4.89E-05</td></tr><tr><td>Z4466</td><td>tdce</td><td>Formate C-acetyltransferase</td><td>4.54</td><td>3.69E-05</td><td>NDE</td><td></td><td>5.05</td><td>4.69E-06</td><td>2.43</td><td>4.10E-02</td><td>3.79</td><td>8.59E-04</td></tr><tr><td>Z5203</td><td>tnaA</td><td>Tryptophanase</td><td>3.16</td><td>6.15E-03</td><td>NDE</td><td></td><td>4.73</td><td>2.43E-05</td><td>2.48</td><td>3.88E-02</td><td>3.62</td><td>1.69E-03</td></tr><tr><td>Z5204</td><td>tnaB</td><td>Low affinity tryptophan permease</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.28</td><td>2.20E-05</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4444</td><td>uxaA</td><td>Altronate hydrolase</td><td>5.41</td><td>1.05E-13</td><td>3.97</td><td>1.52E-07</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z2184</td><td>uxaB</td><td>Tagaturonate reductase</td><td>4.08</td><td>9.81E-13</td><td>2.19</td><td>3.95E-04</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4445</td><td>uxaC</td><td>Glucuronate isomerase</td><td>6.19</td><td>2.40E-19</td><td>4.89</td><td>4.88E-12</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z5920</td><td>uxuA</td><td>Mannonate hydrolase</td><td>7.50</td><td>3.08E-25</td><td>5.27</td><td>1.83E-12</td><td>2.86</td><td>2.82E-04</td><td>2.46</td><td>1.98E-03</td><td>NDE</td><td></td></tr><tr><td>Z5921</td><td>uxuB</td><td>Mannonate oxidoreductase</td><td>6.68</td><td>2.40E-22</td><td>4.21</td><td>4.35E-09</td><td>2.41</td><td>1.52E-03</td><td>NDE</td><td></td><td>NDE</td><td></td></tr><tr><td>Z4990</td><td>xylA</td><td>Xylose isomerase</td><td>NDE</td><td></td><td>3.45</td><td>4.67E-07</td><td>3.22</td><td>2.57E-06</td><td>4.60</td><td>3.99E-12</td><td>2.04</td><td>4.53E-03</td></tr><tr><td>Z4989</td><td>xylB</td><td>Xylulokinase</td><td>NDE</td><td></td><td>2.02</td><td>5.33E-03</td><td>NDE</td><td></td><td>2.85</td><td>2.43E-05</td><td>NDE</td><td></td></tr><tr><td>Z4991</td><td>xylF</td><td>D-xylose transport system (substrate-binding protein)</td><td>NDE</td><td></td><td>4.83</td><td>3.40E-09</td><td>3.48</td><td>2.93E-05</td><td>5.77</td><td>3.03E-13</td><td>2.54</td><td>3.06E-03</td></tr><tr><td>Z4992</td><td>xylG</td><td>D-xylose transport system (ATP-binding protein)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.08</td><td>4.14E-03</td><td>NDE</td><td></td></tr><tr><td>Z4993</td><td>xylH</td><td>D-xylose transport system (permease protein)</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.03</td><td>1.36E-03</td><td>2.07</td><td>1.01E-03</td><td>NDE</td><td></td></tr><tr><td>Z5717</td><td>yjdE</td><td>Arginine:agmatine antiporter</td><td>2.43</td><td>8.29E-03</td><td>NDE</td><td></td><td>3.76</td><td>1.94E-05</td><td>NDE</td><td></td><td>3.19</td><td>3.64E-04</td></tr><tr><td>Z0461</td><td>z0461</td><td>Fucose transporter (optimal import)</td><td>2.66</td><td>1.23E-03</td><td>3.00</td><td>3.12E-04</td><td>2.33</td><td>6.92E-03</td><td>4.75</td><td>1.65E-09</td><td>NDE</td><td>3.61E-04</td></tr><tr><td>Z1245</td><td>z1245</td><td>Putative transporter</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.02</td><td>1.93E-03</td><td>NDE</td><td></td><td>2.93</td><td>9.54E-05</td></tr><tr><td>Z2575</td><td>z2575</td><td>Tat-targeted selenite reductase subunit YnfE</td><td>2.86</td><td>2.35E-05</td><td>NDE</td><td></td><td>2.97</td><td>1.35E-05</td><td>NDE</td><td></td><td>2.68</td><td>2.24E-04</td></tr><tr><td>Z2576</td><td>z2576</td><td>Tat-targeted selenite reductase subunit YnfF</td><td>2.55</td><td>1.58E-04</td><td>NDE</td><td></td><td>3.31</td><td>7.19E-07</td><td>NDE</td><td>1.65E-09</td><td>2.57</td><td>3.61E-04</td></tr><tr><td>Z2577</td><td>z2577</td><td>Tat-targeted selenite reductase subunit YnfG</td><td>2.75</td><td>2.04E-04</td><td>2.04</td><td>1.04E-02</td><td>3.44</td><td>3.00E-06</td><td>NDE</td><td></td><td>2.78</td><td>5.41E-04</td></tr><tr><td>Z3394</td><td>z3394</td><td>MFS transporter</td><td>NDE</td><td></td><td>2.75</td><td>9.22E-03</td><td>2.42</td><td>2.39E-02</td><td>4.54</td><td></td><td>3.92</td><td>4.73E-03</td></tr><tr><td>Z3715</td><td>z3715</td><td>Ethanolamine utilization cobalamin adenosyltransferase</td><td>NDE</td><td></td><td>NDE</td><td></td><td>3.94</td><td>1.25E-03</td><td>NDE</td><td>6.32E-05</td><td>4.28</td><td>4.33E-11</td></tr><tr><td>Z4207</td><td>z4207</td><td>Xanthine dehydrogenase iron-sulfur-binding subunit</td><td>NDE</td><td></td><td>NDE</td><td></td><td>2.86</td><td>2.14E-06</td><td>NDE</td><td></td><td>2.15</td><td>1.85E-04</td></tr><tr><td>Z4464</td><td>z4464</td><td>L-serine dehydratase</td><td>4.15</td><td>1.51E-04</td><td>NDE</td><td></td><td>4.25</td><td>1.30E-04</td><td>NDE</td><td></td><td>3.22</td><td>1.87E-06</td></tr><tr><td>Z4875</td><td>z4875</td><td>PTS system galactitol-specific IIA component</td><td>2.44</td><td>8.99E-06</td><td>NDE</td><td></td><td>3.84</td><td>2.25E-12</td><td>2.36</td><td>1.45E-05</td><td>3.61</td><td>3.91E-09</td></tr><tr><td>Z4876</td><td>z4876</td><td>PTS system galactitol-specific IIB component</td><td>2.98</td><td>2.87E-07</td><td>NDE</td><td></td><td>3.79</td><td>6.23E-11</td><td>2.68</td><td></td><td>3.43</td><td></td></tr><tr><td>Z4877</td><td>z4877</td><td>PTS system galactitol-specific IIC component</td><td>2.22</td><td>2.94E-05</td><td>NDE</td><td></td><td>2.37</td><td>6.12E-06</td><td>NDE</td><td>3.13E-06</td><td>NDE</td><td>7.53E-05</td></tr></tbody></table>

NDE not differentially expressed

intestine contents (except *mglC* during the late-exponential growth phase). Importantly, the genes encoding the fuco-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 re-
sults 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 rec-
tum content. None of the *eut* genes were found differen-
tially expressed in EHEC during the mid-exponential
growth phase in small intestine and rectum contents.

Genes up-regulated relative to C₄-dicarboxylate transporters

The transport of common C₄-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 C₄-dicarboxylates) was up-regulated
in EHEC during incubation in all the DCs tested, includ-
ing 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 re-
spiratory pathways activated by EHEC EDL933 in the
three DCs was undertaken. As shown in Table 3, com-
mon and distinct respiratory systems were induced by
the bacteria according to the bovine DCs. The genes
fdnG and fdoHI (formate dehydrogenase-N and -I re-
spectively), 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 oxi-
doreductase) was only activated in rectum content ex-
cept 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 intes-
tinal mucus, is a filamentous glycoprotein rich in oligo-
saccharides (up to 80% of the mucin biomass) and
represents a niche for bacterial colonization of the intes-
tine [23, 24]. Indeed, mucus is constantly released into
the mammal GIT and degraded by a subset of com-
mensal 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 com-
partments 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 FusK/R 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–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 mammalian 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 (H₂, 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⁻¹ 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 Na₂CO₃, 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 O₂-free CO₂-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 sensitivity was confirmed by spotting each pooled DC (100 µL) on Sorbitol Mac Conkey (SMAC) agar plates supplemented with 100 µg mL⁻¹ 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 O₂-free CO₂-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⁻¹), 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 ≈ 10⁴ bacteria mL⁻¹ 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⁻¹) 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) quantifica-
tion 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 deter-
mined by gas chromatography. For succinate quantification
in filtered DCs, 200 µL of ZnSO₄ 5% and Ba (OH)₂ 3 M
were added to 800 µL of supernatant before quantification
determined by high-performance liquid chromatography.
All the analysis were performed by AFYREN INVEST-
MENT (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), MgSO₄ (1 mM), CaCl₂ (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 ≈ 10⁸ cells mL⁻¹ before inoculating
rumen content, ≈ 10⁶ cells mL⁻¹ before inoculating small
intestine and rectum contents and ≈ 10⁷ cells mL⁻¹ be-
fore 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 incu-
bation, 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 resus-
pended 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 vol-
ume of 10% sodium dodecyl sulfate, and 3.5 µL
β-mercaptoethanol (PROLABO, France). The lysis of
bacterial cells was performed using a FastPrep®-24 in-
struments (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 (TGGT

TGATGACGAAGAGCTG and GCTCTGGTTCCGGA ATGTAG). 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]. Dif-
ferential gene expression was then performed using
DESeq2 version 1.12.4 [52] with R version 3.3.2 [53] fol-
lowing the standard workflow. The DESeq2 method in-
ternally 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 re-
ported 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 compar-
ing 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 iden-
tify significant differences in gene expression (*p* < 0.05 was
considered as significantly different). The statistical ana-
lysis 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.

**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 assess-
ment, 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)

**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; Rif⁹: Rifampicin-resistant; RNA-seq: RNA sequencing; SCFA: Short chain fatty acid; STEC: Shiga toxin producing *Escherichia coli*

**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.

**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.

**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.

**Author details**

¹Université Clermont Auvergne, INRA, UMR 0454 MEDIS, F-63000 Clermont-Ferrand, France. ²Plateforme Bioinformatique Toulouse, Midi-Pyrénées UBIAC, Auzeville, Castanet-Tolosan, France. ³INRA, US 1426, GeT-PlaGe, Genotoul, Castanet-Tolosan, France. ⁴Present address : Institut National de Police Scientifique - Laboratoire de Police Scientifique de Marseille, Marseille, France. ⁵Present address : Ifremer, UMR 241 EIO, Tahiti, French Polynesia.

Received: 13 July 2018 Accepted: 15 October 2018

Published online: 23 October 2018

**References**

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.

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

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

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.

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.

6. Saeedi P, Yazdanfarast 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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. 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.

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.

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.
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.

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.

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. 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.

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.

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.

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. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet J.* 2011;17(1):10-2.

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.

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

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.

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.

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. 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.

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.

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 part 2:D492-6.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_9.png" /></figure>