Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle

Mitochondrial sequence variants: testing imputation accuracy and their association with dairy cattle milk traits

BACKGROUND

Mitochondrial genomes differ from the nuclear genome and in humans it is known that mitochondrial variants contribute to genetic disorders. Prior to genomics, some livestock studies assessed the role of the mitochondrial genome but these were limited and inconclusive. Modern genome sequencing provides an opportunity to re-evaluate the potential impact of mitochondrial variation on livestock traits. This study first evaluated the empirical accuracy of mitochondrial sequence imputation and then used real and imputed mitochondrial sequence genotypes to study the role of mitochondrial variants on milk production traits of dairy cattle.

RESULTS

The empirical accuracy of imputation from Single Nucleotide Polymorphism (SNP) panels to mitochondrial sequence genotypes was assessed in 516 test animals of Holstein, Jersey and Red breeds using Beagle software and a sequence reference of 1883 animals. The overall accuracy estimated as the Pearson's correlation squared (R2) between all imputed and real genotypes across all animals was 0.454. The low accuracy was attributed partly to the majority of variants having low minor allele frequency (MAF < 0.005) but also due to variants in the hypervariable D-loop region showing poor imputation accuracy. Beagle software provides an internal estimate of imputation accuracy (DR2), and 10 percent of the total 1927 imputed positions showed DR2 greater than 0.9 (N = 201). There were 151 sites with empirical R2 > 0.9 (of 954 variants segregating in the test animals) and 138 of these overlapped the sites with DR2 > 0.9. This suggests that the DR2 statistic is a reasonable proxy to select sites that are imputed with higher accuracy for downstream analyses. Accordingly, in the second part of the study mitochondrial sequence variants were imputed from real mitochondrial SNP panel genotypes of 9515 Australian Holstein, Jersey and Red dairy cattle. Then, using only sites with DR2 > 0.900 and real genotypes, we undertook a genome-wide association study (GWAS) for milk, fat and protein yields. The GWAS mitochondrial SNP effects were not significant.

CONCLUSION

The accuracy of imputation of mitochondrial genotypes from the SNP panel to sequence was generally low. The Beagle DR2 statistic enabled selection of sites imputed with higher empirical accuracy. We recommend building larger reference populations with mitochondrial sequence to improve the accuracy of imputing less common variants and ensuring that SNP panels include common variants in the D-loop region.

Effect of breed and dietary composition on the miRNA profile of beef steers divergent for feed efficiency

Identifying and breeding cattle that are more feed efficient is of great benefit to beef production. Additionally, it is crucial that genes contributing to feed efficiency are robust across varying management settings including dietary source as well as being relevant across contrasting breeds of cattle. The aim of this study was to determine miRNAs that are contributing to the expression of residual feed intake (RFI) across two breeds and dietary sources. miRNA profiling was undertaken in Longissimus dorsi tissue of Charolais and Holstein-Friesian steers divergent for RFI phenotype following two contrasting consecutive diets (high-forage and high-concentrate). Ten miRNA were identified as differentially expressed (adj. P < 0.1) across the breed and diet contrasts examined. Of particular interest was the differential expression of miR-2419-5p and miR-2415-3p, both of which were up-regulated in the Low-RFI Charolais steers across each dietary phase. Pathway analysis of target mRNA genes of differentially expressed miRNA revealed enrichment (P < 0.05) for pathways including metabolic related pathways, insulin receptor signalling, adipogenesis as well as pathways related to skeletal muscle growth. These results provide insight into the skeletal muscle miRNAome of beef cattle and their potential molecular regulatory mechanisms relating to feed efficiency in beef cattle.

Genome-Wide Association Study as an Efficacious Approach to Discover Candidate Genes Associated with Body Linear Type Traits in Dairy Cattle

