# Whole blood transcriptomics is relevant to identify molecular changes in response to genetic selection for feed efficiency and nutritional status in the pig

Maëva Jégou, Florence Gondret, Annie Vincent, Christine Trefeu, Hélène

Gilbert, Isabelle Louveau

► To cite this version:

Maëva Jégou, Florence Gondret, Annie Vincent, Christine Trefeu, Hélène Gilbert, et al. Whole blood transcriptomics is relevant to identify molecular changes in response to genetic selection for feed efficiency and nutritional status in the pig. PLoS ONE, 2016, 11 (1), pp.e0146550. 10.1371/journal.pone.0146550. hal-02639563

HAL Id: hal-02639563

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

Submitted on 28 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>

RESEARCH ARTICLE

# Whole Blood Transcriptomics Is Relevant to Identify Molecular Changes in Response to Genetic Selection for Feed Efficiency and Nutritional Status in the Pig

Maëva Jégou¹,², Florence Gondret¹,², Annie Vincent¹,², Christine Tréfeu¹,², Hélène Gilbert³,⁴,⁵, Isabelle Louveau¹,²*

1 INRA, UMR1348 Pegase, Saint-Gilles, France, 2 Agrocampus Ouest, UMR1348 Pegase, Rennes, France, 3 INRA, UMR1388 GenPhySE, Castanet-Tolosan, France, 4 Université de Toulouse, INP, ENSAT, UMR1388 GenPhySE, Castanet-Tolosan, France, 5 Université de Toulouse, INP, ENVT, UMR1388 GenPhySE, Toulouse, France

* isabelle.louveau@rennes.inra.fr

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

OPEN ACCESS

**Citation:** Jégou M, Gondret F, Vincent A, Tréfeu C, Gilbert H, Louveau I (2016) Whole Blood Transcriptomics Is Relevant to Identify Molecular Changes in Response to Genetic Selection for Feed Efficiency and Nutritional Status in the Pig. PLoS ONE 11(1): e0146550. doi:10.1371/journal.pone.0146550

**Editor:** Ramona Natacha PENA i SUBIRÀ, University of Lleida, SPAIN

**Received:** September 18, 2015

**Accepted:** December 19, 2015

**Published:** January 11, 2016

**Copyright:** © 2016 Jégou et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All microarray data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Subserie accession number GSE70838 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70838).

**Funding:** Maëva Jégou is supported by a PhD grant from INRA and Région Bretagne (France). This study was funded by the French National Research Agency (Agence Nationale de la Recherche ANR, ANR-11-BSV7-0004 FatInteger).

## Abstract

The molecular mechanisms underlying feed efficiency need to be better understood to improve animal efficiency, a research priority to support a competitive and sustainable live-stock production. This study was undertaken to determine whether pig blood transcriptome was affected by differences in feed efficiency and by ingested nutrients. Growing pigs from two lines divergently selected for residual feed intake (RFI) and fed isoproteic and isocaloric diets contrasted in energy source and nutrients were considered. Between 74 and 132 days of age, pigs (n = 12 by diet and by line) received a regular diet rich in cereals and low in fat (LF) or a diet where cereals where partially substituted by lipids and fibers (HF). At the end of the feeding trial, the total number of white blood cells was not affected by the line or by the diet, whereas the red blood cell number was higher (P<0.001) in low RFI than in high RFI pigs. Analysis of the whole blood transcriptome using a porcine microarray reveals a higher number of probes differentially expressed (DE) between RFI lines than between diets (2,154 versus 92 probes DE, P<0.01). This corresponds to 528 overexpressed genes and 477 underexpressed genes in low RFI pigs compared with high RFI pigs, respectively. Overexpressed genes were predominantly associated with translational elongation. Under-expressed genes were mainly involved in the immune response, regulation of inflammatory response, anti-apoptosis process, and cell organization. These findings suggest that selec-tion for RFI has affected the immune status and defense mechanisms of pigs. Genes DE between diets were mainly related to the immune system and lipid metabolism. Altogether, this study demonstrates the usefulness of the blood transcriptome to identify the main bio-logical processes affected by genetic selection and feeding strategies.
# Introduction

Production of efficient animals is an important issue for the livestock industry [1, 2] because this implies the reduction of the amount of feed resources needed to produce meat and contributes to reduce environmental wastes and emissions. Feed efficiency can be improved by breeding and feeding strategies. In recent years, residual feed intake (RFI) defined as the difference between the observed feed intake and the feed intake predicted from growth and maintenance requirements, has been studied as a measure of net feed efficiency in selection experiments [3]. Global performance differs between pigs selected for RFI, involving metabolism differences [4, 5] linked with gene expression variations in muscle [6]. The use of alternative feed resources (generally rich in fibers) that do not compete with food for humans is a challenge for the present and the future feeding strategies. However, interactions between RFI genotypes and these resources can affect body composition [7, 8]. An integrated approach combining genetics and nutrition is then necessary to get a better understanding of the complex biology underlying RFI and associated responses [9].

Peripheral blood is an accessible source of information. Indeed, the collection of blood samples is relatively easy compared to other tissues and does not alter the animal physiological status (anesthesia is unnecessary). It may further allow the investigation of kinetics of changes in different blood parameters. Advances in high-throughput technologies such as transcriptomics offer opportunities to answer complex biological questions. To date, genes expressed in peripheral blood cells have been shown to reflect physiological and pathological events occurring in different tissues [10, 11]. It has also been proposed that gene expression profiles from the whole blood or from peripheral blood mononuclear cells (PBMC) can highlight biological processes related to the regulation of body composition in human [12] and rodent [13]. Therefore, the transcriptome analysis of whole blood is a relevant approach to better understand the molecular mechanisms underlying differences in RFI, to determine the possible relationships between blood traits and production traits, and more generally to identify easily accessible targets to monitor physiological changes in response to factors such as nutrition. The current study was undertaken to determine whether diets with contrasted nutrients and genetic selection for RFI affect the whole blood transcriptome profiles of growing pigs.

# Materials and Methods

## Ethics statement

