# Unraveling the Adipose Tissue Proteome of Transition Cows through Severe Negative Energy Balance

Cristian Piras, Valeria Maria Morittu, Anna Antonella Spina, Alessio Soggiu, Viviana Greco, Christelle Rame, Eric Briant, Namya Mellouk, Bruno Tilocca, Luigi Bonizzi, et al.

► To cite this version:

Cristian Piras, Valeria Maria Morittu, Anna Antonella Spina, Alessio Soggiu, Viviana Greco, et al. Unraveling the Adipose Tissue Proteome of Transition Cows through Severe Negative Energy Balance. Animals, 2019, 9 (12), 10.3390/ani9121013. [DOI: 10.3390/ani9121013](https://doi.org/10.3390/ani9121013)

HAL Id: hal-02622959

https://hal.inrae.fr/hal-02622959

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Unraveling the Adipose Tissue Proteome of Transition Cows through Severe Negative Energy Balance

Cristian Piras <fnref n="1" />,<fnref n="2" />, Valeria Maria Morittu <fnref n="3" />,*, Anna Antonella Spina <fnref n="3" />,*, Alessio Soggiu <fnref n="2" />, Viviana Greco <fnref n="4" />,<fnref n="5" />, Christelle Ramé <fnref n="6" />, Eric Briant <fnref n="6" />, Namya Mellouk <fnref n="6" />, Bruno Tilocca <fnref n="3" />, Luigi Bonizzi <fnref n="7" />, Paola Roncada <fnref n="3" />,* and Joëlle Dupont <fnref n="6" />,*

Received: 22 October 2019; Accepted: 19 November 2019; Published: 21 November 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Simple Summary:** This work described the analysis of differential protein expression of subcutaneous adipose tissue of cows that went under negative energy balance during peripartum. In particular, to the best of our knowledge, it represented an original proteomics study that was able to discriminate cows in negative energy balance up to one month before calving. We believed that our findings would open new perspectives to improve animal welfare during peripartum. To know in advance, the metabolic status of cows would permit to correct the status with appropriate measures, like diet or management.

**Abstract:** Fat mobilization in high-yielding dairy cows during early lactation occurs to overcome negative energy balance (NEB), caused by insufficient feed intake and the concomitant increased nutritional requirements. For this reason, adipose tissue represents an essential organ for healthy and performant lactation. However, only a few data are known about adipose tissue proteome and its metabolic status during peripartum. The aim of this study was to analyze the differential proteomics profiles of subcutaneous adipose tissue belonging to cows with different NEB scores (low NEB and severe NEB). Both groups were analyzed at three different time points (one month before calving, one and sixteen weeks after calving) that were related to different levels and rates of adipose tissue mobilization. The dataset highlighted the differential expression of the same four key proteins (annexin A2, actin-related protein 10, glyceraldehyde-3-phosphate dehydrogenase, and fatty acid-binding protein) involved in lipid metabolism during all time points and of other 22 proteins typical of the other comparisons among remaining time points. The obtained dataset suggested that the individual variability in adipose tissue metabolism/mobilization/energy availability could be linked to the different outcomes in levels of energy balance and related physical complications among dairy cows during peripartum.

<footnote n="1">¹ Department of Chemistry, University of Reading, Reading RG66AH, UK; c.piras@reading.ac.uk</footnote>

<footnote n="2">² Dipartimento di Medicina Veterinaria, Università degli Studi di Milano, 20133 Milano, Italy; alessio.soggiu@unimi.it</footnote>

<footnote n="3">³ Department of Health Sciences, University Magna Graæcia, 88100 Catanzaro, Italy; morittu@unicz.it (V.M.M.); anto.spina90@gmail.com (A.A.S.); tilocca@unicz.it (B.T.)</footnote>

<footnote n="4">⁴ Istituto di Biochimica e Biochimica Clinica, Università Cattolica del Sacro Cuore, 00168 Roma, Italy; viviana.greco@unicatt.it</footnote>

<footnote n="5">5 Fondazione Policlinico Universitario A. Gemelli IRCCS, 00168 Roma, Italy</footnote>

<footnote n="6">⁶ Department of Animal Physiology and Livestock Systems, French National Institute for Agricultural Research—INRA, F-37380 Nouzilly, France; Christelle.rame@inra.fr (C.R.); eric.briant@inra.fr (E.B.); mellouk.namya@nih.gov (N.M.)</footnote>

<footnote n="7">⁷ Dipartimento di Scienze Biomediche, Chirurgiche ed Odontoiatriche, Università degli Studi di Milano, 20133 Milano, Italy; luigi.bonizzi@unimi.it</footnote>

<footnote>* Correspondence: roncada@unicz.it (P.R.); joelle.dupont@inra.fr (J.D.); Tel.: +39-09613694284 (P.R.); +33-0247427789 (J.D.)</footnote>
**Keywords:** energy requirements; peripartum; fat mobilization; proteomics

# 1. Introduction

After parturition, the energy intake of high-yielding dairy cows is often not sufficient to meet the elevated energy requirements for milk production. Body energy reserves, in particular, adipose tissue, are mobilized to meet this energy deficit, and this state is often referred to as negative energy balance (NEB). With the mobilization of 50 to 60 kg of adipose tissue during early lactation, lipolysis is thus the essential metabolic process to support energy supply for all physiological functions [1]. It is well known that prolonged NEB can reduce overall fitness. Declining reproductive fitness is one of the main reasons leading to culling decisions in cattle [2]. Indeed, various measurements of reproduction, including days to first ovulation and days open, have been negatively associated with severe or prolonged NEB during early lactation [3,4]. Additionally, the incidences of metabolic diseases, such as ketosis, displaced abomasum, and reproductive disorders, such as retained placenta and susceptibility to infections, increase during early lactation in animals with severe NEB [5]. In this view, adipose tissue (AT) is a key regulator of metabolism in dairy cows [6]. Traditionally, AT was considered as a passive organ for the storage and mobilization of triglycerides during periods of excess or deficit of energy. However, for several years, AT and, more precisely, white AT is now considered as a dynamic endocrine tissue able to produce and secrete molecules called adipokines or adipocytokines [7,8]. These consist of polypeptides but also non-protein factors that are metabolically active molecules involved in different physiological functions, including immunity, metabolism, cardiovascular system regulation, and angiogenesis. In humans and rodents, adipokines regulate satiety, glucose and lipid metabolism, immune functions, angiogenesis, and reproductive functions [9]. In humans with the increased prevalence of obesity, these new roles of adipokines have spawned an increased interest in white adipose tissue. However, in dairy cows, little is known about the relation between adipose tissue proteins and the metabolic status in peripartum. Previous studies of adipose tissue of cows among parturition showed a differential abundance of proteins related to the insulin resistance phenomenon [10] and provided novel biomarkers for this condition. Using proteomic studies, Zachut et al. also identified novel biomarkers of heat stress in adipose tissue in late-pregnant cows [11]. However, proteomic studies on adipose tissue of dairy cows in relation to the net energy balance have not yet been performed.

Our hypothesis took into consideration the possibility that profound metabolic changes are triggered in the cow adipose tissue at peripartum to overcome the NEB status. A massive energy amount is required to overcome the energy demand, and the mobilization of adipose tissue may play a key role. Individual variability in these mechanisms might be the cause of the different outcomes in terms of NEB status. In the present study, we provided a first proteomic survey of subcutaneous adipose tissue of dairy cows at peripartum. Proteomic profiles of cows with different scores of NEB (severe or low) at different stages of adipose tissue mobilization were comparatively evaluated. The findings of the current study contributed to the steadily increasing knowledge in the field of animal breeding management and opened new avenues in the research for the identification of novel biomarkers predictive of the animal energy status; thus, might imply the improvement of livestock breeding strategies by optimizing diet composition and management practice.

# 2. Materials and Methods

## 2.1. *Animals and Experimental Design*

All experimental protocols were approved by the ethics committee "Comité d'Ethique en Expérimentation Animale Val de Loire" (CEEA VdL, protocol registered as n° 2012-10-4) and were carried out in accordance with the guidelines of the French Council for Animal Care.
From a herd of 39 animals, we selected in the present study two sub-groups of animals either with low negative energy balance (INEB, *n* = 6) or with severe negative energy balance (sNEB, *n* = 7) according to their energy balance (EB) at 1 week postpartum (EB < -13 Mcal/day for sNEB and EB > -2 Mcal/day for INEB). Selected animals were reared in the same environment, and calving took place in the same period. The initial weight (one month before the calving) was 664 ± 15.5 for sNEB and 638 ± 17.9 kg for INEB animals.

## 2.2. Body Weight, Milk Yield, Feeding, and Energy Balance

After each milking, cows were automatically weighted (software RIC version RW1.7, Hokofarm Group, Marknesse, The Netherland). Only the morning live body weight (LBW) was used for weight analyses because the afternoon body weight was more variable. All cows were milked twice daily. At the entrance of the milking parlor, the cows were identified by an electronic collar, and the milk yield of each cow was automatically recorded (software Manufeed 500 pro, vc5 version 2.011.14, Manus-Delaval, Elancourt, France). As LBW is affected by digestive contents, the estimation of empty body weight (EBW) was corrected for the digestive tract content. A change of 4.5 kg of digestive contents per kg of dry matter intake (DMI) was assumed [12]. Variation of EBW (VEBW) was calculated day after day: EBW of the previous day was taken as a reference weight [12]. Live body weight and VEBW as compared to one month before calving were measured from 4 weeks before calving until 16 weeks after calving in females. Dry matter intake was determined from the intake of fresh matter and the dry matter content of each feed of the ration. On average, there was one feeder for two cows. When a cow arrived in front of the feeder, it was recognized by a unique passive transponder attached to her ear tag. If the cow was allowed, the feeder opened, and the quantity of food eaten by the cow was automatically recorded (software RIC version RW1.7, Hokofarm Group, Marknesse, The Netherland). Dry matter intake was calculated daily from calving to week 16 after calving. Energy balance (expressed in Mcal/day) was calculated from calving to week 16 postpartum and corresponded to the difference between net energy intake and net energy needs for body maintenance and lactation.

## 2.3. Subcutaneous Adipose Tissue Thickness

Adipose tissue mobilization was assessed through subcutaneous fat thickness measurements in the sacral region using an ultrasonographic examination with a linear probe (LA 332 3.5/10.0-MHz transducer; Mylab30vet; Esaote, Hospimedi, Saint-Crépin-Ibouvillers, France). Backfat thickness was measured at one month and two weeks before calving, and at 1, 4, 8, 12, and 16 weeks postpartum. The variation of the backfat thickness relatively to week 4 before calving was calculated.

## 2.4. Biopsies of Adipose Tissue

Adipose tissue biopsies were carried out on the 6 animals with INEB and the 7 animals with sNEB at one month before calving, one week after calving, and 16 weeks after calving. Cows had fasted for 12 h before surgery, and anesthesia was induced by intravenous (IV) injections of 12 to 14 mg of xylazine (Rompun®, Bayer, Leverkusen, Germany) and subcutaneous (SC) injections of 20 mg lidocaine (Lurocaïne®, Vetoquinol, Lure, France). Subcutaneous fat was collected from the dewlap under the neck, immediately frozen in liquid nitrogen, and stored at -80 °C until use.

## 2.5. Plasma Non-Esterified Fatty Acids (NEFA) Assays

Blood samples were collected from the tail vein immediately before food distribution, once per week (from 4 weeks before calving until two months after calving) and twice a week (from 8 weeks until 16 weeks after calving). They were centrifuged at 3000× g for 10 min at 4 °C, and plasma was stored at -20 °C until its use for assays. Plasma NEFA was determined using enzymatic colorimetry assay (Wako Chemicals GmbH, Neuss, Germany). The intra- and inter-assay coefficient of variation of plasma NEFA measures were 6% and 7.8%, respectively.
## 2.6. Proteomics

Proteomic profiles of adipose tissue samples derived from cows with a different score in NEB were analyzed. Adipose tissue from three animals with INEB and three animals with sNEB was investigated in technical replicates over three time points for a total of 18 samples, nine belonging to the cows with the highest energetic balance, and nine belonging to the animals with the lowest energetic balance (Figure 1). The samples were collected from the same animal during the three time points that were, respectively, one month before calving, one week after calving, and 16 weeks after calving. Protein separation was performed using 2D SDS-PAGE, and differentially expressed proteins were identified by MALDI-TOF MS.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental design and timeline of sample collection. The adipose tissues of two experimental groups (INEB and sNEB) were analyzed in three different time points. Among the samples of each group, three were collected and analyzed for differential protein expression.</figcaption></figure>

## 2.7. Protein Samples Preparation for 2DE

The adipose tissue samples were dissolved in lysis buffer containing 7 M urea, 2 M thiourea, 2% (*w*/*v*) CHAPS, 2% Triton, 1× protease inhibitor cocktails, 1 mM Na₃VO₄, and 1 mM NaF. Thereafter, samples were processed by agitation with magnets for 2 h. Afterward, protein samples were centrifuged at 14,000 rpm for 30 min to remove insoluble materials. Thus, protein extracts were precipitated in Protein Bind tubes (Eppendorf) following a procedure adapted from Wesse and Fugge [13].

## 2.8. 2D SDS-PAGE

Protein concentration in the samples was determined using Bradford assay, with BioRad protein assay stain and 2 µg/µL concentration of BSA as a standard. Optical density was measured using a spectrophotometer (Gene Quant 100, GE Healthcare, Chicago, Il, USA) at 595 nm. Proteins were separated using 2D SDS-PAGE. For isoelectric focusing (IEF) step, immobilized pH gradient (IPG) polyacrylamide gel strips (GE Healthcare, 7 cm, pH 3–10 NL) and Protean IEF Cell (Bio-Rad, Hercules, Ca, USA) were utilized. Prior to IEF, 100 µg of protein sample was dissolved in a solution containing 7 M urea, 2 M thiourea, 2% *w/v* CHAPS, 2% Triton, 30 mM DTT, 0.5% *w/v* ampholine (pH 3.5–10.0), and 1% *w/v* bromophenol blue. IPG strips were first actively rehydrated in the presence of the sample at 50 V and 20 °C for 17 h. After the rehydration step, paper wicks soaked in water were placed between cathode, anode, and gel strip for preventing high voltage to causing the burning of the strips.
The voltage was gradually increased according to the following protocol: 100 V (1 h), 500 V (1 h), 1000 V (1 h), 2500 V (1 h), 5000 V until the cumulative voltage reached 50 kVh. A limitation of current up to 50 µA per gel strip was set. Following IEF, each strip was reduced for 15 min in 5 mL of solution containing 6 M urea, 2% *w/v* SDS, 50 mM Tris-HCl buffer, pH 8.8, and 20% *v/v* glycerol with 1% *w/v* DTT added, and then alkylated in 5 mL of the same solution with 2.5% *w/v* of iodoacetamide (IAA). IPG strips were then washed shortly in 1× running buffer (250 mM Tris-HCl, pH 8.8, 1920 mM glycine, 1% *w/v* SDS, and MilliQ water), loaded onto 12% *w/v* polyacrylamide resolving gels along with the protein ladder and fixed with 0.5% *w/v* agarose gel. The second dimension was carried out in a Mini-Protean Tetra system (Bio-Rad, Hercules, CA, USA). In the first step of electrophoresis, until the bromophenol blue front line entered the resolving gel, 8 mA per gel for 15 min was applied. In the second step, 16 mA per gel was applied until the bromophenol blue front line reached the bottom of the gel. Gels were then removed from the plates, washed three times for 5 min in 100 mL of deionized water, and left overnight to stain in 100 mL of preheated coomassie brilliant blue G-250 (Sigma-Aldrich, Milano, Italy).

Gel images were acquired using a flatbed scanner (ImageScanner III, GE Healthcare, Uppsala, Sweden) with a resolution of 600 dpi. Before scanning, gels were washed for 20 s in 70% *v/v* ethanol and then for 2 min in 100 mL of deionized water.

## 2.9. MALDI-TOF MS Analysis

According to Piras et al. [14,15], MALDI-TOF MS analysis was performed by Ultraflex III MALDI-TOF/TOF spectrometer (Bruker Daltonics, Macerata, Italy) in positive reflectron mode. Briefly, for the external calibration, the standard peptide mixture calibration was used (Bruker Daltonics, Macerata, Italy): *m/z*: 1046.5418, 1296.6848, 1347.7354, 1619.8223, 2093.0862, 2465.1983, 3147.4710), and to select monoisotopic peptide, MS spectra were analyzed by FlexAnalysis 3.3 software (Bruker Daltonics, Macerata, Italy).

After an internal calibration (known autolysis peaks of trypsin, *m/z*: 842.509 and 2211.104) and exclusion of contaminant ions (known matrix and human keratin peaks), the peak lists were analyzed by MASCOT v.2.4.1 algorithm (www.matrixscience.com) against SwissProt database released 2019_09 restricted to *Bos taurus* taxonomy (561,176 sequences).

For the peptide mass fingerprinting analysis (PMF), the query for database searching was set with these established parameters: carbamidomethylation of cysteines and oxidation on methionines as fixed modification and variable modification, respectively; one missed cleavage site for trypsin; maximal tolerance at 50 ppm. For protein identification assignment, only Mascot scores higher than 54 were considered significant (*p* < 0.05).

To confirm PMF identifications, the instrument was switched in LIFT mode, and MS/MS spectra were acquired with 4–8 × 10³ laser shots using the instrument calibration file. For the fragmentation, precursor ions were manually selected, and the precursor mass window was automatically set. Using Flex-Analysis 3.3 software, each MS/MS spectra was processed by spectra baseline subtraction, smoothing (Savitsky–Golay), and centroiding. For database search analysis, the following parameters were used: carbamidomethylation of cysteines and oxidation on methionine among fixed and variable modifications, respectively; maximum of one missed cleavage; mass tolerance at 50 ppm for precursor ions and at a maximum at 0.4 Da for fragments. The database-dependent search was performed against the SwissProt database released 2019_09 (561,176 sequences) restricted to other Mammalia (13,206 sequences) and *Bos taurus* taxonomy. The confidence interval for protein identification was set to 95% (*p* < 0.05), and only peptides with an individual ion score above the identity threshold were considered correctly identified.

## 2.10. Statistical Analyses

Statistical analyses for dry matter intake, milk yield, live body weight, a variation of empty body weight, a variation of backfat thickness, energy balance, and plasma NEFA were performed with SAS® software version 9.4. Unless differently specified, data were analyzed using the MIXED procedure
for linear mixed models. A repeated effect of time (the week before and after calving) within animals was tested. The residuals from the observations generated from the mixed models were tested for normal distribution.

The model used was:

$$Y_{ijk} = \mu + NEB_i + Week_j + NEB_i \times Week_j + e_{ijk}$$

where $Y_{ijk}$ is the dependent variable (dry matter intake, milk yield, live body weight, variation of empty body weight, variation of backfat thickness, energy balance, plasma), $\mu$ is the overall mean, NEB$_i$ is the fixed effect of NEB $i$ ($i$ = INEB, hNEB), Week$_j$ is the fixed effect of week $j$ ($j$ = 17 classes, 19 classes), NEB$_i \times Week_j$ is the interaction between NEB$_i$ and Week$_j$, and $e_{ijk}$ is the residual error.

Least square means (LSmeans, ± standard error of the mean, SEM) estimated by the models were adjusted using the Scheffe adjustment for multiple-post ANOVA comparisons and compared.

Differences with corresponding $p$-values, $p < 0.05$, were considered as significant.

Variations in protein expression between INEB cows and sNEB cows were analyzed using the Progenesis SameSpots software (Nonlinear Dynamics, Newcastle upon Tyne, UK), Version 4.6. After evaluating the quality of the images, the module for 2D gel analysis was used to align the images, subtract background, detect, normalize, and match spots. All spots were then manually reviewed and selected for excision.

Statistical analysis was performed using the Progenesis Stats module on the log-normalized volumes for all spots. Stats module performed automatically a one-way ANOVA on each spot to evaluate the $p$-value between different groups, and the $p$-values under 0.05 were considered statistically significant.

## 2.11. Bioinformatics Analysis

The functional classification of the identified proteins was performed using the PANTHER classification system [16]. The Gene Ontology (GO) slim biological processes related to the overrepresented and underrepresented proteins were obtained through the gene list analysis module.

High confidence (score 0.700) protein-protein interactions and relations between proteins were investigated via the STRING database of protein-protein interactions [17] using the accession numbers of all differentially expressed proteins.

## 3. Results

Adipose tissue samples were collected from cows in INEB and cows in sNEB. The classification of samples based on energy balance was available only at 1 and 16 weeks postpartum when the EB level could be measured according to milk yield (Figures 2 and 3). Therefore, samples of one month before calving were collected as well and stored at -80 °C until the proper classification was known. Three samples of each group (three with the highest positive INEB and three with the lowest negative sNEB) at each time point were chosen for the comparative proteomics analysis, as described in Figure 1.

Live body weight, a variation of empty body weight, dry matter intake, milk yield, a variation of backfat thickness, energy balance, and plasma NEFA are depicted in Figure 2, Figure 3, and resumed in Table 1. As shown in Figure 2 and Table 1, live body weight (Figure 2A) and milk yield (Figure 2B) were changing during the time points ($p < 0.0001$), but there were no differences among EB groups from calving to week 16 postpartum. However, significant differences among least-square means of sNEB and INEB groups were found in the samples collected at week 4 before calving and the samples collected at week 16 after calving. They consisted of the variation of empty body weight (for NEB effect, difference hNEB-INEB: -6.47 kg/day, $p = 0.041$, Figure 2C), dry matter intake (for NEB effect, difference hNEB-INEB: -3.36 kg/day, $p = 0.002$, Figure 2D), backfat thickness (for NEB effect, difference hNEB-INEB: -18.51%, $p < 0.001$, Figure 2E), and plasma NEFA (for NEB effect, difference hNEB-INEB: 0.42 mmol/L, $p = 0.042$, Figure 3B).
Two-dimensional electrophoresis experiments were performed on biological triplicates for each experimental group (Figure 1) and allowed the detection and measurement of around 700 proteoforms. Spot selection based on a changing 2DE pattern among sample groups resulted in 26 statistically significant proteins (Figure 4) and the final identification of 22 proteins. A representation of the selected gel spots along with a table summarizing the obtained MS data is provided as Figure 4 and Table 2, respectively. Identified proteins were also sorted in a Venn diagram, representing all the proteins that were differentially expressed and identified between the peripartum periods and between groups of INEB and sNEB (one-way ANOVA, *p*-value < 0.05, Figure 5). All differentially expressed proteins commonly shared among comparisons followed the same trend (Table 3), confirming their effective role in the NEB response.

<figure><img src="image_5.png" /><figcaption>B.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>C.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>D.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>E.</figcaption></figure>

<figure><img src="image_9.png" /></figure>

**Figure 2.** Live body weight (A), milk yield (B), variation of empty body weight (C), dry matter intake (D), and variation of backfat thickness (E) in Holstein cows with INEB (*n* = 6) and sNEB (*n* = 7) at week 4 before calving and week 16 after calving. Results are presented as LSmeans ± SEM. Asterisks (*) are representative of statistically significant (*p* < 0.05) differences among LSmeans of the two groups. LSmeans, least square means; SEM, standard error mean.
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Energy balance (A) and plasma NEFA (non-esterified fatty acids) (mmol/L) (B) in Holstein cows with INEB (n = 6) or sNEB (n = 7) at week 4 before calving and week 16 after calving. Results are presented as LSmeans ± SEM.</figcaption></figure>

<table><caption>**Table 1.** Management and nutritional parameters for the two groups of animals (severe Negative Energy Balance (n = 7), and low Negative Energy Balance (n = 6)). Results are presented as LSmeans. *p*-values of the effects of NEB, week peripartum, interaction between Negative Energy Balance, and week peripartum.</caption><thead><tr><th rowspan="2">Group</th><th colspan="2">INEB</th><th colspan="2">sNEB</th><th colspan="3">p-Values</th></tr><tr><th>Lsmeans</th><th>SEM</th><th>Lsmeans</th><th>SEM</th><th>NEB</th><th>Week</th><th>NEB × Week</th></tr></thead><tbody><tr><td>Live body weight, kg/day</td><td>560.56</td><td>0.66</td><td>540.26</td><td>0.78</td><td>0.649</td><td>&lt;0.0001</td><td>0.822</td></tr><tr><td>Variation of empty body weight, kg/day</td><td>-12.19</td><td>0.10</td><td>-18.66</td><td>0.11</td><td>0.041</td><td>&lt;0.0001</td><td>0.774</td></tr><tr><td>Milk yield, kg/day</td><td>22.92</td><td>0.11</td><td>24.29</td><td>0.07</td><td>0.364</td><td>&lt;0.0001</td><td>0.004</td></tr><tr><td>Dry matter intake, kg/day</td><td>16.10</td><td>0.07</td><td>12.74</td><td>0.12</td><td>0.002</td><td>&lt;0.0001</td><td>0.004</td></tr><tr><td>Energy balance, Mcal/day</td><td>-0.26</td><td>0.08</td><td>-7.92</td><td>0.11</td><td>&lt;0.0001</td><td>&lt;0.0001</td><td>0.189</td></tr><tr><td>Variation of backfat thickness (%/week 4 before calving)</td><td>-26.75</td><td>1.04</td><td>-45.26</td><td>0.85</td><td>&lt;0.001</td><td>&lt;0.0001</td><td>0.494</td></tr><tr><td>Non Esterified Fatty Acids, mmol/L</td><td>0.86</td><td>0.02</td><td>1.28</td><td>0.04</td><td>0.042</td><td>&lt;0.0001</td><td>0.485</td></tr></tbody></table>
<figure><img src="image_11.png" /><figcaption>**Figure 4.** Selected gel spots from a representative 2DE gel.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 5.** Venn diagram representing the distribution of differentially expressed proteins. This Venn diagram represents the differential expression of the data represented in Table 3 of the differential protein expression of the sNEB group vs. INEB group.</figcaption></figure>
<table><caption>**Table 2.** Summary of the MS-data for protein identification. <sup>a</sup> Name of identified proteins; <sup>b</sup> Accession No. according to Swiss-Prot database; <sup>c</sup> Peptide mass fingerprinting (PMF) score calculated by MASCOT 2.4.1 algorithm (www.matrixscience.com) after database search; <sup>d</sup> Sequence coverage; <sup>e</sup> Number of experimental peptides matched versus searched peptides; <sup>f</sup> Aminoacidic sequence of the peptides identified by MS/MS analysis, and <sup>g</sup> Related score calculated by MASCOT 2.4.1 algorithm. The underlined letter indicates oxidation of methionine residue; <sup>h</sup> Monoisotopic masses of the parent ions used for MS/MS analysis.</caption><thead><tr><th rowspan="2">Spot No.</th><th rowspan="2">Protein Description <sup>a</sup></th><th rowspan="2">Accession No.<sup>b</sup></th><th rowspan="2">Theoretical Mr(kDa)/pI</th><th colspan="3">PMF</th><th colspan="3">Tandem Mass Spectrometry</th></tr><tr><th>Mascot Score <sup>c</sup></th><th>% Seq. Coverage <sup>d</sup></th><th>Matched Peptides <sup>e</sup></th><th>Peptide Sequence <sup>f</sup></th><th>Mascot Score <sup>g</sup></th><th>m/z <sup>h</sup></th></tr></thead><tbody><tr><td>103</td><td>SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A containing DEAD/H box 1 OS = Bos taurus Phospholipase D2 OS = Bos taurus GN = PLD2 PE = 2 SV = 1</td><td>E1B7X9</td><td>118/5.5</td><td>74</td><td>14</td><td>Nov-81</td><td>801-EMSQMLLK-808<br>763-SINNMEKNTEMCNVMMQLR-781</td><td>86<br>979.49</td><td>78<br>2302.007</td></tr><tr><td>191</td><td></td><td>Q0V8L6</td><td>106.6/7.8</td><td>72</td><td>12</td><td>14/41</td><td>558-HFIQRWNFTK-567<br>419-ALMLLHPNIKVMR-431<br>333-SLNLDPSK-340<br>77-VGIPKETPAITINR-90<br>154-GLTEKMAVSIVDYR-167</td><td>76<br>82<br>96<br>98<br>102</td><td>1376.73<br>1551.96<br>873.43<br>1508.8<br>1597.82</td></tr><tr><td>340</td><td>3-ketoacyl-CoA thiolase, mitochondrial OS = Bos taurus</td><td>Q3T0R7</td><td>42.6/8.06</td><td>64</td><td>29</td><td>26-Oct</td><td>40-LNINTATEEELMTLPGVTR-58<br>109-TYPAHQNR-116</td><td>86<br>92</td><td>2118.81<br>986.5</td></tr><tr><td>448</td><td>Endonuclease/exonuclease/phosphatase family domain-containing protein 1 OS = Bos taurus</td><td>Q3MHJ7</td><td>63.2/8.6</td><td>71</td><td>24</td><td>Oct-61</td><td>104-MNFIKTYPAHQNR-116<br>240-VVIGMDVAASEFYR-253</td><td>106<br>126</td><td>1619.7<br>1556.8</td></tr><tr><td>485</td><td>WD repeat and FYVE domain-containing protein 1 OS = Bos taurus GN = WDFY1 PE = 2 SV = 1</td><td>Q2KIY3</td><td>47.29/6.9</td><td>76</td><td>35</td><td>24-Sep</td><td>163-LAMQEFMILPVGAENFR-179<br>197-GYSFTTTAER-206<br>360-QEYDESGPSIVHR-372<br>154-GLTEKMAVSIVDYR-167</td><td>104<br>86<br>74<br>102</td><td>1965.98<br>1132.96<br>1516.66<br>1597.82</td></tr><tr><td>492</td><td>Alpha-enolase OS = Bos taurus GN = ENO1 PE = 1 SV = 4</td><td>Q9XSJ4</td><td>47.64/6.37</td><td>132</td><td>35</td><td>13/39</td><td>40-LNINTATEEELMTLPGVTR-58<br>119-DAEELRAAMK-129<br>154-GLTEKMAVSIVDYR-167</td><td>86<br>76<br>78</td><td>2118.81<br>1071.55<br>1597.82</td></tr><tr><td>530</td><td>Actin, cytoplasmic 1 OS = Bos taurus GN = ACTB PE = 1 SV = 1</td><td>P60712</td><td>42.05/5.29</td><td>94</td><td>28</td><td>Oct-47</td><td>48-TKGVDEVTIVNILTNR-63<br>314-SLYYYIQDTKGDYQK-329<br>404-NQPPLMKR-411</td><td>122<br>102<br>82</td><td>1771.86<br>2012.91<br>999.61</td></tr><tr><td>622</td><td>Annexin A1 OS = Bos taurus GN = ANXA1 PE = 2 SV = 2</td><td>P46193</td><td>63.2/8.6</td><td>98</td><td>46</td><td>15/74</td><td>279-SVATLILDSLMQCPIDTR-296<br>220-LGNAQLKAPR-229<br>198-AQSLQLSLGDSPLKVR-213<br>68-LVFGILNGR-76</td><td>70<br>78<br>96<br>92</td><td>2012.91<br>1067.68<br>1712.08<br>988.59</td></tr><tr><td>627</td><td>Annexin A1 OS = Bos taurus GN = ANXA1 PE = 2 SV = 2</td><td>P46193</td><td>63.2/8.6</td><td>66</td><td>30</td><td>13/74</td><td>40-DHINLPGFSGENPLR-58<br>1-MATNFLVHEK-10</td><td>86<br>82</td><td>2118.81<br>1189.691</td></tr><tr><td>637</td><td>Annexin A2 OS = Bos taurus GN = ANXA2 PE = 1 SV = 2</td><td>P04272</td><td>38.8/6.2</td><td>106</td><td>49</td><td>Nov-81</td><td>259-SHQVEEHVQSVDIAAFNKI-277</td><td>80</td><td>2151.12</td></tr><tr><td>648</td><td>Actin-related protein 10 OS = Bos taurus GN = ACTR10 PE = 2 SV = 1</td><td>Q3ZBD2</td><td>46.84/7.06</td><td>86</td><td>36</td><td>26/109</td><td></td><td></td><td></td></tr><tr><td>713</td><td>Coiled-coil domain-containing protein 71 OS = Bos taurus GN = CCDC71 PE = 2 SV = 1</td><td>Q2HJ91</td><td>48.08/5.13</td><td>76</td><td>31</td><td>Nov-87</td><td></td><td></td><td></td></tr><tr><td>747</td><td>Purine nucleoside phosphorylase OS = Bos taurus GN = PNP PE = 1 SV = 3</td><td>P55859</td><td>32.24/5.92</td><td>90</td><td>28</td><td>20-Aug</td><td></td><td></td><td></td></tr><tr><td>765</td><td>Elongation factor 1-delta OS = Bos taurus</td><td>A5D989</td><td>30.97/5.42</td><td>82</td><td>36</td><td>Jul-57</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th rowspan="2">Spot No.</th><th rowspan="2">Protein Description <sup>a</sup></th><th rowspan="2">Accession No.<sup>b</sup></th><th rowspan="2">Theoretical Mr(kDa)/pI</th><th colspan="3">PMF</th><th colspan="3">Tandem Mass Spectrometry</th></tr><tr><th>Mascot Score <sup>c</sup></th><th>% Seq. Coverage <sup>d</sup></th><th>Matched Peptides <sup>e</sup></th><th>Peptide Sequence <sup>f</sup></th><th>Mascot Score <sup>g</sup></th><th>m/z <sup>h</sup></th></tr></thead><tbody><tr><td>766</td><td>Glyceraldehyde-3-phosphate dehydrogenase OS = Bos taurus GN = GAPDH PE = 1 SV = 4</td><td>P10096</td><td>36.07/8.5</td><td>66</td><td>32</td><td>May-57</td><td><sup>131</sup>-MGVNHEKYN-140</td><td>66</td><td>929.27</td></tr><tr><td>780</td><td>E3 ubiquitin-protein ligase RNF114 OS = Bos taurus</td><td>Q4U5R4</td><td>26.65/6.5</td><td>96</td><td>36</td><td>Oct-64</td><td><sup>171</sup>-GLMTTVHAIT ATQKTVDGPS-190</td><td>88</td><td>1615.88</td></tr><tr><td>789</td><td>N(G),N(G)-dimethylarginine dimethylaminohydrolase 1 OS = Bos taurus GN = DDAH1 PE = 1 SV = 3</td><td>P56965</td><td>31.6/5.7</td><td>84</td><td>58</td><td>12-Jul</td><td><sup>141</sup>-YTFPCPYCPEK-151</td><td>76</td><td>1463.61</td></tr><tr><td>848</td><td>Succinate dehydrogenase cytochrome b560 subunit, mitochondrial OS = Bos taurus GN = SDHC PE = 1 SV = 2</td><td>P35720</td><td>18.7/9.8</td><td>76</td><td>44</td><td>14-Jun</td><td><sup>11</sup>-DGGAQLAGPAAEADPLGR-28</td><td>98</td><td>1665.74</td></tr><tr><td>850</td><td>General transcription factor IIF subunit 2 OS = Bos taurus GN = GTF2F2 PE = 2 SV = 1</td><td>Q2T9L9</td><td>28.49/9.24</td><td>83</td><td>29</td><td>30-Aug</td><td><sup>150</sup>-GAEILADTFK-159</td><td>78</td><td>1063.55</td></tr><tr><td>855</td><td>Fatty acid-binding protein, adipocyte OS = Bos taurus GN = FABP4 PE = 1 SV = 4</td><td>P48035</td><td>15.35/7.57</td><td>96</td><td>39</td><td>Dec-60</td><td><sup>231</sup>-GHVLLHRTPEEYPESAK-247</td><td>84</td><td>1963.08</td></tr><tr><td>869</td><td>Apolipoprotein A-I OS = Bos taurus GN = APOA1 PE = 1 SV = 3</td><td>P15497</td><td>30.25/5.71</td><td>118</td><td>39</td><td>14/54</td><td><sup>8</sup>-HVGRHCLR-15</td><td>67</td><td>1034.52</td></tr><tr><td>882</td><td>Heat shock protein beta-1 OS = Bos taurus GN = HSPB1 PE = 2 SV = 1</td><td>Q3T149</td><td>22.43/5.98</td><td>90</td><td>36</td><td>23-Jul</td><td><sup>18</sup>-NLGSNRPLSPHITIYR-36</td><td>97</td><td>1836.99</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>26</sup>-YLSQQWAK-33</td><td>86</td><td>1023.57</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>208</sup>-QPVSYLKDLK-218</td><td>92</td><td>1303.8</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>11</sup>-WRLVESK-17</td><td>76</td><td>917.5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>111</sup>-LEDGKLVVVCVMNNVTCTR-129</td><td>98</td><td>2225.11</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>142</sup>-VAPLGEEFR-150</td><td>102</td><td>1017.53</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>131</sup>-WHEEVEIYR-159</td><td>98</td><td>1260.57</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>29</sup>-LFDQAFGLPR-38</td><td>96</td><td>1163.97</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td><sup>13</sup>-GPSWDPFRDWYPAHSR-28</td><td>92</td><td>1973.91</td></tr></tbody></table>

OS: Organism; GN: gene name; PE: Protein existence; SV: Sequence version.
<table><caption>**Table 3.** Differentially expressed proteins in severe negative energy balance and low negative energy balance.</caption><thead><tr><th>Uniprot Name/Accession Number</th><th>String Name</th><th>Same Spots Coding Number</th><th>Protein Name</th><th>One Month before Calving: Trend/Fold Change/p-Value</th><th>One Week after Calving: Trend/Fold Change/p-Value</th><th>16 Weeks after Calving: Trend/Fold Change/p-Value</th></tr></thead><tbody><tr><td>E1B7X9</td><td>SMARCAD1</td><td>103</td><td>SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A containing DEAD/H box 1 OS = Bos taurus</td><td>↓/1.9/0.0193</td><td></td><td></td></tr><tr><td>Q0V8L6</td><td>PLD2</td><td>191</td><td>Phospholipase D2 OS = Bos taurus GN = PLD2 PE = 2 SV = 1</td><td>↓/2.1/0.004</td><td></td><td>↓/2.3/0.030</td></tr><tr><td>Q3T0R7</td><td>ACAA2</td><td>340</td><td>3-ketoacyl-CoA thiolase, mitochondrial OS = Bos taurus</td><td></td><td></td><td>↓/3.3/0.028</td></tr><tr><td>Q3MHJ7</td><td>EEPD1</td><td>448</td><td>Endonuclease/exonuclease/phosphatase family domain-containing protein 1 OS = Bos taurus</td><td>↑/1.6/0.0086</td><td>↑/1.5/0.0072</td><td></td></tr><tr><td>Q2KIY3</td><td>WDFY1</td><td>485</td><td>WD repeat and FYVE domain-containing protein 1 OS = Bos taurus GN = WDFY1 PE = 2 SV = 1</td><td></td><td>↓/1.6/0.0207</td><td></td></tr><tr><td>Q9XSJ4</td><td>ENO1</td><td>492</td><td>Alpha-enolase OS = Bos taurus GN = ENO1 PE = 1 SV = 4</td><td>↓/1.8/0.0004</td><td>↓/1.4/0.050</td><td></td></tr><tr><td>P60712</td><td>NA</td><td>530</td><td>Actin, cytoplasmic 1 OS = Bos taurus GN = ACTB PE = 1 SV = 1</td><td></td><td>↑/1.7/0.0297</td><td>↑/2/0.0080</td></tr><tr><td>P46193</td><td>ANXA1</td><td>622</td><td>Annexin A1 OS = Bos taurus GN = ANXA1 PE = 2 SV = 2</td><td>↓/2/0.00186</td><td></td><td></td></tr><tr><td>P46193</td><td>ANXA1</td><td>627</td><td>Annexin A1 OS = Bos taurus GN = ANXA1 PE = 2 SV = 2</td><td>↓/2.9/0.0429</td><td>↓/1.8/0.0502</td><td></td></tr><tr><td>P04272</td><td>ANXA2</td><td>637</td><td>Annexin A2 OS = Bos taurus GN = ANXA2 PE = 1 SV = 2</td><td>↓/2.7/0.0118</td><td>↓/1.7/0.0085</td><td>↓/1.8/0.0098</td></tr><tr><td>Q3ZBD2</td><td>ACTR10</td><td>648</td><td>Actin-related protein 10 OS = Bos taurus GN = ACTR10 PE = 2 SV = 1</td><td>↓/2.1/0.0050</td><td>↓/1.7/0.0085</td><td>↓/1.4/0.038</td></tr><tr><th>Uniprot Name/Accession Number</th><th>String Name</th><th>Same Spots Coding Number</th><th>Protein Name</th><th>One Month before Calving: Trend/Fold Change/p-Value</th><th>One Week after Calving: Trend/Fold Change/p-Value</th><th>16 Weeks after Calving: Trend/Fold Change/p-Value</th></tr><tr><td>Q2HJ91</td><td>CCDC71</td><td>713</td><td>Coiled-coil domain-containing protein 71 OS = Bos taurus GN = CCDC71 PE = 2 SV = 1</td><td>↓/1.2/0.0026</td><td></td><td></td></tr><tr><td>P55859</td><td>PNP</td><td>747</td><td>Purine nucleoside phosphorylase OS = Bos taurus GN = PNP PE = 1 SV = 3</td><td>↓/1.4/0.0439</td><td>↓/1.3/0.0165</td><td></td></tr><tr><td>A5D989</td><td>EEF1D</td><td>765</td><td>Elongation factor 1-delta</td><td>↓/1.6/0.00141</td><td></td><td></td></tr><tr><td>P10096</td><td>G3PDH</td><td>766</td><td>Glyceraldehyde-3-phosphate dehydrogenase OS = Bos taurus GN = GAPDH PE = 1 SV = 4</td><td>↓/1.8/0.0201</td><td>↓/1.6/0.0260</td><td>↓/1.2/0.040</td></tr><tr><td>Q4U5R4</td><td>RNF114</td><td>780</td><td>E3 ubiquitin-protein ligase RNF114 OS = Bos taurus</td><td>↓/1.3/0.024</td><td></td><td></td></tr><tr><td>P56965</td><td>DDAH1</td><td>789</td><td>N(G),N(G)-dimethylarginine dimethylaminohydrolase 1 OS = Bos taurus GN = DDAH1 PE = 1 SV = 3</td><td>↓/2.2/8.84 × 10<sup>-5</sup></td><td>↓/1.8/0.0011</td><td></td></tr><tr><td>P35720</td><td>SDHC</td><td>848</td><td>Succinate dehydrogenase cytochrome b560 subunit, mitochondrial OS = Bos taurus GN = SDHC PE = 1 SV = 2</td><td></td><td>↑/1.8/0.0059</td><td></td></tr><tr><td>Q2T9L9</td><td>GTF2F2</td><td>850</td><td>General transcription factor IIF subunit 2 OS = Bos taurus GN = GTF2F2 PE = 2 SV = 1</td><td></td><td></td><td>↑/3/0.0104</td></tr><tr><td>P48035</td><td>NA</td><td>855</td><td>Fatty acid-binding protein, adipocyte OS = Bos taurus GN = FABP4 PE = 1 SV = 4</td><td>↑/2.2/0.0053</td><td>↑/2.7/0.0231</td><td>↑/3.3/0.0127</td></tr><tr><td>P15497</td><td>APOA1</td><td>869</td><td>Apolipoprotein A-I OS = Bos taurus GN = APOA1 PE = 1 SV = 3</td><td></td><td>↑/3/0.0028</td><td>↑/2.6/0.0343</td></tr><tr><td>Q3T149</td><td>HSPB1</td><td>882</td><td>Heat shock protein beta-1 OS = Bos taurus GN = HSPB1 PE = 2 SV = 1</td><td>↓/2.8/0.0122</td><td>↓/2.7/0.005</td><td></td></tr></tbody></table>

Table of differentially expressed proteins. ↑: proteins overrepresented in high negative energy balance. ↓: proteins underrepresented in high negative energy balance. OS: Organism; GN: gene name; PE: Protein existence; SV: Sequence version.
Among these, four key proteins, common to all investigated time points, were shown to be differentially expressed, regardless of the time points, between both INEB and sNEB groups, namely annexin, actin-related protein 10, glyceraldehyde-3-phosphate dehydrogenase, and fatty acid-binding protein. Two different isoforms of annexin A1 were underrepresented one month before calving with a fold change of at least two. The underexpression of one of these isoforms was maintained one week after calving, and it returned to a steady state 16 weeks after calving. Differential expression of annexin A2 was the same among all comparisons with a consistent decreasing fold change between one month before calving and the other two time points.

The comparison between INEB and sNEB one month before calving highlighted 16 proteins differentially expressed; among them, two were overrepresented and 14 underrepresented.

The comparison between INEB and sNEB of the samples collected one week after calving highlighted a total of 14 proteins with a different abundance profile, among which five overrepresented and nine underrepresented. The same comparison revealed nine differentially expressed proteins between INEB and sNEB 16 weeks after calving. Among these, four were overrepresented, and five were underrepresented.

To increase the depth of analysis and to be able to illustrate better the biology of the rearrangement of adipose tissue in the peripartum, we also analyzed all class and all proteins through PANTHER GO slim analysis. Moreover, their possible interaction was analyzed through the STRING protein-protein interaction network (Figure 6).

<figure><img src="image_13.png" /><figcaption>**Figure 6.** String protein-protein interaction analysis. The figure highlights the connections between differentially represented proteins.</figcaption></figure>

The interaction analysis showed the link between eight of the differentially expressed proteins, and seven out of eight were less abundant in this biological system. Moreover, all the proteins connected by STRING were differentially expressed in at least two of the three time points, except for the succinate dehydrogenase cytochrome b560 subunit mitochondrial that was highly abundant just one week after calving.

Panther GO analysis highlighted the most representative biological processes in which the differentially expressed proteins were involved. The number of underrepresented proteins was higher than the overrepresented ones. The GO analysis highlighted that overrepresented proteins in sNEB
were involved in a greater variety of metabolic processes, including primary metabolic processes (Figures 7 and 8).

<figure><img src="image_14.png" /><figcaption>**Figure 7.** Functional classification of overrepresented proteins among all time points between sNEB and INEB with the enrichment of most represented processes up to the primary metabolic process. (a) Biological process, (b) metabolic process, (c) primary metabolic process.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 8.** Functional classification of underrepresented proteins of all time points between sNEB and INEB with the enrichment of most represented processes up to the primary metabolic process. (a) Biological process, (b) metabolic process, (c) primary metabolic process (d) generation of precursor metabolites and energy.</figcaption></figure>
# 4. Discussion

Our data suggested that there was a consistent difference in the metabolic profile among lNEB and sNEB experimental groups. By comparing the differentially expressed proteins between both groups, four proteins were common in all time points and constant in their differential abundance during the peripartum, suggesting annexin A2, actin-related protein 10, glyceraldehyde-3-phosphate dehydrogenase, and fatty acid-binding (FABP) protein, together with rates of fat tissue mobilization, as important players for the different response to negative energy balance among experimental groups. Interestingly, the FABP family is the master switch of the PPAR signaling pathway [18], and its constant overrepresentation in sNEB during all sampled periods indicated a pivotal role of the PPAR signaling pathway in the metabolic changes occurring in the sNEB phase.

## 4.1. *Biological Processes of Highly Abundant Proteins*

Among all differentially expressed proteins, six were found to be overrepresented in the comparisons over all the time points. One of them (fatty acid-binding protein) was highly abundant in sNEB in the time point lapse ranging from one month before calving up to 16 weeks after calving. Apolipoprotein A-I was strongly overrepresented in the period after calving (one week to 16 weeks). It is the major component of high-density lipoprotein(HDL), and its major role is related to the transport of fat molecules from the cells. Its overrepresentation in sNEB cows after calving follows the trend of FABP4, suggesting both proteins as active players in resuming the impaired function of adipose tissue of sNEB cows. According to previously published experimental evidence on a mouse model [19], ApoA-I overexpression seems to be positively linked to energy expenditure, the faster reduction of white fat mass, and improved insulin sensitivity. Around calving, in order to compensate for the lack of energy, fat is mobilized from the adipose tissue in the form of non-esterified fatty acids to be transported to different organs. Accordingly, a previous study of Folnožić and colleagues on Holstein-Friesian dairy cows reported an increased lipid mobilization after calving [20].

The non-esterified fatty acid can be again converted into triacylglycerol (TAG) from the liver and stored in hepatocytes. This is mainly the reason why up to 40% of dairy cows develop fatty liver disease and could explain as well such an increased amount of ApoA-1 in adipose tissue. Based on the evidences described by Karavia and colleagues [21] on a mouse model and on the assumption that ApoA-1 is one of the major protein components of HDL, it is consistent to assume that its overexpression is due to the higher amount of NEFA in blood, exactly like in these animals, as well as to a required faster clearance of fatty acids from the liver. In this view, Turk et al. already reported a markedly affected concentration of triglycerides, total cholesterol, HDL, beta-hydroxybutyrate, free fatty acids, and paraoxonase-1 activity in heifers throughout transition period [22,23].

Fatty acid-binding protein's main function is related to lipid transport and has the capability to bind with high-affinity hydrophobic ligands as saturated and unsaturated long-chain fatty acids and eicosanoids (leukotrienes and prostaglandins) [24]. Its strong overrepresentation during all the time points around parturition period could reflect the overall condition of lipid mobilization and faster lipid metabolism, in the attempt to provide energy that attenuates the sNEB.

Fatty acid-binding protein 4 (FABP4) is mainly expressed in adipocytes and macrophages. It is responsible for the development of insulin resistance and inflammation [25], it can reversibly bind to hydrophobic ligands, such as saturated and unsaturated long-chain fatty acids (FAs), and transport FAs to specific organelles in the cell [26]. Normally, the amount of FABP in cells is proportional to the rates of FA metabolism [27]. Its overexpression in our model is consistent and correlates with NEB, and it is indeed overexpressed in the sNEB group. A higher metabolism of adipose tissue is also in agreement with the theory of Vries and Veerkamp that describes how a decrease in fat deposit percentage during early lactation might serve as an indicator of energy balance [28].
## 4.2. Biological Processes of Less Abundant Proteins

Most of the underrepresented proteins were assigned to a variety of metabolic processes. Expanding the section of the lipid metabolic process, it was possible to observe that the underrepresented proteins involved in this biological process were annexin A1 and annexin A2. Relevant experimental evidence documents the key role of these proteins in the regulation of fat tissue metabolism. One of its major roles is linked to fat accumulation in cows [29]. We found a positive correlation between the annexin A1 levels and the backfat accumulation. Likewise, annexin A2 deficiency has been linked to white adipose tissue hypotrophy due to reduced fatty acid uptake by endothelium and adipocytes [30]. In our experimental design, annexins A1 and A2 were both strongly underrepresented in severe negative energy balance cows. Considering the previously cited experimental evidence, it seems that annexins A1 and A2 are positively involved in fat storage and adipose tissue formation. Their reduced abundance in sNEB cows might be due to the reduced capability to fulfill the energy requirements necessary during pregnancy and lactation and undertake a fast mobilization of adipose reserves that are not supported by the expression of annexins.

Insulin resistance could represent another important feature to take into account. The differential protein expression of adipose tissue in cows has been evaluated in relation to insulin resistance [10]. We evaluated the differential protein expression of adipose tissue from cows grouped according to insulin resistance or sensitivity following the evaluation of protein kinase B phosphorylation after insulin stimulation. The outcome highlighted that annexin A1 was highly abundant in insulin-resistant adipose tissue and, other experimental shreds of evidence, highlighted as well, its positive relation with backfat thicknesses [29]. These two experimental pieces of evidence could explain why sNEB adipose tissue is faster in lipid metabolism and less efficient in energy storage.

The generation of precursor metabolites and energy highlighted in the panther pie chart (Figure 8b) were composed of two main proteins, and both were underrepresented and part of glycolysis metabolism (Figure 8d). The first one was alpha-enolase, whose overexpression was particularly lower one month before calving, and then the difference decreased up to become non-significant 16 weeks after calving. The second one was the glyceraldehyde-3-phosphate dehydrogenase, whose reduced abundance was accentuated in sNEB in all the time points, but less relevant in the period close or after parturition. The reduced abundance of these two proteins seemed to demonstrate the decreasing amount of the glycolytic process. It was in agreement with the decreased expression of annexin 2, which is also involved in the translocation of GLUT 4 to the membrane for glucose intake [31]. All this evidence highlighted a strong decrease in glucose metabolism in the adipose tissue of sNEB.

Actin-related protein 10 was another protein that was found to be less abundant among all the time points in sNEB. Its role is mainly related to intracellular trafficking and microtubule-based movement. It was again in agreement with the previously mentioned theory about an sNEB adipose tissue more oriented to lipid transport in other organs and then to energy storage.

## 5. Conclusions

The results obtained showed several differences in the protein profile involved in biological processes and pathways, as glycolysis and lipid metabolism and transport. The most important differences highlighted that glycolysis was dramatically underrepresented in adipocytes of sNEB cows, and, on the contrary, lipid transport out of adipocytes was highly increased (Figure 3B). This observation could be explained by the impairment in lipid biosynthesis and in the fatty acid mobilization that was consistent with the level of NEB. The hypothesized mechanism behind this phenomenon was resumed in Figure 9, representing the combination of obtained results and previous experimental evidence, according to Kuhla and Metges [32].

The major requirement of energy is due to the necessity of lipids necessary for milk production. The adipose tissue of sNEB cows seems to be more efficient in lipid transport and mobilization and less efficient in the glycolytic pathway. This may produce a faster metabolism of adipose tissue of sNEB cows that are as well faster in depleting the fat deposit [33]. This different efficiency in adipose tissue
metabolism might be relevant for the complications due to the NEB among calving and play a key role
in the timing necessary for recovering.

<figure><img src="image_16.png" /><figcaption>**Figure 9.** Resuming scheme of different NEB levels in the physiology of lactating cows. The represented arrows represent the trend of the de-regulated metabolisms, pathways, and metabolites in sNEB cows in comparison to INEB cows. Up arrow (↑) and down arrow (↓) refer to increased and decreased concentrations/activity, respectively.</figcaption></figure>

**Author Contributions:** Conceptualization, P.R. and J.D.; methodology, C.P., A.S., V.G., C.R., E.B., N.M., V.M.M. and A.A.S.; software, C.P., A.S., V.G. and B.T.; formal analysis, C.P., A.S., V.G. and V.M.M.; investigation, C.P., A.S., V.G., C.R., E.B., N.M., V.M.M. and A.A.S.; resources, J.D., P.R.; data curation, C.P., A.S., V.G., C.R., E.B., N.M., P.R., J.D., B.T., L.B., V.M.M. and A.A.S.; writing—original draft preparation, C.P., P.R. and J.D.; writing—review and editing, L.B., B.T., A.S., V.G., P.R., J.D. and V.M.M.; supervision, C.P., A.S., V.G., P.R., B.T., J.D. and L.B.; project administration, P.R. and J.D.; funding acquisition, P.R. and J.D.

**Funding:** This research was funded by European Union Seventh Framework Programs (FP7:2007-2013) under the grant agreement n°311776. (JD, PR).

**Acknowledgments:** We thank Mickael Dupont, Mickael Delanoue, Ludovic Métivier, and Christophe Mouaze of the Experimental Unit UEPAO for animal management and involvement in the experiment.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. Schröder, U.J.; Staufenbiel, R. Invited review: Methods to determine body fat reserves in the dairy cow with special regard to ultrasonographic measurement of backfat thickness. *J. Dairy Sci.* **2006**, *89*, 1–14. [CrossRef]

2. Gröhn, Y.; Rajala-Schultz, P.; Allore, H.; DeLorenzo, M.; Hertla, J.; Galliganc, D. Optimizing replacement of dairy cows: Modeling the effects of diseases. *Prev. Vet. Med.* **2003**, *61*, 27–43. [CrossRef]

3. Patton, J.; Kenny, D.; McNamara, S.; Mee, J.; O'Mara, F.; Diskin, M.; Murphy, J. Relationships among milk production, energy balance, plasma analytes, and reproduction in Holstein-Friesian cows. *J. Dairy Sci.* **2007**, *90*, 649–658. [CrossRef]

4. Veerkamp, R.; Oldenbroek, J.; Van Der Gaast, H.; Van Der Werf, J. Genetic correlation between days until start of luteal activity and milk yield, energy balance, and live weights. *J. Dairy Sci.* **2000**, *83*, 577–583. [CrossRef]

5. Grummer, R.R.; Mashek, D.G.; Hayirli, A. Dry matter intake and energy balance in the transition period. *Vet. Clin. Food Anim. Pract.* **2004**, *20*, 447–470. [CrossRef]

6. McNamara, J. Regulation of Adipose Tissue Metabolism in Support of Lactation. *J. Dairy Sci.* **1991**, *74*, 706–719. [CrossRef]
7. Mohamed-Ali, V.; Pinkney, J.; Coppack, S. Adipose tissue as an endocrine and paracrine organ. *Int. J. Obes. Relat. Metab. Disord.* **1998**, *22*, 1145–1158. [CrossRef]

8. Trayhurn, P.; Beattie, J.H. Physiological role of adipose tissue: White adipose tissue as an endocrine and secretory organ. *Proc. Nutr. Soc.* **2001**, *60*, 329–339. [CrossRef]

9. Guerre-Millo, M. Adipose tissue hormones. *J. Endocrinol. Investig.* **2002**, *25*, 855–861. [CrossRef]

10. Zachut, M. Defining the adipose tissue proteome of dairy cows to reveal biomarkers related to peripartum insulin resistance and metabolic status. *J. Proteome Res.* **2015**, *14*, 2863–2871. [CrossRef]

11. Zachut, M.; Kra, G.; Livshitz, L.; Portnick, Y.; Yakoby, S.; Friedlander, G.; Levin, Y. Proteome dataset of subcutaneous adipose tissue obtained from late pregnant dairy cows during summer heat stress and winter seasons. *Data Brief* **2017**, *12*, 535–539. [CrossRef] [PubMed]

12. Coyral-Castel, S.; Faverdin, P.; Rame, C.; Freret, S.; Guillaume, D.; Fritz, S.; Dupont, J. Significant differences in fertility between dairy cows selected for one QTL located on bovine chromosome 3 are not attributable to energy balance, although eating behaviour is affected. *Animal* **2013**, *7*, 610–617. [CrossRef] [PubMed]

13. Wessel, D.; Flügge, U. A method for the quantitative recovery of protein in dilute solution in the presence of detergents and lipids. *Anal. Biochem.* **1984**, *138*, 141–143. [CrossRef]

14. Piras, C.; Soggiu, A.; Bonizzi, L.; Greco, V.; Ricchi, M.; Arrigoni, N.; Bassols, A.; Urbani, A.; Roncada, P. Identification of immunoreactive proteins of *Mycobacterium avium* subsp. paratuberculosis. *Proteomics* **2015**, *15*, 813–823. [CrossRef] [PubMed]

15. Piras, C.; Guo, Y.; Soggiu, A.; Chanrot, M.; Greco, V.; Urbani, A.; Charpigny, G.; Bonizzi, L.; Roncada, P.; Humblotb, P. Changes in protein expression profiles in bovine endometrial epithelial cells exposed to *E. coli* LPS challenge. *Mol. Biosyst.* **2017**, *13*, 392–405. [CrossRef] [PubMed]

16. Mi, H.; Muruganujan, A.; Casagrande, J.T.; Thomas, P.D. Large-scale gene function analysis with the PANTHER classification system. *Nat. Protoc.* **2013**, *8*, 1551–1566. [CrossRef]

17. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein–protein interaction networks, integrated over the tree of life. *Nucleic Acids Res.* **2014**, *43*, 447–452. [CrossRef]

18. Leonardini, A.; Laviola, L.; Perrini, S.; Natalicchio, A.; Giorgino, F. Cross-Talk between PPAR and Insulin Signaling and Modulation of Insulin Sensitivity. *PPAR Res.* **2010**. [CrossRef]

19. Ruan, X.; Li, Z.; Zhang, Y.; Yang, L.; Pan, Y.; Wang, Z.; Feng, G.-S.; Chen, Y. Apolipoprotein A-I possesses an anti-obesity effect associated with increase of energy expenditure and up-regulation of UCP1 in brown fat. *J. Cell. Mol. Med.* **2011**, *15*, 763–772. [CrossRef]

20. Folnožić, I.; Turk, R.; Duričić, D.; Vince, S.; Pleadin, J.; Flegar-Meštrić, Z.; Valpotić, H.; Dobranić, T.; Gračner, D.; Samardžija, M. Influence of Body Condition on Serum Metabolic Indicators of Lipid Mobilization and Oxidative Stress in Dairy Cows During the Transition Period. *Reprod. Domest. Anim.* **2015**, *50*, 910–917. [CrossRef]

21. Karavia, E.A.; Papachristou, D.J.; Liopeta, K.; Triantaphyllidou, I.-E.; Dimitrakopoulos, O.; Kypreos, K.E. Apolipoprotein AI modulates processes associated with diet-induced nonalcoholic fatty liver disease in mice. *Mol. Med.* **2012**, *18*, 901–912. [CrossRef] [PubMed]

22. Turk, R.; Juretić, D.; Gereš, D.; Svetina, A.; Turk, N.; Flegar-Meštrić, Z. Influence of oxidative stress and metabolic adaptation on PON1 activity and MDA level in transition dairy cows. *Anim. Reprod. Sci.* **2008**, *108*, 98–106. [CrossRef] [PubMed]

23. Turk, R.; Podpečan, O.; Mrkun, J.; Kosec, M.; Flegar-Meštrić, Z.; Perkov, S.; Starič, J.; Robić, M.; Belić, M.; Zrimšek, P. Lipid mobilisation and oxidative stress as metabolic adaptation processes in dairy heifers during transition period. *Anim. Reprod. Sci.* **2013**, *141*, 109–115. [CrossRef] [PubMed]

24. Hertzel, A.V.; Bernlohr, D.A. The mammalian fatty acid-binding protein multigene family: Molecular and genetic insights into function. *Trends Endocrinol. Metab.* **2000**, *11*, 175–180. [CrossRef]

25. Furuhashi, M.; Saitoh, S.; Shimamoto, K.; Miura, T. Fatty acid-binding protein 4 (FABP4): Pathophysiological insights and potent clinical biomarker of metabolic and cardiovascular diseases. *Clin. Med. Insights. Cardiol.* **2014**. [CrossRef]

26. Furuhashi, M.; Hotamisligil, G.S. Fatty acid-binding proteins: Role in metabolic diseases and potential as drug targets. *Nat. Rev. Drug Discov.* **2008**, *7*, 489–503. [CrossRef]

27. Veerkamp, J.; Van Moerkerk, H. *Cellular Fatty Acid-Binding Proteins II*; Springer: New York, NY, USA, 1993; pp. 101–106.
28. De Vries, M.; Veerkamp, R. Energy balance of dairy cattle in relation to milk production variables and fertility. *J. Dairy Sci.* **2000**, *83*, 62–69. [CrossRef]

29. Zhao, Y.M.; Basu, U.; Dodson, M.V.; Basarb, J.A. Proteome differences associated with fat accumulation in bovine subcutaneous adipose tissues. *Proteome Sci.* **2010**, *8*, 14. [CrossRef]

30. Salameh, A.; Daquinag, A.C.; Staquicini, D.I.; An, Z.; Hajjar, K.A.; Pasqualini, R.; Arap, W.; Kolonin, M.G. Prohibitin/annexin 2 interaction regulates fatty acid transport in adipose tissue. *JCI Insight* **2016**, *1*, e86351. [CrossRef]

31. Gerke, V.; Creutz, C.E.; Moss, S.E. Annexins: Linking Ca²⁺ signalling to membrane dynamics. *Nat. Rev. Mol. Cell Biol.* **2005**, *6*, 449–461. [CrossRef]

32. Kuhla, B.; Metges, C. *Energy and Protein Metabolism and Nutrition in Sustainable Animal Production*; Wageningen Academic Publishers: Wageningen, The Netherlands, 2013; pp. 213–221.

33. Xu, W.; Vervoort, J.; Saccenti, E.; van Hoeij, R.; Kemp, B.; van Knegsel, A. Milk Metabolomics Data Reveal the Energy Balance of Individual Dairy Cows in Early Lactation. *Sci. Rep.* **2018**, *8*, 15828. [CrossRef] [PubMed]

<figure><img src="image_17.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).