Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Effects of alpha-(1,2)-fucosyltransferase genotype variants on plasma metabolome, immune responses and gastrointestinal bacterial enumeration of pigs pre- and post-weaning

  • Ann-Sofie Riis Poulsen ,

    Contributed equally to this work with: Ann-Sofie Riis Poulsen, Diana Luise

    Roles Data curation, Formal analysis, Writing – original draft

    Affiliation Aarhus University, Faculty of Science and Technology, Department of Animal Science, Tjele, Denmark

  • Diana Luise ,

    Contributed equally to this work with: Ann-Sofie Riis Poulsen, Diana Luise

    Roles Data curation, Formal analysis, Investigation, Writing – original draft

    Affiliation Department of Agricultural and Food Sciences (DISTAL), Alma Mater Studiorum—University of Bologna, Bologna, Italy

  • Mihai Victor Curtasu,

    Roles Formal analysis

    Affiliation Aarhus University, Faculty of Science and Technology, Department of Animal Science, Tjele, Denmark

  • Sugiharto Sugiharto,

    Roles Formal analysis

    Affiliation Diponegoro University, Faculty of Animal and Agricultural Sciences, Semarang, Central Java, Indonesia

  • Nuria Canibe,

    Roles Conceptualization, Project administration, Resources, Supervision, Writing – review & editing

    Affiliation Aarhus University, Faculty of Science and Technology, Department of Animal Science, Tjele, Denmark

  • Paolo Trevisi,

    Roles Writing – review & editing

    Affiliation Department of Agricultural and Food Sciences (DISTAL), Alma Mater Studiorum—University of Bologna, Bologna, Italy

  • Charlotte Lauridsen

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing

    charlotte.lauridsen@anis.au.dk

    Affiliation Aarhus University, Faculty of Science and Technology, Department of Animal Science, Tjele, Denmark

Abstract

In pigs, the alpha-(1,2) fucosyltransferase (FUT1) gene has been highlighted for its properties in controlling the intestinal expression of enterotoxigenic E. coli (ETEC) F18 receptors; a pathogen causing edema disease and post-weaning diarrhoea. In this study, we hypothesized that pigs with different genotypes (ETEC F18 resistant (FUT1AA) versus susceptible (FUT1AG)) differed in following systemic and enteric responses: growth performance, plasma metabolic profiles, expression of candidate genes for intestinal mucosal homeostasis and immunity, number of selected bacteria and the concentration of short-chain fatty acids (SCFA) in faeces and digesta in piglets pre and post-weaning, and on the ETEC F18 adherence ex vivo. Genotype had the strongest impact on plasma metabolomic profile on day 7 and 28 of age. FUT1AG piglets had higher level of N-methyl-2-pyrrolidinone, hippuric acid, oxindole, and 3-oxo-5-beta-chol-7-en-24-oic acid on day 7, and a higher level of guanosine on day 28 than that in the FUT1AA piglets. FUT1AA piglets had a higher level of betaine on day 7 and 3-methylguanine on day 28. On day 34 of age, the FUT1AA pigs had higher levels of S-2-hydroxyglutarate, L-phenylalanine, tauroursodeoxycholic acid and an undetermined PC/LysoPC, while Ile Glu Phe Gly peptide and genistein 5-O-glucuronide, and PC (18:0/0:0) were at higher levels in the FUT1AG piglets. FUT1 genotype did not affect the growth performance and expression of candidate genes. FUT1AG piglets had a higher number of haemolytic bacteria in faeces and in digesta than that in FUT1AA at 34 days of age. The colonic acetic acid concentration was highest in FUT1AG piglets. FUT1 genotype may influence not only the expression of E. coli F18 receptors but could potentially impact the gut homeostasis and metabotype of piglets pre and post-weaning. Further investigations on the relation between FUT1 genotype and these aspects including the intestinal commensal microbiota will expand the knowledge on factors affecting the intestinal ecosystem.

Introduction

The gut microbiota composition and host genetics background are suggested to play key roles in the development of inflammatory intestinal diseases [1]. In pigs, the alpha-(1,2) fucosyltransferase (FUT1) gene has been highlighted for its properties in controlling the intestinal expression of enterotoxigenic E. coli (ETEC) F18 receptors; a pathogen causing edema disease and post-weaning diarrhoea. In humans, the fucosyltransferase 2 (FUT2) gene is responsible for H-antigen expression in the intestine, and this structure has been shown to be particularly relevant in pathogen adhesion and for the susceptibility to chronic inflammation in the intestine, e.g., inflammatory bowel disease (IBD), Crohn’s disease (CD) and ulcerative colitis (UC) [23]. Indeed, genome wide associate studies have associated the mutation in the FUT2 gene to the risk of developing CD and UC [23]. In pigs, both the FUT1 and FUT2 genes are expressed in the intestines and the FUT1 is considered to be orthologous to the human H gene which catalyses the addition of fucose to a terminal galactose in an alpha 1,2-linkage [4].

Sophisticated genetic tools have extensively been used to genetically select pigs with favourable production traits such as increased growth rate, carcass weight, number of piglets per litter, and improved meat quality, while less focus has been given to animal health and robustness. However, the growing concern regarding bacterial antibiotic resistance has increased the demand for alternative therapeutic or preventive options when dealing with infectious diseases. Molecular genetics enables development of pig selection schemes that can include genetic markers associated with innate disease resistance. Escherichia coli (E. coli) F18 is a highly pathogenic bacterium causing post-weaning diarrhoea (PWD) and edema disease (ED) in young pigs, and until now, successful methods to control this bacterium remains to be developed [4,5]. Breeding F18-resistant piglets has been studied as a potential preventative strategy to limit PWD and ED [6]. E. coli F18 recognizes and adheres to highly specific receptors located on the small intestinal brush border. Previous studies have suggested that the F18 receptor has a glycophospholipid structure and that its expression is genetically controlled by FUT1 gene located near the blood group inhibitor S on chromosome 6 [7,8]. Meijerink and co-authors [9] proposed that the single point mutation G>A at base-pair 307 in the FUT1 gene influences piglet susceptibility to E. coli F18 diarrhoea, classifying piglets with the FUT1AA genotype as resistant animals and piglets with the FUT1AG and FUT1GG genotypes as susceptible. The FUT1AA genotype leads to a lower level of FUT1 and FUT2 genes expression, which encode for galactoside 2-L-fucosyltransferase enzyme (FUT2), and thus reduce the quantity of this enzyme [4]. In humans, the FUT2 enzyme catalyses the addition of terminal alpha (1,2) fucose residues, producing the H type 1 carbohydrate expressed on the surface of epithelial cells and in mucosal secretions of secretor individuals [10]. The fucosylation of carbohydrate structures expressed on cell surfaces of the gut by glycoproteins, glycolipids, and proteoglycans have been associated with biological processes such as inflammation, host–pathogen interactions, and tumor metastasis [11,12]. A recent study by Hesselager et al. (2016) showed that there are significant differences in the O-linked glycans of the intestinal mucosal proteins between FUT1AA and FUT1AG pigs [13]. The FUT1AG lead a higher level of O-glycan structures that were fucosylated to show H-antigens [13]. Accordingly, a previous study in mice [14] showed that the intestinal alpha (1,2) fucosylation in FUT2-null mice, which are characterized by a null expression of FUT2 gene and were obtained by replacing FUT2 with the bacterial reporter gene lacZ [15], had a lower level of fucosylation in the intestinal mucosa. Additional studies in mice showed that the intestinal alpha-(1,2) fucosylation represents a significant energy source for some bacteria, suggesting that alpha-(1,2) fucosylated glycans may contribute to the establishment and maintenance of the commensal microbial community [16].