The care and use of pigs were performed in compliance with the European Union legislation (directive 2010/63/EU). The current protocol was approved by the local Ethics Committee in Animal Experiment of Rennes, France (Comité Rennais d'Ethique en matière d'Expérimentation Animale, CREEA, http://ethique.ipbs.fr/creeapresent.html; agreement N°07-2012). All animals were reared and killed in compliance with national regulations and according to procedures approved by the French veterinary Services at INRA Pegase facilities.

## Animals, diets and slaughtering

The pig lines and feeding trials used in the current study have been previously described in details [14]. Briefly, a total of 48 purebred French Large White castrated male pigs in the course of a divergent selection experiment for RFI (8th generation of selection; n = 24 per line) were recruited. From 74 ± 3 days of age and an average body weight (BW) of 22.6 ± 0.5 kg, pigs were housed in full slatted-floor isolated pens. Within each line, they were randomly assigned to one of the two dietary groups (n = 12 per diet), and fed ad libitum either a low-fat, low-fiber (LF) diet or a high-fat, high-fiber (HF) diet. The two diets were mainly based on cereals (wheat and
barley) and soybean meal. The HF diet was formulated by the partial replacing of cereals in the LF diet by wheat straw (11.5%) and a mixture of rapeseed and soya oils (7.5%). Each diet provided the same crude protein and metabolizable energy contents (Table 1). Growing and finishing formula were successively distributed to pigs during their growth (transition around 112 days of age). All pigs were killed at the same age (132.0 ± 0.5 days of age; average BW of 75.6 ± 1.1 kg). The killing procedure included electronarcosis and jugular exsanguination, and was performed 2 h after the first morning meal to obtain animals in a post-prandial state, in the INRA experimental slaughterhouse (Saint-Gilles, France).

Final blood samples were taken at exsanguination at the end of the feeding trial (132 days of age). Blood was collected into EDTA Venosafe tubes (Laboratoires Terumo, Guyancourt, France). For samples dedicated to RNA extraction, one volume of blood sample was mixed with one volume of lysis buffer from the Nucleospin 8 RNA blood kit (Macherey-Nagel, Lyon, France). The obtained mixture was then stored at -70°C for later analyses.

## Blood cell count

Whole blood cells can be separated into three categories. The first category refers to white blood cells (WBC), represented mainly by lymphocytes, monocytes and granulocytes. The two other categories refer to red blood cells (RBC) and platelet cells. Blood cell counts from those three categories were measured on whole blood samples taken at the end of the feeding trial using a hematology automatic cell counter calibrated for pigs (MS9-3, MELET SCHLOESING Laboratoires, Osny, France).

## RNA extraction

Total RNA were extracted from whole blood samples at both ages using a commercial kit (Nucleospin blood kit, Macherey Nagel, Hoerdt, France) according to the manufacturer's instructions. Then, residual genomic DNA was removed from RNA samples by a DNase treatment (DNA-free kit, Applied Biosystems, Foster City, CA, USA) in the presence of a RNase inhibitor (Thermo Scientific, Illkirch, France). Extracted RNA samples were quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Illkirch, France). The integrity of isolated RNA was assessed using the Agilent RNA 6000 Nano kit with an Agilent 2100

<table><caption>**Table 1.** Composition of low-fat, low-fiber (LF) or high-fat, high-fiber (HF) diets given to pigs during the growing and finishing periods.</caption><thead><tr><th rowspan="2">Composition<sup>a</sup></th><th colspan="2">Growing formula<sup>b</sup></th><th colspan="2">Finishing formula<sup>b</sup></th></tr><tr><th>LF</th><th>HF</th><th>LF</th><th>HF</th></tr></thead><tbody><tr><td>Fat</td><td>2.1</td><td>7.0</td><td>2.2</td><td>7.4</td></tr><tr><td>Starch</td><td>42.6</td><td>30.8</td><td>48.5</td><td>36.2</td></tr><tr><td>Crude protein</td><td>17.4</td><td>17.3</td><td>13.4</td><td>13.1</td></tr><tr><td>Neutral detergent fibers</td><td>12.8</td><td>17.7</td><td>11.5</td><td>18.3</td></tr><tr><td>Acid detergent fibers</td><td>3.8</td><td>8.5</td><td>3.5</td><td>7.9</td></tr><tr><td>Metabolizable energy</td><td>12.9</td><td>12.9</td><td>12.9</td><td>12.9</td></tr></tbody></table>

<footnote><sup>a</sup>Details on diet composition can be found in Gondret et al. [14]. Nutrient content is expressed in g per 100 g of feed (as fed basis). Metabolizable energy content is expressed in MJ per kg of feed.</footnote>

<footnote><sup>b</sup>A growing formula was distributed to pigs during 6 weeks from 76 days of age, while the finishing formula was provided thereafter and until slaughter.</footnote>

<footnote>Initial blood samples were collected from the jugular vein at 74 days of age on living pigs.</footnote>

<footnote>doi:10.1371/journal.pone.0146550.t001</footnote>
Bioanalyzer (Agilent Technologies France, Massy, France). All samples met quality criteria. Ratios of A260/280 and A260/230 were greater than 1.8. Average RNA integrity number was of 9.2 with values ranging from 7.9 to 9.7.

RNA labeling and microarray hybridization

Transcriptomics analyses were performed using a custom porcine microarray (8x60K, GPL16524, Agilent Technologies France, Massy, France) containing 60,306 porcine probes and derived at 71% from the porcine commercial Agilent-026440 microarray (V2, 44K, GPL15007), the remaining 29% correspond to a set of probes enriched with immune system, muscle and adipose tissue genes. Total RNA (100 ng) extracted from each whole blood samples taken at the end of the feeding trial was labelled individually with Cy3, using the One-Color Microarray-Based Gene Expression Analysis kit (Agilent Technologies) and following the manufacturer's instructions. Briefly, fluorescent complementary RNA (cRNA) was generated by a two-step procedure using T7 RNA polymerase, which simultaneously amplified target and incorporated cyanine-labeled CTP. Samples were then purified with an RNeasy mini elute kit (Qiagen, Hilden, Germany). The hybridization reactions were performed for 17 h in Agilent's SureHyb hybridization chambers containing 600 ng of Cy3-labeled cRNA per hybridization using Agilent's Gene Expression Hybridization kit. Slides were disassembled, washed according to manufacturer's instructions, and scanned at 3 µm/pixel resolution using the Agilent DNA Microarray Scanner G2505C, and images were analyzed with Agilent Feature Extraction Software (version 10.7.3.1) using the GE1_107_Sep09 extraction protocol. All microarray data have been deposited in NCBI's Gene Expression Omnibus [15] and are accessible through GEO Subserie accession number GSE70838 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE70838).

Microarray data analysis

All analyses were performed using the R software version 3.0.2 [16]. Raw spot intensities were first submitted to quality filtration based on four criteria: background intensity value, diameter, saturation and uniformity of the spot. Positives and negative controls probes were those available in the Agilent-026440 microarray. Intensities of filtered spots were log₂ transformed and median-centered. Altogether, 37,113 spots were finally retained for statistical analyses. Data were submitted to an analysis of variance considering the fixed effects of line, diet and their interaction. Data were then submitted to Benjamini and Hochberg (BH) multiple testing correction procedure. For the diet effect, no genes were found to be differentially expressed (DE) after application of the BH correction. Therefore, DE probes with a fold change (FC) cutoff higher than |1.1| between compared groups, and an uncorrected P-value below 0.01 after statistical analysis [17] were then selected for further functional analysis.

The lists of DE probes were investigated by an enrichment analysis of specific Gene Ontology (GO) terms for Biological Processes (BP), using the functional annotation clustering of Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics resources (http://david.abcc.ncifcrf.gov) [18, 19]. Due to the low number of genes found to be DE between the two dietary treatments, the functional annotation clustering was performed only for genes DE between RFI lines. The list of DE genes in low RFI *versus* high RFI pigs was divided into two lists according to their over- or under-expression. The two lists were then uploaded using the corresponding official gene symbol, when applicable, and obtained with DAVID Gene Accession Conversion Tool. The GO terms_FAT were selected, to filter the broadest terms without overshadowing the more specific ones. The P-values for enrichment were computed by a modified Fisher's exact test, using the *Homo sapiens* repository as
background. The main GOBP term was examined with an enrichment score >1.3 and P<0.1 after Benjamini Hochberg correction. For DE genes exhibiting the greatest overexpression or underexpression between RFI lines (FC > |2|) and for DE genes involved in the response to diets, a manual editing of GOBP terms was also performed using the QuickGo web-based tool [20].

## Quantitative Real-Time PCR (qPCR)

Expression of genes was further evaluated by qPCR. Final blood samples were used to validate transcriptomic results. Initial blood was used to investigate early expression of these genes. Complementary DNA was synthesized from 1 µg of total RNA, using a High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, CA, USA). Primers (S1 Table) were designed from porcine sequences using Primer Express software 3.0 (Applied Biosystems). For each primer pair, the amplification efficiency of qPCR reaction was identified using calibration curves generated with seven decreasing concentrations of cDNA from pooled RNA blood samples (obtained from 5.18 to 1E-3 ng RNA). Amplification reactions were performed in duplicate in 12.5 µL with 1 ng of reverse-transcribed RNA and both forward and reverse primers (5 µM each) in 1X PCR buffer (Fast SYBR® Green Master Mix, Applied Biosystems). A StepOnePlus Real Time PCR system (Applied Biosystems) was used. Thermal cycling conditions were as follows: 50°C for 2 min, 95°C for 20 s, followed by 40 cycles of denaturation at 95°C for 3 s, and annealing at 60°C for 30 s. Specificity of the amplification products was checked by dissociation curves analysis. This allows verifying that a single PCR product was produced. DNA topoisomerase 2-beta (TOP2B) and TATA-box-binding protein 1 (TBP1) genes were used as reference genes to calculate a normalization factor (NF) using geNorm algorithm [21]. For each gene, the normalized expression level N was calculated according to the formula developed by Pfaffl [22]:

$$N = E^{-\Delta Cq(\text{sample}-\text{calibrator})/NF}$$

where E is calculated from the slope of calibration curve and Cq is the quantification cycle and calibrator is a pool of all blood samples. E was between 1.82 and 2.10 for all studied genes.

## Statistical analyses for blood cell counts

Analysis of variance was used to determine the effects of line, diet and their interaction with the R software (version 3.0.2). A P<0.05 was retained for statistical significance, and a P<0.1 were considered as a tendency.

## Results

### Performance and blood cell count traits

Phenotypic performances have been described in detail in our previous study [14] and are briefly summarized here. Irrespective of diet, the gain-to-feed ratio was higher in low RFI pigs than in high RFI pigs; the sum of the main fatty pieces in the carcass (backfat and belly) was lower while the sum of lean pieces was higher in low RFI pigs compared with high RFI pigs. Irrespective of line, no significant difference was elicited in gain-to-feed ratio between pigs fed the HF or LF diet. Pigs fed the HF diet had a reduced BW at slaughter and their carcass fat content was lower (around 26%) than in pigs fed the LF diet. At slaughter, the total number of white blood cells was similar (P>0.1) in the different groups of pigs (Table 2). The numbers of monocytes and neutrophil granulocytes (neutrophils) did not differ between groups. The number of lymphocytes was similar in low and high RFI pigs but tended (P<0.1) to be higher in
<table><caption>Table 2. Blood cell count in pigs with low or high residual feed intake (RFI) fed a low-fat, low-fiber (LF) or high-fat, high-fiber (HF) diet.</caption><thead><tr><th rowspan="2">Traits</th><th colspan="2">Low RFI</th><th colspan="2">High RFI</th><th colspan="2">Effects</th></tr><tr><th>LF</th><th>HF</th><th>LF</th><th>HF</th><th>Line</th><th>Diet</th></tr></thead><tbody><tr><td colspan="7"><strong>WBC<sup>a</sup> (x1000.mm<sup>-3</sup>)</strong></td></tr><tr><td>Total WBC<sup>a</sup></td><td>23.8 ± 1.2</td><td>24.9 ± 1.2</td><td>23.6 ± 0.7</td><td>25.7±1.0</td><td>NS</td><td>NS</td></tr><tr><td>Neutrophils</td><td>11.1 ± 0.8</td><td>10.6 ± 0.8</td><td>10.6 ± 0.4</td><td>12.0±0.8</td><td>NS</td><td>NS</td></tr><tr><td>Lymphocytes</td><td>12.1±0.6</td><td>13.7 ± 0.7</td><td>12.3 ± 0.5</td><td>13.1±0.6</td><td>NS</td><td>&lt;0.1</td></tr><tr><td>Monocytes</td><td>0.56 ± 0.03</td><td>0.60 ± 0.04</td><td>0.63 ± 0.05</td><td>0.61±0.03</td><td>NS</td><td>NS</td></tr><tr><td colspan="7"><strong>RBC<sup>b</sup> (x10<sup>6</sup>.mm<sup>-3</sup>) and associated parameters</strong></td></tr><tr><td>Total RBC<sup>b</sup></td><td>8.14 ± 0.16</td><td>7.96 ± 0.19</td><td>7.24 ± 0.14</td><td>7.38±0.16</td><td>&lt;0.001</td><td>NS</td></tr><tr><td>Hemoglobin, g/dL</td><td>12.8 ± 0.3</td><td>12.0 ± 0.3</td><td>11.4 ± 0.3</td><td>11.4±0.2</td><td>&lt;0.001</td><td>NS</td></tr><tr><td>Hematocrit, %</td><td>47.1 ± 1.0</td><td>44.5 ± 1.1</td><td>42.8 ± 0.7</td><td>42.4±0.5</td><td>&lt;0.001</td><td>&lt;0.1</td></tr></tbody></table>

doi:10.1371/journal.pone.0146550.t002

pigs fed the HF diet than in pigs fed the LF diet. The number of red blood cells and associated parameters (hemoglobin and hematocrit) were markedly higher (P<0.001) in low RFI pigs than in high RFI pigs. These traits did not differ between diets, except a tendency for low RFI pigs fed the HF diet having a lower hematocrit than low RFI pigs fed the LF diet.

## RFI and diet effects on whole blood transcriptome

In whole blood taken at the end of the feeding trial, 2,075 probes corresponding to 1,005 genes were declared as DE (P < 0.01) between RFI lines and 82 probes corresponding to 45 genes were DE between pigs fed the HF and LF diets (Table 3). An interaction between line and diet (P < 0.01) was observed for 106 probes corresponding to 74 annotated genes. The detailed list of these corresponding genes can be found in supplementary file S2 Table.

**RFI effects on whole blood transcriptome.** Among the 2,075 DE probes (P < 0.01; adjusted P value ≤ 0.1 for 88% of the probes and ≤ 0.17 for the other probes) in response to RFI selection, 982 probes corresponding to 528 annotated unique genes were overexpressed whereas 1093 probes corresponding to 477 annotated genes were underexpressed in low RFI pigs compared with high RFI pigs. A large proportion of genes (83 to 90%) exhibited fold

<table><caption>Table 3. Number of genes differentially expressed in blood of pigs divergently selected for low or high residual feed intake (RFI) fed a low-fat, low-fiber (LF) or high-fat, high-fiber (HF) diet.</caption><thead><tr><th rowspan="2">DE genes<sup>c</sup></th><th colspan="2">Low RFI / High RFI<sup>a</sup></th><th colspan="2">HF / LF<sup>b</sup></th></tr><tr><th>Overexpressed</th><th>Underexpressed</th><th>Overexpressed</th><th>Underexpressed</th></tr></thead><tbody><tr><td>Total number</td><td>528</td><td>477</td><td>39</td><td>6</td></tr><tr><td>FC<sup>d</sup> ≥ 1.5</td><td>44</td><td>65</td><td>8</td><td>-</td></tr><tr><td>FC<sup>d</sup> ≥ 2</td><td>10</td><td>17</td><td>2</td><td>-</td></tr></tbody></table>

doi:10.1371/journal.pone.0146550.t003

<footnote><sup>a</sup>WBC: white blood cell count.</footnote>

<footnote><sup>b</sup>RBC: red blood cell count.</footnote>

<footnote>Values are means ± SEM (n = 12 pigs per diet and per line). P-value obtained from analysis of variance for the effects of line and diet. There was no significant line × diet interaction (P > 0.1). NS, not significant (P > 0.1).</footnote>

<footnote><sup>a</sup>Low RFI versus high RFI ratio, overexpressed and underexpressed in low RFI compared with high RFI pigs.</footnote>

<footnote><sup>b</sup>HF versus LF ratio, overexpressed and underexpressed in pigs fed the HF diet compared with pigs fed the LF diet.</footnote>

<footnote><sup>c</sup>DE: Differentially expressed; P < 0.01.</footnote>

<footnote><sup>d</sup>FC: Fold change ratio.</footnote>
changes in expression below |1.5| between the two lines (Table 3). Functional analysis revealed an overrepresentation of genes related to translational elongation among the overexpressed genes in the low RFI pigs compared with high RFI pigs (Table 4). Conversely, genes that were underexpressed in low RFI pigs compared with high RFI pigs shared GO terms associated to 8 different pathways: defense response, leukocyte activation, regulation of inflammatory response, negative regulation of molecular function, antigen processing and presentation, anti-apoptosis, positive regulation of immune system process and regulation of cell adhesion (Table 4).

Differentially expressed genes exhibiting the greatest overexpression or underexpression (FC > |2|) in low RFI pigs compared with high RFI lines (Table 5) were also considered to extend this global functional analysis. Among these eight genes having the greatest overexpression in the low versus high RFI lines, the majority had a documented role in the immune system (IFITM1, SLPI, IL6ST, and TRAF6); the other were mainly related to DNA organization (NPM2), translation (EIF1B), cellular homeostasis (SLCO2B1) and peptidase regulation (WFDC2). With respect to genes exhibiting the greatest underexpression in low RFI pigs

<table><caption>Table 4. Relevant GO biological processes in blood as affected by selection for RFI.</caption><thead><tr><th>GO term<sup>a</sup></th><th>Genes<sup>b</sup></th></tr></thead><tbody><tr><td colspan="2"><strong>Overexpressed genes in low RFI line compared with high RFI line</strong></td></tr><tr><td>GO:0006414 Translational elongation</td><td>RPSA, EEF1A1, RPL14, RPL15, RPL35, RPL23A, SELT, RPS3, RPS25, RPS18, RPL6, RPS13, RPL10, RPL10A, RPS21, UBA52, RPS23</td></tr><tr><td colspan="2"><strong>Underexpressed genes in low RFI line compared with high RFI line</strong></td></tr><tr><td>GO:0006952 Defense response</td><td>KYNU, NMI, FGR, LY86, CCR1, CLU, TLR1, TLR2, CXCR2, PRDX1, IL10, IL10RB, BCL2, TAP1, CSF3R, NFATC4, CD24, THBS1, BLNK, F12, TLR10, PTGER3, LY96, BECN1, HCK, IL1RN, HLA-C, C4BPB, HLA-B, CD40, C4BPA, STAT3, CD84, CD83, CD55, CD19, HIST2H2BE, CLEC7A, CLEC5A, CD14</td></tr><tr><td>GO:0045321 Leukocyte activation</td><td>ZBTB32, BST2, TLR1, TLR2, CXCR2, CD40, TPD52, SKAP2, IL10, CBLB, PKNOX1, PSEN1, ULBP1, BCL2, BCL11A, MS4A1, ADAM17, BCL6, CLEC7A, CD24, SYK, BLNK, ADAM9</td></tr><tr><td>GO:0050727 Regulation of inflammatory response</td><td>F12, ADRB2, PTGS2, SERPINF1, TGM2, CMA1, BCL6, JAK2, CD24, IL10</td></tr><tr><td>GO:0044092 Negative regulation of molecular function</td><td>NF1, PKIG, NPR3, PKIA, PROX1, PDCD4, PSMB8, PSMB9, SH3BP5, SPRY2, ADRB2, PSEN1, PSME2, CDKN2D, RGS4, PSMA3, UBC, SORT1, JAK2, UBE2D1, BUB3, DHCR24</td></tr><tr><td>GO:0019882 Antigen processing and presentation</td><td>HLA-H, ULBP1, HLA-DRB3, HLA-A, ERAP1, HLA-C, CD1A, HLA-B, HLA-DOB, PSMB8, PSMB9</td></tr><tr><td>GO:0006916 Anti-apoptosis</td><td>BECN1, SPHK1, CLU, BIRC3, ANXA4, IL10, TNFSF13B, PSEN1, BCL2, CDKN2D, TGM2, UBC, ADAM17, THBS1, DHCR24</td></tr><tr><td>GO:0002684 Positive regulation of immune system process</td><td>CLU, TLR2, C4BPB, CD40, C4BPA, CD83, CBLB, CD55, CD19, TNFSF13B, PSEN1, ADAM17, BCL6, CD24, CLEC7A, THBS1, SYK</td></tr><tr><td>GO:0030155 Regulation of cell adhesion</td><td>VAV3, GSN, BCL2, NF1, CCDC80, TGM2, BCL6, JAK2, CD24, THBS1, TPM1, ADAM9</td></tr></tbody></table>

doi:10.1371/journal.pone.0146550.t004

<footnote><sup>a</sup>Gene ontology (GO) identification number and term of the biological process. Benjamini-Hochberg adjusted P-value varied from 0.1 and 8.10<sup>-05</sup>.</footnote>

<footnote><sup>b</sup>Unique genes included in each pathway.</footnote>
<table><caption>Table 5. Top-ranked genes with fold changes in expression greater than |2| identified in the whole blood of pigs divergently selected for low or high residual feed intake (RFI).</caption><thead><tr><th>Gene symbol</th><th>Main biological process<sup>a</sup></th><th>FC<sup>b</sup></th><th>P-value<sup>c</sup></th></tr></thead><tbody><tr><td colspan="4"><strong>Overexpressed genes in low RFI line compared with high RFI line</strong></td></tr><tr><td>IFITM1</td><td>GO:0009607 response to biotic stimulus</td><td>4.34</td><td>2.6.10<sup>-04</sup></td></tr><tr><td>NPM2</td><td>GO:0006338 chromatin remodeling</td><td>3.11</td><td>3.5.10<sup>-04</sup></td></tr><tr><td>EIF1B</td><td>GO:0006412 translation</td><td>2.51</td><td>4.8.10<sup>-13</sup></td></tr><tr><td>SLPI</td><td>GO:0045071 negative regulation of viral genome replication</td><td>2.33</td><td>6.1.10<sup>-05</sup></td></tr><tr><td>IL6ST</td><td>GO:0019221 cytokine-mediated signaling pathway</td><td>2.20</td><td>1.3.10<sup>-03</sup></td></tr><tr><td>WFDC2</td><td>GO:0010466 negative regulation of peptidase activity</td><td>2.16</td><td>4.5.10<sup>-08</sup></td></tr><tr><td>TRAF6</td><td>GO:0002726 positive regulation of T cell cytokine production</td><td>2.08</td><td>1.4.10<sup>-06</sup></td></tr><tr><td>SLCO2B1</td><td>GO:0006811 ion transport</td><td>2.04</td><td>6.9.10<sup>-06</sup></td></tr><tr><td colspan="4"><strong>Underexpressed genes in low RFI line compared with high RFI line</strong></td></tr><tr><td>GPX3</td><td>GO:0006979 response to oxidative stress</td><td>-6.07</td><td>1.3.10<sup>-06</sup></td></tr><tr><td>HMG20A</td><td>GO:0006338 chromatin remodeling</td><td>-4.69</td><td>1.2.10<sup>-08</sup></td></tr><tr><td>CAPNS1</td><td>GO:0006508 proteolysis</td><td>-4.09</td><td>3.1.10<sup>-07</sup></td></tr><tr><td>OAZ3</td><td>GO:0043086 negative regulation of catalytic activity</td><td>-4.00</td><td>1.1.10<sup>-06</sup></td></tr><tr><td>NMI</td><td>GO:0045824 negative regulation of innate immune response</td><td>-3.49</td><td>1.6.10<sup>-07</sup></td></tr><tr><td>SERPINF1</td><td>GO:0016525 negative regulation of angiogenesis</td><td>-3.28</td><td>3.8.10<sup>-08</sup></td></tr><tr><td>HTRA1</td><td>GO:0030512 negative regulation of transforming growth factor beta receptor signaling pathway</td><td>-2.68</td><td>1.1.10<sup>-09</sup></td></tr><tr><td>NPR3</td><td>GO:0008217 regulation of blood pressure</td><td>-2.17</td><td>1.4.10<sup>-07</sup></td></tr><tr><td>CHIT1</td><td>GO:0006030 chitin metabolic process</td><td>-2.12</td><td>4.6.10<sup>-07</sup></td></tr><tr><td>SLA-DOA</td><td>GO:0019882 antigen processing and presentation</td><td>-2.12</td><td>1.9.10<sup>-11</sup></td></tr></tbody></table>

doi:10.1371/journal.pone.0146550.t005

compared with high RFI pigs, their roles were also associated with immune system (NMI, CHIT1, SLA-DOA), oxidative process (GPX3, CAPNS1, OAZ3), DNA organization (HMG20A), angiogenesis and the control of blood volume (SERPINF1, HTRA1, NPR3).

**Diet effects on whole blood transcriptome.** The number of DE probes in response to diets was markedly lower than that observed in response to genetic selection (Table 3). Among the 82 DE probes between diets, 62 probes corresponding to 39 annotated genes were overexpressed, whereas 20 probes corresponding to 6 annotated genes were underexpressed in pigs fed the HF diet compared with pigs fed the LF diet. The majority of these genes (78%) had fold change in expression between the two diets below |1.5| (Table 3). As indicated in Table 6, among the 16 genes overexpressed in pigs fed the HF diet compared with pigs fed the LF diet, six genes were related to the immune system (PR39, PMAP-23, CD5L, TCN1, PGLYRP1 and HP), three genes were involved in cellular metabolism such as collagen fiber assembly or cytoskeletal arrangement (BGN, PLEK2, RASL11B), two genes were involved in lipid storage (CPT1A, LCN2), two genes had a role in erythropoietic system (EPOR, KLF1), two genes had a role in DNA synthesis (NOCL4; REXO2) and one gene was involved in embryogenesis (SHISA2). Conversely, three genes were underexpressed in pigs fed the HF diet and were involved in lipid metabolism (PSAP), RNA localization (RAE1), and nucleotides phosphorylation (NTPCR).

<footnote><sup>a</sup>The gene ontology (GO) term for biological process was manually obtained from the QuickGO web-based tool for each gene name.</footnote>

<footnote><sup>b</sup>Fold change value is expressed as the expression ratio of low RFI pigs versus high RFI pigs. Ratio was inversed and preceded by a minus sign for value less than 1 (i.e., a ratio of 0.5 is expressed as -2).</footnote>

<footnote><sup>c</sup>P-value obtained from analysis of variance for the effect of line. The highest P-value is reported when several probes are differentially expressed for a unique gene. Adjusted P-value varied from 0.05 and 1.10<sup>-09</sup>.</footnote>
<table><caption>Table 6. Differentially expressed genes in pigs fed diets with contrasted energy source and nutrients.</caption><thead><tr><th>Gene symbol</th><th>Main biological process<sup>a</sup></th><th>FC<sup>b</sup></th><th>P-value<sup>c</sup></th></tr></thead><tbody><tr><td colspan="4"><strong>Overexpressed genes in HF pigs compared with LF pigs</strong></td></tr><tr><td>PR39</td><td>GO:0042742 defense response to bacterium</td><td>3.06</td><td>6.0.10<sup>-03</sup></td></tr><tr><td>PMAP-23</td><td>GO:0042742 defense response to bacterium</td><td>2.84</td><td>7.3.10<sup>-03</sup></td></tr><tr><td>CD5L</td><td>GO:0006898 receptor-mediated endocytosis</td><td>1.72</td><td>3.3.10<sup>-03</sup></td></tr><tr><td>LCN2</td><td>GO:0006810 transport</td><td>1.61</td><td>8.4.10<sup>-03</sup></td></tr><tr><td>TCN1</td><td>GO:0015889 cobalamin transport</td><td>1.60</td><td>9.2.10<sup>-03</sup></td></tr><tr><td>BGN</td><td>GO:0019800 peptide cross-linking via chondroitin 4-sulfate glycosaminoglycan</td><td>1.53</td><td>9.0.10<sup>-03</sup></td></tr><tr><td>PGLYRP1</td><td>GO:0050728 negative regulation of inflammatory response</td><td>1.51</td><td>8.1.10<sup>-03</sup></td></tr><tr><td>EPOR</td><td>GO:0038162 erythropoietin-mediated signaling pathway</td><td>1.46</td><td>5.5.10<sup>-03</sup></td></tr><tr><td>HP</td><td>GO:0006508 proteolysis</td><td>1.46</td><td>9.7.10<sup>-03</sup></td></tr><tr><td>PLEK2</td><td>GO:0035556 intracellular signal transduction</td><td>1.33</td><td>6.6.10<sup>-03</sup></td></tr><tr><td>RASL11B</td><td>GO:0007264 small GTPase mediated signal transduction</td><td>1.29</td><td>1.1.10<sup>-04</sup></td></tr><tr><td>SHISA2</td><td>GO:0007275 multicellular organismal development</td><td>1.28</td><td>8.3.10<sup>-04</sup></td></tr><tr><td>NOC4L</td><td>GO:0042254 ribosome biogenesis</td><td>1.27</td><td>1.0.10<sup>-03</sup></td></tr><tr><td>KLF1</td><td>GO:0030218 erythrocyte differentiation</td><td>1.25</td><td>6.4.10<sup>-03</sup></td></tr><tr><td>CPT1A</td><td>GO:0032000 positive regulation of fatty acid β-oxidation</td><td>1.24</td><td>1.8.10<sup>-03</sup></td></tr><tr><td>REXO2</td><td>GO:0090305 nucleic acid phosphodiester bond hydrolysis</td><td>1.21</td><td>9.4.10<sup>-03</sup></td></tr><tr><td colspan="4"><strong>Underexpressed genes in HF pigs compared with LF pigs</strong></td></tr><tr><td>NTPCR</td><td>GO:0016311 dephosphorylation</td><td>-1.29</td><td>4.9.10<sup>-03</sup></td></tr><tr><td>RAE1</td><td>GO:0071407 cellular response to organic cyclic compound</td><td>-1.24</td><td>5.2.10<sup>-03</sup></td></tr><tr><td>PSAP</td><td>GO:0006629 lipid metabolic process</td><td>-1.20</td><td>7.3.10<sup>-03</sup></td></tr></tbody></table>

doi:10.1371/journal.pone.0146550.t006

## Data validation by qPCR analysis

The expression of eight genes identified as DE between RFI lines by transcriptome analysis was checked by qPCR (Fig 1). With the exception of *TRAF6* (P = 0.16), the range of difference (P < 0.05) in the expression of *CD40*, *GPX3*, *OAZ3*, *DGAT2*, *NMI*, *PSEN1* and *SLPI* genes was consistent with microarray analyses. Expressions of those eight genes were also determined in 74-day-old pigs (Fig 1). The differences in expression levels between RFI lines observed at 132 days of age were detected at the onset of the growing period for *CD40*, *GPX3*, *OAZ3*, *DGAT2* and *SLPI* genes (P < 0.05). For three other genes (*NMI*, *TRAF6*, *PSEN1*), there was no difference (P > 0.1) between low RFI and high RFI pigs at this early stage of growth.

Finally, the expression levels of three genes (*LCN2*, *CPT1A*, *PSAP*) identified as DE between pigs fed the two diets by transcriptome analysis were also determined by qPCR at 132 days of age (end of the feeding trial), and the range of difference was consistent (P < 0.001) with microarray analyses for *LCN2* and *CPT1A*, while *PSAP* expression by qPCR did not significantly differ (P > 0.10) between HF and LF pigs (Fig 2).

## Discussion

The current study provides new evidence that the investigation of the expression of genes in whole blood allowed clear identification of biological pathways and indicator traits involved in

<footnote><sup>a</sup>The gene ontology (GO) term for biological process (BP) was manually obtained from the QuickGO web-based tool for each gene name.</footnote>

<footnote><sup>b</sup>Fold change (FC) value is expressed as the expression ratio of pigs fed the HF (high fiber, high fat) diet versus pigs fed the LF (low fat low fiber) diet. FC is inversed and preceded by a minus sign for value less than 1 (i.e. a ratio of 0.5 is expressed as -2).</footnote>

<footnote><sup>c</sup>P-value obtained from analysis of variance for the effect of diet. The highest P-value is reported when several probes are differentially expressed for a unique gene.</footnote>
<figure><img src="image_4.png" /></figure>

<table><caption>**Fig 1. Comparison of microarray and qPCR data for eight genes in pigs divergently selected for low or high residual feed intake (RFI).**</caption><thead><tr><th></th><th>CD40</th><th>GPX3</th><th>OAZ3</th><th>DGAT2</th><th>NMI</th><th>TRAF6</th><th>SLPI</th><th>PSEN1</th></tr></thead><tbody><tr><td>P-value</td><td>*</td><td>**</td><td>**</td><td>*</td><td>&gt;0.1</td><td>&gt;0.1</td><td>*</td><td>&gt;0.1</td></tr><tr><td>qPCR 74 days of age</td><td>***</td><td>***</td><td>***</td><td>**</td><td>*</td><td>&gt;0.1</td><td>**</td><td>*</td></tr><tr><td>qPCR 132 days of age</td><td>**</td><td>***</td><td>***</td><td>**</td><td>***</td><td>***</td><td>***</td><td>**</td></tr><tr><td>Microarray</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Transcriptomic differences between the two lines were validated by qPCR at 132 and 74 days of age (onset of the growing period). For values related to microarray, the highest P-value is reported when several probes are differentially expressed for a unique gene (*P < 0.05; **P < 0.01; ***P<0.001). Fold change value is expressed as the expression ratio of low RFI to high RFI samples; ratio was inversed and preceded by a minus sign for value less than 1 (i.e., a ratio of 0.5 is expressed as -2). CD40, tumor necrosis factor receptor superfamily member 5; GPX3, glutathione peroxidase 3; OAZ3, ornithine decarboxylase antizyme 3; DGAT2, diacylglycerol O-acyltransferase 2; NMI, N-myc interactor; TRAF6, TNF receptor-associated factor 6; SLPI, secretory leukocyte peptidase inhibitor; PSEN1, presenilin 1.

doi:10.1371/journal.pone.0146550.g001

the response to genetics and nutrition strategies. To date, the majority of studies including a microarray analysis from whole blood or PBMC have been performed in adults exhibiting several disorders [23]. Therefore, the current study provides novel information on the response of blood transcriptome during growth. This study further supports the interest to use RNA extracted from whole blood which is easier to obtain compared with RNA extracted from PBMC which requires supplementary steps for cell isolation. In the current study, there was no significant difference in the numbers of lymphocytes, monocytes and neutrophils between the
<figure><img src="image_5.png" /><figcaption>**Fig 2. Comparison of microarray and qPCR data for three genes in pigs fed the HF or LF diet.** For values related to microarray, the highest P-value is reported when several probes are differentially expressed for a unique gene (*P* < 0.05; **P* < 0.01; ***P* < 0.001). Fold change value is expressed as the expression ratio of HF (high fiber high fat, n = 24) to LF (low fiber low fat, n = 24) diets; ratio was inversed and preceded by a minus sign for value less than 1 (i.e., a ratio of 0.5 is expressed as -2). *LCN2*, lipocalin 2; *CPT1A*, Carnitine palmitoyltransferase 1A; *PSAP*, prosaposin.</figcaption></figure>

doi:10.1371/journal.pone.0146550.g002

experimental groups, so that identified DE genes after microarray analysis were not related to differences in the number of WBC. These observations differ from a recent study reporting a lower number of WBC, especially lymphocytes and monocytes, in Yorkshire pigs selected for low RFI compared with pigs selected for high RFI and considered at younger age (35-42 days) [24]. With respect to RBC, their total number was higher in low RFI than in high RFI pigs in the current study, while this number did not differ between the divergent RFI lines in the study of Mpetile et al. [24]. It is important to note that other parameters related to RBC were similarly affected by RFI selection in both studies, with higher hemoglobin plasma concentration and higher hematocrit percentage in low RFI pigs than in high RFI pigs. The current findings may reflect differences in intensity of selection, age of pigs considered for the measurement of blood cell traits or environmental stimuli during pig housing.

The current study clearly indicates that genetics has a greater impact on blood gene expression profiles than dietary treatment. In a rat model, genetic background has been shown to have a much larger impact on PBMC transcriptome than on heart transcriptome [25]. With
the current findings, it can be hypothesized that the identified DE genes in whole blood may reflect changes that occur in a similar manner in several tissues of the organism in response to genetic selection for RFI. Indeed, it has been shown that the expression of a large number of genes was shared among all tissues including blood cells [10].

The most important pathways shared by genes DE between RFI lines in the blood concern immunity and defense mechanisms. Most of these genes were underexpressed in low RFI pigs compared with high RFI pigs, although some of them (IFITM1, SLPI, and IL6ST) were rather overexpressed in low RFI pigs. For instance, the CD40 gene coding for a TNF receptor super-family member required for the B-cell function was found to exhibit a lower expression in the low RFI pigs compared with the high RFI pigs at both examined ages. The expression level of PSEN1, a gene participating to T cell activation, was decreased only in 132-day-old pigs. Conversely, the SLPI gene, coding for an antimicrobial protein having inhibitory effects that contribute to the immune response by protecting epithelial surfaces, was overexpressed in low RFI pigs at both ages. An overexpression of the IL6ST gene was also observed in low RFI pigs, a gene encoding a signal transducer shared by many cytokines including interleukin 6. Finding many genes related to immunity in the whole blood transcriptome is expected considering the fact that blood cells constitute one of the first lines of immune defense [10]. A recent study in pigs showing that peripheral blood transcriptome is a relevant source to identify genes related to the immune function responses and to predict the efficiency of individual's immune response further supports our data [26]. From the current data, it is however not possible to predict this efficiency of immune response in RFI pigs. It has been suggested that livestock animals selected for high production efficiency traits may be more vulnerable to diseases and stressors [27, 28]. Nevertheless, the few available data dealing with the inflammatory response of the two RFI lines do not allow a clear conclusion [29-30]. A lower level of basal inflammation, as indicated by reduced concentration of the acute phase protein haptologin, has been recently reported in low RFI pigs compared with high RFI pigs [29]. In response to an immune challenge induced by infection with the porcine reproductive and respiratory syndrome virus [31], pigs selected for a low RFI compared with high RFI pigs had a greater immune response and their growth rate was less affected. Further studies are thus needed to confirm the higher capacity of low RFI pigs and to get a better understanding of the immune capacity and inflammatory response of the two RFI lines.

Genes associated with translational elongation were overexpressed in low RFI pigs compared with high RFI pigs, including *eiF1B* gene involved in translational initiation as one of the top overexpressed genes in the former pigs. Detectable expression of *eiF1B* gene has been previously demonstrated after microarray analyses in normal whole blood from Human patients [32]. An increase in the expression of genes encoding initiation and elongation factor subunits has been also recently reported in skeletal muscle of low RFI pigs of the sixth and seventh generations of the same selection experiment [6]. These findings may reflect a greater protein synthesis in low RFI pigs than in high RFI pigs, to support increased muscle gain in the former pigs during the growing period [5]. This remains to be further investigated. Indeed, translation initiation signaling proteins did not differ in skeletal muscle between low and high RFI lines [4, 33].

The expression of genes encoding glutathione peroxidase 3 (GPX3) and ornithine decarboxylase antizyme 3 (OAZ3) was also found as DE in low *versus* high RFI pigs at both ages. The same genes have been previously reported to be DE in skeletal muscle between pigs of the former generations of selection [6]. Those two genes have a role in oxidative metabolism. GPX3, encoding an enzyme involved in the detoxification of hydrogen peroxide when reactive oxygen species were produced, accounts for the major antioxidant activity in the plasma. OAZ3, encoding an enzyme which is an inhibitor of ornithine decarboxylase, itself involved in the
anti-oxidant defense, converts ornithine in putrescine, the first step in the synthesis of polyamines that have antioxidant properties. These findings argue for a lower oxidative stress in low RFI pigs [34] and the expression of these genes may be considered as relevant indicators of RFI status in the pig.

Only few genes were found DE in whole blood samples collected from pigs fed the HF or LF diets, despite a large impact of the dietary treatment on growth and body composition of pigs from the two RFI lines [14, 35]. Between normal weight and diet-induced obese (cafeteria diet) rats, exhibiting large differences in body weight and adiposity index, the number of genes affected by the diet in PBMC (566 genes; [36]) was much higher than in the present study. The difference in the extent of the gene response between our study and the above-mentioned ones are unknown. This may be related to differences in nutrient composition or energy level of the diet. Alternatively, the difference may be related to the fact that an increase in adipose tissue mass occurred in those studies, whereas there was a decrease in the fat mass of pigs fed the HF diet compared with pigs fed the LF diet.

Despite the low number of genes DE in response to the HF compared with the LF diet, the current study highlighted differences in the expression of genes related to the immune system and lipid metabolism as shown reported in PBMC of miniature pigs fed a hyperlipidic diet *versus* hyperlipidic and hyperglucidic diet for 27 weeks [37]. Overexpressed genes in pigs fed the HF diet were mainly involved in host defense mechanism (*PR39*, *PMAP-23*, *CD5L*, *PGLYRP1*). With respect to lipid metabolism, two genes were also identified as DE in whole blood of pigs fed the two diets. The *CPT1A* (carnitine palmitoyl transferase 1A) gene, encoding an enzyme which is localized in the inner membrane of the mitochondria and participates to the entrance of fatty acids in the mitochondria for β-oxidation, was overexpressed in HF pigs compared with LF pigs. This may reflect a higher fatty acid oxidative capacity in blood cells of HF pigs, as observed in rats [38]. Other studies have shown that *CPT1A* gene expression in PBMC did not vary between overweight and control rats [39], or between obese and normal weight subjects in human [40]. Conversely, the expression of *CPT1A* was higher in overweight children than in normal weight children [41]. In our study, the fact that HF pigs were leaner and lighter than LF pigs [14] suggests no direct relationships between *CPT1A* expression in peripheral blood and changes in adipose tissue mass.

The gene encoding lipocalin 2 (LCN2) was also overexpressed in pigs fed the HF diet compared with pigs fed the LF diet. Lipocalin 2 also named neutrophil gelatinase-associated lipocalin (NGAL) has been recently characterized as an adipokine playing putative roles in glucose and lipid metabolism [42, 43]. This adipokine shares structural similarities with fatty acid binding proteins and retinol-binding proteins. Recent studies have shown a higher LCN2 expression level in PBMC of adults exhibiting higher adiposity [44] and (or) higher body weight and body mass index [45]. With this contrasted findings, further studies are needed to determine whether there is a possible relationship between peripheral blood LCN2 expression and adipose tissue mass variation.

In summary, the current study based on the investigation of pigs selected for feed efficiency and fed two diets contrasted in energy sources indicate that transcriptomics analyses is an additional tool to study variations of phenotypes in a dynamic way throughout the life of the animal. It also supports the potential use of blood transcriptome to highlight biomarkers for future selection process in pig or in other species.

**Supporting Information**

S1 Table. Primer sequences used for analysis of gene expression by qPCR.
(DOCX)
S2 Table. Probes expressed in whole blood of pigs selected for low or high residual feed intake (RFI) fed a low-fat, low-fiber (LF) or high-fat, high-fiber (HF) diet, with significant interaction between line and diet.

(XLSX)

## Acknowledgments

The authors thank Y. Billon and A. Priet (INRA, UE1372 GenESI, Le Magneraud, F-17700 Sur-gères) for line selection, B. Carrisant, H. Demay and J. Liger (INRA, UMR1348 Pegase) for blood sampling, and G. Guillemois (INRA, UMR1348 Pegase) for diet preparation. They are also grateful to E. Baéza (INRA, UR83 Recherches Avicoles, F-37380 Nouzilly) for her critical reading of the manuscript.

## Author Contributions

Conceived and designed the experiments: FG IL. Performed the experiments: MJ FG AV CT HG IL. Analyzed the data: MJ AV. Contributed reagents/materials/analysis tools: MJ AV CT. Wrote the paper: MJ FG HG IL.

## References

1. Niemann H, Kuhla B, Flachowsky G. Perspectives for feed-efficient animal production. J Anim Sci. 2011; 89: 4344-4363. doi: 10.2527/jas.2011-4235 PMID: 21705634

2. Animal Task Force. Research & Innovation for a Sustainable Livestock Sector in Europe. An Animal Task Force White Paper. 2013 [27th of July 2015]; Available from: http://www.animaltaskforce.eu/Portals/0/ATF/documents%20for%20scare/ATF%20white%20paper%20Research%20priorities%20for%20a%20sustainable%20livestock%20sector%20in%20Europe.pdf. Assessed the 15th of September 2015.

3. Gilbert H, Bidanel JP, Gruand J, Caritez JC, Billon Y, Guillouet P, et al. Genetic parameters for residual feed intake in growing pigs, with emphasis on genetic relationships with carcass and meat quality traits. J Anim Sci. 2007; 85: 3182-3188. doi: 10.2527/jas.2006-590 PMID: 17785600

4. Le Naou T, Le Floc'h N, Louveau I, Gilbert H, Gondret F. Metabolic changes and tissue responses to selection on residual feed intake in growing pigs. J Anim Sci. 2012; 90: 4771-4780. doi: 10.2527/jas.2012-5226 PMID: 22871936

5. Faure J, Lefaucheur L, Bonhomme N, Ecolan P, Meteau K, Coustard SM, et al. Consequences of diver-gent selection for residual feed intake in pigs on muscle energy metabolism and meat quality. Meat Sci. 2013; 93:37-45. doi: 10.1016/j.meatsci.2012.07.006 PMID: 22910803

6. Vincent A, Louveau I, Gondret F, Tréfeu C, Gilbert H, Lefaucheur L. Divergent selection for residual feed intake affects the transcriptomic and proteomic profiles of pig skeletal muscle. J Anim Sci. 2015; 93: 2745-2758. doi: 10.2527/jas2015-8928 PMID: 26115262

7. Montagne L, Loisel F, Le Naou T, Gondret F, Gilbert H, Le Gall M. Difference in short-term responses to a high-fiber diet in pigs divergently selected for residual feed intake. J Anim Sci. 2014; 92: 1512-1523. doi: 10.2527/jas.2013-6623 PMID: 24496835

8. Arkfeld EK, Young JM, Johnson RC, Fedler CA, Prusa K, Patience JF, et al. Composition and quality characteristics of carcasses from pigs divergently selected for residual feed intake on high- or low-energy diets. J Anim Sci. 2015; 93: 2530-2545. doi: 10.2527/jas2014-8546 PMID: 26020348

9. Herd RM, Arthur PF. Physiological basis for residual feed intake. J Anim Sci. 2009; 87 (14 Suppl): E64-E71. doi: 10.2527/jas.2008-1345 PMID: 19028857

10. Liew CC, Ma J, Tang HC, Zheng R, Dempsey AA. The peripheral blood transcriptome dynamically reflects system wide biology: a potential diagnostic tool. J Lab Clin Med. 2006; 147: 126-132. doi: 10.1016/j.lab.2005.10.005 PMID: 16503242

11. Mohr S, Liew CC. The peripheral-blood transcriptome: new insights into disease and risk assessment. Trends Mol Med. 2007; 13: 422-432. doi: 10.1016/j.molmed.2007.08.003 PMID: 17919976

12. Ghosh S, Dent R, Harper ME, Gorman SA, Stuart JS, McPherson R. Gene expression profiling in whole blood identifies distinct biological pathways associated with obesity. BMC Med Genomics. 2010; 3: 56. doi: 10.1186/1755-8794-3-56 PMID: 21122113
13. Diaz-Rua R, Keijer J, Caimari A, van Schothorst EM, Palou A, Oliver P. Peripheral blood mononuclear cells as a source to detect markers of homeostatic alterations caused by the intake of diets with an unbalanced macronutrient composition. J Nutr Biochem. 2015; 26: 398-407. doi: 10.1016/j.jnutbio. 2014.11.013 PMID: 25660595

14. Gondret F, Louveau I, Mourot J, Duclos MJ, Lagarrigue S, Gilbert H, et al. Dietary energy sources affect the partition of body lipids and the hierarchy of energy metabolic pathways in growing pigs differing in feed efficiency. J Anim Sci. 2014; 92: 4865-4877. doi: 10.2527/jas.2014-7995 PMID: 25253805

15. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30: 207-210. PMID: 11752295

16. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria 2008; Available from: http://www.R-project.org.

17. Fu WJ, Stromberg AJ, Viele K, Carroll RJ, Wu G. Statistics and bioinformatics in nutritional sciences: analysis of complex data in the era of systems biology. J Nutr Biochem. 2010; 21: 561-572. doi: 10.1016/j.jnutbio.2009.11.007 PMID: 20233650

18. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Protocols. 2009; 4: 44-57. doi: 10.1038/nprot.2008.211 PMID: 19131956

19. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009; 37: 1-13. doi: 10.1093/nar/gkn923 PMID: 19033363

20. Binns D, Dimmer E, Huntley R, Barrell D, O'Donovan C, Apweiler R. QuickGO: a web-based tool for Gene Ontology searching. Bioinformatics. 2009; 25: 3045-3046. doi: 10.1093/bioinformatics/btp536 PMID: 19744993

21. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002; 3: 18.

22. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001; 29: e45. PMID: 11328886

23. Pascual V, Chaussabel D, Banchereau J. A genomic approach to human autoimmune diseases. Annu Rev Immunol. 2010; 28: 535-571. doi: 10.1146/annurev-immunol-030409-101221 PMID: 20192809

24. Mpetile Z, Young JM, Gabler NK, Dekkers JC, Tuggle CK. Assessing peripheral blood cell profile of Yorkshire pigs divergently selected for residual feed intake. J Anim Sci. 2015; 93: 892-899. doi: 10.2527/jas2014-8132 PMID: 26020867

25. Deans KJ, Minneci PC, Chen H, Kern SJ, Logun C, Alsaaty S, et al. Impact of animal strain on gene expression in a rat model of acute cardiac rejection. BMC Genomics. 2009; 10: 1471-2164. doi: 10.1186/1471-2164-10-280

26. Mach N, Gao Y, Lemonnier G, Lecardonnel J, Oswald IP, Estelle J, et al. The peripheral blood transcriptome reflects variations in immunity traits in swine: towards the identification of biomarkers. BMC Genomics. 2013; 14:894. doi: 10.1186/1471-2164-14-894 PMID: 24341289

27. Rauw WM, Kanis E, Noordhuizen-Stassen EN, Grommers FJ. Undesirable side effects of selection for high production efficiency in farm animals: a review. Livest Prod Sci. 1998; 56: 15-33.

28. Prunier A, Heinonen M, Quesnel H. High physiological demands in intensively raised pigs: impact on health and welfare. Animal. 2010; 4: 886-898. doi: 10.1017/S175173111000008X PMID: 22444261

29. Mani V, Harris AJ, Keating AF, Weber TE, Dekkers JC, Gabler NK. Intestinal integrity, endotoxin transport and detoxification in pigs divergently selected for residual feed intake. J Anim Sci. 2013; 91: 2141-2150. doi: 10.2527/jas2012-6053 PMID: 23463550

30. Labussière E, Dubois S, Gilbert H, Thibault JN, Le Floc'h N, Noblet J, et al. Effect of inflammation stimulation on energy and nutrient utilization in piglets selected for low and high residual feed intake. Animal. 2015. doi: 10.1017/S1751731115000932

31. Dunkelberger JR, Boddicker NJ, Serão NVL, Young JM, Rowland RRR, Dekkers JCM. Response of pigs divergently selected for residual feed intake to experimental infection with the PRRS virus. Livest Sci. 2015; 177: 132-41.

32. Demirci H, Reed D, Elner VM. Tissue-based microarray expression of genes predictive of metastasis in uveal melanoma and differentially expressed in metastatic uveal melanoma. J Ophthalmic Vis Res. 2013; 8: 303-307. PMID: 24653816

33. Cruzen SM, Harris AJ, Hollinger K, Punt RM, Grubbs JK, Selsby JT, et al. Evidence of decreased muscle protein turnover in gilts selected for low residual feed intake. J Anim Sci. 2013; 91: 4007-4016. doi: 10.2527/jas2013-6413 PMID: 23739790
34. Grubbs JK, Huff-Lonergan E, Gabler NK, JC MD, Lonergan SM. Liver and skeletal muscle mitochondria proteomes are altered in pigs divergently selected for residual feed intake. J Anim Sci. 2014; 92: 1995–2007. doi: 10.2527/jas2013-7391 PMID: 24671593

35. Jégou M, Gondret F, Lalande-Martin J, Tea I, Baeza E, Louveau I. NMR-based metabolomics highlights differences in plasma metabolites in pigs exhibiting diet-induced differences in adiposity. Eur J Nutr. 2015 May 22. doi: 10.1007/s00394-015-0932-z

36. Caimari A, Oliver P, Rodenburg W, Keijer J, Palou A. Slc27a2 expression in peripheral blood mononuclear cells as a molecular marker for overweight development. Int J Obes (Lond). 2010; 34: 831–839. doi: 10.1038/ijo.2010.17

37. Takahashi J, Waki S, Matsumoto R, Odake J, Miyaji T, Tottori J, et al. Oligonucleotide microarray analysis of dietary-induced hyperlipidemia gene expression profiles in miniature pigs. PLoS One. 2012; 7: e37581. doi: 10.1371/journal.pone.0037581 PMID: 22662175

38. Konieczna J, Sanchez J, van Schothorst EM, Torrens JM, Bunschoten A, Palou M, et al. Identification of early transcriptome-based biomarkers related to lipid metabolism in peripheral blood mononuclear cells of rats nutritionally programmed for improved metabolic health. Genes Nutr. 2014; 9: 366. doi: 10.1007/s12263-013-0366-2 PMID: 24343050

39. Oliver P, Reynes B, Caimari A, Palou A. Peripheral blood mononuclear cells: a potential source of homeostatic imbalance markers associated with obesity development. Pflugers Archiv: Eur J Physiol. 2013; 465: 459–468. doi: 10.1007/s00424-013-1246-8

40. Telle-Hansen VH, Halvorsen B, Dalen KT, Narverud I, Wesseltoft-Rao N, Granlund L, et al. Altered expression of genes involved in lipid metabolism in obese subjects with unfavourable phenotype. Genes Nutr. 2013; 8: 425–434. doi: 10.1007/s12263-012-0329-z PMID: 23296345

41. Sanchez J, Priego T, Pico C, Ahrens W, De Henauw S, Fraterman A, et al. Blood cells as a source of transcriptional biomarkers of childhood obesity and its related metabolic alterations: results of the IDE-FICS study. J Clin Endocrinol Metab. 2012; 4: E648–E652.

42. Guo H, Bazuine M, Jin D, Huang MM, Cushman SW, Chen X. Evidence for the regulatory role of lipocalin 2 in high-fat diet-induced adipose tissue remodeling in male mice. Endocrinology. 2013; 154: 3525–3538. doi: 10.1210/en.2013-1289 PMID: 23885013

43. Wu G, Li H, Zhou M, Fang Q, Bao Y, Xu A, et al. Mechanism and clinical evidence of lipocalin-2 and adipocyte fatty acid-binding protein linking obesity and atherosclerosis. Diabetes Metab Res Rev. 2014; 30: 447–456. doi: 10.1002/dmrr.2493 PMID: 24214285

44. Catalan V, Gomez-Ambrosi J, Rodriguez A, Ramirez B, Valenti V, Moncada R, et al. Peripheral mononuclear blood cells contribute to the obesity-associated inflammatory state independently of glycemic status: involvement of the novel proinflammatory adipokines chemerin, chitinase-3-like protein 1, lipocalin-2 and osteopontin. Genes Nutr. 2015; 10: 460. doi: 10.1007/s12263-015-0460-8 PMID: 25869413

45. Berisha SZ, Serre D, Schauer P, Kashyap SR, Smith JD. Changes in whole blood gene expression in obese subjects with type 2 diabetes following bariatric surgery: a pilot study. PLoS One. 2011; 6: e16729. doi: 10.1371/journal.pone.0016729 PMID: 21423737