Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Changes in bacterial community composition of Escherichia coli O157:H7 super-shedder cattle occur in the lower intestine

  • Rahat Zaheer ,

    Contributed equally to this work with: Rahat Zaheer, Eric Dugat-Bony, Devon Holman

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Eric Dugat-Bony ,

    Contributed equally to this work with: Rahat Zaheer, Eric Dugat-Bony, Devon Holman

    Affiliation UMR Génie et Microbiologie des Procédés Alimentaires, AgroParisTech, INRA, Université Paris-Saclay, Thiverval-Grignon, France

    ⨯
  • Devon Holman ,

    Contributed equally to this work with: Rahat Zaheer, Eric Dugat-Bony, Devon Holman

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Elodie Cousteix,

    Affiliation Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, Canada

    ⨯
  • Yong Xu,

    Affiliation Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, Canada

    ⨯
  • Krysty Munns,

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Lorna J. Selinger,

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Rutn Barbieri,

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Trevor Alexander,

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • Tim A. McAllister ,

    tim.mcallister@agr.gc.ca (TAM); selibl@uleth.ca (LBS)

    Affiliation Agriculture and Agri-Food Canada, Lethbridge Research and Development Centre, Lethbridge, Alberta, Canada

    ⨯
  • L. Brent Selinger

    tim.mcallister@agr.gc.ca (TAM); selibl@uleth.ca (LBS)

    Affiliation Department of Biological Sciences, University of Lethbridge, Lethbridge, Alberta, Canada

    ⨯

Correction

13 Feb 2017: Zaheer R, Dugat-Bony E, Holman DB, Cousteix E, Xu Y, et al. (2017) Correction: Changes in bacterial community composition of Escherichia coli O157:H7 super-shedder cattle occur in the lower intestine. PLOS ONE 12(2): e0172387. https://doi.org/10.1371/journal.pone.0172387 View correction

Abstract

Escherichia coli O157:H7 is a foodborne pathogen that colonizes ruminants. Cattle are considered the primary reservoir of E. coli O157:H7 with super-shedders, defined as individuals excreting > 104 E. coli O157:H7 CFU g-1 feces. The mechanisms leading to the super-shedding condition are largely unknown. Here, we used 16S rRNA gene pyrosequencing to examine the composition of the fecal bacterial community in order to investigate changes in the bacterial microbiota at several locations along the digestive tract (from the duodenum to the rectal-anal junction) in 5 steers previously identified as super-shedders and 5 non-shedders. The overall bacterial community structure did not differ by E. coli O157:H7 shedding status; but several differences in the relative abundance of taxa and OTUs were noted between the two groups. The genus Prevotella was most enriched in the non-shedders while the genus Ruminococcus and the Bacteroidetes phylum were notably enriched in the super-shedders. There was greater bacterial diversity and richness in samples collected from the lower- as compared to the upper gastrointestinal tract (GI). The spiral colon was the only GI location that differed in terms of bacterial diversity between super-shedders and non-shedders. These findings reinforced linkages between E. coli O157:H7 colonization in cattle and the nature of the microbial community inhabiting the digestive tract of super-shedders.

Citation: Zaheer R, Dugat-Bony E, Holman D, Cousteix E, Xu Y, Munns K, et al. (2017) Changes in bacterial community composition of Escherichia coli O157:H7 super-shedder cattle occur in the lower intestine. PLoS ONE 12(1): e0170050. https://doi.org/10.1371/journal.pone.0170050

Editor: Hauke Smidt, Wageningen University, NETHERLANDS

Received: June 5, 2016; Accepted: December 28, 2016; Published: January 31, 2017

Copyright: © 2017 Zaheer et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All raw sequence data described in this study are available through NCBI׳s short-read archives (BioProject PRJNA258277, SRA identifier SRP047129). Additional data has also been provided in the form of Supporting Information files.

Funding: This research work was conducted with the funding from the 'Alberta Livestock and Meat Agency Ltd., Genome Alberta as well as the Abase program of Agriculture and Agri-Food Canada. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Escherichia coli O157:H7 is the serotype responsible for the majority of human enterohemorrhagic E. coli (EHEC) infections in most industrial countries [1]. Cattle are considered the primary reservoir of this pathogen although other animals may be carriers [2]. A number of epidemiological studies have shown that up to 30% of feedlot cattle shed E. coli O157:H7 in their feces [3, 4, 5, 6]. Due to public health concerns and the role of cattle in disseminating E. coli O157:H7, this microorganism has a significant negative impact on the beef production industry.

Heterogeneous shedding among cattle has been observed, where most animals shed only transiently low levels (< 100 CFU g-1 feces) of E. coli O157:H7 while a small number of cattle excrete large quantities of this pathogen or shed lower levels over longer periods of time [7]. The term ‘super-shedder’ was first introduced by Matthews et al., [8, 9] and defined by Chase-Topping et al., [10] in 2008 as an animal that excretes >104 CFU of E. coli O157:H7 per g of feces. Despite their low prevalence within a feedlot (typically <1% of animals), super-shedders are believed to be responsible for the majority of E. coli O157:H7 shed in the feces of cattle and the largest transmitters of this pathogen in feedlots [8, 9, 11].

The primary site of E. coli O157:H7 colonization in cattle has been shown to be the rectal-anal junction (RAJ) [12, 13] and strong associations between RAJ colonization and super-shedder and/or long-term shedder status have been established [14, 15]. The gallbladder has also been suggested as a possible niche for E. coli O157:H7 in cattle [16]. Since cattle are generally asymptomatic carriers of E. coli O157:H7, this bacterium has been described as a commensal. However, the fact that E. coli O157:H7 can cause intestinal inflammation and induce small mucosal haemorrhages and focal petechiae in the cattle intestine as well as induce innate and adaptive immune responses suggests that it could also be viewed as an opportunistic bovine pathogen [17]. Nonetheless, little is known about the etiology of this infection in cattle. Some factors, such as diet, the presence of specific bacteriophages, and other pathogenic strains have been shown to be associated with E. coli O157:H7 infection in cattle [18, 19, 20, 21]. However, these factors alone do not explain the mechanisms that lead to the development of super-shedders and the differences in carriage status among animals.