Nevertheless, so far there exists only limited knowledge of FUT1 polymorphism and pig-gut microbiota interactions. Both, genotype and the gut microbiota, are reported to be intrinsic contributors in regulating the gut mucosal maturation, the intestinal changes occurring during the weaning transition, and more generally, the overall health of the host [1719]. In the present study we, therefore, developed the hypothesis that FUT1 genotype (E. coli F18 susceptible versus resistant) would affect the host metabolism, immune response, and the gastrointestinal microbiota composition in pigs during the suckling and post-weaning period. The dual purposes of the study were to characterise the metabolomic profile and immunological parameters of the gut epithelium of FUT1AA and FUT1AG piglets and to investigate how FUT1 genotype influences the establishment of the gastrointestinal microbiota, using classical culture techniques, in 5 to 34 days old piglets.

Methods

The experiment was carried out at the experimental facility at the Department of Animal Science (Foulum, Aarhus University). Animals were reared following Danish guidelines for animal care and use. The protocol was approved by Danish Animal Experiments Inspectorate, Ministry of Food, Agriculture and Fisheries, Danish Veterinary and Food Administration (Protocol Number: 2012-15-2934-00125).

Animals, study design, and FUT1 genotyping

A total of 17 piglets were included in the study. The animals were obtained by breeding two FUT1AG sows with the same FUT1AA boar ((Danish Landrace x Yorkshire breed) x Duroc), providing FUTAG and FUT1AA, but no FUT1GG, piglets. FUT1AG and FUT1GG piglets are considered equally susceptible to enterotoxigenic Escherichia coli F18 [5]. An ear tissue sample was collected shortly after birth and used for FUT1 genotyping. In brief, the genotyping procedure included a PCR amplification of the FUT1 gene product harbouring the mutation of interest. Subsequently, the PCR fragments were digested with the HinP1I restriction enzyme and run on a 1.5% agarose gel according to Bao et al. [20]. Genetic testing revealed that ten piglets were of the sensitive genotype FUT1AG and seven piglets were of the resistant genotype FUT1AA. (3 FUT1AG and 3 FUT1GG in one litter; 7 FUT1AG and 4 FUT1GG in the other litter).

The piglets were raised with their dam until the day of weaning (day 28), and had free access to water and creep feed (without added zinc oxide) from day 21 of age. After weaning, piglets were housed together in pens according to litter and fed a standard post-weaning diet (without added zinc oxide). The temperature in the weaning unit was 23°C and piglets had free access to rooting material. All animals were individually weighed once a week from day 7 of age, and the average daily weight gain (ADWG) was calculated. On day 34 after birth animals were sacrificed by a trained staff using a penetrating captive bolt followed by exsanguination at our facility, to avoid stress due to transportation.

Sample collection

Blood samples were collected in EDTA-containing vacutainers (Vacuette, Greiner Bio-One GmbH, Kremsmünster, Austria) from the jugular vein on day 7, 28, and 34 after birth. The blood samples were immediately centrifuged at 3,000 × g for 10 min, and the plasma was stored at –80°C for further analyses of the metabolic profile. Faecal samples were collected directly from the rectum at 5, 7, 14, 21, 28, and 34 days of age. Subsamples were stored at -20°C for organic acid analysis. Bacterial enumerations by plating were performed on a subsample of fresh faeces. Due to the small size of the piglets, adequate quantities of faeces for microbial analyses and of blood for metabolomics analysis could not be sampled from all piglets at all ages.

All piglets were sacrificed at 34 days of age. The abdomen was incised and the gastrointestinal tract removed. Luminal contents (digesta) from the stomach, last third of distal small intestine, caecum, and middle part of the colon were sampled immediately and subsamples were stored at -20°C for organic acid analysis. Bacterial enumerations by plating were performed on fresh subsamples. After removal of the lumen content, the epithelium was rinsed and mucosal samples were obtained from four gastrointestinal segments (stomach, ileum, caecum, and mid colon). The samples were stored at -80°C for further candidate gene expression analyses. In addition, 15 cm of distal small intestine were collected in order to investigate the adherence of E. coli O138:F18 to the mucosa using the Porcine Intestinal Organ Culture (PIOC) procedure [21].

LC-MS metabolomics determination

The LC-MS metabolomic profile of the plasma samples was analysed using a Dionex UltiMate 3000 (Dionex, Sunnyvale, CA, USA) ultra-performance liquid chromatography (UPLC) system, coupled to a MaXis impact Quadrupole Time-of-Flight (QTOF) mass spectrometer (Bruker Daltonics GmbH, Bremen, Germany) operating in positive and negative electrospray ionization mode (ESI+, ESI-). The analytical column was an Acquity UPLC HSS T3 (2.1 x 100 mm, 1.8 μm; Waters, Milford, MA, USA), kept at 30°C. Sample preparation and the LC-MS analysis were performed as described by Ingerslev et al. [22]. In order to evaluate the performance of the analytical system, blank samples (100% acetonitrile) and pooled plasma samples from the pigs at each time point were used as quality controls and reinjected for every six samples.

Data acquisition and control of UPLC/MS was done using oTOF control v.3.4 and HyStar v.3.2 (Bruker Daltonics GmbH, Bremen, Germany). Data pre-processing was performed by calibrating the mass spectrum using Data Analysis v.4.0 (Bruker Daltonics GmbH, Bremen, Germany) followed by the Find Molecular Features (FMF) tool. Parameters for compound detection with FMF were set to 5 for S/N threshold, 8 spectra for compound length, and a smoothing width of 2 for the peaks. Other alterations of the data included shortening of the chromatogram retention times to exclude noise and carry-over compounds, the mass range was kept between m/z 50 and 1000, and the most common adducts in positive and negative mode were subtracted.

Metabolite identification

