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

Multiomics reveals blood differential metabolites and differential genes in the early onset of ketosis in dairy cows

Ketosis-a metabolic state characterized by elevated levels of ketone bodies in the blood or urine-reduces the performance and health of dairy cows and causes substantial economic losses for the dairy industry. Currently, beta-hydroxybutyric acid is the gold standard for determining ketosis in cows; however, as this method is only applicable postpartum, it is not conducive to the early intervention of ketosis in dairy cows. In this study, the sera of dry, periparturient, postpartum ketotic, and healthy cows were analyzed by both transcriptomics and metabolomics techniques. Moreover, changes of gene expression and metabolites were observed, and serum physiological and biochemical indexes were detected by ELISA. The purpose was to screen biomarkers that can be used to detect the incidence of dry or periparturient ketosis in cows. The results showed that ketotic cows had increased levels of glycolipid metabolism indexes, oxidizing factors, and inflammatory factors during dry periods and liver damage, which could be used as early biomarkers to predict the onset of ketosis. Transcriptomic results yielded 20 differentially expressed genes (DEGs) between ketotic and healthy cows during dry, peripartum, and postpartum periods. GO and KEGG enrichment analyses indicated that these DEGs were involved in amino acid metabolism, energy metabolism, and disease-related signaling pathways. The metabolomics sequencing results showed that ketotic cows mainly showed enrichment in tricarboxylic acid cycle, butyric acid metabolism, carbon metabolism, lysine degradation, fatty acid degradation, and other signaling pathways. Metabolites differed between ketotic and healthy cows in dry, pre-parturition, and post-parturition periods. Combined transcriptomics and metabolomics analyses identified significant enrichment in the glucagon signaling pathway and the lysine degradation signaling pathway in dry, periparturient, and postpartum ketotic cows. PRKAB2 and SETMAR-key DEGs of the glucagon signaling pathway and lysine degradation signaling pathway, respectively-can be used as key marker genes for determining the early onset of ketosis in dairy cows.

Multi-Omics Reveals Disrupted Immunometabolic Homeostasis and Oxidative Stress in Adipose Tissue of Dairy Cows with Subclinical Ketosis: A Sphingolipid-Centric Perspective

Ketosis, especially its subclinical form, is frequently observed in high-yielding dairy cows and is linked to various diseases during the transition period. Although adipose tissue plays a significant role in the development of metabolic disorders, its exact impact on the emergence of subclinical ketosis (SCK) is still poorly understood. The objectives of this study were to characterize and compare the profiling of transcriptome and lipidome of blood and adipose tissue between SCK and healthy cows and investigate the potential correlation between metabolic disorders and lipid metabolism. We obtained blood and adipose tissue samples from healthy cows (CON, n = 8, β-hydroxybutyric acid concentration < 1.2 mmol/L) and subclinical ketotic cows (SCK, n = 8, β-hydroxybutyric acid concentration = 1.2-3.0 mmol/L) for analyzing biochemical parameters, transcriptome, and lipidome. We found that serum levels of nonesterified fatty acids, malonaldehyde, serum amyloid A protein, IL-1β, and IL-6 were higher in SCK cows than in CON cows. Levels of adiponectin and total antioxidant capacity were higher in serum and adipose tissue from SCK cows than in CON cows. The top enriched pathways in whole blood and adipose tissue were associated with immune and inflammatory responses and sphingolipid metabolism, respectively. The accumulation of ceramide and sphingomyelin in adipose tissue was paralleled by an increase in genes related to ceramide biosynthesis, lipolysis, and inflammation and a decrease in genes related to ceramide catabolism, lipogenesis, adiponectin production, and antioxidant enzyme systems. Increased ceramide concentrations in blood and adipose tissue correlated with reduced insulin sensitivity. The current results indicate that the lipid profile of blood and adipose tissue is altered with SCK and that certain ceramide species correlate with metabolic health. Our research suggests that disruptions in ceramide metabolism could be crucial in the progression of SCK, exacerbating conditions such as insulin resistance, increased lipolysis, inflammation, and oxidative stress, providing a potential biomarker of SCK and a novel target for nutritional manipulation and pharmacological therapy.

Changes in blood biochemical parameters in highly productive cows with ketosis

Background and Aim

Biochemical blood testing is the main diagnostic indicator of the clinical condition of highly productive animals and a method of determining changes in metabolic disorders. This study focuses on metabolic changes (ketosis), which are of the utmost importance in the assessment of the health status of animals, as well as differences in intergroup characteristics. The main focus of this study is to demonstrate the influence of subclinical ketosis in highly productive cows on changes in biochemical blood parameters during different physiological periods to further prevent this disease, adjust feeding rations, and prevent premature culling of animals. This study aimed to evaluate and establish changes in the biochemical status dynamics of highly productive cows with metabolic disorders in an industrial livestock complex.

