Genomics of metabolic adaptations in the peripartal cow

Genotype-by-environment interactions in beef and dairy cattle populations: A review of methodologies and perspectives on research and applications

Modern livestock production systems are characterized by a greater focus on intensification, involving managing larger numbers of animals to achieve higher productive efficiency and animal health and welfare within herds. Therefore, animal breeding programs need to be strategically designed to select animals that can effectively enhance production performance and animal welfare across a range of environmental conditions. Thus, this review summarizes the main methodologies used for assessing the levels of genotype-by-environment interaction (G × E) in cattle populations. In addition, we explored the importance of integrating genomic and phenotypic information to quantify and account for G × E in breeding programs. An overview of the structure of cattle breeding programs is provided to give insights into the potential outcomes and challenges faced when considering G × E to optimize genetic gains in breeding programs. The role of nutrigenomics and its impact on gene expression related to metabolism in cattle are also discussed, along with an examination of current research findings and their potential implications for future research and practical applications. Out of the 116 studies examined, 60 and 56 focused on beef and dairy cattle, respectively. A total of 83.62% of these studies reported genetic correlations across environmental gradients below 0.80, indicating the presence of G × E. For beef cattle, 69.33%, 24%, 2.67%, 2.67%, and 1.33% of the studies evaluated growth, reproduction, carcass and meat quality, survival, and feed efficiency traits, respectively. By contrast, G × E research in dairy cattle populations predominantly focused on milk yield and milk composition (79.36% of the studies), followed by reproduction and fertility (19.05%), and survival (1.59%) traits. The importance of G × E becomes particularly evident when considering complex traits such as heat tolerance, disease resistance, reproductive performance, and feed efficiency, as highlighted in this review. Genomic models provide a valuable avenue for studying these traits in greater depth, allowing for the identification of candidate genes and metabolic pathways associated with animal fitness, adaptation, and environmental efficiency. Nutrigenetics and nutrigenomics are emerging fields that require extensive investigation to maximize our understanding of gene-nutrient interactions. By studying various transcription factors, we can potentially improve animal metabolism, improving performance, health, and quality of products such as meat and milk.

Positive impacts of Nannochloropsis oculata supplementation on gene expression of immune and antioxidant markers and metabolic profile of Barki sheep in the transition period and lipogenic effects on progeny

Nannochloropsis species should be given priority when it comes to microalgae that should be added to feed since they are suitable for intense culture and have a high concentration of PUFAs (especially EPA), antioxidants, and certain vitamins. This study investigated the possible immune and antioxidant impacts of Nannochloropsis supplementation on Barki ewes during transition period and their newly born lambs. Three weeks prior to the expected time of lambing, the researched ewes were divided into two equal groups of thirty ewes each. The second group, on the other hand, was fed the same base diet as the first group plus 10 g of commercially available Nannochloropsis powder per kg of concentrate, given daily to each ewe's concentrate. Findings revealed that supplementation of ewes with Nannochloropsis significantly up-regulated the expression pattern of immune (NFKB, RANTES, HMGB1, TNF-α, IRF4, TLR7, CLA-DRB3.2, IL1B, IL6, CXCL8, S-LZ, and Cathelicidin), and antioxidant (SOD1, CAT, GPX1, GST, ATOX1, Nrf2 and AhpC/TSA) markers in ewes post-lambing and their newly born lambs. Additionally, mRNA levels of lipogenic (ACACA, FASN SCD, LPL, and BTN1A) markers were significantly up-regulated in lambs from supplemented ewes than control ones. There was a significant increase in the WBCs, Hb, RBc count, serum level of glucose, total protein, triacylglycerol and total cholesterol, GPx, catalase, IL1α and IL6 with significantly decreased serum level of TNF-α and MDA in supplemented ewes after lambing as compared with control ones. There was also a significant increase in WBCs, Hb, RBc count, birth weight and body temperature with significantly decreased in the serum levels of TNF-α and stillbirth of newly born lambs from supplemented ewes as compared to other lambs from control ones.

Effect of β-hydroxybutyrate acid on gene expression levels of antioxidant biomarkers and growth hormone-related genes in liver cell culture

Introduction

In dairy cattle, oxidative stress is a predominant problem associated with diseases and reproductive health issues. This study aimed to detect the variation in the antioxidant biomarkers by adding different concentrations of β-hydroxybutyric acid (BHBA) and sought to elucidate its effects on the gene expression levels of growth hormone (GH) and antioxidant biomarkers in bovine hepatocytes.

Material and Methods

Four antioxidant biomarkers, namely malondialdehyde (MDA), superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GSH Px) were evaluated using commercially available bovine ELISA kits. The expression levels of the bovine GH, its receptor (GHR), insulin-like growth factor (IGF), IGF-1, IGF-1 receptor, CAT, SOD, GSH-Px and β-actin (as a reference) genes in liver cell culture were determined by reverse transcriptase-PCR assay.

Results

With the increase of BHBA concentration and culture time, the activities of SOD, CAT, and GSH Px biomarkers in hepatocytes decreased. However, the content of MDA in hepatocytes increased gradually with the increase of hepatocyte culture time and BHBA concentration. The qPCR results revealed that after adding BHBA, gene expression levels of GSH-Px, SOD and IGF biomarkers in hepatocytes began to differ in the culture groups at 12 h, whereas the gene expression level of the CAT and GHR biomarkers in hepatocytes began to differ at 6 h.

Conclusion

Quantitative PCR results showed that the BHBA significantly downregulated the expression levels of the GHR gene and CAT, GSH Px and SOD antioxidant biomarker genes.

Effects of antioxidant vitamins (A, D, E) and trace elements (Cu, Mn, Se, Zn) administration on gene expression, metabolic, antioxidants and immunological profiles during transition period in dromedary camels

BACKGROUND

Nutrition has a primary role for optimum expression of genetic potential, and most of the farmers have limited resources of green fodder. Hence, a fat-soluble vitamin, especially vitamin A and E and trace elements remained most critical in the animal's ration and affects their productive and reproductive performance adversely. Animals cannot be able to produce these vitamins in their bodies; hence, an exogenous regular supply is needed to fulfil the physiological needs and to maintain high production performance. This study elucidated effects of antioxidant vitamins (A, D, E) and trace elements (Cu, Mn, Se, Zn) administration on gene expression, metabolic, antioxidants and immunological parameters in dromedary camels during transition period.

RESULTS

At 0 day, there were no appreciable differences in the expression patterns of the metabolic (IGF-I, ACACA, SCD, FASN, LPL, and BTN1A1) genes between the control and treatment groups, despite lower levels. A substantial variation in the mRNA levels of SOD1, SOD3, PRDX2, PRDX3, PRDX4, PRDX6, and AhpC/TSA was observed between the control and treatment groups, according to the antioxidant markers. In comparison to the control group, the treatment group displayed a significant up-regulation at 0 and 21 days. The treatment and control groups exhibited substantial differences in the mRNA values of IL-1α, IL-1β, IL-6, and TNFα, as indicated by immunological markers. In comparison to the control group, there was a noticeable down-regulation in the treatment group at 0 and + 21 days. But IL10 produced the opposite pattern. No significant difference was observed in glucose, cholesterol, triglyceride, HDL, total protein, NEFA, BHBA, cortisol and IGF-1 levels between control and treatment group. The activity of serum GPx, SOD and TAC was significantly affected by time and treatment x time in supplemented groups as compared with control group. IL-1, IL-1, IL-6, and TNF were noticeably greater in the control group and lower in the treatment group. Additionally, in all groups, the concentration of all pro-inflammatory cytokines peaked on the day of delivery and its lowest levels showed on day 21 following calving. The IL-10 level was at its peak 21 days prior to calving and was lowest on calving day.

CONCLUSION

The results demonstrated a beneficial effect of antioxidant vitamins and trace elements on the metabolic, antioxidant and immunological markers in dromedary camels throughout their transition period.

Identification and validation of putative biomarkers by in silico analysis, mRNA expression and oxidative stress indicators for negative energy balance in buffaloes during transition period

OBJECTIVE

Transition period is considered from 3 weeks prepartum to 3 weeks postpartum, characterized with dramatic events (endocrine, metabolic, and physiological) leading to occurrence of production diseases (negative energy balance/ketosis, milk fever etc). The objectives of our study were to analyze the periodic concentration of serum beta-hydroxy butyric acid (BHBA), glucose and oxidative markers along with identification, and validation of the putative markers of negative energy balance in buffaloes using in-silico and quantitative real time-polymerase chain reaction (qRT-PCR) assay.

METHODS

Out of 20 potential markers of ketosis identified by in-silico analysis, two were selected and analyzed by qRT-PCR technique (upregulated; acetyl serotonin o-methyl transferase like and down regulated; guanylate cyclase activator 1B). Additional two sets of genes (carnitine palmotyl transferase A; upregulated and Insulin growth factor; downregulated) that have a role of hepatic fatty acid oxidation to maintain energy demands via gluconeogenesis were also validated. Extracted cDNA (complementary deoxyribonucleic acid) from the blood of the buffaloes were used for validation of selected genes via qRTPCR. Concentrations of BHBA, glucose and oxidative stress markers were identified with their respective optimized protocols.

RESULTS

The analysis of qRT-PCR gave similar trends as shown by in-silico analysis throughout the transition period. Significant changes (p<0.05) in the levels of BHBA, glucose and oxidative stress markers throughout this period were observed. This study provides validation from in-silico and qRT-PCR assays for potential markers to be used for earliest diagnosis of negative energy balance in buffaloes.

CONCLUSION

Apart from conventional diagnostic methods, this study improves the understanding of putative biomarkers at the molecular level which helps to unfold their role in normal immune function, fat synthesis/metabolism and oxidative stress pathways. Therefore, provides an opportunity to discover more accurate and sensitive diagnostic aids.