Metabolites resulting statistically different between the experimental groups were identified using METLIN (http://masspec.scripps.edu/), the Human Metabolome Database (HMDB, http://www.hmdb.ca/), and LIPID MAPS database (http://www.lipidmaps.org/) searches. Confirmation of compounds was done by comparison of retention times and MS/MS fragmentation patterns to those found in METLIN or HMDB, or to commercial standard compounds when available. Metabolites were reported as described by Sumner et al. [23] using the four level identification system: 1. Identified compounds; 2. Putatively annotated compounds; 3. Putatively characterized compound classes; 4. Unknown compounds.

Gene expression analysis

Total RNA was extracted from the samples of the four tissues using the NucleoSpin RNA kit (Macherey-Nagel) according to the manufacturer’s protocols. The quality and quantity of RNA were evaluated using agarose gel electrophoresis and a spectrophotometer (NanoDrop ND-1000; Saveen Werner), respectively. The sample concentration was adjusted to 200 ng/μl, and the RNA then converted into complementary DNA using the High Capacity cDNA Reverse Transcription Kit (Invitrogen) according to the manufacturer’s instructions. The qPCR of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Interleukin 10 (IL10); Tumor Necrosis Factor Alpha (TNF-α), cyclo-oxygenase-2 (COX-2), Zona Occludens 1 (ZO-1), and Occludin (OCLN) genes were run in 384-well plates by mixing 2 μl of complementary DNA with 8μl of a mix containing 5 pmol of each primer, 1pmol of the probe and 5μl TaqMan Master Mix (catalogue no. 4324018; Applied Biosystems). The sequence (5’-3’) of primers and probes for GAPDH, TNF-α, IL10, and COX-2 genes are reported in Sugiharto et al. [24], while for the ZO-1 and OCLN genes, the primers and probe from 20 x TaqMan Gene Expression Assay, catalogue no. 4351372 (Applied Biosystems), were used. The analysis was performed in duplicates for the target genes and in triplicates for the housekeeping gene (GAPDH) under standard amplification conditions determined for the Applied Biosystem ViiA 7 Real-Time PCR system (Life Technologies). The data were normalized to the expression of the housekeeping gene. Relative quantitation of gene expression was calculated using the 2-ΔΔCt method.

Dry matter and organic acid analysis

Dry matter content of digesta was determined by freeze-drying (ScanVac Coolsafe 55, Labogene ApS, Lynge, Denmark). The concentration of short-chain fatty acids and lactic acid in faeces and digesta were measured as previously described by Canibe et al. [25].

Microbiological enumerations by plating

Faecal samples (approximately 1 g) were transferred to plastic bags and 5 or 10 ml pre-reduced salt medium was added. The content was homogenised in a Smasher paddle blender (bioMérieux Industry, USA) for 2 minutes. Digesta samples (approximately 5 g) were transferred to flasks containing 50 ml pre-reduced salt medium [26]. The flask content was transferred to a CO2 flushed bag and homogenized in a stomacher blender for 2 minutes. 1 ml homogenate was transferred to a Hungate tube containing 9 ml pre-reduced salt medium and 10-fold dilutions were prepared using the technique previously described by Miller and Wolin [27]. Each sample was plated on selective (and indicative) and non-selective agar plates.

Enterobacteriaceae were enumerated on MacConkey agar (Merck 1.05465) after aerobic incubation for 1 day. Haemolytic bacteria were enumerated on blood agar (Oxoid Pb5039A) after aerobic incubation for 1 day. Clostridium perfringens were enumerated using the pour-plate technique on Tryptose Sulfite Cycloserine agar (Merck 1.11972 and 1.00888) after anaerobic incubation for 1 day. Lactic acid bacteria were enumerated on de Man, Rogosa and Sharp agar (Merck 1.10660) after anaerobic incubation for 2 days. Total anaerobic bacteria were enumerated in roll tubes containing pig colon fluid-glucose-cellobiose agar and incubated for 7 days [26]. All plates and roll-tubes were incubated at 37˚C.

Porcine Intestinal Organ Culture (PIOC)

The PIOC procedure was included to investigate the adherence of E. coli O138:F18 to the mucosa of FUT1AA and FUT1AG piglets following the procedure previously described by Sugiharto et al. [21]. In short, a distal small intestinal segment (15 cm) was sampled and incubated with a solution containing E. coli O138:F18 for 1 hour. The tissue was then homogenized, diluted, and E. coli enumerated on MacConkey (Merck 1.05465) agar plates after aerobic incubation at 37°C overnight.

Statistical analysis

Weekly body weight (BW) and average daily weight gain (ADWG) were analysed using an ANOVA with a factor design including the FUT1 genotype and litter as fixed effects, and the initial body weight as a covariate.

With regard to the metabolomics data, the initial step encompassed multi-level PCA (MLPCA) carried out on the full dataset (day 7, 28, 34) and separately on the individual time datasets in order to highlight possible sub-structures in the data. The data was Pareto-scaled and potential outliers were removed after initial visual inspection of PCA plots. To evaluate metabolite differences between FUT1 genotypes, potential confounding factors were evaluated by using a linear regression model. The linear regression model included sex, litter, and weight of the animals at the three different time points (7, 28 and 34) as fixed factors. Confounding effects were removed by computing the residuals of model considering the factors described in the linear regression model. All further analyses were based on these residuals. A sparse Multilevel Partial Least Squares Discriminant Analysis (sMLPLS-DA) [28] was used to investigate differences in the plasma metabolic profile between the two FUT1 genotypes on each of the three-time point. sMLPLS-DA was applied using a dummy variable distinguishing the two genotype (AA and AG) as a response variable and the scaled subject variations as predictors. This version involves the automatically choice, using an internal 5-fold cross-validation procedure of a sparseness coefficient that controls the trade-off between the goodness-of-fit of the model and its complexity (i.e., the number of selected metabolites).

Furthermore, in order to evaluate both the stability and the effect size of the metabolites contributing with non-null coefficients to discriminant dimensions, a validation test of 1000 permutation test (permutation of predictors) coupled with a Leave One Out (LOO) procedure was applied. Metabolites having a P < 0.1 both for effect size (Psi) and stability (Pst), simultaneously, were considered stable and significant.

All statistical analyses on metabolomics data was performed in R v. 3.0.3 (R Development Core Team, 2008) using the “lme4” package [29] for the computation of metabolite residuals, the “spls” package [30] (function cv.splsda and splsda) for the sMLPLS-DA analysis and the function “prcomp” for PCA.

The impact of genetic variation FUT1AA versus FUT1AG on the expression of candidate pro-inflammatory and intestinal barrier regulation genes among the different intestinal tissues was determined using a linear model in which the FUT1 genotype, intestinal tissue, and their interaction were included as fixed effects while litter was included as random effects, and pigs as a repeated measurement.

The effects of FUT1 genotype and age on bacterial composition, organic acid concentration, and pH, were analysed by fitting the data to a linear mixed model in which genotype and age/intestinal segment were included as fixed effects, while litter was included as random effect, and pig as repeated measurement (by including random intercept terms) to account for multiple observations within the same litter or within the same pig. The fixed effects were tested using an F-test with Kenward-Roger approximation, where the reduced model was tested against the full model. When age was found to have a significant effect, a post-hoc test was performed, including Bonferroni adjustment to correct for multiple comparisons [31]. Effects were considered significant when P < 0.05, whereas when P > 0.05 but ≤ 0.10 differences were considered to indicate a trend towards a significant effect.

Statistical analysis was performed with R statistical software (R Development Core Team, 2008) using the “lm” function of “stats” package version 3.3.0 for gene expression and piglet performance data. The “lmer” function from the “lme4” package [29], the “KRmodcomp” function in the “pbkrtest” package [31], and the “multcomp" package [32] were applied for the pH, bacterial and organic acid parameters.

Results

All pigs remained healthy during the trial. Analyses of the data according to the described statistical model revealed that genotype had no significant effect on BW and ADWG during the entire period of the trial. The ADWG from day 7 to day 34 of age was 0.196 kg/day for the FUT1AG and 0.153 kg/day for the FUT1AA genotype (SEM = 0.017 kg/day).

Plasma metabolomics profile

The score plot of the Principal Component Analysis (PCA) conducted on plasma samples is shown in Fig 1. Clusters according to piglet age were clearly visible both during the suckling period (day 7 versus day 28) and after weaning (day 28 versus day 34). After outlier’s filtration (2 outliers on day 7; 1 outlier on day 28 and day 34), the datasets were composed of 12 pigs (7 FUT1AG and 5 FUT1AA) and 789 metabolites on day 7, and 15 animals (10 FUT1AG and 5 FUT1AA) with 727 metabolites on day 28, and 707 compounds on day 34. The litter and genotype effects in the single time datasets are presented in the PCA score plot (Fig 2). At 7 days of age, the effects of both litter and FUT1 genotype were evident on the metabolic profile by exploring the first (PC1) and second principal component (PC2), which explained 18.7% and 13.3% of the total variance, respectively (Fig 2A). The Fig 2B shows that FUT1AA and FUT1AG clustered into two groups along the PC2 (explained 11.7% of the total variance). The PCA showed a clear separation according to litter along the PC1 (explained 15.5% of total variance) on day 28. Litter and genotype did not have an effect on the metabolic profile on day 34. Discriminant identified compounds (P ≤ 0.01) for stability and effect size obtained from sMLPLS-DA coupled with LOO validation model in plasma metabolic profiles of FUT1AG and FUT1AA animals on day 7, 28, and 34 are listed in Table 1, while a complete list, including the unidentified compounds, is shown in S1 Table (S1 Table). Considering the discriminant identified compounds on day 7, the plasma metabolic profile of FUT1AG pigs was characterised by a higher level of N-Methyl-2-pyrrolidinone (Psi = 0.01; Pst = 0.01), hippuric acid (Psi = 0.01; Pst = 0.08), Oxindole (Psi = 0.02; Pst = 0.07) and 3-Oxo-5-beta-chol-7-en-24-oic Acid (Psi = 0.02; Pst = 0.08) compared to the profile of FUT1AA piglets, which had a higher level of Betaine (Psi = 0.03; Pst = 0.08) (Table 1). At weaning (day 28 of age), two identified compounds discriminated the FUT1 genotypes; 3-methylguanine (Psi <0.001; Pst = 0.001) characterised the FUT1AG, and guanosine (Psi<0.001; Pst = 0.02) the FUT1AA piglets. Eight identified compounds discriminated the FUT1AA and FUT1AG pigs based on their metabolome on day 34 of age. (S)-2-hydroxyglutarate (Psi = 0.02; Pst = 0.001), L-phenylalanine (Psi = 0.04; Pst = 0.07), PC 18:0/0:0 (Psi = 0.04; Pst = 0.03), and tauroursodeoxycholic acid (Psi = 0.02; Pst = 0.06 positive mode; Psi = 0.03; Pst = 0.08 negative mode) were higher in plasma of pigs with the FUT1AA genotype, while the FUT1AG piglets had higher values of an unidentified PC/ LysoPC (Psi = 0.03; Pst = 0.001), Ile Glu Phe Gly peptide (Psi = 0.04; Pst = 0.01) and genistein 5-O-glucuronide (Baicalin). (Psi = 0.02; Pst = 0.07) (Table 1).

thumbnail
Fig 1. PCA score plot of the piglets plasma metabolomic profiles at 7, 28, and 34 days of age.

The variances accounted by the principal components are shown in the axes.

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

thumbnail
Fig 2. PCA Score plot of FUT1 genotype and litter differences during suckling, weaning, and after weaning.

A): Score plot of plasma metabolites on day 7. PC1 and PC2 explain 18.67% and 13.31%, respectively; B) Metabolomics score plot on day 28 (weaning), PC1 = 15.53% and PC2 = 11.66%. C) PCA score plot after weaning (day 34), where PC1 and PC2 explain 16.68% and 12.84% of total variance, respectively.

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

