# Spontaneous intake of essential oils after a negative postnatal experience has long-term effects on blood transcriptome in chickens

Aline Foury, Anne Collin, Jean-Christophe Helbling, Christine Leterrier, Marie-Pierre Moisan, Laurence L.A. Guilloteau

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

## ► To cite this version:

Aline Foury, Anne Collin, Jean-Christophe Helbling, Christine Leterrier, Marie-Pierre Moisan, et al. Spontaneous intake of essential oils after a negative postnatal experience has long-term effects on blood transcriptome in chickens. Scientific Reports, 2020, 10 (1), pp.1-14. 10.1038/s41598-020-77732-5. hal-03179256

HAL Id: hal-03179256

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

Submitted on 24 Mar 2021

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>

Check for updates

OPEN

# Spontaneous intake of essential oils after a negative postnatal experience has long-term effects on blood transcriptome in chickens

Aline Foury<fnref n="1" />, Anne Collin², Jean-Christophe Helbling<fnref n="1" />, Christine Leterrier³, Marie-Pierre Moisan<fnref n="1" />,⁴ & Laurence A. Guilloteau²,⁴✉

Chicks subjected to early stressful factors could develop long-lasting effects on their performances, welfare and health. Free access to essential oils (EO) in poultry farming could mitigate these effects and potentially reduce use of antimicrobial drugs. This study on chicken analyzed long-lasting effects of post-hatch adverse conditions (Delayed group), and the impact of EO intake on blood physiological parameters and transcriptome. Half of the Control and Delayed groups had free access to EO, while the other half had only water for the first 13 days post-hatching. Blood analyses of metabolites, inflammation and oxidative stress biomarkers, and mRNA expression showed sex differences. Long-lasting effects of postnatal experience and EO intake persisted in blood transcriptome at D34. The early adverse conditions modified 68 genes in males and 83 genes in females. In Delayed males six transcription factors were over-represented (NFE2L2, MEF2A, FOXI1, Foxd3, Sox2 and TEAD1). In females only one factor was over-represented (PLAG1) and four under-represented (NFIL3, Foxd3, ESR2 and TAL1::TCF3). The genes showing modified expression are involved in oxidative stress, growth, bone metabolism and reproduction. Remarkably, spontaneous EO intake restored the expression levels of some genes affected by the postnatal adverse conditions suggesting a mitigating effect of EO intake.

In poultry production systems, broiler chicks can be exposed to various stressful factors in the hatchery, during transportation to rearing houses and in their first days of life. These early-life stresses have long-lasting effects on their performances, but also on welfare and health<fnref n="1" />,² often requiring the use of antibiotics. Alternatives to antibiotics treatment used in animal production is a field of intense research, driven by the emergence of resistance to bacteria and the societal demand for more sustainable rearing practices and veterinary medicines. An appealing alternative to antibiotics/chemical antimicrobial drugs is the use of specific plants or substrates with medicinal properties to control or to prevent diseases. In particular, essential oils (EO) extracted from aromatic plants are known to offer multi-functional medicinal properties including antimicrobial, antioxidant, anti-inflammatory, immune and nervous regulatory properties³⁻⁶. These regulatory properties are linked to the terpenoids (monoterpenes, sesquiterpenes) and aromatic compounds present in EO. Phenols, alcohols and ketones have also been reported to have an antibacterial action⁷. In chickens, EO have been employed as feed additives to promote growth and health⁸,⁹. However, in most studies, the EO were included in the feed of the animal and their ability to mitigate long-term effects of early-life stresses was never analyzed.

In a previous study, we developed an experimental model reproducing adverse perinatal conditions in chicks including the deprivation of feed and water, unpredictable shaking in transportation boxes and temperature changes 24 h after hatching, summarized as “negative postnatal experience”. These conditions resulted in reduced growth up to 34 days of age¹⁰, fecal-specific odors detectable by rats¹¹ and a fecal metabolomic signature highlighting persistent differences in adaptive response, energy metabolism and microbiota composition compared to optimal conditions¹². Interestingly, we reported that chicks provided with EO independently of their feed were able to self-select which EO they consumed in response to their postnatal experience and that this EO intake mitigated some detrimental effects on growth parameters induced by this negative postnatal experience¹⁰.

<footnote n="1">¹INRAE, Bordeaux INP, NutriNeuro, UMR 1286, University Bordeaux, 33076 Bordeaux, France. ²INRAE, Université de Tours, BOA, 37380 Nouzilly, France. ³INRAE, CNRS, IFCE, Université de Tours, PRC, Nouzilly, France. ⁴These authors contributed equally: Marie-Pierre Moisan and Laurence A. Guilloteau. ✉email: laurence.guilloteau@inrae.fr</footnote>
To complete these data, we examined whether the spontaneous intake of EO would mitigate the long-lasting consequences of negative postnatal experience on health parameters. More precisely, we hypothesized that the negative postnatal experience would alter blood physiological parameters such as metabolites, inflammation and oxidative stress biomarkers and reveal a specific transcriptomic signature in blood cells, but also that EO intake may restore these physiological parameters. We focused on gene expression profiling to assess health status because it is a non-invasive approach that is increasingly used to obtain a general, non-biased view of the biological pathways that are affected by a given condition. This approach has been used to assess the effects of early life stress in humans¹³ and also in animals¹⁴ including laying hens¹⁵ but not in broiler chickens.

Based on our previous study, three EO were chosen for their complementary properties to control microbial infections, reduce stress response, and improve digestive and immune system functions. Among the numerous biological activities of Cardamom EO, antioxidant, antispasmodic, anti-inflammatory, gastroprotective and antibacterial properties have been reported¹⁶⁻¹⁸. Marjoram EO demonstrates a variety of biological activities including a hepatoprotective role¹⁹, and lemon verbena EO has been shown to have analgesic, anti-inflammatory, sedative, and digestive properties²⁰.

Thus, the objectives of the present study were i) to characterize the immediate and long-lasting effects of a negative postnatal experience on blood physiological parameters and the long-lasting effects on whole gene expression in blood cells of fast-growing broiler chickens and ii) to analyze the effects of EO intake on blood metabolites, inflammatory and redox statuses and mRNA expression of a selection of genes in chicks exposed or not to this negative postnatal experience.

## Results

In this section, to facilitate the reading, the postnatal negative experience is called "delayed placement" and the group concerned as the "Delayed group".

**Immediate effects of delayed placement on physiological parameters of blood.** The chicks of the Delayed group demonstrated a significant immediate effect on metabolism shown by a decrease in plasma glucose and triglyceride concentrations (p < 0.0001 for both) and on the redox balance shown by a significant increase in thiobarbituric acid reactive substances (TBARS; p < 0.001) in the liver, and as a tendency there was a rise in blood total antioxidant status (TAS; p = 0.091) (Fig. 1). A significant sex effect (p = 0.009) and interaction (p = 0.011) was found for triglycerides only.