In a recent study by our group using 16S rRNA gene pyrosequencing, the fecal bacterial microbiota of E. coli O157:H7 super-shedders differed significantly from that of non-shedders originating from a feedlot in southern Alberta [22]. This finding supports the idea that microbiota composition may play a critical role in the establishment and/or ecology of E. coli O157:H7 within the intestinal tract of cattle. In the present study, we used 16S rRNA gene pyrosequencing to examine changes in bacterial diversity and community structure along the small and large intestine of cattle identified as super-shedders and non-shedders. The present study is aimed at comparing the composition and structure of the gastrointestinal (GI) microbiota at several different GI locations in E. coli O157:H7 super-shedders and non-shedders to determine if previously observed differences in fecal communities were associated with specific regions of the intestinal tract.

Materials and methods

Super-shedder identification and study design

The identification of E. coli O157:H7 super-shedders within the feedlot has been previously described in detail [22, 23]. Briefly, fecal samples were collected from 400 crossbred Charolaise yearling feedlot beef steers with body weights ranging between 500–525 Kg and housed at a single feedlot (Lethbridge, Alberta, Canada). All cattle were fed a barley-grain based finishing diet. E. coli O157:H7 super-shedders (>104 E. coli O157:H7 CFU/g of fecal material) were identified by enumeration on CT-SMAC (Sorbitol MacConkey with Tellurite and Cefixime; Dalynn Biologicals, Calgary, Alberta, Canada) as described by Niu et al., [20]. Three representative non-sorbitol fermenting colonies from each sampling point were confirmed to be E. coli O157 using the E. coli O157 latex test kit (Oxoid Ltd., Basingstoke, Hampshire, UK). Positive agglutination isolates were confirmed by multiplex PCR to test for the presence of genes specific to the O157:H7 serotype [24]. When E. coli was not detectable by plating, duplicate 1 g subsamples of feces were enriched in 9 mL of modified TSB containing novobiocin (20 mg/L; Sigma-Aldrich Canada Co., Oakville, ON, Canada), bile salts (1.5 g/L; BD—Canada, Mississauga, ON, Canada), dipotassium phosphate (1.5 g/L; Sigma-Aldrich Canada Co.) and TSB (30 g/L; BD—Canada) and incubated for 6 h at 37°C. Enriched samples were then subject to immunomagnetic separation (IMS) using anti-E. coli O157:H7 Dynabeads (Invitrogen, Carlsbad, CA) as per manufacturer's instructions. A 50 μL aliquot of bead-bacteria complex was plated onto CT-SMAC (Dalynn Biologicals) and incubated at 37°C for 18 to 24 h. Three non-sorbitol fermenting clear colonies were randomly selected for latex-test and PCR confirmation as described above. Cattle that were negative both by enumeration and IMS were classified as non-shedders.

Five super-shedders (>104 E. coli O157:H7 CFU/g of feces) from two pens and five non-shedders (negative for IMS) contemporary pen mates were purchased and transported to the Lethbridge Research and Development Centre (LeRDC). Steers were kept on the same barley-grain based finishing diet as in the feedlot. The shedding pattern of all steers was monitored daily until the date of slaughter as described above. Steers were slaughtered over a period of 6 days following the abattoir’s schedule with 2 animals (1 shedder + 1 non-shedder) slaughtered on the next day, 4 animals (2 shedder + 2 non-shedder) on the 4th day and 4 animals (2 shedder + 2 non-shedder) on the 6th day of original identification and arrival at LeRDC. A typical captive bolt procedure according to the animal humane handling and slaughter guidelines by the Canadian Food Inspection Agency was used for pre-slaughter stunning by experienced abattoir staff immediately followed by the jugular vein bleeding. Steers were sequentially slaughtered as it was not possible to slaughter and conduct the necessary intestinal dissections from all cattle within a single day. This study was carried out in strict accordance with the recommendations and guidance of the Canadian Council of Animal Care. The protocol was reviewed and approved by the Lethbridge Research Centre Animal Care Committee.

Digesta and tissue sample collection

Within 10 min of slaughter, the gastrointestinal tract (GIT) of each animal was removed and placed on a clean sheet of plastic on the floor of the abattoir. Eight intestinal sections, approximately 20 cm in length, were collected and aseptically removed from each animal from the following 8 locations in anatomical order: duodenum, proximal jejunum, mid jejunum, distal jejunum, cecum, spiral colon, descending colon and rectum at the RAJ. Bilateral ligatures were applied adjacent to the excision sites to minimize external contamination of the tissues with the digesta. During the excision, the work areas and tools were decontaminated with 70% ethanol and rinsed with sterile water to prevent cross-contamination among samples. Longitudinal incisions were made in each ligated GIT section and the solid and/or semisolid digesta were gently removed (quantities strongly varied between the different GIT sections and between animals), collected in 50 ml falcon tubes, immediately flash frozen in liquid nitrogen and stored at -80°C. After a gentle wash with sterile water, small pieces of the rectum tissue at the RAJ were scraped using a spatula and the samples were collected in 15 ml falcon tubes, immediately flash frozen in liquid nitrogen and stored at -80°C. Samples of the distal ileum were not collected as this tissue was designated as “specified risk material” and could not be legally removed from the abattoir.

Metagenomic DNA extraction

Metagenomic DNA was extracted from each digesta and tissue sample using the Repeated Bead Beating + Column protocol (RBB+C), originally described by Yu and Morrison [25]. Briefly, 0.25 g of each sample was subjected to two rounds of bead-beating in a SDS containing buffer with 15 min of heating at 70°C. DNA was then precipitated with isopropanol, treated with RNase A (Qiagen), and purified using the column provided in the QIAmp DNA Stool Mini Kit (Qiagen). DNA samples were quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific).

Pyrosequencing of the 16S rRNA gene

