# Mammary gland transcriptome and proteome modifications by nutrient restriction in early lactation Holstein cows challenged with intramammary lipopolysaccharide

Karol Pawlowski, José Pires, Yannick Faulconnier, Christophe Chambon, Pierre Germon, Céline Boby, Christine Leroux

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

## To cite this version:

Karol Pawlowski, José Pires, Yannick Faulconnier, Christophe Chambon, Pierre Germon, et al., Mammary gland transcriptome and proteome modifications by nutrient restriction in early lactation Holstein cows challenged with intramammary lipopolysaccharide. International Journal of Molecular Sciences, 2019, 20 (5), 19 p. 10.3390/ijms20051156. [https://doi.org/10.3390/ijms20051156](https://doi.org/10.3390/ijms20051156)

HAL Id: hal-02625031

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

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

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

Article

# Mammary Gland Transcriptome and Proteome Modifications by Nutrient Restriction in Early Lactation Holstein Cows Challenged with Intra-Mammary Lipopolysaccharide

Karol Pawłowski <fnref n="1" />,<fnref n="2" />, José A.A. Pires <fnref n="1" />, Yannick Faulconnier <fnref n="1" />, Christophe Chambon <fnref n="3" />, Pierre Germon <fnref n="4" />, Céline Boby <fnref n="1" /> and Christine Leroux <fnref n="1" />,<fnref n="5" />,*

Received: 19 November 2018; Accepted: 27 February 2019; Published: 6 March 2019

**Abstract:** The objective is to study the effects of nutrient restrictions, which induce a metabolic imbalance on the inflammatory response of the mammary gland in early lactation cows. The aim is to decipher the molecular mechanisms involved, by comparing a control, with a restriction group, a transcriptome and proteome, after an intra-mammary lipopolysaccharide challenge. Multi-parous cows were either allowed *ad libitum* intake of a lactation diet (n = 8), or a ration containing low nutrient density (n = 8; 48% barley straw and dry matter basis) for four days starting at 24 ± 3 days in milk. Three days after the initiation of their treatments, one healthy rear mammary quarter of 12 lactating cows was challenged with 50 µg of lipopolysaccharide (LPS). Transcriptomic and proteomic analyses were performed on mammary biopsies obtained 24 h after the LPS challenge, using bovine 44K microarrays, and nano-LC-MS/MS, respectively. Restriction-induced deficits in energy, led to a marked negative energy balance (41 versus 97 ± 15% of Net Energy for Lactation (NEL) requirements) and metabolic imbalance. A microarray analyses identified 25 differentially expressed genes in response to restriction, suggesting that restriction had modified mammary metabolism, specifically β-oxidation process. Proteomic analyses identified 53 differentially expressed proteins, which suggests that the modification of protein synthesis from mRNA splicing to folding. Under-nutrition influenced mammary gland expression of the genes involved in metabolism, thereby increasing β-oxidation and altering protein synthesis, which may affect the response to inflammation.

**Keywords:** nutrigenomics; transcriptomic analysis; proteomic analysis; mammary gland; cows; inflammation

## 1. Introduction

Milk is synthesized in mammary glands (MG) involving a large number of genes, the expression of which is modulated at a nutritional level [1] and by the health status [2]. Mastitis is the inflammatory response of the mammary gland to pathogens. This pathology is one of the most prevalent disease and has considerable economic impact due to decreased milk production, discarded

<footnote n="1">¹ Université Clermont Auvergne, INRA, VetAgro Sup, UMR Herbivores, F-63122 Saint-Genès-Champanelle, France; karol_pawlowski@sggw.pl (K.P.); jose.pires@inra.fr (J.P.); yannick.faulconnier@inra.fr (Y.F.); celine.boby@inra.fr (C.B.)</footnote>

<footnote n="2">² Department of Pathology and Veterinary Diagnostics, Faculty of Veterinary Medicine, Warsaw University of Life Sciences, 02-776 Warsaw, Poland</footnote>

<footnote n="3">³ INRA, INRA, Plateforme d'Exploration du Métabolisme, composante protéomique PFEMcp), F-63122 Saint-Genès Champanelle, France; christophe.chambon@inra.fr</footnote>

<footnote n="4">⁴ INRA Val de Loire, UMR ISP, F-37380 Nouzilly, France; pierre.germon@inra.fr</footnote>

<footnote n="5">⁵ Department of Food Science and Technology, University of California Davis, Davis, CA 95616, USA</footnote>

<footnote>* Correspondence: christine.leroux@inra.fr</footnote>
milk, cost of veterinary services, and culling [3]. Mastitis is caused by various microorganisms (bacteria, fungus, and viruses). Gram-positive bacteria, such as *Staphylococcus aureus* and *Streptococcus uberis* might cause persistent infections, with pathogens surviving inside host cells [4,5], whereas Gram-negative coliform bacteria, such as *Escherichia coli* most frequently cause an acute inflammation and, eventually, severe mastitis with clinical signs [6]. Cows are particularly susceptible to *E. coli* MG inflammation during the periparturient period, due to altered immune function [7]. Early lactation is often associated with metabolic disorders related to stress, energy deficit and mobilization of the body reserves, hypocalcemia, and metritis, which are likely to influence immune function [8–10]. Negative energy balance (NEB) affects the inflammatory response, which could be due to changes in the metabolic milieu, such as an increased concentration of circulating ketone bodies [8,11]. Undernutrition, however, had minor effects on the response to lipopolysaccharide (LPS) and *S. uberis* challenges in mid-lactation dairy cows [12,13]. The current study is conducted during early lactation Holstein cows, which is characterized by enhanced metabolic deviations in response to NEB that might influence immune system function.

Transcriptome and proteome profiling techniques are available to study inflammation-related changes and enhance the understanding of host-pathogen interactions. Previous research has employed reproducible protocols to challenge lactating cow udders with live pathogens, such as *E. coli*, *S. aureus*, or using LPS to induce an acute inflammatory response [12–15]. High-throughput gene expression technologies, as transcriptomic analyses, have been used to decipher the molecular mechanisms of MG response to inflammation [16], showing that LPS is a strong stimulator of gene expression in an inflammatory response [17]. Transcriptomic analyses showed that MG differentially expressed genes (DEG) 24 h after *E. coli* infection in early-lactating cows, showing the associations of up- and down-regulated genes, respectively with immune response functions and fat metabolism [18]. A study at mid-lactation showed a massive effect of *E. coli* infection on gene expression at 18 h post-infection in infected quarters, but also having an effect on gene expression in neighboring quarters [19]. Furthermore, hierarchic clustering of DEG, showed a sharp separation of the infected and the control group [19]. Similar results were observed in primary bovine mammary epithelial cells after challenging the *E. coli* or *S. aureus* [20–22]. The modification of bovine milk proteome were reported during inflammation, showing an increase in the concentration of proteins of blood serum origin as serum albumin is made up of antimicrobial peptides [23,24]. To the best of our knowledge, the effects of negative energy balance on the responses to acute inflammation have not been studied at a protein expression level in the MG of early lactation cows.

We hypothesized that aggravated undernutrition in early lactating cows modifies the inflammatory response at mRNA and protein levels. Thus, the objective of this study is to evaluate the effects of undernutrition, and the resulting metabolic imbalance on MG transcriptome and proteome in early lactation Holsteins, challenged with intra-mammary LPS. We used complementary approaches to study the effects of undernutrition with a first targeted study of genes involved in the inflammation response, using RT-qPCR, then two global analyses using transcriptomic and proteomic analyses. These two last approaches allowed us to investigate the holistic response of MG. The effects of LPS are already well-described in the literature, thus attention is focused on the effects of restriction during inflammation, by comparisons with MG transcriptome, and the proteome of control versus under-fed early lactating cows. This study was performed with the aim of increasing the knowledge-base on the effects of NEB occurring in early lactation which is the time of high risk of mastitis. A better understanding can help to design news strategies to prevent MG inflammation.

## 2. Results

### 2.1. *Dietary and Inflammatory Challenges Influenced Milk and Blood Composition*

Prior to diet treatments at 24 ± 3 days in milk (DIM), no differences were observed for intake, milk yield, composition and component yield, Net Energy for Lactation (NEL) balance, body score (BS), body weight (BW), plasma metabolite, and insulin concentrations. Feeding the ration containing
48% of straw (restricted group: REST) induced an immediate depression of dry matter intake (DMI) and decreased the energy balance from 5.2 ± 8.9 to -67.2 ± 18.9 MJ/day one day before (corresponding to day 23), and on the last day of restriction (day 27), respectively. The plasma concentrations of glucose and insulin decreased, whereas non-esterified fatty acids (NEFA) and beta-hydroxy butyrate (BHBA) increased dramatically in REST (Supplementary Table S1). The 96 h of nutrient restriction decreased milk yield from 37.9 to 22.4 kg/day (p < 0.001) and milk protein yield from 1.12 to 0.62 kg/day (p < 0.005) in REST. Milk fat percentage increased during feed restriction in REST from 4 to 5.5% (p < 0.001) and returned to pre-restriction concentrations on the same day of re-feeding a regular diet. In CONT cows, milk fat content was greatest only during the 48 h following LPS injection (p < 0.05) compared with all other DIM. These variables were unchanged in the control (CONT) cows [25]. Within 2 to 6 h hours after injection with lipopolysaccharide, we noticed the edema of the challenged quarter of all cows, an increase in rectal temperature up to 39.5 °C (temperature increment was +2.1 ± 0.15 °C). The effect of inflammation also was confirmed by milk somatic cell count (SCC). The day before the LPS challenge, whole udder composite milk (from PM and AM milking) SCC was 78 000/mL and 92 000/mL (p = 0.51) in CONT, and in REST, respectively. The SCC response was greater in REST compared to CONT cows (6919 versus 1956 × 1000 per mL, respectively) in composite milk samples from the two milkings that followed LPS injection. SCC returned to pre-LPS counts, within less than 7 days post-LPS challenge and biopsy. Moreover, quarter milk IL-8, IL1-β, TNF-alpha, and CXCL3 at time zero (before LPS challenge) did not differ between CONT and REST, but their concentrations increased in response to LPS in both groups. This data shows that these indicators of mammary inflammation did not differ between CONT and REST before or after the challenge.

## 2.2. Effects of Under-nutrition on Expression of Genes Involved in Inflammation Response by RT-qPCR Analysis

RT-qPCR analysis was performed to quantify candidate genes (CCL5, LAP, RBP4, IL8, IL1, STAT3, CD36, and TAP) chosen on the basis of their implication in inflammation [26] Expression of these genes in mammary gland (MG) did not differ between REST and CONT (p ≥ 0.1) 24 h after the LPS challenge, except for the defensin Tracheal Antimicrobial Peptide (TAP) gene which tended to decrease in REST (p = 0.07). The expression of INSIG1, and CSN2 genes, which are involved in the biosynthesis of milk components and linked to MG metabolism did not differ between CONT and REST (Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Effects of nutrient restriction and intra-mammary lipopolysaccharide (LPS) challenge on gene mRNA expression quantified by RT-qPCR and presented as ΔCT. Comparison of the gene expression does not show a difference between control (CONT; *n* = 6) and restricted (REST; *n* = 6)</figcaption></figure>
Holstein cows ($p \ge 0.1$). The expression of the *TAP* gene tended to differ ($p = 0.07$). *UXT2*, *CLN3*, and *EIF3K* were used as housekeeping genes.

## 2.3. Microarray Analysis

Mammary gene expression analyzed by a microarray assay allowed the identification of 33 differentially expressed genes, including 25 known genes (corrected $p_{adj} < 0.05$), between CONT and REST, 24 h after the inflammatory challenge by LPS (Table 1). The expression increased for 19 and decreased for 6 genes in REST compared with CONT. All these DEG presented a fold change (FC) greater than 1.4, with two genes (*PDK4* and *SLC25A34*) presenting an FC greater than 4. Gene ontology and function analyses revealed that most DEG are involved in metabolism, including the regulation of fatty acid (FA) oxidation, glucose, and protein metabolism, and in immune responses (Figure 2). The results obtained, using Pathway Studio® software, were consistent with those from Panther software. We focused on the most represented functions, in particular, those involved in metabolism and immune response. We identified DEG in FA and glucose metabolism (*CPT1A*, *PDK4*, *PFKFB4*), carnitine shuttle (*SLC25A20*, *CPT1A*, *SLC25A34*), regulation of cellular ketone metabolic process (*PDK4*), and the key genes in those processes. A number of genes involved autophagy (*PFKFB4*, *DNNED*) and immune function (*PGLYRP3*, *KLF13*, *PLEKHA2*, *WC7*, *TRIB2*, *CXCR7*, and *MBP*) processes also were altered.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Main biological processes of differentially expressed genes in mammary glands (MG) of underfed (REST) versus control (CONT) early lactation cows during an acute inflammation identified by transcriptomic analysis. Each percentage was indicated for each biological processes. Bioinformatics analyses were performed using Panther and Pathway Studio® software.</figcaption></figure>

<table><caption>**Table 1.** Differentially expressed genes (alphabetic classification) in MG of early lactation Holstein cows in response to undernutrition and intra-mammary lipopolysaccharide (LPS) challenge. Transcriptomic comparison of control (CONT; *n* = 6) and restricted diet (REST; *n* = 6). Normalized microarray data were analyzed using GeneSpring software and moderated *t*-tests with Westfall-Young correction. $p_{adj} \le 0.05$ and ranged from 0.01 to 0.04. FC were >1.4.</caption><thead><tr><th>Gene Symbol</th><th>Accession Number</th><th>Gene Name</th><th>Fold change</th></tr></thead><tbody><tr><td>PDK4</td><td>NM_001101883.1</td><td>[Pyruvate Dehydrogenase (Acetyl-Transferring)] Kinase Isozyme 4, Mitochondrial</td><td>6.80</td></tr><tr><td>SLC25A34</td><td>NM_001034497.2</td><td>Solute Carrier Family 25 Member 34</td><td>4.63</td></tr><tr><td>PGLYRP3</td><td>XM_010826801.2</td><td>Peptidoglycan Recognition Protein 3</td><td>3.31</td></tr><tr><td>CPT1A</td><td>NM_001304989.1</td><td>Carnitine O-Palmitoyltransferase 1</td><td>2.77</td></tr><tr><td>NEFL</td><td>NM_174121.1</td><td>Neurofilament Light Polypeptide</td><td>2.76</td></tr></tbody></table>
<table><tr><td>KNTC1</td><td>NM_001192091.</td><td>Kinetochore-associated protein 1</td><td>2.58</td></tr><tr><td>MBP</td><td>NM_001206674.1</td><td>Myelin basic protein</td><td>2.37</td></tr><tr><td>DENND3</td><td>XM_010822660.2</td><td>DENN domain-containing protein 3</td><td>2.33</td></tr><tr><td>TFAP2D</td><td>NM_001192329.1</td><td>Transcription factor AP-2-delta</td><td>2.29</td></tr><tr><td>PFKFB4</td><td>NM_001192835.1</td><td>6-Phosphofructo-2-kinase/Fructose-2,6-bisphosphatase 4</td><td>2.04</td></tr><tr><td>TRIB2</td><td>NM_178317.3</td><td>Tribbles homolog 2</td><td>2.03</td></tr><tr><td>WC-7</td><td>NM_001281912.1</td><td>WC1 isolate CH149</td><td>1.85</td></tr><tr><td>TMEM50B</td><td>NM_001034786.2</td><td>Transmembrane protein 50B</td><td>1.78</td></tr><tr><td>ALAD</td><td>NM_001014895.1</td><td>Delta-aminolevulinic acid dehydratase</td><td>1.77</td></tr><tr><td>LOC517144</td><td>XM_002699133.1</td><td>Putative olfactory receptor 10D3</td><td>1.74</td></tr><tr><td>ESRRA</td><td>NM_001191373.2</td><td>Steroid hormone receptor ERR1</td><td>1.73</td></tr><tr><td>SLC25A20</td><td>NM_001077936.2</td><td>Mitochondrial carnitine/Acylcarnitine carrier protein</td><td>1.70</td></tr><tr><td>PLEKHA2</td><td>NM_001035383.1</td><td>Pleckstrin homology domain-containing family A member 2</td><td>1.69</td></tr><tr><td>KLF13</td><td>NM_001083533.1</td><td>Krueppel-like factor 13</td><td>1.67</td></tr><tr><td>ARMC1</td><td>NM_001015594.2</td><td>Armadillo repeat-containing protein 1</td><td>-1.43</td></tr><tr><td>RPL-37A</td><td>NM_001035008.2</td><td>Similar to 60S ribosomal protein L37a</td><td>-1.76</td></tr><tr><td>BBS9</td><td>NM_001192853.1</td><td>Protein PTHB1 / Bardet-bied syndrom 9</td><td>-1.79</td></tr><tr><td>SLC9A7</td><td>XM_015470441.1</td><td>Sodium/Hydrogen exchanger 7</td><td>-1.99</td></tr><tr><td>CXCR7</td><td>NM_001098381.2</td><td>Atypical chemokine receptor 3</td><td>-2.10</td></tr><tr><td>OOEP</td><td>NM_001077869.2</td><td>Oocyte-expressed protein homolog</td><td>-2.52</td></tr></table>

## 2.4. LC MS/MS Proteomic Analysis

We identified 1475 proteins, 967 of which were validated with more than two peptides and considered for further investigation. Fifty three proteins were differentially expressed proteins, 10 were upregulated and 43 were downregulated in REST (p < 0.05; Table 2; Supplemental File 2). Classification of DEP highlighted proteins involved in the immune process, metabolism (regulation of protein catabolic and carbohydrate metabolic processes) and cell functioning (such as RNA splicing, translation, or regulation of cell adhesion). Proteins involved in apoptosis are also identified (Figure 3). DEP are involved in protein folding and post-translational modifications (GANAB, PDIA3, RPN2, RPN1, CCT4, PDIA4, and PPIB), protein catabolic process (PSMD2 and PPP2CA), carbohydrate metabolism (PAPSS1, RPS27A, GANAB), synthesis of immunoglobulins (F1MH40, F1MLW8), and regulation of inflammatory response (PDIA3, PSMA3, PSMD2, PPP2CA, RPS2, RPL10, RPS15, CASP6, PCBP2, and STAT5A). Many DEP are involved in RNA splicing (HNRNPH1, DHX9, HNRNPC, YBX1, PCBP2, SNRPA1, HNRNPA3, and PPP2CA) or translation processes (RPS27A, RPL10, RPS15, RPS2, EIF3H, RPN2, RPN1, FARSB, and NACA). Additional analysis by Uniprot software revealed similar classifications confirming them.

<figure><img src="image_6.png" /></figure>
**Figure 3.** Main biological processes of differentially expressed proteins in MG of underfed (REST) versus control (CONT) early lactation cows during acute inflammation. Each percentage was indicated for each biological processes. Bioinformatics analyses were performed using Uniprot and Pathway Studio® software.

<table><caption>**Table 2.** The list of differentially expressed proteins (up- and down-regulated then alphabetic classification) 24 h after inflammation challenge by LPS in MG of underfed (REST; *n* = 6) compared with control (CONT; *n* = 6) cows. Proteins were analyzed with Progenesis LC-MS software v4.1 (Nonlinear Dynamics). The minimum mascot score validation for one peptide was 31 with a rate of false discovery <1%.</caption><thead><tr><th>Gene Symbol</th><th>Accession Number</th><th>Protein Name</th><th>Fold Change</th></tr></thead><tbody><tr><td>Unknown</td><td>F1MLW8_BOVIN</td><td>Uncharacterized protein / Immunoglobulin light chain, lambda</td><td>5.5</td></tr><tr><td>SERPINA3</td><td>G8JKW7_BOVIN</td><td>Uncharacterized protein</td><td>2.4</td></tr><tr><td>SERPINA3-5</td><td>SPA35_BOVIN</td><td>Serpin A3-5</td><td>2.0</td></tr><tr><td>Unknown</td><td>Q1RMN8</td><td>TREMBL:Q1RMN8 (Bos taurus) Similar to Immunoglobulin lambda-like polypeptide 1</td><td>1.9</td></tr><tr><td>ALDH18A1</td><td>Q2KJH7_BOVIN</td><td>Aldehyde dehydrogenase 18 family, member A1</td><td>1.7</td></tr><tr><td>Unknown</td><td>F1MH40_BOVIN</td><td>Uncharacterized protein</td><td>1.5</td></tr><tr><td>HBA</td><td>HBA_BOVIN</td><td>Hemoglobin subunit alpha</td><td>1.3</td></tr><tr><td>HBB</td><td>HBB_BOVIN</td><td>Hemoglobin subunit beta</td><td>1.3</td></tr><tr><td>PCBP2</td><td>Q3SYT9_BOVIN</td><td>Poly(RC) binding protein 2</td><td>1.2</td></tr><tr><td>ARPC2</td><td>ARPC2_BOVIN</td><td>Actin-related protein 2/3 complex subunit 2</td><td>1.1</td></tr><tr><td>CCT4</td><td>TCPD_BOVIN</td><td>T-complex protein 1 subunit delta</td><td>-1.2</td></tr><tr><td>DDX17</td><td>A7E307_BOVIN</td><td>DDX17 protein</td><td>-1.2</td></tr><tr><td>GANAB</td><td>F1N6Y1_BOVIN</td><td>Uncharacterized protein</td><td>-1.2</td></tr><tr><td>HNRNPH1</td><td>E1BF20_BOVIN</td><td>Uncharacterized protein</td><td>-1.2</td></tr><tr><td>HNRPC</td><td>Q3SX47_BOVIN</td><td>Heterogeneous nuclear ribonucleoprotein C (C1/C2)</td><td>-1.2</td></tr><tr><td>LMAN2</td><td>A6QP36_BOVIN</td><td>LMAN2 protein</td><td>-1.2</td></tr><tr><td>NACA</td><td>NACA_BOVIN</td><td>Nascent polypeptide-associated complex subunit alpha</td><td>-1.2</td></tr><tr><td>PDIA4</td><td>PDIA4_BOVIN</td><td>Protein disulfide-isomerase A4</td><td>-1.2</td></tr><tr><td>PPP2CA</td><td>PP2AA_BOVIN</td><td>Serine/threonine-protein phosphatase 2A catalytic subunit alpha isoform</td><td>-1.2</td></tr><tr><td>PSMA3</td><td>PSA3_BOVIN</td><td>Proteasome subunit alpha type-3</td><td>-1.2</td></tr><tr><td>RPS2</td><td>RS2_BOVIN</td><td>40S ribosomal protein S2</td><td>-1.2</td></tr><tr><td>RPS27A</td><td>RS27A_BOVIN</td><td>Ubiquitin-40S ribosomal protein S27a</td><td>-1.2</td></tr><tr><td>HNRNPA3</td><td>E1BEG2_BOVIN</td><td>Uncharacterized protein</td><td>-1.1</td></tr><tr><td>PDIA3</td><td>PDIA3_BOVIN</td><td>Protein disulfide-isomerase A3</td><td>-1.1</td></tr><tr><td>PRKCDBP</td><td>PRDBP_BOVIN</td><td>Protein kinase C delta-binding protein</td><td>-1.1</td></tr><tr><td>ACTR1A</td><td>F2Z4F0_BOVIN</td><td>Uncharacterized protein</td><td>-1.3</td></tr><tr><td>DHX9</td><td>DHX9_BOVIN</td><td>ATP-dependent RNA helicase A</td><td>-1.3</td></tr><tr><td>GNB1</td><td>GBB1_BOVIN</td><td>Guanine nucleotide-binding protein G(I)/G(S)/G(T) subunit beta-1</td><td>-1.3</td></tr><tr><td>KIF5B</td><td>F1N1G7_BOVIN</td><td>Kinesin-like protein</td><td>-1.3</td></tr><tr><td>PLBD2</td><td>PLBL2_BOVIN</td><td>Putative phospholipase B-like 2</td><td>-1.3</td></tr><tr><td>RPN1</td><td>A3KN04_BOVIN</td><td>Dolichyl-diphosphooligosaccharide--protein glycosyltransferase subunit 1</td><td>-1.3</td></tr><tr><td>SNRPA1</td><td>A6H788_BOVIN</td><td>SNRPA1 protein</td><td>-1.3</td></tr><tr><td>Unknown</td><td>F6PWD5_BOVIN</td><td>Uncharacterized protein (Fragment)</td><td>-1.3</td></tr><tr><td>CAPZA2</td><td>CAZA2_BOVIN</td><td>F-actin-capping protein subunit alpha-2</td><td>-1.4</td></tr><tr><td>FARSB</td><td>A8E4P2_BOVIN</td><td>FARSB protein</td><td>-1.4</td></tr><tr><td>HSPA4</td><td>E1BBY7_BOVIN</td><td>Uncharacterized protein</td><td>-1.4</td></tr><tr><td>PAPSS1</td><td>Q3T0J0_BOVIN</td><td>3'-phosphoadenosine 5'-phosphosulfate synthase 1</td><td>-1.4</td></tr><tr><td>PPIB</td><td>PPIB_BOVIN</td><td>Peptidyl-prolyl cis-trans isomerase B</td><td>-1.4</td></tr><tr><td>PSMD2</td><td>PSMD2_BOVIN</td><td>26S proteasome non-ATPase regulatory subunit 2</td><td>-1.4</td></tr><tr><td>RPS15</td><td>RS15_BOVIN</td><td>40S ribosomal protein S15</td><td>-1.4</td></tr><tr><td>STAT5A</td><td>STA5A_BOVIN</td><td>Signal transducer and activator of transcription 5A</td><td>-1.4</td></tr><tr><td>CASP6</td><td>CASP6_BOVIN</td><td>Caspase-6</td><td>-1.5</td></tr><tr><td>RPN2</td><td>RPN2_BOVIN</td><td>Dolichyl-diphosphooligosaccharide-protein glycosyltransferase subunit 2</td><td>-1.5</td></tr><tr><td>NIPSNAP3A</td><td>G3X6L8_BOVIN</td><td>Uncharacterized protein</td><td>-1.6</td></tr><tr><td>EIF3H</td><td>EIF3H_BOVIN</td><td>Eukaryotic translation initiation factor 3 subunit H</td><td>-1.7</td></tr><tr><td>YBX1</td><td>YBOX1_BOVIN</td><td>Nuclease-sensitive element-binding protein 1</td><td>-1.7</td></tr><tr><td>COPS7A</td><td>F6QE33_BOVIN</td><td>Uncharacterized protein</td><td>-1.8</td></tr></tbody></table>
<table><tr><td>PDLIM5</td><td>G3MY19_BOVIN</td><td>Uncharacterized protein</td><td>-1.8</td></tr><tr><td>MYBBP1A</td><td>E1BKX3_BOVIN</td><td>Uncharacterized protein</td><td>-1.9</td></tr><tr><td>RPL10</td><td>RL10_BOVIN</td><td>60S ribosomal protein L10</td><td>-1.9</td></tr><tr><td>EEA1</td><td>F1MN61_BOVIN</td><td>Uncharacterized protein (Fragment)</td><td>-2.0</td></tr><tr><td>ARPC1B</td><td>ARC1B_BOVIN</td><td>Actin-related protein 2/3 complex subunit 1B</td><td>-2.2</td></tr><tr><td>C789567</td><td>A6H7H3_BOVIN</td><td>LOC789567</td><td>-2.5</td></tr></table>

### 3. Discussion

This study assesses the effects of undernutrition and the resulting metabolic imbalance on the mammary gland (MG) inflammatory response in early lactation cows using a nutrigenomic approach. The effects of dietary treatments are confirmed by decreased intake, milk yield, energy balance in underfed (REST) cows, and by changes in blood metabolite and insulin concentrations (an increase in plasma NEFA and BHBA and a decrease in insulin and glucose concentrations; Supplementary Table S1; [25]). The inflammatory response to intra-mammary lipopolysaccharide is confirmed by clinical parameters, such as milk SCC, rectal temperature, and other classical clinical symptoms [25]. The effects of nutrient restriction and metabolic imbalance on the inflammatory response, at the RNA and protein levels, were evaluated using transcriptomic and proteomic analyses. These analyses were performed using MG samples were obtained by biopsies performed 24 h after the LPS challenge. We performed a single biopsy to avoid the potential interference of repetitive biopsies on the inflammatory response. Also, the adjacent quarter may not constitute a good control for the LPS challenged quarter, as inflammation cytokines may exert local effects and influence adjacent quarters [19]. However, one limitation of this design is that the present experimental design does not allow a kinetic data to follow the establishment of inflammation.

#### 3.1. Gene Expression Changes at mRNA Level

The study at the mRNA level was performed using two complementary approaches. The first one was a targeted approach to focus our attention on the inflammatory response, then a global analysis was performed using microarray study. We investigated the effects on candidate genes by RT-qPCR, the majority, which are involved in the immune response of interest for inflammation. The expression of candidate genes did not differ between REST and with the control (CONT) group (Figure 1). This result suggests that the expression of genes considered important in the inflammatory response [11,26–28] are not altered by the nutrient restriction in MG 24 h after LPS administration. This experimental design does not allow the ability to evaluate the modification of the expression of these genes during early inflammatory response. A tendency for decreased expression of the *TAP* gene in REST was observed. The product of *TAP* gene is a member of the family of small cationic peptides that have widespread antimicrobial activity; TAP is expressed by bovine mammary epithelial cells [29] and has a broad-spectrum activity against different strains of bacteria, including *E. coli* [30]. The upregulation of *TAP* gene expression in REST may constitute a protective mechanism against pathogens.

To complete the candidate gene analyses, a global gene expression approach, using a bovine microarray was used to assess the molecular mechanisms underlying metabolic and inflammatory MG responses, potentially affected by undernutrition and negative energy balance. Transcriptomic analysis revealed 33 differentially expressed genes in MG 24 h after LPS challenge in REST compared to CONT. The number of DEG detected in our study are small compared to the research assessing the effects of inflammation on MG transcriptome [18,19,31]. This study does not compare normal versus inflamed MG, rather it aims to evaluate the effects of undernutrition during the inflammatory response, and therefore both REST and CONT were challenged with LPS. The DEG were classified in six functional classes. In this discussion, we mainly focus our attention on genes that play a role in the metabolic processes (FA, glucose and protein metabolism) and on those involved in immune function.

##### 3.1.1. DEG Involved in Metabolism
The classification of genes by bioinformatics analyses indicate that metabolic process (FA oxidation, glucose, and protein metabolisms) is the class most altered by undernutrition after an LPS challenge conditions. Among the genes presenting the highest fold change (between 2.7 and 6.8) are *PDK4*, *CPT1*, and *SLC25A34*, which are involved in glucose and FA metabolism. *PDK4* plays a key role by inhibiting the pyruvate dehydrogenase complex. This inhibition prevents the formation of acetyl-coenzyme A from pyruvate [32], resulting in an expected decrease in glucose and an increase in fat utilization in response to prolonged undernutrition [33]. The upregulation of *PDK4* was also found in leucocytes of underfed ewes, however, its expression was downregulated during an intramammary inflammatory challenge [34]. The large increase of *PDK4* expression in REST is in agreement with decreased insulinemia and with the upregulation of *ESRRA*. The gene *ESRRA* upregulates *PDK4* expression [35]. Both genes spare glucose and promote FA β-oxidation, therefore, their upregulation in REST would allow MG to shift the metabolic pathways from glycolysis to β-oxidation. The upregulation of *CPT1* gene expression would also promote the β-oxidation [36,37], since it is a rate-limiting step of FA entry in mitochondria [38]. This is in line with the increased expression of the *CPT1* gene observed in whole blood transcriptome of underfed dairy sheep [34]. The increase in β-oxidation is further supported by an upregulation of *SLC25A20* and *SLC25A34* in REST, which are two members of *SLC25* mitochondrial carrier family. *SLC25A20* transports carnitine and carnitine-FA complexes across the inner mitochondrial membrane. *SLC25A34* is supposed to act in a similar way, but its exact function still is not known totally [39]. MG seems to spare glucose (downregulating glycolysis) and promote FAs as an energy source (upregulating β-oxidation; Figure 4) in order to adapt to underfeeding. The mammary expression of genes involved in lipid metabolism is also modified in comparison with NEB (induced by caloric restriction) and the positive energy balance of cows after the peak in lactation [40]. Interestingly, the downregulation of genes linked to fat metabolism (FA biosynthesis) is observed 24 h after *E. coli* infection in MG of lactating cows [18]. This suggests that inflammation downregulates the FA biosynthesis and may increase the use of preformed FA, derived from other sources, such as from the mobilization of adipose tissue.

These results suggest that energy metabolism modifications, in response to inflammation, are more marked in REST than in CONT, probably due to the expected limited availability of nutrients to support an acute inflammation in REST.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Upregulation of genes detected at mRNA level highlights a potential increase of fatty acid oxidation. Red boxes indicate upregulated genes. Blue and green boxes represent the nucleus, and mitochondria, respectively. *Carnitine O-palmitoyltransferase 1 (CPT1A)* encodes carnitine O-palmitin O-palmitoyltransferase 1, *Mitochondrial Carnitine/Acylcarnitine Carrier Protein (SLC25A20)* and *Solute Carrier Family 25 Member 34 (SLC25A34)* are two members of the solute carrier family 25. *Steroid Hormone Receptor ERR1 (ESRRA)* and *Pyruvate Dehydrogenase Acetyl-Transferring Kinase Isozyme 4 (PDK4)* encode a steroid hormone receptor and a pyruvate dehydrogenase, respectively.</figcaption></figure>
### 3.1.2. DEG Involved in Immune Response

The objective of this study was to identify the effects of an aggravated NEB on MG transcriptomic responses to inflammation. The experiment was not designed to compare gene expression in normal versus LPS challenged MG, which have been previously reported [2,11,18,31]. Microarray analysis shows that, the immune response together with the metabolic processes, are the main biological processes modified by undernutrition after LPS challenge (Figure 2), with modifications of genes different from those selected for RT-qPCR analyses. However, immune response was not the main biological process to be affected by restriction during inflammation. Among the upregulated genes is *PGLYRP3*, which belongs to a family of innate immunity pattern recognition molecules that are activated by LPS and bactericidal and bacteriostatic properties and are activated by LPS [41]. When the invading bacteria survive, neutrophil infiltration is replaced with T and B lymphocytes and monocytes [42]. The upregulation in REST of *KLF13*, *PLEKHA2*, *WC7*, and *MBP*, involved in the immune response by activating T and B lymphocytes, therefore, is in line with the expect recruitment of leukocytes by MG (Figure 5). Additionally, the upregulation of *PLEKHA2* in REST, a gene involved in the cell adhesion process [43], suggests that there is an increased migration of B leucocytes. Together, the upregulation of these genes suggests a different nature or a higher response to LPS stimuli in REST compared with CONT. Nevertheless, the deregulation of *TRIB2* and *CXCR7* suggests that the immune process could be impaired. Indeed, both genes participate in the activation of immune cells and influence IL-8 production, a chemokine upregulated in response to infection [10], where the concentration increased in milk, within 4 h after LPS infusion [25], and within 16–24 h after experimental infection with different strains of *E.coli* or LPS infusion [44,45]. The upregulation of *TRIB2* and the downregulation of *CXCR7* genes, however, suggests a potential IL-8 production alteration in response to inflammation in REST. These results contrast with the expected inflammatory response and could be a sign of deficient immune function under exacerbated NEB. Taken together, differences in gene expression might suggest a modified resolution of inflammation in response to LPS, due to aggravated NEB. During the course of an experiment with LPS, the inflammatory response usually declines within 24 h [44]. The REST cows might have experienced difficulties in restoring the MG homeostasis by 24 h after LPS challenge, due to the metabolic changes inherent in nutrient deficiency. However, this conception needs more detailed investigation, with a kinetic analysis, to be confirmed.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** The potential actions of the differentially expressed genes identified in the comparison of mammary transcriptomes of restricted versus control cows during LPS challenge were from bioinformatics analyses and literature.</figcaption></figure>

### 3.2. Gene Expression at Protein Level
### 3.2.1. Proteins Involved in Protein Synthesis

Among the 53 differentially expressed proteins (DEP) in REST compared to CONT, 43 were downregulated. Most of these are involved in RNA and protein metabolism, with roles that vary from RNA splicing to translation. The downregulation of proteins involved in the splicing process, such as HNRNPH1, HNRPC, HNRNPA3, PCBP2, YBX1, SNRPA1, and DHX9, suggests that splicing is impaired in the MG of REST cows. This could explain in part the reduced synthesis and secretion of milk protein in REST compared with CONT [25]. Moreover, altered splicing and translation mechanisms might have a profound influence on protein biochemical properties and, ultimately, alter immune response to pathogens. Among the four proteins belonging to the HRNPs family (HNRNPH1, HNRPC, HNRNPA3, PCBP2), the first three are RNA binding proteins associated with pre-mRNAs in the nucleus, influencing pre-mRNA processing as well as other aspects of mRNA metabolism and transport. The dysfunction of HRNPs is linked to different proliferative and degenerative diseases [46], but the role of these proteins in the inflammatory response is still not fully understood. Some members of this family are reported to ensure resolution of inflammation [47]. However, the role of HRNPs could be associated with a mammary remodeling due to the restriction. In our study, YBX1 and DHX9 are downregulated in REST and the downregulation of these two genes is associated with impaired inflammatory responses [48]. Additionally, the loss of PP2AC function causes severe immunological disorders in Treg cells [49]. Thus, the downregulation of all these proteins in REST (Figure 6) suggest a modified inflammatory response in underfed early lactation cows.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Identified differentially expressed genes are presented in their potential role in crucial steps of the control of gene expression. Green boxes indicate genes which are downregulated by nutritional restriction. Light and dark blue represent the nucleus and ribosomes, respectively. Grey boxes represent the endoplasmic reticulum. Differentially expressed proteins between restricted and control cows after four days of undernutrition for restricted cows during LPS challenge.</figcaption></figure>

A number of proteins involved in translation are downregulated in REST (Figure 6; RPS27A, RPS15, RPS2, RPL10, EIF3H, RPN2, RPN1, FARSB, and NACA). Four riboprotein family members (3 RPS and 1 RPL) are part of a ribosome. Interestingly, protein building ribosomes alone are shown to affect the other cell processes outside the ribosome like development, apoptosis, and aging during their altered expression levels [50]. Additionally, the decrease of RPN1 and RPN2, which catalyze co-translational N-glycosylation, suggests an impaired post-translational protein modification process in the REST group. This process may play an important role in the immune system by creating the glycans on an immune cell's surface that helps migration of the cell or by glycosylating the various immunoglobulins [51]. Elsewhere, it is reported that this process can be defective during glucose deficit, leading to a reduction of protein glycosylation and harmful accumulation of unfolded
proteins [52]. The decrease in RPN1 and RPN2, observed in the current study, could potentially lead to the creation of misfolded proteins in MG of REST cows.

Additionally, protein folding and its control might be modified in REST, due to the downregulation of chaperone proteins such as PDIA3, PDIA4, and CCT4 (Figure 6). PDIA3 and PDIA4 are part of a larger super-family of a disulfide isomerase family of endoplasmic reticulum proteins that catalyze protein folding [53]. PDIA3 contributes to the correct folding of glycoproteins [54]. The loss of PDI activity and the consequent accumulation of misfolded proteins are associated with chronic inflammation [54,55]. Moreover, PDIA3 is a structural component required for the stable assembly of the peptide-loading complex of the major histocompatibility complex class I pathway. Its activity seems to play a role in lymphocyte T and B function [56]. Added to its role in folding, PDIA4 promotes the immunoglobulin G intermolecular disulfide bonding and antibody assembly in vitro [57]. Because CCT4 assists in the folding of newly translated polypeptides, this function might have been altered in REST [58]. Overall, proteomic data strongly suggest that protein synthesis is impaired by undernutrition at different levels (translation, folding, and post-translation modifications). The modifications of protein metabolism might partially explain the lower milk protein yield from 1.12 to 0.62 kg/ observed during restriction.

### 3.2.2. DEP Involved in Immune Response

Undernutrition downregulated FARSB protein expression. The decrease of this protein is linked with impaired acute inflammation responses in mice [59], suggesting an impairment in immune system function. In contrast, there was an upregulation of proteins such as SERPINA3, SERPINA3-5, F1MLW8, and Q1RMN8. SERPINA is an acute-phase protein, whose concentration can rise during acute and chronic inflammation [60]. F1MLW8 and Q1RMN8 proteins are similar to immunoglobulin lambda and typical for B-cells, and are important for its maturation from pre-B cells to mature ones [61]. The increase of these four proteins in REST MG, suggests that the resolution inflammation process was delayed in REST, compared to CONT 24 h after LPS challenge, whereas it could be considered that, potentially, it has already resolved in CONT. This is in line with the reported peak in SCC at 12 h that declines 24 h after LPS challenge [44].

The decreased translation process and post-translational protein modification (folding and glycosylation), that is observed at the protein level, might alter protein synthesis and activate an unfolded protein response [52]. This role could also be suggested to affect the proteins involved in the immune response.

## 4. Materials and Methods

### 4.1. *Ethics Statement, Treatments and Sampling*

The cows were housed at the Herbivore Research Unit of INRA Research Center of Auvergne-Rhone-Alpes. Animal procedures were performed in compliance with Regional Animal Care Committee guidelines CEMEAA: Auvergne, French Ministry of Agriculture and European Union guidelines for animal research C2EA-02. All procedures were approved by the regional ethics committee on animal experimentation (APAFIS #2018062913565518). The animals were in their second to the fourth days of lactation, with a body condition score (BCS) of 2.0 to 2.2 (0 to 5 scale), a week before feed restricted diet. All animals were observed for uterine disease and did not present any signs of abnormality. Additionally, the health history of each animal was inspected and only those without any health problems, within the last 6 months before calving, were chosen.

At 24 ± 3 days in milk, sixteen multiparous Holstein cows were allowed *ad libitum* intake of a lactation diet CTRL, *n* = 8, 7.1 MJ/kg DM NEL, 17.4% Crude Protein. Their diet constituted of corn (24.2% dry matter), corn silage (29%), grass silage (25.5%), soybean meal (16.9%), and complemented with vitamins and minerals (0.9%). The underfed (REST) group received a ration diluted with barley straw (48% DM) for 96 h (RES, *n* = 8; 5.16 MJ/kg DM NEL, 12.2% CP). Therefore, the ratio of forage to concentrate differed from 58.0/42.0 in control (CONT) group to 79.2/20.8 in REST group [25]. Dry matter intake, milk yield, energy balance, plasma insulin, glucose, non-esterified fatty acids (NEFA)
and BHB (β-hydroxybutyrate) concentrations did not differ between CONT and REST immediately
before underfeeding (21.8 kg/day, 39.0 kg/day, -5.6 MJ/day, 22 µIU/mL, 3.78, 0.415, and 0.66 mM,
respectively, at day -1), but were significantly altered in REST at 72 h of underfeeding
(Supplementary Table S1). Following 72 h of restriction or control diet, one healthy rear mammary
quarter was injected with 50 µg of lipopolysaccharide E. coli 0111:B4; (LPS-EB Ultrapure, InvivoGen,
San Diego, CA, USA) diluted in 10 mL of sterile saline (CDM Lavoisier, Paris, France), containing 0.5
mg/mL BSA cell culture grade, endotoxin free, A9576, (Sigma-Aldrich, St. Louis, MO, USA), using a
sterile disposable syringe fitted with a sterile teat cannula. Mammary biopsies were performed 24 h
after the LPS injection, as previously described [62], corresponding to 96 h of feed restriction or
control diet. Tissue samples were immediately frozen in liquid nitrogen and stored at -80 °C prior to
RNA and protein analyses.

Throughout the study, milk samples were collected at 4 consecutive milkings each week before the beginning of the restriction and just before the LPS challenge and analyzed for SCC. Only healthy cows were included in the study. Additionally cows were screened for mastitis, one week before and immediately before the LPS challenge, using the California Mastitis Test (Neodis, Rambouillet, France) for all quarters, and somatic cell counts of rear quarter milk samples (Galilait, Theix, 63122 Saint Genès-Champanelle, France), one week before and immediately prior to the LPS challenge. Only cows with SCC lower than 100,000 cells/mL, in a rear quarter, were included in the study. Indeed, cows were considered healthy if the quarter SCC was inferior to 100,000 cells/mL and were free of any other signs of health problems [25]. Additionally, foremilk samples were collected from the LPS challenged quarters, immediately before the morning milking that preceded the LPS injection (time 0), and at 4, 6, 10, and 24 h after LPS injection. These quarter milk samples were analyzed for IL-8, IL1-β, TNF-alpha, and CXCL3 using Elisa [25].

4.2. RNA Preparation and Analyses

RNA and protein extractions were performed from the same mammary biopsy samples *n* = 16 animals, (8 CONT and 8 REST). The total RNA was extracted from 50 mg of the mammary gland (MG) by using the mirVana miRNA Isolation Kit (Thermo Fisher Sciences, Waltham, MA USA). The concentration and purity of RNA were estimated by spectrophotometry NanodropTH, (ND-1000, NanoDrop Technologies LLC, Wilmington, DE, USA), and by using the Bioanalyzer 2100 (Agilent Technologies Inc., Santa Clara, CA, USA), respectively. Once these validation steps were completed, only 12 cows (6 RES and 6 CTR) were kept for gene expression analyses at mRNA level, which presented a good and uniform quality of the samples for a microarray experiment.

4.3. RT-PCR Analyses

Reverse transcription (RT) was performed on 2 µg of total RNA using the “High Capacity RNA to cDNA” kit and following the manufacturer’s recommendations (Applied Biosystems, Villebon Sur Yvette, France) in a final volume of 20 µL. In parallel, negative controls were performed without the matrices. The primers are described in Table 3. The genes, *UXT2*, *CLN3*, and *EIF3K* were used as housekeeping genes [63]. Real-time quantitative PCR was performed on the StepOnePlus™ PCR System (Applied Biosystems, Villebon Sur Yvette, France), using 5 µL of 50 fold-diluted single-stranded cDNA and the TFPower SYBRGreen PCR Master Mix, according to the manufacturer’s instructions (Applied Biosystems, Villebon Sur Yvette, France). Subsequent to an initial denaturing step (95 °C for 10 min), the PCR mixture was subjected to the following two-step cycle, which was repeated 40 times: Denaturing for 15 s at 95 °C and annealing and extension for 45 s at 60 or 62 °C. The results were expressed as a fold change of Ct values relative to the control using the ΔCt method [64]. The significance was determined using a t-test with *p* < 0.05 considered as significant.
<table><caption>**Table 3.** Primer sequences and annealing temperatures used in real-time reverse transcription-PCR assays as the size of the amplicons (in bp).</caption><thead><tr><th>Gene Symbol</th><th>Primers Pair</th><th>Amplicon Size (bp)</th><th>T °C annealing</th><th>Reference</th></tr></thead><tbody><tr><td>CCL5</td><td>AGC AGT TGT CTT TAT CAC CAG GA<br>TCC AAA GCG TTG ATG TAC TCT C</td><td>87</td><td>60</td><td>[65]</td></tr><tr><td>CD 36</td><td>ACA GAT GTG GCT TGA GCG TG<br>ACT GGG TCT GTG TTT TGC AGG</td><td>186</td><td>58</td><td>[63]</td></tr><tr><td>CSN2</td><td>CTC AAA CCC CTG TGG TGG TG<br>AAA GGC CTG GAT GGG CAT AT</td><td>332</td><td>60</td><td>[1]</td></tr><tr><td>IL1</td><td>GAA TGG AAA CCC TCT CTC CC<br>GCT GCA GCT ACA TTC TTC CC</td><td>104</td><td>62</td><td>this article</td></tr><tr><td>IL8</td><td>TGG GCC ACA CTG TGA AAA T<br>TCA TGG ATC TTG CTT CTC AGC</td><td>138</td><td>62</td><td>[66]</td></tr><tr><td>LAP</td><td>GAA ATT CTC AAA GCT GCC GTA<br>TCC TCC TGC AGC ATT TTA CTT</td><td>114</td><td>60</td><td>[67]</td></tr><tr><td>RBP4</td><td>CAA CGG TTA CTG TGA TGG<br>GAG GCT GAG TAA GGT TAA TG</td><td>98</td><td>60</td><td>this article</td></tr><tr><td>STAT3</td><td>GTC TAA CAA TGG CAG CCT CTC AGC<br>AAG AGT TTC TCC GCC AGC GTC</td><td>405</td><td>60</td><td>[68]</td></tr><tr><td>TAP</td><td>GCC AGC ATG AGG CTC CAT<br>AAC AGG TGC CAA TCT GT</td><td>166</td><td>60</td><td>[29]</td></tr><tr><td>INSIG1</td><td>CTA GCC TCG AAC TAA AGC CTG ACT<br>TTC CTG TCT CAC CAC ACT TCA TCT</td><td>101</td><td>59</td><td>[69]</td></tr><tr><td colspan="5"><b>Housekeeping genes</b></td></tr><tr><td>UXT</td><td>TGT GGC CCT TGG ATA TGG TT<br>GGT TGT CGC TGA GCT CTG TG</td><td>101</td><td>60</td><td>[63]</td></tr><tr><td>EIF3K</td><td>CCA GGC CCA CCA AGA AGA A<br>TTA TAC CTT CCA GGA GGT CCA TGT</td><td>125</td><td>60</td><td>[63]</td></tr><tr><td>CLN3</td><td>TTC TGA CTC CTT GGG ACA CA<br>CAA CCT GCC CAC CTA TCA GT</td><td>62</td><td>62</td><td>[63]</td></tr></tbody></table>

## 4.4. Microarray Analyses

Microarray analyses were performed on twelve animals (6 RES and 6 CTRL) using 100 ng of total RNA from each MG sample. The total RNA was amplified, fluorescently labeled, and hybridized to the bovine 4 × 44K microarray (Agilent Technologies, Inc. Santa Clara, CA, USA), and all the procedures described below were performed according to the manufacturer instructions (Agilent Technologies, Inc. Santa Clara, CA, USA). Briefly, for each hybridization, the total RNA was linearly amplified and labeled with Cy3 using the one-color Low Input Quick Amp Labeling Kit. Then, 1650 ng of Cy3-labeled cRNA was hybridized on the microarrays using the Gene Expression Hyb Kit. Hybridization was performed for 17 h at 65 °C in a rotating hybridization oven at 10 rpm. Following hybridization, all microarrays were washed and scanned using the Agilent Microarray Scanner G2565A. The resulting TIFF-images (Tagged Image File Format) were processed using Feature Extraction software Version 11 to obtain normalized data. Normalized with 75th percentile shift, the data were analyzed using GeneSpring software. The moderated *t*-test with Westfall-Young familywise error rate (FWER) correction was applied [70]. The differences were considered significant at an adjusted *p* < 0.05. The data were accessible through the GEO series accession number GSE114975. The classification and functional analyses of differentially expressed genes were performed using PANTHER [71] and confirmed using Pathway Studio® software (Elsevier, Nederland).

## 4.5. Protein Preparation and Analyses

The proteins were extracted by homogenizing 80 mg MG tissue (n = 16; 8 RES and 8 CTRL) in 2 mL lysis buffer (8.3 M urea, 2 M thiourea, 2% CHAPS, 1% DDT). Following homogenization, the samples were incubated for 5 min at room temperature and centrifuged at 10,000× g for 30 min at 8 °C. The protein concentrations were measured in supernatant with Quick Start Bradford protein
assay (BioRad, Marnes-La-Coquette, France), aliquoted and then stored at -20 °C, until further preparation. Sample supernatants were mixed with 1 volume of Laemmli buffer and heated at 60 °C for 5 min. Separation, by SDS-PAGE (12% acrylamide), was performed using a Mini-Protean II electrophoresis unit (BioRad, Marnes-La-Coquette, France) and 100 µg protein loaded per lane. To concentrate the samples, the gels were run at 80 V until the dye front reached the bottom of the concentration gel. Gels were stained overnight in Coomassie brilliant blue G-250. Excised lanes were reduced and alkyled before de-staining in 25 mM ammonium bicarbonate with acetonitrile (50/50 v/v). Following dehydration with 100% acetonitrile, gel pieces were dried in a Speed Vacuum and the samples were preserved at -20 °C until LC MS/MS analysis.

## 4.6. LC MS/MS Analysis

The proteins were hydrolyzed overnight at 37 °C, using 800 ng (80 µL) of sequence grade-modified trypsin (Promega, France) per band. Subsequent to extraction by 64 µL of acetonitrile 100% and sonication, the peptides were concentrated in a Speed Vacuum and volume was adjusted to 30 µL with an aqueous solution (99.9% H2O, 0.1% TFA). Peptide mixture (2.5 µL) was injected into the nano HPLC Ultimate 3000, (Thermo Fisher Scientific, Courtaboeuf, France) after a preliminary step of desalting and concentration in the pre-column 300 µm × 5 mm, (ThermoFisher, Courtaboeuf, France) for 6 min, and a second step of separation in an analytical C18 column 75 µm, 25 cm, (Pepmap Thermo Fisher Scientific, Courtaboeuf, France) with a 10-40% gradient (A: 0.1% FA in water, B: 0.1% FA in acetonitrile) at 450 nL/min. The eluate was electrosprayed through the CaptiveSpray ion source into the mass spectrometer QTOF impact II (Bruker, Wissembourg, France) operated in CID Data Dependent mode. Each MS analysis was succeeded by as many MS/MS analysis as possible within 3 sec.

## 4.7. Protein Identification and Label-Free Quantitation

The raw files were loaded, at the end of each LC-MS/MS analysis, into the Progenesis QI software Non-linear Dynamics, v 4.1 (Newcastle upon Tyne, UK) and label-free quantitation was performed using a proprietary workflow alignment, peak picking, normalization, design set up, quantitation, and protein identification.

Regarding protein identification (Supplementary File S2), the Mascot V.2.5, internally licensed version (www.matrixscience.com) was used with uniprot-ref_Bos taurus database 19.840 sequences (07/2015). The following parameters were considered for the searches: Peptide mass tolerance was set to 10 ppm; fragment mass tolerance was set to 0.05 Da and a maximum of two missed cleavages was allowed. Variable modifications were methionine oxidation (M), carbamidomethylation (C) of cysteine and Deamidated (NQ). Protein identification was considered valid, if at least two peptides with a statistically significant Mascot score were assigned, with a false discovery rate (FDR) less than 1%.

Concerning label-free quantitation, all unique validated peptides of an identified protein were included, and the total cumulative abundance was calculated by summing up the abundances of all unique peptides allocated to the respective protein. Statistical analysis was performed, using the "between subject design," and the *p*-values were calculated by an analysis of variance, using the normalized abundances across all runs. Differential proteins were conserved for interpretation if the peptides' individual abundances showed a good correlation with protein abundance. All differential proteins were inspected manually with these correlation criteria. To extract the maximum biological information of differentially expressed proteins, PANTHER [71], Pathway Studio® software (Elsevier, Nederland) and UniProt [72], were used.

# 5. Conclusion

Undernutrition affected multiple aspects of MG function, as demonstrated by modifications of milk secretion, and MG mRNA and protein expression. During this study, expression analyses were performed 24 h post-LPS challenge corresponding to the period of inflammation resolution. The
effects of undernutrition on studied candidate genes, known as major genes relating to the innate immune responses, were weak. Therefore, the transcriptomic and proteomic analyses pointed out modifications of energy metabolism (fatty acid and glucose), and protein metabolism (synthesis and post-translational modification), respectively, but relatively few genes involved in immune response were affected. Our nutrigenomic analyses have suggested that undernutrition of early lactating cows modified the mammary gland metabolism. The holistic analyses of the systemic reaction in the mammary gland expands the knowledge of the effects of NEB and metabolic imbalance occurring in early lactation, during inflammation. These identified genes may be relevant for quantitative trait loci studies and genomic selection.

**Supplementary Materials:** Supplementary materials can be found at www.mdpi.com/xxx/s1. Table S1: Plasma insulin and metabolite concentration at the day of biopsy after dietary treatment and response to LPS challenge. Occurring at day 24 ± 3 of lactation, animals were assigned to a control (CONT, *n* = 8) or restricted (REST, *n* = 8) group. The REST animals received the ration diluted with barley straw (48% DM) for 4 days when cows from CONT were allowed to continue ad libitum intake of a lactation diet (7.1 MJ/kg DM NEL, 17.4% CP). Occurring on day 3, corresponding to the 27th day of lactation, the rear mammary quarter of animals from both groups was injected with 50 µg of LPS. Mammary biopsies were performed 24 h after LPS challenge. *p* < 0.01 for all variables. File S2: Protein analysis report.

**Author contributions:** Conceptualization, J.P., C.L., P.G., and K.P.; funding acquisition, P.G., C.L., and J.P.; transcriptomic analysis, K.P., C.L., Y.F., and C.B.; proteomic analysis, C.C., and K.P.; Bioinformatic analysis, K.P.; supervision, C.L., and J.P.; writing—original draft, K.P.; writing—review, C.L., J.P., Y.F., C.B., C.C., and P.G.

**Funding:** Financial support for this research was provided by GISA meta-program of INRA (Ruminflame project).

**Acknowledgments:** The authors thank the staff at Herbipole INRA (UE1414, Theix, France), Simon Collange for performing biopsies, Sebastien Bes, Arnaud Delavaud, Caroline Soulard, Martine Tourret, Didier Bany, Emile Tixier, Cyril Labonne, and Jacques Rouel for technical assistance, and for sample collection and laboratory analyses. The authors are grateful to the Galilait laboratory (Clermont-Ferrand, France) for the milk component, SCC and microbiological analyses. Special appreciation is extended to Gilles Foucras and Pascal Rainard for their helpful discussion regarding LPS and choice of candidate genes. Transcriptomic analyses were performed at the "PlateForme d'Exploration du Métabolisme" (PFEM, INRA, Theix, France).

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

## References

1. Ollier, S.; Robert-Granié, C.; Bernard, L.; Chilliard, Y.; Leroux, C. Mammary transcriptome analysis of food-deprived lactating goats highlights genes involved in milk secretion and programmed cell death. *J. Nutr.* **2007**, *137*, 560–567, doi:10.1093/jn/137.3.560.

2. Minuti, A.; Zhou, Z.; Graugnard, D.E.; Rodriguez-Zas, S.L.; Palladino, A.R.; Cardoso, F.C.; Trevisi, E.; Loor, J.J. Acute mammary and liver transcriptome responses after an intramammary Escherichia coli lipopolysaccharide challenge in postpartal dairy cows. *Physiol. Rep.* **2015**, *3*, e12388, doi:10.14814/phy2.12388.

3. Seegers, H.; Fourichon, C.; Beauveau, F. Production effects related to mastitis and mastitis economics in dairy cattle herds. *Vet. Res.* **2003**, *34*, 475–491, doi:10.1051/vetres:2003027.

4. Matthews, K.R.; Almeida, R.A.; Oliver, S.P. Bovine mammary epithelial cell invasion by Streptococcus uberis. *Infect. Immun.* **1994**, *62*, 5641–5646.

5. Almeida, R.A.; Matthews, K.R.; Cifrian, E.; Guidry, A.J.; Oliver, S.P. Staphylococcus aureus invasion of bovine mammary epithelial cells. *J. Dairy Sci.* **1996**, *70*, 1021–1026, doi:10.3168/jds.S0022-0302(96)76454-8.

6. Lam, T.J.; Lipman, L.J.; Schukken, Y.H.; Gaastra, W.; Brand, A. Epidemiological characteristics of bovine clinical mastitis caused by Staphylococcus aureus and Escherichia coli studied by DNA fingerprinting. *Am. J. Vet. Res.* **1996**, *571*, 39–42.

7. Mallard, B.A.; Dekkers, J.C.; Ireland, M.J.; Leslie, K.E.; Sharif, S.; Lacey Vankampen, C.; Wagter, L.; Wilkie, B.N. Alteration in immune responsiveness during the periparturient period and its ramification on dairy cow and calf health. *J. Dairy Sci.* **1998**, *81*, 585–595, doi:3168/jds.S0022-0302(98)75612-7.

8. Sordillo, L.M.; Raphael, W. Significance of Metabolic Stress, Lipid Mobilization, and Inflammation on Transition Cow Disorders. *Vet. Clin. Food Anim.* **2013**, *29*, 267–278, doi:10.1016/j.cvfa.2013.03.002.
9. Suriyasathaporn, W.; Heuer, C.; Noordhuizen-Stassen, E.N.; Schukken, Y.H. Hyperketonemia and the impairment of udder defense: A review. *Vet. Res.* **2000**, *31*, 397-412, doi:10.1051/vetres:2000128.

10. Shuster, D.E.; Lee, E.; Kehrli, M.E. Bacterial growth, inflammatory cytokine production, and neutrophil recruitment during coliform mastitis in cows within ten days after calving, compared with cows at midlactation. *Am. J. Vet. Res.* **1996**, *57*, 1569-1575.

11. Bruckmaier, R.M. Gene expression of factors related to the immune reaction in response to intramammary Escherichia coli lipopolysaccharide challenge. *J. Dairy Res.* **2005**, *72*, 120-124, doi:10.1017/S0022029905001159.

12. Moyes, K.M.; Drackley, J.K.; Salak-Johnson, J.L.; Morin, D.E.; Hope, J.C.; Loor, J.J. Dietary-induced negative energy balance has minimal effects on innate immunity during a Streptococcus uberis mastitis challenge in dairy cows during midlactation. *J. Dairy Sci.* **2009**, *92*, 4301-4316, doi:10.3168/jds.2009-2170.

13. Perkins, K.H.; VandeHaar, M.J.; Burton, J.L.; Liesman, J.S.; Erskine, R.J.; Elsasser, T.H. Clinical Responses to Intramammary Endotoxin Infusion in Dairy Cows Subjected to Feed Restriction. *J. Dairy Sci.* **2002**, *85*, 1724-1731, doi:3168/jds.S0022-0302(02)74246-X.

14. Sordillo, L.M.; Shafer-Weaver, K.; DeRosa, D. Immunobiology of the mammary gland. *J. Dairy Sci.* **1997**, *80*, 1851-1865, doi:10.3168/jds.S0022-0302(97)76121-6.

15. Waldron, M.R.; Kulick, A.E.; Bell, A.W.; Overton, T.R. Acute Experimental Mastitis Is Not Causal Toward the Development of Energy-Related Metabolic Disorders in Early Postpartum Dairy Cows. *J. Dairy Sci.* **2006**, *89*, 596-610, doi:10.3168/jds.S0022-0302(06)72123-3.

16. Rinaldi, M.; Li, R.W.; Capuco, A.V. Mastitis associated transcriptomic disruptions in cattle. *Vet. Immunol. Immunopathol.* **2010**, *138*, 267-279, doi:10.1016/j.vetimm.2010.10.005.

17. Loor, J.J.; Moyes, K.M.; Bionaz, M. Functional adaptations of the transcriptome to mastitis-causing pathogens: The mammary gland and beyond. *J. Mammary Gland Biol. Neoplasia* **2011**, *16*, 305-322, doi:10.1007/s10911-011-9232-2.

18. Buitenhuis, B.; Røntved, C.M.; Edwards, S.M.; Ingvartsen, K.L.; Sørensen, P. In depth analysis of genes and pathways of the mammary gland involved in the pathogenesis of bovine Escherichia coli-mastitis. *BMC Genom.* **2011**, *12*, 130, doi:10.1186/1471-2164-12-130.

19. Mitterhuemer, S.; Petzl, W.; Krebs, S.; Mehne, D.; Klanner, A.; Wolf, E.; Zerbe, H.; Blum, H. Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland. *BMC Genom.* **2010**, *11*, 138, doi:10.1186/1471-2164-11-138.

20. Günther, J.; Koczan, D.; Yang, W.; Nurnberg, G.; Repsilber, D.; Schuberth, H.J.; Park, Z.; Maqbool, N.; Molenaar, A.; Seyfert, H.M. Assessment of the immune capacity of mammary epithelial cells: Comparison with mammary tissue after challenge with Escherichia coli. *Vet. Res.* **2009**, *40*, 31, doi:10.1051/vetres/2009014.

21. Günther, J.; Esch, K.; Poschadel, N.; Petzl, W.; Zerbe, H.; Mitterhuemer, S.; Blum, H.; Seyfert, H.M. Comparative kinetics of Escherichia coli- and Staphylococcus aureus-specific activation of key immune pathways in mammary epithelial cells demonstrates that S. aureus elicits a delayed response dominated by interleukin-6 IL-6) but not by IL-1A or tumor necrosis factor alpha. *Infect. Immun.* **2011**, *79*, 695-707, doi:10.1128/IAI.01071-10.

22. Griesbeck-Zilch, B.; Meyer, H.H.D.; Kuhn, C.; Schwerin, M.; Wellnitz, O. Staphylococcus aureus and Escherichia coli cause deviating expression profiles of cytokines and lactoferrin messenger ribonucleic acid in mammary epithelial cells. *J. Dairy Sci.* **2008**, *91*, 2215-2224, doi:10.3168/jds.2007-0752.

23. Hogarth, C.J.; Fitzpatrick, J.L.; Nolan, A.M.; Young, F.J.; Pitt, A.; Eckersall, P.D. Differential protein composition of bovine whey: A comparison of whey from healthy animals and from those with clinical mastitis. *Proteomics* **2004**, *4*, 2094-2100, doi:10.1002/pmic.200300723.

24. Boehmer, J.L.; Bannerman, D.D.; Shefcheck, K.; Ward, J.L. Proteomic analysis of differentially expressed proteins in bovine milk during experimentally induced Escherichia colimastitis. *J. Dairy Sci.* **2008**, *91*, 4206-4218.

25. Pires, J.; Pawlowski, K.; Rouel, J.; Delavaud, C.; Foucras, G.; Rainard, P.; Germon, P.; Leroux, C. Undernutrition modifies metabolic responses to intramammary lipopolysaccharide but has limited effects on inflammation indicators in early lactation cows. *J. Dairy Sci.* **2019**, doi:3168/jds.2018-15446.

26. Rainard, P.; Riotlet, C. Innate immunity of the bovine mammary gland. *Vet. Res.* **2006**, *37*, 369-400.

27. Atabai, K.; Sheppard, D.; Werb, Z. Roles of the Innate Immune System in Mammary Gland Remodeling During Involution. *J. Mammary Gland Biol. Neoplasia* **2007**, *121*, 37-45.
28. Hughes, K.; Watson, C.J. The spectrum of STAT functions in mammary gland development. *Jak-stat* **2012**, *13*, 151–158.

29. López-Meza, J.E.; Gutiérrez-Barroso, A.; Ochoa-Zarzosa, A. Expression of tracheal antimicrobial peptide in bovine mammary epithelial cells. *Res. Vet. Sci.* **2009**, *87*, 59–63, doi:10.1016/j.rvsc.2008.12.005.

30. Yarus, S.; Rosen, J.M.; Cole, A.M.; Diamond, G. Production of active bovine tracheal antimicrobial peptide in milk of transgenic mice. *Proc. Natl. Acad. Sci. USA* **1996**, *93*,14118–14121.

31. Moyes, K.M.; Sørensen, P.; Bionaz, M. The Impact of Intramammary Escherichia coli Challenge on Liver and Mammary Transcriptome and Cross-Talk in Dairy Cows during Early Lactation Using RNAseq. *PLoS ONE* **2016**, *11*, e0157480, doi:1371/journal.pone.0157480.

32. Priestman, D.A.; Mistry, S.C.; Halsall, A.; Randle, P.J. Role of protein synthesis and of fatty acid metabolism in the longer-term regulation of pyruvate dehydrogenase kinase. *Biochem. J.* **1994**, *300*, 659–664.

33. Ibeagha-Awemu, E.M.; Li, R.; Ammah, A.A.; Dudemaine, P.L.; Bissonnette, N.; Benchaar, C.; Zhao, X. Transcriptome adaptation of the bovine mammary gland to diets rich in unsaturated fatty acidFAs shows greater impact of linseed oil over safflower oil on gene expression and metabolic pathways. *BMC Genom.* **2016**, *17*, 104, doi:10.1186/s12864-016-2423-x.

34. Bouvier-Muller, J.; Allain, C.; Tabouret, G.; Enjalbert, F.; Portes, D.; Noirot, C.; Rupp, R.; Foucras, G. Whole blood transcriptome analysis reveals potential competition in metabolic pathways between negative energy balance and response to inflammatory challenge. *Sci. Rep.* **2017**, *7*, 2379, doi:10.1038/s41598-017-02391-y.

35. Lee, J.H.; Kim, E.J.; Kim, D.K.; Lee, J.M.; Park, S.B.; Lee, I.K.; Harris, R.A.; Lee, M.O.; Choi, H.S. Hypoxia Induces PDK4 Gene Expression through Induction of the Orphan Nuclear Receptor ERRc. *PLoS ONE* **2012**, *7*, e46324, doi:10.1371/journal.pone.0046324.

36. Rasmussen, B.B.; Wolfe, R.R. Regulation of fatty acid oxidation in skeletal muscle. *Annu Rev. Nutr.* **1999**, *19*, 463–484, doi:10.1146/annurev.nutr.19.1.463.

37. Akbar, H.; Bionaz, M.; Carlson, D.B.; Rodriguez-Zas, S.L.; Everts, R.E.; Lewin, H.A.; Drackley, J.K.; Loor, J.J. Feed restriction.; but not L-carnitine infusion.; alters the liver transcriptome by inhibiting sterol synthesis and mitochondrial oxidative phosphorylation and increasing gluconeogenesis in mid-lactation dairy cows. *J. Dairy Sci.* **2013**, *96*, 2201–2213, doi:10.3168/jds.2012-6036.

38. Drackley, J.K. Biology of dairy cows during the transition period: The final frontier? *J. Dairy Sci.* **1999**, *82*, 2259–2273, doi:3168/jds.S0022-0302(99)75474-3.

39. Palmieri, F. The mitochondrial transporter family SLC25: Identification, properties and physiopathology. *Mol. Asp. Med.* **2013**, *34*, 465–484, doi:10.1016/j.mam.2012.05.005.

40. Moyes, K.M.; Drackley, J.D.; Morin, D.E.; Rodriguez-Zas, S.L.; Everts, R.E.; Lewin, H.A.; Loor, J.J. Predisposition of cows to mastitis in non-infected mammary glands: Effects of dietary-induced negative energy balance during mid-lactation on immune-related genes. *Funct. Integr. Genom.* **2011**, *11*, 151–156, doi:10.1007/s10142-010-0186-z.

41. Kashyap, D.R.; Rompca, A.; Gaballa, A.; Helmann, J.D.; Chan, J.; Chang, C.J.; Hozo, I.; Gupta, D.; Dziarski, R. Peptidoglycan Recognition Proteins Kill Bacteria by Inducing Oxidative.; Thiol.; and Metal Stress. *PLoS Pathog.* **2014**, *10*, doi:1371/journal.ppat.1004280.

42. Ohtsuka, H.; Koiwa, M.; Fukuda, S.; Satoh, Y.; Hayashi, T.; Hoshi, F.; Yoshino, T.; Kawamura, S. Changes in peripheral leukocyte subsets in dairy cows with inflammatory diseases after calving. *J. Vet. Med. Sci.* **2004**, *66*, 905–909, doi:1292/jvms.66.905.

43. Li, H.; Hou, S.; Wu, X.; Nandagopal, S.; Lin, F.; Kung, S.; Marshall, A.J. The Tandem PH Domain-Containing Protein 2 (TAPP2) Regulates Chemokine-Induced Cytoskeletal; Reorganization and Malignant B Cell Migration. *PLoS ONE* **2013**, *8*, e57809, doi:10.1371/journal.pone.0057809.

44. Lee, J.W.; Paape, M.J.; Elsasser, T.H.; Zhao, X. Elevated milk soluble CD14 in bovine mammary glands challenged with Escherichia coli lipopolysaccharide. *J. Dairy Sci.* **2003**, *86*, 2382–2389, doi:10.3168/jds.S0022-0302(03)73832-6.

45. Bannerman, D.D.; Paape, M.J.; Hare, W.R.; Sohn, E.J. Increased levels of LPS-binding protein in bovine blood and milk following bacterial lipopolysaccharide challenge. *J. Dairy Sci.* **2003**, *86*, 3128–3137, doi:10.3168/jds.S0022-0302(03)73914-9.

46. Han, S.P.; Tang, Y.H.; Smith, R. Functional diversity of the hnRNPs: Past.; present and perspectives. *Biochem. J.* **2010**, *430*, 379–392, doi:10.1042/BJ20100396.
47. Geuens, T.; Bouhy, D.; Timmerman, V. The hnRNP family: Insights into their role in health and disease. *Hum. Genet.* **2016**, *135*, 851–867, doi:10.1007/s00439-016-1683-5.

48. Raffetseder, U.; Liehn, E.A.; Weber, C.; Mertens, P.R. Role of cold shock Y-box protein-1 in inflammation; atherosclerosis and organ transplant rejection. *Eur. J. Cell Biol.* **2012**, *91*, 567–575, doi:10.1016/j.ejcb.2011.07.001.

49. Apostolidis, S.A.; Rodríguez-Rodríguez, N.; Suárez-Fueyo, A.; Dioufa, N.; Ozcan, E.; Crispín, J.C.; Tsokos, M.G.; Tsokos, G.C. Phosphatase PP2A is requisite for the function of regulatory T cells. *Nat. Immunol.* **2016**, *17*, 556–564, doi:10.1038/ni.3390.

50. Bhavsar, R.B.; Makley, L.N.; Tsonis, P.A. The other lives of ribosomal proteins. *Hum. Genom.* **2010**, *4*, 327–344, doi:10.1186/1479-7364-4-5-327.

51. Schachter, H.; Freeze, H.H. Glycosylation diseases: Quo vadis? *Biochim. Biophys. Acta* **2009**, *1792*, 925–930, doi:1016/j.bbadis.2008.11.002.

52. Palorini, R.; Cammarata, F.P.; Balestrieri, C.; Monestiroli, A.; Vasso, M.; Gelfi, C.; Alberghina, L.; Chiaradonna, F. Glucose starvation induces cell death in K-ras-transformed cells by interfering with the hexosamine biosynthesis pathway and activating the unfolded protein response. *Cell Death Dis.* **2013**, *4*, 732, doi:10.1038/cddis.2013.257.

53. Alanen, H.I.; Salo, K.E.; Pirneskoski, A.; Ruddock, L.W. pH dependence of the peptide thiol-disulfide oxidase activity of six members of the human protein disulfide isomerase family. *Antioxid. Redox Signal.* **2006**, *83*, 283–291, doi:10.1089/ars.2006.8.283.

54. Galligan, J.J.; Petersen, D.R. The human protein disulfide isomerase gene family. *Hum. Genom.* **2012**, *6*, 6, doi:10.1186/1479-7364-6-6.

55. Grootjans, J.; Kaser, A.; Kaufman, R.J.; Blumberg, R.S. The unfolded protein response in immunity and inflammation. *Nat. Rev. Immunol.* **2016**, *16*, 469–484, doi:10.1038/nri.2016.62.

56. Garbi, N.; Tanaka, S.; Momburg, F.; Hammerling, G.J. Impaired assembly of the major histocompatibility complex class I peptide-loading complex in mice deficient in the oxidoreductase ERp57. *Nat. Immunol.* **2006**, *7*, 93–102, doi:10.1038/ni1288.

57. Lilie, H.; McLaughlin, S.; Freedman, R.; Buchner, J. Influence of protein disulfide isomerase PDI) on antibody folding in vitro. *J. Biol. Chem.* **1994**, *269*, 14290–14296.

58. Bregier, C.; Kupikowska, B.; Fabczak, H.; Fabczak, S. CCT chaperonins and their cochaperons. *Postepy Biochem.* **2008**, *54*, 64–70.

59. Carneiro, P.; Peters, L.C.; Vorraro, F.; Borrego, A.; Ribeiro, O.G.; Starobinas, N.; Jensen, J.R.; Cabrera, W.H.K.; Ibanez, O.M.; De Franco, M. Gene expression profiles of bone marrow cells from mice phenotype-selected for maximal or minimal acute inflammations: Searching for genes in acute inflammation modifier loci. *Immunology* **2009**, *128*, 562–571, doi:10.1111/j.1365-2567.2008.03032.x.

60. Horvath, A.J.; Irving, J.A.; Rossjohn, J.; Law, R.H.; Bottomley, S.P.; Quinsey, N.S.; Pike, R.N.; Coughlin, P.B.; Whisstock, J.C. The murine orthologue of human antichymotrypsin: A structural paradigm for clade A3 serpins. *J. Biol. Chem.* **2005**, *280*, 43168–43178, doi:10.1074/jbc.M505598200.

61. Nomura, K.; Kanegane, H.; Karasuyama, H.; Tsukada, S.; Agematsu, K.; Murakami, G.; Sakazume, S.; Sako, M.; Tanaka, R.; Kuniya, Y.; et al. Genetic defect in human X-linked agammaglobulinemia impedes a maturational evolution of pro-B cells into a later stage of pre-B cells in the B-cell differentiation pathway. *Blood* **2000**, *96*, 610–617.

62. Bernard, L.; Richard, C.; Gelin, V.; Leroux, C.; Heyman, Y. Milk fatty acid composition and mammary lipogenic genes expression in bovine cloned and control cattle. *Livest. Sci.* **2015**, *176*, 188–195, DOI10.1016/j.livsci.2015.03.016.

63. Bonnet, M.; Bernard, L.; Bes, S.; Leroux, C. Selection of reference genes for quantitative real-time PCR normalisation in adipose tissue.; muscle.; liver and MG from ruminants. *Animal* **2013**, *7*, 1344–1353, doi:10.1017/S1751731113000475.

64. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta C(T)) Method. *Methods* **2001**, *25*, 402–408, doi:10.1006/meth.2001.1262.

