Next Article in Journal
Impact of Fetal Growth Restriction on the Neonatal Microglial Proteome in the Rat
Previous Article in Journal
How Does the Family Influence Adolescent Eating Habits in Terms of Knowledge, Attitudes and Practices? A Global Systematic Review of Qualitative Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Milk Formula Diet Alters Bacterial and Host Protein Profile in Comparison to Human Milk Diet in Neonatal Piglet Model

1
USDA-ARS, Arkansas Children’s Nutrition Center, Little Rock, AR 72202, USA
2
Department of Pediatrics, University of Arkansas for Medical Sciences, Little Rock, AR 72201, USA
3
Department of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, Little Rock, AR 72201, USA
4
Department of Biomedical Informatics, University of Arkansas for Medical Sciences, Little Rock, AR 72201, USA
5
Department of Microbiology and Immunology, University of Arkansas for Medical Sciences, Little Rock, AR 72201, USA
*
Author to whom correspondence should be addressed.
Nutrients 2021, 13(11), 3718; https://doi.org/10.3390/nu13113718
Submission received: 7 September 2021 / Revised: 18 October 2021 / Accepted: 20 October 2021 / Published: 22 October 2021
(This article belongs to the Topic Probiotics, Prebiotics and Postbiotics in Human Health)

Abstract

:
The metaproteome profiling of cecal contents collected from neonatal piglets fed pasteurized human milk (HM) or a dairy-based infant formula (MF) from postnatal day (PND) 2 to 21 were assessed. At PND 21, a subset of piglets from each group (n = 11/group) were euthanized, and cecal contents were collected for further metaproteome analysis. Cecal microbiota composition showed predominantly more Firmicutes phyla and Lachnospiraceae family in the lumen of cecum of HM-fed piglets in comparison to the MF-fed group. Ruminococcus gnavus was the most abundant species from the Firmicutes phyla in the cecal contents of the HM-fed piglets at 21 days of age. A greater number of expressed proteins were identified in the cecal contents of the HM-fed piglets relative to the MF-fed piglets. Greater abundances of proteins potentially expressed by Bacteroides spp. such as glycoside enzymes were noted in the cecal lumen of HM-fed piglets relative to the MF. Additionally, lyases associated with Lachnospiraceae family were abundant in the cecum of the HM group relative to the MF group. Overall, our findings indicate that neonatal diet impacts the gut bacterial taxa and microbial proteins prior to weaning. The metaproteomics data were deposited into PRIDE, PXD025432 and 10.6019/PXD025432.

1. Introduction

Studies have demonstrated positive health outcomes in human milk fed in comparison to formula-fed infants. There is considerable amount of evidence showing that human milk diet minimizes risk of necrotizing enterocolitis in preterm infants [1,2,3,4]. Additionally, the gut microbiota colonization can be influenced by neonatal diet. The literature suggests that breastfed infants have higher abundance of Bifidobacteria and Bacteroides than formula-fed infants [5,6,7,8,9,10]. Previously, we have reported that human milk-fed (HM) piglets had higher fecal Bacteroides abundance relative to a formula-fed group (MF), as well as a stronger immune response by which enhanced T cell proliferation in the mesenteric lymph nodes of HM-fed animals [11]. Gut and immune health are functions of both diet and gut microbes that respond to diet. Several metabolites are known to derive from the microbial metabolism throughout the intestinal regions [12,13]. For instance, indigestible carbohydrates can be fermented by distal gut bacteria (cecum and colon) to short-chain fatty acids [14], and complex human milk oligosaccharides (HMOs) are broken down by microbes in the distal gut, serving as substrates to commensal bacteria among other functions reviewed elsewhere [15]. In addition, derivatives of the tryptophan metabolism (i.e., indoles) [16] and the conversion of primary to secondary bile acids are also metabolized by distal gut microbiota [17]. Metabolomics analysis of the large intestinal contents of these piglets revealed that HM feeding resulted in greater abundance of fatty acids, polyamine derivatives, glutamic acid, and tryptophan metabolites in the distal gut of HM-fed piglets. In contrast, MF-fed piglets had greater abundance of cholesterol, bile acids, and amino acids in the distal colon at 21 days of age relative to the HM-fed group [18]. These findings might be a result from the interaction between neonatal diets and gut microbial activity.
Several approaches have demonstrated that microbiota compositional changes can be altered in response to diet [19,20,21,22,23,24,25,26,27]. However, studies were limited in terms of determining the functional relevance of the microbial changes and which components of microbiota play a role in positive health outcomes observed in human milk-fed infants. Newer technology such as metaproteomics might help to determine the microbial protein presence, abundance, and microbial community. This allows us to understand the functional role of microbiota and their interactions with host and other microbial species in an ecosystem. In addition, host proteins can be identified from the sloughed off cells of the gastrointestinal tract. The proteins provide a measure of the activity of the cells and their abundances provide a phenotype at the molecular level. Metaproteomics was used often to study environmental samples in the 2000s [28]. The first shot gun metaproteomics from human samples was conducted by Verberkmoes et al. in 2009 [29]. They identified that 30% of protein hits were associated with the host, and several microbial pathways related to carbohydrate and energy metabolism were observed. The literature is limited in terms of its understanding of the microbial protein/peptide functions, especially in neonates.
We hypothesized that cecal bacterial community and bacterial proteins act as signaling molecules to promote gut homeostasis and immune function in HM-fed piglets relative to MF-fed piglets. Thus, a metaproteomics approach was used to determine the bacterial protein expression in piglets fed either human milk or cow’s milk formula.

2. Materials and Methods

2.1. Experimental Design

An animal experiment was conducted as per the Institutional Animal Care and Use Committee approval at the University of Arkansas for Medical Sciences (UAMS Institutional Animal Care and committee 3727 and 3471). Diet composition and the experimental design were published previously [11,30]. Piglets were obtained from 4–6 sows. At 2 days of age, White Dutch Landrace Duroc male piglets were randomized into two dietary groups (n = 11/group): pasteurized human milk (HM) provided from the Mother’s Milk Bank of North Texas, or a dairy-based infant formula (MF) (milk formula; Similac Advance powder; Ross products, Abbott Laboratories, Columbus, OH). Piglets were fed to meet the nutrient requirements of growing pigs as per National Research Council (NRC) guidelines [31]. Piglets were euthanized at PND 21 to collect cecal contents. Samples were immediately frozen in liquid nitrogen and transferred to an −80 °C freezer.

2.2. Metaproteome