The 16S rRNA gene in each extracted DNA sample was amplified and sequenced at Molecular Research LP (Shallowater, Texas, USA) using bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP), as originally described by Dowd et al., [26]. Briefly, the HotStarTaq Plus Master Mix Kit, together with primers 27F (5’-GAGTTTGATCNTGGCTCAG-3’) and 519R (5’-GTNTTACNGCGGCKGCTG-3’), was used to amplify the hypervariable regions V1 to V3 of the 16S rRNA gene. The PCR program used to generate these amplicons consisted of: initial denaturation at 94°C for 3 min, followed by 28 cycles of 94°C for 30 s, 53°C for 40 s and 72°C for 1 min, with a final extension step at 72°C for 5 min. To enable multiplexing of samples, each sample was indexed with a unique 8-bp barcode. PCR products from all samples were then mixed in equal concentrations and purified using Agencourt Ampure beads (Agencourt Bioscience Corporation, MA, USA). Sequencing of the 16S rRNA gene amplicons was then carried out using a Roche 454 FLX Titanium system (San Francisco, CA, USA) following manufacturer’s instructions.

16S rRNA gene sequencing analysis

All 16S rRNA gene sequences were processed and analyzed within the QIIME software package v. 1.9.1 [27]. Reads were first demultiplexed with the removal of barcodes and primer sequences. Next, sequences were trimmed and quality filtered with the removal of sequences that were: > 550 bp, < 250, had a Phred score of < 30, or had homopolymer runs of > 6 bp. Chimeric sequences were identified and removed using the UCHIME algorithm [28] implemented in USEARCH v. 6.1544 [29]. The chimera-free sequences were then clustered into OTUs (operational taxonomic units) at 97% similarity using the open reference OTU picking protocol in QIIME and the Greengenes 13_5 database [30]. Sequences that did not match OTUs in the Greengenes database were clustered into OTUs using the de novo approach and USEARCH. The UCLUST consensus taxonomy assigner [29] was used to assign taxonomy to OTUs using Greengenes, with a minimum similarity of 0.9 and max accepts of 3. Representative sequences for the OTUs were aligned using PyNast [31]. Singletons, that is those OTUs found only once in the entire dataset, were also removed prior to analysis.

To account for uneven sequencing depths among samples, all samples were randomly subsampled to 3,100 sequences. This necessitated the loss of five samples from the analysis: one from the rectum and two each from the distal jejunum and cecum. The diversity within each sample (alpha diversity) was calculated within QIIME using the Chao1 estimator [32], observed OTUs, Shannon index [33], and PD whole tree [34] metrics. Diversity metrics for each intestinal section and shedding status group were analyzed using a two-way ANOVA implemented in PROC MIXED in SAS (version 9.4; SAS Institute Inc., Cary, NC, USA) followed by pairwise comparisons using Tukey’s HSD (honest significant difference; P < 0.05). The bacterial community structure (beta diversity) was evaluated using weighted UniFrac distances [35] and visualized using PCoA and the program EMPeror [36]. ANOSIM (analysis of similarities) with 999 permutations was used to compare the weighted UniFrac distances. OTUs that were differentially abundant among the eight different intestinal sections as well as between super-shedders and non-shedders were identified using G-test together with the Benjamini-Hochberg false discovery rate (FDR) correction [37] for multiple comparisons (FDR < 0.05).

Linear discriminant analysis effect size (LEfSe) was used to determine which taxa groups differed according to E.coli O157:H7 shedding status and GI section. LEfSe includes the Kruskal-Wallis test to identify different (P < 0.05) taxa among groups of samples followed by linear discriminant analysis which estimates the size of each of these differences [38].

Results

Escherichia coli O157:H7 shedding pattern at the date of slaughter

At the beginning of the experiment, 5 super-shedders (> 104 CFU/g of feces detected by direct plating on CT-SMAC) and 5 non-shedders (no detection neither by direct plating on CT-SMAC nor by enrichment using immuno-magnetic separation) were identified. The shedding pattern of the five super-shedders was followed twice daily until slaughter as described by Munns et al., [23]. Only one steer, ID 310, met the definition of a super-shedder at the time of slaughter. However, E. coli O157:H7 was consistently isolated in fecal samples until slaughter from all steers originally identified as super-shedders. Remaining steers were slaughtered 4 (IDs 274 and 287) and 6 days (IDs 294 and 299) after arrival at LeRDC. E. coli O157:H7 was not detected in the five non-shedders during the study or at slaughter. The non-shedders were slaughtered along with the shedders on days 2 (ID 165), 4 (IDs 108 and 745) and 6 (IDs 152 and 242).

Bovine gastrointestinal tract microbiota

Eighty-four 16S rRNA gene libraries, corresponding to the 84 digesta and tissue samples, were sequenced using bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). A total of 755,891 sequences with an average length of 458 bp were clustered into 18,542 OTUs after the removal of singletons and following random subsampling of each sample to 3,100 sequences. There were 21 phyla and 120 genera detected among all intestinal samples. Overall, the ten most abundant OTUs accounted for 27.9% of sequences (data not shown). Independent of shedding status or intestinal section sampled, Firmicutes and Bacteoidetes comprised greater than 70% of all sequences (Table A in S1 File). There was considerably more variation among genera as only Turcibacter and Prevotella had an average relative abundance of greater than 4% among all samples (Table B in S1 File). The bacterial diversity in the bovine GI tract was highly variable and was strongly dependent on the GI site that the sample was derived from (Fig 1). There were no OTUs that were shared among all samples, although four OTUs were shared by 90% of the samples (Table C in S1 File). Turcibacter, SMB83, and Clostridium were the only genera that were detected in greater than 90% of all samples. Within the upper GI tract samples, only the genus Mogibacterium was found in each sample, while in the lower GI tract, Turcibacter, Clostridium, and Ruminococcus were common to all samples.

thumbnail
Download:
Fig 1. Box plots displaying the bacterial diversity and richness for E. coli O157:H7 shedder (S) and non-shedder (NS) cattle at each gastrointestinal section (n = 10).