Body shape traits are very important and play a crucial role in the economic development of dairy farming. By improving the accuracy of selection for body size traits, we can enhance economic returns across the dairy industry and on farms, contributing to the future profitability of the dairy sector. Registered body conformation traits are reliable and cost-effective tools for use in national cattle breeding selection programs. These traits are significantly related to the production, longevity, mobility, health, fertility, and environmental adaptation of dairy cows. Therefore, they can be considered indirect indicators of economically important traits in dairy cows. Utilizing efficacious genetic methods, such as genome-wide association studies (GWASs), allows for a deeper understanding of the genetic architecture of complex traits through the identification and application of genetic markers. In the current review, we summarize information on candidate genes and genomic regions associated with body conformation traits in dairy cattle worldwide. The manuscript also reviews the importance of body conformation, the relationship between body conformation traits and other traits, heritability, influencing factors, and the genetics of body conformation traits. The information on candidate genes related to body conformation traits provided in this review may be helpful in selecting potential genetic markers for the genetic improvement of body conformation traits in dairy cattle.

An across breed, diet and tissue analysis reveals the transcription factor NR1H3 as a key mediator of residual feed intake in beef cattle

BACKGROUND

Provision of feed is a major determinant of overall profitability in beef production systems, accounting for up to 75% of the variable costs. Thus, improving cattle feed efficiency, by way of determining the underlying genomic control and subsequently selecting for feed efficient cattle, provides a method through which feed input costs may be reduced. The objective of this study was to undertake gene co-expression network analysis using RNA-Sequence data generated from Longissimus dorsi and liver tissue samples collected from steers of two contrasting breeds (Charolais and Holstein-Friesian) divergent for residual feed intake (RFI), across two consecutive distinct dietary phases (zero-grazed grass and high-concentrate). Categories including differentially expressed genes (DEGs) based on the contrasts of RFI phenotype, breed and dietary source, as well as key transcription factors and proteins secreted in plasma were utilised as nodes of the gene co-expression network.

RESULTS

Of the 2,929 DEGs within the network analysis, 1,604 were reported to have statistically significant correlations (≥ 0.80), resulting in a total of 43,876 significant connections between genes. Pathway analysis of clusters of co-expressed genes revealed enrichment of processes related to lipid metabolism (fatty acid biosynthesis, fatty acid β-oxidation, cholesterol biosynthesis), immune function, (complement cascade, coagulation system, acute phase response signalling), and energy production (oxidative phosphorylation, mitochondrial L-carnitine shuttle pathway) based on genes related to RFI, breed and dietary source contrasts.

CONCLUSIONS

Although similar biological processes were evident across the three factors examined, no one gene node was evident across RFI, breed and diet contrasts in both liver and muscle tissues. However within the liver tissue, the IRX4, NR1H3, HOXA13 and ZNF648 gene nodes, which all encode transcription factors displayed significant connections across the RFI, diet and breed comparisons, indicating a role for these transcription factors towards the RFI phenotype irrespective of diet and breed. Moreover, the NR1H3 gene encodes a protein secreted into plasma from the hepatocytes of the liver, highlighting the potential for this gene to be explored as a robust biomarker for the RFI trait in beef cattle.

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

Combined effect of microbially derived cecal SCFA and host genetics on feed efficiency in broiler chickens

BACKGROUND

Improving feed efficiency is the most important goal for modern animal production. The regulatory mechanisms of controlling feed efficiency traits are extremely complex and include the functions related to host genetics and gut microbiota. Short-chain fatty acids (SCFAs), as significant metabolites of microbiota, could be used to refine the combined effect of host genetics and gut microbiota. However, the association of SCFAs with the gut microbiota and host genetics for regulating feed efficiency is far from understood.

RESULTS