For the comparison of the MF-fed vs. HM-fed piglets, we employed a mass spectrometry-based metaproteomics pipe-line, which was described previously [32,33]. Briefly, steps included (1) protein isolation from cecum samples, (2) fractionation of proteins on SDS-PAGE, (3) in-gel digestion by trypsin, (4) high-resolution mass spectrometry, (5) de novo sequencing, (6) protein inference, (7) organism inference, and (8) differential abundance analyses. To identify peptides and infer proteins, we used PEAKS v.8 (Bioinformatics Solutions, Waterloo, ON, Canada) with the following settings: merge scans—left unchecked; correct precursor—mass only; filter scans—unchecked. The following parameters were used for the de novo sequencing: parent mass error tolerance—5 ppm; fragment mass error tolerance—0.5 Da; enzyme—trypsin; fixed modifications—carbamidomethylation (C); variable modifications—oxidation (M), deamidation (NQ); max variable PTM per peptide—3; report # peptides—5.
Preliminary taxonomy analysis was conducted to evaluate reproducibility of protein extraction. De novo peptide tags obtained by PEAKS for the raw MSMS spectra were filtered using the average local confidence score (ALC ≥ 80), and the filtered peptide lists were supplied to the online metaproteomics tool Unipept [34], with the following settings: equate I and L—checked, filter duplicate peptides—unchecked, advanced missed cleavage handling—checked. The taxonomy information was visualized using a tree diagram provided by Unipept. Bacterial-to-host ratios were calculated using number of de novo sequences matched to each of the two taxonomic kingdoms (Table S1). This step is required the assessment of the quality of protein extraction in order to ensure that sufficient number of bacterial proteins was extracted. After this quality control step, we proceed with the protein identification using the multi-step database strategy as implemented in Peaks Studio.
Multi-step database search strategy for protein identification: The high quality de novo tags (average local confidence score ≥ 50%) were search against a series of protein databases using the multi-step database strategy. The false discovery rate estimation as implemented in PEAKS is compatible with the multi-step searches [35]. Step 1: Uniprot/Tremble protein database (downloaded on 13 April 2020) was searched using Homo sapiens and Sus scrofa taxonomic filter (310,501 entries were searched). Unmatched de novo tags from this step were passed on to Step 2, wherein the Uniprot database was searched using bacteria, archaea, and fungi as taxonomic filters (142,741,860 entries searched). No filters were applied to the search results in these 2 first steps, apart from the de novo quality score (ALC ≥ 50%). All of the identified entries from the first two steps (≈10% estimated F.D.R at this point, 0 unique peptides allowed) were used to compile a sequence database for the final search. Step 3: The de novo tags were re-searched against the final sequence database derived from the results of the previous two steps (172,464 entries), applying stringent FDR criteria to the final result: 1% false discovery rate for peptide-to-spectrum matches (corresponding average −10lgP ≈ 25 across samples) and minimum of 1 unique peptide per protein. One unique peptide hits were further required to have −10lgP = 30 in order to be considered identified. Additional filters were applied at the next step for comparative analysis.
Differential abundance of proteins and bacteria: Spectral counts (number of tandem MS spectra that match to a given protein sequence via the database search) were used to infer differential abundant (DA) proteins and taxonomic units. At the taxonomic unit level, the spectral counts of proteins were grouped using taxonomic information in the sequence database and then were summed to obtain total spectral counts for each species in each sample. If species were not identifiable, higher taxonomic levels were used. Moreover, the identified organism had to be present in at least 4 of the independent biological replicates in either of the two conditions compared. The counts were filtered so that species with less than 10 counts in all samples, but one was removed. Then, counts were normalized to the trimmed mean of M values, a method frequently employed in RNA–Seq analysis [36]. The differential abundance analysis was performed employing Poisson–Tweedie family of distributions using tweeDE package in R [37]. Initially, data analysis for microbiota and microbial and host proteins was conducted by edgeR and DESeq2 methods with different statistical tests (i.e., Wald LRT for DESeq2 and LRT, exactTest for edgeR). Finally, Benjamini–Hochberg correction was used for multiple testing to define differentially abundant proteins and bacterial species (FDR < 0.05).

2.3. Data Accessibility

The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE [38] partner repository with the dataset identifier PXD025432 and 10.6019/PXD025432. Reviewer login details: Username: [email protected]; password: qvFTwXRs.

3. Results

Figures S1 and S2 detail these data as Venn diagrams (Bacteria Venn pairwise and protein Venn pairwise, respectively). Due to zero inflation and overdispersion observed with the data, we used the Poisson–Tweedie method, which enables direct fitting of data with heavy-tails and/or zero-inflation [32]. Heat maps show the most abundant bacterial taxa and proteins altered in cecal lumen due to neonatal feeding (Figure 1 and Figure 2, respectively). The microbial abundance and metaproteome data are listed in Table 1 and Table 2. All the bacteria and proteins identified in this study are presented in the Supplemental Datasets S1 and S2, respectively.

3.1. Microbial Taxonomy Identification in the Cecal Lumen

The bacterial abundance in the luminal cecum of HM- and MF-fed piglets is shown in Table 1. The cecum profile of the HM-fed piglets was predominantly composed of the Firmicutes phylum and of the Lachnospiraceae family, including the species Ruminococcus lactaris, Ruminococcus gnavus, and Lachnospiraceae bacterium, while the cecal lumen of the MF-fed relative to HM-fed piglets had higher abundance of the Bacteroides genera including Bacteroides clarus and Bacteroides stercoris. Additionally, the cecum of MF-fed piglets had greater abundance of the Clostridium clostridioforme (fold-change (FC) = 2.9) compared to the HM-fed group.

3.2. Bacterial Proteins Impacted by Diet Groups in the Lumen of Cecum at PND 21

Bacterial peptide profile of cecal contents of HM- or MF-fed piglets at PND 21 are shown in Table 2. A greater number of bacterial proteins were identified in the HM-fed group relative to the MF piglets. The top 10 bacterial proteins identified in the lumen of cecum of MF group were from the phylum Bacteroidetes, including species from Bacteroides and Phocaeicola genus. Peptides derived from Phocaeicola vulgatus (Bacteroides vulgatus) included RagB/SusD family nutrient uptake outer membrane proteins as well as malate dehydrogenase. In fact, proteins associated with Phocaeicola vulgatus were also identified in the cecal contents of the HM-fed piglets; however, a greater diverse pool of peptides were observed relative to the MF group. For instance, galactose oxidase, sialidase, tetracycline resistance protein, and chaperonin were peptides associated with Phocaeicola vulgatus that had higher abundance in the cecum of the HM group compared to the MF group. Additionally, the Lacl family transcriptional regulator associated with the Firmicutes bacterium was greater in the cecal lumen of HM (FC = 3) relative to the MF group. L-fucose isomerase, D-ribose pyranase, and chaperonin Firmicutes bacterium associated-proteins were greater in the cecal contents of HM compared to MF-fed piglets. The aldehyde-lyase fructose-1,6-bisphosphate aldolase had greater abundance in the cecum of the HM group relative to the MF group. Additionally, this enzyme was associated with different species in the cecum of HM group such as Lachnospiraceae bacterium, Ruminococcus gnavus, and uncultured Ruminococcus sp. The abundance of phosphotransferase acetate kinase was also greater in the cecal contents of HM group, and it was associated with both species Lachnospiraceae bacterium and Clostridium sp. D5.

3.3. Host Proteins Identified in the Cecal Contents at PND 21