**Long-lasting effects of delayed placement and EO intake on physiological parameters of blood.** Metabolite plasma concentrations and redox balance differed significantly between males and females at D34. Triglyceride concentrations were higher in males than in females, whereas uric acid, TAS, ferric reducing ability of plasma (FRAP) and haptoglobin-like activities were higher in females than in males (Table 1).

Data are presented as mean ± SEM. Different letters correspond to significant differences (p < 0.05) between males and females (ANOVA).

At D34 in males, we observed no interaction between EO intake and the delayed placement after hatching and no delayed placement effects. However, EO intake had significant effects on two physiological parameters, namely uric acid and FRAP activity, which were increased with EO intake (Fig. 2A,B). At D34 in females, an interaction between EO intake and the delayed treatment was observed on haptoglobin-like activity (p = 0.034), but without significant differences between groups (Fig. 2D). No effect of the delayed treatment was observed, but triglyceride concentration was decreased with EO intake independently of the postnatal experience (Fig. 2C).

**Long-lasting effects of delayed placement on blood transcriptome in broiler chickens.** Microarray gene expression profiling was performed on blood cells from chickens at D34 in the Delayed and Control groups that had access to water only, to obtain a global biological signature of the long-lasting effects of the delayed placement. Microarrays consisted of 60,000 oligonucleotide probes, representing 12,349 annotated genes, i.e. 80% of the *Gallus gallus* genome. A principal component analysis of gene expression for each individual showed a strong sex effect (Fig. 3A), therefore sexes were analyzed separately thereafter. No genes were found to be expressed differentially when using an adjusted p-value < 0.05. A non-adjusted p-value of p < 0.005 was thus used (Fig. 3B Volcano plots). This analysis enabled 68 differentially regulated genes in males (40 up-regulated and 28 down-regulated in Delayed compared to Control group) and 83 differentially regulated genes to be detected in females (33 up-regulated and 50 down-regulated in Delayed compared to Control group). None of the differentially expressed genes were shared by males and females (Fig. 3C Venn diagram). The lists of differentially expressed genes are provided in Tables S1 (males) and S2 (females) in the supplementary files.

The Gene Ontology (GO) analysis of the differentially expressed genes for each sex indicated modification of broad biological processes (Fig. 4). Interestingly, in Delayed males the GO terms are related to an increase in response to stress and to external stimuli and in defense response and down-regulation of metabolic processes. In Delayed females, the GO terms suggest microtubule-based processes were modified, which would impact cell communication and signaling and decrease in mitochondrial gene expression.

To investigate the regulatory pathways that drive the expression of the differentially expressed genes we used the "single site" analysis of oPOSSUM software to detect over-represented conserved transcription factor binding sites in the set of differentially expressed genes for each sex. The transcription factors significantly over-represented in the Delayed compared to Control groups for each sex together with the corresponding target genes are shown in Fig. 5. In males, in the Delayed group no transcription factors were found significantly under-represented, but six transcription factors were over-represented, Foxd3, FOXI1, MEF2A, NFE2L2, Sox2 and TEAD1. In females, PLAG1 was the only transcription factor over-represented in the Delayed group. The
Glucose

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

D: *p* < 0.0001

Sex: *p* = 0.902

X: *p* = 0.766

Triglycerides

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

D: *p* < 0.0001

Sex: *p* = 0.009

X: *p* = 0.011

TBARS

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

D: *p* < 0.0001

Sex: *p* = 0.661

X: *p* = 0.748

TAS

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

D: *p* = 0.091

Sex: *p* = 0.490

X: *p* = 0.211

**Figure 1.** Effects of delayed placement and sex on metabolic and physiological parameters of plasma from one-day-old chicks: (A) glucose concentration; (B) triglyceride concentration (C) Thiobarbituric acid reactive substances (TBARS) and (D) Total Antioxidant Status (TAS). Bars represent the mean ± SEM of values measured in different groups. Two-way ANOVA results are shown for each parameter. D = delayed effect, X = interaction between delayed placement and sex effects.

other transcription factors NFIL3, ESR2, TAL1::TCF3 and Foxd3 were under-represented in the Delayed group. It should be noted that Foxd3 was over-represented in Delayed males, but under-represented in Delayed females.

**Long-lasting effects of EO intake on blood gene expression in broiler chickens after the delayed placement.**
To evaluate the effect of chick's EO intake on the blood transcriptome according to their postnatal experience, we performed quantitative PCR (qPCR) on a selection of genes based on the microarray results. We focused on the oxidative stress pathway identified in males (represented by NFE2L2 and MEF2A transcription factors) and on the PLAG1 pathway identified in females by selecting gene targets of either NFE2L2 or
<table><caption>**Table 1.** Metabolites and redox balance in blood of 34-day-old broiler chickens.</caption><thead><tr><td>Physiological parameters</td><td>Male</td><td>Female</td><td>p-value</td></tr></thead><tbody><tr><td>Glucose (g/L)</td><td>2.37 ± 0.04<sup>a</sup></td><td>2.36 ± 0.05<sup>a</sup></td><td>0.804</td></tr><tr><td>Triglycerides (mg/L)</td><td>0.80 ± 0.07<sup>a</sup></td><td>0.62 ± 0.03<sup>b</sup></td><td><strong>0.023</strong></td></tr><tr><td>Uric acid (mg/L)</td><td>56.57 ± 2.56<sup>a</sup></td><td>66.34 ± 3.77<sup>b</sup></td><td><strong>0.041</strong></td></tr><tr><td>TAS (mmol/L)</td><td>0.80 ± 0.02<sup>a</sup></td><td>0.96 ± 0.04<sup>b</sup></td><td><strong>0.001</strong></td></tr><tr><td>FRAP (µmol/L)</td><td>705.27 ± 22.53<sup>a</sup></td><td>800.45 ± 24.80<sup>b</sup></td><td><strong>0.007</strong></td></tr><tr><td>SOD (%)</td><td>77.79 ± 0.02<sup>a</sup></td><td>76.01 ± 0.02<sup>a</sup></td><td>0.439</td></tr><tr><td>TBARS (nm MDA/mL)</td><td>0.60 ± 0.02<sup>a</sup></td><td>0.60 ± 0.03<sup>a</sup></td><td>0.952</td></tr><tr><td>Haptoglobin-like (mg/mL)</td><td>0.18 ± 0.01<sup>a</sup></td><td>0.22 ± 0.01<sup>b</sup></td><td><strong>0.002</strong></td></tr></tbody></table>

MEF2A (n = 18) and targets of PLAG1 (n = 5) according to oPOSSUM. To complete the gene selection, eleven genes involved in oxidative stress or inflammation not tested on the microarray were also analyzed and called "oxidative stress" or "inflammation" genes respectively. Additionally, other genes which were significantly differentially expressed between Control and Delayed groups in males (n = 10) and females (n = 4) were selected and called "other responsive genes". The list of these 48 selected genes, their primer sequences used for quantitative real time PCR and their associated signaling pathways are presented in Table S3.

Quantitative PCR were conducted on blood cell cDNA from D34 chickens tested on the microarray analysis (Control and Delayed groups) and also those tested for EO or water intake (EO and W groups).