thumbnail
Table 1. Discriminant plasma metabolites for FUT1AA and FUT1AG obtained from piglets at 7, 28 and 34 days of age by the sMLPLS-DA.

https://doi.org/10.1371/journal.pone.0202970.t001

Gene expression analyses—Intestinal mucosa

No significant differences between FUT1 genotypes were observed in the intestinal expression of analysed genes (S2 Table), but the COX2 expression was higher in mucosa from the mid colon compared to the mucosa from the other gut segments (P = 0.02), and OCLN expression was higher (P = 0.008) in mucosa from the ileum compared to mucosa from the caecum and mid colon. A trend of significance (P < 0.1) was observed for FUT1 genotypes and tissue interaction for ZO-1 expression.

Microbiological enumerations by plating and organic acid concentrations

Faeces.

Higher faecal numbers of haemolytic bacteria at 34 days of age (P = 0.003) and a tendency to higher numbers of lactic acid bacteria (P = 0.06) were observed in FUT1AG piglets (Table 2). The number of haemolytic bacteria was similar between day 5 and 28, whereas a significant increase was detected on day 34 (P < 0.001). An age-dependent decrease of Clostridium perfringens was found with numbers being lowest on day 28 (P < 0.001) and 34 (P < 0.001). The same was observed for lactic acid bacteria, with lower numbers on day 28 (P ≤ 0.03) and day 34 (P ≤ 0.005) compared to day 5 and 7. There was no significant effect of genotype or age on the number of Enterobacteriaceae or total anaerobic bacteria. The interaction between FUT1 genotype and age showed a trend of significance in the number of Enterobacteriaceae (P < 0.1) and a significant effect in the number of haemolytic bacteria (P = 0.03). Furthermore, no difference in the concentration of faecal short chain fatty acids between genotypes was determined, except for a tendency towards a lower concentration of butyric acid in FUT1AA piglets (P < 0.1) (Table 3). The concentration of acetic acid was higher at 34 days of age (one-week post-weaning) compared to day 5, 7, 14, and 21 of age (suckling period) (P < 0.05). The concentration of propionic acid was constant between day 5 and 21, but at 34 days of age, a significant increase in propionic acid was observed (P < 0.001). A tendency to a lower concentration of butyric acid (P < 0.01) with age was measured, whereas no effect of age was detected for valeric, iso-butyric and iso-valeric acid (data not shown).