Host proteins expressed in the cecal contents of HM-fed versus MF-fed piglets at PND 21 is shown in Table S2. Briefly, the human proteins N-sulphoglucosamine sulphohydrolase, epididymis secretory sperm binding protein, alpha-1-antitrypsin, and lactotransferrin were greater (FC > 5) in the cecum of HM-fed piglets compared to the MF group. In contrast, the MF-fed piglets had greater porcine proteins such as secreted folate binding protein, folate_rec domain-containing protein, and transthyretin relative to the HM-fed group.

4. Discussion

This study used a porcine model due to the similarities in the anatomy and physiology of the digestive tract between pigs and humans [39,40]. Previous studies found that different protein sources such as bovine milk, hydrolyzed bovine milk, and soybean formula did not change intestinal trypsin and chymotrypsin and the absorption of nitrogen in the small and large intestine in 3-week-old piglets, similar to the human infants [39]. Furthermore, it has been demonstrated that 3-week-old piglets are suitable for studying parameters of digestion and absorption relative to 3-month-old infants [40]. In our previous study, we observed that MF-fed piglets had an increased microbial diversity and richness across the luminal regions compared to the HM-fed group [26], which is in agreement with microbiota composition findings in infants that have shown higher microbial richness in formula-fed infants [41,42]. Thus, the gut related outcomes from the current study have the potential to be translated to infants consuming human milk or formula.
Metaproteome analysis of gut microbiota are typically conducted with fecal samples, and the latter constitutes a significant amount of microbial biomass in feces, which can reflect the intestinal conditions. However fecal samples are a mixture of microbiota from all intestinal regions, and the piglet model provided the opportunity to measure the specific bioregion of the gut (i.e., cecal contents). In addition, it has been demonstrated that the main microbial fermentation of both carbohydrate and protein occur in the cecum, suggesting a microbiota role in putrefaction [43]; thus, cecal luminal contents were considered for this study. Future studies are needed to determine bioregional differences in bacterial protein expression and its impact on gut health.
Bifidobacterium and Bacteroides are the most abundant genera observed in breastfed infants [24,44], while in formula-fed infants, Bifidobacterium and Bacteroides have been identified in similar levels [9]. Bacteroides vulgatus had persistent abundance from birth up to 4 months of age in the infant gut [45]. Bacteroides vulgatus and Bacteroides dorei abundances have been reported to increase in the feces of infants at 6 months of age [46], while in the adult gut microbiota community of healthy individuals, these species within the Bacteroides genera are the most predominant [47]. Additionally, Bacteroides abundance in the human gut has been associated with the maintenance of a healthy gut [48]. In line with these observations, we previously reported a higher abundance of Bacteroides in the feces of HM-fed piglets relative to the formula-fed group [11]. In the current study, metaproteomic analysis revealed greater abundance of specific bacterial peptides belonging to the Bacteroides vulgatus in the cecal contents of HM-fed piglets relative to MF-fed group at 21 days of age. Interestingly, studies have shown that Bacteroides vulgatus can grow in the presence of human milk oligosaccharides (HMO), as well as metabolize these complex carbohydrates [49,50]. Moreover, proteins associated with Bacteroides vulgatus has been identified in stool samples of breastfed infants at 2–3 months of age [51]. Interestingly, Bacteroides spp. promote Treg cell development [52,53], and it has been shown that infants with decreased allergic colitis had increased Bacteroides spp. in their stool [54] suggesting the role of these species in promoting immune responses and homeostasis in the gut. This further suggests the role of Bacteroides spp. in cell-mediated immunity, but it is yet to be determined how this impacts antibody and humoral immune response.
Recently, Bifidobacterium abundance has been reported to decrease in the feces of infants from 6 to 12 months of age, while Lachnospiraceae abundance increased [46]. Interestingly, in this study, alongside the Bacteroides vulgatus-associated proteins, a greater number of enzymes related to the Lachnospiraceae family were identified in the luminal cecum of HM-fed piglets relative to the MF group. Studies demonstrated that bacteria within the Lachnospiraceae family, in particular Ruminococcus gnavus, has the ability to produce iso-bile acids, and such metabolites can favor the growth of Bacteroides [55,56]. Recently, the PreventADALL cohort study evaluated the microbial composition and the metaproteome of 100 mother–child pairs from Norway and Sweden [57]. Within the Firmicutes phyla, the predominant species identified in the feces of 12 months old infants was Ruminococcus gnavus, and glycoside hydrolases were the enzymes associated with such species [57]. These findings are in agreement with the greater abundance of enzymes involved in the degradation of sialic acids observed in our study, including fructose-1,6-bisphosphate aldolase potentially expressed by Ruminococcus gnavus in the luminal cecum of HM-fed piglets relative to MF group at PND 21. Ruminococcus gnavus and such aldolase enzymes have been reported to metabolize sialic acid to N-acetylmannosamine [58]. Indeed, sialic acids are a family of carbon sugar acids present in human milk in a rich source of oligosaccharide-bound sialic acid [59].
Furthermore, human milk feeding enhanced the expression of glycoside hydrolases, nutrient uptake proteins, and transporters, which are common enzymes involved in the HMO consumption by Bacteroides and Bifidobacteria [60,61]. Therefore, it is plausible to acknowledge that milk glycans such as HMOs promote the growth of beneficial bacteria in HM-fed piglets. Additionally, the different hydrolases identified in the HM-fed piglets potentially expressed by Bacteroides might benefit the immune system, pending mechanistic data.
We acknowledge that in the current study, numbers of microbial taxa and proteins identified were smaller than in usual metagenome analysis. Kleiner et al. described very elegantly the limitations of metaproteomics data acquisition with current mass spectra and how that limits in-depth community analysis [62]. In addition, the metaproteomics analysis requires a well-curated sequence database to assign the proteins to individual microbial species. This potential mismatch between identified proteins and their assignment to bacterial species is evident when one compares hierarchical clustering using bacterial abundance (Figure 1) and protein abundance (Figure 2). The clustering using protein abundance shows clear separation between the HM and MF groups, while clustering using bacterial abundance has one outlier from the HM group (the right-most HM column in the middle of the heatmap, Figure 1). We argue that the limited number of microbial species and hits to microbial proteins in the current study is a result of combination of technology and data acquisition. This is expected to improve in the next few years with better mass spectra technology and metaproteome database tools.
The current study was conducted at controlled environmental (housed at the vivarium) and isocaloric diets for both HM and MF groups. The human milk used in this study was composed of milk samples ranging from 2 to 12 months of lactation, and different components were added to the diet to maintain the growing piglet nutrient requirements. Additionally, the piglets were enrolled in the study at 2 days of age. Thus, this study lacks data on the colostrum intake. These limitations might introduce variation in the milk composition and might affect the luminal microbiota composition and protein expression of the gut microbiota. Moreover, piglets were from 4–6 sows, by which genetic differences could cause some variation.

5. Conclusions