Significant interactions of gene expression were found between the delayed placement and the chickens' intake of EO. Among all the genes tested, a significant interaction or tendency was found in 19 out of 33 genes in males and in three out of 15 genes in females. In males provided with water only, a significant effect of the delayed placement on gene expression was confirmed in five genes involved in oxidative stress (C1QTNF7, SLC16A10, ANKH, HUS1 and SOST), but in only one gene (XDH) in water-only females (Fig. 6A,B). A significant effect (or tendency) of the delayed placement was also confirmed for the expression of four other genes in males and five genes in females, independently of the EO intake (Table 2). These genes were significantly differentially regulated on the microarray as "other responsive genes" (ASPM, BRD9 and PTPLB) in males, or involved in oxidative stress (SOD3) or inflammation (COX2 and PIT54) and in the PLAG1 pathway (MEX3A, TMEM125 and ZSWIM1) in females.

Interestingly, the Delayed group showed a clear favorable effect of EO intake for the SOST gene in males and females, for the SLC16A10 gene in males, and for the RARRES2 gene in females (Fig. 6A,B). Indeed, for these genes, EO intake in the Delayed group significantly restored the gene expression level observed in the Control group supplied with water only. Unexpectedly, the intake of EO significantly modified the gene expression in the Control group for 13 genes in males and for one gene (XDH) in females. These genes were involved in oxidative stress (C1QTNF7, KIAA1217, HUS1, RTN4, TBC1D12, ANKIB1, POLR1D and INOS) or qualified as "other responsive genes" in the microarray analysis (AP3M1, NDFIP1, OTUD3, TPPP, BTBD1 and XDH) (Table S3, Fig. 6A,B). A significant effect of EO intake was also observed in males on the expression of ASPM, COX2 and PIK3AP1 genes, whatever the postnatal experience (Table 2).

Discussion

This study presents immediate and long-lasting effects of a delayed placement of chickens considered as a negative postnatal experience, mimicking early-life stresses encountered in poultry production, and the potential beneficial effects of EO intake on blood physiological parameters and gene expression profiling of chicken blood cells. In a previous study, we showed that chicks spontaneously chose to consume EO or not, according to their postnatal experience when different EO were presented individually¹⁰.

In this study, we detected significant differences between sexes for several blood parameters. Triglyceride concentrations were higher in females than in males after hatching, but this was reversed at D34 indicating lower lipid metabolism in females at that age. At D34, higher anti-oxidative activity (uric acid, TAS and FRAP) and higher haptoglobin-like activity were observed in females than in males.

In both sexes, the negative postnatal experience resulted in an immediate decrease in blood glucose (typically around 2 g/L in control broiler chickens as commonly observed in this species²¹) and triglycerides associated with an increase in blood and liver lipid peroxidation (TBARS). This indicates that chicks exposed to these postnatal conditions adapted metabolically and suggests that they used the nutrients present in the internalized yolk sac²². Such a decrease in glycaemia was not observed in a previous study where chicks were subjected to early post-hatch delayed feeding²³ suggesting that in our study, this decrease could be the result of the combined effects of delayed feeding, variations in ambient temperatures and transport box shaking that may increase chick activity within the box. It could also be linked to a decrease in insulinemia as reported by Bigot et al.²³, not measured in the present experiment. This adaptive metabolic profile was present in fecal metabolome of chicks for at least for 13 days after the delayed placement, suggesting that the adaptive changes in metabolism due to the postnatal experience persisted¹², but no such effect was observed on blood physiological parameters at D34 in the present study.

An in-depth view of the long-lasting effects of the negative postnatal experience on the physiology of the chickens was investigated by analyzing gene expression profiles of whole blood at D34 using microarray and
Males

A

Uric Acid

D: *p* = 0.518

EO: *p* = 0.007

X: *p* = 0.456

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

B

FRAP

D: *p* = 0.462

EO: *p* = 0.026

X: *p* = 0.547

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

Females

C

Triglycerides

D: *p* = 0.898

EO: *p* = 0.045

X: *p* = 0.658

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

D

Haptoglobin-like

D: *p* = 0.892

EO: *p* = 0.781

X: *p* = 0.034

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

Control

Delayed

**Figure 2.** Effects of delayed placement and EO intake on physiological parameters of blood in 34-day-old male and female chickens: (A) urid acid concentration; (B) FRAP activity (C) triglyceride concentration and (D) haptoglobin-like. The bars represent the mean ± SEM of values measured in different groups. Two-way ANOVA results are shown for each parameter. D = delayed effect, EO = essential oil, X = interaction between delayed placement and EO effects.
<figure><img src="image_12.png" /><figcaption>C</figcaption></figure>

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

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

<figure><img src="image_15.png" /><figcaption>**Figure 3.** Microarray analysis of blood cells from 34-day-old chickens after delayed placement. Principal Components Analysis (PCA) conducted on microarray gene expression (normalized data), projection of individuals from each group and sex on the first two dimensions of the PCA (A). Volcano plots in male (blue) and female (orange) using -log10(non-adjusted p-value) on y-axis and log2(Delayed vs. Control chickens Fold Change) on x-axis. The red dotted line corresponds to a non-adjusted p-value = 0.005. Genes selected for further analyses are above this line (B). Venn diagram showing the number of differentially expressed genes between Delayed and Control groups for each sex and the lack of overlap between sexes (zero genes in common). UP = genes upregulated for Delayed chickens, DOWN = genes downregulated for Delayed chickens (C).</figcaption></figure>

qPCR. Very strong sex differences in the response to early external stimuli were observed at D34 in the blood gene expression profiles for which overall sex differences explain 28.7% of the variability in gene expression and from the microarray analysis none of the differentially expressed genes were shared between males and females.

Other studies have reported sex differences in chickens regarding stress response, including early stress, on behaviour, sex hormones and hypothalamic gene expression¹⁵,²⁴,²⁵. The results of the GO analysis were limited by the low number of genes differentially expressed between groups. Nevertheless, it is remarkable that the terms 'response to external stimuli', 'response to stress' and 'defense response' were enriched in Delayed males. The oPOSSUM analysis provided more precise data on the pathway involved. In particular, the negative postnatal experience of chicks induced changes in the expression of genes belonging to the oxidative stress pathway (as represented by the transcription factors NFE2L2 and MEF2A) in males and to a much lower extent in females. Using qPCR, seven gene targets of these two transcription factors were confirmed to be differentially regulated (BRD9, PTPLB, C1QTNF7, SLC16A10, ANKH, HUS1 and SOST) in Delayed males and three genes involved in oxidative stress response (XDH, SOD3 and PIT54) were also differentially regulated in the Delayed group in females. MEF2A is also known to influence muscle development and has been described in chicken²⁶. Thus the MEF2A pathway may be activated to counter the observed decrease in *P. major* muscle yield observed in the Delayed chick group¹⁰.