Materials and Methods

Blood samples were systematically collected from highly productive cows of the Simmental breed (n = 60) and served as the primary material for subsequent analyses. Each methodological step was designed to ensure evaluation of the metabolic changes associated with post-calving adjustments in highly productive dairy cows. This study employed a comprehensive approach integrating clinical assessments, laboratory analyses, biochemical evaluations, instrumental measurements, and statistical analyses.

Results

A biochemical blood test showed that the number of ketone bodies in the experimental group exceeded the norm, varied depending on the physiological state of the animals, and ranged from 0.89 to 1.45 mmol/L. At 10 days after calving, the highest indicator was 1.45 ± 0.05 mmol/L. This indicator was 1.05 mmol/L higher than that in the control group and exceeded the norm by 0.95.

Conclusion

Excess ketone bodies in the blood of animals led to accumulation in urine and milk, indicating a disturbance in metabolic processes in the body and a decrease in the quality of animal husbandry products. The sample size and the focus on a single breed from one geographical location may limit the generalizability of the findings. Further research should explore the mechanistic bases of ketosis development, potentially integrating genomic and proteomic approaches to understand the genetic predispositions and molecular pathways involved.

Alterations in the gut microbiota and its metabolites contribute to metabolic maladaptation in dairy cows during the development of hyperketonemia

Metabolic maladaptation in dairy cows after calving can lead to long-term elevation of ketones, such as β-hydroxybutyrate (BHB), representing the condition known as hyperketonemia, which greatly influences the health and production performance of cows during the lactation period. Although the gut microbiota is known to alter in dairy cows with hyperketonemia, the association of microbial metabolites with development of hyperketonemia remains unknown. In this study, we performed a multi-omics analysis to investigate the associations between fecal microbial community, fecal/plasma metabolites, and serum markers in hyperketonemic dairy cows during the transition period. Dynamic changes in the abundance of the phyla Verrucomicrobiota and Proteobacteria were detected in the gut microbiota of dairy cows, representing an adaptation to enhanced lipolysis and abnormal glucose metabolism after calving. Random forest and univariate analyses indicated that Frisingicoccus is a key bacterial genus in the gut of cows during the development of hyperketonemia, and its abundance was positively correlated with circulating branched-chain amino acid levels and the ketogenesis pathway. Taurodeoxycholic acid, belonging to the microbial metabolite, was strongly correlated with an increase in blood BHB level, and the levels of other secondary bile acid in the feces and plasma were altered in dairy cows prior to the diagnosis of hyperketonemia, which link the gut microbiota and hyperketonemia. Our results suggest that alterations in the gut microbiota and its metabolites contribute to excessive lipolysis and insulin insensitivity during the development of hyperketonemia, providing fundamental knowledge about manipulation of gut microbiome to improve metabolic adaptability in transition dairy cows.IMPORTANCEAccumulating evidence is pointing to an important association between gut microbiota-derived metabolites and metabolic disorders in humans and animals; however, this association in dairy cows from late gestation to early lactation is poorly understood. To address this gap, we integrated longitudinal gut microbial (feces) and metabolic (feces and plasma) profiles to characterize the phenotypic differences between healthy and hyperketonemic dairy cows from late gestation to early lactation. Our results demonstrate that cows underwent excessive lipid mobilization and insulin insensitivity before hyperketonemia was evident. The bile acids are functional readouts that link gut microbiota and host phenotypes in the development of hyperketonemia. Thus, this work provides new insight into the mechanisms involved in metabolic adaptation during the transition period to adjust to the high energy and metabolic demands after calving and during lactation, which can offer new strategies for livestock management involving intervention of the gut microbiome to facilitate metabolic adaptation.

Severe body condition loss lowers hepatic output of IGF1 with adverse effects on the dominant follicle in dairy cows

