- Research article
- Open Access
Unique adaptations in neonatal hepatic transcriptome, nutrient signaling, and one-carbon metabolism in response to feeding ethyl cellulose rumen-protected methionine during late-gestation in Holstein cows
BMC Genomics volume 22, Article number: 280 (2021)
Methionine (Met) supply during late-pregnancy enhances fetal development in utero and leads to greater rates of growth during the neonatal period. Due to its central role in coordinating nutrient and one-carbon metabolism along with immune responses of the newborn, the liver could be a key target of the programming effects induced by dietary methyl donors such as Met. To address this hypothesis, liver biopsies from 4-day old calves (n = 6/group) born to Holstein cows fed a control or the control plus ethyl-cellulose rumen-protected Met for the last 28 days prepartum were used for DNA methylation, transcriptome, metabolome, proteome, and one-carbon metabolism enzyme activities.
Although greater withers and hip height at birth in Met calves indicated better development in utero, there were no differences in plasma systemic physiological indicators. RNA-seq along with bioinformatics and transcription factor regulator analyses revealed broad alterations in ‘Glucose metabolism’, ‘Lipid metabolism, ‘Glutathione’, and ‘Immune System’ metabolism due to enhanced maternal Met supply. Greater insulin sensitivity assessed via proteomics, and efficiency of transsulfuration pathway activity suggested beneficial effects on nutrient metabolism and metabolic-related stress. Maternal Met supply contributed to greater phosphatidylcholine synthesis in calf liver, with a role in very low density lipoprotein secretion as a mechanism to balance metabolic fates of fatty acids arising from the diet or adipose-depot lipolysis. Despite a lack of effect on hepatic amino acid (AA) transport, a reduction in metabolism of essential AA within the liver indicated an AA ‘sparing effect’ induced by maternal Met.
Despite greater global DNA methylation, maternal Met supply resulted in distinct alterations of hepatic transcriptome, proteome, and metabolome profiles after birth. Data underscored an effect on maintenance of calf hepatic Met homeostasis, glutathione, phosphatidylcholine and taurine synthesis along with greater efficiency of nutrient metabolism and immune responses. Transcription regulators such as FOXO1, PPARG, E2F1, and CREB1 appeared central in the coordination of effects induced by maternal Met. Overall, maternal Met supply induced better immunometabolic status of the newborn liver, conferring the calf a physiologic advantage during a period of metabolic stress and suboptimal immunocompetence.
Maternal nutrient and metabolic stresses during pregnancy are important factors that can affect fetal and neonatal growth, development [1, 2], as well as metabolic and inflammatory responses of the offspring . In dairy cattle, the last two months of gestation are the most-critical period for calf fetal growth, since most of the muscle and adipose tissue formation takes place primarily during this time. Around parturition, due to normal decreases in feed intake, cows are exposed to negative energy and essential amino acid (AA) balance, hence, AA such as methionine (Met) become limiting for both cow, fetus, and/or new-born calf [4, 5].
In addition to being a substrate for protein synthesis, Met is a key component of one-carbon metabolism  where it is initially converted into S-adenosylmethionine (SAM), the major biological methyl donor  and where it is involved, through the transsulfuration pathway, in the synthesis of antioxidants glutathione (GSH) and taurine [8, 9]. Physiologically, animals not only obtain Met from the diet, but also from protein breakdown and re-methylation of homocysteine in the Met cycle . Through synthesis of SAM, methyl donors may alter gene transcription in the offspring by methylating DNA and RNA [11, 12]. In non-ruminants, it was demonstrated that maternal methyl donor supplementation (i.e. betaine) led to epigenetic changes that increased expression of genes controlling hepatic gluconeogenesis in the neonate [13, 14].
Epigenetic control of gene expression, also induced by the intrauterine environment, is one of the underlying mechanisms of the ‘Fetal programming’ hypothesis [15, 16] that was first proposed by Barker in 1998 . This concept seeks to explain the effect of maternal nutrition on long-term offspring growth and health [18, 19]; despite its importance, few studies in dairy cattle have addressed the role of nutrient manipulation during late-gestation on fetal and postnatal development. In this regard, for example, Met supplementation is long recognized in dairy cattle as an effective approach to improve productive performance [20, 21], but only recently a promising path in research related to the maternal effect of Met supply on calf health, immune function, and reproductive performance has been highlighted . In particular, it was recently demonstrated that rumen-protected Met (RPM) supplementation during the periparturient period enhanced dry matter intake (DMI) leading to reduced incidence of metabolic disorders and improved overall cow health [23, 24]. Furthermore, enhancing Met supply during late-pregnancy upregulated mRNA abundance of AA and glucose transporters in cow placenta , and was also associated with changes in hepatic one-carbon metabolism and transsulfuration in calf liver . Although the greater DM intake during the last 2–3 wk. prior to parturition that has been consistently reported in cows fed RPM could explain a portion of the greater body mass of the calves at birth [25, 26], other mechanisms potentially encompassing nutrient assimilation efficiency likely play a role.
There are strong associations between Met supplementation during late-pregnancy and body weight and immune response in calves , confirming evidence that AA can affect regulation of metabolic pathways to sustain the immune response against pathogens . The potential role of methyl donors in the early-life innate immune response was recently reported in calves born to cows with high body condition score and after ex vivo lipopolysaccharide challenge . Furthermore, single gene expression studies have suggested that enhancing maternal supply of Met could promote the calf’s ability to quickly adapt to extrauterine life [10, 30, 31]. Lastly, Met supplementation as RPM altered the transcriptome of bovine preimplantation embryos harvested at 70 days postpartum . Although these findings provided some evidence that methyl donors could play a role in nutritional programming in dairy cows, knowledge of the underlying mechanisms between late-gestation methyl donor supply and fetal programing in dairy cattle is still in its infancy.
Since nutritional management of modern dairy cows entails dietary manipulation of energy density and nutrients such as essential AA during the last stages of gestation (~ 4–6 weeks prepartum) , a deeper investigation of the biological outcomes on the neonatal calf is warranted. Particularly considering possible contributions of maternal nutrition to the calf’s immune and sanitary challenges during their first weeks of life . In this context, the use of RNA sequencing technology (RNA-Seq) has already proven to be a promising tool in helping us detect offspring genome-wide alterations in response to maternal post-ruminal Met supply .
In the present work, a subset of calves from a larger cohort [24, 27] was used to investigate the effect of maternal post-ruminal Met supply during late-pregnancy (− 28 ± 2 d to parturition) on changes in plasma systemic physiological indicators, transcriptome profiles, DNA methylation, one-carbon metabolism enzyme activities and protein abundance of nutrient-sensing pathways in the liver of new-born calves. Our general hypothesis was that Met supplementation as RPM during late-gestation improves liver immunometabolic functions in the offspring similar to those observed in the cow , in particular affecting key metabolic and immunological pathways such as one-carbon metabolism and transsulfuration reactions.
Growth performance, blood biomarkers and AA concentrations in plasma
At birth, calves born to dams fed MET had greater hip height (P value = 0.04) and wither height (P value = 0.01). No significant differences were detected for body weight and length, and hip width (P value ≥ 0.10) (Table 1). Calves in MET (from cows fed additional Met) had a tendency for lower concentration of glucose compared with CON (control) calves (P value = 0.08) at day 2, whereas no significant differences were detected for other blood parameters (Table 2). In this regards, it is interesting to note that no significant differences were detected for insulin concentrations (Table 2). At day 2, there was no significant effect of maternal diet for any AA concentration in the plasma; however, there was a tendency for higher concentrations of Phenylalanine (Phe; P value = 0.08) and Taurine (Tau; P value = 0.06) in MET calves compared with CON (Table 3).
Global DNA methylation and western blotting
Maternal supplementation with Met led to greater (P value < 0.05) global DNA methylation compared with CON calves (Fig. 1). Among the proteins measured, the ratio of p-AKT:AKT (AKT Serine/Threonine Kinase) was greater in MET calves (P value < 0.001; Table 4). In contrast, MET calves had a lower ratio of p-S6K1:S6K1 (P value = 0.01; Table 4).
Metabolomics, hepatic enzyme activity and mRNA abundance
At day 4, maternal supplementation with MET led to greater concentrations (P value ≤0.05) of Glycine, Adenosine, Serine, Taurine, Cystathionine, Glutamate, Fumarate, NAD, NADH, Taurocholic acid, Glycocholic acid, Lithocholic acid, and Glycochenodeoxycholic acid (Table 5). In contrast, lower hepatic activity for cystathionine β-synthase (CBS) and 5-methyltetrahydrofolate homocysteine methyltransferase (MTR) was detected in the MET calves (Table 6). Maternal supplementation with MET led to a lower abundance of Phosphate Cytidylyltransferase 1, Choline, Beta (PCYT1B; P value = 0.04), an overall tendency toward lower abundance of Cysteine Dioxygenase Type 1 (CDO1) and greater abundance of Phosphatidylethanolamine N-Methyltransferase (PEMT) and Methionine Adenosyltransferase 1A (MAT1A) respectively (P value ≤0.10; Table 7).
RNA sequencing and gene expression analyses
A summary of sequencing read alignment and mapping is reported in Additional File 1. Overall, samples had approximately 12 million reads of which 11 million (~ 95%) were uniquely mapped and 9.4 million (~ 78%) assigned to genes. Statistical analysis identified 13,867 uniquely annotated (EntrezID) genes. Of these, applying the 0.05 FDR cut-off, 74 genes (36 upregulated and 38 downregulated) were detected as differentially expressed (DEG) comparing MET with CON heifer calves, whereas 568 DEG (273 upregulated and 295 downregulated) were detected at FDR ≤ 0.10 cut-off (Fig. 2). A summary of the top-ten up- and downregulated genes (FDR ≤ 0.10) is reported in Tables 8 and 9. The entire list of DEG is reported in Additional File 2.
KEGG pathway analysis
The Dynamic Impact Approach (DIA) analysis yields the impact and flux of all the manually-curated pathways included in the KEGG database. The term ‘impact’ refers to the biological importance of a given pathway as a function of the change in expression of genes composing the pathway (proportion of DEG and their magnitude) in response to a treatment, condition, or change in physiological state. Consequently, the direction of the impact, or flux, characterizes the average change in expression as up-regulation/activation, down-regulation/inhibition, or no change. Considering DIA results with DEG at FDR ≤ 0.10, a broad effect on the transcriptome due to maternal MET was detected (Fig. 3). All main KEGG categories, both metabolic (‘Metabolism’) and non-metabolic (‘Environmental information processing’, ‘Cellular processes’, and ‘Organismal systems’) were broadly impacted with a negative flux (i.e. downregulated), except for ‘Genetic information processing’ that was markedly upregulated (Fig. 3). In particular, ‘Cell growth and death’, ‘Lipid Metabolism’, ‘Aging’, ‘Metabolism of Cofactors and Vitamins’ and ‘Signaling Molecules and Interaction’ were the most-impacted and downregulated KEGG subcategories (Fig. 5). Along with ‘Genetic Information Processing’ pathways, the ‘Nucleotide Metabolism’, ‘Glycan Biosynthesis and Metabolism’, ‘Environmental adaptation’, ‘Immune System’ and ‘Nervous system’ were the most-impacted subcategories with a positive flux (i.e. upregulation) in MET vs. CON comparison (Fig. 3). Top-10 up- and downregulated single pathways are reported in Figs. 4 and 5, respectively.
Transcription regulator discovery
The transcription factor enrichment analysis with the ChIP-X Enrichment Analysis 3 (ChEA3) tool generated a list of 72 TF significantly-associated (FDR ≤ 0.05) with our DEG at the FDR ≤ 0.10 threshold. The list of top 10-ranked upstream regulators is reported in Table 10, whereas the entire list of TF is reported in Additional File 3. Applying the DIA approach to the ChEA3 results, the TF impact and flux values were predicted (Table 10).
The present findings were broadly consistent with our previous reports investigating the effect of feeding MET during late-pregnancy on cow and calf hepatic function [10, 24, 26, 27, 30, 31, 34, 35]. Briefly, in those previous studies enhanced MET supply improved immunometabolism along with DMI in dairy cows during the peripartal period and through peak lactation [24, 34, 35]. In the larger cohort of calves from this study encompassing birth through the first 9 wk. of life, we reported that maternal MET supply influenced enzyme activity and metabolome in one-carbon metabolism and the tricarboxylic acid (TCA) cycle, with beneficial physiological advantages to calves . Furthermore, previous results on hepatic target-gene transcription have suggested that feeding Met during late-pregnancy was associated with faster maturation of key metabolic pathways involved in the calf’s ability to quickly adapt to extrauterine life [10, 30]. Thus, RNA-seq results alone revealed that feeding Met during late-gestation broadly altered neonatal calf liver transcriptome profiles. In particular, transcriptome profiles confirmed the hypothesis of an enhanced immunometabolic status attributable to the change in expression profiles of several genes mainly involved in ‘Glucose metabolism’, ‘Lipid metabolism, ‘Glutathione’, and ‘Immune System’ metabolism.
Although the focus of the present integrative analyses only encompassed day 4 of age, together, the data supported the idea that feeding Met to enhance post-ruminal supply in the cow during late-gestation primed or programmed the Met cycle in calf liver, hence, contributing to better rates of growth and development . For instance, the greater concentrations of adenosine and serine observed in MET calves at this early age support the hypothesis of a priming effect of the Met cycle. The conversion of Met to SAM is accompanied by ATP consumption in a reaction controlled by Met adenosyltransferase (MAT)  and after demethylation, SAM is converted to SAH, which is subsequently hydrolyzed to homocysteine with the production of adenosine . The tendency for greater MAT1A abundance in MET calves (P value < 0.10; Table 7) supported this idea. In addition, the greater serine concentrations observed in MET calves supported this hypothesis, since it is known that serine supports the Met cycle by providing one-carbon units to regenerate Met from homocysteine and through de novo ATP synthesis . In this context, the downregulation of Thymidylate synthase (TYMS) [FC = − 1.50] was also noteworthy. TYMS converts deoxyuridine monophosphate (dUMP) to deoxythymidine monophosphate (dTMP) in a 5,10-methylene-tetrahydrofolate (THF)-dependent reaction . Its downregulation indirectly suggested a greater concentration of the most reduced form of folate 1C unit, 5-methyl-THF that has a unique cellular fate, the remethylation of homocysteine to form methionine . In light of current and previous observations, we speculate that despite the lower activity of MTR in liver from MET calves, enzymatic efficiency might have been greater. More in-depth discussion on this point is available in Additional File 4.
The greater concentration of serine also indicated the potential for greater metabolic activity through the transsulfuration pathway in MET calves, an idea supported by the greater concentration of a number of metabolites including cystathionine and taurine. The former is the product of the CBS reaction using homocysteine and serine, which is considered rate-limiting in the transsulfuration pathway . Thus, despite the lower CBS activity, based on the greater concentrations of a number of intermediate metabolites (e.g. cystathionine, cysteinesulfinic acid), we speculated that efficiency of the CBS reaction was such that flux through the transsulfuration pathway was overall greater in MET calves. As such, concentrations of the cellular antioxidant GSH  also increased, and could have elicited beneficial effects on antioxidant responses in the calves.
A greater degree of Met metabolism was also indirectly suggested in our experiment by the marked upregulation of the nucleotide metabolism pathway (Fig. 4), since the role played by Met in Purine synthesis is well-documented . Considering that Met plays an essential role in epigenetics via DNA methylation  and that it has been demonstrated that maternal methyl donor supplementation during pregnancy can regulate epigenetics via DNA methylation , the overall upregulation of ‘Genetic information processing’ pathways was compatible with the greater global DNA methylation detected in MET calves (Fig. 1). The evident upregulation of TRMT61A (tRNA Methyltransferase 61A) [FC = 1.69] was consistent with this scenario, as well as the upregulation of POLR3H (RNA Polymerase III Subunit H) [FC = 1.72] and POLR2I (RNA Polymerase II Subunit I) [FC = 1.71] (Additional File 2). In this context, the marked upregulation of histone HIST1H2AC [FC = 2.27] was noteworthy considering that SAM is also the most-important methyl-donor for histone methyltransferases . Recent studies demonstrated how increased SAM and histone methylation levels are associated with genetic ablation of MAT in Saccharomyces pombe . In the general context of nutritional programming, the marked upregulation of Meiotic recombination REC8 [FC = 3.88] and downregulation of cyclin B1 (CCNB1) [FC = − 3.59] and cyclin-dependent kinase 1 (CDK1) [FC = − 3.15] was intriguing. These changes could potentially be biologically-relevant considering the role of these genes in oocyte meiosis [46, 47] and that methionine adenosyltransferase 2B (MAT2B) influences oocyte maturation in mouse by regulating the MAPK pathway . Further studies will have to be performed to ascertain the mechanistic relevance of these changes in the context of hepatic nutritional programming of the calf.
Studies in humans and mice suggest that DNA methylation plays a crucial role in tissue-specific insulin-induced gene expression . Both in vivo and in vitro studies confirmed that AKT is an indicator of hepatic and adipose insulin sensitivity in dairy cows [50, 51]. Thus, the greater phosphorylated AKT ratio in MET calves suggested they were more sensitive to insulin, which is consistent with the tendency for lower plasma glucose concentrations in those calves (P value = 0.08; Table 2) despite the lack of statistical difference in insulin concentrations. In this regard, it is also noteworthy that (at least in non-ruminants) AKT phosphorylates and stimulates sterol regulatory element binding transcription factor 1c (SREPB1c) leading to enhanced liver lipogenesis through the suppression of insulin induced protein 2 (INSIG2). INSIG2 is a protein of the endoplasmic reticulum that blocks the activation of SREBP1c by binding to SREBP cleavage-activating protein (SCAP) and prevents it from escorting SREBPs to the Golgi . This scenario was confirmed in our experiment by the significant downregulation of INSIG2 [FC = − 1.35], and the fact that at this age of the calf most of the supply of fatty acids (FA) are derived from consuming milk replacer.
The lower phosphorylated S6K1 (Ribosomal Protein S6 Kinase) ratio in MET calves was intriguing considering that S6K1 is the predominant regulatory kinase of GSK3 (Glycogen synthase kinase 3)  and that AKT is known to inactivate GSK3β , one of the two isoforms of GSK3. In this context, the significant downregulation of PCK1 (Phosphoenolpyruvate Carboxykinase 1) [FC = − 1.99] is in line with the fact that most of the circulating glucose in newborn calves arises from lactose in the milk replacer fed. The downregulation of FOXO4 (Forkhead Box O4) [FC = − 1.44] also agreed with this idea, since the FOXO complex binds the PCK1 promoter .
It is well-established that lactose intake on its own is not sufficient to meet the newborn glucose demands , thus, calves have to establish endogenous glucose production and gluconeogenesis . In this regard, hepatic glycogen stored during late-gestation represents an important energy source that is quickly exhausted [57, 58] until the gluconeogenic response is better developed. Typically, an increase of glucose-6-phosphatase (G6Pase) activity mirrors the decline of hepatic glycogen storage soon after birth . This scenario was confirmed in our experiment by the significant downregulation of GYS2 (Glycogen Synthase 2) [FC = − 1.39] and the simultaneous upregulation of G6PC3 (Glucose-6-Phosphatase Catalytic Subunit 3) [FC = 1.30].
Mechanistically, greater activity of G6P phosphatase favors conversion of glucose-6-phosphate to glucose and ensures the release of glucose into the circulation . Nevertheless, the downregulation of PYGL (Glycogen Phosphorylase L) [FC = − 1.32], along with upregulation of NR0B2 (Nuclear Receptor Subfamily 0 Group B Member 2) [FC = 2.64] and PPP2R3C (Protein Phosphatase 2 Regulatory Subunit B″ Gamma) [FC = 1.68] was intriguing. These responses seemed to suggest a better capacity of MET calves in the utilization of lactose for meeting glucose requirements with a consequent reduction in depletion of glycogen reserves. The PYGL gene encodes a protein that catalyzes the cleavage of alpha-1,4-glucosidic bonds to release glucose-1-phosphate from liver glycogen stores and plays an important role in maintaining blood glucose homeostasis . NR0B2 encodes a protein that can directly bind to cAMP response element-binding protein (CREB) to block the association with CREB regulated transcription coactivator 2 (CRTC2), leading to an inhibition of hepatic gluconeogenic mRNA abundance [61, 62]. This was in line with the fact that AKT is known to phosphorylate and inhibit CRTC2 whose disruption is known to improve insulin sensitivity . PPP2R3C encodes a subunit of the PP2A (protein phosphatase 2) holoenzyme involved in glucocorticoid receptor (GR) feedback regulation, and particularly in GR-mediated negative regulation . In this regard, it is known that glucocorticoids in neonatal calves through activation of the GR can influence many important liver functions such as gluconeogenesis . Particularly neonatal glucose metabolism when nutrient supply is insufficient . Overall, this idea was also supported by our metabolomics results. Indeed, in absence of differences in feed intake, the greater concentrations of the tricarboxylic acid cycle intermediates fumarate and glutamate along with NAD/NADH in MET calves indicated enhanced rates of energy metabolism .
Lactose intake on its own is not sufficient to meet glucose demands  and mitochondrial FA oxidation increases after birth to cover energy needs . In our experiment, we detected an overall downregulation of ‘Lipid metabolism’ in MET calves (Fig. 3; Additional File 5), and notably of ‘Fatty acid degradation’ genes such as ACOX1 (Acyl-CoA Oxidase 1) [FC = − 1.53], CPT1A (Carnitine Palmitoyltransferase 1A) [FC = − 1.46], and ACADM (Acyl-CoA Dehydrogenase Medium Chain) [FC = − 1.68] all of which are well-known to be crucial for mitochondrial and peroxisomal FA oxidation . There was also a significant downregulation of PPARA (Peroxisome Proliferator Activated Receptor Alpha) [FC = − 1.33] that encodes a nuclear receptor controlling genes required for activation of mitochondrial and peroxisomal FA oxidation . In this context, the downregulation of Gamma-Butyrobetaine Hydroxylase 1 (BBOX1) [FC = − 1.51], which in non-ruminants is regulated by PPARA , was noteworthy. This gene encodes the rate-limiting enzyme in the synthesis of carnitine from butyrobetaine , an essential nutrient for CPT1 activity and mitochondrial FA oxidation [69, 70]. In the absence of differences in hepatic concentration of butyrobetaine (Table 5) and fat intake, and because calves were consuming a high-fat milk replacer (15% crude fat on dry matter basis), the significant upregulation of MCAT (Malonyl-Coa-Acyl Carrier Protein Transacylase) [FC = 1.42] agreed with the downregulation of CPT1A. Indeed, MCAT is a component of the fatty synthase complex  and is specifically involved in the step where malonyl-CoA is loaded into the enzyme complex at the beginning of the FA biosynthesis process . Considering that most of the FA found in the liver would have come as chylomicron remnants or from lipogenesis , it is plausible that lipogenesis from intake of lactose and glucose availability in the liver would have increased malonyl-CoA concentration, an allosteric inhibitor of CPT1A and mitochondrial FA oxidation . Although we did not measure proteins involved in mitochondrial or peroxisomal FA oxidation, the downregulation of the above well-known FA oxidation genes that are transcriptionally-regulated suggested that the most logical alternative for the liver to handle FA coming from the diet in the MET calves might have been through synthesis of very low density lipoproteins (VLDL), phospholipid species, and/or incorporation into cellular membranes. The former hypothesis, formulated previously , was corroborated in our experiment by the tendency for greater mRNA abundance of MAT1A and PEMT (P value < 0.10; Table 7) in MET calves. These genes, indeed, play an important role in hepatic phosphatidylcholine (PC) synthesis [76,77,78], which is known to participate in VLDL synthesis and help prevent fatty liver [79, 80]. Endogenous synthesis of PC from supplemental Met  was already proposed to play a role in the ability of the cow liver to handle influx of FA produced from lipolysis after parturition . More in-depth discussion about possible effects on increased VLDL export through the generation of PC is available in Additional File 4.
The marked downregulation of Fatty Acid Desaturase 2 (FADS2) [FC = − 4.02] was also intriguing. This gene encodes the delta-6 desaturase enzyme, which catalyzes the first step in the desaturation and elongation of the dietary essential FA . Interestingly, a different pool of PC molecular species is synthesized in the liver and is partly dependent on the pathway of synthesis: the cytidine diphosphocholine (CDP)-choline and phosphatidylethanolamine N-methyltransferase (PEMT) pathways . PC synthesized via the PEMT pathway has higher concentrations of long-chain and highly-unsaturated species , whereas PC produced via the CDP-choline pathway has higher monounsaturated and saturated FA . Thus, we speculate that lower FADS2 expression may influence the availability of the preferred substrate required for liver PC synthesis in the PEMT pathway. In this regard, a diminished methylation capacity and hypermethylation silencing of FADS2 mRNA expression has been observed in mice with CBS deficiency .
Amino acid metabolism
Our functional analysis indicated that ‘Amino Acid metabolism’ was overall inhibited (Fig. 3; Additional File 5), due to marked inhibition of the ‘Valine, leucine and isoleucine degradation’ pathway. Leucine, isoleucine and valine are branched-chain amino acids (BCAA) with a number of biological functions beyond protein synthesis . Furthermore, although little is known about the impact of changes in BCAA availability on immune system functionality in ruminants, both human and rodent studies have underscored the relationship between BCAA availability and improved immune function [88,89,90,91], highlighting the need for further studies about BCAA effects on dairy cows during the periparturient period and on newborn calves.
An increased concentration and utilization of essential AA, particularly BCAA, was observed in immortalized bovine mammary epithelial cells when the supply of Met increased . The authors speculated this change was potentially mediated by alterations in AA transporters controlled by mTOR whose activation was observed when Met supply increased [92,93,94]. In this context, the significant upregulation of Late Endosomal/Lysosomal Adaptor, MAPK and MTOR Activator (LAMTOR2) [FC = 1.36] was noteworthy. The protein encoded by this gene is part of the regulatory complex involved in AA sensing and activation of mTORC1 , a signaling complex promoting cell growth in response to growth factors, energy level, and supply of AA .
Although a greater mRNA abundance of AA transporters was observed in the placenta of MET cows , which we speculated was partly due to greater DMI, overall, our results indicated that feeding RPM to enhance post-ruminal supply of MET did not affect hepatic AA transport in calf liver. In this regard, it is important to note that our calves were fed a high-protein milk replacer (28% crude protein on dry matter basis). Without differences in intake, our results suggested a lower utilization of AA in MET calf liver. This idea is supported by the lower p-S6K1, since it is well known the central role played by S6K1 in the regulation of protein translation in mammals . The numerically-lower expression of pRPS6 (P value = 0.15) also partly corroborated this idea. Indeed, RPS6 was the first S6K substrate identified . Considering our scenario of lower AA utilization by the liver, we speculate about the presence of possible AA ‘sparing effects’ in MET calves, i.e. more essential AA available for other organs, e.g. skeletal muscle for growth. This hypothesis agrees with the greater rate of growth already observed in utero and during early postnatal life in calves born to cows fed RPM to enhance post-ruminal Met supply during late gestation [26, 27].
Although no significant differences between CON and MET calves were detected for GSH concentration, upregulation of GSTO1 (Glutathione S-Transferase Omega 1) [FC = 2.85] and GPX1 (Glutathione Peroxidase 1) [FC = 1.56; Additional File 2] both of which are involved in ‘Glutathione metabolism’ underscored the importance of this pathway in the MET neonatal liver. In non-ruminants, the role of GSH (a potent intracellular antioxidant in liver tissue) has been extensively described, in particular its association with methyl-transferase products involved in ‘Cysteine and Methionine metabolism’ . In this regard, the tendency for greater taurine production in MET calves (P value = 0.06; Table 3) suggested possible beneficial effects on antioxidant responses during early postnatal life. An improved antioxidant status has been observed in peripartal dairy cows receiving enhanced post-ruminal Met supply [22, 35], and an increased in GPX activity also has been highlighted .
The idea that endogenous taurine plays a role in the antioxidant response of calves was formulated considering enzyme activity and metabolomics in the bigger cohort of calves from the present study . It was further supported by data from an in vitro study in which taurine supplementation mitigated the inflammatory activation of blood polymorphonuclear leukocytes from neonatal Holstein calves . More in general along with taurine, despite the lower CBS activity, the greater concentration of cystathionine and serine detected in MET calves suggested a greater activity of the transsulfuration pathway in MET calves, which is the main source of the cellular antioxidant GSH . This scenario in our experiment supported the hypothesis of a better antioxidant status in MET calves, also suggested by the marked downregulation of Catalase (CAT) [FC = − 1.77; Additional File 2]. In humans, CAT is the most-adaptive antioxidant enzyme in the presence of oxidative stress and plays a pivotal role in cellular defense against oxidative damage .
The increase of antioxidant enzymes such as CAT can be considered a compensatory response to increased oxidative stress . In non-ruminants, hepatic CAT activity increased during oxidative impairment due to an inhibition of CBS, although no variation in mRNA abundance  and no effect on GHS peroxidase, GHS reductase and GHS S-transferase activities were reported . This latter pattern seemed to be recognizable in our experiment, suggesting that lower CBS activity was not associated with an alteration in epigenetic methylation as reported previously . This is particularly noteworthy considering the role played by CBS in the transition from methylation to transsulfuration. More in-depth discussion about CBS activity is available in Additional File 4.
Overall, our results seemed to support the idea that feeding RPM to enhance post-ruminal supply of Met (along with other nutrients due to the greater DMI) may promote, through the transsulfuration pathway, production of the antioxidants GSH and taurine in calves . The increased potential for antioxidant production may help alleviate oxidative stress and inflammation induced from reactive oxygen metabolite production in the liver [35, 106, 107] particularly at birth when antioxidant mechanisms in the new-born are still immature [108, 109]. Although downregulation of various genes associated with mitochondrial FA oxidation was observed in MET calves, the fact that high intracellular levels of FA are a source of reactive oxygen species  potentially via peroxisomal oxidation underscores the importance for production of antioxidants such as GSH and taurine. This hypothesis was also broadly supported in our experiment by the marked downregulation of ‘p53 signaling pathway’ (Fig. 5), since it is known that p53 activation is induced by multiple stress signals  and notably by oxidative stress .
Furthermore, considering that induction of thermogenic genes was recently observed as a consequence of taurine supplementation , and that liver is one important hub in thermogenesis , the overall upregulation of ‘Enrivonmental Adaptation’ function (Fig. 3; Additional File 5) due to changes in the ‘Thermogenesis’ pathway suggested that adaptive functions also might have been enhanced by feeding RPM to enhance post-ruminal supply of Met in the pregnant cow. The importance of thermoregulation in the newborn is well-establish as well as its linkage with energy sources available at birth , thus, it is plausible that a better energetic and immunometabolic status in MET calves may contribute to the thermoregulation.
New-born calves are immunologically naïve at birth , and particularly vulnerable since they have to face many stressors while their immune system is still immature . Birth is the first direct and critical stressor in the calf’s life that has to be faced through innate immunity, a process that functions less effectively in adults, and passive immunity, i.e. the ingestion of colostrum . In this regard, it is worth noting that feeding RPM in late-pregnancy does not seem to impact colostrum yield, colostral FA or AA profiles, colostral immunoglobulin G (IgG) concentrations, or apparent IgG absorption by the calf [27, 30], suggesting that neither the increase in post-ruminal Met nor total nutrient supply (due to greater DMI) affects the ability of the calf to acquire passive immunity. This is in line with the absence of difference in plasma GGT concentration in MET or CON calves (Table 2). Activity of GGT is a suitable indicator of passive immunity transfer in calves, and is highly- and positively-correlated with plasma IgG . At least in non-ruminants, the liver is a key frontline immune tissue, balanced between an anti-inflammatory and immunotolerant status [118, 119]. In our study the overall upregulation of various pathways involved in ‘Immune system’ (Fig. 3; Additional File 5) and notably of ‘C-type lectin receptor signaling’ and ‘Natural killer cell mediated cytotoxicity’ pathways (Fig. 4) seemed to support our general hypothesis of an enhanced immunometabolic status in response to feeding RPM to enhance post-ruminal supply of Met in late-pregnancy. Indeed, these pathways are well-known to be involved in innate immunity [120,121,122].
The upregulation of the ‘Natural killer cell mediated cytotoxicity’ pathway was especially noteworthy considering that natural killer cells are mainly present in liver tissue and play a key role in the innate immune system [121, 123, 124]. Increasing evidence of cross talk between NK cells and other immune cells suggested that NK cells also play an important role in shaping the adaptive immune response [125,126,127]. On the other hand, C-type lectin receptors (CLR) are powerful pattern-recognition receptors that mediate immune recognition and response , whose role in orchestrating the induction of signaling pathways regulating adaptive immune responses is also well-recognized [129, 130].
The evident upregulation of the ‘Glycan Biosynthesis and Metabolism’ pathway (Fig. 3) agreed with the general scenario of an earlier and more effective immune status. Indeed, the role of glycans in the regulation of both innate and adaptive immune responses is well-documented in non-ruminants , in particular it is well-known that the immune system is highly controlled and fine-tuned by glycosylation through the addition of a diversity of carbohydrate structures (glycans) to the immune cell receptors . This was consistent with the upregulation of the ‘C-type lectin receptor signaling’ pathway in our experiment, since C-type lectins are a class of glycans binding proteins . In addition, considering the role of the phagosome in linking innate and adaptive immunity , the upregulation of the ‘Phagosome’ pathway (Fig. 4) also supported the hypothesis of an enhanced immune status in MET calves.
The general downregulation of ‘Platelet activation’ (Additional File 5) broadly agreed with this scenario, since liver dysfunction is usually reported in combination with an activation of platelets, often considered dynamic sentinels interacting with immune cells [134,135,136]. In this context, the marked upregulation of M-SAA3.2 (Mammary Serum Amyloid A3.2) [FC = 5.70] was noteworthy (Table 8). SAA is an acute-phase protein which is highly-conserved among all vertebrate species . In this regard, it is important to note that, if this response translated into more SAA protein secreted from the liver, it would mean that AA within liver would have been channeled away from protein synthesis, as suggested by the lower p-S6K1 observed in MET calves. Although the functional significance of SAA production remains speculative, several authors suggested its plausible role as an immunological defense molecule providing immediate defense against inflammatory challenges [138, 139], with a protective role during inflammation  and particularly with a possible beneficial function in the well-being and better extrauterine life adaptation in calves . More in-depth discussion on single target genes involved in the general context of enhanced immune status is available in Additional File 4.
Transcription factor enrichment analysis with the ChEA3 tool  allowed us to predict 72 TF significantly associated (FDR ≤ 0.05) with our DEG list (Additional File 3). A similar approach to the DIA analysis  was also applied to predict the potential impact and flux for each highlighted TF. Focusing on the top-10 ranked TF (Table 10) and based on information retrieved from the published literature, many of them are consistent with our general scenario of a higher immunometabolic state, greater insulin sensitivity and lower mitochondrial FA oxidation in MET calves. Among other TF identified were FOXO1 (Forkhead Box O1), PPARG (Peroxisome Proliferator Activated Receptor Gamma), E2F1 (E2F Transcription Factor 1), and CREB1 (CAMP Responsive Element Binding Protein 1). All these have relatively well-established associations with important metabolic and cellular processes in non-ruminants.
At least in non-ruminants, FOXO1 is a key TF associated with coordinating lipid metabolism in the post-absorptive state, integrating insulin signaling with glucose and lipid metabolism . In liver, insulin activates the PI3K/PKB signaling pathway  and results in FOXO1 protein phosphorylation and degradation . FOXO1 binds and promotes transcription of PCK1, which plays a key role in gluconeogenesis . Indeed, hepatic FOXO1 loss-of-function mutant suppresses PCK1 expression determining a decreased hepatic gluconeogenesis . Also the regulation of pyruvate dehydrogenase kinase (PDK) expression by FOXO signaling has been recently investigated and the inhibition of PDK by insulin via phosphorylation of FOXO through PI3K/Akt signaling pathway has been described [147, 148]. The downregulation of PDK1 [FC = − 1.45] in our experiment was in line with this scenario. Furthermore, at least in non-ruminants, FOXO1 regulates transcription of genes involved in hepatic assembly of VLDL  and also plays a fundamental role in development and differentiation of immune cells .
The nuclear receptor PPARG belongs to peroxisome proliferator-activated receptors, members of the superfamily of ligand-activated nuclear transcription factors. Many PPAR-regulated genes encode proteins that regulate FA oxidation (mitochondrial and peroxisomal) and storage, for example, ACOX1 and CPT1A. Furthermore, PPARG regulates pantothenate kinase (PANK) , which catalyzes the rate-controlling step in coenzyme A (CoA) biosynthesis . In our experiment, we detected downregulation of PANK3 [FC = − 1.87], which would be metabolically important because CoA is an essential cofactor supporting a multitude of oxidative and synthetic reactions, including those involved in FA biosynthesis and oxidation . Taking all of this into account, the predicted inhibition of PPARG seemed particularly in line with our general scenario of FA oxidation downregulation.
Among the list of top-10 ranked TF, the presence of E2F1, a well-known TF promoting hepatic gluconeogenesis  and innate immune response  is noteworthy. Similarly, the predicted activation of CREB1 whose role in coordinating hepatic lipid and glucose metabolism through inhibition of PPARG has been documented  along with its role in the immune system through promoting anti-inflammatory responses .
Feeding RPM to enhance maternal post-ruminal Met supply was associated with greater global DNA methylation and distinct hepatic transcriptome, proteome, and metabolome profiles after birth. The hypothesis that enhancing post-ruminal Met supply during late-gestation primed or programmed the Met cycle in calf liver was corroborated. However, in light of the greater DMI in response to feeding RPM (a consistent response across published studies), the increased supply of other nutrients and their potential impact cannot be disregarded. Molecular changes with potential beneficial effects to overall hepatic function were highlighted, especially in terms of metabolism and inflammatory status. This was attributable to the change in expression profiles of several genes mainly involved in ‘Glucose’, ‘Lipid’, ‘Glutathione’, and ‘Immune System’ metabolism. In particular, the results seemed to suggest a greater flux through the TCA cycle along with a better antioxidant status due to taurine and GSH synthesis. Although further studies are required to corroborate this result, the upregulation of ‘Immunity System’ pathways supported the hypothesis that feeding RPM to enhance maternal Met supply may play a role in regulating fetal immunity.
All procedures for this study were conducted in accordance with a protocol approved by the Institutional Animal Care and Use Committee (IACUC) of the University of Illinois (protocol # 14270). Heifer calves in the present study were a subset from a larger cohort of animals born to Holstein cows randomly assigned to receive a basal control (CON) close-up diet (from − 28 ± 2 d to parturition) [1.47 Mcal/kg dry matter (DM) and 15.3% crude protein (CP)] with no added Met or CON plus ethyl cellulose MET (MET, Mepron®, Evonik Nutrition & Care GmbH, Germany) . Diet composition is available in Additional File 6. All management procedures have been described previously . Briefly, during the preliminary period from − 45 to − 29 d relative to parturition all cows received a common low-energy and high-straw far-off diet (1.33 Mcal/kg of DM and 13.9% CP) with no RPM. Cows were individually-fed using Calan gates (American Calan Inc., Northwood, NH). At − 28 days relative to parturition, cows were randomly assigned to CON or MET groups. The Met product was offered at a rate of 0.09% of previous day DMI. This supply of Met was based on studies demonstrating beneficial effects on production performance and health during the prepartum period [23, 158]. Mepron® is a commercial rumen-protected source of DL-Met in the form of small beads containing a minimum of 85% methionine that resists ruminal degradation due to an ethyl-cellulose film coating the methionine core. The intestinal digestibility coefficient of Mepron® is 90%  and its ruminal bypass value is 80% . Feeding this amount of RPM increased Met concentrations in plasma prior to and during the first 3 weeks postpartum .
After parturition, neonatal calves were separated from their dams. Body weight (BW), hip and wither height (HH, WH), hip width (HW) and body length were measured. Heifer calves were kept in the experiment if they fulfilled all the following criteria previously described , 1) single heifer calf; (2) heifer calf birth weight > 36 kg; (3) calving difficulty score < 3; and (4) dam first colostrum volume > 3.8 L. All calves were managed in the same fashion. At birth, the navel was disinfected with 7% tincture of iodine solution (First Priority Inc., Elgin, IL, United States), and neonatal calves were vaccinated with TSV II (Pfizer Inc., New York, NY, United States) via nostril application.
Calves were offered 3.8 L of first milking colostrum from their mother within 6 h. If voluntary colostrum intake had not reached the 3.8 L required, calves were force-fed via an esophageal tube to ensure that all calves consumed the same amount of colostrum. Calves were housed in individual outdoor hutches bedded with straw and fed twice daily (AM and PM) with a milk replacer (Advance Excelerate, Milk Specialties, Carpentersville, IL; 28.5% CP, 15% fat). More in detail, heifer calves received 4.54 kg/day of milk replacer mix (0.59 kg of milk replacer in 3.95 L of water) and had ad libitum access to starter grain mix [Ampli-Calf Starter 20®; 19.9% CP and 13.5% neutral detergent fiber (NDF), Purina Animal Nutrition, Shoreview, MN, United States] fed in the morning. All heifer calves remained clinically healthy during the study.
Per IACUC guidelines, blood samples were collected (n = 6 calves per group) from the jugular vein at 2 d of age using 20-gauge BD Vacutainer needles (Becton Dickinson, Franklin Lakes, NJ). Vacutainer tubes used contained lithium heparin, and plasma was obtained by centrifugation at 1900×g for 15 min at 4 °C and stored at − 80 °C until further analysis. Plasma was used to analyze the concentrations of a number of biomarkers associated with energy metabolism (e.g. NEFA, BHBA), immune function, liver function, and oxidant status according to protocols described in a number of publications from our group [9, 10, 30, 31]. In addition, plasma was used to profile concentrations of free AA and indicators of skeletal muscle protein turnover . Plasma insulin concentrations were determined using a bovine-specific commercial ELISA kit (catalog no. 10–1201-01; Mercodia, Uppsala, Sweden).
Per IACUC guidelines, surgical biopsies of the liver could not be obtained until day 4 of age. At that point, samples were obtained via puncture biopsy under local anesthesia (n = 6 calves per group) using published protocols from our group [26, 31]. Tissue was frozen immediately in liquid nitrogen and stored at − 80 °C until further analysis. Samples from the same calves (n = 6 per treatment) used for DNA methylation and RNA-sequencing were used for protein expression, enzyme activity, and metabolomics analyses. The chosen protein targets are key members of the nutrient signaling cascades that respond to insulin and AA supply, hence, allowed us to evaluate potential effects of maternal MET on key metabolic aspects of neonatal liver function. Since newborn calves do not have a functional rumen, insulin plays an essential role on aspects of hepatic nutrient metabolism including coordination of glucose, FA, and AA metabolism . Similarly, since our previous work with neonatal calves revealed unique adaptations in hepatic one-carbon metabolism in response to feeding RPM in late-pregnancy [10, 30], we elected to perform those assays along with targeted metabolomics analysis.
DNA isolation and global DNA methylation
The DNA was isolated from 50 mg of liver tissue using the Blood and Tissue DNeasy Kit (Qiagen). Methylation of genomic DNA was detected using the Methylamp global DNA methylation quantification ultra-kit (Epigentek) according to manufacturer’s instructions. Details on calculation were reported previously by our group .
The protocol used for western blot was the same as reported previously . Briefly, 50 g of liver tissue was used for total protein extraction using T-PER tissue protein extraction reagent (Cat. 78,510, Thermo Scientific) containing Halt protease and phosphatase inhibitor cocktail (100x, catalog no. 78442; Thermo Fisher Scientific). Concentration of total protein was determined using the Pierce BCA protein assay kit (catalog no. 23227; Thermo Fisher Scientific). Equal amounts of protein were separated using 4–20% Mini-Protean TGX precast gels (Cal. 4,561,096, Bio-rad) and transferred onto a PVDF membrane (Cat. 1,620,261, Bio-rad) using a semidry transfer cell (Serial No. 221BR, Bio-rad). After blocking with 5% non-fat dry milk in trisbuffered saline at room temperature for 2 h, the membrane was washed and incubated with a primary antibody directed against target antibodies. Antibodies raised from rabbit monoclonal antibody were from Cell Signaling Technology and diluted to 1:1000 prior to use [25, 92]. Targets were mTOR, p-mTOR (Ser2448), S6K, p-S6K (Thr389), RPS6, p-RPS6 (Ser235/236), 4EBP1, p-4EBP1 (Thr37/46), Akt, p-Akt (Ser473), eIF2, p-eIF2 (Ser51), eEF2, and p-eEF2 (Thr56) (catalog nos. 2972, 2971, 9202, 9234, 2217, 2211, 9452, 9459, 9272, 9271, 9722, 9721, 2332, and 2331, respectively; Cell Signaling Technology). The secondary antibody was HRP AffiniPure goat anti-rabbit IgG (ab6721, Abcam, diluted to 1:10,000). Visualization of the target proteins was performed using the Clarity Western ECL Substrate (Cat. 170–5060, Bio-rad). Each protein was normalized against GAPDH (Polyclonal antibodies raised in rabbit, ab22555, Abcam, diluted 1:2000). Data are reported as relative expression. The ECL signals were recorded using an imaging system (ChemiDocTM MP, Bio-rad) and analyzed using Image lab 5.2.1 (Bio-Rad, USA).
Approximately 100 mg of frozen tissue was extracted using a 2-step protocol described by Wu et al. . Targeted metabolomics (liquid chromatography–MS) was performed to quantify 31 metabolites related to TCA cycle, one-carbon metabolism, and transsulfuration pathway. Samples were analyzed with the 5500 QTRAP liquid chromatography–tandem MS system (Sciex, Framingham, MA) at the metabolomics core facility of the Roy J. Carver Biotechnology Center, University of Illinois (Urbana). Software Analyst 1.6.2  was used for data acquisition and analysis.
One-carbon metabolism enzyme activity assays
Activity of BHMT, MTR, and CBS was determined according to protocols detailed by Zhou et al. .
RNA extraction and quantitative PCR
All RNA extraction procedures were performed as described previously . Briefly, RNA extracted from liver samples using the miRNeasy kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol. Samples were treated on-column with DNaseI (Qiagen) to remove genomic DNA from the RNA. Quantification was accessed using the NanoDrop ND-1000 (NanoDrop Technologies, Rockland, DE), and RNA quality was measured using an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA). All samples had an RNA integrity number factor greater than 8, RNA was diluted to 100 ng with DNase/RNase-free water for cDNA synthesis by using RT-PCR, and the cDNA was diluted 1:4 with DNase/RNase-free water.
The quantitative PCR (qPCR) performed was SYBER Green based (Quanta Biosciences, Beverly, MA), using a 7-point standard curve obtained from a cDNA pool of all samples. The final data were normalized using the geometric mean of 3 internal control genes: UXT, GAPDH, and RPS9. A total of 24 genes selected related to the Met cycle, DNA methylation, transsulfuration, GSH, and cytidine 5′-diphosphocholine (CDP)–choline pathway were measured. The log2 transformed data were analyzed in order to compare the expression means of CON and MET groups. All evaluated genes and primer information are reported in Additional File 7.
Data were assessed for normality of distribution using the Shapiro-Wilk test. When the normality assumption was rejected, data were log2-transformed before statistical analysis and log2-back transformed after analysis. Differences between the means of the groups were analyzed via one-way ANOVA. Statistical differences were declared significant at P value ≤ 0.05 and tendencies at P value ≤ 0.10. All statistical analyses were performed in R environment .
Sequencing was performed by the High-Throughput Sequencing and Genotyping Unit of the W. M. Keck Biotechnology Center at the University of Illinois at Urbana Champaign (Urbana, IL, USA). A total of 12 mRNA libraries were prepared with Illumina’s ‘TruSeq Stranded Sample Prep kit’ following the manufacturer’s instructions. Libraries were sequenced on one lane for 101 cycles from one end of the fragments on a HiSeq2500 (Illumina Inc.), using HiSeq SBS sequencing kit version 4. Fastq files were generated and demultiplexed with the bcl2fastq (v18.104.22.168) software (Illumina). Approximately a total of 143 million single-reads sequences of 100 nt in length were collected. The quality control of the raw sequences was assessed with FastQC software (v0.11.15) and the reads were mapped to the bovine reference genome (UMD v.3.1.1; Ensembl Genomes website) using the software package STAR (v2.5.3a). The number of reads that mapped to each gene in each sample was calculated using the feature Counts function in the Subread package (v1.5.2). Data are available at the GEO site of NCBI (GSE163644).
Identification of differentially expressed genes
The entire downstream analysis for identification of DEGs was conducted within the R environment using the edgeR pipeline . Non-expressed and weakly-expressed genes were removed , TMM (trimmed mean of M-values) normalization was applied to all samples and the data were log2-transformed. Once dispersion estimates were obtained and a negative binomial generalized linear model was fitted with a design matrix describing the treatment condition (MET and CON groups), quasi-likelihood F-test was used to determine differential expression. Multiplicity correction was performed by applying the Benjamini-Hochberg method on the P values, to control the false discovery rate (FDR). The expression level was deemed to be significantly different between the two groups for FDR ≤ 0.10 .
Dynamic impact approach (DIA)
The DIA software was used for functional analyses . Briefly, DIA uses the KEGG systems information to rank biological pathways, calculating the overall impact (importance of a given pathway) and flux (direction of impact; e.g., upregulation, downregulation, or no change) based on gene fold change (FC) values. For this purpose, the entire DEG data set (FDR ≤ 0.10) with associated statistical P values was imported into DIA and the functional analysis was performed interrogating all KEGG pathway database.
Gene network analysis
The PANEV (Pathway Network Visualizer) tool  was extensively used to visualize the results in a context of gene/pathway networks, and pinpoint candidate genes associated with a subset of pathways of interest. Briefly, PANEV is an R package set for gene/pathway-based network visualization. Based on information available on KEGG, it visualizes genes within a customized network of multiple interconnected pathways, chosen by the user since supposed to be involved with the trait under investigation. The network graph visualization helped us to interpret functional profiles of cluster of genes. In our case, a network of interconnected (i.e. up- and downstream) pathways involved in ‘Glycolysis/gluconeogenesis’, ‘Lipid metabolism’, ‘Amino acid metabolism’ and ‘Immune system’ pathways was created to pinpoint functional candidate genes with coordinate expressions (Additional File 8).
Transcription regulator discovery
ChIP-X Enrichment Analysis 3 (ChEA3) tool  was used to predict the TF associated with the downstream DEG. ChEA3 is a transcription factor enrichment analysis software that ranks TF associated with user-submitted gene sets . The entire list of DEG (FDR ≤ 0.10) was uploaded on ChEA3 and the ‘Literature ChIP-seq’ library, containing ChIP-seq experiments from human, mouse and rat, was interrogated. ChEA3 core analysis was run using default settings and TF with FDR ≤ 0.05 were considered significantly associated with our DEG list. Furthermore, in order to estimate the state (i.e. activation/inhibition) of significantly-predicted TF, the same approach underlying DIA analysis  was applied to the candidate TF list. Briefly, impact and flux value were calculated for each predicted TF according to Bionaz et al. : the impact was obtained by combining the proportion of DEG with the log2 mean fold change and mean –log P value of all TF-target genes. More details about rationale and calculation steps of the analysis are reported in the original publication . The ChEA3 outcomes were matched with PANEV results in order to create a network of possible connections between TF and the downstream DEG involved with the ‘Glycolysis/gluconeogenesis’, ‘Lipid metabolism’, ‘Amino acid metabolism’ and ‘Immune system’ pathways (Additional File 9).
Availability of data and materials
The raw sequencing data of 12 animals in this study are available from the GEO site of NCBI under accession number GSE163644.
Eukaryotic translation initiation factor 4E-binding protein 1
Acyl-CoA dehydrogenase medium chain
Acyl-CoA oxidase 1
Angiopoietin-like protein 3
Analysis of variance
Butyrobetaine hydroxylase 1
Branched-chain amino acids
Cyclin-dependent kinase 1
Complementary deoxyribonucleic acid
Cysteine dioxygenase type 1
ChIP-X enrichment analysis 3
C-type lectin domain containing 6A
C-type lectin receptors
Carnitine palmitoyltransferase 1A
CAMP responsive element binding protein 1
CREB regulated transcription coactivator 2
Differentially expressed gene
Dynamic impact approach
Dry matter intake
E2F transcription factor 1
Eukaryotic elongation factor 2
Eukaryotic initiation factor 2
Fatty acid desaturase 2
False discovery rate
Forkhead transcription factor 1
Forkhead box O4
Glucose-6-phosphatase catalytic subunit 3
Glutamate-cysteine ligase modifier subunit
Glutathione peroxidase 1
Glycogen synthase kinase 3
Glutathione S-transferase omega 1
Glycogen synthase 2
H2A clustered histone 6
Institutional animal care and use committee
Insulin induced protein 2
Kyoto encyclopedia of genes and genomes
Krüppel-like factor 5
Late endosomal/lysosomal adaptor, MAPK and MTOR activator
Mitogen-activated protein kinase
Methionine adenosyltransferase 1A
Methionine adenosyltransferase 2B
Malonyl-CoA-acyl carrier protein transacylase
Maternal methionine supplied
Messenger ribonucleic acid
Mammary serum amyloid A3.2
Mechanistic target of rapamycin kinase 1
5-methyltetrahydrofolate homocysteine methyltransferase
Nicotinamide adenine dinucleotide
Nicotinamide adenine dinucleotide plus hydrogen
Neutral Detergent Fiber
Non-esterified fatty acids
Nuclear receptor subfamily 0 group B member 2
Nuclear receptor subfamily 3 group C member 1
Phosphorylated eukaryotic translation initiation factor 4E-binding protein 1
Phosphorilated serine/threonine-protein kinases
Pathway network visualizer tool
Phosphoenolpyruvate carboxykinase 1
Phosphate cytidylyltransferase 1, choline, beta
Pyruvate dehydrogenase kinase
Phosphorylated eukaryotic elongation factor 2
Phosphorylated eukaryotic initiation factor 2
Protein kinase B
Phospholipid phosphatase 3
Phosphorylated mechanistic target of rapamycin kinase.
RNA polymerase II subunit I
RNA polymerase III subunit H
Protein phosphatase 2
Peroxisome proliferator activated receptor alpha
Peroxisome proliferator sctivated receptor gamma
Protein phosphatase 2 regulatory subunit B″ gamma
Phosphorylated ribosomal protein S6
Phosphorylated ribosomal protein S6 kinase beta-1
Glycogen phosphorylase L
Quantitative polymerase chain reaction
Quadropule ion trap
REC8 meiotic recombination protein
Ribosomal protein S6
Ribosomal protein S9
Real-time polymerase chain reaction
Ribosomal S6 kinase
Ribosomal protein S6 kinase beta-1
SREBP cleavage-activating protein
Solute carrier family 44 member 4
Solute carrier family 51 subunit alpha
Sterol-regulated element-binding protein
Sterol regulatory element binding transcription factor 1c
Tricarboxylic acid cycle
Tumor necrosis factor
tRNA methyltransferase 61A
Ubiquitously expressed prefoldin like chaperone
Very low density lipoproteins
Symonds ME, Sebert SP, Budge H. Nutritional regulation of fetal growth and implications for productive life in ruminants. Animal. 2010;4:1075–83.
Bhutta ZA, Das JK, Rizvi A, Gaffey MF, Walker N, Horton S, et al. Evidence-based interventions for improvement of maternal and child nutrition: what can be done and at what cost? Lancet. 2013;382:452–77.
Ling T, Hernandez-Jover M, Sordillo LM, Abuelo A. Maternal late-gestation metabolic stress is associated with changes in immune and metabolic responses of dairy calves. J Dairy Sci. 2018;101:6568–80.
Hill TM, Bateman HG, Aldrich JM, Schlotterbeck RL, Tanan KG. Optimal concentrations of lysine, methionine, and threonine in milk replacers for calves less than five weeks of age. J Dairy Sci. 2008;91:2433–42.
Schwab CG, Broderick GA. A 100-Year Review: Protein and amino acid nutrition in dairy cows. J Dairy Sci. 2017;100:10094–112.
Kalhan SC, Marczewski SE. Methionine, homocysteine, one carbon metabolism and fetal growth. Rev Endocr Metab Disord. 2012;13:109–19.
McKee SE, Reyes TM. Effect of supplementation with methyl-donor nutrients on neurodevelopment and cognition: considerations for future research. Nutr Rev. 2018;76:497–511.
Atmaca G. Antioxidant effects of sulfur-containing amino acids. Yonsei Med J. 2004;45:776–88.
Coleman DN, Lopreiato V, Alharthi A, Loor JJ. Amino acids and the regulation of oxidative stress and immune function in dairy cattle. J Anim Sci. 2020;98(Supplement_1):S175–93.
Jacometo CB, Zhou Z, Luchini D, Corrêa MN, Loor JJ. Maternal supplementation with rumen-protected methionine increases prepartal plasma methionine concentration and alters hepatic mRNA abundance of 1-carbon, methionine, and transsulfuration pathways in neonatal Holstein calves. J Dairy Sci. 2017;100:3209–19.
Jin Z, Liu Y. DNA methylation in human diseases. Genes Dis. 2018;5:1–8.
Reynolds LP, Borowicz PP, Caton JS, Crouse MS, Dahlen CR, Ward AK. Developmental programming of fetal growth and development. Vet Clin Food Anim Pract. 2019;35:229–47.
Cai D, Jia Y, Song H, Sui S, Lu J, Jiang Z, et al. Betaine supplementation in maternal diet modulates the epigenetic regulation of hepatic gluconeogenic genes in neonatal piglets. PLoS One. 2014;9:e105504.
O’Neill RJ, Vrana PB, Rosenfeld CS. Maternal methyl supplemented diets and effects on offspring health. Front Genet. 2014;5. https://doi.org/10.3389/fgene.2014.00289.
Sookoian S, Gianotti TF, Burgueño AL, Pirola CJ. Fetal metabolic programming and epigenetic modifications: a systems biology approach. Pediatr Res. 2013;73:531–42.
Zhu Z, Cao F, Li X. Epigenetic programming and fetal metabolic programming. Front Endocrinol. 2019;10. https://doi.org/10.3389/fendo.2019.00764.
Barker DJ. In utero programming of chronic disease. Clin Sci. 1998;95:115–28.
Barker DJ. In utero programming of cardiovascular disease. Theriogenology. 2000;53:555–74.
Godfrey KM, Barker DJ. Fetal programming and adult health. Public Health Nutr. 2001;4:611–24.
Papas AM, Vicini JL, Clark JH, Peirce-Sandner S. Effect of rumen-protected methionine on plasma free amino acids and production by dairy cows. J Nutr. 1984;114:2221–7.
Wang C, Liu HY, Wang YM, Yang ZQ, Liu JX, Wu YM, et al. Effects of dietary supplementation of methionine and lysine on milk production and nitrogen utilization in dairy cows. J Dairy Sci. 2010;93:3661–70.
Coleman DN, Alharthi AS, Liang Y, Lopes MG, Lopreiato V, Vailati-Riboni M, et al. Multifaceted role of one-carbon metabolism on immunometabolic control and growth during pregnancy, lactation and the neonatal period in dairy cattle. J Anim Sci Biotechnol. 2021;12:27.
Zhou Z, Vailati-Riboni M, Trevisi E, Drackley JK, Luchini DN, Loor JJ. Better postpartal performance in dairy cows supplemented with rumen-protected methionine compared with choline during the peripartal period. J Dairy Sci. 2016;99:8716–32.
Batistel F, Arroyo JM, Bellingeri A, Wang L, Saremi B, Parys C, et al. Ethyl-cellulose rumen-protected methionine enhances performance during the periparturient period and early lactation in Holstein dairy cows. J Dairy Sci. 2017;100:7455–67.
Batistel F, Alharthi AS, Wang L, Parys C, Pan Y-X, Cardoso FC, et al. Placentome nutrient transporters and mammalian target of rapamycin signaling proteins are altered by the methionine supply during late gestation in dairy cows and are associated with newborn birth weight. J Nutr. 2017;147:1640–7.
Alharthi AS, Coleman DN, Liang Y, Batistel F, Elolimy AA, Yambao RC, et al. Hepatic 1-carbon metabolism enzyme activity, intermediate metabolites, and growth in neonatal Holstein dairy calves are altered by maternal supply of methionine during late pregnancy. J Dairy Sci. 2019;102:10291–303.
Alharthi AS, Batistel F, Abdelmegeid MK, Lascano G, Parys C, Helmbrecht A, et al. Maternal supply of methionine during late-pregnancy enhances rate of Holstein calf development in utero and postnatal growth to a greater extent than colostrum source. J Anim Sci Biotechnol. 2018;9:83.
Li P, Yin Y-L, Li D, Kim SW, Wu G. Amino acids and immune function. Br J Nutr. 2007;98:237–52.
Lopes MG, Alharthi AS, Lopreiato V, Abdel-Hamied E, Liang Y, Coleman DN, et al. Maternal body condition influences neonatal calf whole-blood innate immune molecular responses to ex vivo lipopolysaccharide challenge. J Dairy Sci. 2021;104:2266–79.
Jacometo CB, Zhou Z, Luchini D, Trevisi E, Corrêa MN, Loor JJ. Maternal rumen-protected methionine supplementation and its effect on blood and liver biomarkers of energy metabolism, inflammation, and oxidative stress in neonatal Holstein calves. J Dairy Sci. 2016;99:6753–63.
Jacometo CB, Alharthi AS, Zhou Z, Luchini D, Loor JJ. Maternal supply of methionine during late pregnancy is associated with changes in immune function and abundance of microRNA and mRNA in Holstein calf polymorphonuclear leukocytes. J Dairy Sci. 2018;101:8146–58.
Peñagaricano F, Souza AH, Carvalho PD, Driver AM, Gambra R, Kropp J, et al. Effect of maternal methionine supplementation on the transcriptome of bovine preimplantation embryos. PLoS One. 2013;8:e72302.
Hulbert LE, Moisá SJ. Stress, immunity, and the management of calves. J Dairy Sci. 2016;99:3199–216.
Zhou Z, Bulgari O, Vailati-Riboni M, Trevisi E, Ballou MA, Cardoso FC, et al. Rumen-protected methionine compared with rumen-protected choline improves immunometabolic status in dairy cows during the peripartal period. J Dairy Sci. 2016;99:8956–69.
Osorio JS, Trevisi E, Ji P, Drackley JK, Luchini D, Bertoni G, et al. Biomarkers of inflammation, metabolism, and oxidative stress in blood, liver, and milk reveal a better immunometabolic status in peripartal cows supplemented with Smartamine M or MetaSmart. J Dairy Sci. 2014;97:7437–50.
Mato JM, Martínez-Chantar ML. Lu SC. Methionine metabolism and liver disease. Annu Rev Nutr. 2008;28:273–93.
Lu SC. S-Adenosylmethionine. Int J Biochem Cell Biol. 2000;32:391–5.
Maddocks ODK, Labuschagne CF, Adams PD, Vousden KH. Serine metabolism supports the methionine cycle and DNA/RNA methylation through de novo ATP synthesis in cancer cells. Mol Cell. 2016;61:210–21.
Ducker GS, Rabinowitz JD. One-carbon metabolism in health and disease. Cell Metab. 2017;25:27–42.
Mosharov E, Cranford MR, Banerjee R. The quantitatively important relationship between homocysteine metabolism and glutathione synthesis by the transsulfuration pathway and its regulation by redox changes. Biochemistry. 2000;39:13005–11.
Sanderson SM, Gao X, Dai Z, Locasale JW. Methionine metabolism in health and cancer: a nexus of diet and precision medicine. Nat Rev Cancer. 2019;19:625–37.
Niculescu MD, Zeisel SH. Diet, methyl donors and DNA methylation: interactions between dietary folate, methionine and choline. J Nutr. 2002;132(8 Suppl):2333S–5S.
Sinclair KD, Allegrucci C, Singh R, Gardner DS, Sebastian S, Bispham J, et al. DNA methylation, insulin resistance, and blood pressure in offspring determined by maternal periconceptional B vitamin and methionine status. Proc Natl Acad Sci USA. 2007;104:19351–6.
Serefidou M, Venkatasubramani AV, Imhof A. The impact of one carbon metabolism on histone methylation. Front Genet. 2019;10. https://doi.org/10.3389/fgene.2019.00764.
Hayashi T, Teruya T, Chaleckis R, Morigasaki S, Yanagida M. S-adenosylmethionine synthetase is required for cell growth, maintenance of G0 phase, and termination of quiescence in fission yeast. iScience. 2018;5:38–51.
Kudo NR, Anger M, Peters AHFM, Stemmann O, Theussl H-C, Helmhart W, et al. Role of cleavage by separase of the Rec8 kleisin subunit of cohesin during mammalian meiosis I. J Cell Sci. 2009;122:2686–98.
Li J, Qian W-P, Sun Q-Y. Cyclins regulating oocyte meiotic cell cycle progression. Biol Reprod. 2019;101:878–81.
An Q, Sun H, Zhang J, Lu Z, Peng W, Xu S, et al. Methionine adenosyltransferase 2β participates in mouse oocyte maturation by regulating the MAPK pathway. Reprod Sci. 2020;27:163–71.
Kuroda A, Rauch TA, Todorov I, Ku HT, Al-Abdullah IH, Kandeel F, et al. Insulin Gene Expression Is Regulated by DNA Methylation. PLoS One. 2009;4. https://doi.org/10.1371/journal.pone.0006953.
Zachut M, Honig H, Striem S, Zick Y, Boura-Halfon S, Moallem U. Periparturient dairy cows do not exhibit hepatic insulin resistance, yet adipose-specific insulin resistance occurs in cows prone to high weight loss. J Dairy Sci. 2013;96:5656–69.
Rico JE, Myers WA, Laub DJ, Davis AN, Zeng Q, McFadden JW. Hot topic: Ceramide inhibits insulin sensitivity in primary bovine adipocytes. J Dairy Sci. 2018;101:3428–32.
Yabe D, Brown MS, Goldstein JL. Insig-2, a second endoplasmic reticulum protein that binds SCAP and blocks export of sterol regulatory element-binding proteins. Proc Natl Acad Sci U S A. 2002;99:12753–8.
Zhang HH, Lipovsky AI, Dibble CC, Sahin M, Manning BD. S6K1 regulates GSK3 under conditions of mTOR-dependent feedback inhibition of Akt. Mol Cell. 2006;24:185–97.
Shimura T, Kakuda S, Ochiai Y, Nakagawa H, Kuwahara Y, Takai Y, et al. Acquired radioresistance of human tumor cells by DNA-PK/AKT/GSK3β-mediated cyclin D1 overexpression. Oncogene. 2010;29:4826–37.
Arai T, Kano F, Murata M. Translocation of forkhead box O1 to the nuclear periphery induces histone modifications that regulate transcriptional repression of PCK1 in HepG2 cells. Genes Cells. 2015;20:340–57.
Girard J. Metabolic adaptations to change of nutrition at birth. Biol Neonate. 1990;58(Suppl 1):3–15.
Hammon HM, Steinhoff-Wagner J, Schönhusen U, Metges CC, Blum JW. Energy metabolism in the newborn farm animal with emphasis on the calf: endocrine changes and responses to milk-born and systemic hormones. Domestic Anim Endocrinol. 2012;43:171–85.
Steinhoff-Wagner J, Görs S, Junghans P, Bruckmaier RM, Kanitz E, Metges CC, et al. Maturation of endogenous glucose production in preterm and term calves1. J Dairy Sci. 2011;94:5111–23.
van Schaftingen E, Gerin I. The glucose-6-phosphatase system. Biochem J. 2002;362(Pt 3):513–32.
Gatfield D, Schibler U. Circadian glucose homeostasis requires compensatory interference between brain and liver clocks. Proc Natl Acad Sci U S A. 2008;105:14753–4.
Lee J-M, Seo W-Y, Song K-H, Chanda D, Kim YD, Kim D-K, et al. AMPK-dependent repression of hepatic gluconeogenesis via disruption of CREB.CRTC2 complex by orphan nuclear receptor small heterodimer partner. J Biol Chem. 2010;285:32182–91.
Zhang X, Yang S, Chen J, Su Z. Unraveling the regulation of hepatic gluconeogenesis. Front Endocrinol. 2019;9. https://doi.org/10.3389/fendo.2018.00802.
Wang Y, Inoue H, Ravnskjaer K, Viste K, Miller N, Liu Y, et al. Targeted disruption of the CREB coactivator Crtc2 increases insulin sensitivity. Proc Natl Acad Sci U S A. 2010;107:3087–92.
Govindan MV, Seguin C. Final steps in the feedback regulation of human glucocorticoid receptor gene and role of nuclear protein phosphatase 2A. Mol Biol Genet Eng. 2013;1:2.
Mueller KM, Themanns M, Friedbichler K, Kornfeld J-W, Esterbauer H, Tuckermann JP, et al. Hepatic growth hormone and glucocorticoid receptor signaling in body growth, steatosis and metabolic liver cancer development. Mol Cell Endocrinol. 2012;361:1–11.
Schlegel G, Keller J, Hirche F, Geißler S, Schwarz FJ, Ringseis R, et al. Expression of genes involved in hepatic carnitine synthesis and uptake in dairy cows in the transition period and at different stages of lactation. BMC Vet Res. 2012;8:28.
Jump DB, Botolin D, Wang Y, Xu J, Christian B, Demeure O. Fatty acid regulation of hepatic gene transcription. J Nutr. 2005;135:2503–6.
Rosenbaum S, Ringseis R, Most E, Hillen S, Becker S, Erhardt G, et al. Genes involved in carnitine synthesis and carnitine uptake are up-regulated in the liver of sows during lactation. Acta Vet Scand. 2013;55:24.
Karlic H, Lohninger S, Koeck T, Lohninger A. Dietary l-carnitine stimulates carnitine acyltransferases in the liver of aged rats. J Histochem Cytochem. 2002;50:205–12.
Longo N, Frigeni M, Pasquali M. Carnitine transport and fatty acid oxidation. Biochim Biophys Acta (BBA). 2016;1863:2422–35.
Zhang L, Joshi AK, Smith S. Cloning, expression, characterization, and interaction of two components of a human mitochondrial fatty acid synthase. Malonyltransferase and acyl carrier protein. J Biol Chem. 2003;278:40067–74.
Kremer L, Nampoothiri KM, Lesjean S, Dover LG, Graham S, Betts J, et al. Biochemical characterization of acyl carrier protein (AcpM) and malonyl-CoA:AcpM transacylase (mtFabD), two major components of Mycobacterium tuberculosis fatty acid synthase II. J Biol Chem. 2001;276:27967–74.
Hodson L. Hepatic fatty acid synthesis and partitioning: the effect of metabolic and nutritional state. Proceed Nutr Soc. 2019;78:126–34.
Foster DW. Malonyl-CoA: the regulator of fatty acid synthesis and oxidation. J Clin Invest. 2012;122:1958–9.
Auboiron S, Durand D, Bauchart D, Robert JC, Chapman MJ. Lipoprotein metabolism in the preruminant calf: effect of a high fat diet supplemented with L-methionine. J Dairy Sci. 1994;77:1870–81.
García-Trevijano ER, Latasa MU, Carretero MV, Berasain C, Mato JM, Avila MA. S-adenosylmethionine regulates MAT1A and MAT2A gene expression in cultured rat hepatocytes: a new role for S-adenosylmethionine in the maintenance of the differentiated status of the liver. FASEB J. 2000;14:2511–8.
Martinov MV, Vitvitsky VM, Banerjee R, Ataullakhanov FI. The logic of the hepatic methionine metabolic cycle. Biochim Biophys Acta (BBA). 2010;1804:89–96.
Fernández-Arroyo S, Cuyàs E, Bosch-Barrera J, Alarcón T, Joven J, Menendez JA. Activation of the methylation cycle in cells reprogrammed into a stem cell-like state. Oncoscience. 2016;2:958–67.
Loor JJ, Bionaz M, Drackley JK. Systems physiology in dairy cattle: nutritional genomics and beyond. Annu Rev Anim Biosci. 2013;1:365–92.
Jia H, Li X, Liu G, Loor JJ, Bucktrout R, Sun X, et al. Perilipin 5 promotes hepatic steatosis in dairy cows through increasing lipid synthesis and decreasing very low density lipoprotein assembly. J Dairy Sci. 2019;102:833–45.
Pinotti L, Baldi A, Dell’Orto V. Comparative mammalian choline metabolism with emphasis on the high-yielding dairy cow. Nutr Res Rev. 2002;15:315–32.
Schaeffer L, Gohlke H, Müller M, Heid IM, Palmer LJ, Kompauer I, et al. Common genetic variants of the FADS1 FADS2 gene cluster and their reconstructed haplotypes are associated with the fatty acid composition in phospholipids. Hum Mol Genet. 2006;15:1745–56.
Dushianthan A, Cusack R, Grocott MPW, Postle AD. Abnormal liver phosphatidylcholine synthesis revealed in patients with acute respiratory distress syndrome. J Lipid Res. 2018;59:1034–45.
Ridgway ND, Vance DE. Specificity of rat hepatic phosphatidylethanolamine N-methyltransferase for molecular species of diacyl phosphatidylethanolamine. J Biol Chem. 1988;263:16856–63.
DeLong CJ, Shen YJ, Thomas MJ, Cui Z. Molecular distinction of phosphatidylcholine synthesis between the CDP-choline pathway and phosphatidylethanolamine methylation pathway. J Biol Chem. 1999;274:29683–8.
Devlin AM, Singh R, Wade RE, Innis SM, Bottiglieri T, Lentz SR. Hypermethylation of Fads2 and altered hepatic fatty acid and phospholipid metabolism in mice with hyperhomocysteinemia. J Biol Chem. 2007;282:37082–90.
Nair KS, Short KR. Hormonal and signaling role of branched-chain amino acids. J Nutr. 2005;135:1547S–52S.
Jose DG, Good RA. Quantitative effects of nutritional essential amino acid deficiency upon immune responses to tumors in mice. J Exp Med. 1973;137:1–9.
Petro TM, Bhattacharjee JK. Effect of dietary essential amino acid limitations upon the susceptibility to Salmonella typhimurium and the effect upon humoral and cellular immune responses in mice. Infect Immun. 1981;32:251–9.
Tsukishiro T, Shimizu Y, Higuchi K, Watanabe A. Effect of branched-chain amino acids on the composition and cytolytic activity of liver-associated lymphocytes in rats. J Gastroenterol Hepatol. 2000;15:849–59.
Bassit RA, Sawada LA, Bacurau RFP, Navarro F, Martins E, Santos RVT, et al. Branched-chain amino acid supplementation and the immune response of long-distance athletes. Nutrition. 2002;18:376–9.
Dong X, Zhou Z, Saremi B, Helmbrecht A, Wang Z, Loor JJ. Varying the ratio of Lys:Met while maintaining the ratios of Thr:Phe, Lys:Thr, Lys:His, and Lys:Val alters mammary cellular metabolites, mammalian target of rapamycin signaling, and gene transcription. J Dairy Sci. 2018;101:1708–18.
Nan X, Bu D, Li X, Wang J, Wei H, Hu H, et al. Ratio of lysine to methionine alters expression of genes involved in milk protein transcription and translation and mTOR phosphorylation in bovine mammary cells. Physiol Genomics. 2014;46:268–75.
Ma YF, Batistel F, Xu TL, Han LQ, Bucktrout R, Liang Y, et al. Phosphorylation of AKT serine/threonine kinase and abundance of milk protein synthesis gene networks in mammary tissue in response to supply of methionine in periparturient Holstein cows. J Dairy Sci. 2019;102:4264–74.
Graber TG, Borack MS, Reidy PT, Volpi E, Rasmussen BB. Essential amino acid ingestion alters expression of genes associated with amino acid sensing, transport, and mTORC1 regulation in human skeletal muscle. Nutr Metab. 2017;14:35.
Kim J, Guan K-L. mTOR as a central hub of nutrient signalling and cell growth. Nat Cell Biol. 2019;21:63–71.
Dufner A, Thomas G. Ribosomal S6 kinase signaling and the control of translation. Exp Cell Res. 1999;253:100–9.
Showkat M, Beigh MA, Andrabi KI. mTOR signaling in protein translation regulation: implications in cancer genesis and therapeutic interventions. Mol Biol Int. 2014:686984.
Brosnan JT, Brosnan ME. The sulfur-containing amino acids: an Overview. J Nutr. 2006;136:1636S–40S.
Sun F, Cao Y, Cai C, Li S, Yu C, Yao J. Regulation of nutritional metabolism in transition dairy cows: energy homeostasis and health in response to post-ruminal choline and methionine. PLOS ONE. 2016;11:e0160659.
Abdelmegeid MK, Vailati-Riboni M, Alharthi A, Batistel F, Loor JJ. Supplemental methionine, choline, or taurine alter in vitro gene network expression of polymorphonuclear leukocytes from neonatal Holstein calves. J Dairy Sci. 2017;100:3155–65.
Matés JM, Pérez-Gómez C, Núñez de Castro I. Antioxidant enzymes and human diseases. Clin Biochem. 1999;32:595–603.
Li S, Tan H-Y, Wang N, Zhang Z-J, Lao L, Wong C-W, et al. The role of oxidative stress and antioxidants in liver diseases. Int J Mol Sci. 2015;16:26087–124.
Hamelet J, Seltzer V, Petit E, Noll C, Andreau K, Delabar JM, et al. Cystathionine beta synthase deficiency induces catalase-mediated hydrogen peroxide detoxification in mice liver. Biochim Biophys Acta (BBA). 2008;1782:482–8.
Lee H-O, Wang L, Kuo Y-M, Gupta S, Slifker MJ, Li Y, et al. Lack of global epigenetic methylation defects in CBS deficient mice. J Inherit Metab Dis. 2017;40:113–20.
Zhou Z, Vailati-Riboni M, Luchini DN, Loor JJ. Methionine and choline supply during the periparturient period alter plasma amino acid and one-carbon metabolism profiles to various extents: potential role in hepatic metabolism and antioxidant status. Nutrients. 2016:9.
Batistel F, Arroyo JM, Garces CIM, Trevisi E, Parys C, Ballou MA, et al. Ethyl-cellulose rumen-protected methionine alleviates inflammation and oxidative stress and improves neutrophil function during the periparturient period and early lactation in Holstein dairy cows. J Dairy Sci. 2018;101:480–90.
Ochoa JJ, Ramirez-Tortosa MC, Quiles JL, Palomino N, Robles R, Mataix J, et al. Oxidative stress in erythrocytes from premature and full-term infants during their first 72 h of life. Free Radic Res. 2003;37:317–22.
Gaál T, Ribiczeyné-Szabó P, Stadler K, Jakus J, Reiczigel J, Kövér P, et al. Free radicals, lipid peroxidation and the antioxidant system in the blood of cows and newborn calves around calving. Comp Biochem Physiol Part B. 2006;143:391–6.
Chen Z, Tian R, She Z, Cai J, Li H. Role of oxidative stress in the pathogenesis of nonalcoholic fatty liver disease. Free Rad Biol Med. 2020;152:116–41.
Bieging KT, Mello SS, Attardi LD. Unravelling mechanisms of p53-mediated tumour suppression. Nat Rev Cancer. 2014;14:359–70.
Johnson TM, Yu ZX, Ferrans VJ, Lowenstein RA, Finkel T. Reactive oxygen species are downstream mediators of p53-dependent apoptosis. Proc Natl Acad Sci U S A. 1996;93:11848–52.
Kim KS, Doss HM, Kim H-J, Yang H-I. Taurine stimulates thermoregulatory genes in brown fat tissue and muscle without an influence on inguinal white fat tissue in a high-fat fiet-induced obese mouse model. Foods. 2020;9. https://doi.org/10.3390/foods9060688.
Abumrad NA. The liver as a hub in thermogenesis. Cell Metab. 2017;26:454–5.
Vermorel M, Dardillat C, Vernet J, Null S, Demigne C. Energy metabolism and thermoregulation in the newborn calf. Ann Rech Vet. 1983;14:382–9.
Chase CCL, Hurley DJ, Reber AJ. Neonatal immune development in the calf and its impact on vaccine response. Vet Clin North Am Food Anim Pract. 2008;24:87–104.
Mugnier A, Pecceu K, Schelcher F, Corbiere F. A parallel evaluation of 5 indirect cost-effective methods for assessing failure of passive immunity transfer in neonatal calves. JDS Commun. 2020;1:10–4.
Tiegs G, Lohse AW. Immune tolerance: What is unique about the liver. J Autoimmunity. 2010;34:1–6.
Kubes P, Jenne C. Immune Responses in the Liver. Annu Rev Immunol. 2018;36:247–77.
Kingeter LM, Lin X. C-type lectin receptor-induced NF-κB activation in innate immune and inflammatory responses. Cell Mol Immunol. 2012;9:105–12.
Tian Z, Chen Y, Gao B. Natural killer cells in liver disease. Hepatology. 2013;57:1654–62.
del Fresno C, Iborra S, Saz-Leal P, Martínez-López M, Sancho D. Flexible signaling of myeloid C-type lectin receptors in immunity and inflammation. Front Immunol. 2018;9. https://doi.org/10.3389/fimmu.2018.00804.
Li Z, Diehl AM. Innate immunity in the liver. Curr Opin Gastroenterol. 2003;19:565–71.
Mikulak J, Bruni E, Oriolo F, Di Vito C, Mavilio D. Hepatic natural killer cells: organ-dpecific sentinels of liver Immune homeostasis and physiopathology. Front Immunol. 2019;10. https://doi.org/10.3389/fimmu.2019.00946.
Vivier E, Raulet DH, Moretta A, Caligiuri MA, Zitvogel L, Lanier LL, et al. Innate or adaptive immunity? The example of natural killer cells. Science. 2011;331:44–9.
Crome SQ, Lang PA, Lang KS, Ohashi PS. Natural killer cells regulate diverse T cell responses. Trends Immunol. 2013;34:342–9.
Peng H, Wisse E, Tian Z. Liver natural killer cells: subsets and roles in liver immunity. Cell Mol Immunol. 2016;13:328–36.
Li T-H, Liu L, Hou Y-Y, Shen S-N, Wang T-T. C-type lectin receptor-mediated immune recognition and response of the microbiota in the gut. Gastroenterol Rep (Oxf). 2019;7:312–21.
Geijtenbeek TBH, Gringhuis SI. Signalling through C-type lectin receptors: shaping immune responses. Nat Rev Immunol. 2009;9:465–79.
Brown GD, Willment JA, Whitehead L. C-type lectins in immunity and homeostasis. Nat Rev Immunol. 2018;18:374–89.
Pereira MS, Alves I, Vicente M, Campar A, Silva MC, Padrão NA, et al. Glycans as key checkpoints of T cell activity and function. Front Immunol. 2018;9. https://doi.org/10.3389/fimmu.2018.02754.
Taylor ME, Drickamer K, Schnaar RL, Etzler ME, Varki A. Discovery and classification of glycan-binding proteins. In: Varki A, Cummings RD, Esko JD, Stanley P, Hart GW, Aebi M, et al., editors. Essentials of Glycobiology. 3rd ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 2015. https://www.ncbi.nlm.nih.gov/books/NBK453061/. Accessed 12 Aug 2020.
Kagan JC, Iwasaki A. The phagosome as the organelle linking innate and adaptive immunity. Traffic. 2012;13:1053–61.
Lambert MP. Platelets in liver and renal disease. Hematol Am Soc Hematol Educ Program. 2016:251–5.
Diggs LP, Greten TF. The effects of platelet accumulation in fatty liver disease. Nat Rev Gastroenterol Hepatol. 2019;16:393–4.
Ramadori P, Klag T, Malek NP, Heikenwalder M. Platelets in chronic liver disease, from bench to bedside. JHEP Rep. 2019;1:448–59.
Uhlar CM, Burgess CJ, Sharp PM, Whitehead AS. Evolution of the serum amyloid A (SAA) protein superfamily. Genomics. 1994;19:228–35.
Meek RL, Eriksen N, Benditt EP. Murine serum amyloid A3 is a high density apolipoprotein and is secreted by macrophages. Proc Natl Acad Sci U S A. 1992;89:7949–52.
Uhlar CM, Whitehead AS. Serum amyloid A, the major vertebrate acute-phase reactant. Eur J Biochem. 1999;265:501–23.
Sander LE, Sackett SD, Dierssen U, Beraza N, Linke RP, Müller M, et al. Hepatic acute-phase proteins control innate immune responses during infection by promoting myeloid-derived suppressor cell function. J Exp Med. 2010;207:1453–64.
McDonald TL, Larson MA, Mack DR, Weber A. Elevated extrahepatic expression and secretion of mammary-associated serum amyloid A 3 (M-SAA3) into colostrum. Vet Immunol Immunopathol. 2001;83:203–11.
Keenan AB, Torre D, Lachmann A, Leong AK, Wojciechowicz ML, Utti V, et al. ChEA3: transcription factor enrichment analysis by orthogonal omics integration. Nucleic Acids Res. 2019;47:W212–24.
Bionaz M, Periasamy K, Rodriguez-Zas SL, Hurley WL, Loor JJ. A novel dynamic impact approach (DIA) for functional analysis of time-course omics studies: validation using the bovine mammary transcriptome. PLoS One. 2012;7:e32455.
Lee S, Dong HH. FoxO integration of insulin signaling with glucose and lipid metabolism. J Endocrinol. 2017;233:R67–79.
Niswender KD, Morrison CD, Clegg DJ, Olson R, Baskin DG, Myers MG, et al. Insulin activation of phosphatidylinositol 3-kinase in the hypothalamic arcuate nucleus: a key mediator of insulin-induced anorexia. Diabetes. 2003;52:227–31.
Altomonte J, Richter A, Harbaran S, Suriawinata J, Nakae J, Thung SN, et al. Inhibition of Foxo1 function is associated with improved fasting glycemia in diabetic mice. Am J Physiol-Endocrinol Metab. 2003;285:E718–28.
Jeong JY, Jeoung NH, Park K-G, Lee I-K. Transcriptional regulation of pyruvate dehydrogenase kinase. Diab Metab J. 2012;36:328–35.
Gopal K, Saleme B, Al Batran R, Aburasayn H, Eshreif A, Ho KL, et al. FoxO1 regulates myocardial glucose oxidation rates via transcriptional control of pyruvate dehydrogenase kinase 4 expression. Am J Physiol Heart Circ Physiol. 2017;313:H479–90.
Sparks JD, Dong HH. FoxO1 and hepatic lipid metabolism. Curr Opin Lipidol. 2009;20:217–26.
Cabrera-Ortega AA, Feinberg D, Liang Y, Rossa C, Graves DT. The role of forkhead box 1 (FOXO1) in the immune system: dendritic cells, T cells, B cells, and hematopoietic stem cells. Crit Rev Immunol. 2017;37:1–13.
John E, Wienecke-Baldacchino A, Liivrand M, Heinäniemi M, Carlberg C, Sinkkonen L. Dataset integration identifies transcriptional regulation of microRNA genes by PPARγ in differentiating mouse 3T3-L1 adipocytes. Nucleic Acids Res. 2012;40:4446–60.
Leonardi R, Zhang Y-M, Yun M-K, Zhou R, Zeng F-Y, Lin W, et al. Modulation of pantothenate kinase 3 activity by small molecules that interact with the substrate/allosteric regulatory domain. Chem Biol. 2010;17:892–902.
Leonardi R, Zhang Y-M, Rock CO, Jackowski S, Coenzyme A. Back in action. Progress Lipid Res. 2005;44:125–53.
Giralt A, Denechaud P-D, Lopez-Mejia IC, Delacuisine B, Blanchet E, Bonner C, et al. E2F1 promotes hepatic gluconeogenesis and contributes to hyperglycemia during diabetes. Mol Metab. 2018;11:104–12.
Warg LA, Oakes JL, Burton R, Neidermyer AJ, Rutledge HR, Groshong S, et al. The role of the E2F1 transcription factor in the innate immune response to systemic LPS. Am J Physiol Lung Cell Mol Physiol. 2012;303:L391–400.
Herzig S, Hedrick S, Morantte I, Koo S-H, Galimi F, Montminy M. CREB controls hepatic lipid metabolism through nuclear hormone receptor PPAR-γ. Nature. 2003;426:190–3.
Wen AY, Sakamoto KM, Miller LS. The role of the transcription factor CREB in immune function. J Immunol. 2010;185:6413–9.
Osorio JS, Ji P, Drackley JK, Luchini D, Loor JJ. Supplemental Smartamine M or MetaSmart during the transition period benefits postpartal cow performance and blood neutrophil function. J Dairy Sci. 2013;96:6248–63.
Schwab CG. Protected proteins and amino acids for ruminants. In: Biotechnology in Animal Feeds and Animal Feeding. John Wiley & Sons, Ltd; 1995. p. 115–41. doi:https://doi.org/10.1002/9783527615353.ch7.
Overton TR, LaCount DW, Cicela TM, Clark JH. Evaluation of a ruminally protected methionine product for lactating dairy cows. J Dairy Sci. 1996;79:631–8.
Batistel F, Alharthi AS, Yambao RRC, Elolimy AA, Pan Y-X, Parys C, et al. Methionine supply during late-gestation triggers offspring sex-specific divergent changes in metabolic and epigenetic signatures in bovine placenta. J Nutr. 2019;149:6–17.
Xu T, Alharthi ASM, Batistel F, Helmbrecht A, Parys C, Trevisi E, et al. Hepatic phosphorylation status of serine/threonine kinase 1, mammalian target of rapamycin signaling proteins, and growth rate in Holstein heifer calves in response to maternal supply of methionine. J Dairy Sci. 2018;101:8476–91.
Wu H, Southam AD, Hines A, Viant MR. High-throughput tissue extraction protocol for NMR- and MS-based metabolomics. Anal Biochem. 2008;372:204–12.
Zhou Z, Garrow TA, Dong X, Luchini DN, Loor JJ. Hepatic activity and transcription of betaine-homocysteine methyltransferase, methionine synthase, and cystathionine synthase in periparturient dairy cows are altered to different extents by supply of methionine and choline. J Nutr. 2017;147:11–9.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018. https://www.R-project.org/
Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11:R25.
Anders S, McCarthy DJ, Chen Y, Okoniewski M, Smyth GK, Huber W, et al. Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nat Protoc. 2013;8:1765–86.
Palombo V, Milanesi M, Sferra G, Capomaccio S, Sgorlon S, D’Andrea M. PANEV: an R package for a pathway-based network visualization. BMC Bioinform. 2020;21:1–7.
AA was supported by a fellowship from King Saud University (Riyadh, Saudi Arabia) to perform his PhD studies at the University of Illinois (Urbana). FB was supported by a fellowship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES; Brazil) and National Institute of Food and Agriculture (Grant: ILLU-538-914) to perform her PhD studies at University of Illinois (Urbana).
Evonik Nutrition & Care GmbH had a role in the study design and provided financial support to cover costs of animal use, data collection, and sample analysis.
Ethics approval and consent to participate
All procedures for this study were conducted in accordance with a protocol approved by the Institutional Animal Care and Use Committee (IACUC) of the University of Illinois (protocol # 14270).
Consent for publication
VP, MD, AA, FB, ET, and JJL declare that they have no competing interests. CP and JG are employees of Evonik Nutrition & Care GmbH. This does not alter the authors’ adherence to BMC Genomics policies on sharing data and materials.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of RNA sequencing and alignment for all the samples. mRNA libraries were sequenced on a HiSeq2500 (Illumina Inc.). Quality control metrics were performed on raw sequencing reads using the FASTQC (v0.11.15) application. An index of the reference genome was built and single-end clean reads for each individual were aligned to the reference genome by STAR (v2.5.3a). Reads were mapped and annotated to the Bos Taurus UMD_3.1.1, downloaded from Ensembl Genomes website. Reads aligned were quantified with Subread package (v1.5.2) based on the Refseq gene annotation.
List of all differentially expressed genes (FDR ≤ 0.10) from RNAseq data in liver tissue of 4-d old Holstein calves (n = 6/group) born to cows randomly assigned to receive a basal control (CON) diet from − 28 ± 2 d to parturition [1.47 Mcal/kg dry matter (DM) and 15.3% crude protein (CP)] with no added Met or CON plus ethyl cellulose Met (MET, Mepron®, Evonik Nutrition & Care GmbH, Germany).
List of all transcription factor (TF) obtained with ChEA3 transcription factor enrichment analysis tool and significantly overrepresented (FDR ≤ 0.05) in the differentially expressed gene (DEG) list from RNAseq data in liver tissue of 4-d old Holstein calves (n = 6/group) born to cows randomly assigned to receive a basal control (CON) diet from − 28 ± 2 d to parturition [1.47 Mcal/kg dry matter (DM) and 15.3% crude protein (CP)] with no added Met or CON plus ethyl cellulose Met (MET, Mepron®, Evonik Nutrition & Care GmbH, Germany). For each TF, the ranking (Rank), P value, false discovery rate (FDR), the overlapping genes among the DEG list, and the predicted state (i.e. impact and flux value) calculated using the rationale underlying the Dynamic Impact Approach (DIA) analysis are reported.
Extended discussion on methionine synthase (MTR) and cystathionine beta-synthase (CBS) activity, phosphatidylcholine (PC) and very low density lipoproteins (VLDL), and on single genes involved in immune status in MET calves.
KEGG pathways in ‘Lipid Metabolism’, ‘Amino Acid Metabolism’, and ‘Immune System’ subcategories resulting from the DIA analysis of RNAseq data from liver tissue of 4-d old Holstein calves (n = 6/group) born to cows randomly assigned to receive a basal control (CON) diet from − 28 ± 2 d to parturition [1.47 Mcal/kg dry matter (DM) and 15.3% crude protein (CP)] with no added Met or CON plus ethyl cellulose Met (MET, Mepron®, Evonik Nutrition & Care GmbH, Germany). The columns represent the effect (impact) and flux responses. The blue bars represent the effect value (0 to 300), and the flux columns represent negative (−) and positive (+) flux (− 300 to + 300) based on the direction of the effect. The positive flux (red bars) indicates an upregulation in treated (MET) liver tissue cells compared to control (CON) ones, while the negative flux (green bars) indicates a downregulation.
Ingredient and nutrient composition of far-off (from − 45 to − 29 d) and close-up (from − 28 d to parturition) diets fed to Holstein cows used in the present study.
Extended information about quantitative real-time PCR. Accession number, gene symbol, and forward and reverse primer sequences of genes analyzed in calf liver.
Gene/pathway-based network visualization obtained with Pathway Network Visualizer (PANEV) tool, considering the entire list of differentially expressed genes (FDR ≤ 0.10). The network was created taking into account the main pathways of interest involved in ‘Glycolysis/gluconeogenesis’, ‘Lipid metabolism’, ‘Amino acid metabolism’ and ‘Immune system’, along with the list of their interconnected (i.e. up- and downstream) pathways based on information retrieved on KEGG pathway database.
Transcription Factor/Gene/pathway-based network visualization matching the outcomes from Pathway Network Visualizer (PANEV) tool with the main results obtained by ChIP-X Enrichment Analysis 3 (ChEA3) software.
About this article
Cite this article
Palombo, V., Alharthi, A., Batistel, F. et al. Unique adaptations in neonatal hepatic transcriptome, nutrient signaling, and one-carbon metabolism in response to feeding ethyl cellulose rumen-protected methionine during late-gestation in Holstein cows. BMC Genomics 22, 280 (2021). https://doi.org/10.1186/s12864-021-07538-w
- Methyl donor
- Nutritional programming