A) Chao1, B) number of OTUs, C) PD whole tree, and D) Shannon index (median ± interquartile range).

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

Effect of E. coli O157:H7 shedding status on the bovine GI tract microbiota

Differentially abundant phyla and genera between super-shedders and non-shedders were identified using LEfSe (LDA score [log10] > 3.0; P < 0.05) for pooled data of gastrointestinal samples. The non-shedders were found to have higher relative abundance of Bacteroidetes, Fibrobacteres, and Cyanobacteria phyla (LDA score [log10] > 3.5). At the genus level, Prevotella and Treponema were notably more abundant in non-shedders, while Ruminococcocus, Selenomonas, Campylobacter, and Streptococcus genera were enriched in super-shedders (Fig 2). Of these last four genera, only Ruminococcus had an overall relative abundance of greater than 0.05% (Table B in S1 File). When the GI section that the sample was derived from was included in the analysis, no difference was noted among the five and 20 most relatively abundant phyla and genera, respectively. This most likely reflects the high degree of variability, particularly at the genus-level among individual samples (Tables A and B in S1 File).

thumbnail
Download:
Fig 2. Differentially abundant genera associated with E. coli O157:H7 super-shedders as assessed using linear discriminant analysis effect size (LEfSe).

Yellow bars indicate genera enriched in the non-shedder samples (n = 39) and green bars represent genera enriched in the shedder samples (n = 41). Only genera with a LDA score [log10] > 3.0 are displayed.

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

There were 12 OTUs that were differentially abundant between super-shedders and non-shedders in the lower GI samples (Table D in S1 File). Eight of these OTUs were more abundant in the non-shedder cohort and four in the super-shedder group. Similarly, in the upper GI samples, 12 OTUs were also differentially abundant between super-shedders and non-shedders with eight of these more abundant in the non-shedders (Table E in S1 File). Notably, two OTUs (514059 and 588315) were more abundant in super-shedders in both the upper and lower GI tract. These two OTUs were both classified as Ruminococcus, a genus which was also found to be more abundant in super-shedders using LEfSe. However, at the individual animal level, the enrichment of Ruminococcus within the super-shedding group appeared to be the result of one steer which had a relatively high proportion (> 26%) of this genus within the duodenum, proximal, mid, and distal jejunum, and the descending colon.

Overall, in terms of bacterial diversity, none of the diversity or richness metrics differed significantly based on E. coli O157:H7 shedding status when lower and upper GI tract samples were combined (P > 0.05). However, there were significantly more OTUs and a higher Chao1 value in spiral colon samples from super-shedders compared with non-shedders (Fig 1A and 1B; P < 0.05). Although, not significant, the Shannon index and the phylogenetic diversity (PD whole tree) was also reduced in the non-shedders at this GI location (Fig 1C and 1D). A comparison of only the lower GI tract samples also revealed that super-shedders had significantly greater Chao1 values, OTU richness, and phylogenetic diversity than non-shedders (P < 0.05; data not shown). The bacterial community structure of samples from super-shedders and non-shedders was compared using weighted UniFrac distances. Samples did not cluster by shedding status when all samples were included in the analysis (Fig 3A; P > 0.05; R-value = 0.033).

thumbnail
Download:
Fig 3. Principal coordinate analysis plots (PCoA) plots of the weighted UniFrac distances by E. coli O157:H7 shedding status.

A) all samples and B) lower GI only. The percent variation explained by the principal coordinates is indicated on the axes.

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

When upper GI samples were removed from the analysis, there was a statistically significant difference between shedders and non-shedders (P = 0.014), but the associated R-value was relatively small (0.073) indicating that the two groups of samples were not well separated (Fig 3B). There were also 34 OTUs found in 90% of the lower GI samples from super-shedders (Table F in S1 File) while in the non-shedder lower GI samples, only seven OTUs were shared among 90% of the samples (Table G in S1 File). As the rectum and rectal tissue are most often believed to be the site of E. coli O157:H7 colonization, super-shedder and non-shedder samples taken from these two sites were compared using weighted UniFrac distances. In addition, due to the fact that spiral colon super-shedder and non-shedder samples differed based on alpha diversity metrics, this location was also isolated and analyzed using the same methods. However, the bacterial community structure at these three locations did not differ significantly based on E. coli O157:H7 shedding status (weighted UniFrac distances; data not shown; P > 0.05).

Gut microbiota composition by gastrointestinal location

Whether a sample came from the upper (duodenum, proximal jejunum, mid-jejunum, distal jejunum) or lower (cecum, spiral colon, descending colon, rectal tissue, rectum) GI tract was the largest determinant of bacterial community composition. When samples were grouped by upper vs. lower GI tract, the upper GI tract samples had fewer OTUs, and lower Chao1, Shannon index, and PD whole tree values than the lower GI tract samples (Fig 1; P < 0.05). There were also 43 genera that were differentially abundant at one of the nine GI sections as assessed by LEfSe (Fig A in S2 File; LDA score [log10]> 3.0; P < 0.05). Among the most relatively abundant genera (> 1% of sequences overall), Turicibacter was enriched in the mid jejunum, Prevotella and Lactobacillus in the duodenum, SMB53 in the distal jejunum, Clostridium in the rectum, and CF231 and Oscillospira in the cecum. The number of OTUs shared by 100% of the samples from each GI section varied from only two in the distal jejunum to 45 OTUs in rectal samples (Table H in S1 File).

Not surprisingly, samples from the upper and lower GI tract formed two distinct clusters when analyzed using weighted UniFrac distances and PCoA (Fig B in S2 File; P < 0.001; R-value = 0.57). When the GI section that the sample was derived from was included in the analysis, the samples still clustered largely by GI tract location (upper and lower) rather than by individual GI section, although some duodenum samples appeared to cluster between the lower GI tract samples and the other upper GI tract samples (Fig C(A) in S2 File; P > 0.05). In the lower GI tract, samples clustered together mostly on the basis of host rather than on anatomical location within the GI tract (Fig C(B) in S2 File; P < 0.001; R-value = 0.67).

