Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Environmental Microbiology | Spotlight

Taxonomic and Functional Compositions of the Small Intestinal Microbiome in Neonatal Calves Provide a Framework for Understanding Early Life Gut Health

Nilusha Malmuthuge, Guanxiang Liang, Philip J. Griebel, Le Luo Guan
Andrew J. McBain, Editor
Nilusha Malmuthuge
aDepartment of Agricultural Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Nilusha Malmuthuge
Guanxiang Liang
aDepartment of Agricultural Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Philip J. Griebel
bSchool of Public Health, University of Saskatchewan, Saskatoon, Saskatchewan, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Le Luo Guan
aDepartment of Agricultural Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew J. McBain
University of Manchester
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.02534-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

A lack of information on the intestinal microbiome of neonatal calves prevents the use of microbial intervention strategies to improve calf gut health. This study profiled the taxonomic and functional composition of the small intestinal luminal microbiome of neonatal calves using whole-genome sequencing of the metagenome, aiming to understand the dynamics of microbial establishment during early life. Despite highly individualized microbial communities, we identified two distinct taxonomy-based clusters from the collective luminal microbiomes comprising a high level of either Lactobacillus or Bacteroides. Among the clustered microbiomes, Lactobacillus-dominant ileal microbiomes had significantly lower abundances of Bacteroides, Prevotella, Roseburia, Ruminococcus, and Veillonella compared to the Bacteroides-dominated ileal microbiomes. In addition, the upregulated ileal genes of the Lactobacillus-dominant calves were related to leukocyte and lymphocyte chemotaxis, the cytokine/chemokine-mediated signaling pathway, and inflammatory responses, while the upregulated ileal genes of the Bacteroides-dominant calves were related to cell adhesion, response to stimulus, cell communication and regulation of mitogen-activated protein kinase cascades. The functional profiles of the luminal microbiomes also revealed two distinct clusters consisting of functions related to either high protein metabolism or sulfur metabolism. A lower abundance of Bifidobacterium and a higher abundance of sulfur-reducing bacteria (SRB) were observed in the sulfur metabolism-dominant cluster (0.2% ± 0.1%) compared to the protein metabolism-dominant cluster (12.6% ± 5.7%), suggesting an antagonistic relationship between SRB and Bifidobacterium, which both compete for cysteine. These distinct taxonomic and functional clusters may provide a framework to further analyze interactions between the intestinal microbiome and the immune function and health of neonatal calves.

IMPORTANCE Dietary interventions to manipulate neonatal gut microbiota have been proposed to generate long-term impacts on hosts. Currently, our understanding of the early gut microbiome of neonatal calves is limited to 16S rRNA gene amplicon based microbial profiling, which is a barrier to developing dietary interventions to improve calf gut health. The use of a metagenome sequencing-based approach in the present study revealed high individual animal variation in taxonomic and functional abundance of intestinal microbiome and potential impacts of early microbiome on mucosal immune responses during the preweaning period. During this developmental period, age- and diet-related changes in microbial diversity, richness, density, and the abundance of taxa and functions were observed. A correlation-based approach to further explore the individual animal variation revealed potential enterotypes that can be linked to calf gut health, which may pave the way to developing strategies to manipulate the microbiome and improve calf health.

INTRODUCTION

The gut microbiome consists of millions of nonredundant genes that are necessary for microbial survival in the gastrointestinal tract environment, and the majority of these microbial genes (∼99%) are of bacterial origin (1). In addition to the housekeeping processes universal to microbial survival, various biosynthetic processes are also carried out by the microbial genome (2). Recently, the biosynthetic processes of gut microbiota have been used to explain the mechanisms underlying host-microbial interactions in humans (3). For example, the Firmicutes/Bacteroidetes ratio is high in obese individuals, and a high-energy harvesting capacity of Firmicutes relative to Bacteroidetes has been implicated as a causal factor contributing to the obese phenotype (3). Thus, a thorough understanding of microbial functions and taxonomic composition is vital when elucidating the mechanisms underlying host-microbial interactions.

Metagenome (genetic content of a microbial community) sequencing has been widely used to explore human microbiomes (4). Whole-genome sequencing of metagenomes provides information on the microbial gene composition and abundance that can be used to predict microbial functions as well as to profile taxonomic composition (5). Metagenomic analysis of the human infant gut microbiome has advanced our knowledge of the establishment of the early life gut microbiome (6). However, knowledge of the pre-ruminant intestinal microbiome is limited. Although studies have profiled the preweaned calf gut microbiota using 16S rRNA gene amplicon-based sequencing (7–9), only a few publications have used the whole-genome sequencing of metagenomes to report the preruminant rumen microbiome (10–12). Neonatal dairy calves are highly susceptible to infections by a variety of enteric pathogens that target the mucosal epithelium of the small intestine (13). The lack of information on the microbial functions in the intestinal tracts of neonatal calves, which are colonized by a diverse group of pioneer species (8), prevents the use of microbial intervention strategies to improve neonatal calf gut health. The small intestines of neonatal dairy calves harbor region-specific microbiota (8), and a recent study in mice demonstrated region-specific regulation of the intestinal epithelium by the gut microbiota (14). Linkages between early gut microbial composition and dairy calf health and growth have been reported in a study characterizing the fecal microbiota of preweaned calves (7). A follow-up study reported that oral administration of Faecalibacterium prausnitzii, a bacterium that is associated with calf health and growth (7), lowered the incidence of severe diarrhea and diarrhea-associated calf death compared to control calves (15). These findings together suggest that the early microbiome plays a crucial role in neonatal calf health and that early microbial intervention strategies can be effectively used to improve calf health. However, knowledge of the mechanisms underlying host responses to the early microbiome of preweaned calves is limited. Therefore, it is necessary to characterize the taxonomic and functional composition of the early intestinal microbiome and to elucidate the mechanisms by which the microbiome may alter mucosal immune and barrier functions to improve calf gut health.

In this study, whole-genome sequencing of the luminal metagenome and 16S rRNA gene amplicon sequencing of the mucosal-tissue-attached bacteria were used to characterize the taxonomic and functional composition of the small intestinal microbiomes of neonatal calves. Immune-related genes expressed in the neonatal calf ileum were then extracted from the whole transcriptomes of a previously published work conducted using the same animals (16) to explore whether the variations identified in the microbiome may also be reflected in the host mucosal immune system. The results from these analyses provide the first characterization of functions associated with the intestinal microbiome and provide a framework for investigating the functional consequences of early gut microbiota.

RESULTS

Postnatal changes in the taxonomic composition of the small intestinal microbiota.(i) Luminal microbiota. After quality filtering to remove host DNA and artificial reads from the raw sequences, we generated a 1.82 Gb data set (40,555,690 ± 4,656,914 reads, expressed as the mean ± the standard error of the mean [SEM]) from the small intestinal luminal microbiomes (Data Set S1), which were used to assign taxonomy and functions using the M5NR database and the SEED subsystems within the MG-RAST platform, respectively. Small intestinal luminal microbiomes of neonatal calves consisted of bacteria, archaea, fungi, protozoa, and viruses from the first week of life (Fig. S1; Data Set S2, sheet 1). Further analysis of bacteria, the dominant microbial group, revealed 21 phyla, 46 families, and 510 genera in total from all luminal microbiomes. Four bacterial phyla accounted for nearly 93% of all identified bacteria, with substantial variation in the abundance of phyla among individual microbiomes (Firmicutes, 11 to 80%; Bacteroidetes, 0.5 to 75%; Proteobacteria, 0 to 35%; and Actinobacteria, 1 to 85%) (Fig. 1A). These four main phyla were observed in all intestinal regions of all calves (n = 15), with the exception of Proteobacteria, which was absent from the proximal jejunum of one 6-week-old calf (6W). Individual animal variation in the bacterial composition was highest within the lower taxonomic hierarchy, and no single bacterial family or genus was detected in all of the intestinal samples from all animals.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Postnatal changes in the small intestinal microbiomes of preweaned dairy calves. (A) Individual calf variation of the relative abundance of four main bacterial phyla in the digesta-associated communities collected from three different regions of the small intestine at three different ages. Each data point represents an individual calf, and the bar represents the mean relative abundance of a small intestinal region or an age group. The upper panel presents all three age groups within a small intestinal region (15 data points/phylum, 3 charts for 3 intestinal regions), and the lower panel presents all three small intestinal regions within an age category (15 data points/phylum, 3 charts for 3 age groups). (B) The abundance of the small intestinal subsystems (level 1 microbial functions of the SEED hierarchy, only the abundant subsystems [mean relative abundance >1%] are presented). Numerical values represent individual calf IDs. 1W, 1-week-old calves; 3W, 3-week-old calves; 6W, 6-week-old calves. Each data series represents relative abundance of detected subsystems.