The severe loss of body condition score (BCS) during the early lactation period has been associated with infertility in cows. However, the mechanisms are not fully understood. The aim of this study was to examine the effect of BCS loss on liver health, and ovarian functions in cows during early lactation. Retrospectively multiparous cows from two farms were categorized based on units of BCS (1-5 scale) loss as Moderate (MOD, <0.75 units; n = 11) or Severe (SEV, ≥0.75 units; n = 9) loss groups. From Weeks -3 to 7, relative to calving, MOD and SEV cows lost on average 0.4 and 1.0-unit BCS, respectively. All data except hepatic transcriptomes were analyzed with PROC MIXED procedure of SAS. The plasma concentration of non-esterified fatty acids at Week 0 and 1, ß-hydroxy butyrate at Week 1, and γ-glutamyl transferase at Weeks 1 and 7 relative to calving were higher in SEV cows. Hepatic transcriptome analysis showed that 1 186 genes were differentially expressed in SEV (n = 3) compared to MOD (n = 3) cows at Week 7 after calving. Pathway analysis revealed that significant DEGs in SEV cows enriched in lipid metabolisms including, lipid metabolic process, ether lipid metabolism, fatty acid beta-oxidation, fatty acid biosynthetic process, fatty acid metabolic process, fat digestion and absorption, linoleic acid metabolism, alpha-linolenic acid metabolism. The impaired liver function in SEV cows was associated with 1.5-fold reduction of hepatic IGF1 gene expression and lower serum IGF1 concentrations. At the ovarian level, SEV cows had lower IGF1 concentration in the follicular fluid of the dominant follicle of the synchronized follicular wave compared to that of MOD cows at 7 weeks after calving. Further, the follicular fluid concentration of estradiol-17β was lower in SEV cows along with lower transcript abundance of genes from granulosa cells associated with dominant follicle competence, including CYP19A1, NR5A2, IGF1, and LHCGR. These data show that SEV loss of BCS during early lactation leading up to the planned start of breeding is associated with liver dysfunction, including lower IGF1 secretion, and impaired function of the dominant follicle in the ovary.

Ceramide on the road to insulin resistance and immunometabolic disorders in transition dairy cows: driver or passenger?

Dairy cows must undergo profound metabolic and endocrine adaptations during their transition period to meet the nutrient requirements of the developing fetus, parturition, and the onset of lactation. Insulin resistance in extrahepatic tissues is a critical component of homeorhetic adaptations in periparturient dairy cows. However, due to increased energy demands at calving that are not followed by a concomitant increase in dry matter intake, body stores are mobilized, and the risk of metabolic disorders dramatically increases. Sphingolipid ceramides involved in multiple vital biological processes, such as proliferation, differentiation, apoptosis, and inflammation. Three typical pathways generate ceramide, and many factors contribute to its production as part of the cell's stress response. Based on lipidomic profiling, there has generally been an association between increased ceramide content and various disease outcomes in rodents. Emerging evidence shows that ceramides might play crucial roles in the adaptive metabolic alterations accompanying the initiation of lactation in dairy cows. A series of studies also revealed a negative association between circulating ceramides and systemic insulin sensitivity in dairy cows experiencing severe negative energy balance. Whether ceramide acts as a driver or passenger in the metabolic stress of periparturient dairy cows is an unknown but exciting topic. In the present review, we discuss the potential roles of ceramides in various metabolic dysfunctions and the impacts of their perturbations. We also discuss how this novel class of bioactive sphingolipids has drawn interest in extrahepatic tissue insulin resistance and immunometabolic disorders in transition dairy cows. We also discuss the possible use of ceramide as a new biomarker for predicting metabolic diseases in cows and highlight the remaining problems.

Free Fatty Acids Induce Apoptosis of Mammary Epithelial Cells of Ketotic Dairy Cows via the Mito-ROS/NLRP3 Signaling Pathway

Early lactation increases metabolic stress in ketotic dairy cows, leading to mitochondrial damage, apoptosis, and inflammatory response in mammary epithelial cells. The pyrin domain 3 (NLRP3) pathway involving the mitochondrial reactive oxygen species (Mito-ROS)-induced nucleotide-binding oligomerization domain-like receptor has been recognized as a key mechanism in this inflammatory response and cell apoptosis. This study aimed to elucidate the underlying regulatory mechanism of Mito-ROS-NLRP3 pathway-mediated mammary epithelial cell apoptosis in dairy cows with ketosis. Mitochondrial damage and cellular apoptotic program and NLRP3 inflammasome activation were observed in the mammary gland of ketotic cows. Similar damage was detected in MAC-T cells treated with exogenous fatty acids (FFAs). However, NLRP3 inhibitor MCC950 pretreatment or Mito-ROS scavenging by MitoTEMPO attenuated apoptosis in FFA-induced MAC-T cells by inhibiting the NLRP3 inflammasome pathway. These findings reveal that the Mito-ROS-NLRP3 pathway activation is a potent mechanism underlying mammary epithelial cell apoptosis in response to metabolic stress in ketotic dairy cows, which further contributes to reduced milk yield.

Influence of Cobalt Source, Folic Acid, and Rumen-Protected Methionine on Performance, Metabolism, and Liver Tissue One-Carbon Metabolism Biomarkers in Peripartal Holstein Cows