Besides oxidative stress, a number of the differentially regulated genes in males have interesting properties, such as the up-regulated genes related to bone metabolism: ANKH has a role in osteogenic differentiation of bone marrow stromal cells and is a vascular calcification inhibitor²⁷,²⁸ and SOST is a negative regulator of bone growth²⁹,³⁰. Moreover, the transcription factor TEAD1 over-represented in the Delayed males supports an effect
<figure><img src="image_16.png" /><figcaption>**Figure 4.** Heatmaps displaying 68 and 83 differentially expressed genes among Delayed and Control groups in males and females respectively. Enriched GO terms (non-adjusted p-value < 0.05) and genes are shown on the right for each group. The p-values of GO terms were determined on the WebGestalt website.</figcaption></figure>

of delayed placement on bone metabolism by regulating organ size and skeletal muscle mass³¹. In line with this, SLC16A10 gene over-represented in the Delayed males is a transporter gene belonging to a family of transporters involved in energy metabolism, including thyroid hormone metabolism and possibly skeletal development³².

In females, the PLAG1 pathway was the only pathway over-represented following the negative postnatal experience. PLAG1 is a regulator of growth and reproduction³³ and associated with variability in egg production in laying hens³⁴. Through qPCR, three genes targets of PLAG1 (MEX3A, TMEM215, ZSWIM1) were shown to be up-regulated in the female Delayed group. This is interesting in view of the reduced growth observed in the chicks of the Delayed group. Some other significant pathways were under-represented following the negative postnatal experience and involved gene targets of NFIL3, Foxd3 and ESR2. NFIL3 is involved in the circadian rhythm in the chick's pineal gland and is essential as a regulator of the immune response³⁵. ESR2 plays a central role in folliculogenesis and therefore in reproduction³⁶. Foxd3 is important in embryogenesis being involved in regulating the lineage choice between neural crest-derived glial cells³⁷ and in immunity by regulating IL-10-positive B cells³⁸. Foxd3 under-represented in Delayed females was in contrast over-represented in Delayed males. This may be related to the higher immune resistance of female compared to male chicken as reported in a large study on mortality following infectious ³⁹.

Remarkably, the spontaneous intake of these EO by chicks for 12 days after hatching had long-lasting beneficial effects on the expression of several of these genes modified by the negative postnatal experience. The expression of SOST (both in males and females), SLC16A10 (only in males) and RARRES2 (only in females) genes was restored to the expression observed in the Control group when chickens had access to EO. Interestingly, RARRES2 is an adipokine that plays a major role in the regulation of metabolic and reproductive processes recently studied in the chicken⁴⁰, in line with the list of PLAG1 gene targets, and is also a chemokine for leucocyte populations and participates in antimicrobial activity.

In our study, the beneficial role of EO intake after the negative postnatal experience could be to regulate the growth and bone metabolism in male and female chicken, and also the energy metabolism and reproduction in females.

Surprisingly, EO intake by control chickens modified the expression of 17 genes, 16 in males and one in females. Most of them are involved in oxidative stress, immune regulation and microtubule network integrity.
<figure><img src="image_17.png" /><figcaption>**Figure 5.** Circos plots (GOPlot R package) of over-represented conserved transcription factors (on the right side) and corresponding target genes (left side) enriched in Delayed compared to Control chickens for each sex. Up-regulated and down-regulated genes are indicated by a red and blue square, respectively. LogFC=log2 Fold Change. Significant over- or under-represented transcription factors were detected using oPOSSUM software.</figcaption></figure>

The majority of the genes studied and up-regulated in the Control group were different to those up-regulated in the Delayed group suggesting that the reasons why chicks ingested EO differed between these two placement conditions.

To our knowledge, this is the first study in chickens to show long-lasting effects of spontaneous intake of EO on gene expression, especially changing some blood expression patterns affected by the chick's post-hatch experience. Further studies are needed to improve the understanding of these effects on blood transcriptome, in tissues of the liver and gut for example and on gut microbiota which could be considerably modified by the antimicrobial properties of EO. The self-medication behaviour in domestic animals warrants further investigation, especially regarding infectious challenges. It could be greatly beneficial for animals like chickens raised in large groups to be allowed to manage their health and welfare individually through free access to EO, hence limiting the use of antibiotics in the context of the One Health concept.

## Methods

**Essential oils (EO).** Cardamom (*Elettaria cardamomum*) (1480CQ, batch S12A, Herbes et Traditions, Comines, France), marjoram (*Origanum majorana*) CT thujanol (2507CQ, batch S12D, Herbes et Traditions), and lemon verbena (*Lippia citriodora*) (FLE094, batch H181013MA, Florihana, Caussols, France) were chosen based on the properties of their major components and on our previous study which showed spontaneous intake and long-term effects after a negative postnatal experience in chicks¹⁰. Each EO was diluted in water (0.001%), mixed and shaken vigorously before being made available in a drinking bottle. The main components obtained by gas chromatography coupled to mass spectrometry for each EO are listed in Table 3.