In summary, we observed a 5.4-fold increase in the relative abundance of Ruminococcus gnavus in the cecal microbiota composition of HM-fed piglets relative to the MF-fed piglets at 21 days of age. This bacterial abundance was also associated with the expression of glycolytic enzymes in the cecum lumen of HM-fed piglets compared to the MF group. Furthermore, the greater number of proteins potentially expressed by Bacteroides vulgatus observed in the cecal contents of HM-fed piglets relative to the MF-fed group at 21 days of age might be associated with the ingestion of bioactive components of human milk (i.e., HMOs metabolized by this gut bacteria) and possibly promotes immune function. Overall, our findings highlight the association between gut microbiota composition upon different neonatal diets with the peptides and enzymes originated from this interaction. We believe that further research in the field of metaproteomics might be crucial to understanding the establishment of key gut colonizers and the overall effect on the host metabolism and immune system.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13113718/s1, Figure S1: Venn diagram representing the bacterial pairwise comparisons. Figure S2: Venn diagram representing the proteins pairwise comparisons. Dataset S1: All bacteria identified with the current approach. Dataset S2: All proteins identified in the experimental samples. Table S1: Host/bacteria ratios identified in this study. Table S2: Host protein expression identified in the cecal contents of piglets fed HM or MF at postnatal day 21.

Author Contributions

L.Y. conceptualized the study; B.L.Z., G.V.G., Y.R., S.B. and S.G.M. conducted data analyses and technical aspects of the study; A.K.B. conducted the animal experiment; F.R., B.L.Z., G.V.G. and L.Y. interpreted the data; F.R., B.L.Z. and L.Y. drafted the manuscript and had the primary responsibility for the manuscript revision. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by NIGMS (P20GM121293 to L.Y.) and partly by USDA-ARS (6026-51000-010-06S to L.Y.). In addition, L.Y. was also supported by NIAID (R21AI146521). G.V.G and Y.R. were supported by NIH IDeA Networks of Biomedical Research Excellence (INBRE) grant (P20GM103429). B.L.Z., G.V.G., S.G.M., and S.B. were also supported by NIGMS (P20GM121293).

Institutional Review Board Statement