Among the bacterial genera detected in the neonatal calf intestine, Clostridium and Eubacterium were present in all jejunal microbiomes, whereas Prevotella and Bacteroides were observed in all ileal microbiomes. Bifidobacterium (5.7% ± 2.2%), Prevotella (7.5% ± 1.8%), Lactobacillus (15.4% ± 3.4%), Clostridium (5.6% ± 0.1%), Bacteroides (8.2% ± 2.6%), Streptococcus (4.6% ± 1.2%), and Eubacterium (2.6% ± 0.4%) were present in all sampled calves, although they were not detected in all intestinal regions surveyed (Fig. S2). For example, Eubacterium was detected in the ileum of all but one 1W calf. Significant temporal variations were evident when comparing the relative abundances of the detected taxa among different age groups within a small intestinal region (Table 1 ). For example, the relative abundance of Ruminococcus was significantly higher (P < 0.05) in older calves (3W and 6W) than 1W calves only in the proximal jejunum (Table 1). Similarly, the relative abundance of Lactobacillus was significantly lower (P < 0.05) in older calves than 1W calves in two jejunal regions but not in the ileum (Table 1). Among the observed significantly different bacterial taxa, many were different when comparing 1W calves with either 3W or 6W calves, but not between 3W and 6W calves (Table 1). There were no regional differences observed in the abundance of bacterial taxa along the calf small intestine either when comparing within an age group or when comparing all age groups together.

View this table:
  • View inline
  • View popup
TABLE 1

Temporal variations in the small intestinal digesta-associated metagenome (taxonomy and function) of the preweaned calves within each small intestinal region

Although temporal variations observed in the small intestinal microbiomes regardless of region were not statistically different, there were increasing (higher relative abundance in older calves) and decreasing (lower relative abundance in older calves) trends in the relative abundance of individual bacterial genera as calves aged (Fig. S3). For example, the mean relative abundances of Bifidobacterium (1W, 2.3% ± 1.0%; 3W, 5.0% ± 2.1%; 6W, 9.8% ± 6.9%) and Ruminococcus (1W, 0.8% ± 0.2%; 3W, 1.7% ± 0.1%; 6W, 2.4% ± 0.5%) increased with calf age (Fig. S3 in the supplemental material), whereas the relative abundances of Bacteroides (1W, 17.3% ± 5.6%; 3W, 3.8% ± 0.6%; 6W, 3.5% ± 0.8%) and Lactobacillus (1W, 21.9% ± 2.4%; 3W, 14.9% ± 10.3%; 6W, 9.4% ± 1.2%) decreased with calf age (Fig. S3).

(ii) Mucosa-attached microbiota.Amplicon sequencing of the V1-V3 region of the 16S rRNA gene revealed that the small intestinal mucosa-attached bacterial community comprised 6 phyla, 38 families, and 104 genera. Similar to the luminal bacteria, the relative abundance of each phylum varied greatly among individual calves (Actinobacteria, 0 to 93.3%; Firmicutes, 0 to 94.4%; Proteobacteria, 0 to 99.9%; and Bacteroidetes, 0 to 35.3%) (Data Set S2, sheet 1). No single bacterial genus was observed in all tissue samples analyzed. The predominantly observed mucosa-attached bacterial genera were Propionibacterium (7.3% ± 4.0%), Prevotella (7.0% ± 2.9%), and Lactobacillus (4.8% ± 2.0%) in 1W samples, Propionibacterium (12.2% ± 8.9%) and Prevotella (6.6% ± 3.3%) in 3W samples, and Bifidobacterium (8.5% ± 4.2%), Sharpea (4.6% ± 3.8%), and Prevotella (3.0% ± 1.4%) in 6W samples (Fig. S2). Prevotella was detected in at least one intestinal region of all calves sampled. When the relative abundances of the mucosa-attached bacterial groups were compared, there were no significant differences observed among either the intestinal region or calf age.

Diversity and richness of the small intestinal microbiota.The diversity (Shannon index) of the luminal microbiota increased with calf age and dietary changes, but this increase was not significantly different when comparing among age groups (Table 2). Although the total number of bacterial genera detected was not different among age groups, there was a significantly (P = 0.02) lower number of predominant genera (relative abundance of >1% in at least one sample) in the ileal lumen of 3W and 6W calves than in the age-matched jejunal communities (Table 2). When the bacterial diversity and richness (number of detected genera) were compared for mucosa-attached communities, they were not significantly different among either intestinal regions or age groups (Table 2).

View this table:
  • View inline
  • View popup
TABLE 2

Small intestinal bacterial diversity and richness of neonatal calvesa

Small intestinal bacterial densities in neonatal calves.Real-time PCR-based estimation (using DNA) revealed that the density of total bacteria and the proportion of Lactobacillus were significantly lower (P < 0.01) in the mucosa-attached communities compared to the luminal communities (Table 3). The individual animal variation was high in the estimated total bacterial density, ranging from 2.63 × 108 to 3.22 × 1011 16S rRNA gene copies/g of sample. Similarly, the density of Lactobacillus varied greatly, between 1.49 × 104 and 3.72 × 109 16S rRNA gene copies/g of sample, among individuals. The proportion of Bifidobacterium was affected by the three-way interaction effect of region by age by sample type (P = 0.02). The proportion of Bifidobacterium was significantly higher in the lumen of the proximal jejunum at 1W compared to those of 3W and 6W (Table 4). The proportion of Bifidobacterium in the proximal jejunum at 1W was also higher in the lumen than the respective mucosa-attached community (Table 4). Similar to total bacteria and Lactobacillus, the density of Bifidobacterium varied greatly among individual animals (1.19 × 107 to 2.21 × 1010 16S rRNA gene copies/g of samples). However, the coefficient of variation (CV) of Bifidobacterium (6.5%) was lower than that of total bacteria (18.5%) and Lactobacillus (18.8%).

View this table:
  • View inline
  • View popup
TABLE 3

Bacterial density in the small intestinal communities of preweaned calves

View this table:
  • View inline
  • View popup
TABLE 4

Bifidobacterium density in the small intestinal communities of preweaned calves

Postnatal changes in the functions of the small intestinal luminal microbiome.Functional assignment using MG-RAST revealed 29 subsystems/level 1 functions (Data Set S2, sheet 2) in the SEED hierarchy for the collective small intestinal microbiomes (all gut regions of all calves, from here onward we refer to them as the collective microbiomes). The functions of the collective microbiomes were dominated by genes involved in the metabolism of carbohydrates and proteins (Fig. 1B). The “respiration” subsystem was observed in all microbiomes, and it contained enzymes involved in fermentation, such as anaerobic glycerol-3-phosphate dehydrogenase (EC 1.1.5.3) and formate hydrogenlyase. Although other subsystems were not observed in all microbiomes, they were detected in more than 70% of the intestinal microbiomes generated in the study. The top subsystems identified from the collective microbiomes of the calf small intestine were “carbohydrate metabolism” (range, 0 to 13.7%; mean ± SEM = 10.5% ± 0.6%), “protein metabolism” (range, 0 to 19.5%; mean ± SEM = 9.7% ± 0.6%), “amino acids and derivatives” (range, 0 to 12.6%; mean ± SEM = 6.7% ± 0.4%), “phages, prophages, transposable elements, plasmids” (range, 0 to 95.4%; mean ± SEM = 9.3% ± 3.0%), “DNA metabolism” (range, 0 to 9.6%; mean = 4.8% ± 0.3%), “RNA metabolism” (range, 0 to 13.7%; mean ± SEM = 4.8% ± 0.3%), and “cofactors, vitamins, prosthetic groups, pigments” (range, 0 to 8.1%; mean ± SEM = 4.3% ± 0.3%).

Temporal changes in the relative abundances of the subsystems were observed when pairwise comparisons (1W versus 3W, 1W versus 6W, and 3W versus 6W) were performed within each gut region (Table 1). For example, the abundance of “protein metabolism” was higher in the older calves (3W, 6W) than 1W calves in the two jejunal regions but not in the ileum (Table 1). Similarly, a lower abundance of “phages, prophages, transposable elements, plasmids” in the older calves (compared to 1W calves) was also observed only in the two jejunal regions (Table 1). Once again, the temporal changes in the relative abundances of these subsystems were mainly observed when comparing 1W calves against older (3W and 6W) calves.

When the metagenomes generated from all three small intestinal regions and metagenomes generated from all three age groups were compared, there were no significant temporal or regional changes, respectively. However, the subsystems could be categorized into three patterns depending on the observed temporal trends of the relative abundance (Fig. S4). The abundance of essential functions for all microbiota, such as “protein metabolism,” “carbohydrate metabolism,” and “DNA and RNA metabolism,” as well as biosynthetic processes (“cofactors, vitamins, prosthetic groups, pigments”), increased from the first to the third week and then remained stable from the third week to the sixth week (Fig. S4). In contrast, the abundance of “phages, prophages, transposable elements, plasmids,” “membrane transport,” “virulence, disease, and defense,” “iron acquisition,” “respiration,” and “sulfur and potassium metabolism” decreased from the first to the third week and then remained stable from the third week to the sixth week (Fig. S4). However, the abundances of “motility and chemotaxis,” “cell wall and capsule,” “regulation and cell signaling,” “secondary metabolism,” “phosphorus metabolism,” and “stress response” remained stable throughout the first 6 weeks of life (Fig. S4).