Discussion

Changes in bacterial community structure along the GI tract has been well documented in humans [39, 40, 41] and more recently in other animals such as cats [42], horses [43] and dairy cows [44] using fingerprinting methods or 16S rRNA gene clone libraries. These studies highlight a clear longitudinal effect on dominant bacterial populations along the digestive tract of mammals. de Oliveira et al., [45] reported for the first time the use of pyrosequencing to characterize these changes along the entire gastrointestinal tract of one steer, from the rumen to the feces. The GI tract microbiota of dairy cattle has also recently been described based on high-throughput sequencing of the 16S rRNA gene [46]. In the present study, we used a similar approach to characterize changes in bacterial diversity across eight adjacent intestinal sections of ten feedlot steers, from the duodenum to the rectal-anal junction.

We observed that samples clustered largely based on whether they were derived from the upper or lower GI tract (Fig Ca in S2 File). Bacterial diversity and richness were significantly higher in the lower GI tract than in the upper GI tract where the Firmicutes phylum was largely predominant (> 85% of sequences), with the notable exception of the duodenum (Fig 1; Table A in S1 File). With the exception of samples taken from the spiral colon, a number of genera were enriched within the individual GI sections, including some of the most relatively abundant genera (Table B in S1 File). For example, Prevotella was enriched in the duodenum and Turicibacter in the mid jejunum. Prevotella spp., which are acetate producers, are also typically among the most abundant bacteria in the rumen [47, 48] and as the duodenum is the first site of pH neutrality in the GI tract after the abomasum, it seems likely that rumen digesta serves as a source of this genus. Although Turicibacter has been described as being relatively abundant in the cecum [49], colon [50], and feces [51] of feedlot cattle, it has also been reported to be enriched in the ileum and jejunum of dairy cattle [46].

The environmental conditions in the GI tract such as digesta pH (reduced in the upper GI tract), passage time (faster in the upper GI tract) or nutrient availability and composition, are all factors that may account for differences in microbiota along the GI tract. Furthermore, our results suggest that the bacterial microbiota is relatively stable along the bovine lower GI tract with more variability in the upper intestine. This finding is in agreement with previous observations in both cattle [45], humans [41], and pigs [52].

In a previous study, we determined that super-shedders and non-shedder feedlot steers had a unique fecal microbiota [22]. This suggested that there may be a relationship between changes in the bovine GI tract microbiota and colonization by E. coli O157:H7. This bacterium is preferentially detected in the lower GI tract of cattle, and more precisely within the RAJ [14, 15]. Interestingly, the only significant difference in the microbiota of super-shedders and non-shedders along the length of the GI tract was observed in the spiral colon where super-shedders had greater OTU richness (Fig 1). However, the bacterial community structure was not different between the two groups of cattle at this GI location. Although Ruminococcus appeared to be enriched in the GI microbiota of super-shedders based on LEfSe (Fig 2), analysis of individual steers revealed that this difference was largely the result of one steer having a high proportion of Ruminococcus at several anatomical locations within the GI tract.

Prevotella and Treponema genera were both markedly enriched in the GI tract of the non-shedder cattle. Within this cohort, Prevotella was relatively most abundant in the duodenum, cecum, rectum, and rectal tissue (Table B in S1 File). In addition to being relatively abundant in the rumen of cattle, Prevotella spp. are also significant members of the lower GI tract and fecal microbiota of dairy and feedlot cattle [26, 51, 53]. Certain Treponema spp. have been implicated in diseases in cattle such as bovine digital dermatitis (BDD) and it has recently been suggested that the GI tract serves as a reservoir for BDD-associated Treponema spp. [54]. However, this genus is also ubiquitous and relatively abundant within the rumen and lower GI tract of healthy cattle as evidenced from previous microbiome studies [22, 26, 46, 53].

Recent work aimed at describing and understanding the human gut microbiota has demonstrated the importance of bacterial diversity in the resilience of the gut microbiota and its implication in the prevention of numerous diseases [55]. In particular, it is now generally believed that a more diverse microbiota is more stable and less susceptible to invasion by enteropathogens than one with lower diversity and richness [56, 57]. This is known as the colonization resistance mechanism. Indeed, colonization by pathogens responsible for intestinal inflammation is often associated with dysbiosis of the microbiota in humans. In most cases, bacterial richness is reduced in infected vs. healthy individuals. For example, this was observed in patients suffering from diarrhea caused by Clostridium difficile [58, 59], in mice infected with Salmonella enterica [60] or Campylobacter jejuni [61], in TLR5-KO mice infected with adherent-invasive E. coli [62], and in pigs infected with enterotoxigenic Escherichia coli strain F4 [63].

The loss of colonization resistance leading to the potential overgrowth of enteropathogens and to gut inflammation might be due to several possible mechanisms mediated by complex ecological interactions [57]. One of the proposed mechanisms for this observation suggests that virulence factors enable the pathogen to elicit help from the host inflammatory response to gain a competitive growth advantage over resident microbiota [64]. We observed an increase in microbial diversity and richness in the spiral colon and overall in the lower GI tract of super-shedders (Fig 1), an observation that may reflect the uncertainty of the role of E. coli O157:H7 in intestinal inflammation in cattle [65]. The only similar report concerns helminth colonization of the human gut [66]. Perhaps, because intestinal inflammation in cattle is not as severe as in humans, the bovine immune response is different and as a result does not induce a similar perturbation of the microbiota. Further studies should attempt to define host immune responses in the large intestine of super-shedders, particularly given the recent finding that the expression of genes associated with immunity within the rectal tissue of super-shedders is suppressed [67]. The possibility also exists that differences in the bacterial microbiome within the spiral colon of super-shedders and non-shedders may be a reflection of digesta flow differences from the small intestine and cecum. No ileum samples were available for analysis in this study; however the spiral colon samples of super-shedders appeared to group more closely with their respective cecum samples than in non-shedders. In addition, both the PCoA and the alpha diversity analysis by section of the GI tract indicated more variability among the spiral colon samples in the non-shedders than in the super-shedders.