Vitamin B₁₂ plays a role in the remethylation of homocysteine to Met, which then serves as a substrate for Met adenosyltransferase (MAT) to synthesize S-adenosylmethionine (SAM). We investigated effects of feeding two cobalt sources [Co-glucoheptonate (CoPro) or CoPectin, Zinpro Corp.], an experimental ruminally-available source of folic acid (FOA), and rumen-protected Met (RPM) on performance and hepatic one-carbon metabolism in peripartal Holstein cows. From -30 to 30 d around calving, 72 multiparous cows were randomly allocated to: CoPro, CoPro + FOA, CoPectin + FOA, or CoPectin + FOA + RPM. The Co treatments delivered 1 mg Co/kg of DM (CoPro or CoPectin), each FOA group received 50 mg/d FOA, and RPM was fed at 0.09% of DM intake (DMI). Milk yield and DMI were not affected. Compared with other groups, the percentage of milk protein was greater after the second week of lactation in CoPectin + FOA + RPM. Compared with CoPro or CoPro + FOA, feeding CoPectin + FOA or CoPectin + FOA + RPM led to a greater activity of MAT at 7 to 15 d postcalving. For betaine-homocysteine S-methyltransferase, CoPro together with CoPectin + FOA + RPM cows had greater activity at 7 and 15 d than CoPro + FOA. Overall, supplying FOA with CoPectin or CoPectin plus RPM may enhance S-adenosylmethionine synthesis via MAT in the liver after parturition. As such, these nutrients may impact methylation reactions and liver function.

Stress and inflammatory response of cows and their calves during peripartum and early neonatal period

Stress, inflammatory response, and their relationship were investigated in Simmental cows during the transition period (N = 8; 5 multiparous and 3 primiparous) and in their calves (N = 8; 5 heifers and 3 bulls). From cows, blood was collected at days -21 (±4), 0, +1, +7, and +21 days relative to calving. From calves, blood was collected after birth before colostrum intake (0) and then at 1, 7, and 15 days of age. Cortisol, Interleukin 6 (IL-6) and haptoglobin concentration was assessed by ELISA technique; white blood cells (WBC) were assessed using an ADVIA 2120 Hematology System machine. One-way ANOVA showed an effect of time for all the investigated parameters (P < 0.001) except for lymphocytes in peripartal cows. At calving and 1 d after, cortisol concentration was negatively correlated with levels of IL-6, WBC, and monocytes, whereas levels of IL-6 were positively correlated with WBC, neutrophils, and monocytes count. Cortisol, IL-6, haptoglobin, WBC and all leukocyte populations were affected by the age of neonatal calves (P < 0.001) except for neutrophils. A negative correlation between cortisol and IL-6, neutrophils, monocytes and haptoglobin was found at 15 days of age. A positive correlation between IL-6 and haptoglobin at day 15 of age, and with neutrophils and monocytes at days 7 and 15 of age was found. A positive correlation was obtained between cortisol levels measured in cows around calving and those obtained in calves after birth before colostrum intake (r = 0.83), and between IL-6 concentrations obtained from cows at calving and 1 d after and those obtained in calves at day 1 of age, after the colostrum intake (r = 0.93 and 0.79, respectively). The study suggests that immune function of peripartal cows is in an active state and that, in addition to other well-known factors driving the changes of parameters herein investigated, cortisol could have a role in the immune-modulatory adjustment during peripartum in cows. Furthermore, it can be hypothesized that cortisol is transferred from the cow to newborn calf through the placenta only and not through colostrum, whereas IL-6 levels in calves during the 24 h after birth seem to be influenced by IL-6 values measured in cows around calving due to its transfer through colostrum.

Identification of the metabolic state of surviving cardiomyocytes in the human infarcted heart by spatial single-cell transcriptomics

The metabolic status of surviving cardiomyocytes (CM) in the myocardial tissues of patients who sustained myocardial infarction (MI) is largely unknown. Spatial single-cell RNA-sequencing (scRNA-seq) is a novel tool that enables the unbiased analysis of RNA signatures within intact tissues. We employed this tool to assess the metabolic profiles of surviving CM in the myocardial tissues of patients post-MI.

Methods

A spatial scRNA-seq dataset was used to compare the genetic profiles of CM from patients with MI and control patients; we analyzed the metabolic adaptations of surviving CM within the ischemic niche. A standard pipeline in Seurat was used for data analysis, including normalization, feature selection, and identification of highly variable genes using principal component analysis (PCA). Harmony was used to remove batch effects and integrate the CM samples based on annotations. Uniform manifold approximation and projection (UMAP) was used for dimensional reduction. The Seurat "FindMarkers" function was used to identify differentially expressed genes (DEGs), which were analyzed by the Gene Ontology (GO) enrichment pathway. Finally, the scMetabolism R tool pipeline with parameters method = VISION (Vision is a flexible system that utilizes a high-throughput pipeline and an interactive web-based report to annotate and explore scRNA-seq datasets in a dynamic manner) and metabolism.type = Kyoto Encyclopedia of Genes and Genomes (KEGG) was used to quantify the metabolic activity of each CM.

Results

Analysis of spatial scRNA-seq data showed fewer surviving CM in infarcted hearts than in control hearts. GO analysis revealed repressed pathways in oxidative phosphorylation, cardiac cell development, and activated pathways in response to stimuli and macromolecular metabolic processes. Metabolic analysis showed downregulated energy and amino acid pathways and increased purine, pyrimidine, and one-carbon pool by folate pathways in surviving CM.

Conclusions

Surviving CM within the infarcted myocardium exhibited metabolic adaptations, as evidenced by the downregulation of most pathways linked to oxidative phosphorylation, glucose, fatty acid, and amino acid metabolism. In contrast, pathways linked to purine and pyrimidine metabolism, fatty acid biosynthesis, and one-carbon metabolism were upregulated in surviving CM. These novel findings have implications for the development of effective strategies to improve the survival of hibernating CM within the infarcted heart.

Changes of the liver metabolome following an intravenous lipopolysaccharide injection in Holstein cows supplemented with dietary carnitine

Background

Carnitine facilitates the flux of long-chain fatty acids for hepatic mitochondrial beta-oxidation, which acts to ameliorate the negative energy balance commonly affecting high-yielding dairy cows. Inflammation triggered by lipopolysaccharide (LPS) load can however pose a challenge to the metabolic integrity via the expression of pro-inflammatory mediators, leading to immune system activation and respective metabolic alterations. The effect of enhanced carnitine availability on hepatic metabolome profiles during an inflammatory challenge has not yet been determined in dairy cows. Herein, Holstein cows were supplemented with 25 g/d rumen-protected carnitine from 42 d prepartum until 126 d postpartum (n = 16) or assigned to the control group with no supplementation during the same period (n = 14). We biopsied the liver of the cows before (100 d postpartum) and after (112 d postpartum) an intravenous injection of 0.5 µg/kg LPS. Liver samples were subjected to a targeted metabolomics analysis using the AbsoluteIDQ p180 Kit (Biocrates Life Sciences AG, Innsbruck, Austria). 

Results

Multivariate statistical analyses revealed that hepatic metabolome profiles changed in relation to both the carnitine supplementation and the LPS challenge. Comparing the metabolite profiles on 100 d, carnitine increased the concentration of short- and long-chain acyl-carnitines, which may be explained by an enhanced mitochondrial fatty acid shuttle and hence greater energy availability. The LPS injection affected hepatic metabolite profiles only in the carnitine supplemented group, particularly altering the concentration of biogenic amines.

Conclusions

Our results point to interactions between an acute hepatic inflammatory response and biogenic amine metabolism, depending on energy availability.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00741-z.

Activated autophagy-lysosomal pathway in dairy cows with hyperketonemia is associated with lipolysis of adipose tissues

Activated autophagy-lysosomal pathway (ALP) can degrade virtually all kinds of cellular components, including intracellular lipid droplets, especially during catabolic conditions. Sustained lipolysis and increased plasma fatty acids concentrations are characteristic of dairy cows with hyperketonemia. However, the status of ALP in adipose tissue during this physiological condition is not well known. The present study aimed to ascertain whether lipolysis is associated with activation of ALP in adipose tissues of dairy cows with hyperketonemia and in calf adipocytes. In vivo, blood and subcutaneous adipose tissue (SAT) biopsies were collected from nonhyperketonemic (nonHYK) cows [blood β-hydroxybutyrate (BHB) concentration <1.2 mM, n = 10] and hyperketonemic (HYK) cows (blood BHB concentration 1.2-3.0 mM, n = 10) with similar days in milk (range: 3-9) and parity (range: 2-4). In vitro, calf adipocytes isolated from 5 healthy Holstein calves (1 d old, female, 30-40 kg) were differentiated and used for (1) treatment with lipolysis inducer isoproterenol (ISO, 10 µM, 3 h) or mammalian target of rapamycin inhibitor Torin1 (250 nM, 3 h), and (2) pretreatment with or without the ALP inhibitor leupeptin (10 μg/mL, 4 h) followed by ISO (10 µM, 3 h) treatment. Compared with nonHYK cows, serum concentration of free fatty acids was greater and serum glucose concentration, DMI, and milk yield were lower in HYK cows. In SAT of HYK cows, ratio of phosphorylated hormone-sensitive lipase to hormone-sensitive lipase, and protein abundance of adipose triacylglycerol lipase were greater, but protein abundance of perilipin 1 (PLIN1) and cell death-inducing DNA fragmentation factor-α-like effector c (CIDEC) was lower. In addition, mRNA abundance of autophagy-related 5 (ATG5), autophagy-related 7 (ATG7), and microtubule-associated protein 1 light chain 3 beta (MAP1LC3B), protein abundance of lysosome-associated membrane protein 1, and cathepsin D, and activity of β-N-acetylglucosaminidase were greater, whereas protein abundance of sequestosome-1 (p62) was lower in SAT of HYK cows. In calf adipocytes, treatment with ISO or Torin1 decreased protein abundance of PLIN1, and CIDEC, and triacylglycerol content in calf adipocytes, but increased glycerol content in the supernatant of calf adipocytes. Moreover, the mRNA abundance of ATG5, ATG7, and MAP1LC3B was upregulated, the protein abundance of lysosome-associated membrane protein 1, cathepsin D, and activity of β-N-acetylglucosaminidase were increased, whereas the protein abundance of p62 was decreased in calf adipocytes treated with ISO or Torin1 compared with control group. Compared with treatment with ISO alone, the protein abundance of p62, PLIN1, and CIDEC, and triacylglycerol content in calf adipocytes were higher, but the glycerol content in the supernatant of calf adipocytes was lower in ISO and leupeptin co-treated group. Overall, these data indicated that activated ALP is associated with increased lipolysis in adipose tissues of dairy cows with hyperketonemia and in calf adipocytes.