This study was conducted according to the approval of the Institutional Animal Care and Use Committee at the University of Arkansas for Medical Sciences (UAMS Institutional Animal Care and Committee Protocol# 3727 and 3471).

Informed Consent Statement

Not applicable.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD025432 and 10.6019/PXD025432. Reviewer login details: Username: [email protected] Password: qvFTwXRs.

Acknowledgments

We thank Matt Ferguson, Jessica Besancon, Mallory Jayroe, Bobby Fay, and Trae Pittman for their help with the animal trials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davis, J.A.; Baumgartel, K.; Morowitz, M.J.; Giangrasso, V.; Demirci, J.R. The Role of Human Milk in Decreasing Necrotizing Enterocolitis Through Modulation of the Infant Gut Microbiome: A Scoping Review. J. Hum. Lact 2020, 36, 647–656. [Google Scholar] [CrossRef]
  2. Bering, S.B. Human Milk Oligosaccharides to Prevent Gut Dysfunction and Necrotizing Enterocolitis in Preterm Neonates. Nutrients 2018, 10, 1461. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, B.; Xiu, W.; Dai, Y.; Yang, C. Protective effects of different doses of human milk on neonatal necrotizing enterocolitis. Medicine 2020, 99, e22166. [Google Scholar] [CrossRef]
  4. Cañizo Vázquez, D.; Salas García, S.; Izquierdo Renau, M.; Iglesias-Platas, I. Availability of Donor Milk for Very Preterm Infants Decreased the Risk of Necrotizing Enterocolitis without Adversely Impacting Growth or Rates of Breastfeeding. Nutrients 2019, 11, 1895. [Google Scholar] [CrossRef] [Green Version]
  5. Balmer, S.E.; Wharton, B.A. Diet and faecal flora in the newborn: Breast milk and infant formula. Arch. Dis. Child. 1989, 64, 1672–1677. [Google Scholar] [CrossRef] [Green Version]
  6. Bezirtzoglou, E.; Tsiotsias, A.; Welling, G.W. Microbiota profile in feces of breast- and formula-fed newborns by using fluorescence in situ hybridization (FISH). Anaerobe 2011, 17, 478–482. [Google Scholar] [CrossRef] [PubMed]
  7. Fallani, M.; Amarri, S.; Uusijarvi, A.; Adam, R.; Khanna, S.; Aguilera, M.; Gil, A.; Vieites, J.M.; Norin, E.; Young, D.; et al. Determinants of the human infant intestinal microbiota after the introduction of first complementary foods in infant samples from five European centres. Microbiology 2011, 157, 1385–1392. [Google Scholar] [CrossRef] [Green Version]
  8. Fanaro, S.; Chierici, R.; Guerrini, P.; Vigi, V. Intestinal microflora in early infancy: Composition and development. Acta Paediatr. Suppl. 2003, 91, 48–55. [Google Scholar] [CrossRef] [PubMed]
  9. Harmsen, H.J.; Wildeboer-Veloo, A.C.; Raangs, G.C.; Wagendorp, A.A.; Klijn, N.; Bindels, J.G.; Welling, G.W. Analysis of intestinal flora development in breast-fed and formula-fed infants by using molecular identification and detection methods. J. Pediatr. Gastroenterol. Nutr. 2000, 30, 61–67. [Google Scholar] [CrossRef] [PubMed]
  10. Mackie, R.I.; Sghir, A.; Gaskins, H.R. Developmental microbial ecology of the neonatal gastrointestinal tract. Am. J. Clin. Nutr. 1999, 69, 1035S–1045S. [Google Scholar] [CrossRef]
  11. Miklavcic, J.J.; Badger, T.M.; Bowlin, A.K.; Matazel, K.S.; Cleves, M.A.; LeRoith, T.; Saraf, M.K.; Chintapalli, S.V.; Piccolo, B.D.; Shankar, K.; et al. Human Breast-Milk Feeding Enhances the Humoral and Cell-Mediated Immune Response in Neonatal Piglets. J. Nutr. 2018, 148, 1860–1870. [Google Scholar] [CrossRef] [Green Version]
  12. Davila, A.M.; Blachier, F.; Gotteland, M.; Andriamihaja, M.; Benetti, P.H.; Sanz, Y.; Tomé, D. Intestinal luminal nitrogen metabolism: Role of the gut microbiota and consequences for the host. Pharm. Res. 2013, 68, 95–107. [Google Scholar] [CrossRef] [PubMed]
  13. Neis, E.P.; Dejong, C.H.; Rensen, S.S. The role of microbial amino acid metabolism in host metabolism. Nutrients 2015, 7, 2930–2946. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. den Besten, G.; van Eunen, K.; Groen, A.K.; Venema, K.; Reijngoud, D.J.; Bakker, B.M. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J. Lipid Res. 2013, 54, 2325–2340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Carr, L.E.; Virmani, M.D.; Rosa, F.; Munblit, D.; Matazel, K.S.; Elolimy, A.A.; Yeruva, L. Role of Human Milk Bioactives on Infants’ Gut and Immune Health. Front. Immunol. 2021, 12, 604080. [Google Scholar] [CrossRef] [PubMed]
  16. Zelante, T.; Iannitti, R.G.; Cunha, C.; De Luca, A.; Giovannini, G.; Pieraccini, G.; Zecchi, R.; D’Angelo, C.; Massi-Benedetti, C.; Fallarino, F.; et al. Tryptophan catabolites from microbiota engage aryl hydrocarbon receptor and balance mucosal reactivity via interleukin-22. Immunity 2013, 39, 372–385. [Google Scholar] [CrossRef] [Green Version]
  17. Ridlon, J.M.; Kang, D.J.; Hylemon, P.B. Bile salt biotransformations by human intestinal bacteria. J. Lipid Res. 2006, 47, 241–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Rosa, F.; Matazel, K.S.; Bowlin, A.K.; Williams, K.D.; Elolimy, A.A.; Adams, S.H.; Bode, L.; Yeruva, L. Neonatal Diet Impacts the Large Intestine Luminal Metabolome at Weaning and Post-Weaning in Piglets Fed Formula or Human Milk. Front. Immunol. 2020, 11, 607609. [Google Scholar] [CrossRef]
  19. Abrahamsson, T.R.; Jakobsson, H.E.; Andersson, A.F.; Bjorksten, B.; Engstrand, L.; Jenmalm, M.C. Low gut microbiota diversity in early infancy precedes asthma at school age. Clin. Exp. Allergy 2014, 44, 842–850. [Google Scholar] [CrossRef] [Green Version]
  20. Albenberg, L.G.; Wu, G.D. Diet and the intestinal microbiome: Associations, functions, and implications for health and disease. Gastroenterology 2014, 146, 1564–1572. [Google Scholar] [CrossRef] [Green Version]
  21. Alizadeh, A.; Akbari, P.; Difilippo, E.; Schols, H.A.; Ulfman, L.H.; Schoterman, M.H.; Garssen, J.; Fink-Gremmels, J.; Braber, S. The piglet as a model for studying dietary components in infant diets: Effects of galacto-oligosaccharides on intestinal functions. Br. J. Nutr. 2016, 115, 605–618. [Google Scholar] [CrossRef] [Green Version]
  22. Ardeshir, A.; Narayan, N.R.; Mendez-Lagares, G.; Lu, D.; Rauch, M.; Huang, Y.; Van Rompay, K.K.; Lynch, S.V.; Hartigan-O’Connor, D.J. Breast-fed and bottle-fed infant rhesus macaques develop distinct gut microbiotas and immune systems. Sci. Transl. Med. 2014, 6, 252ra120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Arora, T.; Seyfried, F.; Docherty, N.G.; Tremaroli, V.; le Roux, C.W.; Perkins, R.; Backhed, F. Diabetes-associated microbiota in fa/fa rats is modified by Roux-en-Y gastric bypass. ISME J. 2017, 11, 2035–2046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Avershina, E.; Storro, O.; Oien, T.; Johnsen, R.; Pope, P.; Rudi, K. Major faecal microbiota shifts in composition and diversity with age in a geographically restricted cohort of mothers and their children. FEMS Microbiol. Ecol. 2014, 87, 280–290. [Google Scholar] [CrossRef] [PubMed]
  25. Azad, M.B.; Konya, T.; Maughan, H.; Guttman, D.S.; Field, C.J.; Chari, R.S.; Sears, M.R.; Becker, A.B.; Scott, J.A.; Kozyrskyj, A.L.; et al. Gut microbiota of healthy Canadian infants: Profiles by mode of delivery and infant diet at 4 months. CMAJ 2013, 185, 385–394. [Google Scholar] [CrossRef] [Green Version]
  26. Brink, L.R.; Matazel, K.; Piccolo, B.D.; Bowlin, A.K.; Chintapalli, S.V.; Shankar, K.; Yeruva, L. Neonatal Diet Impacts Bioregional Microbiota Composition in Piglets Fed Human Breast Milk or Infant Formula. J. Nutr. 2019, 149, 2236–2246. [Google Scholar] [CrossRef] [PubMed]
  27. Brink, L.R.; Mercer, K.E.; Piccolo, B.D.; Chintapalli, S.V.; Elolimy, A.; Bowlin, A.K.; Matazel, K.S.; Pack, L.; Adams, S.H.; Shankar, K.; et al. Neonatal diet alters fecal microbiota and metabolome profiles at different ages in infants fed breast milk or formula. Am. J. Clin. Nutr. 2020, 111, 1190–1202. [Google Scholar] [CrossRef]
  28. Wilmes, P.; Bond, P.L. Metaproteomics: Studying functional gene expression in microbial ecosystems. Trends Microbiol. 2006, 14, 92–97. [Google Scholar] [CrossRef]
  29. Verberkmoes, N.C.; Russell, A.L.; Shah, M.; Godzik, A.; Rosenquist, M.; Halfvarson, J.; Lefsrud, M.G.; Apajalahti, J.; Tysk, C.; Hettich, R.L.; et al. Shotgun metaproteomics of the human distal gut microbiota. ISME J. 2009, 3, 179–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Elolimy, A.A.; Washam, C.; Byrum, S.; Chen, C.; Dawson, H.; Bowlin, A.K.; Randolph, C.E.; Saraf, M.K.; Yeruva, L. Formula Diet Alters the Ileal Metagenome and Transcriptome at Weaning and during the Postweaning Period in a Porcine Model. mSystems 2020, 5, e00457-20. [Google Scholar] [CrossRef]
  31. Council, N.R. Nutrient Requirements of Swine: Eleventh Revised Edition; The National Academies Press: Washington, DC, USA, 2012; p. 420. [Google Scholar]
  32. Karaduta, O.; Glazko, G.; Dvanajscak, Z.; Arthur, J.; Mackintosh, S.; Orr, L.; Rahmatallah, Y.; Yeruva, L.; Tackett, A.; Zybailov, B. Resistant starch slows the progression of CKD in the 5/6 nephrectomy mouse model. Physiol. Rep. 2020, 8, e14610. [Google Scholar] [CrossRef] [PubMed]
  33. Zybailov, B.L.; Glazko, G.V.; Rahmatallah, Y.; Andreyev, D.S.; McElroy, T.; Karaduta, O.; Byrum, S.D.; Orr, L.; Tackett, A.J.; Mackintosh, S.G.; et al. Metaproteomics reveals potential mechanisms by which dietary resistant starch supplementation attenuates chronic kidney disease progression in rats. PLoS ONE 2019, 14, e0199274. [Google Scholar] [CrossRef] [Green Version]
  34. Mesuere, B.; Debyser, G.; Aerts, M.; Devreese, B.; Vandamme, P.; Dawyndt, P. The Unipept metaproteomics analysis pipeline. Proteomics 2015, 15, 1437–1442. [Google Scholar] [CrossRef]
  35. Zhang, J.; Xin, L.; Shan, B.; Chen, W.; Xie, M.; Yuen, D.; Zhang, W.; Zhang, Z.; Lajoie, G.A.; Ma, B. PEAKS DB: De novo sequencing assisted database search for sensitive and accurate peptide identification. Mol. Cell Proteom. 2012, 11, M111.010587. [Google Scholar] [CrossRef] [Green Version]
  36. Robinson, M.D.; Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11, R25. [Google Scholar] [CrossRef] [Green Version]
  37. Esnaola, M.; Puig, P.; Gonzalez, D.; Castelo, R.; Gonzalez, J.R. A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments. BMC Bioinform. 2013, 14, 254. [Google Scholar] [CrossRef] [Green Version]
  38. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef] [PubMed]
  39. Miller, E.R.; Ullrey, D.E. The pig as a model for human nutrition. Annu Rev. Nutr. 1987, 7, 361–382. [Google Scholar] [CrossRef]
  40. Moughan, P.J.; Birtles, M.J.; Cranwell, P.D.; Smith, W.C.; Pedraza, M. The piglet as a model animal for studying aspects of digestion and absorption in milk-fed human infants. World Rev. Nutr. Diet. 1992, 67, 40–113. [Google Scholar] [CrossRef] [PubMed]
  41. Fan, W.; Huo, G.; Li, X.; Yang, L.; Duan, C. Impact of diet in shaping gut microbiota revealed by a comparative study in infants during the six months of life. J. Microbiol. Biotechnol. 2014, 24, 133–143. [Google Scholar] [CrossRef]
  42. Azad, M.B.; Konya, T.; Persaud, R.R.; Guttman, D.S.; Chari, R.S.; Field, C.J.; Sears, M.R.; Mandhane, P.J.; Turvey, S.E.; Subbarao, P.; et al. Impact of maternal intrapartum antibiotics, method of birth and breastfeeding on gut microbiota during the first year of life: A prospective cohort study. Bjog 2016, 123, 983–993. [Google Scholar] [CrossRef] [PubMed]
  43. Heinzmann, S.S.; Schmitt-Kopplin, P. Deep metabotyping of the murine gastrointestinal tract for the visualization of digestion and microbial metabolism. J. Proteome Res. 2015, 14, 2267–2277. [Google Scholar] [CrossRef] [PubMed]
  44. Jost, T.; Lacroix, C.; Braegger, C.P.; Chassard, C. New insights in gut microbiota establishment in healthy breast fed neonates. PLoS ONE 2012, 7, e44595. [Google Scholar] [CrossRef] [PubMed]
  45. Ferretti, P.; Pasolli, E.; Tett, A.; Asnicar, F.; Gorfer, V.; Fedi, S.; Armanini, F.; Truong, D.T.; Manara, S.; Zolfo, M.; et al. Mother-to-Infant Microbial Transmission from Different Body Sites Shapes the Developing Infant Gut Microbiome. Cell Host Microbe 2018, 24, 133–145.e135. [Google Scholar] [CrossRef] [PubMed]
  46. Nilsen, M.; Lokmic, A.; Angell, I.L.; Lødrup Carlsen, K.C.; Carlsen, K.H.; Haugen, G.; Hedlin, G.; Jonassen, C.M.; Marsland, B.J.; Nordlund, B.; et al. Fecal Microbiota Nutrient Utilization Potential Suggests Mucins as Drivers for Initial Gut Colonization of Mother-Child-Shared Bacteria. Appl. Environ. Microbiol. 2021, 87, e02201-20. [Google Scholar] [CrossRef]
  47. Arumugam, M.; Raes, J.; Pelletier, E.; Le Paslier, D.; Yamada, T.; Mende, D.R.; Fernandes, G.R.; Tap, J.; Bruls, T.; Batto, J.-M.; et al. Enterotypes of the human gut microbiome. Nature 2011, 473, 174–180. [Google Scholar] [CrossRef]
  48. Wexler, A.G.; Goodman, A.L. An insider’s perspective: Bacteroides as a window into the microbiome. Nat. Microbiol. 2017, 2, 17026. [Google Scholar] [CrossRef] [Green Version]
  49. Ward, R.E.; Niñonuevo, M.; Mills, D.A.; Lebrilla, C.B.; German, J.B. In vitro fermentation of breast milk oligosaccharides by Bifidobacterium infantis and Lactobacillus gasseri. Appl. Environ. Microbiol. 2006, 72, 4497–4499. [Google Scholar] [CrossRef] [Green Version]
  50. Marcobal, A.; Barboza, M.; Froehlich, J.W.; Block, D.E.; German, J.B.; Lebrilla, C.B.; Mills, D.A. Consumption of human milk oligosaccharides by gut-related microbes. J. Agric. Food Chem. 2010, 58, 5334–5340. [Google Scholar] [CrossRef] [Green Version]
  51. Cortes, L.; Wopereis, H.; Tartiere, A.; Piquenot, J.; Gouw, J.W.; Tims, S.; Knol, J.; Chelsky, D. Metaproteomic and 16S rRNA Gene Sequencing Analysis of the Infant Fecal Microbiome. Int. J. Mol. Sci. 2019, 20, 1430. [Google Scholar] [CrossRef] [Green Version]
  52. Round, J.L.; Lee, S.M.; Li, J.; Tran, G.; Jabri, B.; Chatila, T.A.; Mazmanian, S.K. The Toll-like receptor 2 pathway establishes colonization by a commensal of the human microbiota. Science 2011, 332, 974–977. [Google Scholar] [CrossRef] [Green Version]
  53. Round, J.L.; Mazmanian, S.K. Inducible Foxp3+ regulatory T-cell development by a commensal bacterium of the intestinal microbiota. Proc. Natl. Acad. Sci. USA 2010, 107, 12204–12209. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, S.X.; Li, Y.H.; Dai, W.K.; Li, X.S.; Qiu, C.Z.; Ruan, M.L.; Zou, B.; Dong, C.; Liu, Y.H.; He, J.Y.; et al. Fecal microbiota transplantation induces remission of infantile allergic colitis through gut microbiota re-establishment. World J. Gastroenterol. 2017, 23, 8570–8581. [Google Scholar] [CrossRef]
  55. Tailford, L.E.; Owen, C.D.; Walshaw, J.; Crost, E.H.; Hardy-Goddard, J.; Le Gall, G.; de Vos, W.M.; Taylor, G.L.; Juge, N. Discovery of intramolecular trans-sialidases in human gut microbiota suggests novel mechanisms of mucosal adaptation. Nat. Commun. 2015, 6, 7624. [Google Scholar] [CrossRef] [Green Version]
  56. Devlin, A.S.; Fischbach, M.A. A biosynthetic pathway for a prominent class of microbiota-derived bile acids. Nat. Chem. Biol. 2015, 11, 685–690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Nilsen, M.; Madelen Saunders, C.; Leena Angell, I.; Arntzen, M.; Lødrup Carlsen, K.C.; Carlsen, K.H.; Haugen, G.; Hagen, L.H.; Carlsen, M.H.; Hedlin, G.; et al. Butyrate Levels in the Transition from an Infant- to an Adult-Like Gut Microbiota Correlate with Bacterial Networks Associated with Eubacterium Rectale and Ruminococcus Gnavus. Genes 2020, 11, 1245. [Google Scholar] [CrossRef] [PubMed]
  58. Kentache, T.; Thabault, L.; Deumer, G.; Haufroid, V.; Frédérick, R.; Linster, C.L.; Peracchi, A.; Veiga-da-Cunha, M.; Bommer, G.T.; Van Schaftingen, E. The metalloprotein YhcH is an anomerase providing N-acetylneuraminate aldolase with the open form of its substrate. J. Biol. Chem. 2021, 296, 100699. [Google Scholar] [CrossRef] [PubMed]
  59. Nakano, T.; Sugawara, M.; Kawakami, H. Sialic acid in human milk: Composition and functions. Acta Paediatr. Taiwan 2001, 42, 11–17. [Google Scholar]
  60. Marcobal, A.; Sonnenburg, J.L. Human milk oligosaccharide consumption by intestinal microbiota. Clin. Microbiol. Infect. 2012, 18, 12–15. [Google Scholar] [CrossRef] [Green Version]
  61. Kitaoka, M. Bifidobacterial enzymes involved in the metabolism of human milk oligosaccharides. Adv. Nutr. 2012, 3, 422s–429s. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Kleiner, M. Metaproteomics: Much More than Measuring Gene Expression in Microbial Communities. mSystems 2019, 4, e00115-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Metaproteomics revealed distinct microbiome profiles in human milk-fed and milk formula-fed neonatal piglets. Heatmap visualization of the top 31 significantly abundant bacteria is shown. The two phenotypes are indicated by the bottom labels: HM and MF (n = 11/group). Unsupervised clustering of the metaproteomics-derived bacterial abundance resulted in HM and MF piglets being grouped together, except for one outlier from the HM group.