In total, 185 level 2 functions in the SEED hierarchy were detected within the small intestinal microbiomes (proximal jejunum, 179; distal jejunum, 144; ileum, 183), with no single level 2 function shared among all microbiomes. We identified 16 and 5 level 2 functions present in all distal jejunum and ileum microbiomes, respectively, but no level 2 functions were conserved among all of the proximal jejunum microbiomes (Data Set S2, sheet 2). Microbial functions related to “protein biosynthesis,” “folate and pterines,” “DNA repair,” “lysine, threonine, methionine, and cysteine,” and “RNA processing and modification” were conserved in the distal jejunum and ileum. Temporal changes in the abundance of level 2 functions were evident through pairwise comparisons conducted within each small intestinal region (Table 1). When the collective metagenomes were compared, there were no significant temporal or regional differences observed in the abundance of level 2 functions. The abundances of “DNA repair,” “di- and oligosaccharides,” “carbohydrates,” “transcription,” “alanine, serine, and glycine,” and “pyridoxine (vitamin B6)” increased numerically with age (Fig. S4). In contrast, the abundances of “electron-donating reactions,” “protein biosynthesis,” “cell division and cell cycle,” “central carbohydrate metabolism,” “lysine, threonine, methionine, and cysteine,” and “ABC transporters” decreased numerically with age (Fig. S4).

Individual variation in the taxonomic and functional composition of the luminal microbiome.Principal-coordinate analysis (weighted and unweighted UniFrac distance matrices built using operational taxonomic units clustered at 97% similarity) resulted in no segregation of luminal microbiomes based on either small intestinal region or calf age (Fig. S5). We compared coefficients of variance (CV) within an individual animal to those of the relevant age group to understand the individual variation of luminal microbiomes (Fig. S6; Data Set S2, sheet 4). A lower intraindividual variation (variation among three intestinal regions) than the interindividual variation (variation among individuals within an age group) was evident in microbiomes of all age groups at the functional and taxonomic levels. Therefore, the taxonomic and functional compositions of the collective luminal microbiomes (45 metagenomes generated from all luminal samples) were further analyzed using a hierarchical clustering approach.

Use of the 27 most frequently detected bacterial genera (in at least half of the lumen microbiomes or ≥22 of 45 microbiomes) revealed that the majority of luminal microbiomes (24 of 45 microbiomes) were not correlated with one other based on the taxonomic composition. However, we observed two distinct clusters (Fig. 2A) among the highly correlated (R2 ≥ 0.9, P < 0.01) luminal microbiomes. One cluster consisted of microbiomes containing a higher relative abundance of Lactobacillus (ileal profiles of four calves: 1W, calf5 [72.4%]; 3W, calf7 [86.4%] and calf8 [85.8%]; 6W, calf12 [56.4%]; jejunal profiles of six calves: proximal jejunum: 1W, calf2 [40.4%], calf3 [23.0%], and calf5 [36.2%]; 6W, calf12 [28.0%]; distal jejunum: 1W, calf1 [38.7%], calf3 [14.9%], and calf5 [58.1%]; 3W, calf8 [14.9%]; 6W, calf12 [31.2%]) compared to all other microbiomes. Another cluster consisted of microbiomes containing a higher relative abundance of Bacteroides compared to all other microbiomes (0.01 to 5.5%; all three small intestinal regions of one 1W calf [calf4: proximal jejunum, 23.7%; distal jejunum, 41.7%; ileum, 56.3%], the distal jejunum of one 1W calf [calf2, 22.2%] and one 6W calf [calf11, 16.4%], and the ileum of two 1W calves [calf2, 45.8%; calf3, 33.9%], as well as one 3W calf [calf6, 17.9%]). Comparison of the relative abundance of bacterial genera between two clusters revealed that Bacteroides, Prevotella, Roseburia, Ruminococcus, and Veillonella were significantly lower in the Lactobacillus-dominant ileal microbiomes compared to the Bacteroides-dominant ileal microbiomes (Table S1).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Clustering of the collective microbial profiles based on taxonomic/functional similarities. (A) Clustering of calves based on microbial taxonomic composition. (B) Clustering of calves based on microbial functional abundance. Numbers on the x axis represent calf IDs. PJ, proximal jejunum; DJ, distal jejunum; IL, ileum. Abundances of bacterial genera between two function-based clusters are compared using Metastats, and data are presented as means ± the SEM. Clustering is based on the Spearman rank correlation coefficient between two samples generated either based on the relative abundance of 27 bacterial genera identified in at least 50% of samples or based on the relative abundance of subsystems. Each row and column represent an individual microbial profile generated from small intestinal samples (45 microbial profiles).

Hierarchical clustering of the collective metagenome using functional composition revealed a high correlation (R2 ≥ 0.9, P < 0.01) among the majority of luminal microbiomes (34 of 45 microbiomes), and we identified two distinct clusters among these highly correlated microbiomes (Fig. 2B). One cluster contained samples from 3W and 6W calves (cluster 1), whereas the other cluster contained samples from 1W and 3W calves (cluster 2). Comparison of these two clusters revealed that functions related to “protein metabolism” (cluster 1, 12.3% ± 0.5%; cluster 2, 9.0% ± 0.5%; P < 0.01), “amino acids and derivatives” (cluster 1, 7.9% ± 0.4%; cluster 2, 6.4% ± 0.6%; P = 0.06) and “nucleosides and nucleotides” (cluster 1, 5.8% ± 0.3%; cluster 2, 3.8% ± 0.3%; P < 0.01) were significantly enriched in cluster 1 compared to cluster 2. In contrast, functions related to “cell wall and capsule” (cluster 1, 3.1% ± 0.3%; cluster 2, 4.4% ± 0.3%; P = 0.01) and “sulfur metabolism” (cluster 1, 0.3% ± 0.05%; cluster 2, 0.8% ± 0.2%; P = 0.02) were enriched in cluster 2 compared to cluster 1. Comparison of the taxonomic composition of the microbiomes belonging to the two functional clusters revealed that the relative abundances of the predominant genera Bifidobacterium, Bacillus, Streptococcus, Lactococcus, and Corynebacterium were significantly higher in cluster 1 than cluster 2, while the abundance of sulfur-reducing bacteria (Desulfarculus, Dethiosulfovibrio, Desulfatibacillum, Desulfuromonas, Desulfurispirillum, and Desulfotalea) was significantly higher in cluster 2 than cluster 1 (Fig. 2B; Table S2). Furthermore, there were 19 and 16 additional genera that were relatively enriched in cluster 1 and cluster 2, respectively (Table S2).

Linking identified ileal bacterial taxonomy-based clusters to the host transcriptome.To explore whether the variations observed in the ileal microbiomes could also be reflected in the host at the transcriptome level, ∼3,000 transcripts involved in mucosal immune-related functions were extracted from a previously published ileal transcriptome data set generated from the same calves (16) (Data Set S3). Among these immune-related transcripts, 28 were upregulated (>2-fold change) in the calves with a higher abundance of Lactobacillus in the ileum than in those with a higher abundance of Bacteroides in the ileum (Data Set S2, sheet 5). Among these upregulated (enriched) transcripts, six were significantly different (P < 0.05) between the two groups, and five tended to be different (0.05 < P < 0.1). Similarly, 63 immune-related transcripts were upregulated in the calves with a high abundance of Bacteroides in the ileum compared to levels in calves with a high abundance of Lactobacillus in the ileum (Data Set S3, sheet 5). Among these enriched transcripts, 32 were significantly different (P < 0.05) between the two groups, while 15 tended to be different (0.05 < P < 0.1). GO enrichment of the transcripts upregulated in the ileal Lactobacillus-dominant calves revealed that they were involved in “leukocyte and lymphocyte chemotaxis,” “cytokine/chemokine-mediated signaling pathway,” and “inflammatory responses” (Fig. 3A; Data Set S2, sheet 6). Protein-coding genes CXCL9 (P = 0.12), CXCL10 (P < 0.05), CXCL11 (P < 0.05), CCL2 (P = 0.18), and CCL14 (P < 0.1) were mainly annotated to these enriched functions (Fig. 3B). In contrast to the ileal Lactobacillus-dominant calves, the ileal Bacteroides-dominant calves had functions enriched for processes related to “cell adhesion,” “response to stimulus,” “cell communication,” and “regulation of MAPK cascades” (Fig. 3A; Data Set S2, sheet 7). The main annotated protein coding genes for these enriched functions included PLCE1 (P < 0.05), CCL22 (P < 0.05), EGFR (P < 0.05), IL1RL1 (P = 0.05), IL1R2 (P = 0.10), IL17RD (P = 0.53), and ITGB6 (P = 0.15) (Fig. 3B).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Ileal mucosal immune-related genes and functions of calves belong to two taxonomic clusters. (A) Functions of the genes upregulated/enriched in taxonomic clusters. (Upper panel) Functions of the genes enriched in the ileal tissue of the Lactobacillus-dominant calves; (lower panel) functions of the genes enriched in the ileal tissue of the Bacteroides-dominant calves. Functional enrichment was performed using GO enrichment analysis in Gene Ontology Consortium. The P value indicates significance for the enrichment in the data set of the listed GO identifier. A P value close to zero, that is, a higher –log(P), means that the group of genes associated with the particular GO term is more significant and it is less likely the observed annotation of the particular GO term to a group of genes occurs by chance. (B) Main genes annotated to enriched immune functions. Each bar represents the fold change expression of annotated genes. The fold change is calculated by dividing the expression of genes (counts per million [cpm]) in Lactobacillus-dominant calves by the expression of genes (mean cpm) in Bacteroides-dominant calves for the Lactobacillus-dominant cluster and vice versa for the Bacteroides-dominant cluster. CXCL9, chemokine ligand 9; CXCL10, chemokine ligand 10; CXCL11, chemokine ligand 11; PLCE1, phospholipase C epsilon 1; CCL22, C-C motif chemokine ligand 22; EGFR, epidermal growth factor receptor.