Transcriptome analysis reveals disruption of circadian rhythms in late gestation dairy cows may increase risk for fatty liver and reduced mammary remodeling

Circadian disruption increased insulin resistance and decreased mammary development in late gestation, nonlactating (dry) cows. The objective was to measure the effect of circadian disruption on transcriptomes of the liver and mammary gland. At 35 days before expected calving (BEC), multiparous dry cows were assigned to either control (CON) or phase-shifted treatments (PS). CON was exposed to 16-h light and 8-h dark. PS was exposed to 16-h light to 8-h dark, but phase of the light-dark cycle was shifted 6 h every 3 days. On day 21 BEC, liver and mammary were biopsied. RNA was isolated (n = 6 CON, n = 6 PS per tissue), and libraries were prepared and sequenced using paired-end reads. Reads mapping to bovine genome averaged 27 ± 2 million and aligned to 14,222 protein-coding genes in liver and 15,480 in mammary analysis. In the liver, 834 genes, and in the mammary gland, 862 genes were different (nominal P < 0.05) between PS and CON. In the liver, genes upregulated in PS functioned in cholesterol biosynthesis, endoplasmic reticulum stress, wound healing, and inflammation. Genes downregulated in liver function in cholesterol efflux. In the mammary gland, genes upregulated functioned in mRNA processing and transcription and downregulated genes encoded extracellular matrix proteins and proteases, cathepsins and lysosomal proteases, lipid transporters, and regulated oxidative phosphorylation. Increased cholesterol synthesis and decreased efflux suggest that circadian disruption potentially increases the risk of fatty liver in cows. Decreased remodeling and lipid transport in mammary may decrease milk production capacity during lactation.

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.

Aspects of transition cow metabolomics-Part III: Alterations in the metabolome of liver and blood throughout the transition period in cows with different liver metabotypes

The liver plays a central role in the postpartum (PP) energy metabolism of the transition dairy cow; however, studies describing the liver metabolome during this period were lacking. The aim of the presented study was therefore to compare the alterations in the liver and blood metabolome of transition dairy cows. For this purpose, an on-farm trial with 80 German Holstein cows (mean lactation number: 3.9; range: 2-9) was performed, with thorough documentation of clinical traits and clinical chemistry, as well as production data. Liver biopsies and blood samples were collected at d 14 (mean: 12 d, range: 1-26 d) antepartum (AP), d 7 (7, 4-13) and 28 (28, 23-34; mean, earliest-latest) PP for targeted mass spectroscopy-based metabolomics analysis using the AbsoluteIDQ p180 kit (Biocrates Life Sciences). Statistical analysis was performed using multivariate (partial least squares discriminant analysis) as well as univariate methods (linear mixed model). Multivariate data analysis of the liver metabolome revealed 3 different metabotypes (A = medium, B = minor, C = large alterations in the liver metabolome profile between AP and PP). In metabotype C, an increase of almost all acylcarnitines, lysophosphatidylcholines (lysoPC), sphingomyelins, and some phosphatidylcholines (PC, mainly at 7 d PP) was observed after calving. In contrast to metabotype C, the clinical data of the metabotype B animals indicated a higher PP lipomobilization and occurrence of transition cow diseases. The liver metabolome profile of these animals most likely mirrors a failure of adaptation to the PP state. This strong occurrence of metabotypes was much less pronounced in the blood metabolome. Additionally, differences in metabolic patterns were observed across the transition period when comparing liver and blood matrices (e.g., in different biogenic amines, acylcarnitines and sphingolipids). In summary, the blood samples at 7 d PP showed lower acylcarnitines and PC, with minor alterations and a heterogeneous pattern in AA, biogenic amines, and sphingomyelins compared with 14 d AP. In contrast to 7 d PP, the blood samples at 28 PP revealed an increase in several AA, lysoPC, PC, and sphingomyelins in comparison to the AP state, irrespective of the metabotype. In the liver biopsies metabotype B differed from metabotype C animals ante partum by following metabolites: higher α aminoadipic acid, lower AA, serotonin, taurine, and symmetric dimethylarginine levels, lower or higher concentrations of certain acylcarnitines (higher: C2, C3, C5, C4:1; lower: C12:1, C14:1-OH, C16:2), and lower lysoPC (a C16:0, C18:0, C20:3, C20:4) and hexose levels. In blood samples, fewer differences were observed, with lower serotonin, acylcarnitine C16:2, lysoPC (a C16:0, C17:0, C18:0 and C18:1), PC aa C38:0, and PC ae C42:2. The results show that the use of only the blood metabolome to assess liver metabolism may be hampered by the fact that blood profiles are influenced by the metabolism of many organs, and metabolomics analysis from liver biopsies is a more suitable method to identify distinct metabotypes. Future studies should investigate the stability and reproducibility of the metabotype and phenotypes observed, and the possible predictive value of the metabolites already differing AP between metabotype B and C.