65. Magee, P., Pearson, S., Whittingham-Dowd, J., Allen, J. PPARgamma as a molecular target of EPA anti-inflammatory activity during TNF-alpha-impaired skeletal muscle cell differentiation. *The Journal of nutritional biochemistry* **2012**, *23*(11):1440-1448.

66. Zaros, L.G., Bricarello P.A., Amarante, A.F.T., Coutinho L.L. Quantification of bovine cytokine expression using real-time RT-PCR methodology. *Genet. Mol. Biol.* **2007**, *30* (3), 575-579.
67. Swanson, K., Gorodetsky, S., Good, L., Davis, S., Musgrave, D., Stelwagen, K., Farr, V., Molenaar, A. Expression of a beta-defensin mRNA, lingual antimicrobial peptide, in bovine mammary epithelial tissue is induced by mastitis. *Infection and immunity* **2004**, *72*(12):7311-7314.

68. Wang, S. W., Wang, S. S., Wu, D. C., Lin, Y. C., Ku, C. C., Wu, C. C., Chai, C. Y., Lee, J. N., Tsai, E. M., Lin, C. L., Yang, R. C., Ko, Y. C., Yu, H. S., Huo, C., Chuu, C. P., Murayama, Y., Nakamura, Y., Hashimoto, S., Matsushima, K., Jin, C., Eckner, R., Lin, C. S., Saito, S., Yokoyama, K. K.. Androgen receptor-mediated apoptosis in bovine testicular induced pluripotent stem cells in response to phthalate esters. *Cell death & disease* **2013**, *4*:e907. doi:10.1038/cddis.2013.420.

69. Leroux, C., Bernard, L., Faulconnier, Y., Rouel, J., de la Foye, A., Domagalski, J., Chilliard, Y. Bovine Mammary Nutrigenomics and Changes in the Milk Composition due to Rapeseed or Sunflower Oil Supplementation of High-Forage or High-Concentrate Diets. *Journal of nutrigenetics and nutrigenomics* 2016, *9*(2-4):65-82.

70. Westfall, P.H.; Zaykin, D.V.; Young, S.S. Multiple tests for genetic effects in association studies. In *Biostatistical Methods (Methods in Molecular Biology)*; Looney, S., Ed.; Humana Press: Toloway, NJ, USA, 2001; pp. 143–168.

71. Thomas, P.D.; Campbell, M.J.; Kejariwal, A.; Mi, H.; Karlak, B.; Daverman, R.; Diemer, K.; Muruganujan, A.; Narechania, A. PANTHER: A library of protein families and subfamilies indexed by function. *Genome Res.* **2003**, *13*, 2129–2214, doi:10.1101/gr.772403.

72. The UniProt Consortium. UniProt: The universal protein knowledgebase. *Nucleic Acids Res.* **2017**, *45*, D158–D169, doi:10.1093/nar/gkw1099.

<figure><img src="image_10.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/).