DISCUSSION

Although functions of the human gut microbiome have been studied extensively through metagenomics approaches, such knowledge is very limited for ruminants, especially preruminants. The present study profiled the small intestinal microbiomes of healthy neonatal calves to gain a better understanding of the taxonomic and functional composition and temporal dynamics associated with the establishment of the gut microbiome during early life. Calf small intestinal microbiomes were mainly colonized by bacteria, although metagenome sequencing of the luminal content revealed the presence of archaea, fungi, protozoa, and viruses from the first week of life. The genus Methanobrevibacter dominated the small intestine archaeal population, whereas the family Microviridae dominated the viral (phage) population of the neonatal calf intestine (Fig. S1B). These observations are in agreement with human studies reporting a higher abundance of Methanobrevibacter (17) and Microviridae (18) in fecal microbiomes and the presence of Methanobrevibacter in preweaned dairy calves (19, 20). The presence of the fungal phyla Ascomycota and Basidiomycota (Fig. S1B), which were the most abundant fungi identified in the present study, was also reported in a metagenomics-based study conducted using fecal samples of healthy humans (21). The main protozoa observed in the neonatal calf small intestine belonged to the genus Cryptosporidium (Fig. S1B). C. parvum is the main parasitic protozoan that contributes to neonatal calf diarrhea; however, calves may harbor Cryptosporidium even in the absence of clinical signs or symptoms (22). This is the first study to present the composition of all microbial groups in the intestine of neonatal calves.

The increased diversity and richness in the small intestinal bacterial communities (both luminal and tissue) with increasing age suggest a progressive establishment of a complex microbiome during the first 6 weeks of life. Furthermore, total bacterial densities (qPCR-based quantification) of the small intestinal tissue and lumen environments of neonatal calves reached levels similar to those of weaned calves (23) from the first week of life. The calves used in this study were fed whole milk during the first week postpartum and whole milk with ad libitum calf starter supplementation thereafter. Thus, the observed changes in the luminal microbiomes were both age and diet driven. The dietary changes may have contributed substantially to the significant differences in microbial composition observed when comparing 1W versus 3W and 6W calves. The higher similarity between 3W and 6W calves further suggests that the luminal microbiome may have started to stabilize once solid feed was introduced. Regardless, a high individual animal variation in the composition and density of bacteria was evident throughout the neonatal period. In humans, the fecal microbiome begins to stabilize and increase in similarity among individuals around 3 years of age with the intake of a stable solid diet (24). Similarly, high individual animal variation in the rumen microbiota has been reported to decrease with increasing age in cattle (25). Thus, high individual variation in the intestinal microbiome may be a conserved trait in neonates of different species. Host genetics has also been shown to shape gut microbiomes, leading to an individualized gut microbiome in mice (26) and heritable bacterial taxa in humans (27–29). The impact of host genetics may be one of the factors contributing to the high individual variation observed in the gut microbiomes of neonates, whose diet is not yet stabilized.

Principal-coordinate analysis revealed no segregation of luminal microbial profiles in the present study. However, the correlation-based hierarchical clustering approach we used grouped the collective luminal microbiomes based on their taxonomic and functional compositions. Metagenome sequencing of the fecal microbiome has also been used to categorize the human population broadly into three enterotypes based on the individual microbial profile (30). Despite the small number of samples, our approach also resulted in two taxonomic clusters of microbiomes with either high Lactobacillus or high Bacteroides abundances in the ilea of the neonatal calves, suggesting potential enterotyping of neonatal calves. Comparisons of the ileal transcriptomes of the calves belonging to these two clusters further revealed variations in the expression of mucosal immune-related genes and enriched immune functions between them. These findings together suggest a link between potential enterotyping of neonatal calves and their mucosal immune status during early life.

The Lactobacillus-dominant calves in the ileal lumen exhibited a higher expression of CXCL9, CXCL10, and CXCL11 in the ileal tissue than did the Bacteroides-dominant calves in the ileal lumen. These proinflammatory chemokines have been shown to be upregulated in dendritic cells and human intestinal epithelial cells after exposure to Lactobacillus acidophilus (31, 32). However, such upregulation quickly returned to baseline levels one hour after exposure to L. acidophilus (32), suggesting a potential priming of immune responses by Lactobacillus during the first encounter. The priming effect of early Lactobacillus exposure may prepare the host for subsequent anti-inflammatory regulatory responses, which may be a crucial process in neonatal calves with a naive mucosal immune system. In contrast to these studies, infectious models (colitis and pneumovirus infection in mice) have demonstrated that Lactobacillus exerts an anti-inflammatory effect through the downregulation of CXCL10 and CXCL11 (33, 34). Therefore, these studies suggest that the beneficial effects of Lactobacillus may vary depending on the context, with a priming of immune responses upon initial encounter but a reduction of inflammatory responses during infections.

The increased expression of PLCE1, CCL22, and EGFR in the ileal tissue of Bacteroides-dominant calves suggests an increased capacity for T-cell-dependent responses in the guts of these calves, since these mitogen-activated protein kinase cascade-related genes are involved in the regulation of T-cell differentiation (35). Our transcriptome data from calves with a higher abundance of either Lactobacillus or Bacteroides in the ileal lumen revealed differential expression of genes that may be involved in the activation of mucosal immune responses. The identified enriched mucosal immune functions associated with taxonomy-based microbiome clusters may indicate that the calves with a “beneficial-bacterium-enriched microbiome” are resilient to enteric health issues such as diarrhea during early life. Thus, validation of the enriched genes in a larger population of healthy animals and/or use of enteric challenged models with phenotypic measurements will be required to determine whether Lactobacillus and Bacteroides represent enterotypes in the calf small intestine that are associated with specific immune functions and health outcomes.

Although taxonomic profiles of the luminal microbiomes were markedly different among individual calves, such differences were not directly reflected in functional profiles. This may be due to the observed higher abundance of universal microbial functions that are important for microbial survival, such as “protein biosynthesis,” “RNA processing and modification,” “DNA replication,” “central carbohydrate metabolism,” and “cell cycle.” Nevertheless, two distinct clusters based on functional abundance were observed from the collective luminal microbiomes of neonatal calves. Functional cluster 1 was enriched for functions related to the microbial protein metabolism, which has been suggested to play an important role in host health through the regulation of amino acid homeostasis in the gut (36).

The luminal microbiomes of functional cluster 1 had a higher abundance (69-fold) of Bifidobacterium than those in functional cluster 2, while those of functional cluster 2 had a higher abundance of sulfur-reducing bacteria (SRB) than functional cluster 1. A recent study has shown that Bifidobacterium utilizes milk protein effectively and requires cysteine, a sulfide amino acid, for growth (37). SRB also utilize cysteine as an energy source (38). Thus, the observed low abundance or absence of Bifidobacterium in the luminal microbiomes of functional cluster 2 calves may be due to a lack of cysteine for their growth as a result of competition with SRB. Dissimilatory sulfur reduction and cysteine degradation by SRB produce H2S, which can have detrimental effects on intestinal epithelial cells (38). An accumulation of H2S in the colon inhibits the activity of acyl coenzyme A dehydrogenase, which catalyzes the oxidation of butyrate in colonocytes and compromises epithelial barrier permeability (39). Bifidobacterium has been reported to promote intestinal epithelial barrier functions (40, 41), while SRB have been reported to have a negative impact on the epithelial barrier (39). Although these observed differences in bacterial composition might suggest potential differences in epithelial barrier functions between functional clusters 1 and 2, transcriptome data (16) revealed no difference in the expression of tight-junction genes (<2-fold change in the expression of genes belonging to the claudin family, occludin, zonula occludens, and junctional adhesion molecules, with P > 0.1) when comparing animals in these two clusters. Intestinal barrier integrity loss can be measured using changes in the expression of tight-junction proteins, epithelial cell damage, or the translocation of bacterial products (42). Therefore, assessing barrier function by using these alternative indicators may provide a more sensitive measurement of how compositional and functional variations in the gut microbiome impact epithelial barrier function.

A high rate of diarrhea-related neonatal deaths remains a major health problem in the cattle industry (43). Neonatal calf diarrhea is often treated with antibiotics, and often prophylactic antibiotics are given to calves during the first month to prevent colonization of the small intestine by pathogenic Escherichia coli (one of the diarrhea-causing bacteria) (44). With the increasing restrictions on the prophylactic use of antibiotics (44), it is important to explore other microbial manipulation techniques that can limit or prevent pathogen colonization of the calf small intestine. Feeding probiotics has been shown to reduce the frequency of pathogenic E. coli in the calf gastrointestinal tract (45), suggesting that promoting the colonization of beneficial bacteria may exclude pathogens. Limited knowledge of the calf microbiome, especially in the small intestine, is a barrier to developing such strategies to prevent or treat diarrhea. Thus, our findings on gut microbial markers linked with mucosal immune functions during early life may pave the way to developing strategies to manipulate the microbiome and improve calf health.