Figure 1. Metaproteomics revealed distinct microbiome profiles in human milk-fed and milk formula-fed neonatal piglets. Heatmap visualization of the top 31 significantly abundant bacteria is shown. The two phenotypes are indicated by the bottom labels: HM and MF (n = 11/group). Unsupervised clustering of the metaproteomics-derived bacterial abundance resulted in HM and MF piglets being grouped together, except for one outlier from the HM group.
Nutrients 13 03718 g001
Figure 2. Metaproteomics-derived protein profiles uniquely defined in HM and MF-fed neonatal piglets. Heatmap constructed from differentially abundant proteins is shown. The two dietary groups are indicated by the bottom labels: HM and MF (n = 11/group). Unsupervised clustering of bacterial and host protein resulted in HM and MF piglets being grouped together. Host and human proteins are identified with *, red—human, blue—pig.
Figure 2. Metaproteomics-derived protein profiles uniquely defined in HM and MF-fed neonatal piglets. Heatmap constructed from differentially abundant proteins is shown. The two dietary groups are indicated by the bottom labels: HM and MF (n = 11/group). Unsupervised clustering of bacterial and host protein resulted in HM and MF piglets being grouped together. Host and human proteins are identified with *, red—human, blue—pig.
Nutrients 13 03718 g002
Table 1. Cecal bacterial profile of piglets (n = 11/group) fed either with human milk (HM) or milk formula (MF) at PND 21.
Table 1. Cecal bacterial profile of piglets (n = 11/group) fed either with human milk (HM) or milk formula (MF) at PND 21.
Bacterial Organism 1HM 2SDMF 2SDlog2 FC 3p-Value 4
Bacteroides clarus6.927.195.868.4−3.8<0.0001
Bacteroides stercoris CAG:1203.413.344.533.9−3.70.0028
Bacteroides sp. AF35-2239.249.2266.5194.2−2.80.0076
Clostridium clostridioforme CAG:1323.06.221.721.7−2.90.0076
Bacteroides stercoris ATCC 431834.015.735.930.4−3.20.0112
Bacteroides sp. AM44-1957.662.1312.0229.9−2.40.0131
Firmicutes bacterium CAG:1291.86.918.012.8−3.30.0159
Firmicutes bacterium CAG:1372.511.278.982.5−5.00.0164
Firmicutes bacterium CAG:1762.23.817.016.0−3.00.0197
Bacteroides clarus YIT 120563.413.332.030.7−3.20.0375
(Clostridium) clostridioforme 90A76.77.926.927.2−2.00.0375
Clostridium bolteae (strain ATCC BAA-613/WAL 16351)2.56.416.215.0−2.70.0375
(Ruminococcus) lactaris135.1113.03.610.95.2<0.0001
Firmicutes bacterium OM07-1182.151.32.55.45.10.0001
Ruminococcus gnavus303.5202.57.115.45.40.0001
Firmicutes bacterium AM41-5BH199.0119.815.615.93.70.0001
Mediterraneibacter sp. gm00247.739.04.57.73.40.0009
Firmicutes bacterium AM43-11BH275.4158.219.518.63.80.0012
(Ruminococcus) torques ATCC 2775673.267.43.211.24.50.0028
Firmicutes bacterium AM31-12AC209.2169.66.410.85.00.0067
Firmicutes bacterium CAG:212116.857.13.58.35.00.0087
Faecalicatena orotica50.338.16.212.63.00.0131
Selenomonas ruminantium233.3140.541.9106.82.50.0152
Dorea sp. CAG:31713.111.02.05.12.70.0485
Hespellia stercorisuis DSM 1548012.811.70.00.0NA0.0152
Lachnospiraceae bacterium AM40-2BH139.2116.633.051.52.10.0152
Faecalicatena contorta164.5122.533.141.52.30.0174
Lachnospiraceae bacterium 2_1_58FAA88.473.10.30.98.30.0197
Lachnospiraceae bacterium TM07-2AC117.5112.50.00.0NA0.0197
Dorea sp. OM07-544.036.50.00.0NA0.0389
Lachnospiraceae bacterium OF09-33XD18.820.50.00.0NA0.0482
1 The raw spectral counts matching to the identified bacterial species were analyzed using tweeDEseq package in Bioconductor. 2 HM and MF columns indicate mean value of the total spectral counts followed by the standard deviation of the mean (SD). 3 Log2 FC is the log2 of the HM to MF ratio. 4 Benjamini–Hochberg correction for multiple testing was applied to adjust p-values.
Table 2. Bacterial protein profile of cecal content of piglets (n = 11/group) fed either with human milk (HM) or milk formula (MF) at PND 21.
Table 2. Bacterial protein profile of cecal content of piglets (n = 11/group) fed either with human milk (HM) or milk formula (MF) at PND 21.
Uniprot_IDOrganismProtein NameHM 1
(sum)
MF 1
(sum)
HM 2MF 2Log2 FC 3p-Value 4
A0A412VME7Phocaeicola vulgatus (Bacteroides vulgatus)RagB/SusD family nutrient uptake outer membrane protein515733.447.5−3.80.002
A0A412VME7Phocaeicola doreiRagB/SusD family nutrient uptake outer membrane protein515733.445.9−3.80.008
R9H865Bacteroides vulgatus CL09T03C04Uncharacterized protein17831113.432.9−3.30.03
C6Z771Bacteroides sp. 4_3_47FAASusD family protein513533.432.9−3.30.03
A0A1Y3Z4G2Bacteroides clarusPolyribonucleotide nucleotidyltransferase (EC 2.7.7.8) (polynucleotide phosphorylase) (PNPase)443863.735.4−3.20.05
A0A397WJ15Phocaeicola vulgatus (Bacteroides vulgatus)Beta-galactosidase (EC 3.2.1.23)1076148.452.2−2.60.02
A0A415SZ11Phocaeicola vulgatus (Bacteroides vulgatus)Malate dehydrogenase (EC 1.1.1.37)1008199792.3180.3−10.03
C3RBX6Bacteroides vulgatus CL09T03C0460 kDa chaperonin (GroEL protein) (protein Cpn60)32541195245.9118.31.10.03
C3RBX6Bacteroides sp. AM18-960 kDa chaperonin (GroEL protein) (protein Cpn60)32541195245.9118.31.10.03
A0A069SFK4Bacteroides vulgatus str. 3975 RP4TonB-linked outer membrane, SusC/RagA family protein1370468103.237.41.50.02
A0A069SKV9Bacteroides vulgatus str. 3975 RP4Tetracycline resistance protein TetQ37921001302.996.41.70.03
A0A069SKV9Bacteroides sp. AF32-15BHTetracycline resistance protein TetQ37921001302.996.41.70.03
A0A1H6GDE8Selenomonas ruminantiumPhosphoenolpyruvate carboxykinase (ATP) (PCK) (PEP carboxykinase) (PEPCK) (EC 4.1.1.49)2383519238.143.52.50.04
A0A174QRN1Phocaeicola vulgatus (Bacteroides vulgatus)DUF1735 domain-containing protein (galactose oxidase) (EC 3.2.1.18)1712913.122.70.03
I9UCD9Phocaeicola vulgatus (Bacteroides vulgatus)SusC/RagA family TonB-linked outer membrane protein137046932.34.72.80.03
A0A417FRM0Firmicutes bacterium AM31-12ACLacI family transcriptional regulator2513524.52.83.10.02
A0A397WGI7Phocaeicola vulgatus (Bacteroides vulgatus)RagB/SusD family nutrient uptake outer membrane protein (starch-binding associating with outer membrane)108711196.810.53.20.02
A0A316A2W0Faecalicatena contortaPhosphoglycerate kinase (EC 2.7.2.3)114296107.18.73.60.002
R9H865Bacteroides vulgatus dnLKV7Uncharacterized protein1783111128.510.53.60.01
A0A417FRM0Firmicutes bacterium AM43-11BHLacI family transcriptional regulator2513521.61.53.80.05
R5ICG7Firmicutes bacterium CAG:124L-fucose isomerase (FucIase) (EC 5.3.1.25) (6-deoxy-L-galactose isomerase)937084.90NA<0.0001
A0A3C1BZY4Clostridiales bacteriumCarbon monoxide dehydrogenase (EC 1.2.7.4)0535053.7NA0.007
F0Z4P0Clostridium sp. D5Acetate kinase (EC 2.7.2.1) (acetokinase)296025.90NA0.008
A0A374BZR2Firmicutes bacterium AM31-12ACD-ribose pyranase (EC 5.4.99.62)155013.10NA0.02
A0A417BM01Firmicutes bacterium AM43-11BHUPF0210 protein DW928_02850148012.90NA0.02
A0A2N5NQW2Ruminococcus gnavusClass II fructose-1,6-bisphosphate aldolase (EC 4.1.2.13) (fructose-1,6-bisphosphate aldolase, class II)777078.20NA0.02
A0A494ZWH6Ruminococcus sp. B05UPF0210 protein D8Q48_03220148012.90NA0.02
A0A374BZR2Mediterraneibacter sp. gm002D-ribose pyranase (EC 5.4.99.62)155013.10NA0.02
A0A316A2W0Dorea sp. CAG:105Phosphoglycerate kinase (EC 2.7.2.3)11429647.30NA0.02
F0Z4P0Lachnospiraceae bacterium 1_4_56FAAAcetate kinase (EC 2.7.2.1) (acetokinase)296018.60NA0.03
A0A3A6G758Lachnospiraceae bacterium TM07-2ACClass II fructose-1,6-bisphosphate aldolase (EC 4.1.2.13)776078.10NA0.03
A0A374BZR2Ruminococcus sp. AM42-11D-ribose pyranase (EC 5.4.99.62)1550120NA0.04
C3RBX6Firmicutes bacterium AM43-11BH60 kDa chaperonin (GroEL protein) (protein Cpn60)32541195107.20NA0.05
1 The raw spectral counts matching to the identified proteins were analyzed using tweeDEseq package in Bioconductor. 2 HM and MF columns indicate mean value of the total spectral counts. 3 Log2 FC is the log2 of the HM-to-MF ratio. 4 Benjamini–Hochberg correction for multiple testing was applied to adjust p-values.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rosa, F.; Zybailov, B.L.; Glazko, G.V.; Rahmatallah, Y.; Byrum, S.; Mackintosh, S.G.; Bowlin, A.K.; Yeruva, L. Milk Formula Diet Alters Bacterial and Host Protein Profile in Comparison to Human Milk Diet in Neonatal Piglet Model. Nutrients 2021, 13, 3718. https://doi.org/10.3390/nu13113718