Complementary hepatic metabolomics and proteomics reveal the adaptive mechanisms of dairy cows to the transition period

The transition period from late pregnancy to early lactation is a vital time of the lifecycle of dairy cows due to the marked metabolic challenges. Besides, the liver is the pivot point of metabolism in cattle. Nevertheless, the hepatic physiological molecular adaptation during the transition period has not been elucidated, especially from the metabolomics and proteomics view. Therefore, the present study aims to investigate the hepatic metabolic alterations in transition cows by using integrative metabolomics and proteomics methods. Gas chromatography quadrupole-time-of-flight mass spectrometry-based metabolomics and data-independent acquisition-based quantitative proteomics methods were used to analyze liver tissues collected from 8 healthy multiparous Holstein dairy cows 21 d before and after calving. In total, 44 metabolites and 250 proteins were identified as differentially expressed from 233 metabolites and 3,539 proteins detected from the liver biopsies during the transition period. Complementary functional analysis of different metabolites and proteins indicated the upregulated gluconeogenesis, tricarboxylic acid cycles, AA degradation, fatty acid oxidation, AMP-activated protein kinase signaling pathway, peroxisome proliferator-activated receptor signaling pathway, and ribosome proteins in postpartum dairy cows. In terms of the metabolites and proteins, glucose-6-phosphate, fructose-6-phosphate, carnitine palmitoyltransferase 1A, and phosphoenolpyruvate carboxykinase played a significant role in these pathways. The upregulated oxidative status may be accompanied by the pathways mentioned above. In addition, the upregulated glucagon and insulin signaling pathways also indicated the significant requirement for glucose in postpartum dairy cows. These outcomes, from the view of global metabolites and proteins, may present a better comprehension of the biology of the transition period, which can be helpful in further developing nutritional regulation strategies targeting the liver to help cows overcome this metabolically challenging time.

Effects of diacylglycerol O-acyltransferase 1 (DGAT1) on endoplasmic reticulum stress and inflammatory responses in adipose tissue of ketotic dairy cows

Adipose tissue of ketotic dairy cows exhibits greater lipolytic rate and signs of inflammation, which further aggravate the metabolic disorder. In nonruminants, the endoplasmic reticulum (ER) is a key organelle coordinating metabolic adaptations and cellular functions; thus, disturbances known as ER stress lead to inflammation and contribute to metabolic disorders. Enhanced activity of diacylglycerol O-acyltransferase 1 (DGAT1) in murine adipocytes undergoing lipolysis alleviated ER stress and inflammation. The aim of the present study was to investigate the potential role of DGAT1 on ER stress and inflammatory response of bovine adipose tissue in vivo and in vitro. Adipose tissue and blood samples were collected from cows diagnosed as clinically ketotic (n = 15) or healthy (n = 15) following a veterinary evaluation based on clinical symptoms and serum concentrations of β-hydroxybutyrate, which were 4.05 (interquartile range = 0.46) and 0.52 mM (interquartile range = 0.14), respectively. Protein abundance of DGAT1 was greater in adipose tissue of ketotic cows. Among ER stress proteins measured, ratios of phosphorylated PKR-like ER kinase (p-PERK) to PERK and phosphorylated inositol-requiring enzyme 1 (p-IRE1) to IRE1, and protein abundance of cleaved ATF6 protein were greater in adipose tissue of ketotic cows. Furthermore, ratios of phosphorylated RELA subunit of NF-κB (p-RELA) to RELA and phosphorylated c-jun N-terminal kinase (p-JNK) to JNK were greater, whereas protein abundance of NF-κB inhibitor α (NFKBIA) was lower in adipose tissue of ketotic cows. In addition, mRNA abundance of proinflammatory cytokines including TNF and IL-6 was greater in adipose tissue of ketotic cows. To better address mechanistic aspects of these responses, primary bovine adipocytes isolated from the harvested adipose tissue of healthy cows were subjected to lipolysis-stimulating conditions via incubation with 1 μM epinephrine (EPI) for 2 h. In another experiment, adipocytes were cultured with DGAT1 overexpression adenovirus and DGAT1 small interfering RNA for 48 h, respectively, followed by EPI (1 μM) exposure for 2 h. Treatment with EPI led to greater ratios of p-PERK to PERK, p-IRE1 to IRE1, p-RELA to RELA, p-JNK to JNK, and cleaved ATF6 protein, whereas EPI stimulation inhibited protein abundance of NFKBIA. Furthermore, treatment with EPI upregulated the secretion of proinflammatory cytokines into culture medium, including TNF-α and IL-6. Overexpression of DGAT1 in EPI-treated adipocytes attenuated ER stress, the activation of NF-κB and JNK signaling pathways, and the secretion of inflammatory cytokines. In contrast, silencing DGAT1 further aggravated EPI-induced ER stress and inflammatory responses. Overall, these data indicated that activation of DGAT1 may act as an adaptive mechanism to dampen metabolic dysregulation in adipose tissue. As such, it contributes to relief from ER stress and inflammatory responses.

Nutrigenomic Interventions to Address Metabolic Stress and Related Disorders in Transition Cows

For dairy cattle, the period involving a shift from late pregnancy to early lactation termed transition or periparturient is an excruciating phase. Health-related disorders are likely to happen in this time frame. Timely postpartum and metabolic adjustments to this new physical state demands correct management strategies to fulfill the cow's needs for a successful transition to this phase. Among the management strategies, one of the most researched methods for managing transition-related stress is nutritional supplementation. Dietary components directly or indirectly affect the expression of various genes that are believed to be involved in various stress-related responses during this phase. Nutrigenomics, an interdisciplinary approach that combines nutritional science with omics technologies, opens new avenues for studying the genome's complicated interactions with food. This revolutionary technique emphasizes the importance of food-gene interactions on various physiological and metabolic mechanisms. In animal sciences, nutrigenomics aims to promote the welfare of livestock animals and enhance their commercially important qualities through nutritional interventions. To this end, an increasing volume of research shows that nutritional supplementation can be effectively used to manage the metabolic stress dairy cows undergo during the transition period. These nutritional supplements, including polyunsaturated fatty acids, vitamins, dietary amino acids, and phytochemicals, have been shown to modulate energy homeostasis through different pathways, leading to addressing metabolic issues in transition cows.

Comparative transcriptome analysis of goat (Capra hircus) adipose tissue reveals physiological regulation of body reserve recovery after the peak of lactation

Adipose tissue is the energy storage organ providing energy to other tissues, including mammary gland, that supports the achievement of successive lactation cycles. Our objective was to investigate the ability of goats to restore body fat reserves by comparing lipogenic enzyme activities and by transcriptomic RNA-Seq data at two different physiological stages, mid- and post-lactation. Key lipogenic enzyme activities were higher in goat omental adipose tissue during mid-lactation (74 days in milk) than during the post-lactation period (300 days postpartum). RNA-Sequencing analysis revealed 19,271 expressed genes in the omental adipose tissue. The comparison between adipose transcriptome analysis from mid- and post-lactation goats highlighted 252 differentially expressed genes (p_adj < 0.05) between these two physiological stages. The differential expression of 11 genes was confirmed by RT-qPCR. Functional genomic analysis revealed that 31% were involved in metabolic processes among which 38% in lipid metabolism. Most of the genes involved in lipid synthesis and those in lipid transport and storage were upregulated in adipose tissue of mid- compared to post-lactation goats. In addition, adipose tissue plasticity was emphasized by genes involved in cellular signaling and tissue integrity. Network analyses also highlighted three key regulators of lipid metabolism (LEP, APOE and HNF4A) and a key target gene (VCAM1). The greatest lipogenic enzyme activities with the upregulation of genes involved in lipid metabolism highlighted a higher recovery of lipid reserves after the lactation peak than 4 months post-lactation. This study contributes to a better understanding of the molecular mechanisms controlling the body lipid reserves management during the successive lactations.

Consolidating and Upscaling Molecular Research Capacity in Nigeria: On Who's Account?

Molecular research and researchers engage in studies that seek to understand the structures, functions, and interactions of biomolecules as the basis for cellular and systemic effects in living organisms. This research approach was made possible by considerable technological advancements that equip researchers with tools to view biomolecules. Although molecular research holds great promises for improving lives and living, the technological requirements and equipment to undertake molecular research are quite expensive, often requiring a heavy start-up capital or investment. In developing countries such as Nigeria, where the majority of the population lives below the poverty line and research funding is abysmally low, such heavy investments into research that do not provide immediate solutions to societal problems are difficult. This is mostly due to limited resources available to tackle many urgent and pressing needs, and limited perspective and understanding of policymakers, leading to infrastructural and skilled personnel deficit to support molecular research. Despite all these, the field of molecular research continues to grow exponentially globally, hence, funding and investments into this critical life science research area have become imperative. With the rich biodiversity of humans, animals, and plants in Nigeria, and the huge burden of infectious diseases in the country or region, global advances in genomics and proteomics studies will be incomplete without adequate contribution from Nigeria and sub-Saharan Africa region. This paper examines the progression and challenges of undertaking molecular research in Nigeria, and how Nigerian molecular research scientists are tackling these issues, with recommendations for improved molecular research capacity and output in the country or region.

Increased adipose tissue lipolysis in dairy cows with fatty liver is associated with enhanced autophagy activity

