# GWAS analyses reveal QTL in egg layers that differ in response to diet differences

Hélène Romé, Amandine Varenne, Frédéric Herault, Hervé Chapuis, Christophe Alleno, Patrice Dehais, Alain Vignal, Thierry Burlot, Pascale Le Roy

► To cite this version:

Hélène Romé, Amandine Varenne, Frédéric Herault, Hervé Chapuis, Christophe Alleno, et al.. GWAS analyses reveal QTL in egg layers that differ in response to diet differences. Genetics Selection Evolution, 2015, 47 (1), n.p. 10.1186/s12711-015-0160-2 . hal-01230151

HAL Id: hal-01230151

https://hal.science/hal-01230151

Submitted on 17 Nov 2015

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>

# GWAS analyses reveal QTL in egg layers that differ in response to diet differences

Hélène Romé¹,<fnref n="2" />, Amandine Varenne³, Frédéric Hérault¹,<fnref n="2" />, Hervé Chapuis⁴, Christophe Alleno⁵, Patrice Dehais⁶, Alain Vignal⁶, Thierry Burlot³ and Pascale Le Roy¹,<fnref n="2" />*

## Abstract

**Background:** The genetic architecture of egg production and egg quality traits, i.e. the quantitative trait loci (QTL) that influence these traits, is still poorly known. To date, 33 studies have focused on the detection of QTL for laying traits in chickens, but less than 10 genes have been identified. The availability of a high-density SNP (single nucleotide polymorphism) chicken array developed by Affymetrix, i.e. the 600K Affymetrix® Axiom® HD genotyping array offers the possibility to narrow down the localization of previously detected QTL and to detect new QTL. This high-density array is also anticipated to take research beyond the classical hypothesis of additivity of QTL effects or of QTL and environmental effects. The aim of our study was to search for QTL that influence laying traits using the 600K SNP chip and to investigate whether the effects of these QTL differed between diets and age at egg collection.

**Results:** One hundred and thirty-one QTL were detected for 16 laying traits and were spread across all marked chromosomes, except chromosomes 16 and 25. The percentage of variance explained by a QTL varied from 2 to 10 % for the various traits, depending on diet and age at egg collection. Chromosomes 3, 9, 10 and Z were overrepresented, with more than eight QTL on each one. Among the 131 QTL, 60 had a significantly different effect, depending on diet or age at egg collection. For egg production traits, when the QTL × environment interaction was significant, numerous inversions of sign of the SNP effects were observed, whereas for egg quality traits, the QTL × environment interaction was mostly due to a difference of magnitude of the SNP effects.

**Conclusions:** Our results show that numerous QTL influence egg production and egg quality traits and that the genomic regions, which are involved in shaping the ability of layer chickens to adapt to their environment for egg production, vary depending on the environmental conditions. The next question will be to address what the impact of these genotype × environment interactions is on selection.

## Background

Over the last decades, layer chicken lines have been selected and improved for egg production and egg quality performance. However, the genetic architectures of the underlying traits, i.e. the quantitative trait loci (QTL) that influence these traits, are still poorly known. To date, 33 studies have focused on the detection of QTL for laying traits in chickens, but less than 10 genes have been identified [1].

The high-density array for chicken recently developed by Affymetrix, i.e. the 600K Affymetrix® Axiom® HD genotyping array [2], offers the possibility to use high-density genotype data for genomic selection in laying hens. It will also contribute to improve the localization of previously detected QTL and to detect new QTL. This high-density array is also anticipated to take research beyond the classical hypothesis of additivity of QTL effects or of QTL and environmental effects. Indeed, some studies suggest that genotype × environment (G × E) interactions may explain a large part of the phenotypic variance in laying traits in chickens [3]. However, to date, no study has tested the robustness of QTL across environments.

<footnote>*Correspondence: pascale.leroy@rennes.inra.fr</footnote>

<footnote n="2">² Agrocampus Ouest, UMR1348 PEGASE, 65 Rue de Saint Brieuc, 35042 Rennes, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
Therefore, in this study, a genome-wide association study (GWAS) using the 600K Affymetrix® Axiom® HD genotyping array was conducted to detect QTL that influence egg production and egg quality traits in layer chickens. In order to investigate whether the QTL detected differed between environmental conditions, animals were divided into two groups that were fed a different diet.

# Methods

## Animals

The population studied consisted of 438 sires from a commercial pure line that was created and selected by NOVOGEN (Le Foeil, France) and 31,381 of their F1 crossbred female offspring. Hens were hatched in three batches in November 2010, May 2011 and November 2011. At 18 weeks of age, they were housed in a production farm in collective cages that contained 12 half-sisters of the same sire. The hens laid from 18 to 75 weeks of age. Fifty percent of the hens were fed ad libitum a high-energy diet (HE) that supplied 2881 kcal of metabolizable energy (ME) (1342 cages) and 50% were fed ad libitum a low-energy diet (LE) that supplied 2455 kcal of ME (1346 cages).

## Genotyping

Blood was sampled from the brachial veins of the sires, DNA was extracted and hybridized on the 600K Affymetrix® Axiom® HD genotyping array [2] by Ark-Genomics (Edinburgh, UK). In total, 438 sires were genotyped for 580,961 SNPs that were distributed over chromosomes 1 to 28, two linkage groups (LGE22C19W28_E50C23 and LGE64) and the two sex chromosomes, along with a group of 7883 markers of unknown location. Genotypes were filtered in five successive steps: (1) 14 SNPs on chromosome W with a call rate less than 5% were excluded; (2) none of the animals had a call rate less than 95%; (3) 260,945 SNPs with a minor allele frequency less than 0.05 were excluded; (4) 9041 SNPs with a call rate less than 95% were excluded; and (5) 26,318 SNPs that deviated significantly (P < 5%) from Hardy-Weinberg equilibrium were excluded. Finally, 284,643 SNPs remained for analysis and no individual was excluded.

## Measurement of traits

In this paper, traits are named according to Animal Trait Ontology for Livestock [4].

Egg production was recorded daily from week 18 to week 75 and egg production rate (EPR) of one cage (in %) was calculated by dividing the number of eggs produced by the number of hen days in the cage. Based on the laying curve, four periods of production were defined: increasing from week 18 to week 30 (EPR1); plateau from week 31 to week 49 (EPR2); persistence from week 50 to

week 75 (EPR3); and global production from week 18 to week 75 (EPR).

At about 50 weeks of age, and then again at 70 weeks, all eggs produced on the farm were collected and egg quality traits were measured by the Zootests company (Ploufragan, France) on 27,747 week-50 eggs and 25,964 week-70 eggs (Table 1). The first step consisted of measuring the short length of the egg (SLE, in mm) and egg weight (EW, in g), before calculating egg shell shape (ESshape) as: $ESshape = (SLE/10)/(EW/10)^{1/3}$. Second, shell color was measured with a Minolta chromameter and three traits were recorded: redness of egg shell a* (RSS), yellowness of egg shell b* (YSS) and lightness of egg shell L* (LSS). Egg shell color was then calculated as: $ESC = 100 - (L^* - A^* - B^*)$. Third, shell strength was measured using a compression machine to evaluate the static stiffness of the shell. The egg was compressed between two flat plates moving at constant speed and at a constant force of 15 N to record egg shell stiffness (ESS-TIF, in mm). Egg shell strength is the maximum force recorded before fracture of the shell (ESS in N). Then, each egg was broken and albumen height (H) was measured using a tripod. The Haugh unit (HU) measure of albumen firmness was then calculated as: $HU = 100 \log (H - 1.7 EW^{0.37} + 7.57)$ [5]. The yolk was weighed and a yolk index (YOLKIND) was calculated as yolk weight divided by EW. Finally, the egg was scored for blood and meat spots (EMTSP) on a scale of 0 (without spots)–3 (many spots).