AMA Style

Rosa F, Zybailov BL, Glazko GV, Rahmatallah Y, Byrum S, Mackintosh SG, Bowlin AK, Yeruva L. Milk Formula Diet Alters Bacterial and Host Protein Profile in Comparison to Human Milk Diet in Neonatal Piglet Model. Nutrients. 2021; 13(11):3718. https://doi.org/10.3390/nu13113718

Chicago/Turabian Style

Rosa, Fernanda, Boris L. Zybailov, Galina V. Glazko, Yasir Rahmatallah, Stephanie Byrum, Samuel G. Mackintosh, Anne K. Bowlin, and Laxmi Yeruva. 2021. "Milk Formula Diet Alters Bacterial and Host Protein Profile in Comparison to Human Milk Diet in Neonatal Piglet Model" Nutrients 13, no. 11: 3718. https://doi.org/10.3390/nu13113718

APA Style

Rosa, F., Zybailov, B. L., Glazko, G. V., Rahmatallah, Y., Byrum, S., Mackintosh, S. G., Bowlin, A. K., & Yeruva, L. (2021). Milk Formula Diet Alters Bacterial and Host Protein Profile in Comparison to Human Milk Diet in Neonatal Piglet Model. Nutrients, 13(11), 3718. https://doi.org/10.3390/nu13113718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop