Transcriptional changes detected in fecal RNA of neonatal dairy calves undergoing a mild diarrhea are associated with inflammatory biomarkers

After birth, a newborn calf has to adapt to an extrauterine life characterized by several physiological changes. In particular, maturation of the gastrointestinal tract in a new environment loaded with potential pathogens, which can predispose neonatal calves to develop diarrhea, and is a major cause of morbidity and mortality during the first 4 wks of life. We aimed to investigate the inflammatory adaptations at a transcriptomic level in the gastrointestinal (GI) tract to a mild diarrhea in neonatal dairy calves using RNA isolated from fresh fecal samples. Eight newborn Jersey male calves were used from birth to 5 wks of age and housed in individual pens. After birth, calves received 1.9 L of colostrum from their respective dams. Calves had ad-libitum access to water and starter grain (22% CP) and were fed twice daily a total of 5.6 L pasteurized whole milk. Starter intake, body weight (BW), fecal score, withers height (WH), and rectal temperature (RT) were recorded throughout the experiment. Blood samples were collected weekly for metabolic and inflammatory profiling from wk 0 to wk 5. Fresh fecal samples were collected weekly and immediately flash frozen until RNA was extracted using a Trizol-based method, and subsequently, an RT-qPCR analysis was performed. Orthogonal contrasts were used to evaluate linear or quadratic effects over time. Starter intake, BW, and WH increased over time. Fecal score was greatest (2.6 ± 0.3) during wk 2. The concentrations of IL-6, ceruloplasmin, and haptoglobin had a positive quadratic effect with maximal concentrations during wk 2, which corresponded to the maximal fecal score observed during the same time. The concentration of serum amyloid A decreased over time. The mRNA expression of the proinflammatory related genes TLR4, TNFA, IL8, and IL1B had a positive quadratic effect of time. A time effect was observed for the cell membrane sodium-dependent glucose transporter SLC5A1, for the major carbohydrate facilitated transporter SLC2A2, and water transport function AQP3, where SLC5A1 and AQP3 had a negative quadratic effect over time. Our data support the use of the fecal RNA as a noninvasive tool to investigate intestinal transcriptomic profiling of dairy calves experiencing diarrhea, which would be advantageous for future research including nutritional effects and health conditions.


Introduction
After birth, the gastrointestinal tract (GI) of a newborn calf undergoes morphological and functional adaptations, which starts with the ingestion of colostrum affecting the gastrointestinal motility as well as the acquisition of the passive immunity [1,2].In addition to the gastrointestinal maturation, the calf has to adapt to a new environment, which is loaded with pathogens, thermoregulate its own body and transition from pre-ruminant to ruminant during the weaning period.Together, these factors can prone calves to develop diarrhea, which is the leading cause of mortality in neonatal dairy calves in the US during the first 4 weeks of life [3].Fecal score is commonly used by researchers and dairy producers to keep track and identify young animals with diarrhea.Also, more sophisticated methods such as microscopy and bacterial culturing have been used to identifying specific strains of infectious pathogens causing diarrhea.Blood biomarkers of inflammation, oxidative stress, metabolism are useful indicators of overall animal health.However, these analyses along with fecal score do not provide additional information on the actual biological adaptations the GI tract undergoes during the onset of diarrhea.Epithelial cells from the GI tract act as a barrier against pathogens, and the daily shedding of these cells into feces has been used in rodents and humans as a source of RNA to evaluate gene expression in the GI tract [4].The exfoliation of these cells has been associated with intestinal natural turnover process to preserve tissue structure [5].This is a noninvasive method that can be used instead of performing tissue biopsies and postmortem tissue analysis [6].Originally, RNA from fecal samples was used for detecting molecular biomarkers of colon cancer in rats [7,8] and later for intestinal development and function in adult humans [9,10].The basis for using fecal RNA is that, unlike mammalian RNA, bacteria RNA lacks a polyA tail making possible to discriminate between RNA from the host and RNA from bacteria by targeting RNA with the polyA.Certain strains of fungi, which are also considered pathogens, can produce polyA tail [11,12], however, the extent and influence of these strains in the fecal RNA remain unknown.In humans, the RNA isolated from stool samples were used to study nutritional regimens [13], as well as, to quantify intestinal gene expression profiles in developing human infants [14].In fact, in the latter study, using fecal RNA, Chapkin and collaborators [14] were able to differentiate between neonates that were breastfed and formula-fed by their different transcriptome profile.These data underscore the potential of the fecal RNA method to study biological alterations to the GI tract product of nutritional effects.Similarly, such method could be applied to study the effects of inflammatory conditions arise from diarrhea in neonatal dairy calves.
Therefore, we hypothesized that inflammatory signals could be detected in gene expression of fecal RNA in response to adaptations prior to-and during the onset of diarrhea in neonatal calves.Thus, the objective of this study was to use RNA isolated from fecal samples to study the molecular adaptations to a mild-diarrhea via assessment of expression of specific genes involved in inflammation and cell membrane transporters through RT-qPCR and validate the findings using well-established profiling of inflammatory biomarkers in blood.