The taxonomic composition of the small intestinal bacteria observed using metagenome sequencing revealed that the luminal microbiome in the present study was similar to previous studies reporting a higher abundance of Firmicutes in the small intestine (8) and the fecal microbiome (7) of young calves using 16S rRNA amplicon sequencing. However, there were significant differences within the lower taxonomic hierarchy that could be attributable to differences in the microbial profiling tools (whole-genome sequencing versus 16S amplicon sequencing), calf age, and sample type (fecal versus small intestine). Sequencing the V1-V3 region of 16S rRNA genes revealed a higher abundance of Lactobacillus in the ileum and Sharpea in the proximal jejunum at 3 weeks of life (8). The use of metagenomic profiling of luminal content in the present study, however, revealed a higher abundance of Prevotella (17.9% ± 12.8%) in the ileum and Streptococcus (12.0% ± 6.1%) in the proximal jejunum of 3-week-old calves. In contrast to luminal communities, the mucosa-attached communities in the present study differed from those reported in a previous study, even at higher levels of the taxonomic hierarchy (8). In contrast to our previous study (8), the present study revealed a higher abundance of Actinobacteria (41.2% ± 15.1%) in the ileal tissues, followed by Firmicutes (26.8% ± 9.6%), Proteobacteria (20.1% ± 15.5%), and Bacteroidetes (11.8% ± 6.1%), and a higher abundance of Proteobacteria (66.5% ± 23.2%) in the jejunal tissues, followed by Firmicutes (14.0% ± 8.1%), Actinobacteria (11.0% ± 10.5%), and Bacteroidetes (8.0% ± 4.2%), at 3 weeks of age. The use of animals from different herds and feeding calves different starter rations may explain some of the observed differences between these two studies. In addition, the small intestinal tissue-attached bacterial community is closely scrutinized by the host mucosal immune system (46), which may also contribute to taxonomic variation in the microbial communities among different animals.

Conclusion.This study generated a detailed analysis of the taxonomic and functional development of the small intestinal microbiome in healthy preweaned calves using metagenome sequencing of lumen content. During this developmental period, age- and diet-related increases in microbially diverse populations, richness, and density were observed; however, the taxonomic composition (both tissue and lumen) varied widely among individual animals. The taxonomic comparison of luminal microbiomes classified calves as either Lactobacillus dominant or Bacteroides dominant in the ileum, while the functional comparison of luminal microbiomes revealed an antagonism between SRB and Bifidobacterium. These distinct taxonomic and functional clusters may provide a framework to further analyze interactions between the intestinal microbiome and the immune function, health, and growth of neonatal calves to improve calf gut health.

MATERIALS AND METHODS

Animal experiments and sampling.All experimental protocols were approved by the Livestock Care Committee of the University of Alberta (AUP00001012) and were conducted following the guidelines of the Canadian Council on Animal Care. This study used Holstein bull calves (n = 18) reared at the Dairy Research and Technology Center (DRTC), University of Alberta (Edmonton, Alberta, Canada), under the DRTC standard management practices. Calves used in the present study received 4 liters of colostrum/day during the first 3 days postpartum and 4 liters of whole milk/day from the fourth day up to the end of the trial. Calves had ad libitum access to calf starter (canola meal, soybean meal, wheat, barley, corn, peas, and vitamins and minerals; 29.5% crude protein; 27.1% starch [Wetaskiwin Co-Op Country Junction, Wetaskiwin, Alberta, Canada]) from the second to the sixth week of life. Calves used in the study had no records of respiratory or enteric diseases or use of treatments for any disease. Calves were humanely euthanized using a captive bolt gun stunning method and sampling was completed within 30 min after euthanization. Mucosal tissue and lumen (digesta) samples (proximal jejunum, distal jejunum and ileum) were collected from calves at week 1 (1W; n = 6), week 3 (3W; n = 6), and week 6 (6W; n = 6); snap-frozen; and stored at −80°C. First, digesta was collected into 50-ml Falcon tubes, and then mucosal tissue was rinsed three times with sterile phosphate-buffered saline (pH 7.2) before transfer into tubes. Sampling sites were maintained consistently among all animals. Ileal samples were collected 30 cm proximal to the ileocecal junction, distal jejunal samples were collected 30 cm proximal to the collateral branch of the cranial mesenteric artery, and proximal jejunal samples were collected 100 cm distal to the pylorus sphincter. Ten-centimeter intestinal segments were collected from each site in the middle of each segment aligned with the above-mentioned measurements.

Analysis of the small intestinal microbiomes.Total DNA was extracted from tissue and digesta samples using the repeated bead-beating plus column method (47). Briefly, samples were first ground using mortar and pestle with liquid nitrogen and ∼0.5 g of samples was measured into a 2-ml screw-cap tube (BioSpec, Bartlesville, OK) containing 0.5 g of 0.5-mm zirconia beads (BioSpec). The samples were then subjected to bead beating using a Mini-Beadbeater 8 (BioSpec) at 5,000 rpm for 3 min in the presence of a sodium dodecyl sulfate (SDS)-based lysis buffer. Impurities and SDS used during bead-beating steps were removed by ammonium acetate and nucleic acids were precipitated using isopropanol. Finally, recovered nucleic acids were purified using QIAamp Fast DNA stool minikit (Qiagen, Frederick, MD) according to the manufacturer’s instructions. DNA concentration was measured using Qubit 2.0 (Thermo Fisher Scientific, Massachusetts) and used in library preparation.

DNA libraries were prepared for metagenome sequencing (n = 5/age group, 45 small intestinal luminal samples; one animal was removed during the library preparation due to the low DNA amount) using the TruSeq DNA PCR-free library preparation kit (Illumina, CA). The genomic DNA was first normalized with a resuspension buffer to a final volume of 55 µl at 20 ng/µl, and 50 µl of the solution was transferred into a Covaris microTUBE (Covaris, Inc., MA) for fragmentation using a Covaris S2 focused-ultrasonicator. Then, the cleaned up fragmented DNA was subjected to end repair and size selection (350 bp), followed by the adenylation of the 3′ ends and ligation of the adaptor index. Each metagenomic library was then quantified using a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, MA) and whole-genome sequencing of 100-bp paired-end reads was performed using HiSeq 2000 (Illumina, CA) and TrueSeq cluster kit V3 and TrueSeq SBS kit V3 sequencing kits at Génome Québec (Montréal, Canada).

Raw sequences (6.91 Gb) of metagenome sequencing were first demultiplexed (CASAVA version 1.8; Illumina) and uploaded into the MG-RAST metagenomic analysis server (48), version 3.3.9, and paired ends were joined before submitting for processing. Artificial replicates, host (bovine) DNA, and low quality (Phred score < 25) sequences were removed from the raw data, and the remaining good quality sequences (1.82 Gb; 26% of the generated raw data) were used to assign taxonomy and functions. The subsystems approach in the SEED hierarchy was used to assign microbial functions of the small intestinal microbiomes, and the taxonomic assignments were made using M5 nonredundant protein database (M5NR) within MG-RAST. A maximum cutoff E value of 1e–10, maximum identity of 70% and a maximum alignment length of 50 were used for data selection criteria for both taxonomy and function abundance analyses. Taxonomic and functional abundances of luminal microbial communities obtained from the MG-RAST platform were analyzed using metastats (49) to explore temporal (1 week versus 3 weeks, 1 week versus 6 weeks, and 3 weeks versus 6 weeks) and regional (proximal jejunum [PJ] versus distal jejunum [DJ]; PJ versus ileum [IL]; and DJ versus IL) changes in the abundance of the identified microbial taxa and functions.

DNA extracted from tissue samples was used in 16S rRNA gene-based amplicon sequencing (V1-V3 region) due to the higher abundance of host DNA in the genomic DNA extracted from tissue, which would be sequenced during the whole-genome sequencing. DNA extracted from tissue samples (n = 3/age group, 27 small intestinal tissue samples) was amplified using 27F and 338R primers containing barcodes (8), and purified PCR products (∼400 bp) were subjected to 454 sequencing (Roche GS-FLX Titanium) at Génome Québec. Data were analyzed using the Quantitative Insight into Microbial Ecology (QIIME) toolkit (50), following the removal of low-quality (Phred score < 25, length < 100 bp) and chimeric sequences as described by Malmuthuge et al. (8). Briefly, chimeric sequences were filtered out using usearch61 within the QIIME platform, and the remaining high quality sequences were used to perform operational taxonomic unit clustering at 97% similarity and taxonomic assignment using the Greengenes database (13_8 release).

Estimation of small intestinal bacterial densities using quantitative real-time PCR.Total bacterial, Lactobacillus, and Bifidobacterium densities (n = 6/age group, 54 luminal samples and 54 tissue samples) were estimated using previously published primers (51) and SYBR green chemistry (Fast SYBR Green Master Mix; Applied Biosystems) with a StepOnePlus real-time PCR system (Applied Biosystems, Foster City, CA). The copy number of 16S rRNA gene/g of fresh tissue or luminal content was then calculated using the equation described by Li et al. (52). All bacterial density data were analyzed using the MIXED procedure and repeated measurement experimental design in SAS (SAS 9.4; SAS, Inc., Cary, NC) with small intestinal region as the repeated measurement and animal as the experimental unit. Compound symmetry covariance structure was selected as the best fit by the Bayesian information criteria (BIC). The following statistical model was fitted to test the effect of age, gut region, and sample type on bacterial density: Yijkl = μ + Ai + Rj + Sk + (AR)ij + (AS)ik + (RS)jk + (ARS)ijk + eijkl, where Y is the bacterial density (total bacteria, Lactobacillus, Bifidobacterium), μ is the mean, A is the calf age, R is the small intestinal region, S is the sample type (tissue, content), and e is the residual error.