Lipolysis is increased in adipose tissue of cows with fatty liver during the transition period. Autophagy, a major cellular degradation process, plays a critical role in adipose tissue homeostasis. The objective of this study was to explore the relationship between lipolysis and autophagy in adipose tissue of cows with fatty liver. Using a nested case-control design, we compared blood and adipose tissue samples from 10 control cows [parity: median = 3, range = 2-4; days in milk: median = 8 d, range = 5-10 d; hepatic triacylglycerol content: median = 0.55% liver wt, range = 0.48-0.61% liver wt] and 10 lactation stage-matched cows with fatty liver (parity: median = 3, range = 2-4; days in milk: median = 9 d, range = 5-11 d; hepatic triacylglycerol content: median = 6.28% liver wt, range = 2.86-7.75% liver wt). Data were analyzed using paired t-tests. Serum concentrations of free fatty acids and β-hydroxybutyrate were greater and glucose concentration was lower in cows with fatty liver, which we determined by using commercially-available kits. Furthermore, western blotting showed that increased protein abundance of ATGL (adipose triglyceride lipase), ATG5 (autophagy-related gene 5), and ATG7; ratio of phosphorylated (p)-HSL (hormone-sensitive lipase) to HSL and MAP1LC3 (microtubule-associated protein 1 light chain 3, also called LC3-II) to LC3-I along with decreased abundance of PLIN1 (perilipin 1), SQSTM1 (sequestosome-1, also called p62), and the ratio of p-mTOR (phosphorylated mechanistic target of rapamycin) to mTOR in cows with fatty liver. Quantitative reverse-transcription PCR revealed an increase in abundance of MAP1LC3 mRNA and a decrease in SQSTM1 mRNA in cows with fatty liver. These findings were replicated using an adipocyte model. Primary cultures of calf adipocytes isolated from the adipose tissue of the peritoneal omentum and mesentery were treated with 10 mM 3-methyladenine (3-MA), 5 nM rapamycin, 1 µM isoproterenol (ISO), and 1 µM ISO + 10 mM 3-MA. Comparisons among groups were analyzed using one-way ANOVA. Compared with the control, the 1 µM ISO treatment upregulated the abundance of ATGL, the ratio of p-HSL to HSL and LC3-II to LC3-I, and the glycerol content, whereas it downregulated the abundance of PLIN1 and p62 in calf adipocytes. Compared with the 1 μM ISO treatment group, 1 µM ISO + 10 mM 3-MA downregulated the abundance of ATGL, the ratio of p-HSL to HSL and LC3-II to LC3-I, and the glycerol content, whereas it upregulated the abundance of PLIN1 and p62. Compared with the control, the 5 nM rapamycin treatment upregulated the abundance of ATGL, the ratio of p-HSL to HSL and LC3-II to LC3-I, and the glycerol content, whereas it downregulated the abundance of PLIN1 and p62 in calf adipocytes. Overall, these data indicated that increased lipolysis in adipose tissue of cows with fatty liver was associated with enhanced autophagy. However, the specific molecular mechanisms that link lipolysis and autophagy need to be further investigated.

Metabogenomic analysis to functionally annotate the regulatory role of long non-coding RNAs in the liver of cows with different nutrient partitioning phenotype

Long non-coding RNAs (lncRNAs) hold gene regulatory potential, but require substantial further functional annotation in livestock. Applying two metabogenomic approaches by combining transcriptomic and metabolomic analyses, we aimed to identify lncRNAs with potential regulatory function for divergent nutrient partitioning of lactating crossbred cows and to establish metabogenomic interaction networks comprising metabolites, genes and lncRNAs. Through correlation analysis of lncRNA expression with transcriptomic and metabolomic data, we unraveled lncRNAs that have a putative regulatory role in energy and lipid metabolism, the urea and tricarboxylic acid cycles, and gluconeogenesis. Especially FGF21, which correlated with a plentitude of differentially expressed genes, differentially abundant metabolites, as well as lncRNAs, suggested itself as a key metabolic regulator. Notably, lncRNAs in close physical proximity to coding-genes as well as lncRNAs with natural antisense transcripts appear to perform a fine-tuning function in gene expression involved in metabolic pathways associated with different nutrient partitioning phenotypes.

Novel Facets of the Liver Transcriptome Are Associated with the Susceptibility and Resistance to Lipid-Related Metabolic Disorders in Periparturient Holstein Cows

Simple Summary

Energy and nutrient demands of the early lactation period can result in the development of metabolic disorders, such as ketosis and fatty liver, in dairy cows. Variability in the incidence of these disorders suggests that some cows have an ability to adapt. The objective of this study was to discover differences in liver gene expression that are associated with a cow’s susceptibility (disposition to disorder during typical conditions) or resistance (disposition to disorder onset and severity when presented a challenge) to metabolic disorders. Cows in a control treatment and a ketosis induction protocol treatment were retrospectively grouped into susceptibility and resistance groups, respectively, by a machine learning algorithm using lipid biomarker concentrations. Whole-transcriptome RNA sequencing was performed on liver samples from these cows. More susceptible cows had lower expression of glutathione metabolism genes, while less resistant cows had greater expression of eicosanoid-metabolism-related genes. Additionally, differentially expressed genes suggested a role for immune-response-related genes in conferring susceptibility and resistance to metabolic disorders. The overall inferred metabolism suggests that responses to oxidative stress may determine susceptibility and resistance to metabolic disorders, with novel implications for immunometabolism.

Abstract

Lipid-related metabolic disorders (LRMD) are prevalent in early lactation dairy cows, and have detrimental effects on productivity and health. Our objectives were to identify cows resistant or susceptible to LRMD using a ketosis induction protocol (KIP) to discover differentially expressed liver genes and metabolic pathways associated with disposition. Clustering cows based on postpartum lipid metabolite concentrations within dietary treatments identified cows more or less susceptible (MS vs. LS) to LRMD within the control treatment, and more or less resistant (MR vs. LR) within the KIP treatment. Whole-transcriptome RNA sequencing was performed on liver samples (−28, +1, and +14 days relative to calving) to assess differential gene and pathway expression (LS vs. MS; MR vs. LR; n = 3 cows per cluster). Cows within the MS and LR clusters had evidence of greater blood serum β-hydroxybutyrate concentration and liver triglyceride content than the LS and MR clusters, respectively. The inferred metabolism of differentially expressed genes suggested a role of immune response (i.e., interferon-inducible proteins and major histocompatibility complex molecules). Additionally, unique roles for glutathione metabolism and eicosanoid metabolism in modulating susceptibility and resistance, respectively, were implicated. Overall, this research provides novel insight into the role of immunometabolism in LRMD pathology, and suggests the potential for unique control points for LRMD progression and severity.

Methionine supply during the peripartum period and early lactation alter immunometabolic gene expression in cytological smear and endometrial tissue of holstein cows

The objective of the present study was to evaluate the effect of feeding rumen-protected methionine (RPM) during the peripartal period and early lactation on mRNA gene expression profiles of uterine cytological smear and endometrial samples of Holstein cows (n = 20). Treatments consisted of a supplementation with RPM [MET; n = 11; RPM at a rate of 0.08 % of DM: Lys:Met = 2.8:1, (Smartamine® M Adisseo, Alpharetta, GA, USA)] and no supplementation (CON; n = 9; Lys:Met = 3.5:1). Uterine cytology smears and endometrial samples were collected at 15, 30, and 73 days in milk (DIM) and analyzed for expression of genes related with metabolism, inflammation, and methionine metabolism. Regarding the cytological smear samples, RPM supplementation tended to increase mRNA expression of methionine adenosyltransferase 1 alpha (MAT1A) and increased the mRNA expression of fibroblast growth factor 7 (FGF7), with an effect of time for the latter. On the other hand, RPM decreased mRNA expression for glucose transporter 4 (GLUT4), interleukin 1 beta (IL-1β), interleukin 6 (IL-6), interleukin 8 (IL-8), prostaglandin E synthase 3 (PTGES3), translocator protein 18 kDa (TSPO), mucin 1 (MUC1) and superoxide dismutase (SOD1) in cytological smear samples. There was an effect of time for all variables except MAT1A, with decreasing expression over time. There was a TRT × TIME interaction for GLUT4 mRNA expression, with higher GLUT4 mRNA expression for cows fed CON than for cows fed RPM at time 15 and a tendency to higher expression for cows fed CON on time 30 when compared with cows fed RPM. For uterine tissue samples, feeding RPM increased the mRNA expression of lecithin-cholesterol acyltransferase (LCAT), S-adenosyl-l-homocysteine hydrolase (SAAH), FGF7, GLUT4, and apolipoproteins 3 (APOL3), with an effect of time for APOL3 where its expression increased over time. There was a tendency for cows fed RPM to have decreased IL1β mRNA expression. In conclusion, feeding RPM during transition period and early lactation is beneficial for uterine immune response and metabolism in early lactation as indicated by the favorable expressions of genes affecting the uterine immunometabolism during such a challenging period.

Disruption of endoplasmic reticulum homeostasis exacerbates liver injury in clinically ketotic cows

Disruption of endoplasmic reticulum (ER) homeostasis, a condition termed "ER stress," contributes to the development of liver injury in nonruminants. Because liver injury is a prominent pathological feature associated with overproduction of ketone bodies in dairy cows with ketosis, understanding the ER stress state and its functional consequences on liver injury is of particular interest. Here, 30 multiparous cows (within 3 wk postpartum) classified based on blood β-hydroxybutyrate (BHB) as healthy (n = 15, BHB <0.6 mM) or clinically ketotic (n = 15, BHB >3.0 mM) were used. Compared with healthy cows, ketotic cows had greater levels of serum fatty acids and activities of serum aspartate aminotransferase, alanine aminotransferase, alkaline phosphatase, γ-glutamyl transferase, and glutamate dehydrogenase but lower serum glucose. Furthermore, dairy cows with ketosis had greater protein abundance of ER stress markers in liver tissue, including protein kinase RNA-like ER kinase (PERK), inositol-requiring protein-1α (IRE1α), and cleaved activating transcription factor-6 (ATF6). Cows with ketosis also had higher mRNA levels of hepatic 78-kDa glucose-regulated protein (GRP78) and spliced X-box binding protein 1 (sXBP1). These data confirmed an enhanced ER stress state in clinically ketotic cows. To explore whether enhanced hepatic ER stress was induced by elevated ketone bodies and the possible contribution of ER stress to liver injury, in vitro experiments were then performed using isolated primary calf hepatocytes treated with incremental concentrations of BHB (0, 0.6, 1.2, 3.0, and 4.8 mM) for 12 h with or without overexpression of GRP78 (the master regulator of unfolded protein response). Phosphorylation levels of PERK and IRE1α proteins, level of cleaved ATF6 protein, and mRNA abundance of GRP78 and sXBP1 in hepatocytes increased after treatment with high (3.0 and 4.8 mM) BHB, indicating a mechanistic link between excessive BHB and enhanced hepatic ER stress. Furthermore, treatment with 3.0 and 4.8 mM BHB markedly elevated activities of aspartate aminotransferase and alanine aminotransferase in cell supernatant, indicating exacerbated hepatocyte damage after ER stress was enhanced. Overexpression of GRP78 attenuated both BHB-induced ER stress and the ensuing cellular damage, suggesting that hepatocyte damage caused by excessive BHB can be mediated via enhanced ER stress. Overall, the present study revealed that ER stress may exacerbate liver injury development in clinically ketotic cows, underscoring the biological relevance of this pathway in the context of liver injury.

Feeding a Saccharomyces cerevisiae fermentation product improves udder health and immune response to a Streptococcus uberis mastitis challenge in mid-lactation dairy cows

Background

We aimed to characterize the protective effects and the molecular mechanisms of action of a Saccharomyces cerevisiae fermentation product (NTK) in response to a mastitis challenge. Eighteen mid-lactation multiparous Holstein cows (n = 9/group) were fed the control diet (CON) or CON supplemented with 19 g/d NTK for 45 d (phase 1, P1) and then infected in the right rear quarter with 2500 CFU of Streptococcus uberis (phase 2, P2). After 36-h, mammary gland and liver biopsies were collected and antibiotic treatment started until the end of P2 (9 d post challenge). Cows were then followed until day 75 (phase 3, P3). Milk yield (MY) and dry matter intake (DMI) were recorded daily. Milk samples for somatic cell score were collected, and rectal and udder temperature, heart and respiration rate were recorded during the challenge period (P2) together with blood samples for metabolite and immune function analyses. Data were analyzed by phase using the PROC MIXED procedure in SAS. Biopsies were used for transcriptomic analysis via RNA-sequencing, followed by pathway analysis.

Results

DMI and MY were not affected by diet in P1, but an interaction with time was recorded in P2 indicating a better recovery from the challenge in NTK compared with CON. NTK reduced rectal temperature, somatic cell score, and temperature of the infected quarter during the challenge. Transcriptome data supported these findings, as NTK supplementation upregulated mammary genes related to immune cell antibacterial function (e.g., CATHL4, NOS2), epithelial tissue protection (e.g. IL17C), and anti-inflammatory activity (e.g., ATF3, BAG3, IER3, G-CSF, GRO1, ZFAND2A). Pathway analysis indicated upregulation of tumor necrosis factor α, heat shock protein response, and p21 related pathways in the response to mastitis in NTK cows. Other pathways for detoxification and cytoprotection functions along with the tight junction pathway were also upregulated in NTK-fed cows.

Conclusions

Overall, results highlighted molecular networks involved in the protective effect of NTK prophylactic supplementation on udder health during a subclinical mastitic event.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00560-8.

Glucose and Fatty Acid Metabolism of Dairy Cows in a Total Mixed Ration or Pasture-Based System During Lactation

In this study, we explored mechanisms related to glucose and fatty acid metabolism in Holstein–Friesian multiparous dairy cows during lactation under two feeding strategies. From 0 to 180 days postpartum, cows were fed total mixed ration (TMR) ad libitum (non-grazing group, G0) or grazed Festuca arundinacea or Medicago sativa and were supplemented with 5.4 kg DM/d of an energy-protein concentrate (grazing group, G1). From 180 to 250 days postpartum, all cows grazed F. arundinacea and were supplemented with TMR. Plasma samples and liver biopsies were collected at −14, 35, 60, 110, 180, and 250 days in milk (DIM) for metabolite, hormone, gene expression, and western blot analysis. Our results showed increased levels of negative energy balance markers: plasma non-esterified fatty acids (NEFA), liver triglyceride and plasma β-hydroxybutyrate (BHB) (P &lt; 0.01), triglyceride and β-hydroxybutyrate concentration were especially elevated for G1 cows. Also, hepatic mRNA expression of gluconeogenic enzymes was upregulated during early lactation (P &lt; 0.05). In particular, methymalonyl-CoA mutase expression was increased for G0 cows (P &lt; 0.05) while pyruvate carboxylase (PC) expression was increased for G1 cows (P &lt; 0.05), suggesting differential gluconeogenic precursors for different feeding strategies. Phosphorylation of AMP-activated protein kinase was increased in early lactation vs. late lactation (P &lt; 0.01) and negatively correlated with PC mRNA levels. The positive association of gluconeogenic genes with proliferator-activated receptor gamma coactivator 1-alpha (PPARGC1A) hepatic expression supported the importance of this transcription factor in glucose metabolism. The peroxisome proliferator-activated receptor alpha (PPARA) mRNA was increased during early lactation (P &lt; 0.05), and was positively associated to PPARGC1A, carnitine palmitoyl-transferase 1, and hydroxymethylglutaryl-CoA synthase 2 (HMGCS2) mRNA expression. Alongside, hepatic mRNA expression of FABP was decreased for G1 vs. G0 cows (P &lt; 0.05), possibly linked to impaired fatty acid transport and related to accumulation of liver triglycerides, evidencing G1 cows fail to adapt to the demands of early lactation. In sum, our results showed that metabolic adaptations related to early lactation negative energy balance can be affected by feeding strategy and might be regulated by the metabolic sensors AMPK, SIRT1, and coordinated by transcription factors PPARGC1A and PPARA.

Glucose metabolism and the somatotropic axis in dairy cows after abomasal infusion of essential fatty acids together with conjugated linoleic acid during late gestation and early lactation

Sufficient glucose availability is crucial for exploiting the genetic potential of milk production during early lactation, and endocrine changes are mainly related to repartitioning of nutrient supplies toward the mammary gland. Long-chain fatty acids, such as essential fatty acids (EFA) and conjugated linoleic acid (CLA), have the potential to improve negative energy balance and modify endocrine changes. In the present study, the hypothesis that combined CLA and EFA treatment supports glucose metabolism around the time of calving and stimulates insulin action and the somatotropic axis in cows in an additive manner was tested. Rumen-cannulated German Holstein cows (n = 40) were investigated from wk 9 antepartum (AP) until wk 9 postpartum (PP). The cows were abomasally supplemented with coconut oil (CTRL, 76 g/d); 78 g/d of linseed and 4 g/d of safflower oil (EFA); Lutalin (CLA, isomers cis-9,trans-11 and trans-10,cis-12 CLA, each 10 g/d); or the combination of EFA+CLA. Blood samples were collected several times AP and PP to determine the concentrations of plasma metabolites and hormones related to glucose metabolism and the somatotropic axis. Liver tissue samples were collected several days AP and PP to measure glycogen concentration and the mRNA abundance of genes related to gluconeogenesis and the somatotropic axis. On d 28 AP and 21 PP, endogenous glucose production (eGP) and glucose oxidation (GOx) were measured via tracer technique. The concentration of plasma glucose was higher in CLA than in non-CLA-treated cows, and the plasma β-hydroxybutyrate concentration was higher in EFA than in non-EFA cows on d 21 PP. The eGP increased from AP to PP with elevated eGP in EFA and decreased eGP in CLA-treated cows; GOx was lower in CLA than in CTRL on d 21 PP. The plasma insulin concentration decreased after calving in all groups and was higher in CLA than in non-CLA cows at several time points. Plasma glucagon and cortisol concentrations on d 21 PP were lower in CLA than non-CLA groups. The glucagon/insulin and glucose/insulin ratios were higher in CTRL than in CLA group during the transition period. Plasma IGF-I concentration was lower in EFA than non-EFA cows on d 42 AP and was higher during the dry period and early lactation in CLA than in non-CLA cows. The IGF binding protein (IGFBP)-3/-2 ratio in blood plasma was higher in CLA than in non-CLA cows. Hepatic glycogen concentration on d 28 PP was higher, but the mRNA abundance of PC and IGFBP2 was lower in CLA than non-CLA cows on d 1 PP. The EFA treatment decreased the mRNA abundance of IGFBP3 AP and PCK1, PCK2, G6PC, PCCA, HMGCS2, IGFBP2, and INSR at several time points PP. Results indicated elevated concentrations of plasma glucose and insulin along with the stimulation of the somatotropic axis in cows treated with CLA, whereas EFA treatment stimulated eGP but not mRNA abundance related to eGP PP. The systemic effects of the combined EFA+CLA treatment were very similar to those of CLA treatment, but the effects on hepatic gene expression partially corresponded to those of EFA treatment.

Effects of a combined essential fatty acid and conjugated linoleic acid abomasal infusion on metabolic and endocrine traits, including the somatotropic axis, in dairy cows

The objective of this study was to test the effects of essential fatty acids (EFA), particularly α-linolenic acid (ALA), and conjugated linoleic acid (CLA) supplementation on metabolic and endocrine traits related to energy metabolism, including the somatotropic axis, in mid-lactation dairy cows. Four cows (126 ± 4 d in milk) were used in a dose-escalation study design and were abomasally infused with coconut oil (CTRL; 38.3 g/d; providing saturated fatty acids), linseed and safflower oils (EFA; 39.1 and 1.6 g/d; n-6:n-3 FA ratio = 1:3), Lutalin (CLA; cis-9,trans-11 and trans-10,cis-12 CLA, 4.6 g/d of each), or EFA and CLA (EFA+CLA) for 6 wk. The initial dosage was doubled twice after 2 wk, resulting in 3 dosages (dosages 1, 2, and 3). Each cow received each fat treatment at different times. Cows were fed with a corn silage-based total mixed ration providing a low-fat content and a high n-6:n-3 fatty acid ratio. Plasma concentrations of metabolites and hormones (insulin-like growth factor-binding proteins only on wk 0 and 6) were analyzed at wk 0, 2, 4, and 6 of each treatment period. Liver biopsies were taken before starting the trial and at wk 6 of each treatment period to measure hepatic mRNA abundance of genes linked to glucose, cholesterol and lipid metabolism, and the somatotropic axis. The changes in the milk and blood fatty acid patterns and lactation performance of these cows have already been published in a companion paper. The plasma concentration of total cholesterol increased with dosage in all groups, except CLA, reaching the highest levels in EFA+CLA and CTRL compared with CLA. The high-density lipoprotein cholesterol plasma concentration increased in CTRL and was higher than that in EFA and CLA, whereas the concentration of low-density lipoprotein cholesterol increased in a dose-dependent manner in EFA and EFA+CLA, and was higher than that in CLA. Hepatic mRNA expression of 3-hydroxy-3-methyl-glutaryl-CoA synthase 1 was upregulated in all groups but was highest in EFA+CLA. Expression of sterol regulatory element-binding factor 1 tended to be lowest due to EFA treatment, whereas expression of long chain acyl-CoA-synthetase was lower in EFA than in CTRL. Hepatic mRNA expression of GHR1A tended to be higher in EFA+CLA than in CTRL. The plasma concentration of insulin-like growth factor I increased in CLA, and the plasma IGFBP-2 concentration was lower in EFA+CLA than in CTRL at wk 6. The plasma concentration of adiponectin decreased in EFA+CLA up to dosage 2. Plasma concentrations of albumin and urea were lower in CLA than in CTRL throughout the experimental period. Supplementation with EFA and CLA affected cholesterol and lipid metabolism and their regulation differently, indicating distinct stimulation after the combined EFA and CLA treatment. The decreased IGFBP-2 plasma concentration and upregulated hepatic mRNA abundance of GHR1A in EFA+CLA-supplemented cows indicated the beneficial effect of the combined EFA and CLA treatment on the somatotropic axis in mid-lactation dairy cows. Moreover, supplementation with CLA might affect protein metabolism in dairy cows.

Role of peroxisome proliferator-activated receptor-α on the synthesis of monounsaturated fatty acids in goat mammary epithelial cells

A key member of the nuclear receptor superfamily is the peroxisome proliferator-activated receptor alpha (PPARA) isoform, which in nonruminants is closely associated with fatty acid oxidation. Whether PPARA plays a role in milk fatty acid synthesis in ruminants is unknown. The main objective of the present study was to use primary goat mammary epithelial cells (GMEC) to activate PPARA via the agonist WY-14643 (WY) or to silence it via transfection of small-interfering RNA (siRNA). Three copies of the peroxisome proliferator-activated receptor response element (PPRE) contained in a luciferase reporter vector were transfected into GMEC followed by incubation with WY at 0, 10, 20, 30, 50, or 100 µM. A dose of 50 µM WY was most effective at activating PPRE without influencing PPARA mRNA abundance. Transfecting siRNA targeting PPARA decreased its mRNA abundance to 20% and protein level to 50% of basal levels. Use of WY upregulated FASN, SCD1, ACSL1, DGAT1, FABP4, and CD36 (1.1-, 1.5-, 2-, 1.4-, 1.5-, and 5-fold, respectively), but downregulated DGAT2 and PGC1A (-20% and -40%, respectively) abundance. In contrast, triacylglycerol concentration decreased and the content and desaturation index of C16:1 and C18:1 increased. Thus, activation of PPARA via WY appeared to channel fatty acids away from esterification. Knockdown of PPARA via siRNA downregulated ACACA, SCD1, AGPAT6, CD36, HSL, and SREBF1 (-43%, -67%, -16%, -56%, -26%, and -29%, respectively), but upregulated ACSL1, DGAT2, FABP3, and PGC1A (2-, 1.4-, 1.3-, and 2.5-fold, respectively) mRNA abundance. A decrease in the content and desaturation index of C16:1 and C18:1 coupled with an increase in triacylglycerol content accompanied those effects at the mRNA level. Overall, data suggest that PPARA could promote the synthesis of MUFA in GMEC through its effects on mRNA abundance of genes related to fatty acid synthesis, oxidation, transport, and triacylglycerol synthesis.

Body condition alters glutathione and nuclear factor erythroid 2-like 2 (NFE2L2)-related antioxidant network abundance in subcutaneous adipose tissue of periparturient Holstein cows

Dairy cows with high body condition score (BCS) in late prepartum are more susceptible to oxidative stress (OS). Nuclear factor erythroid 2-like 2 (NFE2L2) is a major antioxidant transcription factor. We investigated the effect of precalving BCS on blood biomarkers associated with OS, inflammation, and liver function, along with mRNA and protein abundance of targets related to NFE2L2 and glutathione (GSH) metabolism in s.c. adipose tissue (SAT) of periparturient dairy cows. Twenty-two multiparous Holstein cows were retrospectively classified into a high BCS (HBCS; n = 11, BCS ≥3.5) or normal BCS (NBCS; n = 11, BCS ≤3.17) on d 28 before parturition. Cows were fed a corn silage- and wheat straw-based total mixed ration during late prepartum, and a corn silage- and alfalfa hay-based total mixed ration postpartum. Blood samples obtained at -10, 7, 15, and 30 d relative to parturition were used for analyses of biomarkers associated with inflammation, including albumin, ceruloplasmin, haptoglobin, and myeloperoxidase, as well as OS, including ferric reducing ability of plasma (FRAP), reactive oxygen species (ROS), and β-carotene. Adipose biopsies harvested at -15, 7, and 30 d relative to parturition were analyzed for mRNA (real-time quantitative PCR) and protein abundance (Western blotting) of targets associated with the antioxidant transcription regulator nuclear factor, NFE2L2, and GSH metabolism pathway. In addition, concentrations of GSH, ROS and malondialdehyde were measured. High BCS cows had lower prepartum dry matter intake expressed as a percentage of body weight along with greater BCS loss between -4 and 4 wk relative to parturition. Plasma concentrations of ROS and FRAP increased after parturition regardless of treatment. Compared with NBCS, HBCS cows had greater concentrations of FRAP at d 7 postpartum, which coincided with peak values in those cows. In addition, NBCS cows experienced a marked decrease in plasma ROS after d 7 postpartum, while HBCS cows maintained a constant concentration by d 30 postpartum. Overall, ROS concentrations in SAT were greater in HBCS cows. However, overall mRNA abundance of NFE2L2 was lower and cullin 3 (CUL3), a negative regulator of NFE2L2, was greater in HBCS cows. Although HBCS cows had greater overall total protein abundance of NFE2L2 in SAT, ratio of phosphorylated NFE2L2 to total NFE2L2 was lower, suggesting a decrease in the activity of this antioxidant system. Overall, mRNA abundance of the GSH metabolism-related genes glutathione reductase (GSR), glutathione peroxidase 1 (GPX1), and transaldolase 1 (TALDO1), along with protein abundance of glutathione S-transferase mu 1 (GSTM1), were greater in HBCS cows. Data suggest that HBCS cows might experience greater systemic OS after parturition, while increased abundance of mRNA and protein components of the GSH metabolism pathway in SAT might help alleviate tissue oxidant status. Data underscored the importance of antioxidant mechanisms at the tissue level. Thus, targeting these pathways in SAT during the periparturient period via nutrition might help control tissue remodeling while allowing optimal performance.

Short communication: Enhanced autophagy activity in liver tissue of dairy cows with mild fatty liver

During the transition period, dairy cows are challenged by increased energy demands and decreased dry matter intake, which can induce a variety of metabolic disorders, especially fatty liver. Dairy cows suffering from mild or moderate fatty liver in this period show no distinct clinical symptoms, indicating the occurrence of adaptive processes. The process of autophagy (an adaptive response) leads to degradation of intracellular components to generate energy and maintains cellular homeostasis during negative nutrient status. Whether autophagy is involved in metabolic adaptations of the pathological course of mild fatty liver is unclear. Thus, the aim of this study was to determine hepatic autophagy status in dairy cows with mild fatty liver. Liver samples were collected from healthy cows (n = 15), defined as having hepatic triglyceride (TG) content <1% on a wet weight basis, and cows with mild fatty liver (n = 15), defined as having hepatic TG content between 1 and 5%. The abundance of the ubiquitinated proteins, microtubule-associated protein 1 light chain 3 (MAP1LC3, also called LC3-II) and sequestosome-1 (SQSTM1, also called p62) was lower, whereas the mRNA abundance of MAP1LC3 and SQSTM1 was greater in cows with mild fatty liver. The hepatic mRNA abundance of autophagy-related (ATG) genes ATG5 and ATG7 was greater in response to fatty liver. However, the protein abundance of ATG5 and ATG7 did not differ between healthy and mild fatty liver cows. Together, these data indicate that the formation and degradation of autophagosomes is enhanced in the liver of cows with mild fatty liver. Besides, these results are conducive to define the adaptation mechanisms of dairy cows during the transition period.

Hepatic metabolomics and transcriptomics to study susceptibility to ketosis in response to prepartal nutritional management

Background

Ketosis in dairy cows is associated with body fat mobilization during the peripartal period. Sub-clinical and clinical ketosis arise more frequently in cows that are overfed energy during the entire dry (last 50 to 45 days prior to parturition) or close-up period (last ~ 28 days prepartum).

Methods