## Statistical analysis

First, egg measurements were adjusted for environmental effects. Separately for each hatch, the covariates and fixed effects were tested using the SAS® 9.2 GLM procedure based on the following linear models: $Y_{ijlmn} = s_i + d_j + b_l + c_m + f_n + E_{ijlmn}$, for EPR, EPR1, EPR2 and EPR3 and $Y_{ijklmno} = s_i + d_j + a_k + b_l + c_m + f_n + e_o + \beta_1 W_{ijklmno} + \beta_2 R_{ijklmno} + E_{ijklmno}$, for SLE, EW, ESshape, RSS, YSS, LSS, ESC, ESSTIF, ESS, HU, YOLKIND and EMTSP, where $Y_{ijlmn}$ and $Y_{ijklmno}$ are trait values, $s_i$ is the fixed effect of sire i (438 levels); $d_j$ is the fixed effect of diet j (two levels: HE or LE); $a_k$ is the fixed effect of age class k (two levels: 50 or 70 weeks); $b_l$,$c_m$ and $f_n$ represent the location of the cage in the building, respectively the fixed effect of battery l (four levels), the fixed effect of column m (two levels: middle or edges of the battery) and the fixed effect of floor n (two levels); $e_o$ is the fixed effect of the person who made the measurement (eight to 10 levels according to quality trait), $W_{ijklmno}$ and $R_{ijklmno}$ are the waiting time between sample and egg measurement (in days) covariate and the age of the hen (in days) covariate; $E_{ijlmn}$ and $E_{ijklmno}$ are random residual variables.
<table><caption>Table 1 Summary statistics on phenotypic data for the high energy (HE) and low energy (LE) diets</caption><thead><tr><th rowspan="2">Trait<br>Diet</th><th colspan="2">Number of records<sup>a</sup></th><th colspan="2">Mean</th><th colspan="2">Standard deviation</th></tr><tr><th>HE</th><th>LE</th><th>HE</th><th>LE</th><th>HE</th><th>LE</th></tr></thead><tbody><tr><td>Egg production traits</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>EPR (%)</td><td>1342</td><td>1340</td><td>84.22</td><td>83.26</td><td>6.12</td><td>7.17</td></tr><tr><td>EPR1 (%)</td><td>1344</td><td>1342</td><td>72.12</td><td>76.44</td><td>9.40</td><td>9.35</td></tr><tr><td>EPR2 (%)</td><td>1342</td><td>1338</td><td>91.14</td><td>89.28</td><td>6.06</td><td>8.16</td></tr><tr><td>EPR3 (%)</td><td>1344</td><td>1342</td><td>85.32</td><td>81.97</td><td>6.35</td><td>7.62</td></tr><tr><td>Egg quality traits at 50 weeks of age</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ESshape</td><td>12,375</td><td>12,403</td><td>1.1</td><td>1.1</td><td>0.02</td><td>0.02</td></tr><tr><td>EW (g)</td><td>13,547</td><td>13,603</td><td>61.11</td><td>60.59</td><td>4.78</td><td>4.74</td></tr><tr><td>SLE (mm)</td><td>12,378</td><td>12,409</td><td>43.37</td><td>43.32</td><td>1.29</td><td>1.28</td></tr><tr><td>ESC</td><td>13,528</td><td>13,580</td><td>24.78</td><td>25.31</td><td>9.66</td><td>10.09</td></tr><tr><td>LSS</td><td>13,548</td><td>13,596</td><td>66.77</td><td>66.91</td><td>4.57</td><td>4.72</td></tr><tr><td>RSS</td><td>13,546</td><td>13,605</td><td>13.72</td><td>13.42</td><td>3.16</td><td>3.27</td></tr><tr><td>YSS</td><td>13,445</td><td>13,493</td><td>28.3</td><td>28.22</td><td>2.65</td><td>2.81</td></tr><tr><td>ESS (N)</td><td>12,379</td><td>12,410</td><td>39.15</td><td>39.58</td><td>7.51</td><td>7.33</td></tr><tr><td>ESSTIF (mm)</td><td>12,388</td><td>12,418</td><td>193.42</td><td>191.49</td><td>29.53</td><td>28.57</td></tr><tr><td>EMTSP</td><td>6959</td><td>6283</td><td>1.68</td><td>1.55</td><td>0.74</td><td>0.72</td></tr><tr><td>HU</td><td>6938</td><td>6257</td><td>73.03</td><td>74.45</td><td>7.69</td><td>8.02</td></tr><tr><td>YOLKIND</td><td>6157</td><td>6808</td><td>0.27</td><td>0.27</td><td>0.02</td><td>0.02</td></tr><tr><td>Egg quality traits at 70 weeks of age</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ESshape</td><td>11,315</td><td>10,983</td><td>1.1</td><td>1.1</td><td>0.02</td><td>0.02</td></tr><tr><td>EW (g)</td><td>12,945</td><td>12,804</td><td>61.05</td><td>60.61</td><td>5.04</td><td>4.95</td></tr><tr><td>SLE (mm)</td><td>11,329</td><td>11,004</td><td>43.41</td><td>43.36</td><td>1.37</td><td>1.35</td></tr><tr><td>ESC</td><td>12,929</td><td>12,259</td><td>25.4</td><td>25.48</td><td>10.18</td><td>10.64</td></tr><tr><td>LSS</td><td>12,937</td><td>12,272</td><td>67.09</td><td>66.99</td><td>4.92</td><td>5.05</td></tr><tr><td>RSS</td><td>12,942</td><td>12,274</td><td>13.4</td><td>13.3</td><td>3.27</td><td>3.42</td></tr><tr><td>YSS</td><td>12,833</td><td>12,188</td><td>28.35</td><td>28.23</td><td>2.65</td><td>2.86</td></tr><tr><td>ESS (N)</td><td>11,327</td><td>11,007</td><td>37.4</td><td>37.82</td><td>7.53</td><td>7.75</td></tr><tr><td>ESSTIF (mm)</td><td>11,332</td><td>11,004</td><td>192.83</td><td>190.46</td><td>29.98</td><td>30.53</td></tr><tr><td>EMTSP</td><td>6459</td><td>6705</td><td>1.63</td><td>1.63</td><td>0.72</td><td>0.71</td></tr><tr><td>HU</td><td>6437</td><td>6683</td><td>64.88</td><td>67.28</td><td>9.31</td><td>9.46</td></tr><tr><td>YOLKIND</td><td>5729</td><td>5491</td><td>0.27</td><td>0.27</td><td>0.02</td><td>0.02</td></tr></tbody></table>

For each trait, a sub-model that took only effects exceeding the significance level (P < 0.2) into account was retained. Raw data were then adjusted using the estimates of all effects in this model, except the sire effect. Distributions of the adjusted data were tested for each trait and extreme individual values, i.e. values that were more than four phenotypic standard deviations from the mean, were discarded. Finally, for each trait, the “performance” of one sire was calculated as the mean of its daughters’ adjusted performances.

To examine the genetic architecture of a trait, all data were considered together, i.e. one mean per sire across the two diets and ages, whereas to examine the QTL effect for each condition, data were considered separately

for each condition, i.e. two means per sire (HE and LE diets) for egg production traits and four means per sire (for egg collection at 50 and 70 weeks of age, and HE and LE diets) for egg quality traits.

## GWAS analysis

According to the recommendations made by Teysse-dre et al. [6], data were analyzed using a mixed model that takes pedigree kinship into account. For each SNP (1 to 284,643) and each trait (1 to 16), the following mixed model was applied: $Z = 1\mu + X\alpha + E$, with $V(E) = A\sigma_g^2 + I\sigma_e^2$, where Z is the vector of the sire performances; μ is the general mean; X is the incidence matrix of genotypes for the SNP evaluated; α is the allele substitution