Clustering analysis of luminal microbiomes.We used a hierarchical clustering approach to obtain the associations among the luminal microbiomes using metagenomes generated from all small intestinal samples (45 samples), regardless of calf age or small intestinal region. The 27 most frequently detected (in at least half of the luminal samples or ≥22 of 45 samples) bacterial genera were used for the taxonomic composition clustering and 29 subsystems (level 1 functions) were used for the functional composition. First, a Spearman rank correlation analysis was performed using the cor() function in R to obtain a correlation (ρ) between each pair of small intestinal microbiomes. Then, a distance matrix containing the square root of ρ values was used to perform hierarchical clustering using the hclust() function within R package (V 3.4.1). The observed taxonomic and functional clusterings were then validated by assessing the clustering tendency using the get_clust_tendency() function in the R packages cluster and factoextra to calculate Hopkins statistics. Approximately unbiased (AU) probability values (P values) of the validated clusters were then calculated using a multiscale bootstrap resampling (1,000 bootstrap replications) in the R package pvclust (53) to identify distinct clusters among the highly correlated luminal microbiomes. The relative abundances of bacterial genera and microbial functions between two taxonomic clusters and two functional clusters were then compared using a nonparametric t test (Mann-Whitney-Wilcoxon) in R (3.3.1v) and a multiple test correction was performed according to the method described by Benjamini and Hochberg in 1995 (54) to identify the differences associated with microbial clusters.

Clustering tendency validation.Validation of clustering tendency revealed that Lactobacillus-dominant and Bacteroides-dominant clusters are significantly different clusters (Hopkins statistic, 0.1898534; Hopkins value close to 0 [far less than 0.5] means that data are significantly clusterable). AU probability values (P values) by multiscale bootstrap resampling was 100% for Lactobacillus-dominant cluster and 98% AU for Bacteroides-dominant cluster, indicating that these two clusters were strongly supported by the taxonomic composition data and that they were not randomly formed clusters. Similarly, the observed lower Hopkins statistic (0.1367969) for function-based clusters also suggested that observed two function-based clusters are significantly different. The AU for the protein metabolism-dominant cluster was 98%, and it was 97% for the sulfur metabolism-dominant cluster, implying that these two significantly different clusters were supported by the functional composition data.