A retrospective analysis was performed on 12 cows from a larger cohort that were fed a higher-energy diet [1.54 Mcal/kg of dry matter (DM); 35.9% of DM corn silage and 13% of DM ground corn] during the close-up dry period, of which 6 did not develop clinical ketosis (OVE, 0.83 mmol/L plasma hydroxybutyrate; BHB) and 6 were diagnosed with clinical ketosis (KET, 1.4 mmol/L BHB) during the first week postpartum. A whole-transcriptome bovine microarray (Agilent Technologies) and metabolomics (GC-MS, LC-MS; Metabolon® Inc.) were used to perform transcript and metabolite profiling of liver tissue harvested at − 10 days relative to parturition which allowed to establish potential associations between prepartal transcriptome/metabolome profiles and susceptibility to clinical ketosis postpartum.

Results

Cows in KET had greater (P = 0.01) overall body weight between − 2 and 1 week around parturition, but similar body condition score than OVE. Although dry matter intake (DMI) did not differ prepartum, KET cows had lower (P < 0.01) DMI and similar milk yield as OVE cows during the first week postpartum. Transcriptome analysis revealed a total of 3065 differentially expressed genes (DEG; P ≤ 0.05) in KET. Metabolomics identified 15 out of 313 total biochemical compounds significantly affected (P ≤ 0.10) in KET. Among those, greater concentrations (P ≤ 0.06, + 2.3-fold) of glycochenodeoxycholate in KET cows also have been detected in humans developing non-alcoholic fatty liver disease. Bioinformatics analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database and the DEG revealed that, among the top 20 most-impacted metabolic pathway categories in KET, 65% were overall downregulated. Those included ‘Metabolism of cofactors and vitamins’, ‘Biosynthesis of other secondary metabolites’, ‘Lipid’, ‘Carbohydrate’, and ‘Glycan biosynthesis and metabolism’. The lower relative concentration of glucose-6-phosphate and marked downregulation of fructose-1,6-bisphosphatase 2 and pyruvate dehydrogenase kinase 4 support a strong impairment in gluconeogenesis in prepartal liver of cows developing KET postpartum. Among the top 20 most-impacted non-metabolic pathways, 85% were downregulated. Pathways such as ‘mTOR signalling’ and ‘Insulin signalling’ were among those. ‘Ribosome’, ‘Nucleotide excision repair’, and ‘Adherens junctions’ were the only upregulated pathways in cows with KET.

Conclusions

The combined data analyses revealed more extensive alterations of the prepartal liver transcriptome than metabolome in cows overfed energy and developing ketosis postpartum. The causative link between these tissue-level adaptations and onset of clinical ketosis needs to be studied further.

Genome-wide association studies for heat stress response in Bos taurus × Bos indicus crossbred cattle

Heat stress is an important issue in the global dairy industry. In tropical areas, an alternative to overcome heat stress is the use of crossbred animals or synthetic breeds, such as the Girolando. In this study, we performed a genome-wide association study (GWAS) and post-GWAS analyses for heat stress in an experimental Gir × Holstein F₂ population. Rectal temperature (RT) was measured in heat-stressed F₂ animals, and the variation between 2 consecutive RT measurements (ΔRT) was used as the dependent variable. Illumina BovineSNP50v1 BeadChip (Illumina Inc., San Diego, CA) and single-SNP approach were used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene-transcription factor (TF) networks, generated from enriched TF. The breed origin of marker alleles in the F₂ population was assigned using the breed of origin of alleles (BOA) approach. Heritability and repeatability estimates (± standard error) for ΔRT were 0.13 ± 0.08 and 0.29 ± 0.06, respectively. Association analysis revealed 6 SNP significantly associated with ΔRT. Genes involved with biological processes in response to heat stress effects (LIF, OSM, TXNRD2, and DGCR8) were identified as putative candidate genes. After performing the BOA approach, the 10% of F₂ animals with the lowest breeding values for ΔRT were classified as low-ΔRT, and the 10% with the highest breeding values for ΔRT were classified as high-ΔRT. On average, 49.4% of low-ΔRT animals had 2 alleles from the Holstein breed (HH), and 39% had both alleles from the Gir breed (GG). In high-ΔRT animals, the average proportion of animals for HH and GG were 1.4 and 50.2%, respectively. This study allowed the identification of candidate genes for ΔRT in Gir × Holstein crossbred animals. According to the BOA approach, Holstein breed alleles could be associated with better response to heat stress effects, which could be explained by the fact that Holstein animals are more affected by heat stress than Gir animals and thus require a genetic architecture to defend the body from the deleterious effects of heat stress. Future studies can provide further knowledge to uncover the genetic architecture underlying heat stress in crossbred cattle.

Liver fatty acid metabolism associations with reproductive performance of dairy cattle

The peri-calving period is characterized by a negative energy balance, which leads to lipid mobilization. Thus, during this period, the liver has important functions related to optimizing milk yield, preventing metabolic and infectious diseases, and improving fertility. To clarify the relationship between liver fatty acid metabolism and reproductive performance, the present study was conducted to assess the abundance of specific hepatic proteins related to lipid metabolism in both plasma and follicular fluid in dairy cattle with different days to conception (DC). Sixteen animals were grouped according to DC, as more and fewer DC (MDC and FDC, respectively). Blood and liver biopsies were sampled 14 days before the expected calving date and 4, 14 and 28 days after calving. The plasma beta-hydroxybutyric acid (BHBA) concentrations and the liver triacylglycerol (TAG) content were greater in the MDC group (P <  0.05), whereas the protein abundance of carnitine palmitoyl transferase 1 was greater in the FDC group (P < 0.05). Additionally, total bilirubin (TBil) concentration was less in the FDC than MDC group on day 28 (P < 0.05). These results indicate lipid mobilization and liver fatty acid oxidation capacity in dairy cows could contribute to the adaptations and reproductive performance.

Fatty acid-induced endoplasmic reticulum stress promoted lipid accumulation in calf hepatocytes, and endoplasmic reticulum stress existed in the liver of severe fatty liver cows

Disruption of endoplasmic reticulum (ER) homeostasis, often termed ER stress, is intrinsically linked with perturbation of lipid metabolism in humans and mice. Whether ER homeostasis is affected in cows experiencing fatty liver is unknown. The aim of this study was to investigate the potential role of ER stress in hepatic lipid accumulation in calf hepatocytes and ER stress status in dairy cows with severe fatty liver. In vitro experiments were conducted in which hepatocytes were isolated from calves and treated with different concentrations of fatty acids, tauroursodeoxycholic acid (TUDCA; a canonical inhibitor of ER stress), or both. The increase in phosphorylation level of protein kinase RNA-like ER kinase (PERK) and inositol requiring protein-1α (IRE1α) proteins, and the cleavage of activating transcription factor-6 (ATF6) protein in response to increasing doses of fatty acids (which were reversed by TUDCA treatment) in primary hepatocytes underscored a mechanistic link between fatty acids and ER stress. In addition, fatty acid treatment increased the abundance of sterol regulatory element-binding protein 1c, acetyl-CoA carboxylase-α, fatty acid synthase, and diacylglycerol acyltransferase 1, and lipid accumulation in calf primary hepatocytes, whereas inhibition of ER stress by incubating with TUDCA significantly weakened these effects. Overall, results in vitro indicate that inhibition of ER stress in calf hepatocytes alleviates fatty acid-induced lipid accumulation by downregulating the expression of lipogenic genes. In vivo experiments, liver and blood samples were collected from cows diagnosed as healthy (n = 15) or with severe fatty liver (n = 15). The phosphorylation level of PERK and IRE1α, the cleavage of ATF6 protein, and the abundance of several unfolded protein response genes (78 kDa glucose-regulated protein, AMP-dependent transcription factor 4, and spliced X-box binding protein 1) were greater in liver of cows with severe fatty liver. The present in vivo study confirms the occurrence of ER stress in dairy cows with severe fatty liver. Considering the causative role of fatty acid-induced ER stress in hepatic lipid accumulation in calf hepatocytes, the existence of ER stress in the liver of severe fatty liver cows may presage its participation in fatty liver progression in dairy cows. However, the mechanistic relationship between ER stress and fatty liver in dairy cows remain to be determined.

Glutathione metabolism and nuclear factor erythroid 2-like 2 (NFE2L2)-related proteins in adipose tissue are altered by supply of ethyl-cellulose rumen-protected methionine in peripartal Holstein cows

Enhancing the supply of rumen-protected Met (RPM) during the peripartum period alleviates inflammation and oxidative stress status in dairy cows. We tested the hypothesis that RPM could increase abundance of genes and proteins related to glutathione (GSH) metabolism and the antioxidant transcription factor nuclear factor erythroid 2-like 2 (NFE2L2) in subcutaneous adipose tissue. Multiparous Holstein cows were fed a basal diet [control prepartum diet = 1.47 Mcal/kg of dry matter (DM) and 15.3% crude protein; control postpartum diet = 1.67 Mcal/kg of DM and 17.7% crude protein] or the control plus ethyl-cellulose RPM at a rate of 0.09 and 0.10% of DM intake before expected calving and after calving, respectively. Sixty cows were assigned to treatments based on parity, previous 305-d milk yield, and body condition score at 28 d from parturition. Diets were fed from -28 to 30 d. Biopsies of subcutaneous adipose tissue collected on d -10, 10, and 30 relative to parturition from 7 cows in each group were used for measuring concentrations of GSH, reactive oxygen species, superoxide dismutase, malondialdehyde, and mRNA and protein abundance (Western blotting). A repeated-measures ANOVA was used for statistics. The statistical model included the random effect of block and fixed effects of treatment, time, and its interaction. There was a diet × time effect for reactive oxygen species due to lower concentrations in Met versus control cows specifically at d -10. Cows fed Met also had lower concentrations of malondialdehyde in subcutaneous adipose tissue. Compared with controls, overall mRNA abundance of the GSH metabolism-related genes cystathionine-β-synthase (CBS), glutamate-cysteine ligase modifier subunit (GCLM), glutathione reductase (GSR), and glutathione peroxidase 1 (GPX1) was greater in cows fed Met. Furthermore, supply of Met resulted in an overall upregulation of protein abundance of glutathione peroxidase (GPX) 1, GPX3, glutathione S-transferase mu 1 (GSTM1), and glutathione S-transferase α 4 (GSTA4), all related to GSH metabolism. There was a diet × time effect for protein abundance of NFE2L2 and its repressor Kelch-like ECH associated protein 1 (KEAP1) due to lower values at 30 d in cows fed Met versus controls. The abundance of phosphorylated NFE2L2 was lower at 30 d in response to Met. Overall, the data suggest that exogenous Met may play a role in activating GSH metabolism and the antioxidant NFE2L2 pathways in subcutaneous adipose tissue.