Experimental design, animal management
The Institutional Animal Care and Use Committee (IACUC) of the Oregon State University (OSU) approved all procedures for this study (Protocol # 4747).Eight healthy newborn male Jersey calves from the OSU Dairy Center were used in the study by applying the following inclusion criteria: 1) calving difficulty score < 3, 2) single calf, and 3) calf birth weight !34 kg.After birth, calves were supplemented with BO-SE (Merck Animal Health, Germany) and vaccinated with CALF-GUARD 1 (Bovine Rota-Coronavirus vaccine, Zoetis, MI), and received 1.9 L of colostrum from the respective dam in the first two feedings.Calves remained in the study from birth to 5 wk of age (i.e., preweaned period).Calves were housed in individual pens bedded with straw at OSU Dairy Center and fed twice daily with pasteurized whole cow milk.Calves had ad libitum access to water and a starter grain (Ampli-Calf 1 ; 22% CP, Purina Animal Nutrition, Shoreview, MN) throughout the experiment.
Growth performance including BW and withers height (WH) were recorded weekly before morning feeding.Health evaluations including fecal score [scale 1-4, 1: firm, well formed (not hard); 2: soft, pudding-like; 3: runny, pancake batter; 4: liquid, splatters; [15]] and respiratory score [scale 1-5, 1: normal 2: runny nose, 3: heavy breathing, 4: cough moist, 5: cough dry; [15]] were recorded daily throughout the experiment.Based on the above fecal score scale, we classified calves with a mild-diarrhea with fecal score of 3, and diarrhea with fecal score of 4. Rectal temperature was recorded daily until 21 d of age.Fresh fecal samples (~10 g) were collected from the rectum weekly ~30 min after morning feeding from wk 1 to 5 (or 7 to 35 d).Fecal samples were immediately flash-frozen in liquid N 2 and kept at -80˚C until further analysis.

Blood sample collection
Blood samples were collected prior to morning feeding from the jugular vein using a 20-gauge BD vacutainer needles (Becton Dickinson, Franklin Lakes, NJ).Samples were collected into evacuated tubes (10 mL, BD Vacutainer 1 , Becton Dickinson, Franklin Lakes, NJ) containing either serum clot activator or sodium heparin within 30 minutes after the morning feeding starting from the 1 d after calving (i.e., 0 wk) and subsequent samples were collected weekly until wk 5.After blood collection, tubes with sodium heparin were placed on ice (4˚C), and tubes with clot activator were kept at 21˚C (~30 min) until centrifugation.Serum and plasma were obtained by centrifugation of clot activator and sodium heparin tubes, respectively, at 1,200 × g for 15 minutes.Aliquots of serum and plasma were frozen (-80˚C) until further analysis.

RNA isolation, cDNA synthesis, primer design and evaluation and quantitative PCR
Complete details regarding the RNA isolation and cDNA synthesis are presented in S1 File.The RNA isolation from fecal samples was carried using the RNeasy Plus Mini Kit (Qiagen, Cat.No. 74134), following the manufacturer's instructions with some modifications.Briefly, prior to the RNA isolation, the fecal samples (~200 mg) were homogenized using a Bullet Blender Next Advance (Laboratory Instruments, USA) immersed in 1.2 mL of TRIzol 1 Reagent (Cat# 15596018, Ambion, Carlsbad, CA, USA).After homogenization, 240 μl of chloroform (EMD Millipore, Germany, Cat.No. CX1055-6) was added in order to isolate the RNA from the organic phase.The RNA quantity (822.9 ± 491.8 ng/μL; mean ± SD) and purity as 260/280 ratio (1.98 ± 0.3) were determined via spectrophotometry using SpectraMax 1 Plus 384 (Molecular Devices, Sunnyvale, CA, USA).The RNA integrity was assessed using the TapeStation 1 (Agilent Technologies), and the final RIN number was 3.4 ± 1.5 with a range from 1.8 to 7.4.This RIN value can be considered as a limitation for this study since is lower than desirable values (! 5) for RT-qPCR.However, in this study, this effect did not translate in an adverse effect on the PCR amplification curves across samples.Future optimizations to this technique will improve this parameter.

Statistical analysis
Data were analyzed using the PROC MIXED procedure of SAS 9.4 (SAS Institute, Inc., Cary, NC, USA).The model included time as fixed effect and calf as random.The autoregressive (1) covariance structure was used for repeated measures for all parameters analyzed.Blood metabolites and gene expression results were log-scale transformed if needed, in order to comply with a normal distribution of residuals.Daily fecal scores and starter intakes were average weekly before statistical analysis.Orthogonal contrasts were used to determine linear and quadratic effects when a trend (P < 0.10) for a significant time effect was observed.Statistical significance and tendencies were declared at P 0.05 and 0.05 P 0.10, respectively.

Growth performance and health
Main effect of time for performance parameters are shown in Fig 1 .There was a time effect (P < 0.01) observed on starter intake, which increased over time.Similar to starter intake, BW and WH increased (P < 0.01) over time.Fecal score was affected by time (P < 0.001), where a maximal fecal score (2.6 ± 0.3) was observed during wk 2, and returned to initial values by wk 5. Based on our classification of mild-diarrhea or fecal score 3 (i.e., runny or package batter), 7 out of the 8 calves experienced at least a mild-diarrhea for at least 1 day.Similarly, 4 out of the 8 calves experience diarrhea or fecal score 4 for at least 1 day.The descriptive distribution of daily fecal scores per week presented in S1 Fig .indicates a total occurrence of ~40% and 30% for at least a mild-diarrhea events, respectively, during wk 2. The morbidity based on fecal scores was 4.6 ± 2.5 and 5.8 ± 4.3 d, for mild-diarrhea and diarrhea, respectively.Consequently, morbidity increased to 7.9 ± 5.9 d, when evaluating numbers of days calves experience at least a mild-diarrhea (i.e., fecal score !3).The maximal fecal score during wk 2 was reflected in a morbidity of 4.0 ± 3.5 d based on calves having at least a mild-diarrhea.

Blood biomarkers
Main effect of time on energy and muscle mass catabolism metabolites are presented in Fig 2 .Energy metabolites.There was a time effect (P = 0.02) observed for BHBA and a trend (P = 0.06) for glucose.A linear (P 0.04) response was observed for glucose and BHBA,  quadratic effect (P 0.04) was observed in ceruloplasmin, haptoglobin, and IL-6.The trend (P = 0.10) for a time effect on IL-1β was not explained by a linear (P = 0.21) or quadratic effect (P = 0.90), but rather by an overall trend (P 0.08) for an increase at wk 1 in comparison with wk 0, 2, 3, and 4 (Fig 4).
Oxidative stress.Main effect of time on inflammatory biomarkers are presented in Fig 5 .There was a main effect (P < 0.01) on ROM and FRAP.There was a linear (P 0.01) and quadratic (P< 0.01) response on ROM and FRAP.A positive quadratic effect was observed for ROM, and a negative quadratic effect was observed for FRAP (Fig 5).

Gene expression
Immune receptors and proinflammatory cytokines.Main effects for immune receptors related genes are shown in Fig 6 .There was a time effect (P <0.01) on TLR4, TNFA, IL1B, IL8, and IFNG and a trend on NFKB1(P = 0.07) over time.There was a linear response (P< 0.04) for all genes evaluated.A positive quadratic effect (P = 0.05) was observed for TLR4, TNFA, IL1B, and IL8.No quadratic effect (P !0.64) was observed for TLR2, NFKB1 or IFNG, however, the expression of TLR2 and IFNG decreased (P 0.03) over time while NFKB1 increased (P = 0.04).There was a time effect (P < 0.04) observed for SLC5A1, AQP3, and SLC2A2.A linear response was observed (P < 0.01) for SLC2A2 and with no linear response for AQP3 and SLC5A1.A negative quadratic effect (P 0.05) was observed on AQP3 and SLC5A1.

Animal performance
Although there was an evident increase in scours during 2 and 3 wk (Fig 1B), overall growth was not compromised throughout the study.The starter intake at the end of the trial was 0.193 kg/d, which is in agreement with starter intakes (0.2 kg/d starter intake) reported by others [24].Although ADG reached a nadir level at wk 2, it increased linearly until wk 5 (Fig 1E ), and, in fact, the overall ADG (0.42 vs 0.21 kg/d) was greater than that reported in Jersey calves across treatments in [25].

Blood biomarkers for energy and protein metabolism during early life
The concentration of BHB was likely increased linearly by the increased starter intake due to the stimulation of the ruminal fermentation, which is an indicator of rumen development of young ruminants [26].The BHB levels in this study ranged from 0.12 mmol/L (wk 0) to 0.20 mmol/L (wk 5) which are comparable to those observed by Quigley et al. [26] and Klotz and Heitmann [25], where BHB concentrations increased in neonatal Jersey calves regardless of dietary treatments.The increased concentration of glucose throughout this trial can be associated with a physiological response, where neonates have to develop their gluconeogenesis pathway to generate endogenous energy to reach the normal metabolism [27].In our study, glucose concentration ranged from 6.8 mmol/L at wk 0 to 7.5 mmol/L at wk 5, and was higher than other studies [25,28].
In this study, the blood creatinine was in agreement with Khan et al. [29], reporting creatinine values from 80.4 (wk 2) to 67.2 μmol/L (wk 12) and in both experiments creatinine levels decreased with age.The latter suggests that creatinine levels at birth might be associated with a transfer from maternal creatinine as observed in humans [30].The overall increase in blood urea in our study is associated with the development of the ruminal fermentation, which has been demonstrated that blood urea concentrations in dairy calves increase over time due to rumen development stimulated by solid feed, especially after weaning [31].Also, plasma urea concentration increased in neonatal dairy calves, when calves began to consume greater amounts of grain [32,33] which result in greater production of ruminal ammonia, and consequently more metabolized urea into the bloodstream.

Biomarkers of liver function
Liver function is commonly depressed during inflammatory episodes [17], which results in reduced albumin, PON, and cholesterol in the blood.In contrast, greater blood GGT and GOT are associated with liver cell damage [34].Therefore, in the current study, the linear increase of cholesterol, albumin, and PON in the blood of calves does not suggest a liver dysfunction even during 2 wk of age when the fecal score was maximal.In addition, the clearance of total bilirubin in blood is associated with a progressive increase in liver function (Fig 3).The high GGT concentration at wk 0 can be associated with colostrum intake, as this enzyme is present in high levels in colostrum (Gungor et al., 2004).Then, the decreased GGT concentration over time is in accordance with other studies in neonatal dairy calves [35,36].

Biomarkers of oxidative stress
Oxidative stress is driven by the imbalance between the production of reactive oxygen metabolites (ROM) and antioxidants.The process of phagocytosis by which leukocytes such as neutrophils and macrophages can destroy invading pathogens can also lead to increased production of ROM [37].The fact that ROM reached a maximum level in blood at 2 wk of age (19 mg H 2 O 2 /100 mL), together with the decrease in FRAP concentration (antioxidant parameter; Fig 5) during the same time which calves experienced an increase in the fecal score.Taken these data together, suggest that this is a physiological response during a period of stress such as diarrhea.Therefore, increased ROM could be associated with the increase in diarrhea events leading to a potential increase phagocytosis, arachidonic acid metabolism, oxidation of tissue reserves, among other mechanisms to counteract this condition [38].

Inflammatory response in fecal RNA and blood biomarkers
Toll-like receptors (TLRs) present in the cell membrane can identify an invading pathogen, and start an intracellular signaling pathway with the purpose to activate the transcription factor NFKB1 which promotes the transcription of genes encoding proinflammatory proteins (e.g., cytokines and chemokines), such as IL1B, TNFA, IL6 and IL8 [39].Among TLRs, TLR2 and TLR4 are perhaps the most widely studied, which recognize components of the bacterial cell walls such as lipoteichoic acids (LTA) and lipopolysaccharides (LPS) on Gram-positive and Gram-negative bacteria, respectively [40].Therefore, it is plausible that the upregulation of TLR4 mRNA expression observed in fecal samples was triggered by the presence of gramnegative bacteria (i.e., E. coli).Similar responses have been observed in other tissues such as the mammary gland where an upregulation of TLR4 was observed in dairy cows experiencing an inflammatory condition due to LPS [41,42] or E. coli challenge [43].Among the limited data on TLR4 expression in GI tract in dairy calves, Malmuthuge et al. [44] did not observe differences in TLR4 expression across several GI tract tissues (i.e., rumen, jejunum, ileum, cecum, and colon) in 3 wk-old calves.Although an upregulation of TLR4 was observed with a maximal expression at 4 wk, this does not precisely overlap with the increase in fecal scores at 1 wk.However, the latter could be associated with an alternative pathogen and consequently mechanism, other than TLR4 pathway, responsible for triggering the increased fecal scores at 1 wk.Certainly, there still much to know regarding the GI tract of neonatal calves, for instance, Ontsouka et al. [45] described the importance in the maturation of the GI tract in neonatal calves and how bioactive substance in colostrum can affect this maturation.Therefore, the latter could be another potential explanation for the delay observed in TLR4 with respect of the increased fecal scores in 1 wk observed in this study.
The concomitant downregulation over time observed in this study for TLR2 and IFNG can be associated with early life exposure to Mycobacterium species in colostrum [46], for which TLR2 has been reported as the main TLR to recognized these pathogens [47], and IFNγ has a significant role during the activation of TLR2 signaling pathway [48].
Because of the downregulation of TLR2, it is likely that most of the inflammatory response observed in fecal RNA was originated through the upregulation of TLR4 at 4 wk.For instance, the simultaneous upregulation of NFKB1 in fecal RNA at 4 wk and its target genes TNFA, IL1B, and IL8 can be associated with an overall activation of the TLR4 pathway.A similar effect was observed by Santos et al. [49], where ileal sections of calves were inoculated with Salmonella typhimurium, and the overall results demonstrated an upregulation of IL8 and IL1B.
Serum cytokines such as IL-6 have been used as a pro-inflammatory [50,51] or anti-inflammatory [52,53] biomarker, and it is one of the main cytokines responsible for initiating the hepatic acute phase response (APR) [54].During inflammation, IL-6 can be synthesized by the intestinal mucosa and enterocyte by which stimulates enterocytes production of +APP [55,56].Upon an inflammatory condition, the APR increases hepatic synthesis of positive acute phase proteins (+APP) such as SAA, ceruloplasmin, and haptoglobin, and the concomitant decrease in the production of negative acute phase proteins (-APP), such as albumin and apolipoproteins [34,57].The results of this study showed a simultaneous increase in IL-6 in blood and fecal scores, which indicates that calves underwent an inflammatory response.In fact, Fischer et al. [58] reported that serum IL-6 is a reliable biomarker for the inflammatory response to diarrhea in calves.The latter was further confirmed by increased concentrations of +APP such as ceruloplasmin, haptoglobin, and SAA which has been observed previously in neonatal dairy calves infected with Salmonella [56,59,60] and inoculated with bovine viral diarrhea virus [61].As expected inoculation with Salmonella causes diarrhea resulting in a fecal score of 2.5 [60], which is in agreement with our study.Overall, in our study, the pattern of fecal score and pro-inflammatory blood biomarkers (peaking both between 2 and 3 wk of age) indicate a pro-inflammatory response during diarrhea in neonatal dairy calves.The inflammatory response observed in fecal RNA did not overlap with that observed in fecal scores and blood biomarkers at 2 and 3 wk.However, the fact that inflammation-related genes in fecal RNA followed a similar pattern as fecal scores and blood biomarkers of inflammation provide support for fecal RNA being a promising method to study intestinal response to inflammation and its variation over time in neonatal dairy calves.Future improvements to the fecal RNA method will consolidate its reliability to study gut health and response to inflammation, dietary conditions, and growth among others.

Cell membrane transporters in fecal RNA
The absorption of glucose and galactose from the intestinal lumen into the enterocytes is regulated by the sodium-dependent glucose transporter SGLT1 (gene symbol SLC5A1) present in the apical membrane [62].After glucose, galactose, and fructose are transported into the enterocytes, these are further translocated out of the cell into the portal venous system through the facilitated transporter GLUT2 present in the basolateral membrane [63].Previous work in colonic human cells and rabbits showed a decrease in the mRNA expression of SLC5A1 in response to an inflammatory condition induced by synthetic TNFα [64,65].Similarly, the downregulation of SLC5A1 in fecal RNA can be partially related to the time of increased inflammatory response observed in blood biomarkers and fecal score.Thus, the decreased expression of SLC5A1 around the same time of the diarrhea could be suggestive of a physiological response to an inflammatory condition.
As mentioned above GLUT2 is a major transporter of carbohydrates and is encoded by SLC2A2, which was detected in gut segments (i.e., duodenum and jejunum) of 4 d old dairy calves [66].The hepatic mRNA expression of SLC2A2 was downregulated overtime in transition calves from preruminant to ruminants [67], indicating a decreased in the flux of glucose in the portal vein.Therefore, the downregulation overtime in SLC2A2 expression in the fecal samples can be associated with the transition from preruminant to ruminant in those calves, where SLC2A2 expression decreased as the rumen functionality increased.
Another family of cell membrane transporters is the aquaporins, and the Aquaporin-3 (encoded by AQP3) transports water, glycerol, urea, and other small uncharged solutes across the plasma membrane [68,69].Previous studies have shown that AQP3 is down-regulated during inflammatory bowel disease in humans [70] and rats [71].Indeed, it was reported that during bacterial-induced diarrhea in mice the flux of water through the enterocytes is altered compared to a healthy gut environment, where water flux in bacterial-induced diarrhea is blocked from entering or exiting the lumen of the intestine into the epithelium, mainly due to the mislocalization of Aquaporin-3 within the enterocytes [72].Thus, the reduced AQP3 expression observed in the fecal RNA during diarrhea in the present study might be an indication of restricted water exchange between the intestinal lumen and epithelium as a physiological response to diarrhea.

Conclusions
The greater concentration of biomarkers of inflammation in the blood coupled with the upregulation of inflammatory-related genes in fecal RNA from neonatal dairy calves undergoing a mild diarrhea is indicative of an inflammatory response.Such a response confirmed our hypothesis, where physiological changes such as inflammation in response to diarrhea could be detected through transcriptional alterations in fresh fecal samples from neonatal dairy calves.Although further optimization to this method is required in the future, results from this study suggest that RNA isolated from fecal samples has a potential to be a reliable tool to evaluate physiological changes not only due to diarrhea but also to other factors such as dietary treatments.

Fig 1 .
Fig 1. Performance parameters during the neonatal period of Jersey calves.Starter intake (A), fecal scores (B), withers height (C), body weight (D), and average daily gain (ADG) (E).The P values for main effect of time are shown in each plot.The time effect (P< 0.10) was further analyzed in starter intake and fecal scores through orthogonal contrast in order to determine linear (L) and quadratic (Q) effects over time.Values are means, with standard errors represented by vertical bars.https://doi.org/10.1371/journal.pone.0191599.g001

Fig 2 .
Fig 2. Profiles of energy (Glucose and BHBA) and muscle mass catabolism (Creatinine and urea) metabolites in blood during the neonatal period in Jersey dairy calves.The time effect (P< 0.10) was further analyzed through orthogonal contrast in order to determine linear (L) and quadratic (Q) effects over time.Values are means, with standard errors represented by vertical bars.https://doi.org/10.1371/journal.pone.0191599.g002

Fig 4 .
Fig 4. Profiles of inflammatory response (haptoglobin, ceruloplasmin, IL-6, SAA, and IL-1β) biomarkers in blood during the neonatal period in Jersey dairy calves.The time effect (P< 0.10) were further analyzed through orthogonal contrast in order to determine linear (L) and quadratic (Q) effects over time.An overall trend for an increase (P 0.08) in IL-1β at wk 1 in comparison to wk 0, 2, 3, and 4 was denoted by an asterisk ( Ã ).Values are means, with standard errors represented by vertical bars.https://doi.org/10.1371/journal.pone.0191599.g004

Fig 5 .
Fig 5. Profiles of oxidative stress [reactive oxygen metabolites (ROM) and ferric reducing antioxidant power (FRAP)] biomarkers in blood during the neonatal period in Jersey calves.The time effect (P < 0.10) were further analyzed through orthogonal contrast in order to determine linear (L) and quadratic (Q) effects over time.Values are means, with standard errors represented by vertical bars.https://doi.org/10.1371/journal.pone.0191599.g005

Fig 6 .Fig 7 .
Fig 6.Expression of genes related to the activation of the innate immune system (TLR2 and TLR4) and proinflammatory response (NFKB1, TNFA, IL1B, IL8, and IFNG) in fecal samples collected during the neonatal period of Jersey calves.The time effect (P< 0.10) were further analyzed through orthogonal contrast in order to determine linear (L) and quadratic (Q) effects over time.Values are means, with standard errors represented by vertical bars.https://doi.org/10.1371/journal.pone.0191599.g006

Table 1 . Target genes and their biological function. Symbol Name Biological function and process Source Inflammatory
response TLR2 Toll-like receptor 4 Participates in the innate immune response to gram positive microbial agents [16] TLR4 Toll-like receptor 2 Participates in the innate immune response to gram negative microbial agents [16] NFKB1 Nuclear factor of kappa light polypeptide gene enhancer in B-cells 1It is involved in many biological processes such as inflammation, immunity, differentiation, cell growth, tumorigenesis, and apoptosis