In this study, 464 broilers were housed for RFI measuring and examining the host genome sequence. And 300 broilers were examined for cecal microbial data and SCFA concentration. Genome-wide association studies (GWAS) showed that four out of seven SCFAs had significant associations with genome variants. One locus (chr4: 29414391-29417189), located near or inside the genes MAML3, SETD7, and MGST2, was significantly associated with propionate and had a modest effect on feed efficiency traits and the microbiota. The genetic effect of the top SNP explained 8.43% variance of propionate. Individuals with genotype AA had significantly different propionate concentrations (0.074 vs. 0.131 μg/mg), feed efficiency (FCR: 1.658 vs. 1.685), and relative abundance of 14 taxa compared to those with the GG genotype. Christensenellaceae and Christensenellaceae_R-7_group were associated with feed efficiency, propionate concentration, the top SNP genotypes, and lipid metabolism. Individuals with a higher cecal abundance of these taxa showed better feed efficiency and lower concentrations of caecal SCFAs.

CONCLUSION

Our study provides strong evidence of the pathway that host genome variants affect the cecal SCFA by influencing caecal microbiota and then regulating feed efficiency. The cecal taxa Christensenellaceae and Christensenellaceae_R-7_group were identified as representative taxa contributing to the combined effect of host genetics and SCFAs on chicken feed efficiency. These findings provided strong evidence of the combined effect of host genetics and gut microbial SCFAs in regulating feed efficiency traits. Video Abstract.

Plasma and milk metabolomics in lactating sheep divergent for feed efficiency

Enhancing the ability of animals to convert feed into meat or milk by optimizing feed efficiency (FE) has become a priority in livestock research. Although untargeted metabolomics is increasingly used in this field and may improve our understanding of FE, no information in this regard is available in dairy ewes. This study was conducted to (1) discriminate sheep divergent for FE and (2) provide insights into the physiological mechanisms contributing to FE through high-throughput metabolomics. The ultra-high-performance liquid chromatography-quadrupole time-of-flight mass spectrometry (UHPLC-Q/TOF-MS) technique was applied to easily accessible animal fluids (plasma and milk) to assess whether their metabolome differs between high- and low-feed efficient lactating ewes (H-FE and L-FE groups, respectively; 8 animals/group). Blood and milk samples were collected on the last day of the 3-wk period used for FE estimation. A total of 793 features were detected in plasma and 334 in milk, with 100 and 38 of them, respectively, showing differences between H-FE and L-FE. The partial least-squares discriminant analysis separated both groups of animals regardless of the type of sample. Plasma allowed the detection of a greater number of differential features; however, results also supported the usefulness of milk, more easily accessible, to discriminate dairy sheep divergent for FE. Regarding pathway analysis, nitrogen metabolism (either anabolism or catabolism) seemed to play a central role in FE, with plasma and milk consistently indicating a great impact of AA metabolism. A potential influence of pathways related to energy/lipid metabolism on FE was also observed. The variable importance in the projection plot revealed 15 differential features in each matrix that contributed the most for the separation in H-FE and L-FE, such as l-proline and phosphatidylcholine 20:4e in plasma or l-pipecolic acid and phosphatidylethanolamine (18:2) in milk. Overall, untargeted metabolomics provided valuable information into metabolic pathways that may underlie FE in dairy ewes, with a special relevance of AA metabolism in determining this complex phenotype in the ovine. Further research is warranted to validate these findings.

Natural Antioxidant By-Product Mixture Counteracts the Effects of Aflatoxin B1 and Ochratoxin A Exposure of Piglets after Weaning: A Proteomic Survey on Liver Microsomal Fraction