Invited review: Metabolic challenges and adaptation during different functional stages of the mammary gland in dairy cows: Perspectives for sustainable milk production

Milk production of dairy cows has increased markedly during recent decades and continues to increase further. The evolutionarily conserved direction of nutrients to the mammary gland immediately after calving provided the basis for successful selective breeding toward higher performance. Considerable variation in adaptive responses toward energy and nutrient shortages exists; however, this variation in adaptability recently gained interest for identifying more metabolically robust dairy cows. Metabolic challenges during periods of high milk production considerably affect the immune system, reproductive performance, and product quality as well as animal welfare. Moreover, growing consumer concerns need to be taken into consideration because the public perception of industrialized dairy cow farming, the high dependency on feed sources suitable for human nutrition, and the apparently abundant use of antibiotics may affect the sales of dairy products. Breeding for high yield continues, but the metabolic challenges increasingly come close to the adaptational limits of meeting the mammary gland's requirements. The aim of the present review is to elucidate metabolic challenges and adaptational limitations at different functional stages of the mammary gland in dairy cows. From the challenges and adaptational limitations, we derive perspectives for sustainable milk production. Based on previous research, we highlight the importance of metabolic plasticity in adaptation mechanisms at different functional stages of the mammary gland. Metabolic adaptation and plasticity change among developing, nonlactating, remodeling, and lactational stages of the mammary gland. A higher metabolic plasticity in early-lactating dairy cows could be indicative of resilience, and a high performance level without an extraordinary occurrence of health disorders can be achieved.

Methyl donor supplementation suppresses the progression of liver lipid accumulation while modifying the plasma triacylglycerol lipidome in periparturient Holstein dairy cows

Co-supplementation of methyl donors may lower hepatic lipid content in transition cows. To define the ability of methyl donor supplementation (MDS) to reduce hepatic lipid content and modify the plasma lipidome, 30 multiparous Holstein cows (2.04 ± 0.69 lactations; 689 ± 58 kg of body weight; 3.48 ± 0.10 units of body condition score) were fed a ration with or without rumen-protected methyl donors (22 g/d of Met, 10 g/d of choline chloride, 3 g/d of betaine, 96 mg/d of riboflavin, and 1.4 mg/d of vitamin B₁₂) from d -28 before expected calving through d 14 postpartum. Cows were randomly enrolled based on predefined selection criteria (body condition score and parity). Base diets without MDS were formulated for gestation (15.4% crude protein with a predicted Lys-to-Met ratio of 3.25; 1.44 Mcal of net energy for lactation/kg of dry matter) and lactation (16.6% crude protein with a predicted Lys-to-Met ratio of 3.36; 1.64 Mcal of net energy for lactation/kg of dry matter). Blood sampling occurred from d -28 relative to expected calving through d 14 postpartum. Liver tissue was biopsied at d -28 relative to expected calving and on d 5 and 14 postpartum. In addition to routine analyses, serum AA concentrations on d 10 and 12 were quantified using mass spectrometry. Plasma triacylglycerol (TAG) and cholesteryl esters (CE) were qualitatively measured using time-of-flight mass spectrometry. Data were analyzed using a mixed model with repeated measures. Dry matter intake and milk yield were not modified by MDS. The transition from d -28 relative to expected parturition to d 14 postpartum was characterized by increased plasma fatty acid (0.15 to 0.71 mmol/L) and β-hydroxybutyrate (0.34 to 0.43 mmol/L) levels and liver lipid content (3.91 to 9.16%). Methyl donor supplementation increased the serum Met level by 26% and decreased the serum Lys-to-Met ratio by 21% on d 10 and 12, respectively. Moreover, the increase in hepatic lipid content from d 5 through 14 postpartum was suppressed with MDS relative to control (3.57 vs. -0.29%). Dietary MDS modified the TAG and CE lipidome. For example, MDS increased plasma TAG 46:3 (carbon number:double bond) by 116% relative to control cows on d 5 postpartum. Moreover, MDS tended to increase plasma CE 34:6. In contrast, MDS lowered plasma TAG 54:8 by 39% relative to control cows on d 5 postpartum. We concluded that in the absence of gains in dry matter intake and milk and milk protein yields, dietary MDS slows the progression of hepatic lipid accumulation and modifies the plasma TAG lipidome in transition cows.

Genome-wide association study using haplotype alleles for the evaluation of reproductive traits in Nelore cattle

Zebu cattle (Bos taurus indicus) are highly adapted to tropical regions. However, females reach puberty after taurine heifers, which affects the economic efficiency of beef cattle breeding in the tropical regions. The aims of this study were to establish associations between the haplotype alleles of the bovine genome and age at first calving (AFC) in the Nelore cattle, and to identify the genes and quantitative trait loci (QTL) related to this phenotype. A total of 2,273 Nelore cattle (995 males and 1,278 females) genotyped using the Illumina BovineHD BeadChip were used in the current study. The association analysis included females with valid first calving records as well as open heifers. Linkage disequilibrium (LD) analysis among the markers was performed using blocks of 5, 10, and 15 markers, which were determined by sliding windows shifting one marker at a time. Then, the haplotype block size to be used in the association study was chosen based on the highest r² average among the SNPs in the block. The five HapAlleles most strongly associated with the trait (top five) were considered as significant associations. The results of the analysis revealed four genomic regions related to AFC, which overlapped with 20 QTL of the reproductive traits reported previously. Furthermore, there were 19 genes related to reproduction in those regions. In conclusion, the use of haplotypes allowed the detection of chromosomal regions associated with AFC in Nelore cattle, and provided the basis for elucidating the mechanisms underlying this trait.

Different mitochondrial DNA copy number in liver and mammary gland of lactating cows with divergent genetic background for milk production

Adequate metabolic adaptation of key tissues playing an essential role for bioenergetic homeostasis and lactogenesis is critical in cows to adapt to changes in energy requirements and physiological processes during the lactation period. Mitochondria are recognized as central to meet energy needs and maintaining of metabolic homeostasis because mitochondrial DNA (mtDNA) is template for several polypeptides of the respiratory chain complexes essential for ATP generation. The quantity of mtDNA in a cell has been widely used as a surrogate marker for the capacity of cells for energy generation. In our study we analyzed the mtDNA copy number and the mRNA expression of important nuclear encoded genes controlling mitochondrial biogenesis in liver and mammary gland. We compared cows with a nuclear genome dairy × beef crossbred make-up to purebred German Holstein dairy cows. The study revealed tissue-specific variations of mtDNA copy number and expression levels of nuclear genes involved in mitochondrial biogenesis when comparing lactating cows with different genetic predisposition regarding milk performance. This may reflect nuclear genome-determined genetic differences between the cow groups in coping with metabolic demands and physiological changes during lactation. The results indicate that mitochondrial biogenesis processes in the liver and mammary gland appear to be impaired in high lactating dairy cows, which consequently, would point to a disturbed energy adaptation. The results provide a basis to further elucidate the adaptive and regulatory modulation of the mitochondrial biogenesis in response to lactation-associated metabolic challenges in lactating cows.

Regulation of carnitine status in ruminants and efficacy of carnitine supplementation on performance and health aspects of ruminant livestock: a review

Carnitine has long been known to play a critical role for energy metabolism. Due to this, a large number of studies have been carried out to investigate the potential of supplemental carnitine in improving performance of livestock animals including ruminants, with however largely inconsistent results. An important issue that has to be considered when using carnitine as a feed additive is that the efficacy of supplemental carnitine is probably dependent on the animal's carnitine status, which is affected by endogenous carnitine synthesis, carnitine uptake from the gastrointestinal tract and carnitine excretion. The present review aims to summarise the current knowledge of the regulation of carnitine status and carnitine homeostasis in ruminants, and comprehensively evaluate the efficacy of carnitine supplementation on performance and/or health in ruminant livestock by comparing the outcomes of studies with carnitine supplementation in dairy cattle, growing and finishing cattle and sheep. While most of the studies show that supplemental carnitine, even in ruminally unprotected form, is bioavailable in ruminants, its effect on either milk or growth performance is largely disappointing. However, supplemental carnitine appears to be a useful strategy to offer protection against ammonia toxicity caused by consumption of high levels of non-protein N or forages with high levels of soluble N both, in cattle and sheep.

STAT6, PBX2, and PBRM1 Emerge as Predicted Regulators of 452 Differentially Expressed Genes Associated With Puberty in Brahman Heifers

The liver plays a central role in metabolism and produces important hormones. Hepatic estrogen receptors and the release of insulin-like growth factor 1 (IGF1) are critical links between liver function and the reproductive system. However, the role of liver in pubertal development is not fully understood. To explore this question, we applied transcriptomic analyses to liver samples of pre- and post-pubertal Brahman heifers and identified differentially expressed (DE) genes and genes encoding transcription factors (TFs). Differential expression of genes suggests potential biological mechanisms and pathways linking liver function to puberty. The analyses identified 452 DE genes and 82 TF with significant contribution to differential gene expression by using a regulatory impact factor metric. Brain-derived neurotrophic factor was observed as the most down-regulated gene (P = 0.003) in post-pubertal heifers and we propose this gene influences pubertal development in Brahman heifers. Additionally, co-expression network analysis provided evidence for three TF as key regulators of liver function during pubertal development: the signal transducer and activator of transcription 6, PBX homeobox 2, and polybromo 1. Pathway enrichment analysis identified transforming growth factor-beta and Wnt signaling pathways as significant annotation terms for the list of DE genes and TF in the co-expression network. Molecular information regarding genes and pathways described in this work are important to further our understanding of puberty onset in Brahman heifers.