Effect of overconditioning on the hepatic global gene expression pattern of dairy cows at the end of pregnancy

Overconditioning is a risk factor for upregulated pre- and postpartum fat mobilization. Therefore, we hypothesized that overconditioning at the end of pregnancy leads to the accumulation of lipids in the liver and modifications of the hepatic gene expression pattern. The aim of this study was to evaluate the effect of normal- versus overconditioning on the hepatic transcriptomic profile of dairy cows at the end of pregnancy. Ten dry multiparous Holstein cows were killed 2 wk before expected calving. Body condition score (BCS) and backfat thickness (BFT) were evaluated, and blood samples for nonesterified fatty acids (NEFA) were taken before cows were killed. After cows were killed, liver biopsy samples were collected for further assessment of total lipids and RNA sequencing. Five cows were classified as normal-conditioned (median BCS = 3, range 2.75-3.5) and 5 as overconditioned (median BCS = 4, range 4-5). Regression models confirmed that normal-conditioned cows had lower BFT (1.29 ± 0.29 cm; least squares means ± standard error) and serum NEFA (0.16 ± 0.04 mmol/L) in comparison to overconditioned cows (3.14 ± 0.43 cm and 0.38 ± 0.07 mmol/L for BFT and NEFA, respectively). Total liver lipid percentage tended to be lower in normal- versus overconditioned cows (4.63 ± 0.40% and 6.06 ± 0.44%, respectively). In comparison to the mean liver lipid percentage of the normal- and overconditioned cows, 1 overconditioned cow had a relatively low (5.21%) and 1 normal-conditioned cow had a relatively high (6.07%) liver lipid percentage. Differentially expressed genes analysis (edgeR quasi-likelihood method) showed that normal-conditioned cows presented 11 upregulated and 12 downregulated genes in comparison to overconditioned cows. Linear discriminant analysis effects size revealed 133 differentially expressed genes between normal- versus overconditioned cows. Notably, the liver of normal-conditioned cows had upregulated genes associated with liver functionality (ALB, SELENOP, IGF1, and IGF2). On the other hand, overconditioned cows had upregulated genes associated with the acute-phase response (C3, HPX, and, LBP). High basal lipolysis in overconditioned cows at the end of pregnancy increased liver lipid content, and this may alter the hepatic gene expression pattern to a pro-inflammatory state.

Greater liver PNPLA3 protein abundance in vivo and in vitro supports lower triglyceride accumulation in dairy cows

Fatty liver syndrome is a prevalent metabolic disorder in peripartum dairy cows that unfavorably impacts lactation performance and health. Patatin-like phospholipase domain-containing protein 3 (PNPLA3) is a lipase that plays a central role in human non-alcoholic fatty liver disease etiology but has received limited attention in bovine fatty liver research. Thus, we investigated the relationship between tissue PNPLA3 expression and liver triglyceride accumulation in vivo via a ketosis induction protocol in multiparous dairy cows peripartum, as well as in vitro via small interfering RNA knockdown of PNPLA3 mRNA expression in bovine primary hepatocytes. Results demonstrated a negative association (P = 0.04) between liver PNPLA3 protein abundance and liver triglyceride content in peripartum dairy cows, while adipose PNPLA3 protein abundance was not associated with liver triglyceride content or blood fatty acid concentration. Knockdown of PNPLA3 mRNA resulted in reduced PNPLA3 protein abundance (P < 0.01) and greater liver triglyceride content (P < 0.01). Together, these results suggest greater liver PNPLA3 protein abundance may directly limit liver triglyceride accumulation peripartum, potentially preventing bovine fatty liver or accelerating recovery from fatty liver syndrome.