thumbnail
Table 2. Enumeration (log cfu/g sample) of selected microbial groups in faeces from piglets at 5, 7, 14, 21, 28 and 34 days of age1.

https://doi.org/10.1371/journal.pone.0202970.t002

thumbnail
Table 3. Short-chain fatty acid concentration (mmol/kg sample) in faeces from FUT1AG and FUT1AA piglets at 5, 7, 14, 21, 28, and 34 days of age1.

https://doi.org/10.1371/journal.pone.0202970.t003

Digesta.

Genotype did not influence dry matter content or pH (S3 Table) of gastrointestinal digesta. FUT1AG piglets had a higher number of Enterobacteriaceae (P ≤ 0.02) in digesta from the distal small intestine, caecum, and mid colon; of haemolytic bacteria (P = 0.02) in the stomach, small intestine, caecum, and colon; and of total anaerobic bacteria in the distal small intestine compared to FUT1AA piglets (P = 0.004) (Table 4). The concentration of acetic acid was higher in the mid colon of FUT1AG piglets (P = 0.01) (Table 5). Digesta concentrations of propionic, butyric, lactic, valeric, and the sum of iso-butyric and iso-valeric acid were similar in both genotypes.

thumbnail
Table 4. Enumeration (log cfu/g sample) of selected microbial groups in digesta from the gastrointestinal tract of 34 days old piglets (one week post-weaning)1 from the FUT1AG and FUT1AA groups.

https://doi.org/10.1371/journal.pone.0202970.t004

thumbnail
Table 5. Organic acid concentrations in digesta (mmol/kg wet sample) from the gastrointestinal tract of 34 days old piglets (one week post-weaning) belong to FUT1AG and FUT1AA groups.

https://doi.org/10.1371/journal.pone.0202970.t005

Quantification of E. coli attachment to distal small intestinal tissue

The E. coli numbers in the small intestinal tissue were similar in both genotypes (P = 0.41). However, FUT1AG piglets had numerically higher E. coli numbers before (FUT1AG 6.3 log cfu/g versus FUT1AA 6.0 log cfu/g), and after E. coli O138:F18 inoculation (FUT1AG 8.1 log cfu/g versus FUT1AA 7.8 log cfu/g) than the FUT1AA piglets.

Discussion

In this study, a metabolomics approach was applied to evaluate the influence of FUT1 genotype on system-biological responses in healthy piglets during the suckling and post-weaning period. Bacterial counts for important bacteria groups involved in the microbial balance such as lactic acid bacteria, Enterobacteriaceae and haemolytic bacteria, gene expression analysis for genes associated with intestinal epithelial barrier (ZO-1 and OCL) [33] and mucosal gut inflammation such as the TNF-α, and COX-2, which are upregulated upon antigens stimulation, and IL-10, which encodes for an anti-inflammatory cytokine that plays a key role in immune response [34], were analysed to further add information on the impact of genotype on the host’s gastrointestinal homeostatic environment.

Age of the piglets influenced the plasma metabolomics profile and faecal bacterial enumerations, which is ascribed to the piglets physiological and immunological maturation during the initial phase of life [35,36], and to the dietary change that occurred at weaning. In our study, we found a progressive decrease in the number of lactic acid bacteria and Clostridium perfringens in the period from 5 to 34 days of age. Lactic acid bacteria, especially Lactobacillus spp., are important in maintaining intestinal homeostasis as they are capable of reducing pH by lactic acid production, compete with potential pathogens, and produce bacteriocins [37,38]. However, the number of Lactobacillus has been observed to decrease with the age and specially after the weaning transition period [3942]. The process of weaning the piglet from the sow and thereby the switch from milk to solid feed significantly increased in the number of haemolytic bacteria in faeces between 28 and 34 days of age. Indeed, the gastrointestinal microbiota is strongly influenced by diet, both by the change in milk composition during suckling period [3942], by the switch from the dam’s milk to the solid feed, and generally an increase in microbial richness and complexity occurs between birth and the post-weaning period [40,41]. The observed effect of age of piglets on the metabolomics profile and bacterial population will not be discussed in depth here, as it was not the primary focus of the present study. However, it is important to note that an interaction between FUT1 genotype and age was detected in the number of haemolytic bacteria in the faeces on day 34, suggesting that the interaction of factors can affect the number of possible pathogenic bacteria.