In conclusion, this study provides new insight into the total gut microbiota of cattle shedding E. coli O157:H7. Although the overall bacterial community structure was not altered by E. coli O157:H7 shedding status, several major bacterial genera, such as Prevotella and Treponema, as well as the phylum Bacteroidetes, were differentially abundant in the two groups of cattle. In addition, spiral colon samples from super-shedders had greater bacterial richness than non-shedders. Despite the fact that super-shedding status is usually defined by E. coli O157:H7 shedding concentrations in the rectum, the microbiota of this GI section and that of the rectal tissue samples was not different between super-shedders and non-shedders.

Supporting information

S1 File. Tables supporting information.

Table A in S1 File. Comparisons of the five most relatively abundant phyla by gastrointestinal section and E. coli O157:H7 shedding status.

Table B in S1 File. Comparisons of the 20 most relatively abundant genera by gastrointestinal section and E. coli O157:H7 shedding status.

Table C in S1 File. OTUs found in 90% of all samples (n = 80).

Table D in S1 File. Differentially abundant OTUs between super-shedders (n = 21) and non-shedders (n = 21) in lower GI samples.

Table E in S1 File. Differentially abundant OTUs between super-shedders (n = 20) and non-shedders (n = 18) in upper GI samples.

Table F in S1 File. OTUs found in 90% of all lower GI samples from super-shedding cattle (n = 21).

Table G in S1 File. OTUs found in 90% of all lower GI samples from non-shedding cattle (n = 21).

Table H in S1 File. OTUs found in 100% of the samples taken from each GI section. n refers to the number of samples from each GI section.

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

(PDF)

S2 File. Figures supporting information.

Fig A in S2 File. Differentially abundant genera in each gastrointestinal section as assessed using LEfSE.

Fig B in S2 File. PCoA plot of the weighted UniFrac distances for lower and upper gastrointestinal samples.

Fig C (A) in S2 File. PCoA plots of the weighted UniFrac distances for each gastrointestinal section.

Fig C (B) in S2 File. PCoA plots of the weighted UniFrac distances for each animal for lower GI samples only

https://doi.org/10.1371/journal.pone.0170050.s002

(PDF)

Author Contributions

  1. Conceptualization: TAM LBS.
  2. Data curation: RZ DH.
  3. Formal analysis: RZ DH EDB YX.
  4. Funding acquisition: TAM LBS.
  5. Investigation: RZ DH EDB YX RB LJS TAM LBS.
  6. Methodology: RZ EDB EC KM LJS RB.
  7. Project administration: LBS TAM.
  8. Resources: LBS TAM.
  9. Supervision: TAM LBS TA.
  10. Writing – original draft: EDB RZ TAM.
  11. Writing – review & editing: RZ DH TAM.