Metabolomic and Proteomic Profiles Associated With Ketosis in Dairy Cows

Ketosis is a common metabolic disease in dairy cows during early lactation. However, information about the metabolomic and proteomic profiles associated with the incidence and progression of ketosis is still limited. In this study, an integrated metabolomics and proteomics approach was performed on blood serum sampled from cows diagnosed with clinical ketosis (case, ≥ 2.60 mmol/L plasma β-hydroxybutyrate; BHBA) and healthy controls (control, < 1.0 mmol/L BHBA). Samples were taken 2 weeks before parturition and 2 weeks after parturition from 19 animals (nine cases, 10 controls). All serum samples (n = 38) were subjected to Liquid Chromatography-Mass Spectrometry (LC-MS) based metabolomic analysis, and 20 samples underwent Data-Independent Acquisition (DIA) LC-MS based proteomic analysis. A total of 97 metabolites and 540 proteins were successfully identified, and multivariate analysis revealed significant differences in both metabolomic and proteomic profiles between cases and controls. We investigated clinical ketosis-associated metabolomic and proteomic changes using statistical analyses. Correlation analysis of statistically significant metabolites and proteins showed 78 strong correlations (correlation coefficient, R ≥ 0.7) between 38 metabolites and 25 proteins, which were then mapped to pathways using IMPaLA. Results showed that ketosis altered a wide range of metabolic pathways, such as metabolism, metabolism of proteins, gene expression and post-translational protein modification, vitamin metabolism, signaling, and disease related pathways. Findings presented here are relevant for identifying molecular targets for ketosis and biomarkers for ketosis detection during the transition period.

Effect of functional single nucleotide polymorphism g.-572 A > G of apolipoprotein A1 gene on resistance to ketosis in Chinese Holstein cows

The ketosis has negative effects on the high-yielding dairy cows during early lactation. Apolipoprotein A1 (APOA1) is a component of high-density lipoprotein. However, the association of APOA1 gene with ketosis, and the molecular mechanisms of expression of APOA1 gene are not fully understood in dairy cows. In this study, expression of APOA1 in the liver and blood was investigated by RT-qPCR and immunohistochemistry, and genetic variation in the 5'-flanking region of the AOPA1 gene was also screened and identified. In addition, correlation of the single nucleotide polymorphisms (SNPs) of APOA1 gene with blood ketone characters, and activity of APOA1 promoter were analyzed in dairy cows. The results showed that ApoA1 protein was expressed in the liver, and the mRNA level of APOA1 was significantly higher in the cows with ketosis comparing to the healthy cows. In addition, a novel SNP (g.-572 A > G) in the core promoter of the APOA1 gene was identified between base g.-714 and g.-68 through transient transfection in both HepG2 cell and FFb cell, and luciferase report assay. Moreover, there was lower concentration of blood β-hydroxybutyrate in cows with genotype GG comparing to the cows with genotypes AA and AG. This study reported for the first time that the genetic variant g.-572 A > G in the core promoter region of APOA1 gene was associated with the ketosis in Chinese Holstein cows, and g.-572 A > G may be used as a genetic marker for ketosis prevention.

The Bovine Metabolome

From an animal health perspective, relatively little is known about the typical or healthy ranges of concentrations for many metabolites in bovine biofluids and tissues. Here, we describe the results of a comprehensive, quantitative metabolomic characterization of six bovine biofluids and tissues, including serum, ruminal fluid, liver, Longissimus thoracis (LT) muscle, semimembranosus (SM) muscle, and testis tissues. Using nuclear magnetic resonance (NMR) spectroscopy, liquid chromatography–tandem mass spectrometry (LC–MS/MS), and inductively coupled plasma–mass spectrometry (ICP–MS), we were able to identify and quantify more than 145 metabolites in each of these biofluids/tissues. Combining these results with previous work done by our team on other bovine biofluids, as well as previously published literature values for other bovine tissues and biofluids, we were able to generate quantitative reference concentration data for 2100 unique metabolites across five different bovine biofluids and seven different tissues. These experimental data were combined with computer-aided, genome-scale metabolite inference techniques to add another 48,628 unique metabolites that are biochemically expected to be in bovine tissues or biofluids. Altogether, 51,801 unique metabolites were identified in this study. Detailed information on these 51,801 unique metabolites has been placed in a publicly available database called the Bovine Metabolome Database.