Impact of post-ruminally infused macronutrients on bovine mammary gland expression of genes involved in fatty acid synthesis, energy metabolism, and protein synthesis measured in RNA isolated from milk fat
Journal of Animal Science and Biotechnology volume 11, Article number: 53 (2020)
Characterising the regulation of milk component synthesis in response to macronutrient supply is critical for understanding the implications of nutritional interventions on milk production. Gene expression in mammary gland secretory cells was measured using RNA isolated from milk fat globules from 6 Holstein-Friesian cows receiving 5-d abomasal infusions of saline, essential amino acids (AA), or glucose (GG) or palm olein (LG) without (LAA) or with (HAA) essential AA, according to a 6 × 6 Latin square design. RNA was isolated from milk fat samples collected on d 5 of infusion and subjected to real-time quantitative PCR. We hypothesised that mRNA expression of genes involved in de novo milk fatty acid (FA) synthesis would be differently affected by GG and LG, and that expression of genes regulating transfer of tricarboxylic acid cycle intermediates would increase at the HAA level. We also hypothesised that the HAA level would affect genes regulating endoplasmic reticulum (ER) homeostasis but would not affect genes related to the mechanistic target of rapamycin complex 1 (mTORC1) or the integrated stress response (ISR) network.
Infusion of GG did not affect de novo milk FA yield but decreased expression of FA synthase (FASN). Infusion of LG decreased de novo FA yield and tended to decrease expression of acetyl-CoA carboxylase 1 (ACC1). The HAA level increased both de novo FA yield and expression of ACC1, and tended to decrease expression of mitochondrial phosphoenolpyruvate carboxykinase (PCK2). mRNA expression of mTORC1 signaling participants was not affected by GG, LG, or AA level. Expression of the ε subunit of the ISR constituent eukaryotic translation initiation factor 2B (EIF2B5) tended to increase at the HAA level, but only in the presence of LG. X-box binding protein 1 (XBP1) mRNA was activated in response to LG and the HAA level.
Results show that expression of genes involved in de novo FA synthesis responded to glucogenic, lipogenic, and aminogenic substrates, whereas genes regulating intermediate flux through the tricarboxylic acid cycle were not majorly affected. Results also suggest that after 5 d of AA supplementation, milk protein synthesis is supported by enhanced ER biogenesis instead of signaling through the mTORC1 or ISR networks.
Lactating mammary gland secretory cells respond to glucose, fatty acids (FA), and amino acids (AA) through regulated catabolic and anabolic pathways facilitating biosynthesis of milk lactose, fat, and protein. Lactose synthesis is the major anabolic fate of glucose in bovine mammary cells, while glucose oxidation generates energy for FA and protein biosynthesis [1, 2]. The glucose moiety of lactose arises entirely from intramammary glucose, but hexose phosphate intermediates can contribute to the galactose moiety of lactose [3, 4]. Carbon skeletons of catabolised intramammary AA contribute as intermediates in the tricarboxylic acid (TCA) cycle, with fates including non-essential AA (NEAA) synthesis, FA synthesis, hexoneogenesis, and oxidation [5, 6]. Acetate is the primary precursor for de novo FA synthesis in ruminant mammary glands, but presence of long-chain FA (LCFA) for mammary uptake can downregulate this pathway in favour of incorporation of preformed LCFA into milk fat [7, 8]. Intersections between glycolysis, hexoneogenesis, AA catabolism, protein synthesis, and lipogenesis allow mammary cells to be flexible in their use of available substrates in response to the nutritional status of the animal. However, expression of genes related to carbon transfer and energy generation in bovine mammary secretory cells have gone largely unexplored in vivo, particularly in response to changes in energy and protein supply.
Understanding regulation of milk protein synthesis is crucial when aiming to optimize absorptive protein supply and AA profile in dairy cattle diets to minimize excess AA catabolism and nitrogen (N) excretion [9,10,11,12]. Intracellular signaling cascades of the mechanistic target of rapamycin complex 1 (mTORC1) and the integrated stress response (ISR) have been implicated in regulation of protein synthesis in dairy cattle mammary glands during postruminal infusions of starch, glucose, casein, or AA mixtures lasting up to 36 h [13, 14]. However, several multi-day (≥ 5 d) studies with glucose and essential AA (EAA) infusions have revealed that signaling through these networks, assessed at the level of protein phosphorylation [10, 11, 15] or mRNA expression , does not entirely account for persistent milk protein production responses to extra energy and AA supply. Instead, studies of mRNA expression of constituent proteins of the unfolded-protein response (UPR) in mammary glands of dairy cattle suggest that endoplasmic reticulum (ER) biogenesis and secretory cell differentiation could be activated in response to 5-d increases in protein supply [16, 17]. The UPR defines a 3-arm signaling cascade that copes with perturbations in ER homeostasis by regulating protein synthesis rate, stimulating ER biogenesis, and initiating cell apoptosis .
We recently reported results on milk production and composition and mammary gland metabolite utilization from cows post-ruminally supplemented for 5 d with glucose or palm olein (major FA constituents are palmitic, oleic, and linoleic acid) at low and high metabolizable protein (MP) levels (created through supplementation of EAA in a casein profile) [19, 20]. Independent of MP level, glucose infusion stimulated body energy retention and reduced milk energy output, increased milk N efficiency, reduced arterial concentration of branched-chain AA, Arg, and Lys, and increased plasma flow to the mammary gland. Palm olein infusion increased milk energy output independent of MP level, did not affect body energy retention, plasma flow to the mammary gland, and milk N efficiency, and did not elicit any major effects on mammary gland AA metabolism. The high MP level increased milk energy output, did not affect body energy retention and plasma flow to the mammary gland, decreased milk N efficiency, and increased arterial concentration and mammary gland uptake of EAA. Considering these responses in milk production and mammary gland AA utilization in response to glucose, palm olein, and EAA, we investigated mammary gland expression of select genes involved in de novo FA synthesis, the TCA cycle, protein synthesis, and the UPR using RNA isolated from milk fat as a source of genetic material from mammary secretory cells. We hypothesised that after 5 d, mRNA expression of genes involved in de novo FA synthesis would be differently affected by glucose and palm olein, and that mRNA expression of proteins regulating transfer of TCA cycle intermediates would increase in response to the infused EAA. Furthermore, we hypothesised that EAA infusion would affect genes related to ER homeostasis but would not affect genes related to mTORC1 signaling or the ISR network.
Details of the experimental design, animal housing, ration composition and preparation, and feed chemical analyses have been described by Nichols et al. . Briefly, 6 rumen-fistulated, second lactation Holstein-Friesian dairy cows (97 ± 13 d in milk) were randomly assigned to a 6 × 6 Latin square design to test the effects of energy from glucogenic (GG) or lipogenic (LG) substrates at low (LAA) and high (HAA) EAA levels. Individual cows were housed in tie stalls within identical climate respiration chambers, and were allowed 5 d of adaptation to the housing environment before the first experimental period began. Each experimental period consisted of 5 d of continuous abomasal infusion followed by 2 d of no infusion. Cows were fed a total mixed ration (14% crude protein) consisting of 42% corn silage, 31% grass silage, and 27% concentrate on a dry matter basis, formulated to meet 100% and 83% of net energy (NEL) and MP requirements for lactation , respectively. Daily feed intake of individual cows during the entire experiment was restricted by 10% based on their daily individual ad libitum intake measured during an adaptation period. Fresh feed was allocated twice daily, with the exception of a 34-h window over d 4 and 5 of each period where feed was dispensed in equal portions every 2 h to promote metabolic steady-state conditions in preparation for a blood sampling protocol described by Nichols et al. . Cows were milked twice daily at 05:30 and 15:30 h and milk volume was recorded by weight at each milking time. Sampling and analytical procedures for milk fat, FA, protein, and lactose composition have been described by Nichols et al. .
Infusion lines were placed in the abomasum via the rumen cannula 7 d before the first experimental period and were checked daily for patency and position. The 6 abomasal infusion treatments were: 1) 0.9% saline (LAA-C; 90% NEL, 75% MP), 2) 1319 g/d glucose (LAA-GG; 100% NEL, 75% MP), 3) 676 g/d palm olein (LAA-LG; 100% NEL, 75% MP), 4) 844 g/d of a complete EAA mixture in the same profile and amount as found in casein (HAA-C; 100% NEL, 120% MP), 5) 1319 g/d glucose + 844 g/d EAA (HAA-GG; 110% NEL, 120% MP), and 6) 676 g/d palm olein + 844 g/d EAA (HAA-LG; 110% NEL, 120% MP), where MP and NEL reflect the proportion of requirements met by the restricted feeding level of the diet plus the infusate, and are expressed relative to target requirements for cows consuming 20 kg dry matter/d and producing 30 kg/d of milk containing 40 g/kg fat and 34 g/kg protein. Treatment solutions (prepared in 10-L batches) were infused via multi-channel peristaltic pumps at a rate of 6.95 mL/min continuously for 5 d. Daily infusion dosages (g/d) of LAA-GG, LAA-LG, and HAA-C were designed to be isoenergetic based on the gross energy content of the EAA infusion. Details of the composition of the EAA mixture and palm olein FA composition are described by Nichols et al. .
Milk fat collection, RNA extraction, and real-time quantitative PCR
Procedures for milk fat collection, RNA extraction, and real-time quantitative PCR were performed according to those described by Nichols et al. , with some modifications. During morning milking on d 5 of each period, 10 mL of milk was collected by hand from individual cows within 10 min from the time the milking machine was removed. This approach allowed each sample to be processed individually, thus minimizing the amount of time between alveolar ejection of fresh milk and placing the collected milk fat into RNase inhibitor (TRIzol Reagent; Invitrogen, ThermoFisher Scientific, Waltham, MA, USA). Samples were immediately centrifuged at 2000×g for 10 min at room temperature to facilitate separation of the fat fraction. This centrifugation process of whole milk is expected to render the contribution of other cellular material (immune cells, sloughed epithelial cells, stem cells) negligible in the milk fat layer [22, 23]. Approximately 1 g of the supernatant cream layer was transferred into 6 mL TRIzol Reagent, mixed vigorously, snap frozen in liquid N2, and stored at − 80 °C until RNA extraction. Total RNA was isolated from milk fat according to TRIzol manufacturer’s instructions for handling samples with a high fat content. Total RNA concentrations and purity were determined by optical density measurement using a Nano-Drop ND-1000 (ThermoFisher Scientific). The average 260/280 absorbance ratio of total RNA samples was adequate at 1.89 ± 0.117 (mean ± SD) and the average RNA concentration was 115 ± 64.1 ng/μL (mean ± SD). The mean RNA integrity number (RIN) was 4.5 ± 2.30 (mean ± SD). The low concentration of RNA found with this procedure may have impacted RIN measurement . This RIN is low by conventional standards , but was consistent across all samples, suggesting it was inherent to the sample type and not due to degradation of only a few samples. An aliquot of 100 ng total RNA was reverse transcribed in a total volume of 20 μL with Superscript III (100 U; ThermoFisher Scientific) in the presence of random hexamers (250 ng; Roche, Almere, the Netherlands), dithiothreitol (0.5 mmol/L; ThermoFisher Scientific) and dNTPs (0.5 mmol/L; Roche) at 25 °C for 5 min and 50 °C for 1 h. After inactivation of the enzyme (70 °C, 15 min), cDNA was stored at − 80 °C until further analysis.
Specific primers (Eurogentec, Maastricht, the Netherlands) were intron-spanning and designed to yield amplicons in the range of 100 and 180 bp (Table 1) with efficiencies of 90–100%. Templates were amplified after a preincubation of 10 min at 95 °C, followed by amplification for 40 cycles (10 s at 95 °C, 5 s at 60 °C, and 5 s at 72 °C) on a 7500 Fast Real-Time PCR System (Applied Biosystems Deutschland GmbH, Darmstadt, Germany) by using the SensiMix SYBR Low-ROX kit (Bioline UK Ltd., London, UK). All melting curves confirmed that a single amplicon was produced. As a control, a standard cDNA was included on each qPCR plate, resulting in a single and similar amplicon for the standard and cDNA samples as evidenced by melting curve analysis. As internal standards, expression of housekeeping genes glyceraldehyde 3-phosphate dehydrogenase (GAPDH), ribosomal protein S9 (RPS9), and ubiquitously expressed transcript isoform 2 (UXT) were analysed. NormFinder  identified RPS9 as the most stable gene across experimental period and treatment. Fold changes in gene expression relative to LAA-C were calculated by the 2−ΔΔCt method  using RPS9 expression as the reference gene.
Variances in milk and component production and gene expression data were analysed using the MIXED procedure of SAS (SAS Institute Inc., Cary, NC). The model contained main and interaction effects of infusion treatment factors (GG, LG, and AA level) and experimental period as fixed effects and cow as a random effect. No carryover effects on milk and milk component production were observed between periods, as assessed by testing for an effect of the previous treatment in the ANOVA. Differences were considered significant at P ≤ 0.05 and tendencies at 0.05 < P ≤ 0.10. Multiple comparisons between treatment means were made using the Tukey-Kramer method when GG × AA or LG × AA interactions were detected at P ≤ 0.10.
Production of milk protein, fat, lactose, and fatty acids
Infusing GG or LG at the HAA level did not affect yield of total milk or milk components differently than at the LAA level (no significant GG × AA or LG × AA interactions; P > 0.28; Table 2). Total milk yield was unaffected by GG and LG (P > 0.60) and increased at the HAA level (P < 0.01). Milk fat yield decreased in response to GG (P < 0.01), increased in response to LG (P < 0.01), and increased at the HAA level (P = 0.02). De novo FA (< 16 carbons) yield decreased in response to LG (P = 0.02), increased at the HAA level (P < 0.01), and was unaffected by GG (P = 0.88). Yield of preformed FA (> 16 carbons) decreased in response to GG (P < 0.01), increased in response to LG (P < 0.01), and was not affected by AA level (P = 0.25). Yield of mixed FA (16 carbons) was not affected by GG, LG, or AA level (P ≥ 0.14). Yields of milk protein and lactose were unaffected by GG and LG and increased at the HAA level (P < 0.01).
Expression of CSN2 and genes related to fatty acid synthesis and energy-yielding pathways
Infusing GG or LG at the HAA level did not affect the expression of CSN2 or genes associated with energy-yielding pathways differently than at the LAA level (no GG × AA or LG × AA interactions; P ≥ 0.17; Table 3). Infusion of GG decreased CSN2 and FASN expression 11% and 38%, respectively, compared with infusions with no GG (P = 0.02). Infusion of LG tended to decrease ACC1 expression (P = 0.08), and decreased CSN2 expression 22% compared with infusions with no LG (P < 0.01). The HAA level increased ACC1 expression 76% over infusions at the LAA level (P = 0.01), and tended to decrease PCK2 expression (P = 0.06).
Expression of genes related to the mechanistic target of rapamycin complex 1 and the integrated stress response network
Expression of mTORC1- and ISR-related genes are presented in Table 4. Expression of translation initiation factor EIF2B5 tended to be affected by a LG × AA interaction, where EIF2B5 expression was numerically higher (55%) at the HAA level than at the LAA level (P = 0.07). Gene expression of EIF2A, EIF2AK3, EIF4E, PRKAA1 and RPS6KB1 was not affected by treatment (P ≥ 0.12).
Expression of genes related to endoplasmic reticulum homeostasis and biogenesis, and cell differentiation and turnover
Genes related to ER homeostasis and biogenesis, and cell differentiation and turnover are presented in Table 5. mRNA expression of the serine/threonine kinase CDC42BPA was affected by a GG × AA interaction, where expression was 40% higher with GG infusion at the LAA level compared with the HAA level (P = 0.05). Infusion of GG tended to increase the expression of cell growth regulator MYC (P = 0.10), and tended to decrease expression of unspliced (u) ER stress-related transcription factor XBP1 (P = 0.08). Infusion of LG decreased the expression of XBP1u (P = 0.02) and increased the proportion of spliced (s) XBP1 to XBP1u (P = 0.04). The HAA level tended to increase MYC expression (P = 0.08), and increased the expression of XBP1s and the proportion of XBP1s to XBP1u (P ≤ 0.01).
Milk fat RNA technique
RNA isolated from milk fat was used as the source of genetic material from mammary secretory cells to measure gene expression in this study. A small fraction of cytoplasmic material from mammary secretory cells is entrapped during secretion of milk fat globules [29, 30]. Milk fat globules are an unconventional non-tissue source of extracellular mRNA , but representation of the mammary secretory cell transcriptome in the genomic material found in cytoplasmic crescents of milk fat globules has been validated in primates [23, 32] and ruminants [33, 34]. Cytoplasmic crescents in bovine milk are small and are present in a low proportion of milk fat globules . However, the relatively high fat content of ruminant milk allows the abundance of cytoplasmic material in milk fat to be adequate for RNA isolation [33, 35]. Compared with the common technique of mammary tissue biopsy, milk fat collection to obtain RNA from secretory cells for in vivo signaling analysis is noninvasive and enables a larger number of animals to be sampled at a higher frequency, and with relative ease compared to tissue biopsy. In this way, the technique allows mRNA expression of genes related to mammary gland metabolism to be characterized under a variety of experimental conditions. RNA isolated from milk fat has been previously used to assess the effects of nutritional intervention on mammary cell metabolic pathways in studies with dairy cattle [17, 35]. We employed this technique in the current study to assess changes in mRNA expression of a select set of genes related to FA synthesis, energy metabolism, and protein synthesis in response to increased absorptive supply of glucose, palm olein, and EAA. The genes of interest were selected based on results of previous investigations in response to energy and protein supply using RNA isolated from biopsied mammary tissue  and milk fat globules , and identification of their associated pathways via microarray analysis .
Fatty acid synthesis
Acetyl-CoA carboxylase 1 (ACC1), encoded by ACC1, and fatty acid synthase (FAS), encoded by FASN, are key participants in de novo biosynthesis of milk lipids in bovine mammary cells. The first step of de novo FA synthesis, the carboxylation of acetyl-CoA to form malonyl-CoA, is catalysed by ACC1 . Next, FAS catalyses the condensation of malonyl-CoA and acetyl-CoA and forms FA of ≤16 carbons . Coordinated regulation of these enzymes in dairy cattle mammary glands has been predominantly investigated under conditions of milk fat depression, a phenotype characterised by reduced yield of de novo-synthesized FA (mainly < 16 carbon FA) and total milk fat [39, 40]. In the current study, GG infusion did not affect de novo FA yield or ACC1 expression, but decreased FASN expression by 38%. Considering the key role of FAS in de novo FA synthesis, reduced FASN expression does not align with our observation of no change in de novo FA yield. However, GG infusion decreased total milk fat yield by 16%. This was driven by the stimulatory effect of increased circulating insulin on lipogenesis in extra-mammary tissues, leading to the reduction of circulating LCFA (> 16 carbons) for incorporation into milk fat [19, 20]. Lipogenic genes in bovine mammary cells are under the control of transcription factor families, such as the sterol regulatory element binding proteins (SREBP) [41, 42]. FASN is a transcriptional target of SREBP1 , and both constituents are downregulated when total milk fat yield is decreased during diet-induced milk fat depression or feed restriction [35, 44]. The relatively lower abundance of LCFA in mammary cells for incorporation into milk triacyglycerides (TAG), as suggested by reduced arterial LCFA concentration during GG infusion [19, 20], may have resulted in a coordinated downregulation of the expression of rate-limiting enzymes in the glycerol 3-phosphate pathway, which are also controlled by SREBP . Taken together, regulation of FAS at the transcriptional level after 5 d of glucose infusion did not impact downstream changes required to affect de novo FA yield.
Infusion of LG decreased de novo FA yield by 8%, tended to decrease ACC1 expression, but did not affect FASN expression. High intracellular concentrations of LCFA or their acyl-CoA esters inhibit ACC1 in lipogenic tissues [37, 45]. This inhibitory effect could have arisen in the presence of LCFA from palm olein (predominantly C16:0 and C18:1), and may have had a relatively larger impact on ACC1 than FAS, thus impacting their transcripts accordingly. A similar pattern was observed in response to EAA infusion, where de novo milk FA yield and ACC1 expression both increased at the HAA level, but FASN expression was not affected. Overall, coordinated changes in mRNA expression of ACC1 and FASN were not observed in response to glucose, palm olein, or EAA infusions in this study. Our results suggest that when mammary cells are supplied with glucose, palm olein, or EAA for 5 d, a distinctly coupled increase or decrease in both enzymes at the transcriptional level is not required to produce phenotypic changes in the FA composition of milk fat. In line with this observation, previous studies have observed variations between milk FA synthesis and lipogenic gene expression. In lactating cows and goats, supplements of corn oil or sunflower oil plus wheat starch reduced milk FA content and yield, but little or no variation in mRNA levels of mammary lipogenic genes was observed after 25 or 27 d of feeding [46, 47]. In the study of Invernizzi et al. , mammary gland expression of ACC1 and FASN and de novo milk FA synthesis from cows was increased after 7 d of dietary supplementation with saturated fat but not fish oil. By d 21 of this supplementation, differences in gene expression between saturated fat and fish oil had disappeared but the differences in de novo milk FA synthesis remained.
Tricarboxylic acid cycle
Intermediates arising from AA catabolism in mammary cells contribute to TCA cycle anaplerosis [5, 6]. The mitochondrial isoforms of malic enzyme (ME), encoded by ME2, and phosphoenolpyruvate carboxykinase (PEPCK), encoded by PCK2, channel TCA cycle intermediates derived from NEAA and certain EAA (Ile, Thr, Val) to pyruvate [49, 50]. PEPCK catalyses the GTP-driven transformation of oxaloacetate to phosphoenolpyruvate (PEP) [49, 51]. Malic enzyme converts malate to pyruvate, independent of glycolytic flux [50, 52]. The mitochondrial forms of both enzymes support oxidative flux through pyruvate. Therefore, these channels are important in mammary secretory cells where glucose is used for biosynthesis and is not the main fuel for the TCA cycle , and they may be upregulated when the rate of intramammary AA catabolism increases in response to MP supplementation . In the current study, both the intramammary catabolism of Ile and Val and the calculated intramammary glucose deficit for lactose and fat synthesis increased in response to infused EAA at the high MP level . Therefore, we hypothesised that mRNA expression of PCK2 and ME2 would increase, reflecting increased flux of carbon skeletons of AA (AA-C) through the TCA cycle to spare glucose in support of the increase in lactose and fat yield during EAA infusion.
In contrast to our hypothesis, PCK2 expression tended to decrease at the HAA level, and ME2 expression was not affected. Alternative to sparing glucose from oxidation in the TCA cycle, intracellular AA-C could support lactose synthesis through contribution to galactose synthesis [6, 20, 54]. Upregulation of the cytosolic isoform of PEPCK would favour this process over the mitochondrial form by incorporating AA-C into cytosolic PEP . Cytosolic PEP can enter the pentose phosphate cycle, producing hexose phosphate intermediates contributing to galactose [3, 4], whereas AA-C incorporated into mitochondrial PEP favours NEAA synthesis and oxidation [1, 49]. Data on gene expression or enzyme activity of the cytosolic or mitochondrial isoforms of PEPCK and ME in bovine mammary glands in response to energy or protein supply are scarce. We recently reported that extra energy from protein or fat supplementation did not affect PCK2 expression in mammary cells, and protein supplementation decreased expression of ME2, but only in combination with low fat diets . However, such an interaction between aminogenic and lipogenic nutrient supply on ME2 did not occur in the present experiment. Further investigation at the transcriptional and translational level of enzymes involved in intramammary carbon transfer would allow a clearer picture of energy balance and AA-C metabolism in lactating bovine mammary glands in response to glucogenic, lipogenic, and aminogenic substrates.
Previous reports where expression of milk protein genes are not indicative of observed milk protein yields were substantiated by the current study [16, 17, 55]. Milk protein yield was not affected by GG or LG infusion, but expression of CSN2, encoding β-casein, decreased 11% and 22% in response to GG and LG, respectively. The HAA level increased milk protein yield 28% but did not affect CSN2 expression. Similarly, in a previous study, the same dose and profile of abomasally infused EAA increased milk protein yield 30% and did not affect mammary CSN2 expression .
As hypothesised, and in agreement with others [10, 11, 16], we observed minimal effects on mRNA expression of mTORC1 and ISR participants in response to the 5-d infusion of EAA, despite the 0.2 kg/d increase in milk protein yield at the HAA level. Expression of EIF2B5, encoding the ε subunit of eukaryotic translation initiation factor 2B (eIF2B), tended to increase at the HAA level, but more so in the presence of LG. eIF2B is a 5-subunit guanine nucleotide exchange factor responsible for the formation of active eIF2-GTP complexes that deliver initiator Met-tRNA to ribosomes during each round of translation initiation . To our knowledge, there are no previous reports investigating mRNA expression of eIF2B in response to high fat conditions in mammary cells. In heart muscle of rats, free FA increased the nucleotide exchange activity of eIF2B . Greater EIF2B5 expression has the potential to facilitate protein translation. Although this observation is in alignment with the increase in milk protein yield at the HAA level, EIF2B5 expression increased more at the HAA level in the presence of LG, while the response in milk protein yield was independent of LG. This, and our finding that mRNA expression of no other mTORC1 or ISR constituents were affected by GG, LG, or AA level, supports our hypothesis that other mechanisms may have greater influence on phenotypic milk protein yield responses after 5 d of nutrient infusion.
Previous investigation at the mRNA level from biopsied mammary tissue suggested that ER biogenesis in mammary cells is responsive to EAA supply, potentially through regulation on separate UPR signaling arms . Initiation of the UPR in mammalian cells occurs by dissociation of the ER chaperone binding immunoglobulin protein (BiP/GRP78) onto 3 ER transmembrane proteins, namely protein kinase R-like ER kinase (PERK), activating transcription factor 6 (ATF6), and inositol-requiring enzyme 1 (IRE1), which initiate their respective UPR arms . The elevated proportion of the active spliced XBP1 available for translation relative to the inactive unspliced XBP1 mRNA (spliced/unspliced) at the HAA level suggests induction of the IRE1-mediated arm of the UPR. Splicing of XBP1 mRNA is needed to generate the active transcription factor (XBP1s), which has been implicated in lactating bovine mammary glands in response to increased protein supply [16, 17], and at the onset of lactation . The translation product of XBP1s regulates genes involved in ER formation and expansion, and secretory vesicle maturation in a variety of differentiated secretory cells, such as antibody-secreting plasma cells , gastric zymogenic cells , pancreatic endocrine and exocrine cells [61, 62], and mammary secretory cells [63, 64].
The increased ratio of XBP1s to XBP1u in response to EAA infusion suggests potential stimulation of ER biogenesis in support of greater mammary cell protein secretory capacity, in line with our hypothesis. Similarly, Nichols et al.  observed an increased XBP1s to XBP1u ratio in response to 5-d abomasal EAA infusion in the same AA profile and dose, and Nichols et al.  observed a similar increase in response to rumen-protected protein supplementation for 27 d, but only in low-fat diets. The ratio of XBP1s to XBP1u also increased in response to LG infusion in the current study. Increased XBP1 splicing in response to high fat diets has been reported in murine liver and skeletal muscle [65, 66]. Moreover, there could be a link between increased XBP1 transcript splicing in response to LG and increased EIF2B5 expression at the HAA level in the presence of LG. The IRE-1-mediated arm of the UPR is protective in its effort to regulate protein synthetic capacity of cells by reducing ER protein load and inhibiting apoptosis . Likewise, activation of eIF2Bε upon toll-like receptor stimulation was protective against ER-stress-mediated cell death in macrophages . Together with the increased XBP1s to XBP1u ratio, the tendency for EIF2B5 expression to be increased relatively more at the HAA level in the presence of LG may arise from a response whereby FA promote adaptation in mammary secretory cells in support of persistent milk protein synthesis when EAA are in abundant supply. Overall, adaptive ER biogenesis in bovine mammary cells in response to EAA supplementation and the presence of LCFA, suggested here by the increased XBP1s to XBP1u ratio, is a potentially attractive explanation for observed increases in milk protein synthesis. However, continued investigation of mechanisms involving changes in rate of cell turnover, secretory cell differentiation, and secretory cell capacity under supplementation with energy and protein substrates is necessary to characterize the adaptive mechanisms of lactating mammary glands in response to nutritional interventions. In particular, the temporal characterisation of transcriptional and translational changes in participants of cell signaling networks, such as the link between mRNA expression of UPR genes and ER biogenesis, should be further elucidated in lactating bovine mammary glands from both conventional tissue samples and from application of the herein described milk fat RNA technique.
Altered secretory cell number in response to increased nutrient supply could be contributing to faster protein synthesis in response to EAA supplementation. To test this hypothesis, we analysed mRNA expression of select genes encoding proliferation and apoptosis pathways. Genes related to cytoskeletal dynamics and cell proliferation, CDC42BPA and MYC, were both affected by glucose and EAA infusion. Expression of CDC42BPA tended to increase 40% when GG was infused at the LAA level compared with HAA-GG. CDC42BPA encodes myotonic dystrophy kinase-related CDC42-binding kinase (MRCK) α, a downstream effector of the cell division control protein 42 homolog (CDC42) involved in cytoskeleton assembly and organization. This gene has been characterised in mammary cells primarily in view of breast cancer pathogenesis , and in cattle in association with cell reorganization during bacterial infection from Mycobacterium avium ssp. paratuberculosis . The fact that no GG × AA interaction was observed for milk protein yield suggests that the increase in CDC42BPA expression with GG infusion at the LAA level may be unrelated to synthesis of export proteins in mammary cells. The c-Myc transcription factor, encoded by MYC, regulates genes required for cell proliferation, growth, differentiation, and apoptosis, but is also associated with ribosome biogenesis and protein synthesis required for cell growth [70, 71]. Expression of MYC tended to increase independently in response to GG and the HAA level, but milk protein yield was only affected at the HAA level. Nichols et al.  observed increased mammary MYC expression but no extra stimulation of milk protein synthesis when glucose was added to EAA infusions. In murine pancreatic islet cells, mRNA expression of MYC was increased in vitro under supraphysiological glucose concentrations , and during hyperglycaemia in vivo . In hepatocytes, growth and proliferation in response to protein and energy restriction and re-feeding schemes altered expression of MYC as cells progressed through the development cycle [74,75,76]. In bovine mammary glands, MYC expression was associated with cell proliferation in response to milking frequency , and during the dry period and early lactation . Taken together, these findings suggest that mammary cell hyperplasia or hypertrophy could be affected by transcription factors in response to glucogenic and aminogenic substrates, but further study is required to determine the biological importance with respect to persistent milk component synthesis throughout lactation.
Expression of genes involved in de novo FA synthesis in mammary gland secretory cells is sensitive to 5-d infusions of glucose, palm olein, and EAA, but a coupled increase or decrease in ACC1 and FASN mRNA expression was not required to produce phenotypic changes in de novo FA yield. Reduced FASN expression in response to glucose infusion may have been related to downregulation of TAG synthesis. Possible inhibition of ACC1 by LCFA during palm olein infusion may have decreased ACC1 expression, in line with increased preformed FA in milk fat. Stimulation of de novo FA synthesis by EAA infusion was associated with increased ACC1 expression. Expression of PCK2 tended to decrease during EAA infusion, and ME2 expression was not affected, suggesting that increased intramammary AA catabolism did not impact these enzymes at the mRNA level, or that AA-C shuttling occurred through different channels of the TCA cycle. We found no effects of 5-d glucose, palm olein, or EAA infusions on mRNA expression of participants in the mTORC1 signaling network. Expression of the ISR constituent EIF2B5 tended to increase with EAA infusion, but only in the presence of palm olein, and this did not coincide with extra milk protein synthesis when EAA and palm olein were infused together. Instead, splicing of the XBP1 transcript at the HAA level and the absence of an effect on CSN2 expression suggests increased milk protein synthesis in response to EAA infusion potentially occurred through stimulation of ER biogenesis in support of greater cell secretory capacity. Palm olein infusion also increased the proportion of spliced XBP1 relative to unspliced, but this was not associated with an increase in milk protein yield. XBP1 splicing in response to palm olein and the tendency for EIF2B5 expression to be increased at the HAA level only in the presence of palm olein may indicate a protective effect of FA in mammary cells in support of persistent milk protein synthesis when EAA are in abundant supply.
Availability of data and materials
The datasets generated and analysed during the current study are available from the corresponding author upon reasonable request.
Amino acid carbon skeleton
Acetyl-CoA carboxylase 1
Activating transcription factor 6
Binding immunoglobulin protein
Cell division control protein 42 homolog
Essential amino acid
Eukaryotic translation initiation factor 2B
Fatty acid synthase
Inositol-requiring enzyme 1
Integrated stress response
Long-chain fatty acid
Myotonic dystrophy kinase-related CDC42-binding kinase
Mechanistic target of rapamycin complex 1
- NEL :
Net energy for lactation
Non-essential amino acid
Protein kinase R-like ER kinase
Sterol regulatory element binding proteins
Tricarboxylic acid cycle
Unfolded protein response
X-box binding protein 1
Scott RA, Bauman DE, Clark JH. Cellular gluconeogenesis by lactating bovine mammary tissue. J Dairy Sci. 1976;59:50–6.
Smith GH, Crabtree B, Smith R. Energy metabolism in the mammary gland. In: Mepham TB, editor. Biochemistry of lactation. Amsterdam: Elsevier; 1983. p. 121–41.
Wood HG, Joffe S, Gillespie R, Hansen RG, Hardenbrook H. Lactose synthesis. IV. The synthesis of milk constituents after unilateral injection of glycerol-1,3-C14 into the pudic artery. J Biol Chem. 1958;233:1264–70.
Wood HG, Peeters GJ, Verbeke R, Lauryssens M, Jacobson B. Estimation of the pentose cycle in the perfused cow's udder. Biochem J. 1965;96:607–15.
Owen OE, Kalhan SC, Hanson RW. The key role of anaplerosis and cataplerosis for citric acid cycle function. J Biol Chem. 2002;277:30409–12.
Bequette BJ, Sunny NE, El-Kadi SW, Owens SL. Application of stable isotopes and mass isotopomer distribution analysis to the study of intermediary metabolism of nutrients. J Anim Sci. 2006;84(Suppl):E50–9.
Bauman DE, Brown RE, Davis CL. 1970. Pathways of fatty acid synthesis and reducing equivalent generation in mammary gland of rat, sow, and cow. Arch Biochem Biophys. 1970;140:237–44.
Chilliard Y. Dietary fat and adipose tissue metabolism in ruminants, pigs, and rodents: a review. J Dairy Sci. 1993;76:3897–931.
Arriola Apelo SI, Singer LM, Lin XY, McGilliard ML, St-Pierre NR, Hanigan MD. Isoleucine, leucine, methionine, and threonine effects on mammalian target of rapamycin signaling in mammary tissue. J Dairy Sci. 2014;97:1047–56.
Doelman J, Curtis RV, Carson M, Kim JJM, Metcalf JA, Cant JP. Essential amino acid infusions stimulate mammary expression of eukaryotic initiation factor 2Bε but milk protein yield is not increased during an imbalance. J Dairy Sci. 2015;98:4499–508.
Doelman J, Kim JJM, Carson M, Metcalf JA, Cant JP. Branched-chain amino acid and lysine deficiencies exert different effects on mammary translational regulation. J Dairy Sci. 2015;98:7846–55.
Cant JP, Kim JJM, Cieslar SRL, Doelman J. Symposium review: amino acid uptake by the mammary glands: where does the control lie? J Dairy Sci. 2018;101:5655–66.
Rius AG, Appuhamy JADRN, Cyriac J, Kirovski D, Becvar O, Escobar J, et al. Regulation of protein synthesis in mammary glands of lactating dairy cows by starch and amino acids. J Dairy Sci. 2010;93:3114–27.
Toerien CA, Trout DR, Cant JP. Nutritional stimulation of milk protein yield of cows is associated with changes in phosphorylation of mammary eukaryotic initiation factor 2 and ribosomal S6 kinase 1. J Nutr. 2010;140:285–92.
Curtis RV, Kim JJM, Bajramaj DL, Doelman J, Osborne VR, Cant JP. Decline in mammary translational capacity during intravenous glucose infusion into lactating dairy cows. J Dairy Sci. 2014;97:430–8.
Nichols K, Doelman J, Kim JJM, Carson M, Metcalf JA, Cant JP. Exogenous essential amino acids stimulate an adaptive unfolded protein response in the mammary glands of lactating cows. J Dairy Sci. 2017;100:5909–21.
Nichols K, Dijkstra J, van Laar H, Kim JJM, Cant JP, Bannink A. Expression of genes related to energy metabolism and the unfolded protein response in dairy cow mammary cells is affected differently during dietary supplementation with energy from protein and fat. J Dairy Sci. 2019;102:6603–13.
Rutkowski DT, Hegde RS. Regulation of basal cellular physiology by the homeostatic unfolded protein response. J Cell Biol. 2010;189:783–94.
Nichols K, Dijkstra J, van Laar H, Pacheco S, van Valenberg HJ, Bannink A. Energy and nitrogen partitioning in dairy cows at low or high metabolizable protein levels is affected differently by postrumen glucogenic and lipogenic substrates. J Dairy Sci. 2019;102:395–412.
Nichols K, Bannink A, Doelman J, Dijkstra J. Mammary gland metabolite utilization in response to exogenous glucose or long-chain fatty acids at low and high metabolizable protein levels. J Dairy Sci. 2019;102:7150–67.
CVB (Centraal Veevoederbureau). CVB Table Ruminants 2008. Series nr. 43. CVB, The Hague, The Netherlands; 2008.
Patton S, Hutson GE. A method for isolation of milk fat globules. Lipids. 1986;21:170–4.
Lemay DG, Hovey RC, Hartono SR, Hinde K, Smilowitz JT, Ventimiglia F, et al. Sequencing the transcriptome of milk production: milk trumps mammary tissue. BMC Genomics. 2013;14:872.
Mueller O, Lightfoot S, Schroeder A. RNA integrity number (RIN)–standardization of RNA quality control. Agilent Application Note. 2004;1:1–8.
Fleige S, Pfaffl MW. RNA integrity and the effect on the real-time qRT-PCR performance. Mol Asp Med. 2006;27:126–39.
Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245–50.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25:402–8.
Glasser F, Doreau M, Ferlay A, Chilliard Y. Technical note: estimation of milk fatty acid yield from milk fat data. J Dairy Sci. 2007;90:2302–4.
Hutson GE, Patton S. Factors related to the formation of cytoplasmic crescents on milk fat globules. J Dairy Sci. 1990;73:2061–6.
Heid HW, Keenan TW. Intracellular origin and secretion of milk fat globules. Eur J Cell Biol. 2005;84:245–58.
Lemay DG, Ballard OA, Hughes MA, Morrow AL, Horseman ND, Nommsen-Rivers LA. RNA sequencing of the human milk fat layer transcriptome reveals distinct gene expression profiles at three stages of lactation. PLoS One. 2013;8:e67531.
Maningat PD, Sen P, Sunehag AL, Hadsell DL, Haymond MW. Regulation of gene expression in human mammary epithelium: effect of breast pumping. J Endocrinol. 2007;195:503–11.
Brenaut P, Bangera R, Bevilacqua C, Rebours E, Cebo C, Martin P. Validation of RNA isolated from milk fat globules to profile mammary epithelial cell expression during lactation and transcriptional response to bacterial infection. J Dairy Sci. 2012;95:6130–44.
Cánovas A, Rincón G, Bevilacqua C, Islas-Trejo A, Brenaut P, Hovey RC, et al. Comparison of five different RNA sources to examine the lactating bovine mammary gland transcriptome using RNA-sequencing. Sci Rep. 2014;4:5297.
Abdelatty AM, Iwaniuk ME, Garcia M, Moyes KM, Teter BB, Delmonte P, et al. Effect of short-term feed restriction on temporal changes in milk components and mammary lipogenic gene expression in mid-lactation Holstein dairy cows. J Dairy Sci. 2017;100:4000–13.
Bionaz M, Periasamy K, Rodriguez-Zas SL, Everts RE, Lewin HA, Hurley WL, et al. Old and new stories: revelations from functional analysis of the bovine mammary transcriptome during the lactating cycle. PLoS One. 2012;3:e33268.
Brownsey RW, Boone AN, Elliott JE, Kulpa JE, Lee WM. Regulation of acetyl-CoA carboxylase. Biochem Soc T. 2006;34:223–7.
Smith S, Witkowski A, Joshi AK. Structural and functional organization of the animal fatty acid synthase. Prog Lipid Res. 2003;42:289–317.
Bauman DE, Griinari JM. Nutritional regulation of milk fat synthesis. Annu Rev Nutr. 2003;23:203–27.
Harvatine KJ, Boisclair YR, Bauman DE. Recent advances in the regulation of milk fat synthesis. Animal. 2009;3:40–54.
Eberlé D, Hegarty B, Bossard P, Ferré P, Foufelle F. SREBP transcription factors: master regulators of lipid homeostasis. Biochimie. 2004;86:839–48.
Bionaz M, Loor JJ. Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics. 2008;9:366.
Ma L, Corl BA. Transcriptional regulation of lipid synthesis in bovine mammary epithelial cells by sterol regulatory element binding protein-1. J Dairy Sci. 2012;95:3743–55.
Harvatine KJ, Bauman DE. SREBP1 and thyroid hormone responsive spot 14 (S14) are involved in the regulation of bovine mammary lipid synthesis during diet-induced milk fat depression and treatment with CLA. J Nutr. 2006;136:2468–74.
Hansen HO, Knudsen J. Effect of exogenous long-chain fatty acids on lipid biosynthesis in dispersed ruminant mammary gland epithelial cells: esterification of long-chain exogenous fatty acids. J Dairy Sci. 1987;70:1344–9.
Bernard L, Toral PG, Chilliard Y. Comparison of mammary lipid metabolism in dairy cows and goats fed diets supplemented with starch, plant oil, or fish oil. J Dairy Sci. 2017;100:9338–51.
Fougère H, Bernard L. Effect of diets supplemented with starch and corn oil, marine algae, or hydrogenated palm oil on mammary lipogenic gene expression in cows and goats: a comparative study. J Dairy Sci. 2019;102:768–79.
Invernizzi G, Thering BJ, McGuire MA, Savoini G, Loor JJ. Sustained upregulation of stearoyl-CoA desaturase in bovine mammary tissue with contrasting changes in milk fat synthesis and lipogenic gene networks caused by lipid supplements. Funct Integr Genomics. 2010;10:561–75.
Agca C, Greenfield RB, Hartwell JR, Donkin SS. Cloning and characterization of bovine cytosolic and mitochondrial PEPCK during transition to lactation. Physiol Genomics. 2002;11:53–63.
Pongratz RL, Kibbey RG, Shulman GI, Cline GW. Cytosolic and mitochondrial malic enzyme isoforms differentially control insulin secretion. J Biol Chem. 2007;282:200–7.
Méndez-Lucas A, Hyroššová P, Novellasdemunt L, Viñals F, Perales JC. Mitochondrial phosphoenolpyruvate carboxykinase (PEPCK-M) is a pro-survival, endoplasmic reticulum (ER) stress response gene involved in tumor cell adaptation to nutrient availability. J Biol Chem. 2014;289:22090–102.
Mandella RD, Sauer LA. The mitochondrial malic enzymes. I. Submitochondrial localization and purification and properties of the NAD (P)+−dependent enzyme from adrenal cortex. J Biol Chem. 1975;250:5877–84.
Chang G, Tong L. Structure and function of malic enzymes, a new class of oxidative decarboxylases. Biochemistry. 2003;42:12721–33.
Maxin G, Ouellet DR, Lapierre H. Contribution of amino acids to glucose and lactose synthesis in lactating dairy cows. In: Oltjen JW, Kebreab E, Lapierre H, editors. Energy and protein metabolism and nutrition in sustainable animal production, vol. 134. Wageningen: Wageningen Academic Publishers; 2013. p. 443–4.
Boutinaud M, Ben Chedly MH, Delamaire E, Guinard-Flament J. Milking and feed restriction regulate transcripts of mammary epithelial cells purified from milk. J Dairy Sci. 2008;91:988–98.
Baird TD, Wek RC. Eukaryotic initiation factor 2 phosphorylation and translational control in metabolism. Adv Nutr. 2012;3:307–21.
Crozier SJ, Anthony JC, Schworer CM, Reiter AK, Anthony TG, Kimball SR, et al. Tissue-specific regulation of protein synthesis by insulin and free fatty acids. Am J Physiol Endocrinol Metab. 2003;285:E754–62.
Yonekura S, Tsuchiya M, Tokutake Y, Mizusawa M, Nakano M, Miyaji M, et al. The unfolded protein response is involved in both differentiation and apoptosis of bovine mammary epithelial cells. J Dairy Sci. 2018;101:3568–78.
Reimold AM, Iwakoshi NN, Manis J, Vallabhajosyula P, Szomolanyi-Tsuda E, Gravellese EM, et al. Plasma cell differentiation requires the transcription factor XBP-1. Nature. 2001;412:300–7.
Huh WJ, Esen E, Geahlen JH, Bredemeyer AJ, Lee AH, Shi G, et al. XBP1 controls maturation of gastric zymogenic cells by induction of MIST1 and expansion of the rough endoplasmic reticulum. Gastroenterol. 2010;139:2038–49.
Harding HP, Zeng H, Zhang Y, Jungries R, Chung P, Plesken H, et al. Diabetes mellitus and exocrine pancreatic dysfunction in Perk −/− mice reveals a role for translational control in secretory cell survival. Mol Cell. 2001;7:1153–63.
Lee A, Chu GC, Iwakoshi NN, Glimcher LH. XBP-1 is required for biogenesis of cellular secretory machinery of exocrine glands. EMBO J. 2005;24:4368–80.
Davis KR, Giesy SL, Long Q, Krumm CS, Harvatine KJ, Boisclair YR. XBP1 regulates the biosynthetic capacity of the mammary gland during lactation by controlling epithelial expansion and endoplasmic reticulum formation. Endocrinology. 2016;157:417–28.
Tsuchiya M, Koizumi Y, Hayashi S, Hanaoka M, Tokutake Y, Yonekura S. The role of unfolded protein response in differentiation of mammary epithelial cells. Biochem Bioph Res Co. 2017;484:903–8.
Wang D, Wei Y, Pagliassotti MJ. Saturated fatty acids promote endoplasmic reticulum stress and liver injury in rats with hepatic steatosis. Endocrinology. 2006;147:943–51.
Deldicque L, Cani PD, Philp A, Raymackers J, Meakin PJ, Ashford MLJ, et al. The unfolded protein response is activated in skeletal muscle by high-fat feeding: role in the downregulation of protein synthesis. Am J Physiol Endocrinol Metab. 2010;299:E695–705.
Woo CW, Kutzler L, Kimball SR, Tabas I. Toll-like receptor activation suppresses ER stress factor CHOP and translation inhibition through activation of eIF2B. Nat Cell Biol. 2012;14:192–200.
Collins KAL, Stuhlmiller TJ, Zawistowski JS, East MP, Pham TT, Hall CR, et al. Proteomic analysis defines kinase taxonomies specific for subtypes of breast cancer. Oncotarget. 2018;9:15480–97.
Bannantine JP, Bermudez LE. No holes barred: invasion of the intestinal mucosa by Mycobacterium avium subsp. paratuberculosis. Infect Immun. 2013;81:3960–5.
Kim S, Li Q, Dang CV, Lee LA. Induction of ribosomal genes and hepatocyte hypertrophy by adenovirus-mediated expression of c-Myc in vivo. PNAS. 2000;97:11198–202.
van Riggelen J, Yetil A, Felsher DW. MYC as a regulator of ribosome biogenesis and protein synthesis. Nat Rev Cancer. 2010;10:301–7.
Elouil H, Cardozo AK, Eizirik DL, Jonas JC. High glucose and hydrogenperoxide increase c-Myc and haeme-oxygenase 1 mRNA levels in rat pancreatic islets without activating NFκB. Diabetologia. 2005;48:496–505.
Jonas J, Laybutt DR, Steil GM, Trivedi N, Pertusa JG, Van de Casteele M, et al. High glucose stimulates early response gene c-Myc expression in rat pancreatic β cells. J Biol Chem. 2001;276:35375–81.
Horikawa S, Sakata K, Hatanaka M, Tsukada K. Expression of c-myc oncogene in rat liver by a dietary manipulation. Biochem Bioph Res Co. 1986;140:574–80.
Kanamoto R, Yokota T, Hayashi S. Expressions of c-myc and insulin-like growth factor-1 mRNA in the liver of growing rats vary reciprocally in response to changes in dietary protein. J Nutr. 1994;124:2329–34.
Kim HH, Park CS. A compensatory nutrition regimen during gestation stimulates mammary development and lactation potential in rats. J Nutr. 2004;134:756–61.
Bernier-Dodier P, Delbecchi L, Wagner GF, Talbot BG, Lacasse P. Effect of milking frequency on lactation persistency and mammary gland remodeling in mid-lactation cows. J Dairy Sci. 2010;93:555–64.
Nørgaard JV, Thiel PK, Sørensen MT, Sejrsen K. Cellular mechanisms in regulating mammary cell turnover during lactation and dry period in dairy cows. J Dairy Sci. 2008;91:2319–27.
The authors gratefully acknowledge technical assistance from Sven Alferink, Marcel Heetkamp, Tamme Zandstra, and the animal caretakers of CARUS (Wageningen University and Research, Wageningen, the Netherlands), and from Inge van Bakel and Koen van Helvoort (students of Wageningen University, Wageningen, the Netherlands).
This research was conducted by Wageningen University and Research (Wageningen Livestock Research, Wageningen, the Netherlands), commissioned and funded by the Ministry of Agriculture, Nature and Food Quality (The Hague, the Netherlands) within the framework of Policy Support Research theme ‘Feed4Foodure’ (BO-31.03-005-001; TKI-AF12039), and by the Vereniging Diervoederonderzoek Nederland (Rijswijk, the Netherlands).
Ethics approval and consent to participate
All experimental procedures were conducted under the Dutch Law on Animal Experiments in accordance with EU Directive 2010/63.
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Nichols, K., Bannink, A., van Baal, J. et al. Impact of post-ruminally infused macronutrients on bovine mammary gland expression of genes involved in fatty acid synthesis, energy metabolism, and protein synthesis measured in RNA isolated from milk fat. J Animal Sci Biotechnol 11, 53 (2020). https://doi.org/10.1186/s40104-020-00456-z
- Cytoplasmic crescent
- Endoplasmic reticulum biogenesis
- Mammary cell
- Milk fat globule
- Milk synthesis
- Tricarboxylic acid cycle