The presented metabolomics results shed light on the influence of FUT1 genotype to specific architecture in the pig metabolome. The PCA score plot showed relatively consistent clusters attributable to FUT1 genotype when considering each of the single time points separately, i.e., day 7, 28, and 34. The influence of genotype was clearest in the metabolome profile of young animals (day 7 and 28) compared to the plasma metabolic profile after weaning (day 34). In general, no particular pathway in the plasma metabolome was shown to be responsible for the differentiation of the two FUT1 genotypes. Nevertheless, discriminant compounds were detected and validated thought the LOO procedure in the metabolomics data from day 7, 28, and 34. On day 7, the FUT1AG piglets had higher plasma values of metabolites originating from the metabolism of phenylalanine, hippuric acid, tryptophan, and oxindole. Hippuric acid, constituted by the conjugation of benzoic and acyl glycine, is synthesised in the intestine, liver, and spleen [43], and it accumulates in the blood, causes ammoniagenesis stimulation, and it is generally associated with an increased feed consumption of phenylalanin or phenolic compounds, and with the activity of the gastrointestinal microbiota [4446]. We speculate that the higher hippuric acid amount in FUT1AG piglets could be ascribed to their higher feed intake since the sMLPLS-DA was performed on the metabolomics data corrected for BW of the animals on day 7, and further that the diet (i.e. sow milk and post weaning diet) was not different between FUT1AG and FUT1AA pigs. However, the microbial activity in the gut is also associated with the concentration of hippuric acid in the blood [45], and in humans lower urinary hippurate levels in Crohn’s disease patients have been associated with an altered gut microbial metabolism [46]. No significant difference in the number of various microbial groups in faeces of the two groups was found on day 7, though the analysis of digesta at this age might have provided different results. The oxindole, also named 1,3-Dihydroindol-2-one (KEEG C12312), is a putative tryptophan metabolite, which is involved in the benzoxazinoid biosynthesis, it is a secondary metabolite that serves as an important factor in the host’s resistance against microbial pathogens. No further information has been found regarding this compound in porcine biological systems. On the other hand, the animals with the FUT1AA genotype showed a higher intensity in the signal of betaine, a small N-trimethylated amino acid that originates from the oxidation of choline, which functions as a methyl donor. Donation of methyl groups is important for an optimal liver function, cellular replication, and detoxification reactions [47]. Furthermore, a decreased betaine level has been associated with stress due to weaning and exposure to E. coli [48]. Our results identified two nitrogenous bases, guanosine and 3-methylguanine level, at weaning (day 28) that were affected by genotype, FUT1AA versus FUT1AG. The higher level of guanosine in FUT1AG piglets probably indicated that these animals could had suffered of an inflammatory process. Indeed, it has been shown that extracellular guanosine tended to increase and accumulate in response to injury as oxygene and glucose deprivation in rats [49,50] and it promotes a glioprotective mechanism against the inflammatory and oxidative damage induced by lipopolysaccharide (LPS) exposure in primary cultures of hippocampal astrocytes by the activation of heme oxygenase-1 (HO-1) pathway [51]. Piglets with the FUT1AA genotype were characterized by a higher 3-methylguanine level, a metabolite involved in the DNA repair process [52]. However, we are not aware of any animal studies available on the role of this compound in plasma, and our work does therefore not allow us to associate this metabolite with animal health or biological processes.

Regarding the post-weaning period, the metabolome profile of plasma samples obtained on day 34 showed that the influence of FUT1 genotype on the plasma metabolome had decreased when compared to day 7 and 28. This result may be ascribed to a range of different factors perturbing the piglet’s physiological homeostasis after weaning. The effect of FUT1 genotype on the metabolic profile discriminations was attributed to compounds such as tauroursodeoxycholic acid, a bile acid compound, and L-phenylalanine, an amino acid precursor. Phenylalanine is an essential amino acid, derived by the protein catabolism and which is mainly converted to tyrosine in the kidneys [53], while tauroursodeoxycholic acid is involved in the digestion and intestinal absorption of hydrophobic nutrients. A higher level of these compounds was found in plasma of in plasma of FUT1AA piglets, and may hence be relevant for their nutrient absorption capability. On the other hand, in our study, FUT1AG piglets have been found to have a higher level of genistein 5-o-glucuronide, a polyphenol metabolite involved in the intestinal barrier regulation [54]. Our results on the expression of candidate genes involved in the intestinal barrier function (ZO-1 and OCLN) and immunity response (IL10, TNF-α and COX-2) showed no differences between FUT1AA and FUT1AA piglets, however, it could be possible that the differences due to FUT1 variants were not sufficient to demonstrate differences in these gene markers that are generally stimulated after infection, indeed our animals stayed healthy until the end of the trial.

In order to assess the microbiota metabolism along the gastrointestinal tract of the two FUT1 groups, organic acid analysis was carried out. Acetic acid, which is the major short-chain fatty acid in the gut, was found in higher concentrations in mid-colon digesta of FUT1AG. As the short-chain fatty acid profile reflects bacterial metabolism [55], the higher colonic concentration of acetic acid in FUT1AG piglets most likely reflects bacterial community differences between genotypes. Bacterial culture showed a higher number of haemolytic bacteria, both in faeces and digesta, and of Enterobacteriaceae in digesta content from FUT1AG piglets. E. coli belongs to this family and has been shown to produce higher amounts of acetic acid at high growth rates [56]. This alone is unlikely to explain the difference in acetic acid concentrations as we have a lack of information about the detailed commensal bacteria present in the intestinal tract of the animals, but all these results together may indicate that FUT1 genotype can influence the microbiota balance, the microbial metabolism and thus the host physiological condition.

We were not able to detect significantly higher E. coli F18 numbers in the small intestinal tissue using an ex vivo porcine intestinal model. Further, we were not able to detect weight gain differences, although some plasma metabolomic differences between genotypes could be ascribed to feed intake. The lack of significance between the genotypes regarding epithelial adherence of E. coli as obtained in the ex vivo challenge model, and the pig body weight responses might be due to the low number of replicates in the present study.

Conclusion

The results obtained in the present study contribute to the existing knowledge on FUT1 polymorphisms effect in pigs immune responses and gastrointestinal bacteria. Our investigation showed that FUT1 genotype influenced the number of haemolytic bacteria and Enterobacteriaceae in the intestinal content. In addition, FUT1 genotype influenced the plasma metabolic profile, especially during the pre-weaning period, and specific plasma metabolites related with gut microbial metabolism, inflammation and nutrient absorption capability. Further study on host FUT1 genotype are desirable to improve the knowledge on the relationship between FUT1 genotype, the mucosal bacterial community, and the host immune response in order to elucidate the perspective of e.g. breeding strategies on gut health parameters others than diarrhoea. Moreover, our work encourages future studies that evaluate the role and contribution of Fucosyltransferase genes family on the mechanisms involved in the pathophysiology of intestinal disease and inflammation.

Supporting information

S1 Table. Complete list of discriminant plasma metabolites for FUT1AG and FUT1AA obtained from piglets at 7, 28 and 34 days of age by the sMLPLS-DA.

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

(DOCX)

S2 Table. Relative expression of selected genes in mucosa from the stomach, ileum, caecum and mid colon of piglets 34 days of age (one-week post-weaning) belong to FUT1AG and FUT1AA groups.

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

(DOCX)

S3 Table. Dry matter content (%) and pH of digesta from the gastrointestinal tract of 34 days old piglets (one-week post-weaning)1.

https://doi.org/10.1371/journal.pone.0202970.s003

(DOCX)

Acknowledgments

This study is based upon work from COST Action FA1401, supported by COST (European Cooperation in Science and Technology). The authors would like to thank laboratory technicians Trine Poulsen and Karin Durup, technicians Mette Lykkegaard and Karin Johansen and the staff of SB32 at the Department of Animal Science for their skilful technical assistance during the course of the experiment. The authors want to thank Scientist Lone Bruhn Madsen, and research professor Christian Bendixen at Department of Molecular Biology and Genetics for FUT1 genotyping the animals.