Identifying host transcriptome differences in response to the microbiome.A total of 3,306 mucosal immune-related genes (counts per million >1 in at least 50% of the samples) extracted from a whole ileal transcriptome generated using the same calves (16) (Data Set S3) were used to explore the differences in the mucosal immune systems of the calves belonging to the two taxonomy-based clusters. Calves with either Lactobacillus-dominant ileal microbiomes or Bacteroides-dominant ileal microbiomes were chosen to compare the expression of immunity-related genes using a nonparametric t test (Mann-Whitney-Wilcoxon) in R (3.3.1V). Multiple test correction was performed according to the method described in Benjamini and Hochberg (54). Differentially expressed genes were declared at fold change >2 and a corrected P value of <0.05. All transcripts with >2-fold change between the two groups were then subjected to a GO enrichment analysis (http://www.geneontology.org/page/go-enrichment-analysis) to identify the enriched mucosal immune responses in the calves belonging to the two taxonomic clusters.

Data availability.Sequence data have been deposited in the NCBI Sequence Read Archive (SRA) under accession numbers SRP097207 and SRP050950.

ACKNOWLEDGMENTS

L.L.G. received funding support from Alberta Livestock and Meat Agency, Ltd., Edmonton, Canada (project 2011F129R), and an NSERC discovery grant. N.M. holds a Banting Postdoctoral Fellowship and received funding support from the Alberta Innovates Doctoral Graduate Student Scholarship. G.L. received funding support from the Alberta Innovates Doctoral Graduate Student Scholarship. P.J.G. holds a Tier I CRC funded by the Canadian Institute for Health research (CIHR).

We thank the staff at the Dairy Research Technology Center (University of Alberta), M. Zhou, X. Sun, Y. Chen, E. Hernandez-Sanabria, and S. Urrutia for their assistance with sample collection and X. Sun and Y. Chen for their assistance with metagenomics library preparation and real-time PCR.

FOOTNOTES

    • Received 17 October 2018.
    • Accepted 3 January 2019.
    • Accepted manuscript posted online 18 January 2019.
  • Supplemental material for this article may be found at https://doi.org/10.1128/AEM.02534-18.

  • Copyright © 2019 American Society for Microbiology.

All Rights Reserved.

REFERENCES

  1. 1.↵
    1. Qin J,
    2. Li R,
    3. Raes J,
    4. Arumugam M,
    5. Burgdorf KS,
    6. Manichanh C,
    7. Burgdorf KS,
    8. Manichanh C,
    9. Nielsen T,
    10. Pons N,
    11. Levenez F,
    12. Yamada T,
    13. Mende DR,
    14. Li J,
    15. Xu J,
    16. Li S,
    17. Li D,
    18. Cao J,
    19. Wang B,
    20. Liang H,
    21. Zheng H,
    22. Xie Y,
    23. Tap J,
    24. Lepage P,
    25. Bertalan M,
    26. Batto J-M,
    27. Hansen T,
    28. Paslier DL,
    29. Linneberg A,
    30. Nielsen HB,
    31. Pelletier E,
    32. Renault P,
    33. Sicheritz-Ponten T,
    34. Turner K,
    35. Zhu H,
    36. Yu C,
    37. Li S,
    38. Jian M,
    39. Zhou Y,
    40. Li Y,
    41. Zhang X,
    42. Li S,
    43. Qin N,
    44. Yang H,
    45. Wang J,
    46. Brunak S,
    47. Dore J,
    48. Guarner F,
    49. Kristiansen K,
    50. Pedersen O,
    51. Parkhill J,
    52. Weissenbach J,
    53. MetaHIT Consortium, Bork P,
    54. Ehrlich SD,
    55. Wang J
    . 2010. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464:59–65. doi:10.1038/nature08821.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Shafquat A,
    2. Joice R,
    3. Simmons SL,
    4. Huttenhower C
    . 2014. Functional and phylogenetic assembly of microbial communities in the human microbiome. Trends Microbiol 22:261–266. doi:10.1016/j.tim.2014.01.011.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    1. Turnbaugh PJ,
    2. Hamady M,
    3. Yatsunenko T,
    4. Cantarel BL,
    5. Duncan A,
    6. Ley RE,
    7. Sogin ML,
    8. Jones WJ,
    9. Roe BA,
    10. Affourtit JP,
    11. Egholm M,
    12. Henrissat B,
    13. Heath AC,
    14. Knight R,
    15. Gordon JI
    . 2009. A core gut microbiome in obese and lean twins. Nature 457:480–484. doi:10.1038/nature07540.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    The Human Microbiome Project Consortium. 2012. A framework for human microbiome research. Nature 486:215–221. doi:10.1038/nature11209.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Petrosino JF,
    2. Highlander S,
    3. Luna RA,
    4. Gibbs RA,
    5. Versalovic J
    . 2009. Metagenomic pyrosequencing and microbial identification. Clin Chem 55:856–866. doi:10.1373/clinchem.2008.107565.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Clemente JC,
    2. Ursell LK,
    3. Parfrey LW,
    4. Knight R
    . 2012. The impact of the gut microbiota on human health: an integrative view. Cell 148:1258–1270. doi:10.1016/j.cell.2012.01.035.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    1. Oikonomou G,
    2. Teixeira AGV,
    3. Foditsch C,
    4. Bicalho ML,
    5. Machado VS,
    6. Bicalho RC
    . 2013. Fecal microbial diversity in pre-weaned dairy calves as described by pyrosequencing of metagenomic 16S rDNA: associations of Faecalibacterium species with health and growth. PLoS One 8:e63157. doi:10.1371/journal.pone.0063157.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Malmuthuge N,
    2. Griebel PJ,
    3. Guan LL
    . 2014. Taxonomic identification of commensal bacteria associated with the mucosa and digesta throughout the gastrointestinal tracts of pre-weaned calves. Appl Environ Microbiol 80:2012–2028. doi:10.1128/AEM.03864-13.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Dill-McFarland KA,
    2. Breaker JD,
    3. Suen G
    . 2017. Microbial succession in the gastrointestinal tract of dairy cows from 2 weeks to first lactation. Sci Rep 7:40864. doi:10.1038/srep40864.
    OpenUrlCrossRef
  10. 10.↵
    1. Li RW,
    2. Connor EE,
    3. Li C,
    4. Baldwin Vi RL,
    5. Sparks ME
    . 2012. Characterization of the rumen microbiota of pre-ruminant calves using metagenomic tools. Environ Microbiol 14:129–139. doi:10.1111/j.1462-2920.2011.02543.x.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Li W,
    2. Han Y,
    3. Yuan X,
    4. Wang G,
    5. Wang Z,
    6. Pan Q,
    7. Gao Y,
    8. Qu Y
    . 2017. Metagenomic analysis reveals the influences of milk containing antibiotics on the rumen microbes of calves. Arch Microbiol 199:433–443. doi:10.1007/s00203-016-1311-8.
    OpenUrlCrossRef
  12. 12.↵
    1. Austin K,
    2. Cunningham HC,
    3. Powell SR,
    4. Carpenter KT,
    5. Cammack KM
    . 2018. Investigation of maternal breed and rearing type on the calf rumen microbiome from day 28 through weaning. Transl Anim Sci 2:S125–S129. doi:10.1093/tas/txy034.
    OpenUrlCrossRef
  13. 13.↵
    1. Cho Y,
    2. Yoon K-J
    . 2014. An overview of calf diarrhea: infectious etiology, diagnosis, and intervention. J Vet Sci 15:1–17. doi:10.4142/jvs.2014.15.1.1.
    OpenUrlCrossRef
  14. 14.↵
    1. Sommer F,
    2. Nookaew I,
    3. Sommer N,
    4. Fogelstrand P,
    5. Backhed F
    . 2015. Site-specific programming of the host epithelial transcriptome by the gut microbiota. Genome Biol 16:62. doi:10.1186/s13059-015-0614-4.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Foditsch C,
    2. Pereira RVV,
    3. Ganda EK,
    4. Gomez MS,
    5. Marques EC,
    6. Santin T,
    7. Bicalho RC
    . 2015. Oral administration of Faecalibacterium prausnitzii decreased the incidence of severe diarrhea and related mortality rate and increased weight gain in preweaned dairy heifers. PLoS One 10:e0145485. doi:10.1371/journal.pone.0145485.
    OpenUrlCrossRef
  16. 16.↵
    1. Liang G,
    2. Malmuthuge N,
    3. Bao H,
    4. Stothard P,
    5. Griebel PJ,
    6. Guan LL
    . 2016. Transcriptome analysis reveals regional and temporal differences in mucosal immune system development in the small intestine of neonatal calves. BMC Genomics 17:602. doi:10.1186/s12864-016-2957-y.
    OpenUrlCrossRef
  17. 17.↵
    1. Lozupone CA,
    2. Stombaugh JI,
    3. Gordon JI,
    4. Jansson JK,
    5. Knight R
    . 2012. Diversity, stability, and resilience of the human gut microbiota. Nature 489:220–230. doi:10.1038/nature11550.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Minot S,
    2. Bryson A,
    3. Chehoud C,
    4. Wu GD,
    5. Lewis JD,
    6. Bushman FD
    . 2013. Rapid evolution of the human gut virome. Proc Natl Acad Sci U S A 110:12450–12455. doi:10.1073/pnas.1300833110.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Yeoman CJ,
    2. Ishaq SL,
    3. Bichi E,
    4. Olivo SK,
    5. Lowe J,
    6. Aldridge BM
    . 2018. Biogeographical differences in the influence of maternal microbial sources on the early successional development of the bovine neonatal gastrointestinal tract. Sci Rep 8:3197. doi:10.1038/s41598-018-21440-8.
    OpenUrlCrossRef
  20. 20.↵
    1. Zhou M,
    2. Chen Y,
    3. Griebel PJ,
    4. Guan LL
    . 2014. Methanogen prevalence throughout the gastrointestinal tract of pre-weaned dairy calves. Gut Microbes 5:628–638. doi:10.4161/19490976.2014.969649.
    OpenUrlCrossRef
  21. 21.↵
    1. Hoffmann C,
    2. Dollive S,
    3. Grunberg S,
    4. Chen J,
    5. Li H,
    6. Wu GD,
    7. Lewis JD,
    8. Bushman FD
    . 2013. Archaea and fungi of the human gut microbiome: correlations with diet and bacterial residents. PLoS One 8:e66019. doi:10.1371/journal.pone.0066019.
    OpenUrlCrossRef
  22. 22.↵
    1. De Waele V,
    2. Van den Broeck F,
    3. Huyse T,
    4. McGrath G,
    5. Higgins I,
    6. Speybroeck N,
    7. Berzano M,
    8. Raleigh P,
    9. Mulcahy GM,
    10. Murphy TM
    . 2013. Panmictic structure of the Cryptosporidium parvum population in Irish calves: influence of prevalence and host movement. Appl Environ Microbiol 79:2534–2541. doi:10.1128/AEM.03613-12.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Malmuthuge N,
    2. Li M,
    3. Chen Y,
    4. Fries P,
    5. Griebel PJ,
    6. Baurhoo B,
    7. Zhao X,
    8. Guan LL
    . 2012. Distinct commensal bacteria associated with ingesta and mucosal epithelium in the gastrointestinal tracts of calves and chickens. FEMS Microbiol Ecol 79:337–347. doi:10.1111/j.1574-6941.2011.01220.x.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Yatsunenko T,
    2. Rey FE,
    3. Manary MJ,
    4. Trehan I,
    5. Dominguez-Bello MG,
    6. Contreras M,
    7. Magris M,
    8. Hidalgo G,
    9. Baldassano RN,
    10. Anokhin AP,
    11. Heath AC,
    12. Warner B,
    13. Reeder J,
    14. Kuczynski J,
    15. Caporaso JG,
    16. Lozupone CA,
    17. Lauber C,
    18. Clemente JC,
    19. Knights D,
    20. Knight R,
    21. Gordon JI
    . 2012. Human gut microbiome viewed across age and geography. Nature 486:222–227. doi:10.1038/nature11053.
    OpenUrlCrossRefPubMedWeb of Science
  25. 25.↵
    1. Jami E,
    2. Israel A,
    3. Kotser A,
    4. Mizrahi I
    . 2013. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J 7:1069–1079. doi:10.1038/ismej.2013.2.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Benson AK,
    2. Kelly SA,
    3. Legge R,
    4. Ma F,
    5. Low SJ,
    6. Kim J,
    7. Zhang M,
    8. Oh PL,
    9. Nehrenberg D,
    10. Hua K,
    11. Kachman SD,
    12. Moriyama EN,
    13. Walter J,
    14. Peterson DA,
    15. Pomp D
    . 2010. Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc Natl Acad Sci U S A 107:18933–18938. doi:10.1073/pnas.1007028107.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Goodrich JK,
    2. Waters JL,
    3. Poole AC,
    4. Sutter JL,
    5. Koren O,
    6. Blekhman R,
    7. Beaumont M,
    8. Van Treuren W,
    9. Knight R,
    10. Bell JT,
    11. Spector TD,
    12. Clark AG,
    13. Ley RE
    . 2014. Human genetics shape the gut microbiome. Cell 159:789–799. doi:10.1016/j.cell.2014.09.053.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    1. Goodrich JK,
    2. Davenport ER,
    3. Beaumont M,
    4. Jackson MA,
    5. Knight R,
    6. Ober C,
    7. Spector TD,
    8. Bell JT,
    9. Clark AG,
    10. Ley RE
    . 2016. Genetic determinants of the gut microbiome in UK twins. Cell Host Microbe 19:731–743. doi:10.1016/j.chom.2016.04.017.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Turpin W,
    2. Espin-Garcia O,
    3. Xu W,
    4. Silverberg MS,
    5. Kevans D,
    6. Smith MI,
    7. Guttman DS,
    8. Griffiths A,
    9. Panaccione R,
    10. Otley A,
    11. Xu L,
    12. Shestopaloff K,
    13. Moreno-Hagelsieb G, GEM Project Research Consortium,
    14. Paterson AD,
    15. Croitoru K
    . 2016. Association of host genome with intestinal microbial composition in a large healthy cohort. Nat Genet 48:1413–1417. doi:10.1038/ng.3693.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Arumugam M,
    2. Raes J,
    3. Pelletier E,
    4. Le Paslier D,
    5. Yamada T,
    6. Mende DR,
    7. Fernandes GR,
    8. Tap J,
    9. Bruls T,
    10. Batto J-M,
    11. Bertalan M,
    12. Borruel N,
    13. Casellas F,
    14. Fernandez L,
    15. Gautier L,
    16. Hansen T,
    17. Hattori M,
    18. Hayashi T,
    19. Kleerebezem M,
    20. Kurokawa K,
    21. Leclerc M,
    22. Levenez F,
    23. Manichanh C,
    24. Nielsen HB,
    25. Nielsen T,
    26. Pons N,
    27. Poulain J,
    28. Qin J,
    29. Sicheritz-Ponten T,
    30. Tims S,
    31. Torrents D,
    32. Ugarte E,
    33. Zoetendal EG,
    34. Wang J,
    35. Guarner F,
    36. Pedersen O,
    37. de Vos WM,
    38. Brunak S,
    39. Doré J,
    40. Antolín M,
    41. Artiguenave F,
    42. Blottiere HM,
    43. Almeida M,
    44. Brechot C,
    45. Cara C,
    46. Chervaux C,
    47. Cultrone A,
    48. Delorme C,
    49. Denariaz G,
    50. Dervyn R,
    51. Foerstner KU,
    52. Friss C,
    53. van de Guchte M,
    54. Guedon E,
    55. Haimet F,
    56. Huber W,
    57. van Hylckama-Vlieg J,
    58. Jamet A,
    59. Juste C,
    60. Kaci G,
    61. Knol J,
    62. Lakhdari O,
    63. Layec S,
    64. Le Roux K,
    65. Maguin E,
    66. Mérieux A,
    67. Melo Minardi R,
    68. M’rini C,
    69. Muller J,
    70. Oozeer R,
    71. Parkhill J,
    72. Renault P,
    73. Rescigno M,
    74. Sanchez N,
    75. Sunagawa S,
    76. Torrejon A,
    77. Turner K,
    78. Vandemeulebrouck G,
    79. Varela E,
    80. Winogradsky Y,
    81. Zeller G,
    82. Weissenbach J,
    83. Ehrlich SD,
    84. Bork P
    . 2011. Enterotypes of the human gut microbiome. Nature 473:174–180. doi:10.1038/nature09944.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Weiss G,
    2. Rasmussen S,
    3. Nielsen Fink L,
    4. Jarmer H,
    5. Nøhr Nielsen B,
    6. Frøkiær H
    . 2010. Bifidobacterium bifidum actively changes the gene expression profile induced by Lactobacillus acidophilus in murine dendritic cells. PLoS One 2010. 5:e11065. doi:10.1371/journal.pone.0011065.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. O’Flaherty S,
    2. Klaenhammer TR
    . 2012. Influence of exposure time on gene expression by human intestinal epithelial cells exposed to Lactobacillus acidophilus. Appl Environ Microbiol 78:5028–5032. doi:10.1128/AEM.00504-12.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Gabryszewski SJ,
    2. Bachar O,
    3. Dyer KD,
    4. Percopo CM,
    5. Killoran KE,
    6. Domachowske JB,
    7. Rosenberg HF
    . 2011. Lactobacillus-mediated priming of the respiratory mucosa protects against lethal pneumovirus infection. J Immunol 186:1151–1161. doi:10.4049/jimmunol.1001751.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. von Schillde MA,
    2. Hormannsperger G,
    3. Weiher M,
    4. Alpert CA,
    5. Hahne H,
    6. Bauerl C,
    7. van Huynegem K,
    8. Steidler L,
    9. Hrncir T,
    10. Pérez-Martínez G,
    11. Kuster B,
    12. Haller D
    . 2012. Lactocepin secreted by Lactobacillus exerts anti-inflammatory effects by selectively degrading proinflammatory chemokines. Cell Host Microbe 11:387–396. doi:10.1016/j.chom.2012.02.006.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Rincon M,
    2. Flavell RA,
    3. Davis RJ
    . 2001. Signal transduction by MAP kinases in T lymphocytes. Oncogene 20:2490–2497. doi:10.1038/sj.onc.1204382.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Neis EP,
    2. Dejong CH,
    3. Rensen SS
    . 2015. The role of microbial amino acid metabolism in host metabolism. Nutrients 7:2930–2946. doi:10.3390/nu7042930.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Ferrario C,
    2. Duranti S,
    3. Milani C,
    4. Mancabelli L,
    5. Lugli GA,
    6. Turroni F,
    7. Mangifesta M,
    8. Viappiani A,
    9. Ossiprandi MC,
    10. van Sinderen D,
    11. Ventura M
    . 2015. Exploring amino acid auxotrophy in Bifidobacterium bifidum PRL2010. Front Microbiol 6:1331. doi:10.3389/fmicb.2015.01331.
    OpenUrlCrossRef
  38. 38.↵
    1. Carbonero F,
    2. Benefiel AC,
    3. Alizadeh-Ghamsari AH,
    4. Gaskins HR
    . 2012. Microbial pathways in colonic sulfur metabolism and links with health and disease. Front Physiol 3:448. doi:10.3389/fphys.2012.00448.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Babidge W,
    2. Millard S,
    3. Roediger W
    . 1998. Sulfides impair short chain fatty acid beta-oxidation at acyl-CoA dehydrogenase level in colonocytes: implications for ulcerative colitis. Mol Cell Biochem 181:117–124.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    1. Ulluwishewa D,
    2. Anderson RC,
    3. McNabb WC,
    4. Moughan PJ,
    5. Wells JM,
    6. Roy NC
    . 2011. Regulation of tight junction permeability by intestinal bacteria and dietary components. J Nutr 141:769–776. doi:10.3945/jn.110.135657.
    OpenUrlAbstract/FREE Full Text
  41. 41.↵
    1. Hsieh CY,
    2. Osaka T,
    3. Moriyama E,
    4. Date Y,
    5. Kikuchi J,
    6. Tsuneda S
    . 2015. Strengthening of the intestinal epithelial tight junction by Bifidobacterium bifidum. Physiol Rep 3:e12327. doi:10.14814/phy2.12327.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Grootjans J,
    2. Thuijls G,
    3. Verdam F,
    4. Derikx JPM,
    5. Lenaerts K,
    6. Buurman WA
    . 2010. Non-invasive assessment of barrier integrity and function of the human gut. World J Gastrointest Surg 2:61–69. doi:10.4240/wjgs.v2.i3.61.
    OpenUrlCrossRefPubMed
  43. 43.↵
    USDA. 2010. Dairy 2007, heifer calf health and management practices on US dairy operations, p 122–131. USDA/APHIS/VS/CEAH, Fort Collins, CO.
  44. 44.↵
    1. Smith G
    . 2015. Antimicrobial decision making for enteric diseases of cattle. Vet Clin North Am Food Anim Pract 31:47–60. doi:10.1016/j.cvfa.2014.11.004.
    OpenUrlCrossRef
  45. 45.↵
    1. Zhao T,
    2. Doyle MP,
    3. Harmon BG,
    4. Brown CA,
    5. Mueller POE,
    6. Parks AH
    . 1998. Reduction of carriage of enterohemorrhagic Escherichia coli O157:H7 in cattle by inoculation with probiotic bacteria. J Clin Microbiol 36:641–647.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Hansen CHF,
    2. Nielsen DS,
    3. Kverka M,
    4. Zakostelska Z,
    5. Klimesova K,
    6. Hudcovic T,
    7. Tlaskalova-Hogenova H,
    8. Hansen AK
    . 2012. Patterns of early gut colonization shape future immune responses of the host. PLoS One 7:e34043. doi:10.1371/journal.pone.0034043.
    OpenUrlCrossRef
  47. 47.↵
    1. Yu Z,
    2. Morrison M
    . 2004. Improved extraction of PCR-quality community DNA from digesta and fecal samples. Biotechniques 36:808–812. doi:10.2144/04365ST04.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    1. Meyer F,
    2. Paarmann D,
    3. D’Souza M,
    4. Olson R,
    5. Glass EM,
    6. Kubal M,
    7. Paczian T,
    8. Rodriguez A,
    9. Stevens R,
    10. Wilke A,
    11. Wilkening J,
    12. Edwards RA
    . 2008. The metagenomics RAST server: a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi:10.1186/1471-2105-9-386.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. White JR,
    2. Nagarajan N,
    3. Pop M
    . 2009. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol 5:e1000352. doi:10.1371/journal.pcbi.1000352.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Peña AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Knights D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. Sevinsky JR,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Liang G,
    2. Malmuthuge N,
    3. McFadden TB,
    4. Bao H,
    5. Griebel PJ,
    6. Stothard P,
    7. Guan LL
    . 2014. Potential regulatory role of microRNAs in the development of bovine gastrointestinal tract during early life. PLoS One 9:e92592. doi:10.1371/journal.pone.0092592.
    OpenUrlCrossRef
  52. 52.↵
    1. Li M,
    2. Penner GB,
    3. Hernandez-Sanabria E,
    4. Oba M,
    5. Guan LL
    . 2009. Effects of sampling location and time, and host animal on assessment of bacterial diversity and fermentation parameters in the bovine rumen. J Appl Microbiol 107:1924–1934. doi:10.1111/j.1365-2672.2009.04376.x.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Suzuki R,
    2. Shimodaira H
    . 2006. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22:1540–1542. doi:10.1093/bioinformatics/btl117.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Benjamini Y,
    2. Hochberg Y
    . 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc 57:289–300. doi:10.1111/j.2517-6161.1995.tb02031.x.
    OpenUrlCrossRef
PreviousNext
Back to top
Download PDF
Citation Tools
Taxonomic and Functional Compositions of the Small Intestinal Microbiome in Neonatal Calves Provide a Framework for Understanding Early Life Gut Health
Nilusha Malmuthuge, Guanxiang Liang, Philip J. Griebel, Le Luo Guan
Applied and Environmental Microbiology Mar 2019, 85 (6) e02534-18; DOI: 10.1128/AEM.02534-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Taxonomic and Functional Compositions of the Small Intestinal Microbiome in Neonatal Calves Provide a Framework for Understanding Early Life Gut Health
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Taxonomic and Functional Compositions of the Small Intestinal Microbiome in Neonatal Calves Provide a Framework for Understanding Early Life Gut Health
Nilusha Malmuthuge, Guanxiang Liang, Philip J. Griebel, Le Luo Guan
Applied and Environmental Microbiology Mar 2019, 85 (6) e02534-18; DOI: 10.1128/AEM.02534-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

gut microbiome
metagenomics
mucosal immune system
neonatal calves

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336