<footnote><sup>a</sup> Number of cages for production traits and number of eggs for quality traits</footnote>
effect of the SNP for the trait [7]; E is the vector of residuals with variance-covariance matrix V(E); A is the pedigree kinship matrix; I is the identity matrix; σ²g is the genetic variance and σ²e is the environmental variance.

This model was fitted using the BLUPF90 program [8]. Chromosome-wide thresholds (P < 1 %) to test the H0 hypothesis, i.e. no effect of the SNP on the trait, were estimated using the method of Müller et al. [9] and Müller's software [10], which offsets the heavy computations to parallel implementation. Genome-wide thresholds (P < 5 %) were calculated according to the Bonferroni correction as chromosome-wide thresholds at P < 0.0015 (0.05/32 chromosomes). For the group of SNPs with unknown locations, the threshold was set at P < 1e-5 (0.05/3681 SNPs), according to the Bonferroni correction.

### QTL detection

For each chromosome, a QTL was detected when one SNP was genome-wide significant (P < 5 %). The confidence interval of this QTL was defined with the adjacent SNPs that were chromosome-wide significant (P < 1 %). For each QTL, the SNP with the highest estimated effect (α) was called the “top SNP”. Based on the estimate of the allele substitution effect for the top SNP, α, the genetic variance explained by the top SNP was calculated as 2p(1 – p)α², where p is the minor allele frequency. The percentage of variance explained by the top SNP was then calculated for each trait by dividing by the variance of the sire performances, and called “% of explained variance”.

### QTL × environment interaction

For each trait, the GWAS analysis was repeated separately, as above, for the four datasets for each condition (egg collections at 50 and 70 weeks of age, and HE and LE diets) in order to estimate the effect of the top SNP for each QTL in the two environments, α₁ and α₂, i.e. egg collection at 50 vs. 70 weeks of age and HE vs. LE diet. QTL × age and QTL × diet interactions were tested by comparing the two corresponding estimates using a Z test statistic. Variances explained by the top SNPs for each condition were calculated, and the residual variances for the different conditions were compared with a Fisher test to verify their equality (which was always accepted). The Z test statistic was calculated as follows:

$$Z = \frac{|\alpha_1 - \alpha_2|}{\sigma} \times \sqrt{n},$$

where σ² is the average residual variance and n is the average number of sires over the two conditions. Z followed a normal distribution because n was large. However, if the interaction effect is large, some QTL may not be detected by analyzing the whole dataset. Indeed,

additional QTL, i.e. QTL that were not detected when analyzing the whole dataset, were identified by analyzing the four within-condition datasets separately. In such cases, we also tested the QTL × environment interaction. The QTL were defined using the results of the within-condition GWAS analyses, and the estimated effects of the top SNP for each QTL were compared using the Z test statistic described previously.

### Results and discussion

Number of observations, means and standard deviations of traits are in Table 1 for the raw data. Differences in average performance between the HE and LE diets ranged from 0 to 0.5 standard deviations. Differences in performance between the two age groups (egg collections at 50 and 70 weeks of age) ranged from 0 to 0.9 standard deviations. Overall, these differences remained small, even if, in general, they were significant.

Correlations between the sires' “performances” depending on diet and age at collection are in Table 2. They are much lower for production traits than for egg quality traits, except for EMTSP. For quality traits, they are lower for alternate diets than for alternate ages and they vary more between traits. At a broad level, these correlations confirm the putative presence of interaction effects.

<table><caption>Table 2 Phenotypic correlations between the “performances” of sires for alternate diets and alternate ages</caption><thead><tr><th></th><th>Diet<sup>a</sup></th><th>Age<sup>b</sup></th></tr></thead><tbody><tr><td>Egg production traits</td><td></td><td></td></tr><tr><td>EPR</td><td>0.42</td><td>—</td></tr><tr><td>EPR1</td><td>0.25</td><td>—</td></tr><tr><td>EPR2</td><td>0.28</td><td>—</td></tr><tr><td>EPR3</td><td>0.10</td><td>—</td></tr><tr><td>Egg quality traits</td><td></td><td></td></tr><tr><td>ESshape</td><td>0.67</td><td>0.78</td></tr><tr><td>EW</td><td>0.78</td><td>0.92</td></tr><tr><td>SLE</td><td>0.76</td><td>0.87</td></tr><tr><td>ESC</td><td>0.73</td><td>0.83</td></tr><tr><td>LSS</td><td>0.74</td><td>0.84</td></tr><tr><td>RSS</td><td>0.72</td><td>0.82</td></tr><tr><td>YSS</td><td>0.56</td><td>0.62</td></tr><tr><td>ESS</td><td>0.67</td><td>0.77</td></tr><tr><td>ESSTIF</td><td>0.71</td><td>0.77</td></tr><tr><td>EMTSP</td><td>0.42</td><td>0.47</td></tr><tr><td>HU</td><td>0.50</td><td>0.57</td></tr><tr><td>YOLKIND</td><td>0.56</td><td>0.67</td></tr></tbody></table>

<footnote>a Correlations between the sires' performances depending on the diet</footnote>

<footnote>b Correlations between the sires' performances depending on the age of measurement</footnote>
GWAS analyses

For egg production traits, the GWAS analysis on the whole dataset resulted in 1202 significant tests at the 1 % chromosome significance level and identified 861 SNPs with an effect on at least one egg production trait. For egg quality traits, it resulted in 8116 significant tests at the 1 % chromosome significance level and identified 5384 SNPs with an effect on at least one trait.

Analyses of the four subsets estimated the effect of each SNP for each of the four conditions investigated (HE diet, LE diet, egg collections at 50 and 70 weeks of age). Correlations between estimates of SNP effects based on HE and LE diets ranged from 0.16 to 0.80 (Table 3). They were lower for egg production traits than for egg quality traits. Figure 1 shows an example of the correlations between estimates of SNP effect for the HE and LE diets for egg shell color and EPR1. Estimates of the effects of each SNP for each diet were plotted and SNPs within a QTL that had an interaction with diet were investigated, in order to determine the direction of the interaction. For egg quality traits, in most cases, the sign of the allele substitution effect was not reversed between diets. Instead, interactions were due to deviations of the magnitude of the SNP effects, except for two traits with a low correlation between diets, YOLKIND and EMTSP. However, for

<table><caption>Table 3 Correlations between the effects of SNPs for alternate diets and alternate ages</caption><thead><tr><th></th><th>Diet<sup>a</sup></th><th>Age<sup>b</sup></th></tr></thead><tbody><tr><td>Egg production traits</td><td></td><td></td></tr><tr><td>EPR</td><td>0.35</td><td>-</td></tr><tr><td>EPR1</td><td>0.16</td><td>-</td></tr><tr><td>EPR2</td><td>0.27</td><td>-</td></tr><tr><td>EPR3</td><td>0.27</td><td>-</td></tr><tr><td>Egg quality traits</td><td></td><td></td></tr><tr><td>ESshape</td><td>0.69</td><td>0.8</td></tr><tr><td>EW</td><td>0.8</td><td>0.93</td></tr><tr><td>SLE</td><td>0.78</td><td>0.88</td></tr><tr><td>ESC</td><td>0.75</td><td>0.86</td></tr><tr><td>LSS</td><td>0.8</td><td>0.87</td></tr><tr><td>RSS</td><td>0.76</td><td>0.86</td></tr><tr><td>YSS</td><td>0.61</td><td>0.74</td></tr><tr><td>ESS</td><td>0.66</td><td>0.78</td></tr><tr><td>ESSTIF</td><td>0.7</td><td>0.76</td></tr><tr><td>EMTSP</td><td>0.46</td><td>0.54</td></tr><tr><td>HU</td><td>0.54</td><td>0.62</td></tr><tr><td>YOLKIND</td><td>0.6</td><td>0.72</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Examples of correlations between estimates of SNP effects for egg shell color (a) and egg production rate 1 (b) for the high-energy diet (x-axis) and the low-energy diet (y-axis). Chromosome-wide significant SNPs (P < 1 %) located in QTL that show an interaction with diet are plotted in black; correlations between estimates of SNPs effects are equal to 0.75 for egg shell color and to 0.16 for egg production rate</figcaption></figure>

egg production traits, in most cases, the sign of the allele substitution effect was reversed. Correlations between estimates of the SNP effects at 50 and 70 weeks of age ranged from 0.54 to 0.93.

<footnote><sup>a</sup> Correlations between the effects of SNPs depending on the diet</footnote>

<footnote><sup>b</sup> Correlations between the effects of SNPs depending on the age of measurement</footnote>
## QTL detected

To analyze the genetic architecture of each trait, all the data were analyzed together, whereas to determine which QTL were specific to performances for egg collection at 50 and 70 weeks of age or with HE and LE diets, their detection was carried out by using the data divided into four subsets. To calculate the number of QTL detected and to give them a name, we considered the QTL according to their chromosomal localization, regardless of which traits they influenced. For QTL that were identified by analyzing the whole dataset, 11 had a significant effect on at least one egg production trait and 58 QTL had an effect on at least one egg quality trait (see Additional file 1). When GWAS was applied separately to the four subsets, 19 additional QTL with a significant effect on at least one egg production trait and 46 additional QTL with an effect on at least one egg quality trait were detected (see Additional files 2, 3). Three QTL were common to egg production traits and egg quality traits.

QTL localizations are in Additional files 1, 2 and 3, along with estimates of the variance explained by each QTL for each trait for each condition (HE diet, LE diet, egg collections at 50 and 70 weeks of age) are in Table 4.

### Egg production (EPR, EPR1, EPR2, EPR3) traits

For egg production traits, analysis of the whole dataset identified 11 QTL that had an effect on at least one

trait and separate analysis of the four subsets detected 21 additional QTL. Previously, several QTL for egg production have been reported [11, 12] that are located close to QTL65 and QTL67 identified here.

For EPR1, 23 QTL were detected on 15 chromosomes, which indicates that it has a polygenic determination. Together, these QTL explained 34 % of the variance of EPR1, regardless of environmental condition. For EPR and EPR2, nine and six QTL were detected, which explained 26.8 and 24.1 % of the variance, respectively. Although fewer QTL were detected for EPR and EPR2, they had a larger effect when the data are considered together. QTL7 had an effect on the four traits (EPR, EPR1, EPR2 and EPR3) but its effect was diet-dependent. Generally, QTL influencing EPR3 also had an effect on EPR1.

### Egg weight (EW) and egg form (SLE and ESshape) traits

Analysis of the whole dataset identified 13 QTL that had an effect on at least one of the following traits EW, SLE and ESshape and analysis of the four subsets separately detected 14 additional QTL. Together, these QTL explained 65.4 % of the variance of ESshape, 55 % of the variance of EW but only 28.7 % of the variance of SLE.

QTL106 was common to all three egg weight and shape traits and was located on chromosome 26, between 1.2 and 1.3 Mb. Two QTL for EW have been reported on

<table><caption>Table 4 QTL detected for each trait: number of QTL detected, number of QTL that show an interaction, percentage of the variance explained across conditions and percentage of the variance explained within each condition</caption><thead><tr><th>Trait</th><th>Nb QTL</th><th>QTL diet</th><th>QTL age</th><th>Var (%)</th><th>Var LE (%)</th><th>Var HE (%)</th><th>Var 50 (%)</th><th>Var 70 (%)</th></tr></thead><tbody><tr><td colspan="9"><strong>Egg production traits</strong></td></tr><tr><td>EPR</td><td>9</td><td>7</td><td>-</td><td>26.8</td><td>34.07</td><td>14.61</td><td>NA</td><td>NA</td></tr><tr><td>EPR1</td><td>23</td><td>19</td><td>-</td><td>33.99</td><td>85.31</td><td>44.61</td><td>NA</td><td>NA</td></tr><tr><td>EPR2</td><td>6</td><td>3</td><td>-</td><td>24.07</td><td>22.03</td><td>10.57</td><td>NA</td><td>NA</td></tr><tr><td>EPR3</td><td>7</td><td>7</td><td>-</td><td>16.25</td><td>34.04</td><td>4.97</td><td>NA</td><td>NA</td></tr><tr><td colspan="9"><strong>Egg quality traits</strong></td></tr><tr><td>ESshape</td><td>15</td><td>3</td><td>1</td><td>65.36</td><td>45.19</td><td>64.42</td><td>59.55</td><td>59.21</td></tr><tr><td>EW</td><td>12</td><td>1</td><td>0</td><td>55</td><td>49.5</td><td>50.58</td><td>58.48</td><td>47.96</td></tr><tr><td>SLE</td><td>7</td><td>0</td><td>0</td><td>28.7</td><td>26.37</td><td>24.23</td><td>29.05</td><td>25.2</td></tr><tr><td>ESC</td><td>27</td><td>3</td><td>1</td><td>105.35</td><td>96</td><td>91.37</td><td>93.76</td><td>101.89</td></tr><tr><td>LSS</td><td>24</td><td>0</td><td>0</td><td>102.11</td><td>90.47</td><td>90.02</td><td>94.44</td><td>95.84</td></tr><tr><td>RSS</td><td>19</td><td>3</td><td>1</td><td>85.97</td><td>81.79</td><td>68.54</td><td>76.94</td><td>81.68</td></tr><tr><td>YSS</td><td>15</td><td>6</td><td>0</td><td>61.1</td><td>58.58</td><td>43.54</td><td>51.71</td><td>51.02</td></tr><tr><td>ESS</td><td>8</td><td>3</td><td>3</td><td>27.16</td><td>22.7</td><td>26.37</td><td>21.88</td><td>27.47</td></tr><tr><td>ESSTIF</td><td>8</td><td>2</td><td>1</td><td>32.88</td><td>33.82</td><td>23.73</td><td>29.36</td><td>30.16</td></tr><tr><td>EMTSP</td><td>9</td><td>3</td><td>2</td><td>42.55</td><td>24.87</td><td>31.16</td><td>33.93</td><td>27.09</td></tr><tr><td>HU</td><td>11</td><td>2</td><td>7</td><td>39.26</td><td>24.8</td><td>34.79</td><td>36.08</td><td>28.62</td></tr><tr><td>YOLKIND</td><td>16</td><td>2</td><td>2</td><td>62.45</td><td>72.49</td><td>40.71</td><td>56.95</td><td>47.64</td></tr></tbody></table>

Nb QTL number of QTL detected for each trait, QTL diet number of QTL that show an interaction with diet; QTL age number of QTL that show an interaction with age, var LE variance explained for the low energy diet (%), var HE variance explained for the high energy diet, in percent; var 50 variance explained at 50 weeks of age (%), var 70 variance explained at 70 weeks of age (%)

## QTL detected
chromosome 26 at 2 and 3 Mb [13]. QTL 106 explained about 4 % of the variance of EW and SLE, but only 1.6 % of the variance of ESshape.

Three other QTL common to EW and SLE were detected (see Additional file 1), of which two, QTL121 and QTL122, were located on chromosome Z. This agrees with studies that have shown that chromosome Z plays a role in the regulation of EW [14, 15].

### Egg shell color (ESC, RSS, LSS and YSS) traits

For egg shell color traits, analysis of the whole dataset identified 22 QTL that had an effect on at least one trait and analysis of the four subsets separately detected 17 additional QTL. Several QTL have been reported for egg shell color [16, 17] that are localized close to QTL73 and QTL10. Together, these QTL explained more than 100 % of the variance of ESC and LSS, which, although their effects may be overestimated, suggests the existence of other interactions such as dominance or epistasis, which were not considered in the model.

The number of QTL detected and their distribution along the genome suggest that egg shell color has a polygenic determinism. Most of the QTL for ESC also had an effect on RSS and YSS or on YSS, but only a few had effects on RSS, LSS and YSS. Indeed, only two QTL (QTL20 and QTL119) had an effect on all four traits.

### Egg shell solidity (ESSTIF and ESS) traits

For egg shell solidity traits, analysis of the whole dataset identified seven QTL that had an effect on at least one trait, i.e. five for ESSTIF and two for ESS but no QTL was common to both traits. Analysis of the four subsets separately detected seven additional QTL that had a significant effect on at least one of these traits, i.e. two for ESSTIF and five for ESS. Overall, the same number of QTL was found for ESS and ESSTIF, but the total variance explained by these QTL was larger for ESSTIF than for ESS. QTL103 was common to both traits and had a pleiotropic effect.

### Internal egg quality (YOLKIND, HU and EMTSP) traits

For internal egg quality traits, analysis of the whole dataset identified 18 QTL and analysis of the four subsets separately detected 15 additional QTL that had a significant effect on YOLKIND and/or HU. One of these QTL (QTL123) was previously reported by [18]. YOLKIND, for which 16 QTL were identified, is the trait that had the largest proportion of variance explained, followed by EMTSP for which nine QTL were detected.

Among the 33 QTL detected for the three internal egg quality traits, QTL33 was common to HU and YOLKIND and QTL100 was common to HU and EMTSP. For YOLKIND, QTL88 and QTL89 were localized on

chromosome 20 between 11 and 12 Mb, nearby a QTL detected for EMTSP. Previously, a QTL for HU, on chromosome 20 at 11 Mb was reported by [13]. Together, these results suggest that this region has a pleiotropic effect.

### QTL x environment interactions

Our main aim was to determine whether diet and age had an effect on the genetic architecture of egg production and egg quality traits. As described in the Methods section, we estimated the interaction of all QTL detected with diet or age at egg collection. Forty-three QTL showed a significant interaction with diet, 12 QTL showed a significant interaction with age at egg collection and five QTL showed a significant interaction with both diet and age at egg collection, according to the Z test. A same QTL, i.e. the same location on the genome, could happen to be in interaction or not, depending on the trait. This is mostly due to the fact that the top SNP of a QTL can differ between traits. For example, QTL84 was detected for both ESshape and ESC but two different top SNPs were identified, i.e. AX-75934213, which showed an interaction with diet, and AX-80996918, which showed an interaction with age at egg collection.

### Egg production (EPR, EPR1, EPR2, EPR3) traits

Among the 32 QTL detected for egg production traits, 26 had a significant interaction with diet (Table 5). Among the QTL that showed an interaction with diet, 21 explained a larger proportion of variance for the LE diet and five explained a larger proportion of variance for the HE diet. For EPR1, 22 QTL were detected of which 19 showed an interaction with diet. This suggests that EPR1 has a polygenic determinism and that it is highly subject to interaction with diet. For EPR2 and ERP3, only three and seven of the detected QTL, respectively, showed an interaction with diet. For these two traits, the QTL that showed an interaction with diet explained a larger proportion of variance for the LE diet than for HE diet. Finally, seven QTL were identified for EPR, of which QTL68 explained a larger proportion of variance for the HE diet.

### Egg weight (EW) and egg form (SLE and ESshape) traits

Among the 27 QTL previously detected for egg weight and egg form, four showed an interaction with diet and one showed a significant interaction with age. All the QTL that showed an interaction with diet accounted for a larger proportion of variance explained for the HE diet than for the LE diet. QTL59 was detected for EW on chromosome 10 and had a greater effect for the HE diet than for the LE diet (5.0 % of variance explained vs. 2.4 %). The three other QTL which had an interaction
<table><caption>**Table 5 QTL with significantly different effects between diets**</caption><thead><tr><th>QTL</th><th>Trait</th><th>Chr</th><th>Var LE (%)</th><th>Var HE (%)</th><th>Z Diet</th></tr></thead><tbody><tr><td>SNP5</td><td>EPR</td><td>0</td><td>5.06</td><td>1.16</td><td>4.41***</td></tr><tr><td>QTL7</td><td>EPR</td><td>2</td><td>5.44</td><td>2.61</td><td>2.75*</td></tr><tr><td>QTL23</td><td>EPR</td><td>4</td><td>4.8</td><td>0.54</td><td>9.27***</td></tr><tr><td>QTL47</td><td>EPR</td><td>9</td><td>4.92</td><td>0.53</td><td>7.14***</td></tr><tr><td>QTL68</td><td>EPR</td><td>11</td><td>0.02</td><td>4.41</td><td>7.78***</td></tr><tr><td>QTL79</td><td>EPR</td><td>17</td><td>2.39</td><td>0.83</td><td>14.23***</td></tr><tr><td>QTL90</td><td>EPR</td><td>20</td><td>4.89</td><td>0.01</td><td>8.48***</td></tr><tr><td>QTL7</td><td>EPR1</td><td>2</td><td>3.86</td><td>1.79</td><td>22.38***</td></tr><tr><td>QTL14</td><td>EPR1</td><td>3</td><td>5.98</td><td>2.1</td><td>20.68***</td></tr><tr><td>QTL28</td><td>EPR1</td><td>5</td><td>3.48</td><td>1.01</td><td>11.56***</td></tr><tr><td>QTL38</td><td>EPR1</td><td>7</td><td>2.64</td><td>1.12</td><td>18.8***</td></tr><tr><td>QTL39</td><td>EPR1</td><td>7</td><td>3.19</td><td>2.27</td><td>13.72***</td></tr><tr><td>QTL40</td><td>EPR1</td><td>7</td><td>3.01</td><td>1.37</td><td>13.76***</td></tr><tr><td>QTL42</td><td>EPR1</td><td>8</td><td>4.15</td><td>6.22</td><td>14.25***</td></tr><tr><td>QTL44</td><td>EPR1</td><td>9</td><td>1.42</td><td>3.57</td><td>3.74***</td></tr><tr><td>QTL46</td><td>EPR1</td><td>9</td><td>0.01</td><td>4.18</td><td>7.8***</td></tr><tr><td>QTL47</td><td>EPR1</td><td>9</td><td>5.52</td><td>0.46</td><td>7.28***</td></tr><tr><td>QTL51</td><td>EPR1</td><td>9</td><td>4.26</td><td>2.66</td><td>22.76***</td></tr><tr><td>QTL71</td><td>EPR1</td><td>12</td><td>0.02</td><td>4.26</td><td>7.04***</td></tr><tr><td>QTL75</td><td>EPR1</td><td>13</td><td>3.27</td><td>2.21</td><td>19.83***</td></tr><tr><td>QTL79</td><td>EPR1</td><td>17</td><td>5.85</td><td>1.68</td><td>21.41***</td></tr><tr><td>QTL90</td><td>EPR1</td><td>20</td><td>7.09</td><td>0.21</td><td>11.12***</td></tr><tr><td>QTL91</td><td>EPR1</td><td>21</td><td>4.2</td><td>1.05</td><td>12.59***</td></tr><tr><td>QTL99</td><td>EPR1</td><td>23</td><td>3.92</td><td>0.24</td><td>8.92***</td></tr><tr><td>QTL110</td><td>EPR1</td><td>27</td><td>3.48</td><td>1.19</td><td>14.93***</td></tr><tr><td>QTL116</td><td>EPR1</td><td>28</td><td>9.68</td><td>0</td><td>13.09***</td></tr><tr><td>QTL7</td><td>ESshape</td><td>2</td><td>4.9</td><td>1.84</td><td>3.66***</td></tr><tr><td>QTL67</td><td>ESshape</td><td>11</td><td>5.07</td><td>0.89</td><td>4.55***</td></tr><tr><td>QTL87</td><td>ESshape</td><td>20</td><td>4.46</td><td>0.68</td><td>6.46***</td></tr><tr><td>SNP4</td><td>EPR3</td><td>0</td><td>6.04</td><td>0.95</td><td>5.77***</td></tr><tr><td>SNP5</td><td>EPR3</td><td>0</td><td>6.87</td><td>0.9</td><td>6.39***</td></tr><tr><td>QTL7</td><td>EPR3</td><td>2</td><td>5.65</td><td>0.94</td><td>5.42***</td></tr><tr><td>QTL26</td><td>EPR3</td><td>4</td><td>3.19</td><td>0.07</td><td>15.03***</td></tr><tr><td>QTL51</td><td>EPR3</td><td>9</td><td>4.57</td><td>0.78</td><td>18.88***</td></tr><tr><td>QTL75</td><td>EPR3</td><td>13</td><td>3.6</td><td>1.3</td><td>18.47***</td></tr><tr><td>QTL90</td><td>EPR3</td><td>20</td><td>4.12</td><td>0.03</td><td>8.24***</td></tr><tr><td>QTL1</td><td>ESshape</td><td>1</td><td>1.88</td><td>5.4</td><td>3.58***</td></tr><tr><td>QTL45</td><td>ESshape</td><td>9</td><td>1.58</td><td>4.7</td><td>3.44***</td></tr><tr><td>QTL84</td><td>ESshape</td><td>19</td><td>1.29</td><td>4.11</td><td>5.9***</td></tr><tr><td>QTL59</td><td>EW</td><td>10</td><td>2.37</td><td>4.98</td><td>2.74*</td></tr><tr><td>QTL15</td><td>ESC</td><td>3</td><td>4.97</td><td>1.7</td><td>3.16**</td></tr><tr><td>QTL74</td><td>ESC</td><td>13</td><td>2.13</td><td>4.55</td><td>2.96**</td></tr><tr><td>QTL83</td><td>ESC</td><td>19</td><td>4.05</td><td>1.41</td><td>3.13**</td></tr><tr><td>QTL69</td><td>RSS</td><td>12</td><td>5.62</td><td>2.46</td><td>3.82***</td></tr><tr><td>QTL83</td><td>RSS</td><td>19</td><td>3.95</td><td>1.09</td><td>3.37**</td></tr><tr><td>QTL117</td><td>RSS</td><td>LGE22</td><td>0.47</td><td>2.41</td><td>4.97***</td></tr><tr><td>QTL13</td><td>YSS</td><td>3</td><td>4.88</td><td>2.12</td><td>2.87*</td></tr><tr><td>QTL15</td><td>YSS</td><td>3</td><td>5.84</td><td>1.61</td><td>3.98***</td></tr></tbody></table>

<table><caption>**Table 5 continued**</caption><thead><tr><th>QTL</th><th>Trait</th><th>Chr</th><th>Var LE (%)</th><th>Var HE (%)</th><th>Z Diet</th></tr></thead><tbody><tr><td>QTL16</td><td>YSS</td><td>3</td><td>4.96</td><td>1.72</td><td>3.45***</td></tr><tr><td>QTL52</td><td>YSS</td><td>9</td><td>4.34</td><td>0.64</td><td>4.43***</td></tr><tr><td>QTL56</td><td>YSS</td><td>10</td><td>4.33</td><td>0.9</td><td>7.3***</td></tr><tr><td>QTL82</td><td>YSS</td><td>19</td><td>4.26</td><td>0.26</td><td>5.93***</td></tr><tr><td>QTL11</td><td>ESS</td><td>2</td><td>1.81</td><td>5.71</td><td>3.84***</td></tr><tr><td>QTL102</td><td>ESS</td><td>24</td><td>4.89</td><td>1.3</td><td>3.01**</td></tr><tr><td>QTL117</td><td>ESS</td><td>LGE22</td><td>0.06</td><td>4.46</td><td>6.14***</td></tr><tr><td>QTL4</td><td>ESSTIF</td><td>1</td><td>7.56</td><td>2.94</td><td>3.45***</td></tr><tr><td>QTL85</td><td>ESSTIF</td><td>20</td><td>3.02</td><td>1.19</td><td>4.38***</td></tr><tr><td>QTL46</td><td>HU</td><td>9</td><td>1.22</td><td>5.17</td><td>3.95***</td></tr><tr><td>QTL110</td><td>HU</td><td>27</td><td>1.73</td><td>4.49</td><td>3.55***</td></tr><tr><td>QTL19</td><td>EMTSP</td><td>3</td><td>2.75</td><td>4.92</td><td>2.77*</td></tr><tr><td>QTL87</td><td>EMTSP</td><td>20</td><td>1.24</td><td>3.37</td><td>4.04***</td></tr><tr><td>QTL100</td><td>EMTSP</td><td>24</td><td>1.46</td><td>3.91</td><td>2.84*</td></tr><tr><td>QTL62</td><td>YOLKIND</td><td>10</td><td>5.14</td><td>0.89</td><td>3.03**</td></tr><tr><td>QTL105</td><td>YOLKIND</td><td>26</td><td>4.19</td><td>0.64</td><td>3.55***</td></tr></tbody></table>

Chr chromosome that carries the QTL, var LE variance explained for the low energy diet, var HE variance explained for the high energy diet, LGE22 linkage group LGE22C19W28_E50C23, Z Diet Z test statistics for QTL × diet interaction

* P < 0.1, ** P < 0.05, *** P < 0.01

with diet, QTL1, QTL45 and QTL84 were detected for ESshape. QTL101, which showed an interaction with age at egg collection, explained a larger proportion of variance for ESshape for egg collection at 50 weeks of age (4.2 vs. 1.3 %) than at 70 weeks of age. None of the QTL detected for SLE showed an interaction with age at egg collection or diet. Thus, this trait does not seem to be sensitive to environment, or at least it was metabolically stable regardless of age at egg collection or diet.

### Egg shell color (ESC, RSS, LSS and YSS) traits

Among the 39 QTL detected for egg shell color, 10 showed an interaction with diet, one an interaction with age at egg collection and one an interaction with both age at egg collection and diet (see Additional files 1, 2, 3). The QTL that showed an interaction with diet were located on chromosomes 3, 9, 10, 12, 13 and 19, and nine of these 10 QTL had a greater effect for the LE diet than for the HE diet (Table 5). QTL117, which was located on LGE22C19W28_E50C23 and had AX-80867378 as the top SNP, showed an interaction with both diet and age at egg collection. It accounted for a larger proportion of variance for egg collection at 70 weeks of age and for the HE diet (Table 5). QTL84 was identified on chromosome 19 around 5 Mb and its effect varied with age at egg collection. For LSS, no QTL showed an interaction with either age at egg collection or diet, which suggests that it is not sensitive to variations in environment.
### Egg shell solidity (ESSTIF and ESS) traits

Among the 15 QTL detected for egg shell solidity, four showed an interaction with diet, three with age at egg collection and one with both diet and age at egg collection (Table 5). Among the QTL that showed an interaction with diet, QTL4 for ESSTIF and QTL102 for ESS accounted for a larger proportion of variance for the LE diet than for the HE diet, and QTL11 and QTL117 for ESS accounted for a larger proportion of variance for the HE diet than for the LE diet. The three QTL that showed an interaction with age at egg collection (QTL18, QTL72 and QTL103) were related to ESS and accounted for a larger proportion of variance for egg collection at 70 weeks of age. QTL85, detected for ESSTIF, showed an interaction with both diet and age and accounted for a

larger proportion of variance for the LE diet and for egg collection at 50 weeks of age. The same number of QTL was detected for ESSTIF and ESS but, for ESS, more QTL showed an interaction with diet or age, which suggests that ESS is more sensitive to environment than ESSTIF.

### Internal egg quality (YOLKIND, HU and EMTSP) traits

Among the 33 QTL detected for internal egg quality traits, five showed an interaction with diet, nine with age at egg collection and two with both diet and age at egg collection. Among the QTL that showed an interaction with diet only, three (QTL46, QTL87 and QTL110) accounted for a larger proportion of variance for the HE diet than for the LE diet. The two other QTL (QTL62 and QTL105), which were detected for YOLKIND, accounted

<table><caption>Table 6 Number and names of QTL detected on each chromosome and traits affected</caption><thead><tr><th>Chr</th><th>Nb QTL</th><th>QTL</th><th>Trait</th></tr></thead><tbody><tr><td>0</td><td>5</td><td>SNP5, SNP4, SNP2, SNP1, SNP3</td><td>EPR, EPR3, ESC, ESshape, EW, LSS, RSS</td></tr><tr><td>1</td><td>5</td><td>QTL1, QTL2, QTL4, QTL3, QTL5</td><td>ESshape, ESSTIF, YOLKIND</td></tr><tr><td>2</td><td>7</td><td>QTL7, QTL9, QTL10, QTL11, QTL8, QTL12, QTL6</td><td>EPR, EPR1, EPR2, EPR3, ESC, ESS, ESSTIF, LSS, RSS, YSS</td></tr><tr><td>3</td><td>10</td><td>QTL19, QTL14, QTL15, QTL20, QTL18, QTL17, QTL21, QTL13, QTL16, QTL22</td><td>EMTSP, EPR1, ESC, ESS, HU, LSS, RSS, YSS</td></tr><tr><td>4</td><td>4</td><td>QTL23, QTL26, QTL24, QTL25</td><td>EPR, EPR3, EW, HU</td></tr><tr><td>5</td><td>6</td><td>QTL28, QTL30, QTL31, QTL27, QTL32, QTL29</td><td>EPR1, ESC, ESshape, LSS, RSS</td></tr><tr><td>6</td><td>4</td><td>QTL35, QTL33, QTL34, QTL36</td><td>EW, HU, YOLKIND</td></tr><tr><td>7</td><td>4</td><td>QTL38, QTL39, QTL40, QTL37</td><td>EPR1, ESSTIF</td></tr><tr><td>8</td><td>3</td><td>QTL43, QTL42, QTL41</td><td>EMTSP, EPR1, ESSTIF</td></tr><tr><td>9</td><td>11</td><td>QTL44, QTL47, QTL46, QTL51, QTL48, QTL49, QTL53, QTL54, QTL50, QTL45, QTL52</td><td>EPR, EPR1, EPR2, EPR3, ESC, ESS, ESshape, HU, LSS, RSS, YSS</td></tr><tr><td>10</td><td>10</td><td>QTL60, QTL61, QTL58, QTL60-61, QTL59, QTL63, QTL55, QTL62, QTL56, QTL57</td><td>ESC, ESshape, EW, LSS, YOLKIND, YSS</td></tr><tr><td>11</td><td>5</td><td>QTL65, QTL68, QTL67, QTL66, QTL64</td><td>EPR, EPR2, EW, RSS, SLE</td></tr><tr><td>12</td><td>5</td><td>QTL70, QTL71, QTL73, QTL72, QTL69</td><td>EMTSP, EPR1, ESC, ESS, LSS, RSS</td></tr><tr><td>13</td><td>2</td><td>QTL75, QTL74</td><td>EPR1, EPR3, ESC</td></tr><tr><td>14</td><td>1</td><td>QTL76</td><td>EPR1</td></tr><tr><td>15</td><td>2</td><td>QTL77, QTL78</td><td>HU, YOLKIND</td></tr><tr><td>17</td><td>2</td><td>QTL79, QTL80</td><td>EPR, EPR1, ESC, LSS, RSS</td></tr><tr><td>18</td><td>1</td><td>QTL81</td><td>ESC, LSS, YSS</td></tr><tr><td>19</td><td>3</td><td>QTL82, QTL83, QTL84</td><td>ESC, ESshape, LSS, RSS, YSS</td></tr><tr><td>20</td><td>6</td><td>QTL87, QTL90, QTL85, QTL86, QTL88, QTL89</td><td>EMTSP, EPR, EPR1, EPR2, EPR3, ESSTIF, YOLKIND</td></tr><tr><td>21</td><td>4</td><td>QTL91, QTL94, QTL93, QTL92</td><td>EPR1, EW, HU, SLE, YOLKIND</td></tr><tr><td>22</td><td>2</td><td>QTL96, QTL95</td><td>ESSTIF, LSS</td></tr><tr><td>23</td><td>3</td><td>QTL99, QTL97, QTL98</td><td>EPR1, ESC, ESshape, LSS</td></tr><tr><td>24</td><td>5</td><td>QTL100, QTL102, QTL103, QTL101, QTL104</td><td>EMTSP, ESS, ESshape, ESSTIF, HU, SLE</td></tr><tr><td>26</td><td>5</td><td>QTL106, QTL109, QTL108, QTL107, QTL105</td><td>ESshape, EW, HU, SLE, YOLKIND</td></tr><tr><td>27</td><td>2</td><td>QTL110, QTL111</td><td>EPR1, HU</td></tr><tr><td>28</td><td>6</td><td>QTL115, QTL116, QTL114, QTL113-114, QTL112, QTL113</td><td>EPR1, EPR2, ESS, YOLKIND</td></tr><tr><td>LGE22</td><td>1</td><td>QTL117</td><td>ESS, ESshape, RSS, YOLKIND</td></tr><tr><td>LGE64</td><td>1</td><td>QTL118</td><td>ESshape, YSS</td></tr><tr><td>Z</td><td>8</td><td>QTL120, QTL123, QTL124, QTL125, QTL119, QTL122, QTL126, QTL121</td><td>EMTSP, ESC, ESshape, EW, LSS, RSS, SLE, YSS</td></tr></tbody></table>

Chr chromosome that carries the QTL, Nb QTL number of QTL detected on each chromosome
for a larger proportion of variance for the LE diet than for the HE diet. Among the QTL that showed an interaction with age, four (QTL21, QTL25, QTL111 and QTL113) accounted for a larger proportion of variance for egg collection at 50 weeks of age. The five other QTL (QTL17, QTL33, QTL78, QTL93 and QTL100) accounted for a larger proportion of variance for egg collection at 70 weeks of age. The two QTL that showed an interaction with both diet and age were detected for EMTSP and accounted for a larger proportion of variance for the HE diet and for egg collection at 50 weeks of age.

## Conclusions

This study detected 131 QTL for egg production and egg quality traits. These QTL were distributed across 27 chromosomes, two linkage groups (LGE22C19W28_E50C23 and LGE64) and a group of unassigned SNPs (Table 6). Among the 131 QTL detected, 60 showed a significant interaction with age at egg collection and/or diet (Table 5) although the average phenotypic performance varied only slightly with diet and age at collection (Table 1). This shows that laying hens have an in-built ability to adapt to their environment that probably involves different genetic pathways. These complex G × E interactions could have an effect on genetic selection, since the best candidates may differ depending on the environmental conditions in which the hens are reared. Performance of crossbred daughters was used to characterize pure line sires. This study pinpoints the existence of “unrobust” QTL, which raises the question which QTL are expressed in the commercial hybrids. To answer this, QTL detection based on the genotypes of commercial hybrids would be useful.

## Additional files

**Additional file 1.** QTL detected using the whole dataset to determine the genetic architecture of egg production and egg quality traits. This file gives the position of all the QTL detected using the whole dataset, with the top SNP corresponding to the SNP with the highest effect in the QTL region. The QTL is defined by the first (left) and last (right) SNPs that are 1 % significant at the chromosome level, respectively. Var (%) is the percentage of variance explained by the top SNP in the analysis with the whole dataset. Var LE(%) is the percentage of variance explained by the top SNP in the analysis with data for the low-energy diet only. Var HE(%) is the percentage of variance explained by the top SNP in the analysis with data for the high-energy diet only. Var 50(%) is the percentage of variance explained by the top SNP in the analysis with data for egg collection at 50 weeks only. Var 70(%) is the percentage of variance explained by the top SNP in the analysis with data for egg collection at 70 weeks only. Z Diet is the Z test statistics used to compare the two estimates calculated from the data for LE and HE diets. Z Age is the Z test statistics used to compare the two estimates calculated from the data for egg collection at 50 and 70 weeks of age. The difference was significant when P < 0.1 *, P < 0.05 ** and P < 0.01***.

**Additional file 2.** Additional QTL detected using the data for egg production and egg quality traits within diet. This file gives the position of the QTL detected with the dataset divided by diet, with the top SNP corresponding to the SNP with the highest effect in the QTL region. The QTL is defined by the first (left) and last (right) SNPs that are 1 % significant at the chromosome level, respectively. Var (%) is the percentage of variance explained by the top SNP in the analysis with the whole dataset. Var LE(%) is the percentage of variance explained by the top SNP in the analysis with data for the low-energy diet only. Var HE(%) is the percentage of variance explained by the top SNP in the analysis with data for the high-energy diet only. Z Diet is the Z test statistic used to compare the two estimates calculated from the data for LE and HE diets. The difference was significant when P < 0.1 *, P < 0.05 ** and P < 0.01***.

**Additional file 3.** Additional QTL detected using the data for egg quality traits according to age. This file gives the position of the QTL detected with the dataset divided by age, with the top SNP corresponding to the SNP with the highest effect in the QTL region. The QTL is defined by the first (left) and last (right) SNPs that are 1 % significant at the chromosome level, respectively. Var (%) is the percentage of variance explained by the top SNP in the analysis with the whole dataset. Var 50(%) is the percentage of variance explained by the top SNP in the analysis with data for egg collection at 50 weeks only. Var 70(%) is the percentage of variance explained by the top SNP in the analysis with data for egg collection at 70 weeks only. Z Age is the Z test statistic used to compare the two estimates calculated from the data for egg collection at 50 and 70 weeks of age. The difference was significant when P < 0.1 *, P < 0.05 ** and P < 0.01***.

## Authors' contributions

HR performed the QTL analyses. FH filtered the genotype data and contributed to QTL analyses. AVa and TB supervised animal management and production. CA supervised the phenotyping of animals. AVi supervised the genome scan. PD and HC assisted with management and analysis of genotyping data. HR drafted the manuscript. PLR supervised the QTL analyses and participated in drafting the manuscript. PLR and TB conceived the study. All authors contributed to the ideas and methods. All authors read and approved the final manuscript.

## Author details

¹ INRA, UMR1348 PEGASE, Domaine de La Prise, 35590 Saint-Gilles, France.

² Agrocampus Ouest, UMR1348 PEGASE, 65 Rue de Saint Brieuc, 35042 Rennes, France. ³ Novogen, Mauguérand, 22800 Le Foeil, France. ⁴ SYSAAF, INRA UR83 Recherches Avicoles, 37380 Nouzilly, France. ⁵ Zootests, Parc Technologique Du Zoopôle, 5 Rue Gabriel Calloet Kerbrat, 22440 Ploufragan, France. ⁶ INRA, UMR1388 GenPhySe, Auzeville BP52627, 31326 Castanet-Tolosan, France.

## Acknowledgements

This research project was supported by the French national research agency "ANR" within the framework of project ANR-10-GENOM_BTV-015 UtOplGe. HR is a PhD fellow supported by the Brittany region (France) and the INRA's Animal Genetics division.

## Competing interests

The authors declare that they have no competing interests.

Received: 7 April 2015 Accepted: 6 October 2015

Published online: 19 October 2015

## References

1. Chicken QTL database. Iowa State University 2003. [http://www.animal-genome.org/cgi-bin/QTldb/GG/index]. Accessed 13 Mar 2015.

2. Kranis A, Gheyas AA, Boschiero C, Turner F, Yu L, Smith S, et al. Development of a high density 600K SNP genotyping array for chicken. BMC Genomics. 2013;14:59.
3. El-Soda M, Malosetti M, Zwaan BJ, Koornneef M, Aarts MGM. Genotype × environment interaction QTL mapping in plants: lessons from Arabidopsis. *Trends Plant Sci.* 2014;19:390–8.

4. Atol Ontology. INRA 2012. [http://www.atol-ontology.com]. Accessed 11 Mar 2015.

5. Eisen EJ, Bohren BB, McKean HE. The Haugh unit as a measure of egg albumen quality. *Poult Sci.* 1962;41:1461–8.

6. Teyssèdre S, Elsen JM, Ricard A. Statistical distributions of test statistics used for quantitative trait association mapping in structured populations. *Genet Sel Evol.* 2012;44:32.

7. Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Harlow: Pearson, Prentice Hall; 2009.

8. Misztal I, Tsuruta S, Strabel T, Auvery B, Druet T, Lee DH. BLUPF90 and related programs (BGF90). In Proceedings of the 7th World Congress on Genetics Applied to Livestock Production. Montpellier; 2002. p. 19–23.

9. Müller BU, Stich B, Piepho H-P. A general method for controlling the genome-wide type I error rate in linkage and association mapping experiments in plants. *Heredity* (Edinb). 2011;106:825–31.

10. Ricard A, Filangi O, Elsen JM. GWAS Muller—Guide de l'utilisateur. Version 3.0.0. 2013. Accessed 29 Sept 2015.

11. Atzmon G, Blum S, Feldman M, Lavi U, Hillel J. Detection of agriculturally important QTLs in chickens and analysis of the factors affecting genotyping strategy. *Cytogenet Genome Res.* 2007;117:327–37.

12. Tuiskula-Haavisto M, de Koning DJ, Honkatukia M, Schulman NF, Mäki-Tanila A, Vilkki J. Quantitative trait loci with parent-of-origin effects in chicken. *Genet Res.* 2004;84:57–66.

13. Liu W, Li D, Liu J, Chen S, Qu L, Zheng J, et al. A genome-wide SNP scan reveals novel loci for egg production and quality traits in white leghorn and brown-egg dwarf layers. *PLoS One*. 2011;6:e28600.

14. Tuiskula-Haavisto M, Honkatukia M, Vilkki J, de Koning DJ, Schulman NF, Mäki-Tanila A. Mapping of quantitative trait loci affecting quality and production traits in egg layers. *Poult Sci.* 2002;81:919–27.

15. Wolc A, Arango J, Settar P, Fulton JE, O'Sullivan NP, Preisinger R, et al. Genome-wide association analysis and genetic architecture of egg weight and egg uniformity in layer chickens. *Anim Genet.* 2012;43:87–96.

16. Schreiweis MA, Hester PY, Settar P, Moody DE. Identification of quantitative trait loci associated with egg quality, egg production, and body weight in an F2 resource population of chickens. *Anim Genet.* 2006;37:106–12.

17. Wolc A, Arango J, Jankowski T, Dunn I, Settar P, Fulton JE, et al. Genome-wide association study for egg production and quality in layer chickens. *J Anim Breed Genet.* 2014;131:173–82.

18. Honkatukia M, Tuiskula-Haavisto M, Ahola V, Uimari P, Schmutz M, Preisinger RD, et al. Mapping of QTL affecting incidence of blood and meat inclusions in egg layers. *BMC Genet.* 2011;12:55.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

( ) **BioMed** Central