Mycotoxins are toxic compounds produced by certain strains of fungi that can contaminate raw feed materials. Once ingested, even in small doses, they cause multiple health issues for animals and, downstream, for people consuming meat. It was proposed that inclusion of antioxidant-rich plant-derived feed might diminish the harmful effects of mycotoxins, maintaining the farm animals' health and meat quality for human consumption. This work investigates the large scale proteomic effects on piglets' liver of aflatoxin B1 and ochratoxin A mycotoxins and the potential compensatory effects of grapeseed and sea buckthorn meal administration as dietary byproduct antioxidants against mycotoxins' damage. Forty cross-bred TOPIGS-40 hybrid piglets after weaning were assigned to three (n = 10) experimental groups (A, M, AM) and one control group (C) and fed with experimental diets for 30 days. After 4 weeks, liver samples were collected, and the microsomal fraction was isolated. Unbiased label-free, library-free, data-independent acquisition (DIA) mass spectrometry SWATH methods were able to relatively quantify 1878 proteins from piglets' liver microsomes, confirming previously reported effects on metabolism of xenobiotics by cytochrome P450, TCA cycle, glutathione synthesis and use, and oxidative phosphorylation. Pathways enrichment revealed that fatty acid metabolism, steroid biosynthesis, regulation of actin cytoskeleton, regulation of gene expression by spliceosomes, membrane trafficking, peroxisome, thermogenesis, retinol, pyruvate, and amino acids metabolism pathways are also affected by the mycotoxins. Antioxidants restored expression level of proteins PRDX3, AGL, PYGL, fatty acids biosynthesis, endoplasmic reticulum, peroxisome, amino acid synthesis pathways, and, partially, OXPHOS mitochondrial subunits. However, excess of antioxidants might cause significant changes in CYP2C301, PPP4R4, COL18A1, UBASH3A, and other proteins expression levels. Future analysis of proteomics data corelated to animals growing performance and meat quality studies are necessary.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Full-lactation performance of multiparous dairy cows with differing residual feed intake

Residual feed intake (RFI) is an efficiency trait underpinning profitability and environmental sustainability in dairy production. This study compared performance during a complete lactation of 36 multiparous dairy cows divided into three equal-sized groups with high (HRFI), intermediate (IRFI) or low RFI (LRFI). Residual feed intake was determined by two different equations. Residual feed intake according to the NorFor system was calculated as (RFI_NorFor) = (NE_intake)–(NE_maintenance + NE_gestation + NE_milk—NE_mobilisation + NE_deposition). Residual feed intake according to the USA National Research Council (NRC) (RFI_NRC) was calculated as: RFI = DMI − predicted DMI where predicted_s DMI = [(0.372× ECM)+(0.0968×BW^0.75)]×(1−e^{−0.192×(DIM/7+3.67)}). Cows in the HRFI_NorFor group showed higher daily CH₄ production, CH₄/ECM and CH₄ yield (g/kg DMI) than IRFI_NorFor and LRFI_NorFor cows. Cows characterized by high efficiency (LRFI_NorFor) according to the NorFor system had lower body weight. Dry matter intake and apparent dry matter digestibility were not affected by efficiency group but milk yield was lower in the low efficiency, HRFI_NorFor, group. Cows characterized by high efficiency according to the NRC system (LRFI_NRC) had lower dry matter intake while yield of CH₄ was higher. Daily CH₄ production and CH₄ g/kg ECM did not differ between RFI_NRC groups. Dairy cows characterized by high efficiency (both LRFI_NorFor and LRFI_NRC cows) over a complete lactation mobilized more of their body reserves in early lactation as well as during the complete lactation. The results also indicated great phenotypic variation in RFI between different stages the lactation.

Proteomic analysis of adipose tissue revealing differentially abundant proteins in highly efficient mid-lactating dairy cows

The improvement of nutrient utilization efficiency in dairy cows represents an important task in view of the current rising demand for animal products and sustainable resource usage. In this perspective, the identification of appropriate markers to identify the most efficient animals for dairy production becomes a crucial factor. Residual feed intake (RFI), which represents the difference between predicted and actual intake, is used to define the efficiency of cows. In this study, subcutaneous adipose tissue (AT) was collected from five high efficient (HEF) and five low efficient (LEF) mid-lactation Holstein dairy cows, that represented subgroups of the 20% lowest RFI values (HEF) and highest 20% RFI values (LEF), out of a cohort of 155 cows that were examined for feed efficiency at the individual dairy barn at Volcani Institute, Israel. Adipose samples were examined for proteomic analysis by nano-LC/MS–MS and gene expression by RT-PCR. A total of 101 differential proteins (P ≤ 0.05 and fold change ± 1.5) and two protein networks related to feed efficiency were found between HEF and LEF cows. Among the enriched top canonical pathways, FAT10 signaling, EIF2 signaling, Sirtuin signaling, Acute phase response signaling, Protein ubiquitination and mTOR signaling pathways were related to feed efficiency in AT. Furthermore, abundance of transferrin (TF; FC = 78.35, P = 0.02) enriched pathways, including mTOR signaling, LXR/RXR and FXR/RXR activation was found in AT of HEF cows. Relative mRNA expression of RBM39, which is involved in energy metabolism, was decreased in AT of HEF versus LEF. The relationship found between the AT proteins and/or metabolic pathways and the feed efficiency demonstrates that AT may reflect metabolic adaptations to high efficiency, and suggests that these proteins together with their metabolic mechanisms are suitable candidates as biomarkers to identify efficient cows for dairy production.

