Advertisement
Research| Volume 106, ISSUE 5, P3465-3476, May 2023

Download started.

Ok

Serum metabolomics assessment of etiological processes predisposing ketosis in water buffalo during early lactation

Open AccessPublished:March 17, 2023DOI:https://doi.org/10.3168/jds.2022-22209

      ABSTRACT

      Metabolic disorders as ketosis are manifestations of the animal's inability to manage the increase in energy requirement during early lactation. Generally, buffaloes show a different response to higher metabolic demands than other ruminants with a lower incidence of metabolic problems, although ketosis is one of the major diseases that may decrease the productivity in buffaloes. The aim of this study was to characterize the metabolic profile of Mediterranean buffaloes (MB) associated with 2 different levels of β-hydroxybutyrate (BHB). Sixty-two MB within 50 days in milk (DIM) were enrolled and divided into 2 groups according to serum BHB concentration: healthy group (37 MB; BHB <0.70 mmol/L; body condition score: 5.00; parity: 3.78; and DIM: 30.70) and group at risk of hyperketonemia (25 MB; BHB ≥0.70 mmol/L; body condition score: 4.50; parity: 3.76; and DIM: 33.20). The statistical analysis was conducted by one-way ANOVA and unpaired 2-sample Wilcoxon tests. Fifty-seven metabolites were identified and among them, 12 were significant or tended to be significant. These metabolites were related to different metabolic changes such as mobilization of body resources, ruminal fermentations, urea cycle, thyroid hormone synthesis, inflammation, and oxidative stress status. These findings are suggestive of metabolic changes related to subclinical ketosis status that should be further investigated to better characterize this disease in the MB.

      Key words

      INTRODUCTION

      The transition period is defined as the period between 3 wk before to 3 wk after parturition, and it is critically important to the health and profitability of dairy cows as well as in buffaloes (
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      ). This phase is characterized by major physiological, nutritional, metabolic, and immunological changes (
      • Raphael W.
      • Sordillo L.M.
      Dietary polyunsaturated fatty acids and inflammation: The role of phospholipid biosynthesis.
      ;
      • Lisuzzo A.
      • Laghi L.
      • Fiore F.
      • Harvatine K.
      • Mazzotta E.
      • Faillace V.
      • Spissu N.
      • Zhu C.
      • Moscati L.
      • Fiore E.
      Evaluation of the metabolomic profile through 1H-NMR spectroscopy in ewes affected by postpartum hyperketonemia.
      ). Buffaloes adjust their metabolism to deal with the considerable increase of energy and nutrient requirements needed for milk production which makes them susceptible to negative energy balance (NEB;
      • Purohit G.N.
      • Gaur M.
      • Saraswat C.S.
      • Bihani D.K.
      Metabolic disorders in the parturient buffalo.
      ). The metabolic adaptation to NEB requires interactions with different energy resources, and its failure may occur in various tissues such as the liver, adipose tissue, and others (
      • Herdt T.H.
      Ruminant adaptation to negative energy balance: Influences on the etiology of ketosis and fatty liver.
      ).
      Metabolic disorders in clinical or subclinical form are manifestations of the animal's inability to manage the greater metabolic demands (
      • Sundrum A.
      Metabolic disorders in the transition period indicate that the dairy cows' ability to adapt is overstressed.
      ). The metabolic response to lactation of buffaloes shows a different pattern compared with other ruminants, as demonstrated by the low incidence of metabolic disorders (
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      ). However, as reported for different buffalo species, a NEB is still one of the major concerns that may decrease the productivity in these ruminants and predispose to other pathologies and fertility disorders (
      • Ghanem M.M.
      • El-Deeb W.M.
      Lecithin cholesterol acyltransferase (LCAT) activity as a predictor for ketosis and parturient haemoglobinuria in Egyptian water buffaloes.
      ;
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ;
      • Sundrum A.
      Metabolic disorders in the transition period indicate that the dairy cows' ability to adapt is overstressed.
      ). The energy deficit is characterized by elevated concentrations of the ketone bodies BHB, acetoacetate, and acetone in blood (hyperketonemia), urine, and milk. Similar to cattle, the disorder could have a clinical and subclinical exhibition of diseases or production decreases in dairy buffaloes (
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ) and represents the inadequate metabolic adaptation to contribute in the development of ketosis (
      • Herdt T.H.
      Ruminant adaptation to negative energy balance: Influences on the etiology of ketosis and fatty liver.
      ). The gold-standard test to diagnose an energy deficit is the measurement of serum BHB concentration. Nevertheless, a specific BHB threshold for dairy buffaloes has not been established and dairy cow reference ranges are often used (
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ;
      • Purohit G.N.
      • Gaur M.
      • Saraswat C.S.
      • Bihani D.K.
      Metabolic disorders in the parturient buffalo.
      ). The early diagnosis of metabolic disorders can be performed by different metabolism analyses, and it is essential to properly treat ongoing disorders (
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ;
      • Gianesella M.
      • Fiore E.
      • Arfuso F.
      • Vecchio D.
      • Curone G.
      • Morgante M.
      • Mazzotta E.
      • Badon T.
      • Rossi P.
      • Bedin S.
      • Zumbo A.
      • Piccione G.
      Serum haptoglobin and protein electrophoretic fraction modifications in buffaloes (Bubalus bubalis) around calving and during early lactation.
      ).
      The metabolic processes can be investigated using the metabolomics approach, which reflects the animals health status (
      • Sun L.
      • Guo Y.
      • Fan Y.
      • Nie H.
      • Wang R.
      • Wang F.
      Metabolic profiling of stages of healthy pregnancy in Hu sheep using nuclear magnetic resonance (NMR).
      ). Metabolomics is an important branch of system biology that studies endogenous metabolism stimulated by internal and external factors (
      • Nicholson J.K.
      • Lindon J.C.
      • Holmes E.
      ‘Metabonomics': Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data.
      ). In recent years, the metabolomic approach has been used in other ruminants as cows and sheep to study metabolic alterations associated with ketosis. The differences found between healthy and diseased animals concern different metabolites belonging to diverse metabolic classes with differences depending on the type of technique used (
      • Zhang G.
      • Dervishi E.
      • Dunn S.M.
      • Mandal R.
      • Liu P.
      • Han B.
      • Wishart D.S.
      • Ametaj B.N.
      Metabotyping reveals distinct metabolic alterations in ketotic cows and identifies early predictive serum biomarkers for the risk of disease.
      ;
      • Fiore E.
      • Lisuzzo A.
      • Tessari R.
      • Spissu N.
      • Moscati L.
      • Morgante M.
      • Gianesella M.
      • Badon T.
      • Mazzotta E.
      • Berlanda M.
      • Contiero B.
      • Fiore F.
      Milk fatty acids composition changes according to β-hydroxybutyrate concentrations in ewes during early lactation.
      ;
      • Lisuzzo A.
      • Laghi L.
      • Faillace V.
      • Zhu C.
      • Contiero B.
      • Morgante M.
      • Mazzotta E.
      • Gianesella M.
      • Fiore E.
      Differences in the serum metabolome profile of dairy cows according to the BHB concentration revealed by proton nuclear magnetic resonance spectroscopy (1H-NMR).
      ). The 1H-nuclear magnetic resonance (1H-NMR) spectroscopy is one of the main platforms of metabolomics because the very simple sample preparation and highly reproducible molecule quantification counter-balance a sensitivity lower than the one granted by other platforms such as mass spectrometry (
      • Jones O.A.H.
      • Cheung V.L.
      An introduction to metabolomics and its potential application in veterinary science.
      ).
      The serum metabolic profiling of buffaloes related to ketone bodies has not been investigated, to the best of the authors' knowledge. For all these reasons, the investigation has been performed in this dairy species, hypothesizing that different etiological processes predispose buffalo to ketosis at the start of lactation compared with other ruminant species. Based on the statements, the aim of the current study was to use 1H-NMR to assess the metabolomic profile of Mediterranean buffaloes (MB) at milking resumption to investigate the metabolic changes associated with different levels of energy deficit.

      MATERIALS AND METHODS

      Animals and Farm

      The current investigation received an institutional approval by the Ethical Animal Care and Use Committee of the University of Naples Federico II (n.PG/2017/0099607). All the clinical procedures were performed by the investigators abide by the common good clinical practices (
      • European Medicines Association
      VICH Topic GL49: Studies to Evaluate the Metabolism and Residues Kinetics of Veterinary Drugs in Human Food-Producing Animals: Validation of -Analytical Methods Used in Residue Depletion Studies - Revision at Step 9 – for Implementation.
      ). Moreover, the farmer was previously informed and in agreement with purposes as well as methods used. Finally, the protocol of this study was carried out according to the standards recommended by the Guide for the Care and Use of Laboratory Animals and Directive 2010/63/EU.
      Sixty-two Italian MB were selected from an artificially induced seasonal calving herd (late winter–springtime) consisting of 400 dairy MB and located in Caserta district (Campania, Italy), between January 2019 and April 2019. The sample size was calculated according to
      • Friedman H.
      Simplified determinations of statistical power, magnitude of effect, and research sample sizes.
      : assuming an effect size of 0.40, a correlation analysis with 0.90 power level, and a 2-tailed significant level of 0.05. All the animals were randomly selected within the entire group of fresh buffaloes (<50 DIM) available during the considered period. All MB were characterized by an average: BCS of 4.87 ± 1.03 points (9-point scale), parity of 3.77 ± 2.13, daily milk production of 14.50 ± 3.27 L, and 31.69 ± 11.90 DIM.
      Regarding the farm, it was casually extracted within a group of 10 (n°5 placed in Caserta district and n°5 in Salerno one) regularly requesting consultancy services at the Veterinary Teaching Hospital – Didactic Mobile Clinic Service of the Department of Veterinary Medicine and Animal Production of Napoli (Italy). Moreover, the group of farms met the following eligibility criteria: (1) a similar herd size (∼400 buffaloes, consistent along the year); (2) a feeding system TMR based, given 2 times/d; (3) the absence of a regular monitoring program for metabolic diseases; (4) housing and overall management system respecting the minimum welfare standard for buffalo (
      • De Rosa G.
      • Grasso F.
      • Winckler C.
      • Bilancione A.
      • Pacelli C.
      • Masucci F.
      • Napolitano F.
      Application of the Welfare Quality protocol to dairy buffalo farms: Prevalence and reliability of selected measures.
      ).
      On the farm, milking buffaloes were kept in a separate group up to 60 DIM, so all trial subjects were selected from one group. The animals were housed in a roofed area consisting of solid-grooved concrete floors in the walking and feeding alleys. The lying area was represented by elevated cubicles covered with rubber mattresses. As reported by
      • De Rosa G.
      • Grasso F.
      • Winckler C.
      • Bilancione A.
      • Pacelli C.
      • Masucci F.
      • Napolitano F.
      Application of the Welfare Quality protocol to dairy buffalo farms: Prevalence and reliability of selected measures.
      , a mean space allowance >16 m2/head was guaranteed, as well as a minimum space of 0.75 m/head, and 0.08 m/head at manger and drinking frontage, respectively. Milking MB were milked twice a day by means of a herringbone parlor. Animals enrolled were fed 2 times/d with a TMR including as ingredients: dry-hay, ryegrass silage (plastic-wrapped baled) and corn silage, buffalo cake (Stick-Florido, Fusco, characterized by 23% of CP originating mainly from legumes and cereals, crude fat 5.5%, crude fiber 7.5%, ash 6.9%, sodium 0.32%), and sodium bicarbonate.
      The chemical composition of the diets used during dry and early lactation periods are reported in Table 1. The composition of TMR was determined using a portable analyzer based on near infrared reflectance spectroscopy (AgriNIR Analyzer, Dinamica Generale s.p.a.).
      Table 1Feed chemical composition of TMR used for dry and milking Mediterranean buffaloes (MB)
      Results obtained by portable analyzer based on near-infrared reflectance spectroscopy.
      ItemMilking MBDried-off MB
      DM (%)5770
      Percentage of DM
       Starch22.67.50
       NFC37.018.0
       CP13.07.50
       ADF22.028.0
       NDF39.064.0
       Ash6.508.20
       Ether extract3.802.40
      1 Results obtained by portable analyzer based on near-infrared reflectance spectroscopy.

      Clinical Procedures

      This study used a cross-sectional experimental design. Mediterranean buffalo selected for enrollment were individually confined to a trimming chute and submitted to complete a clinical examination to rule out organic or systemic diseases, as well as foot disorders recognized as causing a reduction of feed intake in MB (
      • Guccione J.
      • Carcasole C.
      • Alsaaod M.
      • D'Andrea L.
      • Di Loria A.
      • De Rosa A.
      • Ciaramella P.
      • Steiner A.
      Assessment of foot health and animal welfare: Clinical findings in 229 dairy Mediterranean Buffaloes (Bubalus bubalis) affected by foot disorders.
      ,
      • Guccione J.
      • Della Valle G.
      • Carcasole C.
      • Kuhnert P.
      • Alsaaod M.
      Detection of Treponema pedis associated with digital dermatitis in Mediterranean buffalo (Bubalus bubalis).
      ). A complete examination of the gastrointestinal systems (including forestomaches, small intestine, and large intestine) was performed according to
      • Jackson P.
      • Cockcroft P.
      The general clinical examination of cattle.
      , to assess the presence of ongoing clinically obvious pathologies, potentially resulting from a NEB, and the udder was examined for clinical mastitis. During the exam, the BCS of the animals was performed using a 9-point scoring system (score 4.5 = ideal BCS), according to
      • De Rosa G.
      • Napolitano F.
      • Grasso F.
      • Pacelli C.
      • Bordi A.
      On the development of a monitoring scheme of buffalo welfare at farm level.
      and slightly modified by the authors, using visual inspection and manual palpation of 4 areas of the body where MB store fat: ribs, spine, hips, and base of the tail as reported also by
      • Guccione J.
      • Carcasole C.
      • Alsaaod M.
      • D'Andrea L.
      • Di Loria A.
      • De Rosa A.
      • Ciaramella P.
      • Steiner A.
      Assessment of foot health and animal welfare: Clinical findings in 229 dairy Mediterranean Buffaloes (Bubalus bubalis) affected by foot disorders.
      . At the end of the examination procedures, blood sampling was performed by jugular venipuncture with a 10-mL monouse syringe (Becton Dickinson Hypodermic Syringe equipped with 21-gauge needle). Some drops of blood (obtained directly by the syringe) were immediately used for a BHB test in field (FreeStyle, Abbott), while the remaining amount of blood was placed into serum clot activator tubes (Vacutainer, Becton and Dickinson) and centrifuged (908 × g for 15 min at room temperature; centrifuge model DMO412, Giorgio-Bormac s.r.l.) to obtain the serum in field.

      Experimental Design and Biochemical Analysis

      Samples were placed in a cool box (4°C) and brought at the same temperature to the reference laboratory of the University of Naples within 1 h of collection for further investigations. The obtained sera were immediately transferred to Eppendorf tubes (1mL of serum/tube) to obtain 2 aliquots for each animal. These aliquots were immediately sent on dry ice to the Department of Animal Medicine, Production, and Health at the University of Padua (Italy) arriving within 24 h. One aliquot of serum was stored at −18°C until biochemical analysis, and the other one was sent at the same temperature using a portable freezer (CoolFreeze CFX65 W professional, Dometic; minimum temperature −22°C) to the University of Bologna where it was stored at −80°C until metabolomic analysis by 1H-NMR.
      Biochemical analysis was performed using an automatic analyzer (BT3500 Biotecnica Instruments s.p.a.). β-hydroxybutyrate concentration was determined by a RANBUT RX Monza test (Randox). Nonesterified fatty acid (NEFA) concentration was determined by a colorimetric method (NEFA RX Monza test; Randox). Glucose concentration was determined by Glucose Monoreagent, LR (Gesan s.r.l.).
      Considering that there is no specific BHB cut-off for buffaloes to identify animals as healthy and hyperketonemic, we selected a subjective cut-off. Specifically, we used the mean BHB value of 0.4 mmol/L plus 3 standard deviations of 0.1 described in previous work of healthy buffaloes around 30 DIM (
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      ) to establish a subjective BHB cut-off of 0.7 mmol/L. According to serum BHB concentration obtained in the laboratory, MB were divided into 2 groups: healthy group (group H) enrolled 37 MB with level of BHB <0.70 mmol/L and group at risk of hyperketonemia (group K) enrolled 25 MB with level BHB ≥0.70 mmol/L.

      Metabolomics Analysis

      An NMR analysis solution was created with 3-(trimethylsilyl)-propionic-2,2,3,3-d4 acid sodium salt 10 mM in non-deuterium oxide, set at pH 7.00 ± 0.02 by means of 1 M phosphate buffer, also containing 10 μL of NaN3 2 mM. We employed 3-(trimethylsilyl)-propionic-2,2,3,3-d4 acid sodium salt as the NMR chemical-shift reference, while NaN3 avoided microbial proliferation as suggested by
      • Zhu C.
      • Laghi L.
      • Zhang Z.
      • He Y.
      • Wu D.
      • Zhang H.
      • Huang Y.
      • Li C.
      • Zou L.
      First steps toward the giant panda metabolome database: Untargeted metabolomics of feces, urine, serum, and saliva by 1H NMR.
      . Serum samples were prepared for 1H-NMR by thawing them. This was done by placing the still-sealed Eppendorf tubes, each containing 1 mL of serum, under a gentle flux of air at room temperature for 10 min. After centrifuging the same Eppendorf tubes for 15 min at 18,630 × g and 4°C, 700 μL of supernatant was transferred to a new tube where 100 μL of NMR analysis solution was added. Finally, each of the so-obtained samples were centrifuged again at the above conditions (15 min at 18,630 × g and 4°C) immediately before analysis.
      The 1H-NMR spectra were recorded at 298 K with an AVANCE III spectrometer (Bruker) operating at a frequency of 600.13 MHz, equipped with the software Topspin 3.5. Following
      • Zhu C.
      • Petracci M.
      • Li C.
      • Fiore E.
      • Laghi L.
      An untargeted metabolomics investigation of Jiulong Yak (Bos grunniens) meat by 1H-NMR.
      , the signals from broad resonances originating from large molecules were suppressed by a CPMG-filter composed by 400 echoes with a τ of 400 μs and a 180° pulse of 24 μs, for a total filter of 330 ms. The oxide of deuterium and protium residual signal was suppressed by means of presaturation. This was done by employing the cpmgpr1d sequence, part of the standard pulse sequence library. Each spectrum was acquired by summing up 256 transients using 32 K data points over a 7,184 Hz spectral window, with an acquisition time of 2.28 s.
      Differences in water content among samples were taken into consideration by probabilistic quotient normalization (
      • Dieterle F.
      • Ross A.
      • Schlotterbeck G.
      • Senn H.
      Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures: Application in 1H NMR metabonomics.
      ), more reliable than the once more common normalization on creatinine. Spectra phase was manually adjusted in Topspin, whereas the subsequent adjustments were performed in R computational language by means of script developed in-house (
      • R Core Team
      R: A language and environment for statistical computing.
      ). After the removal of the residual water signal, 1H-NMR spectra were baseline-corrected by means of peak detection, according to the “rolling ball” principle (
      • Kneen M.A.
      • Annegarn H.J.
      Algorithm for fitting XRF, SEM, and PIXE X-ray spectra backgrounds.
      ), implemented in the baseline R package (
      • Liland K.H.
      • Almøy T.
      • Mevik B.H.
      Optimal choice of baseline correction for multivariate calibration of spectra.
      ). The signals were assigned by comparing their chemical shift and multiplicity with the Chenomx software library (ver. 8.3, Chenomx Inc.). According to the Metabolomics Standard Initiative for metabolites annotation (
      • Salek R.M.
      • Steinbeck C.
      • Viant M.R.
      • Goodacre R.
      • Dunn W.B.
      The role of reporting standards for metabolite annotation and identification in metabolomic studies.
      ), this allowed the confidence in the identification of each metabolite to be of level 1. To apply NMR as a quantitative technique (
      • Zhu C.
      • Li C.
      • Wang Y.
      • Laghi L.
      Characterization of yak common biofluids metabolome by means of proton nuclear magnetic resonance spectroscopy.
      ), the recycle delay was set to 5 s, by considering the relaxation time of the protons under investigation. The molecules of the first serum sample analyzed were quantified by means of an external standard, by taking advantage of the principle of reciprocity (
      • Hoult D.I.
      The principle of reciprocity.
      ).
      Molecules' quantification was performed by means of rectangular integration, considering one of the corresponding signals, free from interferences (
      • Foschi C.
      • Laghi L.
      • D'Antuono A.
      • Gaspari V.
      • Zhu C.
      • Dellarosa N.
      • Salvo M.
      • Marangoni A.
      Urine metabolome in women with Chlamydia trachomatis infection.
      ).

      Statistical Analysis

      Statistical analysis was performed with R software ver. 4.0.3 (
      • R Core Team
      R: A language and environment for statistical computing.
      ). Normal distributions of data were assessed using the Shapiro-Wilk test. The comparison between the 2 groups was conducted by one-way ANOVA for normally distributed data and the unpaired 2-samples Wilcoxon test for not normally distributed data. Data were expressed as least squares means and standard error of the mean. The significance threshold was set at P-value ≤ 0.05. Parameters that presented a P-value between 0.05 and 0.1 were considered as a trend to significance.
      The MetaboAnalyst 5.0 software (https://www.metaboanalyst.ca) was used to assess metabolite fold change expressed as the ratio between group K and group H concentration of each metabolite. A volcano plot was generated using the fold changes and P-values of all identified metabolites. A principal component analysis (PCA), a partial least squares-discriminant analysis (PLS-DA), and an orthogonal partial least squares-discriminant analysis (OPLS-DA) were generated using autoscaling data (mean-centered and divided by the standard deviation of each variable) to highlight the trends of serum metabolome according to BHB value. The variable importance in projection (VIP) scores was applied to PLS-DA and OPLS-DA analyses to identify the metabolites contributing the most to variance between groups. A hierarchical clustering heatmap was then generated to identify trends of significant metabolites between groups. Using the MetaboAnalyst 5.0 software, the website of PubChem (https://pubchem.ncbi.nlm.nih.gov/), Human Metabolome Database (https://hmdb.ca/metabolites/), and Kyoto Encyclopedia of Genes and Genomes (https://www.genome.jp/kegg/) were consulted to assess the function of significant and not significant metabolites (
      • Shi W.
      • Yuan X.
      • Cui K.
      • Li H.
      • Fu P.
      • Rehman S.
      • Shi D.
      • Liu Q.
      • Li Z.
      LC-MS/MS based metabolomics reveal candidate biomarkers and metabolic changes in different buffalo species.
      ). The software's function “Enrichment Analysis” was used to assess the metabolic pathways influenced by the increase of BHB concentration.

      RESULTS

      No MB showed clinical signs of diseases at examination time. Mean animal demographic and biochemical data are provided in Table 2 for healthy animals and animals at risk of hyperketonemia. The only significant differences between groups were BHB (P-value < 0.0001) and aspartate transaminase (AST; P-value = 0.034) values.
      Table 2Characterization of Mediterranean buffaloes (MB) categorized as healthy with low BHB (group H) or at risk of hyperketonemia group with high BHB (group K; BHB ≥0.70 mmol/L)
      Results obtained by automatic analyzer.
      Parameter
      NEFA = nonesterified fatty acids; CHO = cholesterol; TGR = triacylglycerol; AST = aspartate amino transferase; ALT = alanine amino transferase; GGT = γ-glutamyl-transferase.
      Group H (37 MB)Group K (25 MB)SEMP-value
      BHB (mmol/L)0.470.750.02<0.0001
      BCS5.004.500.180.112
      Parity3.783.760.390.966
      DIM30.733.22.180.447
      Milk yield (L/d)14.114.60.600.892
      Glucose (mg/dL)64.762.11.250.154
      NEFA (mmol/L)0.250.240.020.460
      CHO (mg/dL)77.283.05.770.486
      TGR (mg/dL)9.4510.00.490.435
      AST (units/L)142.0164.06.930.034
      ALT (units/L)48.348.92.080.837
      GGT (units/L)20.021.51.070.333
      1 Results obtained by automatic analyzer.
      2 NEFA = nonesterified fatty acids; CHO = cholesterol; TGR = triacylglycerol; AST = aspartate amino transferase; ALT = alanine amino transferase; GGT = γ-glutamyl-transferase.
      A total of 57 molecules were characterized in MB serum samples: 27 AA and derivates, 9 organic acids, 5 alcohols, 4 carbohydrates, 3 amines and derivates, 2 fatty acids, 2 ketone bodies, 1 sulfone, 1 vitamin, 1 imidazole, 1 nucleoside, and 1 guanidine. Six of the quantified metabolites were significantly different between groups: glycerol, taurine, creatinine, acetone, acetate, and 3-hydroxybutyrate. Six of the quantified metabolites tended to be significant: methanol, formate, citrate, Glu, Pro, and Gly. The metabolite concentrations that differed between the groups or tended to differ and the fold changes are listed in Table 3, whereas the concentrations and fold changes of the remaining metabolites that did not differ are shown in Supplemental Table S1 (https://data.mendeley.com/datasets/c5h6wpv856;
      • Lisuzzo A.
      • Fiore E.
      • Laghi L.
      • Harvatine K.
      • Mazzotta E.
      • Alterisio M.C.
      • Ciaramella P.
      • Contiero B.
      • Guccione J.
      Serum metabolomics assessment of etiological processes predisposing ketosis in water buffalo during early lactation - Supplementary Table 1. Mendeley Data, V1.
      ). The volcano plot displayed the association between the base-2 logarithm of fold change (x-axis) and the base-10 negative logarithm of P-value (y-axis; Supplemental Figure S1, https://data.mendeley.com/datasets/crx6hygbf8;
      • Lisuzzo A.
      • Fiore E.
      • Laghi L.
      • Harvatine K.
      • Mazzotta E.
      • Alterisio M.C.
      • Ciaramella P.
      • Contiero B.
      • Guccione J.
      Serum metabolomics assessment of etiological processes predisposing ketosis in water buffalo during early lactation - Supplementary Figure 1. Mendeley Data, V1.
      ). The PCA analysis was first conducted as an overview with an unsupervised method to identify the difference between groups. However, the groups were not well clustered (Figure 1A). Then, PLS-DA and OPLS-DA analyses were performed to maximize the separation between groups (Figure 1B and 1D). According to VIP score >1.5, the most important metabolites for this separation were acetate, BHB, acetone, and glycerol (Figure 1C and 1E). Two heatmaps were generated: the first one was generated for each MB enrolled (Figure 2A), whereas the second one reflects the trends within the group for significant metabolites (Figure 2B).
      Table 3Representative metabolites of Mediterranean buffaloes (MB) expressed as μmol/L of serum categorized as healthy with low BHB (group H) or at risk of hyperketonemia group with high BHB (group K; BHB ≥0.70 mmol/L)
      Results obtained by 1H-nuclear magnetic resonance spectroscopy.
      ClassMetaboliteGroup HGroup KSEMFC
      Fold change expressed as ratio between group K on group H.
      Log2FC
      Base-2 logarithm of fold change.
      P-value
      AA and derivates (5 of 27)Creatinine12.611.00.520.91−0.140.047
      Glycine138.0127.04.130.92−0.120.067
      Glutamate86.591.12.741.050.080.058
      Proline29.127.80.610.96−0.060.068
      Taurine10.79.560.370.90−0.160.040
      Organic acids (3 of 9)Acetate198.0310.014.11.570.65<0.0001
      Citrate27.030.41.261.130.170.063
      Formate7.528.200.251.090.120.064
      Alcohols (2 of 5)Glycerol17.715.60.630.88−0.180.022
      Methanol3.252.860.150.88−0.190.096
      Ketone bodies (2 of 2)3-Hydroxybutyrate58.376.02.231.300.38<0.0001
      Acetone2.763.650.211.320.410.003
      1 Results obtained by 1H-nuclear magnetic resonance spectroscopy.
      2 Fold change expressed as ratio between group K on group H.
      3 Base-2 logarithm of fold change.
      Figure thumbnail gr1
      Figure 1Results obtained by the Chemometrics Analysis performed with MetaboAnalyst 5.0. (A) Scores plot of principal component analysis between principal component (PC) 1 and 2 applied to group H (Healthy) and group K (Risk). (B and D) Scores plot of partial least squares of discriminant analysis (B; PLS-DA) and of orthogonal partial least squares of discriminant analysis (D; OPLS-DA) applied to the groups. The variance displayed in the plot is the explained variance for X. (C and E) Variable importance in projection (VIP) scores of the most important metabolites that differ between groups for the PLS-DA (C) and for OPLS-DA (E). The colored boxes on the right indicate the relative concentrations of the corresponding metabolite in the 2 groups.
      Figure thumbnail gr2
      Figure 2Results obtained by the Cluster Analysis performed with MetaboAnalyst 5.0. (A) Heatmap of all metabolite concentrations on 62 Mediterranean buffaloes (MB). (B) Heatmap of significantly different metabolite concentrations in group H (Healthy) and group K (Risk). The top right shows the colorimetric scale: if the color tends to dark red, the metabolite concentration was increased; if the color tends to dark blue, the metabolite concentration was decreased. On top right is also present the class of group H (Healthy; red) and group K (Risk; green), which is represented in the first line of each heatmap.
      The identified metabolites were used to perform a quantitative enrichment analysis (Figure 3) to understand the metabolic pathways influenced by the increase of BHB concentration. Five metabolic pathways were influenced in group K: glyoxylate and dicarboxylate metabolism, pyruvate metabolism, glycolysis or gluconeogenesis, glycerolipid metabolism, and taurine and hypotaurine metabolism (Table 4).
      Figure thumbnail gr3
      Figure 3Dot plot of metabolic pathways influenced by the increase of BHB concentration (group K). Color gradient and symbol size represent significant metabolite changes in the corresponding pathway. Results obtained by the enrichment analysis performed with MetaboAnalyst 5.0.
      Table 4Metabolic pathways influenced by the increase of BHB concentration (group K; BHB ≥0.70 mmol/L)
      Results obtained by the enrichment analysis performed with MetaboAnalyst 5.0.
      Metabolic pathwayTotal
      Total number of metabolites related to metabolic pathway.
      Hits
      Number of identified metabolites related to metabolic pathway.
      Metabolites
      Name of identified metabolites related to metabolic pathway.
      P-value
      Glyoxylate and dicarboxylate metabolism328Acetate, citrate, formate, glycine, glutamate, glutamine, pyruvate, and serine<0.0001
      Pyruvate metabolism224Pyruvate, lactate, acetate, and fumarate<0.0001
      Glycolysis/gluconeogenesis264Ethanol, pyruvate, lactate, and acetate<0.0001
      Glycerolipid metabolism161Glycerol0.022
      Taurine and hypotaurine metabolism81Taurine0.040
      1 Results obtained by the enrichment analysis performed with MetaboAnalyst 5.0.
      2 Total number of metabolites related to metabolic pathway.
      3 Number of identified metabolites related to metabolic pathway.
      4 Name of identified metabolites related to metabolic pathway.

      DISCUSSION

      A specific BHB cut-point for hyperketonemia in dairy buffaloes has not been established and dairy cow cut-points (BHB ≥1.0–1.4 mmol/L) are often used (
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ;
      • Purohit G.N.
      • Gaur M.
      • Saraswat C.S.
      • Bihani D.K.
      Metabolic disorders in the parturient buffalo.
      ;
      • Lisuzzo A.
      • Laghi L.
      • Faillace V.
      • Zhu C.
      • Contiero B.
      • Morgante M.
      • Mazzotta E.
      • Gianesella M.
      • Fiore E.
      Differences in the serum metabolome profile of dairy cows according to the BHB concentration revealed by proton nuclear magnetic resonance spectroscopy (1H-NMR).
      ). According to the study of
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      , BHB concentration has a tendency to increase from the prepartum to postpartum period in dairy buffaloes. In the current study, buffaloes had a BHB concentration under 0.50 mmol/L after 20 DIM. In this study, the BHB concentrations of both groups were under the threshold value for dairy cows. For this reason, animals enrolled in group H were considered healthy (BHB = 0.47 mmol/L; DIM = 31), whereas animals enrolled in group K were considered at risk of hyperketonemia (BHB = 0.75 mmol/L; DIM = 33).

      Evaluation of Biochemical Analysis

      Ketosis is associated with lower glycemia and milk production (
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ;
      • Fiore E.
      • Arfuso F.
      • Gianesella M.
      • Vecchio D.
      • Morgante M.
      • Mazzotta E.
      • Badon T.
      • Rossi P.
      • Bedin S.
      • Piccione G.
      Metabolic and hormonal adaptation in Bubalus bubalis around calving and early lactation.
      ). The increase of the energy requirement and fat for milk yield around the peak of lactation around 20–30 DIM can lead to a reduction of triacylglycerols which normally range between 14 and 25 mg/dL (
      • Ciaramella P.
      Semeiologia Funzionale indagini collaterali -Allegato 6.
      ;
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      ). We detected no significant difference between groups in triacylglycerols concentration, though both groups had values below the normal reference range (9.45 and 10.0 mg/dL in group H and K, respectively) suggesting a state of energy deficiency at the peak of lactation. The peripartum period is generally associated with higher lipolysis rate due to higher energy demands, which leads to an increase in NEFA concentrations. Nonesterified fatty acids are metabolized in hepatic tissue via complete oxidation to energetic production or partial oxidation to generate ketone bodies. Moreover, NEFA increase is often associate with BHB increases (
      • Fiore E.
      • Giambelluca S.
      • Morgante M.
      • Contiero B.
      • Mazzotta E.
      • Vecchio D.
      • Vazzana I.
      • Rossi P.
      • Arfuso F.
      • Piccione G.
      • Gianesella M.
      Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
      ;
      • Lisuzzo A.
      • Bonelli F.
      • Sgorbini M.
      • Nocera I.
      • Cento G.
      • Mazzotta E.
      • Turini L.
      • Martini M.
      • Salari F.
      • Morgante M.
      • Badon T.
      • Fiore E.
      Differences of the plasma total lipid fraction from pre-foaling to post-foaling period in donkeys.
      ). In this study, NEFA concentration was not different between the groups which may be due to the NEFA increase occurring before study measurements began (
      • McCarthy M.M.
      • Mann S.
      • Nydam D.V.
      • Overton T.R.
      • McArt J.A.A.
      Short communication: Concentrations of nonesterified fatty acids and β-hydroxybutyrate in dairy cows are not well correlated during the transition period.
      ). Liver injury may lead to an increase of hepatic enzymes such as AST. According to the literature, increased AST of 122 units/L has been previously reported in ketotic buffaloes due to fat infiltration in the liver (
      • Russell K.E.
      • Roussel A.J.
      Evaluation of the ruminant serum chemistry profile.
      ;
      • Youssef M.A.
      • El-Khodery S.A.
      • El-deeb W.M.
      • El-Amaiem W.E.E.A.
      Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
      ). In this study, AST concentrations were different between groups with a concentration of 142 and 164 units/L in group H and K, respectively. The higher values of AST may suggest a potential state of hepatic lipidosis in both groups, with a possibly worse state in group K. However, the normal range of cholesterol concentration suggested the absence of fatty liver suggesting that the increment of AST is due to another type of liver or muscle injury (
      • Van Saun R.J.
      Pregnancy toxemia in a flock of sheep.
      ;
      • Russell K.E.
      • Roussel A.J.
      Evaluation of the ruminant serum chemistry profile.
      ).

      Assessment of Metabolome Trends

      The PCA, PLS-DA, and OPLS-DA analyses are multivariate statistical methods that summarize and transform hundreds of metabolite features into few key components (
      • Xia J.
      • Broadhurst D.I.
      • Wilson M.
      • Wishart D.S.
      Translational biomarker discovery in clinical metabolomics: An introductory tutorial.
      ). Among them, the PCA analysis is generally used as an exploratory clustering technique, whereas the PLS-DA is a supervised method which maximize group separation. Furthermore, the OPLS-DA removes the systematic variation that is not correlated with classes improving the interpretation (
      • Sundekilde U.K.
      • Larsen L.B.
      • Bertram H.C.
      NMR-based milk metabolomics.
      ). Our PCA analysis failed to generate a net division between groups. On the contrary, the PLS-DA and OPLS-DA analysis and its VIP plots showed a separation between groups principally due to acetate, 3-hydroxybutyrate, acetone, and glycerol concentrations. The first 3 metabolites showed greater concentrations in group K, while glycerol concentration was greater in group H.

      Metabolomic Analysis: The Lipid Mobilization

      The mobilization of fat stores due to high energy requirement led to an increase of lipolysis in ruminants with a NEB status (
      • de Vries M.J.
      • Veerkamp R.F.
      Energy balance of dairy cattle in relation to milk production variables and fertility.
      ;
      • Li Y.
      • Xu C.
      • Xia C.
      • Zhang H.
      • Sun L.
      • Gao Y.
      Plasma metabolic profiling of dairy cows affected with clinical ketosis using LC/MS technology.
      ). During lipolysis, the triacylglycerols were catabolize in their component, glycerol and NEFA, and released into blood stream (
      • Contreras G.A.
      • Sordillo L.M.
      Lipid mobilization and inflammatory responses during the transition period of dairy cows.
      ). Glycerol represent one of the glucose precursors that animals can use for gluconeogenesis (
      • Drackley J.K.
      • Dann H.M.
      • Douglas N.
      • Guretzky N.A.J.
      • Litherland N.B.
      • Underwood J.P.
      • Loor J.J.
      Physiological and pathological adaptations in dairy cows that may increase susceptibility to periparturient diseases and disorders.
      ;
      • Zhang G.
      • Ametaj B.N.
      Ketosis an old story under a new approach.
      ). The glycerol concentration of our MB showed a reduction in group K. This reduction may be associated with its utilization for gluconeogenesis, possibly indicating a higher energy requirement supplied by glucose precursors. This context may suggest that group K presented a change of metabolism with early ketosis state. Furthermore, ketone bodies are products of fat metabolism that commonly increase during the higher energy requirements and the NEB status (
      • Ceciliani F.
      • Lecchi C.
      • Urh C.
      • Sauerwein H.
      Proteomics and metabolomics characterizing the pathophysiology of adaptive reactions to the metabolic challenges during the transition from late pregnancy to early lactation in dairy cows.
      ;
      • Puppel K.
      • Gołębiewski M.
      • Solarczyk P.
      • Grodkowski G.
      • Slósarz J.
      • Kunowska-Slósarz M.
      • Balcerak M.
      • Przysucha T.
      • Kalińska A.
      • Kuczyńska B.
      The relationship between plasma β-hydroxybutyric acid and conjugated linoleic acid in milk as a biomarker for early diagnosis of ketosis in postpartum Polish Holstein-Friesian cows.
      ). Two ketone bodies (acetone and BHB) were also identified by metabolomic analysis. These metabolites increased in the K group although with a different absolute value than the biochemical analysis probably related to the different sensitivity of analysis.

      Metabolomic Analysis: Muscle Mobilization

      Muscular metabolism is used to supply the energy requirement (
      • Drackley J.K.
      • Dann H.M.
      • Douglas N.
      • Guretzky N.A.J.
      • Litherland N.B.
      • Underwood J.P.
      • Loor J.J.
      Physiological and pathological adaptations in dairy cows that may increase susceptibility to periparturient diseases and disorders.
      ). The mobilization of this tissue is related to a reduction of creatinine (
      • Grasso F.
      • Maria Terzano G.
      • De Rosa G.
      • Tripaldi C.
      • Napolitano F.
      Influence of housing conditions and calving distance on blood metabolites in water buffalo cows.
      ). Creatinine is a spontaneous product of creatine or creatine phosphate breakdown and they are related to total muscular mass and cell energy metabolism (
      • Kohlmeier M.
      Amino acids and nitrogen compounds.
      ;
      • Megahed A.A.
      • Hiew M.W.H.
      • Ragland D.
      • Constable P.D.
      Changes in skeletal muscle thickness and echogenicity and plasma creatinine concentration as indicators of protein and intramuscular fat mobilization in periparturient dairy cows.
      ;
      • Yanibada B.
      • Hohenester U.
      • Pétéra M.
      • Canlet C.
      • Durand S.
      • Jourdan F.
      • Boccard J.
      • Martin C.
      • Eugène M.
      • Morgavi D.P.
      • Boudra H.
      Inhibition of enteric methanogenesis in dairy cows induces changes in plasma metabolome highlighting metabolic shifts and potential markers of emission.
      ). In this study, creatinine concentration was reduced in group K, likely due to increased muscle metabolism, in agreement with previous reports (
      • Nozad S.
      • Ramin A.-G.
      • Moghadam G.
      • Asri-Rezaei S.
      • Babapour A.
      • Ramin S.
      Relationship between blood urea, protein, creatinine, triglycerides, and macro-mineral concentrations with the quality and quantity of milk in dairy Holstein cows.
      ;
      • Zhang G.
      • Ametaj B.N.
      Ketosis an old story under a new approach.
      ).

      Metabolomic Analysis: Ruminal Fermentations

      Acetate is a short-chain fatty acid produced during ruminal fermentation, particularly by fibrous components of the diet (
      • Drackley J.K.
      • Dann H.M.
      • Douglas N.
      • Guretzky N.A.J.
      • Litherland N.B.
      • Underwood J.P.
      • Loor J.J.
      Physiological and pathological adaptations in dairy cows that may increase susceptibility to periparturient diseases and disorders.
      ;
      • Zhu C.
      • Li C.
      • Wang Y.
      • Laghi L.
      Characterization of yak common biofluids metabolome by means of proton nuclear magnetic resonance spectroscopy.
      ). In this study, acetate concentration increased in group K in agreement with previous studies (
      • Zhang G.
      • Ametaj B.N.
      Ketosis an old story under a new approach.
      ;
      • Lisuzzo A.
      • Laghi L.
      • Faillace V.
      • Zhu C.
      • Contiero B.
      • Morgante M.
      • Mazzotta E.
      • Gianesella M.
      • Fiore E.
      Differences in the serum metabolome profile of dairy cows according to the BHB concentration revealed by proton nuclear magnetic resonance spectroscopy (1H-NMR).
      ). Formate is another metabolite produced during ruminal fermentation from methanethiol (
      • Saleem F.
      • Bouatra S.
      • Guo A.C.
      • Psychogios N.
      • Mandal R.
      • Dunn S.M.
      • Ametaj B.N.
      • Wishart D.S.
      The bovine ruminal fluid metabolome.
      ;
      • Yanibada B.
      • Hohenester U.
      • Pétéra M.
      • Canlet C.
      • Durand S.
      • Jourdan F.
      • Boccard J.
      • Martin C.
      • Eugène M.
      • Morgavi D.P.
      • Boudra H.
      Inhibition of enteric methanogenesis in dairy cows induces changes in plasma metabolome highlighting metabolic shifts and potential markers of emission.
      ). The conversion of methanethiol to formate produces hydrogen sulfides and consumes hydrogen. This conversion is performed by anaerobic microbes that compete with methanogens, and for this reason formate level is negatively related to methane (
      • Yanibada B.
      • Hohenester U.
      • Pétéra M.
      • Canlet C.
      • Durand S.
      • Jourdan F.
      • Boccard J.
      • Martin C.
      • Eugène M.
      • Morgavi D.P.
      • Boudra H.
      Inhibition of enteric methanogenesis in dairy cows induces changes in plasma metabolome highlighting metabolic shifts and potential markers of emission.
      ). Methanogens microbes can synthetize methane from methanol (
      • Hook S.E.
      • Wright A.D.G.
      • McBride B.W.
      Methanogens: Methane producers of the rumen and mitigation strategies.
      ). In this study, formate increased in group K while methanol was reduced in the same group. These findings associated with the greater concentration of acetate may suggest an alteration of ruminal microbial population and fermentation according to BHB concentration in this study.

      Metabolomic Analysis: AA Metabolism

      In this study, 4 AA (Gly, taurine, Glu, and Pro) were different between groups. Glycine, taurine, and Pro concentrations were reduced in group K, while Glu increased in the same group. Changes in these AA may suggest links with the urea cycle, tricarboxylic acid (TCA) cycle, oxidative stress, and thyroid hormone synthesis. However, further more specific studies are needed to assess their possible influence.

      Effect on Urea Cycle

      Glutamic acid and Pro are linked to each other, because Pro is a precursor of Glu trough pyrroline-5-carboxylate (
      • Kohlmeier M.
      Amino acids and nitrogen compounds.
      ;
      • Albaugh V.L.
      • Mukherjee K.
      • Barbul A.
      Proline precursors and collagen synthesis: Biochemical challenges of nutrient supplementation and wound healing.
      ). Furthermore, pyrroline-5-carboxylate is a precursor of Orn, an intermediate of the urea cycle (
      • Zhang H.
      • Wu L.
      • Xu C.
      • Xia C.
      • Sun L.
      • Shu S.
      Plasma metabolomic profiling of dairy cows affected with ketosis using gas chromatography/mass spectrometry.
      ;
      • Albaugh V.L.
      • Mukherjee K.
      • Barbul A.
      Proline precursors and collagen synthesis: Biochemical challenges of nutrient supplementation and wound healing.
      ). Orn concentration was slightly reduced in group K. However, it was not significant in our study with a P-value of 0.12. The reduced concentration of Orn and consequent influence of the urea cycle were found in other studies related to ketosis (
      • Zhang H.
      • Wu L.
      • Xu C.
      • Xia C.
      • Sun L.
      • Shu S.
      Plasma metabolomic profiling of dairy cows affected with ketosis using gas chromatography/mass spectrometry.
      ;
      • Guo C.
      • Xue Y.
      • Seddik H.E.
      • Yin Y.
      • Hu F.
      • Mao S.
      Dynamic changes of plasma metabolome in response to severe feed restriction in pregnant ewes.
      ). The reduction of only Pro may suggest an initial influence of the urea cycle before reaching subclinical ketosis state.

      Effect on the Oxidative Stress Status

      The oxidation of Pro to pyrroline-5-carboxylate generates a reactive oxygen species and creates a state of cellular oxidative stress (
      • Krishnan N.
      • Dickman M.B.
      • Becker D.F.
      Proline modulates the intracellular redox environment and protects mammalian cells against oxidative stress.
      ;
      • Kohlmeier M.
      Amino acids and nitrogen compounds.
      ). The presence of free radicals can be controlled by antioxidants such as glutathione (
      • McPherson P.A.C.
      • McEneny J.
      The biochemistry of ketogenesis and its role in weight management, neurological disease, and oxidative stress.
      ). Glutamic acid and Gly are 2 of the 3 AA involved in glutathione synthesis. Moreover, hypotaurine, the reduced form of taurine, is reported to be protective for the cell's oxygen free radicals damages (
      • Kohlmeier M.
      Amino acids and nitrogen compounds.
      ). The reduction concentrations of Pro, Gly, and taurine may suggest an increase of reactive oxygen species production and the use of antioxidant systems represented by glutathione and hypotaurine. This state may indicate an initial oxidative stress status according to the increase of BHB concentration.

      Relationships with Thyroid Hormone Synthesis

      Glutamic acid and taurine can form the glutaurine, a metabolite involved into the regulation of thyroid hormone synthesis (
      • Bittner S.
      • Win T.
      • Gupta R.
      γ-L-glutamyltaurine.
      ;
      • Kohlmeier M.
      Amino acids and nitrogen compounds.
      ). The glutaurine concentrations were reported as inversely linked to triiodothyronine (T3) levels (
      • Bittner S.
      • Win T.
      • Gupta R.
      γ-L-glutamyltaurine.
      ). According to the study of
      • Fiore E.
      • Arfuso F.
      • Gianesella M.
      • Vecchio D.
      • Morgante M.
      • Mazzotta E.
      • Badon T.
      • Rossi P.
      • Bedin S.
      • Piccione G.
      Metabolic and hormonal adaptation in Bubalus bubalis around calving and early lactation.
      on dairy buffaloes, the thyroid hormone levels tended to decrease during the postpartum period. The reduction of taurine may indicate an influence on glutaurine concentrations and a consequent effect on thyroid hormone levels.

      Relationships with the TCA Cycle

      Amino acids are glucogenic substrates that can influence the TCA, with a consequent risk of developing ketosis. Glycine and Glu are precursors of pyruvate and α-ketoglutarate, respectively (
      • Xue Y.
      • Guo C.
      • Hu F.
      • Liu J.
      • Mao S.
      Hepatic metabolic profile reveals the adaptive mechanisms of ewes to severe undernutrition during late gestation.
      ;
      • Lisuzzo A.
      • Laghi L.
      • Faillace V.
      • Zhu C.
      • Contiero B.
      • Morgante M.
      • Mazzotta E.
      • Gianesella M.
      • Fiore E.
      Differences in the serum metabolome profile of dairy cows according to the BHB concentration revealed by proton nuclear magnetic resonance spectroscopy (1H-NMR).
      ). Pyruvate represents the starting point for TCA, by its conversion in acetyl-CoA, and gluconeogenesis. The Gly reduction may also suggest a potential reduction of pyruvate. However, it is the only precursor that presented a difference between the 2 groups. Furthermore, the joining of acetyl-CoA with oxaloacetate generates citrate, an intermediate of TCA (
      • Duffield T.
      Subclinical ketosis in lactating dairy cattle.
      ;
      • Xu W.
      • Van Knegsel A.
      • Saccenti E.
      • Van Hoeij R.
      • Kemp B.
      • Vervoort J.
      Metabolomics of milk reflects a negative energy balance in cows.
      ). In this study, citrate concentrations showed a tendency toward significative increment in group K. This finding associated with the higher concentration of Glu may suggest that the TCA cycle was not influenced in the group at risk of hyperketonemia.

      CONCLUSIONS

      The higher levels of BHB in our group identified as being at risk of hyperketonemia may suggest metabolic changes related to ketosis such as (1) initial mobilization of body resources, (2) a potential state of oxidative stress, (3) possible changes in rumen microbial populations and fermentations, (4) an initial influence on the urea cycle, and (5) thyroid hormone synthesis. These results thus suggest the need to further refine BHB levels in buffaloes by identifying a specific threshold value for this species to clearly define animals as healthy or ketotic.

      ACKNOWLEDGMENTS

      This study was funded and supported by University of Padua (Padua, Italy) in the Bovine OMICS Project (SID Fiore-protocol B.I.RD.-195883/19). Author contributions were as follows: conception of the work: E.F., A.L., J.G., K.J.H., and P.C.; acquisition of data: J.C., M.C.A, and P.C.; analysis of data: E.F., A.L., L.L., V.F., B.C., and C.Z.; interpretation of data: E.F., J.G., M.C.A, and A.L.; drafting the work: E.F., J.G., M.C.A., A.L., V.F., and P.C.; revising the original draft: E.F., A.L., L.L., K.J.H., E.M., M.C.A., C.Z., B.C., V.F., and J.G.; and final approval of original draft: E.F., A.L., L.L., K.J.H., E.M., M.C.A., C.Z., B.C., V.F., J.G., and P.C. The authors have not stated any conflicts of interest.

      REFERENCES

        • Albaugh V.L.
        • Mukherjee K.
        • Barbul A.
        Proline precursors and collagen synthesis: Biochemical challenges of nutrient supplementation and wound healing.
        J. Nutr. 2017; 147 (28978679): 2011-2017
        • Bittner S.
        • Win T.
        • Gupta R.
        γ-L-glutamyltaurine.
        Amino Acids. 2005; 28 (15838590): 343-356
        • Ceciliani F.
        • Lecchi C.
        • Urh C.
        • Sauerwein H.
        Proteomics and metabolomics characterizing the pathophysiology of adaptive reactions to the metabolic challenges during the transition from late pregnancy to early lactation in dairy cows.
        J. Proteomics. 2018; 178 (29055723): 92-106
        • Ciaramella P.
        Semeiologia Funzionale indagini collaterali -Allegato 6.
        in: Semeiologia Clinica Veterinaria. 1st ed. Poletto editore srl, Vermazzo (MI)2014: 618
        • Contreras G.A.
        • Sordillo L.M.
        Lipid mobilization and inflammatory responses during the transition period of dairy cows.
        Comp. Immunol. Microbiol. Infect. Dis. 2011; 34 (21316109): 281-289
        • De Rosa G.
        • Grasso F.
        • Winckler C.
        • Bilancione A.
        • Pacelli C.
        • Masucci F.
        • Napolitano F.
        Application of the Welfare Quality protocol to dairy buffalo farms: Prevalence and reliability of selected measures.
        J. Dairy Sci. 2015; 98 (26233458): 6886-6896
        • De Rosa G.
        • Napolitano F.
        • Grasso F.
        • Pacelli C.
        • Bordi A.
        On the development of a monitoring scheme of buffalo welfare at farm level.
        Ital. J. Anim. Sci. 2005; 4: 115-125
        • de Vries M.J.
        • Veerkamp R.F.
        Energy balance of dairy cattle in relation to milk production variables and fertility.
        J. Dairy Sci. 2000; 83 (10659965): 62-69
        • Dieterle F.
        • Ross A.
        • Schlotterbeck G.
        • Senn H.
        Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures: Application in 1H NMR metabonomics.
        Anal. Chem. 2006; 78 (16808434): 4281-4290
        • Drackley J.K.
        • Dann H.M.
        • Douglas N.
        • Guretzky N.A.J.
        • Litherland N.B.
        • Underwood J.P.
        • Loor J.J.
        Physiological and pathological adaptations in dairy cows that may increase susceptibility to periparturient diseases and disorders.
        Ital. J. Anim. Sci. 2005; 4: 323-344
        • Duffield T.
        Subclinical ketosis in lactating dairy cattle.
        Vet. Clin. North Am. Food Anim. Pract. 2000; 16 (11022338): 231-253
        • European Medicines Association
        VICH Topic GL49: Studies to Evaluate the Metabolism and Residues Kinetics of Veterinary Drugs in Human Food-Producing Animals: Validation of -Analytical Methods Used in Residue Depletion Studies - Revision at Step 9 – for Implementation.
        http://www.vichsec.org
        Date: 2015
        Date accessed: October 28, 2018
        • Fiore E.
        • Arfuso F.
        • Gianesella M.
        • Vecchio D.
        • Morgante M.
        • Mazzotta E.
        • Badon T.
        • Rossi P.
        • Bedin S.
        • Piccione G.
        Metabolic and hormonal adaptation in Bubalus bubalis around calving and early lactation.
        PLoS One. 2018; 13 (29617370)0193803
        • Fiore E.
        • Giambelluca S.
        • Morgante M.
        • Contiero B.
        • Mazzotta E.
        • Vecchio D.
        • Vazzana I.
        • Rossi P.
        • Arfuso F.
        • Piccione G.
        • Gianesella M.
        Changes in some blood parameters, milk composition, and yield of buffaloes (Bubalus bubalis) during the transition period.
        Anim. Sci. J. 2017; 88 (28776872): 2025-2032
        • Fiore E.
        • Lisuzzo A.
        • Tessari R.
        • Spissu N.
        • Moscati L.
        • Morgante M.
        • Gianesella M.
        • Badon T.
        • Mazzotta E.
        • Berlanda M.
        • Contiero B.
        • Fiore F.
        Milk fatty acids composition changes according to β-hydroxybutyrate concentrations in ewes during early lactation.
        Animals (Basel). 2021; 11 (34065915)1371
        • Foschi C.
        • Laghi L.
        • D'Antuono A.
        • Gaspari V.
        • Zhu C.
        • Dellarosa N.
        • Salvo M.
        • Marangoni A.
        Urine metabolome in women with Chlamydia trachomatis infection.
        PLoS One. 2018; 13 (29566085)0194827
        • Friedman H.
        Simplified determinations of statistical power, magnitude of effect, and research sample sizes.
        Educ. Psychol. Meas. 1982; 42: 521-526
        • Ghanem M.M.
        • El-Deeb W.M.
        Lecithin cholesterol acyltransferase (LCAT) activity as a predictor for ketosis and parturient haemoglobinuria in Egyptian water buffaloes.
        Res. Vet. Sci. 2010; 88 (19682715): 20-25
        • Gianesella M.
        • Fiore E.
        • Arfuso F.
        • Vecchio D.
        • Curone G.
        • Morgante M.
        • Mazzotta E.
        • Badon T.
        • Rossi P.
        • Bedin S.
        • Zumbo A.
        • Piccione G.
        Serum haptoglobin and protein electrophoretic fraction modifications in buffaloes (Bubalus bubalis) around calving and during early lactation.
        J. Dairy Res. 2019; 86 (31292012): 291-295
        • Grasso F.
        • Maria Terzano G.
        • De Rosa G.
        • Tripaldi C.
        • Napolitano F.
        Influence of housing conditions and calving distance on blood metabolites in water buffalo cows.
        Ital. J. Anim. Sci. 2004; 3: 275-282
        • Guccione J.
        • Carcasole C.
        • Alsaaod M.
        • D'Andrea L.
        • Di Loria A.
        • De Rosa A.
        • Ciaramella P.
        • Steiner A.
        Assessment of foot health and animal welfare: Clinical findings in 229 dairy Mediterranean Buffaloes (Bubalus bubalis) affected by foot disorders.
        BMC Vet. Res. 2016; 12 (27297174): 107
        • Guccione J.
        • Della Valle G.
        • Carcasole C.
        • Kuhnert P.
        • Alsaaod M.
        Detection of Treponema pedis associated with digital dermatitis in Mediterranean buffalo (Bubalus bubalis).
        Vet. Rec. Case Rep. 2018; 6: 1-6
        • Guo C.
        • Xue Y.
        • Seddik H.E.
        • Yin Y.
        • Hu F.
        • Mao S.
        Dynamic changes of plasma metabolome in response to severe feed restriction in pregnant ewes.
        Metabolites. 2019; 9 (31185597): 112
        • Herdt T.H.
        Ruminant adaptation to negative energy balance: Influences on the etiology of ketosis and fatty liver.
        Vet. Clin. North Am. Food Anim. Pract. 2000; 16 (11022337): 215-230
        • Hook S.E.
        • Wright A.D.G.
        • McBride B.W.
        Methanogens: Methane producers of the rumen and mitigation strategies.
        Archaea. 2010; 2010 (21253540)945785
        • Hoult D.I.
        The principle of reciprocity.
        J. Magn. Reson. 2011; 213 (21889377): 344-346
        • Jackson P.
        • Cockcroft P.
        The general clinical examination of cattle.
        in: Clinical Examination of Farm Animals. 1st ed. John Wiley & Sons, 2002: 9-11
        • Jones O.A.H.
        • Cheung V.L.
        An introduction to metabolomics and its potential application in veterinary science.
        Comp. Med. 2007; 57 (17974125): 436-442
        • Kneen M.A.
        • Annegarn H.J.
        Algorithm for fitting XRF, SEM, and PIXE X-ray spectra backgrounds.
        Nucl. Instrum. Methods Phys. Res. B. 1996; 109–110: 209-213
        • Kohlmeier M.
        Amino acids and nitrogen compounds.
        in: Kohlmeier M. Nutrient Metabolism. 2nd ed. Academic Press, 2015: 265-467
        • Krishnan N.
        • Dickman M.B.
        • Becker D.F.
        Proline modulates the intracellular redox environment and protects mammalian cells against oxidative stress.
        Free Radic. Biol. Med. 2008; 44 (18036351): 671-681
        • Li Y.
        • Xu C.
        • Xia C.
        • Zhang H.
        • Sun L.
        • Gao Y.
        Plasma metabolic profiling of dairy cows affected with clinical ketosis using LC/MS technology.
        Vet. Q. 2014; 34 (25299384): 152-158
        • Lisuzzo A.
        • Fiore E.
        • Laghi L.
        • Harvatine K.
        • Mazzotta E.
        • Alterisio M.C.
        • Ciaramella P.
        • Contiero B.
        • Guccione J.
        Serum metabolomics assessment of etiological processes predisposing ketosis in water buffalo during early lactation - Supplementary Table 1. Mendeley Data, V1.
        • Lisuzzo A.
        • Fiore E.
        • Laghi L.
        • Harvatine K.
        • Mazzotta E.
        • Alterisio M.C.
        • Ciaramella P.
        • Contiero B.
        • Guccione J.
        Serum metabolomics assessment of etiological processes predisposing ketosis in water buffalo during early lactation - Supplementary Figure 1. Mendeley Data, V1.
        • Liland K.H.
        • Almøy T.
        • Mevik B.H.
        Optimal choice of baseline correction for multivariate calibration of spectra.
        Appl. Spectrosc. 2010; 64 (20828437): 1007-1016
        • Lisuzzo A.
        • Bonelli F.
        • Sgorbini M.
        • Nocera I.
        • Cento G.
        • Mazzotta E.
        • Turini L.
        • Martini M.
        • Salari F.
        • Morgante M.
        • Badon T.
        • Fiore E.
        Differences of the plasma total lipid fraction from pre-foaling to post-foaling period in donkeys.
        Animals (Basel). 2022; 12 (35158628): 304
        • Lisuzzo A.
        • Laghi L.
        • Faillace V.
        • Zhu C.
        • Contiero B.
        • Morgante M.
        • Mazzotta E.
        • Gianesella M.
        • Fiore E.
        Differences in the serum metabolome profile of dairy cows according to the BHB concentration revealed by proton nuclear magnetic resonance spectroscopy (1H-NMR).
        Sci. Rep. 2022; 12 (35169190)2525
        • Lisuzzo A.
        • Laghi L.
        • Fiore F.
        • Harvatine K.
        • Mazzotta E.
        • Faillace V.
        • Spissu N.
        • Zhu C.
        • Moscati L.
        • Fiore E.
        Evaluation of the metabolomic profile through 1H-NMR spectroscopy in ewes affected by postpartum hyperketonemia.
        Sci. Rep. 2022; 12 (36183000)16463
        • McCarthy M.M.
        • Mann S.
        • Nydam D.V.
        • Overton T.R.
        • McArt J.A.A.
        Short communication: Concentrations of nonesterified fatty acids and β-hydroxybutyrate in dairy cows are not well correlated during the transition period.
        J. Dairy Sci. 2015; 98 (26188575): 6284-6290
        • McPherson P.A.C.
        • McEneny J.
        The biochemistry of ketogenesis and its role in weight management, neurological disease, and oxidative stress.
        J. Physiol. Biochem. 2012; 68 (21983804): 141-151
        • Megahed A.A.
        • Hiew M.W.H.
        • Ragland D.
        • Constable P.D.
        Changes in skeletal muscle thickness and echogenicity and plasma creatinine concentration as indicators of protein and intramuscular fat mobilization in periparturient dairy cows.
        J. Dairy Sci. 2019; 102 (30954258): 5550-5565
        • Nicholson J.K.
        • Lindon J.C.
        • Holmes E.
        ‘Metabonomics': Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data.
        Xenobiotica. 1999; 29 (10598751): 1181-1189
        • Nozad S.
        • Ramin A.-G.
        • Moghadam G.
        • Asri-Rezaei S.
        • Babapour A.
        • Ramin S.
        Relationship between blood urea, protein, creatinine, triglycerides, and macro-mineral concentrations with the quality and quantity of milk in dairy Holstein cows.
        Vet. Res. Forum. 2012; 3 (25653747): 55-59
        • Puppel K.
        • Gołębiewski M.
        • Solarczyk P.
        • Grodkowski G.
        • Slósarz J.
        • Kunowska-Slósarz M.
        • Balcerak M.
        • Przysucha T.
        • Kalińska A.
        • Kuczyńska B.
        The relationship between plasma β-hydroxybutyric acid and conjugated linoleic acid in milk as a biomarker for early diagnosis of ketosis in postpartum Polish Holstein-Friesian cows.
        BMC Vet. Res. 2019; 15 (31653264): 367
        • Purohit G.N.
        • Gaur M.
        • Saraswat C.S.
        • Bihani D.K.
        Metabolic disorders in the parturient buffalo.
        in: Bubaline Theiogenology. International Veterinary Information Service, 2013: 229-237
        • R Core Team
        R: A language and environment for statistical computing.
        R Foundation for Statistical Computing, Vienna, Austria2018
        • Raphael W.
        • Sordillo L.M.
        Dietary polyunsaturated fatty acids and inflammation: The role of phospholipid biosynthesis.
        Int. J. Mol. Sci. 2013; 14 (24152446): 21167-21188
        • Russell K.E.
        • Roussel A.J.
        Evaluation of the ruminant serum chemistry profile.
        Vet. Clin. North Am. Food Anim. Pract. 2007; 23 (17920455): 403-426
        • Saleem F.
        • Bouatra S.
        • Guo A.C.
        • Psychogios N.
        • Mandal R.
        • Dunn S.M.
        • Ametaj B.N.
        • Wishart D.S.
        The bovine ruminal fluid metabolome.
        Metabolomics. 2013; 9: 360-378
        • Salek R.M.
        • Steinbeck C.
        • Viant M.R.
        • Goodacre R.
        • Dunn W.B.
        The role of reporting standards for metabolite annotation and identification in metabolomic studies.
        Gigascience. 2013; 2 (24131531): 13
        • Shi W.
        • Yuan X.
        • Cui K.
        • Li H.
        • Fu P.
        • Rehman S.
        • Shi D.
        • Liu Q.
        • Li Z.
        LC-MS/MS based metabolomics reveal candidate biomarkers and metabolic changes in different buffalo species.
        Animals (Basel). 2021; 11 (33672725): 560
        • Sun L.
        • Guo Y.
        • Fan Y.
        • Nie H.
        • Wang R.
        • Wang F.
        Metabolic profiling of stages of healthy pregnancy in Hu sheep using nuclear magnetic resonance (NMR).
        Theriogenology. 2017; 92 (28237326): 121-128
        • Sundekilde U.K.
        • Larsen L.B.
        • Bertram H.C.
        NMR-based milk metabolomics.
        Metabolites. 2013; 3 (24957988): 204-222
        • Sundrum A.
        Metabolic disorders in the transition period indicate that the dairy cows' ability to adapt is overstressed.
        Animals (Basel). 2015; 5 (26479480): 978-1020
        • Van Saun R.J.
        Pregnancy toxemia in a flock of sheep.
        J. Am. Vet. Med. Assoc. 2000; 217 (11128547): 1536-1539
        • Xia J.
        • Broadhurst D.I.
        • Wilson M.
        • Wishart D.S.
        Translational biomarker discovery in clinical metabolomics: An introductory tutorial.
        Metabolomics. 2013; 9 (23543913): 280-299
        • Xu W.
        • Van Knegsel A.
        • Saccenti E.
        • Van Hoeij R.
        • Kemp B.
        • Vervoort J.
        Metabolomics of milk reflects a negative energy balance in cows.
        J. Proteome Res. 2020; 19 (32633519): 2942-2949
        • Xue Y.
        • Guo C.
        • Hu F.
        • Liu J.
        • Mao S.
        Hepatic metabolic profile reveals the adaptive mechanisms of ewes to severe undernutrition during late gestation.
        Metabolites. 2018; 8 (30486444): 85
        • Yanibada B.
        • Hohenester U.
        • Pétéra M.
        • Canlet C.
        • Durand S.
        • Jourdan F.
        • Boccard J.
        • Martin C.
        • Eugène M.
        • Morgavi D.P.
        • Boudra H.
        Inhibition of enteric methanogenesis in dairy cows induces changes in plasma metabolome highlighting metabolic shifts and potential markers of emission.
        Sci. Rep. 2020; 10 (32973203)15591
        • Youssef M.A.
        • El-Khodery S.A.
        • El-deeb W.M.
        • El-Amaiem W.E.E.A.
        Ketosis in buffalo (Bubalus bubalis): Clinical findings and the associated oxidative stress level.
        Trop. Anim. Health Prod. 2010; 42 (20582718): 1771-1777
        • Zhang G.
        • Ametaj B.N.
        Ketosis an old story under a new approach.
        Dairy. 2020; 1: 42-60
        • Zhang G.
        • Dervishi E.
        • Dunn S.M.
        • Mandal R.
        • Liu P.
        • Han B.
        • Wishart D.S.
        • Ametaj B.N.
        Metabotyping reveals distinct metabolic alterations in ketotic cows and identifies early predictive serum biomarkers for the risk of disease.
        Metabolomics. 2017; 13: 43
        • Zhang H.
        • Wu L.
        • Xu C.
        • Xia C.
        • Sun L.
        • Shu S.
        Plasma metabolomic profiling of dairy cows affected with ketosis using gas chromatography/mass spectrometry.
        BMC Vet. Res. 2013; 9 (24070026): 186
        • Zhu C.
        • Laghi L.
        • Zhang Z.
        • He Y.
        • Wu D.
        • Zhang H.
        • Huang Y.
        • Li C.
        • Zou L.
        First steps toward the giant panda metabolome database: Untargeted metabolomics of feces, urine, serum, and saliva by 1H NMR.
        J. Proteome Res. 2020; 19 (31994893): 1052-1059
        • Zhu C.
        • Li C.
        • Wang Y.
        • Laghi L.
        Characterization of yak common biofluids metabolome by means of proton nuclear magnetic resonance spectroscopy.
        Metabolites. 2019; 9 (30832316): 41
        • Zhu C.
        • Petracci M.
        • Li C.
        • Fiore E.
        • Laghi L.
        An untargeted metabolomics investigation of Jiulong Yak (Bos grunniens) meat by 1H-NMR.
        Foods. 2020; 9 (32290528): 481