References

  1. 1. Armstrong GL, Hollingsworth J, Morris JG Jr. Emerging foodborne pathogens: Escherichia coli O157:H7 as a model of entry of a new pathogen into the food supply of the developed world. Epidemiol Rev. 1996; 18: 29–51. pmid:8877329
  2. 2. Duffy G. Verocytoxigenic Escherichia coli in animal faeces, manures and slurries. J Appl Microbiol. 2003; 94: 94S–103S. pmid:12675941
  3. 3. Callaway TR, Edrington TS, Brabban AD, Keen JE, Anderson RC, Rossman ML, et al. Fecal prevalence of Escherichia coli O157, Salmonella, Listeria, and Bacteriophage Infecting E. coli O157:H7 in feedlot cattle in the Southern Plains region of the United States. Foodborne Pathog Dis. 2006; 3: 234–244. pmid:16972771
  4. 4. Reinstein S, Fox JT, Shi X, Nagaraja TG. Prevalence of Escherichia coli O157:H7 in gallbladders of beef cattle. Appl Environ Microbiol. 2007; 73: 1002–1004. pmid:17158626
  5. 5. Stanford K, Bach SJ, Marx TH, Jones S, Hansen JR, Wallins GL, et al. Monitoring Escherichia coli O157:H7 in inoculated and naturally colonized feedlot cattle and their environment. J Food Prot. 2005; 68: 26–33. pmid:15690800
  6. 6. Stanford K, Croy D, Bach SJ, Wallins GL, Zahiroddini H, McAllister TA. Ecology of Escherichia coli O157:H7 in commercial dairies in southern Alberta. J Dairy Sci. 2005; 88: 4441–4451. pmid:16291636
  7. 7. Ferens WA, Hovde CJ. Escherichia coli O157:H7: animal reservoir and sources of human infection. Foodborne Pathog Dis. 2011; 8: 465–487. pmid:21117940
  8. 8. Matthews L, Low JC, Gally DL, Pearce MC, Mellor DJ, Heesterbeek JA, et al. Heterogeneous shedding of Escherichia coli O157 in cattle and its implications for control. Proc Natl Acad Sci USA 2006; 103: 547–552. pmid:16407143
  9. 9. Matthews L, McKendrick IJ, Ternent H, Gunn GJ, Synge B, Woolhouse ME. Super-shedding cattle and the transmission dynamics of Escherichia coli O157. Epidemiol Infect. 2006; 134: 131–142. pmid:16409660
  10. 10. Chase-Topping M, Gally D, Low C, Matthews L, Woolhouse M. Super-shedding and the link between human infection and livestock carriage of Escherichia coli O157. Nat Rev Microbiol. 2008; 6: 904–912. pmid:19008890
  11. 11. Arthur TM, Keen JE, Bosilevac JM, Brichta-Harhay DM, Kalchayanand N, Shackelford SD. Longitudinal study of Escherichia coli O157:H7 in a beef cattle feedlot and role of high-level shedders in hide contamination. Appl Environ Microbiol. 2009; 75: 6515–6523. pmid:19684164
  12. 12. 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: 93–97. pmid:15640175
  13. 13. 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: 1505–1512. pmid:12595469
  14. 14. 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: 1380–1382. pmid:17189448
  15. 15. Cobbold RN, Hancock DD, Rice DH, Berg J, Stilborn R, Hovde CJ, et al. Rectoanal junction colonization of feedlot cattle by Escherichia coli O157:H7 and its association with supershedders and excretion dynamics. Appl Environ Microbiol. 2007; 73: 1563–1568. pmid:17220263
  16. 16. Stoffregen WC, Pohlenz JF, Dean-Nystrom EA. Escherichia coli 0157:H7 in the gallbladders of experimentally infected calves. J Vet Diagn Invest. 2004; 16: 79–83. pmid:14974852
  17. 17. Vande Walle K, Vanrompay D, Cox E. Bovine innate and adaptive immune responses against Escherichia coli O157:H7 and vaccination strategies to reduce faecal shedding in ruminants. Vet Immunol Immunopathol. 2013: 152: 109–20. pmid:23084625
  18. 18. Callaway TR, Carr MA, Edrington TS, Anderson RC, Nisbet DJ. Diet, Escherichia coli O157:H7, and cattle: a review after 10 years. Curr Issues Mol Biol. 2009; 11: 67–79. pmid:19351974
  19. 19. Chase-Topping ME, McKendrick IJ, Pearce MC, MacDonald P, Matthews L, Halliday J, et al. Risk factors for the presence of high-level shedders of Escherichia coli O157 on Scottish farms. J Clin Microbiol. 2007; 45: 1594–1603. pmid:17360845
  20. 20. Niu YD, McAllister TA, Xu Y, Johnson RP, Stephens TP, Stanford K. Prevalence and impact of bacteriophages on the presence of Escherichia coli O157:H7 in feedlot cattle and their environment. Appl Environ Microbiol. 2009: 75: 1271–1278. pmid:19139243
  21. 21. Carlson BA, Nightingale KK, Mason GL, Ruby JR, Choat WT, Loneragan GH, et al. Escherichia coli O157:H7 strains that persist in feedlot cattle are genetically related and demonstrate an enhanced ability to adhere to intestinal epithelial cells. Appl Environ Microbiol. 2009: 75; 5927–5937. pmid:19617387
  22. 22. Xu Y, Dugat-Bony E, Zaheer R, Selinger L, Barbieri R, Munns K, et al. Escherichia coli O157:H7 Super-shedder and non-shedder feedlot steers harbour distinct fecal bacterial communities. PLoS One 2014: 9; e98115. pmid:24858731
  23. 23. Munns KD, Selinger L, Stanford K, Selinger LB, McAllister TA. Are super-shedder feedlot cattle really super? Foodborne Path Dis. 2014; 11: 329–331. pmid:24484291
  24. 24. Gannon VP, D'Souza S, Graham T, King RK, Rahn K, Read S. Use of the flagellar H7 gene as a target in multiplex PCR assays and improved specificity in identification of enterohemorrhagic Escherichia coli strains. J Clin Microbiol. 1997; 35: 656–662. pmid:9041407
  25. 25. Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques 2004; 36: 808–812. pmid:15152600
  26. 26. Dowd SE, Callaway TR, Wolcott RD, Sun Y, McKeehan T, Hagevoort RG, et al. Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol. 2008; 8: 125. pmid:18652685
  27. 27. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010; 7: 335–336. pmid:20383131
  28. 28. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011; 27: 2194–2200. pmid:21700674
  29. 29. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010; 26: 2460–2461. pmid:20709691
  30. 30. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006; 72: 5069–5072. pmid:16820507
  31. 31. Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 2010; 26: 266–267. pmid:19914921
  32. 32. Chao A. Nonparametric estimation of the number of classess in a population. Scand J Statist. 1984; 11: 265–270.
  33. 33. Shannon C.E. A mathematical theory of communication. Bell System Tech J. 1948; 27: 379–423.
  34. 34. Faith DP. Conservation evaluation and phylogenetic diversity. Biol Con. 1992; 61: 1–10.
  35. 35. Lozupone C, Knight R. UniFrac: a New Phylogenetic Method for Comparing Microbial Communities. Appl Environ Microbiol. 2005; 71: 8228–8235. pmid:16332807
  36. 36. Vázquez-Baeza Y, Pirrung M, Gonzalez A, Knight R. EMPeror: a tool for visualizing high-throughput microbial community data. Gigascience 2013; 2:16. pmid:24280061
  37. 37. Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Royal Stat Soc. Series B (Methodological). 1995; 57: 289–300.
  38. 38. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011; 12: R60. pmid:21702898
  39. 39. Hayashi H, Takahashi R, Nishi T, Sakamoto M, Benno Y. Molecular analysis of jejunal, ileal, caecal and recto-sigmoidal human colonic microbiota using 16S rRNA gene libraries and terminal restriction fragment length polymorphism. J Med Microbiol. 2005; 54: 1093–1101. pmid:16192442
  40. 40. Marteau P, Pochart P, Dore J, Bera-Maillet C, Bernalier A, Corthier G. Comparative study of bacterial groups within the human cecal and fecal microbiota. Appl Environ Microbiol. 2001; 67; 4939–4942. pmid:11571208
  41. 41. Wang M, Ahrne S, Jeppsson B, Molin G. Comparison of bacterial diversity along the human intestinal tract by direct cloning and sequencing of 16S rRNA genes. FEMS Microbiol Ecol. 2005; 54: 219–231. pmid:16332321
  42. 42. Ritchie LE, Steiner JM, Suchodolski JS. Assessment of microbial diversity along the feline intestinal tract using 16S rRNA gene analysis. FEMS Microbiol Ecol. 2008; 66; 590–598. pmid:19049654
  43. 43. Schoster A, Arroyo LG, Staempfli HR, Weese JS. Comparison of microbial populations in the small intestine, large intestine and feces of healthy horses using terminal restriction fragment length polymorphism. BMC Res Notes 2013; 6: 91. pmid:23497580
  44. 44. Frey JC, Pell AN, Berthiaume R, Lapierre H, Lee S, Ha JK, et al. Comparative studies of microbial populations in the rumen, duodenum, ileum and faeces of lactating dairy cows. J Appl Microbiol. 2010; 108; 1982–1993. pmid:19863686
  45. 45. de Oliveira MN, Jewell KA, Freitas FS, Benjamin LA, Totola MR, Borges AC et al. Characterizing the microbiota across the gastrointestinal tract of a Brazilian Nelore steer. Vet Microbiol. 2013; 164: 307–314. pmid:23490556
  46. 46. Mao S, Zhang M, Liu J, Zhu W. Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: membership and potential function. Sci Rep. 2015; 5: 16116. pmid:26527325
  47. 47. Jami E, Mizrahi I. Composition and Similarity of Bovine Rumen Microbiota across Individual Animals. PLoS One 2012; 7: e33306. pmid:22432013
  48. 48. Thoetkiattikul H, Mhuantong W, Laothanachareon T, Tangphatsornruang S, Pattarajinda V, Eurwilaichitr L, et al. Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing. Curr Microbiol. 2013; 67: 130–137. pmid:23471692
  49. 49. Myer PR, Wells JE, Smith TP, Kuehn LA, Freetly HC. Cecum microbial communities from steers differing in feed efficiency. J Anim Sci. 2015; 93: 5327–40. pmid:26641052
  50. 50. Myer PR, Wells JE, Smith TP, Kuehn LA, Freetly HC. Microbial community profiles of the colon from steers differing in feed efficiency. Springerplus 2015; 4: 454. pmid:26322260
  51. 51. Kim M, Kim J, Kuehn LA, Bono JL, Berry ED, Kalchayanand N, et al. Investigation of bacterial diversity in the feces of cattle fed different diets. J Anim Sci. 2014; 92: 683–694. pmid:24352967
  52. 52. Looft T, Allen HK, Cantarel BL, Levine UY, Bayles DO, Alt DP, et al. Bacteria, phages and pigs: the effects of in-feed antibiotics on the microbiome at different gut locations. ISME J. 2014; 8: 1566–1576. pmid:24522263
  53. 53. Callaway TR, Dowd SE, Edrington TS, Anderson RC, Krueger N, Bauer N, et al. Evaluation of bacterial diversity in the rumen and feces of cattle fed different levels of dried distillers grains plus solubles using bacterial tag-encoded FLX amplicon pyrosequencing. J Anim Sci. 1010; 88: 3977–83. pmid:20729286
  54. 54. Sullivan LE, Carter SD, Duncan JS, Grove-White DH, Angell JW, Evans NJ, et al. The gastrointestinal tract as a potential infection reservoir of digital dermatitis-associated treponemes in beef cattle and sheep. Appl Environ Microbiol 2015; 81: 7460–7469. pmid:26276110
  55. 55. Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature 2012; 489: 220–230. pmid:22972295
  56. 56. Pham TA, Lawley TD. Emerging insights on intestinal dysbiosis during bacterial infections. Curr Opin Microbiol. 2014; 17: 67–74. pmid:24581695
  57. 57. Stecher B, Hardt WD. The role of microbiota in infectious disease. Trends Microbiol. 2008; 16: 107–114. pmid:18280160
  58. 58. Antharam VC, Li EC, Ishmael A, Sharma A, Mai V, Rand KH, et al. Intestinal dysbiosis and depletion of butyrogenic bacteria in Clostridium difficile infection and nosocomial diarrhea. J Clin Microbiol. 2013; 51: 2884–92. pmid:23804381
  59. 59. Chang JY, Antonopoulos DA, Kalra A, Tonelli A, Khalife WT, Schmidt TM, et al. Decreased diversity of the fecal microbiome in recurrent Clostridium difficile-associated diarrhea. J Infect Dis. 2008; 197: 435–438. pmid:18199029
  60. 60. Stecher B, Robbiani R, Walker AW, Westendorf AM, Barthel M, Kremer M, et al. Salmonella enterica serovar typhimurium exploits inflammation to compete with the intestinal microbiota. PLoS Biol, 2007; 5: 2177–2189. pmid:17760501
  61. 61. Lone AG, Selinger LB, Uwiera RRE, Xu Y, Inglis GD. Campylobacter jejuni colonization is associated with a dysbiosis in the cecal microbiota of mice in the absence of prominent inflammation. PLoS ONE. 2013; 38: e75325.
  62. 62. Chassaing B, Koren O, Carvalho FA, Ley RE, Gewirtz AT. AIEC pathobiont instigates chronic colitis in susceptible hosts by altering microbiota composition. Gut 2013; 63: 1069–1080. pmid:23896971
  63. 63. Messori S, Trevisi P, Simongiovanni A, Priori D, Bosi P. Effect of susceptibility to enterotoxigenic Escherichia coli F4 and of dietary tryptophan on gut microbiota diversity observed in healthy young pigs. Vet Microbiol. 2013; 162: 173–179. pmid:23021862
  64. 64. Thiennimitr P, Winter SE, Baumler AJ. Salmonella, the host and its microbiota. Curr Opin Microbiol. 2012; 15: 108–114. pmid:22030447
  65. 65. Nart P, Naylor SW, Huntley JF, McKendrick IJ, Gally DL, Low JC. Responses of cattle to gastrointestinal colonization by Escherichia coli O157:H7. Infect Immun. 2008; 76: 5366–5372. pmid:18765741
  66. 66. Lee SC, Tang MS, Lim YAL, Choy SH, Kurtz ZD, Cox LM, et al. Helminth colonization is associated with increased diversity of the gut microbiota. PLoS Negl. Trop. Dis. 2014; 8: e2880. pmid:24851867
  67. 67. Wang O, Liang G, McAllister TA, Plastow G, Stanford K, Selinger LB, et al. Comparative transcriptomic analysis of rectal tissue from beef steers revealed reduced host immunity in Escherichia coli O157:H7 super shedders. PLoS One 2016; 11: e0151284. pmid:26959367