Genetic and genomic analysis of oxygen consumption in mice

We estimated genetic parameters for oxygen consumption (OC), OC per metabolic body weight (OCMBW) and body weight at three through 8 weeks of age in divergently selected mice populations, with an animal model considering maternal genetic, common litter environmental and cytoplasmic inheritance effects. Cytoplasmic inheritance was considered based on maternal lineage information. With respect to OC, estimated direct heritability was moderate (0.32) and the estimated proportion of the variance of cytoplasmic inheritance effects to the phenotypic variance was very low (0.01), implying that causal genes for OC could be located on autosomes. To assess this hypothesis, we attempted to identify possible candidate causal genes through selective signature detection with the results of pooled whole-genome resequencing using pooled DNA samples from high and low OC mice. We made a list of possible candidate causal genes for OC, including those relating to electron transport chain and ATP-binding proteins (Ndufa12, Sdhc, Atp10b, etc.), Prr16 encoding Largen protein, Cry1 encoding a key component of the circadian core oscillator and so on. The results, although careful interpretation must be required, could contribute to elucidate the genetic mechanism of OC, an indicator for maintenance energy requirement, and therefore feed efficiency.

Characterization of the liver proteome in dairy cows experiencing negative energy balance at early lactation

Negative energy balance (NEB) is associated with metabolic disorders in early lactation dairy cows. Therefore, our objective was to characterize the liver proteome in cows experiencing either NEB or positive energy balance (PEB). Forty-two multiparous Holstein dairy cows were milked either 2 or 3 times daily for the first 30 days in milk (DIM) to alter EB, and were classified retrospectively as NEB (n = 18) or PEB (n = 22). Liver biopsies were collected from 10 cows (n = 5 from each milking frequency) at 17 ± 3 DIM (NEB, n = 6; PEB, n = 4). The liver proteome was characterized using label-free quantitative shotgun proteomics and Ingenuity Pathway Analysis used to identify key affected canonical pathways. Overall, 2741 proteins were identified, and 68 of those were differentially abundant (P ≤ 0.05 and FC ± 1.5). ENO3 (FC = 10.3, P < 0.01) and FABP5 (FC = -12.5, P = 0.045) were the most dramatically upregulated and downregulated proteins, respectively, in NEB cows. Numerous mitochondrial proteins (NDUFA5, NDUFS3, NDUFA6, COX7A2L, COX6C, and COA5) were differentially abundant. Canonical pathways associated with NEB were LPS/IL-1 mediated inhibition of RXR function, oxidative phosphorylation, and mitochondrial dysfunction. Additionally, cows experiencing NEB had less hepatic IL10 transcript abundance than PEB. Together, NEB was associated with altered hepatic inflammatory status, likely due to oxidative stress from mitochondrial dysfunction. SIGNIFICANCE: Our manuscript describes the associations of negative energy balance with the liver proteome in early lactation dairy cows, when metabolic stress and the incidence of diseases is increased. Specifically, we found associations of negative energy balance with shifts in hepatic protein abundance involved in fatty acid uptake, impaired anti-inflammatory responses, and mitochondrial dysfunction. Moving forward, differentially abundant proteins found in this study may be useful as either biological markers for disease or therapeutic targets to improve metabolic adaptations to lactation in postpartum dairy cattle.