Differential transcript abundance in adipose tissue of mature beef cows during feed restriction and realimentation

Feed costs account for over 70% of the annual expenditures in cow/calf production. During the production year the cow uses nutrients to support conceptus growth, milk production, work (grazing and locomotion), and maintenance requirements. The majority of the nutrients are used to support maintenance. Substrate cycling has been identified as one of the major contributors toward energy expenditure associated with maintenance in mature cows. The objective of this study was to determine whether beef cows that differ in the efficiency of weight gain differ in the relative abundance of transcripts for metabolic regulation in adipose tissue. Mature beef cows were subjected to feed restriction followed by ad libitum feed. Adipose tissue from twelve cows with high (n = 6) and low (n = 6) gain based on growth performance during the ad libitum feeding period was evaluated for transcriptome expression differences. A total of 496 genes were differentially expressed and passed Bonferroni correction for the animals with greater gain between restriction and realimentation and 491 genes were differentially expressed among animals with lesser gains between feed restriction and realimentation. Of these two differentially expressed gene lists, 144 genes were common between animals with greater and those with lesser gain. Enriched biological processes included the TCA cycle, oxidative phosphorylation, respiratory electron transport chain and fatty acid metabolic processes. Specific to adipose tissue of low gaining animals was glycolysis and to high gain animals was coenzyme, steroid, cellular amino acid, nitrogen compound metabolic processes, and sensory perception. The oxidative phosphorylation, mitochondrial dysfunction and cholesterol biosynthesis pathways were commonly associated with the high gain animals between feed restriction and realimentation, as well as with the low gaining animals between the two time points. Unique to the high gaining animals were valine degradation and LPS/IL-1 mediated inhibition of RXR function pathways. In this discovery study, genes involved in lipid metabolism, mitochondrial respiration and oxidative phosphorylation pathways appear to be critical to mature cows during times of abundant feed after feed restriction.

Integrated regulatory network reveals novel candidate regulators in the development of negative energy balance in cattle

Negative energy balance (NEB) is an altered metabolic state in modern high-yielding dairy cows. This metabolic state occurs in the early postpartum period when energy demands for milk production and maintenance exceed that of energy intake. Negative energy balance or poor adaptation to this metabolic state has important effects on the liver and can lead to metabolic disorders and reduced fertility. The roles of regulatory factors, including transcription factors (TFs) and micro RNAs (miRNAs) have often been separately studied for evaluating of NEB. However, adaptive response to NEB is controlled by complex gene networks and still not fully understood. In this study, we aimed to discover the integrated gene regulatory networks involved in NEB development in liver tissue. We downloaded data sets including mRNA and miRNA expression profiles related to three and four cows with severe and moderate NEB, respectively. Our method integrated two independent types of information: module inference network by TFs, miRNAs and mRNA expression profiles (RNA-seq data) and computational target predictions. In total, 176 modules were predicted by using gene expression data and 64 miRNAs and 63 TFs were assigned to these modules. By using our integrated computational approach, we identified 13 TF-module and 19 miRNA-module interactions. Most of these modules were associated with liver metabolic processes as well as immune and stress responses, which might play crucial roles in NEB development. Literature survey results also showed that several regulators and gene targets have already been characterized as important factors in liver metabolic processes. These results provided novel insights into regulatory mechanisms at the TF and miRNA levels during NEB. In addition, the method described in this study seems to be applicable to construct integrated regulatory networks for different diseases or disorders.

Expression of hepatic genes related to energy metabolism during the transition period of Holstein and F1 Holstein-Gir cows

The aim of this study was to investigate the expression of genes encoding enzymes and other factors involved with carbohydrate and lipid metabolism in the liver of 2 genetic groups of dairy cows during the transition period. We analyzed the expression of glucose-6-phosphatase (G6PC), cytosolic phosphoenolpyruvate carboxykinase (PEPCK-C), methylmalonyl-CoA mutase (MUT), β-hydroxybutyrate dehydrogenase-2 (BDH2), acetyl-CoA carboxylase (ACC), carnitine palmitoyltransferase-2 (CPT2), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), glucose transporter-2 (SLC2A2), and the transcription factor peroxisome proliferator-activated receptor α (PPARA). Blood concentrations of glucose, nonesterified fatty acids, and β-hydroxybutyrate were also determined. Liver biopsies and blood samples were taken at d 15 prepartum and at d 6, 21, 36, 51, and 66 postpartum from Holsteins (n = 6) and F₁ Holstein-Gir (n = 6) cows. Cows were kept under the same prepartum and postpartum management conditions. The results showed that the expression of G6PC, PEPCK-C, BDH2, ACC, CPT2, HMGCR, SLC2A2, and PPARA genes did not differ between genetic groups. Except for PEPCK-C, no interaction between genetic groups and the experimental period was observed. Within both groups of cows, G6PC and PEPCK-C gene expression decreased when comparing prepartum gene expression with 21 and 36 DIM, and increased in d 51 postpartum. MUT mRNA levels differed between the 2 genetic groups and displayed a significant increase after d 36 postpartum, whereas mRNA levels of HMGCR tended to increase when comparing d 21 and 36 to d 51 postpartum. Glucose concentrations also differed between genetic groups, being significantly higher in the plasma of F₁ Holstein-Gir cows than in Holstein cows, but no differences were found within each group during the analysis period. β-Hydroxybutyrate and nonesterified fatty acid concentrations did not differ between genetic groups, but displayed increased levels from prepartum to d 6 and 21 postpartum. Our results indicated that expression in the liver of genes involved with glucose and fatty acid metabolism were similar in both groups of cows and significant differences were observed between the 2 groups in the expression of MUT, a gene involved in propionate metabolism.

Analysis of hepatic transcript profile and plasma lipid profile in early lactating dairy cows fed grape seed and grape marc meal extract

Background

It was recently reported that dairy cows fed a polyphenol-rich grape seed and grape marc meal extract (GSGME) during the transition period had an increased milk yield, but the underlying reasons remained unclear. As polyphenols exert a broad spectrum of metabolic effects, we hypothesized that feeding of GSGME influences metabolic pathways in the liver which could account for the positive effects of GSGME in dairy cows. In order to identify these pathways, we performed genome-wide transcript profiling in the liver and lipid profiling in plasma of dairy cows fed GSGME during the transition period at 1 week postpartum.

Results

Transcriptomic analysis of the liver revealed 207 differentially expressed transcripts, from which 156 were up- and 51 were down-regulated, between cows fed GSGME and control cows. Gene set enrichment analysis of the 155 up-regulated mRNAs showed that the most enriched gene ontology (GO) biological process terms were dealing with cell cycle regulation and the most enriched Kyoto Encyclopedia of Genes and Genomes pathways were p53 signaling and cell cycle. Functional analysis of the 43 down-regulated mRNAs revealed that a great part of these genes are involved in endoplasmic reticulum (ER) stress-induced unfolded protein response (UPR) and inflammatory processes. Accordingly, protein folding, response to unfolded protein, unfolded protein binding, chemokine activity and heat shock protein binding were identified as one of the most enriched GO biological process and molecular function terms assigned to the down-regulated genes. In line with the transcriptomics data the plasma concentrations of the acute phase proteins serum amyloid A (SAA) and haptoglobin were reduced in cows fed GSGME compared to control cows. Lipidomic analysis of plasma revealed no differences in the concentrations of individual species of major and minor lipid classes between cows fed GSGME and control cows.

Conclusions

Analysis of hepatic transcript profile in cows fed GSGME during the transition period at 1 week postpartum indicates that polyphenol-rich feed components are able to inhibit ER stress-induced UPR and inflammatory processes, both of which are considered to contribute to liver-associated diseases and to impair milk performance in dairy cows, in the liver of dairy cows during early lactation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3638-1) contains supplementary material, which is available to authorized users.

Transcriptomics of liver and muscle in Holstein cows genetically divergent for fertility highlight differences in nutrient partitioning and inflammation processes

Background

The transition between pregnancy and lactation is a major physiological change for dairy cows. Complex systemic and local processes involving regulation of energy balance, galactopoiesis, utilisation of body reserves, insulin resistance, resumption of oestrous cyclicity and involution of the uterus can affect animal productivity and hence farm profitability. Here we used an established Holstein dairy cow model of fertility that displayed genetic and phenotypic divergence in calving interval. Cows had similar genetic merit for milk production traits, but either very good genetic merit for fertility traits (‘Fert+’; n = 8) or very poor genetic merit for fertility traits (‘Fert-’; n = 8). We used RNA sequencing to investigate gene expression profiles in both liver and muscle tissue biopsies at three distinct time-points: late pregnancy, early lactation and mid lactation (-18, 1 and 147 days relative to parturition, respectively).

Results

We found 807 and 815 unique genes to be differentially expressed in at least one time-point in liver and muscle respectively, of which 79 % and 83 % were only found in a single time-point; 40 and 41 genes were found differentially expressed at every time-point indicating possible systemic or chronic dysregulation. Functional annotation of all differentially expressed genes highlighted two physiological processes that were impacted at every time-point in the study, These were immune and inflammation, and metabolic, lipid and carbohydrate-binding.

Conclusion

These pathways have previously been identified by other researchers. We show that several specific genes which are differentially regulated, including IGF-1, might impact dairy fertility. We postulate that an increased burden of reactive oxidation species, coupled with a chronic inflammatory state, might reduce dairy cow fertility in our model.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2938-1) contains supplementary material, which is available to authorized users.