References

  1. 1. Frank DN, Robertson CE, Hamm CM, Kpadeh Z, Zhang T, Chen H, et al. Disease phenotype and genotype are associated with shifts in intestinal-associated microbiota in inflammatory bowel diseases. Inflamm Bowel Dis. 2011;17:179–184. pmid:20839241
  2. 2. McGovern DPB, Jones MR, Taylor KD, Marciante K, Yan X, Dubinsky M, et al. Fucosyltransferase 2 (FUT2) non-secretor status is associated with Crohn’s disease. Hum Mol Genet. 2010;19:3468–3476. pmid:20570966
  3. 3. Jostins L, Ripke S, Weersma RK, Duerr RH, Dermot P. Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature. 2013;491:119–124.
  4. 4. Bao W Bin, Ye L, Zi C, Su XM, Pan ZY, Zhu J, et al. Study on the age-dependent tissue expression of FUT1 gene in porcine and its relationship to E. coli F18 receptor. Gene. 2012;497:336–9. pmid:22305985
  5. 5. Meijerink E, Neuenschwander S, Fries R, Dinter A, Bertschinger HU, Stranzinger G, et al. A DNA polymorphism influencing (1,2)fucosyltransferase activity of the pig FUT1 enzyme determines susceptibility of small intestinal epithelium to Escherichia coli F18 adhesion. Immunogenetics. 2000;52:129–136. pmid:11132149
  6. 6. Frydendahl K, Jensen TK, Andersen JS, Fredholm M, Evans G. Association between the porcine Escherichia coli F18 receptor genotype and phenotype and susceptibility to colonisation and postweaning diarrhoea caused by E. coli O138:F18. Vet Microbiol. 2003;93:39–51. pmid:12591205
  7. 7. Francis DH. Enterotoxigenic Escherichia coli infection in pigs and its diagnosis. J Swine Heal Prod. 2002;10:171–175.
  8. 8. Vogeli P, Bertschinger HU, Stamm M, Stricker C, Hagger C, Fries R, et al. Genes specifying receptors for F18 fimbriated Escherichia coli, causing oedema disease and postweaning diarrhoea in pigs, map to chromosome 6. Anim Genet. 1996;27:321–328. pmid:8930072
  9. 9. Meijerink E, Fries R, Vogeli P, Masabanda J, Wigger G, Stricker C, et al. Two α (1,2) fucosyltransferase genes on porcine chromosome 6q11 are closely linked to the blood group inhibitor (S) and Escherichia coli F18 receptor (ECF18R) loci. Mamm Genome. 1997;8:736–741. pmid:9321466
  10. 10. Henry S, Mollicone R, Fernandez P, Samuelsson B, Oriol R, Larson G. Molecular basis for erythrocyte Le(a+ b+) and salivary ABH partial-secretor phenotypes: expression of a FUT2 secretor allele with an A—>T mutation at nucleotide 385 correlates with reduced alpha(1,2) fucosyltransferase activity. Glycoconj J. 1996;13:985–993. pmid:8981090
  11. 11. Pickard JM, Maurice CF, Kinnebrew MA, Abt MC, Schenten D, Golovkina TV, et al. Rapid fucosylation of intestinal epithelium sustains host–commensal symbiosis in sickness. Nature. 2014;514:638–641. pmid:25274297
  12. 12. Hollingsworth MA, Swanson BJ. Mucins in cancer: protection and control of the cell surface. Nat Rev Cancer. 2004;4:45–60. pmid:14681689
  13. 13. Hesselager MO, Everest-Dass A V, Thaysen-Andersen M, Bendixen E, Packer NH. FUT1 genetic variants impact protein glycosylation of porcine intestinal mucosa. Glycobiology. 2016;26:607–622. pmid:26858341
  14. 14. Hasson CJ, Caldwell GE, Emmerik REA Van. Gastrointestinal Mucins of Fut2-Null Mice Lack Terminal Fucosylation without Affecting Colonization by Candida albicans. Glycobiology. 2005;15:1002–1007. pmid:15958416
  15. 15. Domino SE, Zhang L, Gillespie PJ, Saunders TL, Lowe JB. Deficiency of Reproductive Tract α (1, 2) Fucosylated Glycans and Normal Fertility in Mice with Targeted Deletions of the FUT1 or FUT2 α (1, 2) Fucosyltransferase Molecular and Cellular Biology. 2001;21:8336–8345.
  16. 16. Nanthakumar NN, Meng D, Newburg DS. Glucocorticoids and microbiota regulate ontogeny of intestinal fucosyltransferase 2 requisite for gut homeostasis. Glycobiology. 2013;23:1131–1141. pmid:23887940
  17. 17. McKnite AM, Perez-Munoz ME, Lu L, Williams EG, Brewer S, Andreux PA, et al. Murine gut microbiota is defined by host genetics and modulates variation of metabolic traits. PLoS One. 2012;7: e39191. pmid:22723961
  18. 18. Sommer F, Backhed F. The gut microbiota—masters of host development and physiology. Nat Rev Micro. 2013;11:227–238.
  19. 19. Schokker D, Veninga G, Vastenhouw SA, Bossers A, de Bree FM, Kaal-Lansbergen LMTE, et al. Early life microbial colonization of the gut and intestinal development differ between genetically divergent broiler lines. BMC Genomics. 2015;16:418. pmid:26017153
  20. 20. Bao WB, Wu SL, Musa HH, Zhu GQ, Chen GH. Genetic variation at the alpha-1-fucosyltransferase (FUT1) gene in Asian wild boar and Chinese and Western commercial pig breeds. J Anim Breed Genet. 2008;125 (6):427–430. pmid:19134079
  21. 21. Sugiharto B, Jensen B, Lauridsen C. Development of an ex vivo model for investigating the bacterial association to the gut epithelium of pigs. J Anim Sci. 2012;90:397–399. pmid:23365391
  22. 22. Ingerslev AK, Karaman I, Bağcıoğlu M, Kohler A, Theil PK, Bach Knudsen KE, et al. Whole Grain Consumption Increases Gastrointestinal Content of Sulfate-Conjugated Oxylipins in Pigs − A Multicompartmental Metabolomics Study. J Proteome Res. 2015;14:3095–3110. pmid:26133160
  23. 23. Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics: Official journal of the Metabolomic Society. 2007;3:211–221.
  24. 24. Sugiharto S, Poulsen A- SR, Canibe N, Lauridsen C. Effect of bovine colostrum feeding in comparison with milk replacer and natural feeding on the immune responses and colonisation of enterotoxigenic Escherichia coli in the intestinal tissue of piglets. Br J Nutr. 2015;113:923–934. pmid:25743486
  25. 25. Canibe N, Miettinen H, Jensen BB. Effect of adding Lactobacillus plantarum or a formic acid containing-product to fermented liquid feed on gastrointestinal ecology and growth performance of piglets. Livest Sci. 2008;114:251–262.
  26. 26. Holdeman LV, Cato EP, Moore EC. Anaerobe laboratory manual. 1977. Virginia Polytechnic Institute and State University: Blacksburg
  27. 27. Miller TL, Wolin MJ. A serum bottle modification of the Hungate technique for cultivating obligate anaerobes. Appl Microbiol. 1974;27:985–987. pmid:4598231
  28. 28. Westerhuis JA, van Velzen EJJ, Hoefsloot HCJ, Smilde AK. Multivariate paired data analysis: Multilevel PLSDA versus OPLSDA. Metabolomics. 2010;6:119–128. pmid:20339442
  29. 29. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. arXiv preprint arXiv 2014: 1406.5823.
  30. 30. Chun H, Keleş S. Sparse partial least squares regression for simultaneous dimension reduction and variable selection. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2010 (72): 3–25. pmid:20107611
  31. 31. Halekoh U, Højsgaard S. Kenward-Roger A. Approximation and Parametric Bootstrap Methods for Tests in Linear Mixed Models—The R Package pbkrtest. J Stat Softw. 2014;59:1–32.
  32. 32. Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008;50:346–363. pmid:18481363
  33. 33. Suzuki T. Regulation of intestinal epithelial permeability by tight junctions. Cell. Mol. Life Sci. 20313; 70, 631–659. pmid:22782113
  34. 34. Plaza-Diaz J, Gomez-Llorente C, Fontana L, Gil A. Modulation of immunity and inflammatory gene expression in the gut, in inflammatory diseases of the gut and in the liver by probiotics. World Journal of Gastroenterology: WJG. 2014;20(42):15632–15649. pmid:25400447
  35. 35. Lallès J- P, Bosi P, Smidt H, Stokes CR. Nutritional management of gut health in pigs around weaning. Proc Nutr Soc. 2007;66:260–268. pmid:17466106
  36. 36. Stokes CR, Bailey M, Haverson K, Harris C, Jones P, Inman C, et al. Postnatal development of intestinal immune system in piglets: implications for the process of weaning. Anim Res. 2004;53:325–334.
  37. 37. Hou C, Zeng X, Yang F, Liu H, Qiao S. Study and use of the probiotic Lactobacillus reuteri in pigs: a review. J Anim Sci Biotechnol. 2015;6:1–8.
  38. 38. Daly K, Darby AC, Hall N, Nau A, Bravo D, Shirazi-Beechey SP. Dietary supplementation with lactose or artificial sweetener enhances swine gut Lactobacillus population abundance. Br J Nutr. 2014;111:30–35. pmid:24382146
  39. 39. Bian G, Ma S, Zhu Z, Su Y, Zoetendal EG, Mackie R, et al. Age, introduction of solid feed and weaning are more important determinants of gut bacterial succession in piglets than breed and nursing mother as revealed by a reciprocal cross-fostering model. Environ Microbiol. 2016;18:1566–1577. pmid:26940746
  40. 40. Chen L, Xu Y, Chen X, Fang C, Zhao L, Chen F. The Maturing Development of Gut Microbiota in Commercial Piglets during the Weaning Transition. Front. Microbiol. 2017;8:1–13.
  41. 41. Mach N, Berri M, Estellé J, Levenez F, Lemonnier G, Denis C, et al. Early-life establishment of the swine gut microbiome and impact on host phenotypes. Environ Microbiol Rep. 2015;7:554–569. pmid:25727666
  42. 42. Slifierz MJ, Friendship RM, Weese JS. Longitudinal study of the early-life fecal and nasal microbiotasof the domestic pig. Bmc Microbiol. 2015;15:12.
  43. 43. Poon K, Pang KS. Benzoic acid glycine conjugation in the isolated perfused rat kidney. Drug Metab Dispos. 1995;23:255–260. pmid:7736921
  44. 44. Armstrong MD, Chao FC, Parker VJ, Wall PE. Endogenous formation of hippuric acid. Proc Soc Exp Biol Med. 1955;90:675–679. pmid:13289873
  45. 45. Wikoff WR, Anfora AT, Liu J, Schultz PG, Scott A, Peters EC, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Biological Sciences—Biochemistry: PNAS 2009 106 (10) 3698–3703. pmid:19234110
  46. 46. Williams HR, Cox IJ, Walker DG, Cobbold JF, Taylor-Robinson SD, Marshall SE, et al. Differences in gut microbial metabolism are responsible for reduced hippurate synthesis in Crohn’s disease. BMC Gastroenterol. 2010;10:108. pmid:20849615
  47. 47. Wang J, Wu Z, Li D, Li N, Dindot S V, Satterfield MC, et al. Nutrition, epigenetics, and metabolic syndrome. Antioxidants Redox Signal. 2012;17:282–301. pmid:22044276
  48. 48. Sugiharto S, Hedemann MS, Lauridsen C. Plasma metabolomic profiles and immune responses of piglets after weaning and challenge with E. coli. J Anim Sci Biotechnol. 2014;5:17. pmid:24624922
  49. 49. Ciccarelli R, Di Iorio P, Giuliani P, D’Alimonte I, Ballerini P, Caciagli F, et al. Rat cultured astrocytes release guanine-based purines in basal conditions and after hypoxia/hypoglycemia. Glia. 1999;25:93–98. pmid:9888301
  50. 50. Dal-Cim T, Ludka FK, Martins WC, Reginato C, Parada E, Egea J, et al. Guanosine controls inflammatory pathways to afford neuroprotection of hippocampal slices under oxygen and glucose deprivation conditions. J Neurochem. 2013;126:437–450. pmid:23713463
  51. 51. Bellaver B, Souza DG, Bobermin LD, Goncalves CA, Souza DO, Quincozes-Santos A. Guanosine inhibits LPS-induced pro-inflammatory response and oxidative stress in hippocampal astrocytes through the heme oxygenase-1 pathway. Purinergic Signal. 2015;11:571–580. pmid:26431832
  52. 52. Curtis CD, Thorngren DL, Nardulli AM. Immunohistochemical analysis of oxidative stress and DNA repair proteins in normal mammary and breast cancer tissues. BMC Cancer. 2010;10:9. pmid:20064251
  53. 53. Bröer S, Bröer A. Amino acid homeostasis and signalling in mammalian cells and organisms. Biochemical Journal. 2017;474(12):1935–1963. pmid:28546457
  54. 54. Suzuki T. Regulation of intestinal epithelial permeability by tight junctions. Cell Mol Life Sci. 2013;70:631–659. pmid:22782113
  55. 55. Williams BA, Bosch MW, Awati A, Konstantinov SR, Smidt H, Akkermans ADL, et al. The effects of prestorage incubation and length of storage of broiler breeder eggs on hatchability and subsequent growth performance of progeny. Anim Res. 2005;54:191–201.
  56. 56. Han K, Lim HC, Hong J. Acetic acid formation in escherichia coli fermentation. Biotechnol Bioeng. 1992;39:663–671. pmid:18600996