**Experimental design.** A model of postnatal negative experience was previously developed to analyze the consequences of this experience over the whole growing period of broiler chickens¹⁰. This model involved reproducing suboptimal transportation conditions for chicks after hatching. The chicks (Hubbard Classic, Quintin, France) were either placed immediately in pens in the rearing facility after their removal from the incubator (Control group, C, n=192) or they were removed and subjected to 24 h of negative experience before being placed in pens (Delayed group, D, n=192). This delayed group was deprived of both food and water and subjected to irregular movement and various room temperatures¹⁰. Half of the chicks (six pens each for C and D groups) had access to water only in four bottles (W-C and W-D groups). Besides feed and water supplies, the other half of the chicks had ad libitum access to the three EO (EO-C or EO-D) from D1 until D12 post-hatching. One bottle containing water, and the three others each containing one of the EO, were placed in each pen. Chicks were reared at the Experimental Poultry Facility (PEAT, INRAE, 2018, 37,380 Nouzilly, France, https://doi.org/10.15454/1.5572326250887292E12) under standard temperature and light conditions with ad libitum access to water and with a wire mesh platform and a perch for environmental enrichment. At D13, the chickens
A

Males

<figure><img src="image_18.png" /><figcaption>**Figure 6.** Effects of EO intake on blood-cell gene expression in 34-day-old male (A) and female (B) chickens after the delayed placement (qPCR analyses). Bars represent the mean ± SEM of values measured in different groups. Different letters between bar-plots indicate significant or tendency differences (p < 0.05 or p < 0.10) between groups detected by post-hoc analyses following two-way ANOVA.</figcaption></figure>
A

Males

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

B

Females

<figure><img src="image_20.png" /><figcaption>**Figure 6.** (continued)</figcaption></figure>

were transferred to another poultry building for the growth phase until D34. They had ad libitum access to feed without anticoccidial drugs. They were fed with a standard starting diet (metabolizable energy = 12.8 MJ/kg, crude protein = 22%) until D19 and then a rearing diet from D19 to D34.

Physiological analyses were performed at D1 (first experiment), and D34 post hatching (second experiment) on plasma and liver samples from 12 chickens (6 males and 6 females) per condition (2 conditions × 2 experiments × 12 chickens = 48 chickens). Transcriptomic analysis was performed at D34 on blood samples from the same chickens as those sampled for physiological parameters. All procedures used in these experiments were approved by the local ethics committee (Comité d'Ethique en Expérimentation Animale Val de Loire, Tours, France; permission no 01730.02 and 2015070815347034v2 (APAFIS#1082) and carried out in accordance with current European legislation (EU Directive 2010/63/EU).

**Physiological parameters.** Different metabolic, antioxidant and oxidative statuses and inflammation parameters were measured at D1 and D34 after hatching. Commercial kits (THERMO FISHER DIAGNOSTICS SAS, Courtaboeuf, France) were used to determine plasma glucose (g/L) (MG981780), uric acid (mg/L) (MG981788) and triglycerides (mg/L) (MG981786). Total plasma antioxidant activity was determined through
<table><caption>**Table 2.** Placement and EO effects on blood RNA expression in 34-day-old chickens after a negative postnatal experience. q PCR data (mean ± SEM) were only shown for significant gene expression without significant interaction between the two factors (ANOVA). *p* < 0.05 considered significant and † 0.05 < *p* < 0.1 considered to be a tendency.</caption><thead><tr><th rowspan="2">Gene</th><th colspan="3">Placement</th><th colspan="3">Treatment</th><th rowspan="2">Placement<br>x Treatment<br>interaction p-value</th></tr><tr><th>Control</th><th>Delayed</th><th>Placement p-value</th><th>Water</th><th>EO</th><th>Treatment p-value</th></tr></thead><tbody><tr><td colspan="8"><strong>Males</strong></td></tr><tr><td>BRD9</td><td>1.26 ± 0.14</td><td>0.89 ± 0.11</td><td>0.050<sup>†</sup></td><td>0.96 ± 0.12</td><td>1.19 ± 0.14</td><td>0.195</td><td>0.388</td></tr><tr><td>PTPLB</td><td>1.13 ± 0.18</td><td>0.64 ± 0.13</td><td>0.055<sup>†</sup></td><td>1.02 ± 0.13</td><td>0.81 ± 0.19</td><td>0.493</td><td>0.364</td></tr><tr><td>ASPM</td><td>1.50 ± 0.19</td><td>1.00 ± 0.11</td><td><strong>0.023*</strong></td><td>1.05 ± 0.14</td><td>1.50 ± 0.18</td><td><strong>0.038*</strong></td><td>0.174</td></tr><tr><td>COX2</td><td>0.82 ± 0.12</td><td>0.52 ± 0.11</td><td>0.052<sup>†</sup></td><td>0.97 ± 0.08</td><td>0.44 ± 0.10</td><td><strong>0.001*</strong></td><td>0.580</td></tr><tr><td>PIK3AP1</td><td>1.34 ± 0.15</td><td>1.11 ± 0.09</td><td>0.135</td><td>1.01 ± 0.08</td><td>1.44 ± 0.14</td><td><strong>0.011*</strong></td><td>0.255</td></tr><tr><td>DUS4L</td><td>1.46 ± 0.18</td><td>1.34 ± 0.16</td><td>0.633</td><td>1.18 ± 0.19</td><td>1.63 ± 0.13</td><td>0.086<sup>†</sup></td><td>0.615</td></tr><tr><td colspan="8"><strong>Females</strong></td></tr><tr><td>MEX3A</td><td>3.10 ± 0.86</td><td>12.30 ± 3.95</td><td><strong>0.038*</strong></td><td>7.74 ± 2.23</td><td>7.65 ± 4.01</td><td>0.673</td><td>0.336</td></tr><tr><td>TMEM125</td><td>1.05 ± 0.31</td><td>3.63 ± 0.92</td><td><strong>0.010*</strong></td><td>2.05 ± 0.57</td><td>2.81 ± 1.06</td><td>0.225</td><td>0.073</td></tr><tr><td>ZSWIM1</td><td>1.51 ± 0.16</td><td>2.34 ± 0.35</td><td><strong>0.035*</strong></td><td>1.67 ± 0.19</td><td>2.19 ± 0.37</td><td>0.127</td><td>0.668</td></tr><tr><td>PIT54</td><td>1.18 ± 0.63</td><td>4.58 ± 1.11</td><td>0.052<sup>†</sup></td><td>4.23 ± 1.19</td><td>1.91 ± 0.91</td><td>0.268</td><td>0.911</td></tr><tr><td>SOD3</td><td>2.70 ± 1.25</td><td>10.29 ± 3.01</td><td>0.080<sup>†</sup></td><td>8.61 ± 2.87</td><td>5.10 ± 1.57</td><td>0.417</td><td>0.764</td></tr></tbody></table>

<table><caption>**Table 3.** Essential oils composition (from [10]).<sup>a</sup>1480CQ, batch S12A, Herbes et Traditions, Comines, France.<sup>b</sup>2507CQ, batch S12D, Herbes et Traditions.<sup>c</sup>FLE094, batch H181013MA, Florihana, Caussols, France.</caption><thead><tr><th rowspan="2">Compound</th><th colspan="2">Cardamom<br>Elettaria cardamomum</th><th colspan="2">Marjoram<br>Origanum majorana<br>CT thujanol</th><th colspan="2">Verbena<br>Lippia citriodora</th></tr><tr><th>Specification (%)</th><th>Relative content (%)<sup>a</sup></th><th>Specification (%)</th><th>Relative content (%)<sup>b</sup></th><th>Specification (%)</th><th>Relative content (%)<sup>c</sup></th></tr></thead><tbody><tr><td>Monoterpenes</td><td>6—12</td><td>13</td><td>30</td><td>40</td><td>5—15</td><td>29</td></tr><tr><td>Sesquiterpenes</td><td></td><td></td><td></td><td>3</td><td>18—26</td><td>24.5</td></tr><tr><td>Monoterpenols</td><td>3—6</td><td>5</td><td>40—50 (20<br>thujanol)</td><td>50 (25 thujanol)</td><td>3—15</td><td>2</td></tr><tr><td>Esters</td><td>39—51</td><td>36</td><td></td><td>2</td><td></td><td></td></tr><tr><td>Oxides</td><td>27—35</td><td>34</td><td></td><td></td><td>&lt;7</td><td>5</td></tr><tr><td>Aldehydes</td><td></td><td></td><td></td><td></td><td>20—40</td><td>24</td></tr></tbody></table>

Total antioxidant status (TAS) (mmol/L) (NX 2332, RANDOX LABORATOIRES, Roissy, France). Protocols were used in accordance with supplier instructions and adapted to the automated Thermo Scientific Arena 20XT photometric analyzer (THERMO FISHER DIAGNOSTICS). The Ferric reducing/antioxidant power (FRAP) was determined as described by Benzie and Strain<sup>41</sup> and results expressed as µmol Trolox/L. Plasma superoxide dismutase (SOD) activity was measured with a commercial kit (19,160, Sigma-Aldrich Chemie GmbH, Buchs, Switzerland), using a microplate reader (TECAN infinite 200, Tecan Group Ltd, Männedorf, Switzerland) and expressed as an inhibition activity (%). Haptoglobin-like activity described as PIT54 in chicken<sup>42</sup> was evaluated in plasma using a colorimetric assay (TP-801, Tridelta Development Ltd, Maynooth, Ireland) measuring inhibition of hemoglobin peroxidase activity and expressed as mg/mL. Lipid peroxidation was determined using spectrophotometry of thiobarbituric acid reacting substances (TBARS) as previously described by Lynch and Frei for the liver<sup>43</sup>, and adapted from Lin et al.<sup>44</sup> for plasma TBARS.

**RNA extraction and gene expression analyses.** RNA extraction. Blood was collected from chickens at D34 using 5 ml EDTA vacutainer tubes. Blood (v = 100 µL) was immediately suspended in 1 mL of Invitrogen TRIzol reagent (FISHER SCIENTIFIC, Illkirch-Graffenstaden, France) and vigorously shaken for 5 min on ice. The samples were kept at -80 °C until RNA extraction. Total RNA was extracted using TRIzol (FISHER SCIENTIFIC) following the manufacturer's protocol modified in accordance with Desert et al.<sup>45</sup>. Ten µl of acetic acid 5 N were added with chloroform to reduce DNA contamination. The quality of total RNA was assessed using RNA Nano chips on a 2100 Bioanalyzer Instrument (Agilent, Waldbronn, Germany). All the samples had an RNA Integrity Number (RIN) score of > 8.0.

**Microarrays.** Gene expression profiles were performed at the GeT-TRiX facility (GenoToul, Génopole Toulouse, Toulouse, Midi-Pyrénées) using Agilent SurePrint G3 gallus_exp_microarray_SLagarrigue_8×60k_
V2_july2012 microarrays (8×60 K, design 042,004) following the manufacturer's instructions. For each sample, Cyanine-3 (Cy3) labeled cRNA was prepared from 200 ng of total RNA using the One-Color Quick Amp Labeling kit (Agilent) in line with the manufacturer's instructions, followed by RNA clean-up using Agencourt RNAClean XP (Agencourt Bioscience Corporation, Beverly, Massachusetts). Dye incorporation and cRNA yield were checked using Dropsense 96 UV/VIS droplet reader (Trinean, Ghent, Belgium). Aliquots of 600 ng of Cy3-labeled cRNA were hybridized on the microarray slides following the manufacturer's instructions. Immediately after washing, the slides were scanned on an Agilent G2505C Microarray Scanner using Agilent Scan Control A.8.5.1 software and the fluorescence signal was extracted using Agilent Feature Extraction software v10.10.1.1 with default parameters. Microarray data and the experimental details are available in the NCBI's Gene Expression Omnibus⁴⁶ and are accessible through GEO Series accession number GSE102358 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE102358).

**Real-time quantitative PCR (qPCR).** To avoid genomic DNA amplification, primer pairs were designed in two different exons (thus spanning an intron) using the Primer Express software (Applied Biosystems, THERMO FISHER DIAGNOSTICS). The sequences of primers used are provided in Table S3. The specificity of the PCR reaction was validated according to MIQE (Minimum Information for publication of Quantitative real time PCR Experiments) guidelines⁴⁷. An aliquot of 2,000 ng of total RNA was reverse-transcribed in cDNA with Superscript III (Invitrogen, THERMO FISHER DIAGNOSTICS) and random hexamers in accordance with the manufacturer's protocol. High throughput real time quantitative PCR was performed on the Biomark HD System (Fluidigm, Les Ullis, France) following the manufacturer's protocol. The chip was placed into the IFC Controller, where 6.3 nl of Sample Mix and 0.7 nl of Assay Mix were mixed. Real-time PCR was performed on the Biomark System as follows: Thermal Mix at 50 °C, 2 min; 70 °C, 30 min; 25 °C, 10 min, UNG at 50 °C, 2 min, Hot Start at 95 °C, 10 min, PCR Cycle of 35 cycles at 95 °C, 15 s; 60 °C, 60 s and Melting curves (from 60 °C to 95 °C). Results were analyzed using Fluidigm Real-Time PCR Analysis software v4.1.3. (Fluidigm) to control specific amplification for each primer, then the raw results of the qPCR were analyzed using GenEx software (MultiD analyses AB, Göteborg, Sweden) in order to choose the best reference gene to normalize mRNA expression and measure the relative expression of each gene between groups. HPRT was found to be the best reference gene in this experiment and was thus used for the normalization of gene expression.

**Statistical analyses.** The effects of a negative postnatal experience and EO intake and their interaction on physiological responses and qPCR analyses were determined using 2-way analysis of variance (ANOVA) after having checked the normality of data distribution. When there was an interaction between variables, a Fisher (LSD) test was used to determine the statistical significance of the difference. Differences were considered significant when p-values < 0.05 and a tendency for p < 0.10. Analyses were performed using XLSTAT software (version 2015, Addinsoft, Paris, France).

Microarray data were analyzed using R⁴⁸ and Bioconductor packages (www.bioconductor.org, v 3.0)⁴⁹ as described in GEO accession GSE102358. Briefly, raw data (median of pixel intensity) were filtered, log2 transformed, corrected for bath effects (washing and labeling serials) and normalized using the quantile method⁵⁰. Principal Component Analysis (PCA) of the normalized microarray data was used to create a factor map of individual global expression using R package FactoMineR⁵¹.

A model was fitted using the limma lmFit function⁵² considering sex as a blocking factor. A correction for multiple testing was then applied using the Benjamini-Hochberg procedure (BH)⁵³ for the false discovery rate (FDR). Probes with FDR ≤ 0.05 were considered to be differentially expressed between conditions. For bioinformatic analyses we selected probes that displayed a non-adjusted p value of p < 0.005 as too few probes were found with a FDR ≤ 0.05. All differentially expressed genes, screened on conditions of showing a ≥ 20% difference in mean expression levels between samples from delayed and control conditions, and a non-adjusted p-value < 0.005, were subjected to GO classification [GO-BP (biological process)] to assign these genes to relevant GO Terms using WebGestalt website⁵⁴. oPOSSUM website was used to determine the over-representation of transcription factor binding sites within the differentially expressed genes⁵⁵.

Received: 20 January 2020; Accepted: 10 November 2020

Published online: 26 November 2020

**References**

1. Ericsson, M. & Jensen, P. Domestication and ontogeny effects on the stress response in young chickens (*Gallus gallus*). *Sci. Rep.-Uk.* https://doi.org/10.1038/srep35818 (2016).

2. Bergoug, H. *et al*. Effect of transportation duration of 1-day-old chicks on postplacement production performances and pododermatitis of broilers up to slaughter age. *Poult. Sci.* **92**, 3300–3309. https://doi.org/10.3382/ps.2013-03118 (2013).

3. Bakkali, F., Averbeck, S., Averbeck, D. & Waomar, M. Biological effects of essential oils: a review. *Food Chem. Toxicol.* **46**, 446–475. https://doi.org/10.1016/j.fct.2007.09.106 (2008).

4. Adorjan, B. & Buchbauer, G. Biological properties of essential oils: an updated review. *Flavour. Frag. J.* **25**, 407–426. https://doi.org/10.1002/ffj.2024 (2010).

5. Dobetsberger, C. & Buchbauer, G. Actions of essential oils on the central nervous system: an updated review. *Flavour. Frag. J.* **26**, 300–316. https://doi.org/10.1002/ffj.2045 (2011).

6. de Cassia da Silveira, E. S. R., Andrade, L. N., Dos Reis Barreto de Oliveira, R. & de Sousa, D. P. A review on anti-inflammatory activity of phenylpropanoids found in essential oils. *Molecules* **19**, 1459–1480. https://doi.org/10.3390/molecules19021459 (2014).

7. Nazzaro F, F. F., De Martino L, Coppola R, De Feo V. Effects of essential oils on pathogenic bacteria. *Pharmaceuticals* **6**, 1451–1474 (2013).
8. Brenes, A. & Roura, E. Essential oils in poultry nutrition: main effects and modes of action. *Anim. Feed Sci. Technol.* **158**, 1–14. https://doi.org/10.1016/j.anifeedsci.2010.03.007 (2010).

9. Adaszynska-Skwirzynska, M. & Szczersbinska, D. Use of essential oils in broiler chicken production a review. *Ann. Anim. Sci.* **17**, 317–335. https://doi.org/10.1515/aoas-2016-0046 (2017).

10. Guilloteau, L. A., Collin, A., Koch, A. & Letterrier, C. Spontaneous intake and long-term effects of essential oils after a negative postnatal experience in chicks. *Front. Vet. Sci.* **6**, 72. https://doi.org/10.3389/fvets.2019.00072 (2019).

11. Bombail, V. *et al*. In search of stress odours across species: behavioural responses of rats to faeces from chickens and rats subjected to various types of stressful events. *Appl. Anim. Behav. Sci.* **205**, 216–226. https://doi.org/10.1016/j.applanim.2017.10.013 (2018).

12. Beaucercq, S. *et al*. A multiplatform metabolomic approach to characterize fecal signatures of negative postnatal events in chicks: a pilot study. *J. Anim. Sci. Biotechnol.* **10**, 21. https://doi.org/10.1186/s40104-019-0335-8 (2019).

13. Slavich, G. M. & Cole, S. W. The emerging field of human social genomics. *Clin. Psychol. Sci.* **1**, 331–348. https://doi.org/10.1177/2167702613478594 (2013).

14. Lansade, L. *et al*. Behavioral and Transcriptomic fingerprints of an enriched environment in horses (*Equus caballus*). *PLoS ONE* **9**, e114384. https://doi.org/10.1371/journal.pone.0114384 (2014).

15. Hedlund, L., Whittle, R. & Jensen, P. Effects of commercial hatchery processing on short- and long-term stress responses in laying hens. *Sci. Rep.* **9**, 2367. https://doi.org/10.1038/s41598-019-38817-y (2019).

16. al-Zuhair, H., el-Sayeh, B., Ameen, H. A. & al-Shoora, H. Pharmacological studies of cardamom oil in animals. *Pharmacol. Res.* **34**, 79–82. https://doi.org/10.1006/phrs.1996.0067 (1996).

17. Jamal, A., Javed, K., Aslam, M. & Jafri, M. A. Gastroprotective effect of cardamom, Elettaria cardamomum Maton. fruits in rats. *J. Ethnopharmacol.* **103**, 149–153. https://doi.org/10.1016/j.jep.2005.07.016 (2006).

18. Ashokkumar, K., Murugan, M., Dhanya, M. K. & Warkentin, T. D. Botany, traditional uses, phytochemistry and biological activities of cardamom [Elettaria cardamomum (L.) Maton]: a critical review. *J. Ethnopharmacol.* **246**, 112244. https://doi.org/10.1016/j.jep.2019.112244 (2020).

19. Bina, F. & Rahimi, R. Sweet marjoram: a review of ethnopharmacology, phytochemistry, and biological activities. *J. Evid. Based Complement. Altern. Med.* **22**, 175–185. https://doi.org/10.1177/2156587216650793 (2017).

20. Pascual ME, S. K., Carretero E, Sanchez Mata D, Villar A. Lippia: traditional uses, chemistry and pharmacology: a review. *J. Ethnopharmacol.* **76**, 201–214 (2001).

21. Simon, J. Chicken as a useful species for the comprehension of insulin action. *Crit. Rev. Poult. Biol.* **2**, 121–148 (1989).

22. Yadgary, L., Cahaner, A., Kedar, O. & Uni, Z. Yolk sac nutrient composition and fat uptake in late-term embryos in eggs from young and old broiler breeder hens. *Poult. Sci.* **89**, 2441–2452. https://doi.org/10.3382/ps.2010-00681 (2010).

23. Bigot, K., Mignon-Grasteau, S., Picard, M. & Tesseraud, S. Effects of delayed feed intake on body, intestine, and muscle development in neonate broilers. *Poult. Sci.* **82**, 781–788 (2003).

24. Elfwing, M. *et al*. Early stress causes sex-specific, life-long changes in behaviour, levels of gonadal hormones, and gene expression in chickens. *PLoS ONE* **10**, e0125808. https://doi.org/10.1371/journal.pone.0125808 (2015).

25. Madison, F. N., Jurkevich, A. & Kuenzel, W. J. Sex differences in plasma corticosterone release in undisturbed chickens (*Gallus gallus*) in response to arginine vasotocin and corticotropin releasing hormone. *Gen. Comp. Endocrinol.* **155**, 566–573. https://doi.org/10.1016/j.ygcen.2007.08.014 (2008).

26. Zhou, Y. *et al*. Polymorphism of chicken myocyte-specific enhancer-binding factor 2A gene and its association with chicken carcass traits. *Mol. Biol. Rep.* **37**, 587–594. https://doi.org/10.1007/s11033-009-9838-2 (2010).

27. Minashima, T., Quirno, M., Lee, Y. J. & Kirsch, T. The role of the progressive ankylosis protein (ANK) in adipogenic/osteogenic fate decision of precursor cells. *Bone* **98**, 38–46. https://doi.org/10.1016/j.bone.2017.03.003 (2017).

28. Chen, B. *et al*. Wnt1 inhibits vascular smooth muscle cell calcification by promoting ANKH expression. *J. Mol. Cell. Cardiol.* **135**, 10–21 (2019).

29. Loots, G. G. *et al*. Genomic deletion of a long-range bone enhancer misregulates sclerostin in Van Buchem disease. *Genome Res.* **15**, 928–935. https://doi.org/10.1101/gr.3437105 (2005).

30. Sebastian, A. & Loots, G. G. Transcriptional control of Sost in bone. *Bone* **96**, 76–84. https://doi.org/10.1016/j.bone.2016.10.009 (2017).

31. Watt, K. I., Goodman, C. A., Hornberger, T. A. & Gregorevic, P. The Hippo signaling pathway in the regulation of skeletal muscle mass and function. *Exerc. Sport Sci. Rev.* **46**, 92–96. https://doi.org/10.1249/JES.0000000000000142 (2018).

32. Mariotta, L. *et al*. T-type amino acid transporter TAT1 (Slc16a10) is essential for extracellular aromatic amino acid homeostasis control. *J. Physiol.* **590**, 6413–6424. https://doi.org/10.1113/jphysiol.2012.239574 (2012).

33. Juma, A. R., Damdimopoulou, P. E., Grommen, S. V., Van de Ven, W. J. & De Groef, B. Emerging role of PLAG1 as a regulator of growth and reproduction. *J. Endocrinol.* **228**, R45–R56. https://doi.org/10.1530/JOE-15-0449 (2016).

34. Chen, C. *et al*. Laying traits and underlying transcripts, expressed in the hypothalamus and pituitary gland, that were associated with egg production variability in chickens. *Theriogenology* **68**, 1305–1315 (2007).

35. Male, V., Nisoli, I., Gascoyne, D. & Brady, H. E4BP4: an unexpected player in the immune response. *Trends Immunol.* **33**, 98–102 (2012).

36. Niu, X. *et al*. Sequence variations in estrogen receptor 1 and 2 genes and their association with egg production traits in Chinese Dagu chickens. *J. Vet. Med. Sci.* **79**, 927–934. https://doi.org/10.1292/jvms.17-0014 (2017).

37. Thomas, A. & Erickson, C. FOXD3 regulates the lineage switch between neural crest-derived glial cells and pigment cells by repressing MITF through a non-canonical mechanism. *Development* **136**, 1849–1858 (2009).

38. Zhang, Y. *et al*. Foxd3 suppresses interleukin-10 expression in B cells. *Immunology* **150**, 478–488. https://doi.org/10.1111/imm.12701 (2017).

39. Leitner, G., Heller, E. D. & Friedman, A. Sex-related differences in immune response and survival rate of broiler chickens. *Vet. Immunol. Immunopathol.* **21**, 249–260. https://doi.org/10.1016/0165-2427(89)90035-4 (1989).

40. Mellouk, N. *et al*. Chicken is a useful model to investigate the role of Adipokines in metabolic and reproductive diseases. *Int. J. Endocrinol.* **2018**, 4579734. https://doi.org/10.1155/2018/4579734 (2018).

41. Benzie, I. F. & Strain, J. J. The ferric reducing ability of plasma (FRAP) as a measure of "antioxidant power": the FRAP assay. *Anal. Biochem.* **239**, 70–76. https://doi.org/10.1006/abio.1996.0292 (1996).

42. Wicher, K. B. & Fries, E. Haptoglobin, a hemoglobin-binding plasma protein, is present in bony fish and mammals but not in frog and chicken. *Proc. Natl. Acad. Sci. USA* **103**, 4168–4173. https://doi.org/10.1073/pnas.0508723103 (2006).

43. Lynch, S. M. & Frei, B. Mechanisms of copper- and iron-dependent oxidative modification of human low density lipoprotein. *J. Lipid Res.* **34**, 1745–1753 (1993).

44. Lin, H., Decuypere, E. & Buyse, J. Oxidative stress induced by corticosterone administration in broiler chickens (*Gallus gal- lus* domesticus) 2. Short-term effect. *Comp. Biochem. Physiol. B Biochem. Mol. Biol.* **139**, 745–751. https://doi.org/10.1016/j. cbpc.2004.09.014 (2004).

45. Desert, C. *et al*. Transcriptomes of whole blood and PBMC in chickens. *Comp. Biochem. Physiol. Part D Genomics Proteomics* **20**, 1–9. https://doi.org/10.1016/j.cbd.2016.06.008 (2016).

46. Edgar, R., Domrachev, M. & Lash, A. E. Gene expression omnibus: NCBI gene expression and hybridization array data repository. *Nucleic Acids Res.* **30**, 207–210 (2002).
47. Bustin, S. A. et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. *Clin. Chem.* **55**, 611–622. https://doi.org/10.1373/clinchem.2008.112797 (2009).

48. R Core Team (2008). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna. URL http://www.R-project.org/.

49. Gentleman, R. C. et al. Bioconductor: open software development for computational biology and bioinformatics. *Genome Biol.* **5**, R80. https://doi.org/10.1186/gb-2004-5-10-r80 (2004).

50. Bolstad, B. M., Irizarry, R. A., Astrand, M. & Speed, T. P. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. *Bioinformatics* **19**, 185–193 (2003).

51. Lé, S., Josse, J. & Husson, F. FactoMineR: an R package for multivariate analysis. *J. Stat. Softw.* **25**, 1–18 (2008).

52. Smyth, G. K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.* https://doi.org/10.2202/1544-6115.1027 (2004).

53. Benjamini, Y. H. Y. Controlling the false discovery rate: a pratical and powerful approach to multiple testing. *J. R. Stat. Soc.* **57**, 289–300 (1995).

54. Liao, Y., Wang, J., Jaehnig, E., Shi, Z. & Zhang, B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. *Nucleic Acids Res.* **47**, W199–W205. https://doi.org/10.1093/nar/gkz401 (2019).

55. Kwon, A. T., Arenillas, D. J., Worsley Hunt, R. & Wasserman, W. W. oPOSSUM-3: advanced analysis of regulatory motif over-representation across genes or ChIP-Seq datasets. *G3 (Bethesda)* **2**, 987–1002. https://doi.org/10.1534/g3.112.003202 (2012).

## Acknowledgements

We would like to thank the staff of the Experimental Poultry Facility (PEAT, INRAE, 2018, 37380 Nouzilly, France. https://doi.org/10.15454/1.5572326250887292E12) for producing the chickens and for their assistance in experiments, S. Crochet, E. Cailleau-Audouin and P. Chartrin from the Avian Biology and Poultry Science Unit (UMR BOA, INRAE, Université de Tours, 37380 Nouzilly, France) and P. Constantin (UMR PRC, INRAE, CNRS, IFCE, Université de Tours, 37380 Nouzilly, France) for their help in experiments and biological analyses, and A. Koch and I. Colina Moreno, the students involved in the study. We also thank S. Lagarrigue and C. Désert (Pegase, Agrocampus Ouest, INRAE, 35590 Saint-Gilles, France) for their advice in blood RNA extraction method, M. Derval and C. Ingraham for their expertise in aromatherapy and self-medication practices with livestock. Lastly, we thank C. Naylies and Y. Lippi for their help in acquiring microarray fingerprints and microarray data analysis carried out at GeT-TRiX Genopole Toulouse Midi-Pyrénées facility. The English was reviewed by Susan Edrich (Inter-connect LC).

## Author contributions

L.G. and C.L. conceived and designed the experiments; L.G., C.L. and A.C. performed the experiments; L.G., M.-P.M., A.F., J.-C.H. analyzed the data; L.G. and M.-P.M. wrote the paper. A.C. and C.L. contributed to the manuscript revision. All the authors read and approved the submitted version.

## Funding

This work was supported by a grant from the Integrated Management of Animal Health Metaprogram of INRAE for the GISA-WHELP project (www.gisa.inra.fr/en).

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-77732-5.

**Correspondence** and requests for materials should be addressed to L.A.G.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020