Skip to main content
Download PDF
  • Research Article
  • Open access
  • Published:

Identification of quantitative trait loci for body temperature, body weight, breast yield, and digestibility in an advanced intercross line of chickens under heat stress

Genetics Selection Evolution volume 47, Article number: 96 (2015) Cite this article

Abstract

Background

Losses in poultry production due to heat stress have considerable negative economic consequences. Previous studies in poultry have elucidated a genetic influence on response to heat. Using a unique chicken genetic resource, we identified genomic regions associated with body temperature (BT), body weight (BW), breast yield, and digestibility measured during heat stress. Identifying genes associated with a favorable response during high ambient temperature can facilitate genetic selection of heat-resilient chickens.

Methods

Generations F18 and F19 of a broiler (heat-susceptible) × Fayoumi (heat-resistant) advanced intercross line (AIL) were used to fine-map quantitative trait loci (QTL). Six hundred and thirty-one birds were exposed to daily heat cycles from 22 to 28 days of age, and phenotypes were measured before heat treatment, on the 1st day and after 1 week of heat treatment. BT was measured at these three phases and BW at pre-heat treatment and after 1 week of heat treatment. Breast muscle yield was calculated as the percentage of BW at day 28. Ileal feed digestibility was assayed from digesta collected from the ileum at day 28. Four hundred and sixty-eight AIL were genotyped using the 600 K Affymetrix chicken SNP (single nucleotide polymorphism) array. Trait heritabilities were estimated using an animal model. A genome-wide association study (GWAS) for these traits and changes in BT and BW was conducted using Bayesian analyses. Candidate genes were identified within 200-kb regions around SNPs with significant association signals.

Results

Heritabilities were low to moderate (0.03 to 0.35). We identified QTL for BT on Gallus gallus chromosome (GGA)14, 15, 26, and 27; BW on GGA1 to 8, 10, 14, and 21; dry matter digestibility on GGA19, 20 and 21; and QTL of very large effect for breast muscle yield on GGA1, 15, and 22 with a single 1-Mb window on GGA1 explaining more than 15 % of the genetic variation.

Conclusions

This is the first study to estimate heritabilities and perform GWAS using this AIL for traits measured during heat stress. Significant QTL as well as low to moderate heritabilities were found for each trait, and these QTL may facilitate selection for improved animal performance in hot climatic conditions.

Background

The climate is becoming increasingly warmer, according to the Intergovernmental Panel on Climate Change, and the global average temperature will continue to increase by 0.2 °C per decade. Heat stress in poultry impacts animal production and welfare and, in the poultry industry, it causes an estimated economic loss of $125 to 165 million in the US, with the broiler industry alone accounting for $58.1 million [1]. In 2007, an extreme heat wave in California resulted in more than 700,000 deaths in poultry [2] and in 2009, over 1.5 million layer hens died during a summer heat wave (National Oceanic and Atmospheric Association).

Production losses due to heat stress may result from mortality, reduced body weight, reduced egg production, reduced feed intake, and higher feed to gain ratio. A recent study on broilers that were exposed to chronic heat stress from 1 to 42 days of age showed a reduced body weight (32.6 %), increased feed conversion ratio (25.6 %), and reduced feed intake (16.4 %) [3], and another study using shorter periods of heat stress on younger birds, from 2 to 4 weeks of age showed a reduced feed intake by 14 % [4]. In a paired feed study, genetically lean broilers that were exposed to chronic heat stress, from hatch to 9 weeks of age showed increased weight gain and feed efficiency compared to less lean counterparts, which supports the hypothesis that increased fat accretion is inversely related to thermo-regulation [5]. Previous studies in poultry have elucidated a genetic influence on response to heat. Layers that were divergently selected for heat tolerance displayed differences in survivability during increased heat conditions [6]. Significant differences in production traits have been found between a commercially fast growing chicken line and a local chicken breed from China during heat stress [7]. Microsatellites were used to identify quantitative trait loci (QTL) for several traits measured during heat stress in a Japanese quail F2 intercross including body weight (BW), feed intake, and body temperature (BT) [8]. Because in poultry, response to heat stress involves a genetic component, it is possible to use genomic selection for heat tolerance, which will increase accuracies and response to selection [9]. To increase our understanding of the genetic influence on response to heat stress in chickens, we used the F18 and F19 generations of a broiler (heat-susceptible) × Fayoumi (heat-resistant) advanced intercross line (AIL) and an environmentally-controlled experiment to identify genomic markers related to response to high ambient temperatures.

Chickens of this AIL were exposed to high ambient temperatures for 7 days during which BW, BT, breast yield, and digestibility were measured. These traits, as well as the changes in BT and BW due to heat treatment, were used for genome-wide association studies (GWAS). The genes and markers associated with thermal tolerance can help elucidate the genetic architecture of traits involved in heat stress and, subsequently, be used to breed more heat-resilient chickens.

Methods

Chicken lines

All animal experiments were approved by the Institutional Animal Care and Use Committee at Iowa State University: Log #4-11-7128-G. We used the F18 and F19 generations of an AIL between two highly divergent chicken lines for thermo-tolerance, i.e. it was created by crossing a single broiler sire to six highly inbred Fayoumi dams [10]. Although this population has limited variability due to the initial mating, the broiler sire was characterized by the offsprings’ phenotypic means and variances of body composition phenotypes, which showed that it was representative of the entire broiler population [10]. We hypothesize that the highly inbred Fayoumi breed became fixed for alleles that had the highest frequency in the founder line. Thus, this population is a powerful resource to identify QTL. The broiler breed has been commercially selected for muscle accretion, whereas the Fayoumi breed has not undergone commercial selection. The Fayoumi breed originated in the Fayoum region of Egypt that is characterized by a high-temperature climate and thus this breed has undergone natural selection for tolerance to heat. Birds were reared on floor pens with wood shavings and had ad libitum access to water and corn-soy feed that met all NRC requirements for the duration of the study [11].

Heat stress experimental design

We used birds from two generations with each generation producing two hatches. Six hundred and thirty-one birds from the four hatches were used for independent heat stress experiments (four replicates). Six hundred animals from 17 sire families were phenotyped for BT, BW, and breast yield. Digestibility measurements were available for 461 animals from 14 sire families. At 17 days of age, birds were transferred to environmentally-controlled chambers and acclimated for 5 days. There were four chambers, each containing six pens, per replicate. From days 22 to 28 of age, the chambers were heated to 35 °C for 7 h per day and remained at 25 °C at all other times.

Phenotypic measurements

Cloacal BT was measured by inserting a digital thermometer approximately 2.5 cm into the cloaca on days 20, 22, and 28 of age. The precision of the digital thermometer was 0.1 °C. BW was measured using a digital scale on days 21 and 28 of age. Breast yield (%) was determined by weighing one half of the pectoralis major and minor muscle, multiplying this value by two, and then dividing it by the total BW on day 28. Dry matter digestibility was measured as described in [12]. Briefly, dry matter was determined by drying ileal and feed samples for 24 h at 110 °C. Titanium dioxide was used as a marker for both ileal and feed samples and was analysed as described in [13]. Dry matter digestibility was calculated by the following equation [12]:

$$\left [ {\frac{{{\rm \%\, Diet \,DM} - ( {\frac{{\rm \%\, Fecal \,DM \; \times \% \,Diet \,Ti}}{{\rm \% \,Fecal \,Ti}}} )}}{{\rm \% \,Diet \,DM}}}\right ]\; \times \;100,$$

where DM is dry matter and Ti is titanium dioxide. Dry matter digestibility was log-transformed to obtain a normal distribution of the data and the transformed data was used for all downstream analyses.

DNA isolation and genotyping

Blood was collected from the wing vein by using an EDTA-coated syringe and needle, and then stored at −20 °C. DNA was extracted using a salting out method. Briefly, whole blood was incubated with lysis buffer containing proteinase K. Proteins were precipitated out using 5 M NaCl and the supernatant was recovered. 70 % ethanol was added to the supernatant to precipitate DNA. DNA isolated from 468 AIL, six broiler, and six Fayoumi chickens was genotyped using the Affymetrix 600 K chicken SNP (single nucleotide polymorphism) axiom array [14] by GeneSeek Inc., Lincoln, NE. SNP chromosomal locations were based on the Gallus_gallus_4.0 assembly through Ensembl.

Statistical analyses

Means, standard errors, fixed effects, and covariates for the GWAS analyses were calculated based on ANOVA (analysis of variance) estimates, and significant terms were included as fixed effects with a P value less than 0.05 using JMP statistical software [15]. Heritability was estimated using a single-trait animal model in ASReml [16]. For all phenotypes, sex was fitted as a fixed effect, while chamber nested within replicate and within animal were fitted as random effects. To estimate BW21 heritability, dam was fitted as a random effect. For the GWAS for BT, the closest BW measurement in days was fitted as a covariate.

Genotyping console (Affymetrix) software was used to obtain genotyping calls and to perform quality control based on whole animal DishQC score ≥0.7. The SNPolisher (Affymetrix) R package was used to perform quality control of individual SNPs for all the animals that passed the DishQC criterion. For SNP genotypes to be included in the analysis, SNP call rate had to be greater than 95 % and minor allele frequency (MAF) higher than 5 %.

GWAS for phenotypic traits with SNP genotypes was done using GenSel software [17]. Bayes B, which fits all SNPs simultaneously as random effects, was used for the analysis. The following mixed model was used for the GWAS:

$${\mathbf{y = Xb}}{ + }\sum\limits_{\text{j}}^{\text{k}} {{\mathbf{z}}_{\text{j}} \upalpha_{\text{j}} \delta_{\text{j}} { + } {\upvarepsilon} },$$

where y is a vector of phenotypes, X is an incidence matrix to account for fixed effects on phenotypes, b is a vector of fixed effects, z j is a vector of genotypes for SNP j based on the number of B alleles (−10, 0, +10, or the average of the genotypes at SNP j), αj is the allele substitution effect for SNP j, δj is a parameter that indicates whether SNP j was included in the Markov chain Monte Carlo (MCMC) chain, and ε is the error associated with the analysis. For one analysis per trait, a total of 41,000 MCMC iterations were completed and the first 1000 iterations were discarded.

SNPs were split into 1001 non-overlapping 1-Mb windows across the genome. Thus, the windows that have the SNP, which is most frequently included in the MCMC iterations (post-burn-in), are predicted to have an effect on the phenotype. δj was set so that π = 0.9978 to avoid fitting more SNPs than the number of animals in a given iteration. Using a true infinitesimal model, each window is expected to explain 0.1 % (100 %/1001) of the genetic variation; therefore, a 1-Mb window was considered significant if it explained more than 0.5 % of the total genetic variation.

Candidate genes

For each trait, the window explaining the largest percentage of genetic variation was investigated, and, within this window, the SNP that was most frequently included in the model was identified. Then, all annotated genes within 200 kb (100 kb upstream and 100 kb downstream) of that SNP were identified using ENSEMBL biomart [18]. We chose a 200-kb window based on the average linkage disequilibrium (LD) in commercial broiler populations i.e., less than 1 cM on average [19], and on the fact that the chicken genome contains 250 kb per cM on average [20]. In the F18 and F19 AIL chickens, LD was expected to cover a shorter distance than in the commercial broiler population because of their unique population structure and the large number of generations in which recombination occurred.

Results

Phenotypic measurements and heritabilities

Phenotypic means, standard errors, ranges, and heritabilities are in Table 1. BT measurements had low heritabilities that ranged from 0.03 to 0.11. Changes in BT after acute heat (1 day) and chronic heat (7 days) treatments had low heritabilities of 0.03 with large standard errors and were not statistically different from 0. BW measurements had moderate heritabilities that ranged from 0.15 to 0.35. Breast yield and dry matter digestibility, both measured on day 28, had moderate heritabilities of 0.15 and 0.33, respectively.

Table 1 Phenotypic means and heritabilities (h2)
Full size table

Genotypes

Of the 480 birds that were genotyped, 458 AIL and all 12 parental line birds passed the whole animal DishQC criterion. Of the 580,961 SNPs on the array, filtering based on a SNP call rate greater than 95 % removed a small proportion, i.e. 59,789 SNPs, while filtering based on MAF removed a much larger proportion, i.e. 311,055 SNPs, thus 210,117 SNPs remained for subsequent analyses.

GWAS

The detailed results for each window that explained a significant percentage of the genetic variation (>0.5 %) and the SNP within each window that had the highest effect on the trait are in Table 2. To increase clarity, significant consecutive windows were designated as a single QTL. In total, 35 QTL were identified across all traits and measurement phases.

Table 2 Identified windows that explain a significant percentage (≥0.5 %) of the genetic variance
Full size table

Eight QTL were identified for BT phenotypes (Fig. 1a–d), i.e. (1) four for BT20 with two on GGA1 (GGA for Gallus gallus chromosome), one on GGA14, and one on GGA15; (2) two for BT28 with one on GGA15 and one on GGA26; and (3) two for BT28-20 (change in BT measured at pre-heat treatment and after 1 week of heat treatment) with one on GGA14 and one on GGA27. No QTL was found for BT22-20 (change in BT measured at pre-heat treatment and after 1 day of heat treatment). QTL co-localizations were identified on GGA14 for BT20 and BT28-20 and on GGA15 for BT20 and BT28.

Fig. 1
figure 1

Genome-wide plot of percentage of genetic variance for traits measured during heat stress. Traits were measured before heat treatment (day 20 or 21), during acute heat treatment (day 22) and chronic heat treatment (day 28), and the differentials between trait measurements due to 1 day of heat treatment (day 22–20) and to 7 days of heat treatment (day 28–20 or 28–21) were calculated. Only traits that reached significance in the GWAS (≥0.05 % of the genetic variation) are displayed. Plots for body temperature (BT) measured on days 20, 22, 28, and the differential 28-20 (a–d); plots for body weight (BW) measured on days 21, 28, and the differential 28–21 (e–g); plot for digestibility, measured from ileal content, on day 28 (H); plot for  % of breast weight (i), calculated from  % of total body weight, and measured on day 28. Results show the percentage of genetic variance that is explained by each non-overlapping 1-Mb window, labelled by the index number of the windows coloured and ordered by chromosome (1–27, and Z)

Full size image

Twenty-one QTL for BW phenotypes were identified, i.e. (1) seven for BW21 with one each on GGA1, 2, 4, 6, 7, 8, and 14; (2) nine for BW28 with two on GGA1, two on GGA5, and one each on GGA2, 3, 4, 6, and 26; and (3) five for BW28-21 (change in BW measured at pre-heat treatment and after 1 week of heat treatment) with one each on GGA1, 4, 6, 10, and 21 (Fig. 1e–g). QTL for all BW phenotypes co-localized on GGA6 in a region containing four adjacent 1-Mb windows. QTL also co-localized on GGA1 for BW21 and BW28 and on GGA4 for BW28 and BW28-21.

Three QTL were identified for digestibility with one each on GGA19, 20 and 21 (Fig. 1h). QTL co-localized for digestibility and BW28-21 on GGA21.

For breast yield, an economically important trait, three QTL were identified with one each on GGA1, 15, and 22 (Fig. 1i). The QTL on GGA15 included five adjacent 1-Mb windows and cumulatively accounted for 24.5 % of the genetic variation. The most significant single 1-Mb window in this region accounted for 15.4 % of the genetic variation. QTL co-localization was not identified between breast yield and any of the other traits measured in the current study.

Candidate genes

For each trait, positional candidate genes were identified within a 200 kb region i.e. 100 kb upstream and 100 kb downstream of the SNP with the highest effect within the 1-Mb window that explained the highest percentage of genetic variation. Fifty annotated genes were identified (see Table 3). For BW measurements, five, four, and one positional candidate genes for BW21, BW28, and BW28-21, respectively, were found among the 18, 24, and six genes that were located within the corresponding 1-Mb windows. For BT measurements, seven, seven and 11 positional candidate genes for BT20, BT28, and BT28-20, respectively, were found among the 48, 38, and 29 genes that were located within the corresponding 1-Mb windows. For dry matter digestibility and breast yield, six and four positional candidate genes, respectively, were identified among the 19 and 10 genes that were located within the corresponding 1-Mb windows.

Table 3 Positional candidate genes categorized by function for windows explaining the largest percentage of genetic variation
Full size table

Discussion

The aims of this study were to identify and estimate the effect of QTL, and to identify positional candidate genes, for BT, BW, dry matter digestibility, and breast yield using a novel AIL under heat stress and a 600 K SNP panel for genotyping.

Population used

Previous generations of this AIL were used for several QTL mapping studies and allowed the identification of many QTL including 257 QTL for growth and body composition [21–25], 93 for skeletal integrity [26], 51 for metabolic traits [27], and 12 for response to Salmonella enteritidis challenge [28–30]. Therefore, collectively, a wide range of traits has been associated with a large number of QTL in previous generations of this AIL. The continued erosion of LD in this population over subsequent generations, combined with the availability of more dense SNP panels, creates a unique opportunity to more finely map the location of QTL that are in LD with a causal mutation.

Phenotypes and heritabilities

Phenotypic measurements for BT and BW consisted of repeated measures on individual birds. This allowed us to use both absolute measures and the differences between measures carried out before (pre-heat) and during heat treatments. Since measurements of dry matter digestibility and breast yield required euthanization, they were only performed after 7 days of heat treatment, on day 28.

A previous study reported a significant correlation between BT and survival in chickens during heat stress [31], which suggested that selection for BT during heat stress has potential to reduce mortality. In our study, heritabilities were low (0.03–0.11) for BT and higher (0.10–0.11) for absolute BT than that previously estimated for a broiler line i.e. 0.05 [32]. Heritabilities for changes in BT from pre-heat to acute and chronic heat conditions were both low i.e. 0.03, which could be due to the low precision and large variation of the measurements. This indicates that it will be challenging to genetically select for resistance to BT change during heat stress. More precise methods of BT measurement (e.g., infrared thermography) should be explored.

Heritabilities for BW were low to moderate (0.15–0.35), which agree with the heritabilities of 0.4 to 0.6 previously reported for BW in broiler lines [33]. Heritability for breast yield was moderate i.e. 0.15 and agreed with previously reported estimates [34].

In this study, dry matter digestibility was measured using a titanium oxide marker to calculate dry matter in the feed and ileal contents. We estimated a heritability of 0.33 for digestibility, which is similar to that (0.33 to 0.47) reported for broilers that were fed a wheat-based diet [35]. Many previous studies on feed conversion ratio in chickens reported moderate heritabilities [36–38]. Feed intake and conversion ratios, which are associated with digestibility, are arguably the most costly impacts of heat stress. Because the heritability estimated for digestibility during heat stress is moderate, there may be potential for improvement of this trait via selection.

Generally, heritabilities estimated for most traits in our study are lower than those previously reported, which is likely due to lower genetic variation within the population studied.

Genome-wide association study

To date, seven QTL for BT have been reported on GGA2, 3, 4, 5, 6 and 11 [39–41], but none overlapped with those detected here. This absence of QTL overlap may be due to differences in experimental protocols since, in these previous studies [39–41], the traits that were measured were response to disease challenge and resting BT between lines selected for growth or fat accretion. In addition, these QTL may be population-specific.

The two QTL which co-localized for BT20 and BT28-20 and for BT20 and BT28 were near QTL for hematocrit on GGA14 [42] and GGA15 [39], respectively. One mechanism to regulate BT during periods of heat is to increase blood flow towards the body surface [43], and it has also been well documented that panting behaviour occurs under high temperatures [44]; both mechanisms result in changes in the blood system. Moreover, the co-localized QTL for BT on GGA15 was located near a QTL for corticosterone that was measured in response to manual restraint [45]. This co-localization suggests that it may be a good candidate for further investigation of the pleiotropic response to stress.

Heat stress specific QTL were identified for BT28 and BT28-20 on GGA26 and 27, respectively. Both regions present considerable overlap with previously reported QTL for growth [23, 24, 40, 46–49]. The large overlap between BT and QTL for growth is not surprising given the highly negative correlation between BT and growth during heat stress in chickens [50], although we attempted to account for this relationship by fitting BW as a co-variate in the GWAS analysis.

The co-localized QTL for all BW measurements on GGA6 was near a previously reported QTL for growth in many different chicken populations, including a broiler × layer cross [48, 49, 51, 52], White Plymouth Rock, New Hampshire and White Leghorn chickens [53], a commercial broiler line [54], high and low growth broiler lines [40], white leghorn x red jungle fowl [55], and the F2 broiler × Fayoumi generation used in the current study [23]. The QTL that we detected on GGA6 is confirmed by previously reported QTL in this region in a wide range of chicken populations, which suggests a conserved QTL, and supports our results. This region explained a relatively large percentage of the genetic variation for BW21 (3.0 %), BW28 (8.1 %), and BW28-21 (4.6 %), which confirms the importance of this QTL.

The co-localized QTL for BW21 and BW28 on GGA1 were also localized near QTL for growth that were previously reported in a large range of chicken populations including crosses between Silkie fowl and Cornish broiler [46], White Recessive Rock and Xinghua chicken [56], broiler and layer [49], and Shamo and White Plymouth Rock [57]. As observed for the co-localized QTL on GGA6, this region on GGA1 is probably highly conserved given the variety of populations for which QTL have been reported near this region.

Several QTL for average daily gain are located near the region where QTL for BW28 and BW28-21 co-localize on GGA4, in two chicken populations including a Silky fowl × White Plymouth rock cross [58], and a broiler × layer cross [49]. Again, this overlap between the QTL detected in our study and previously reported QTL supports our results.

Although it is interesting to discuss overlapping QTL for measurements performed at different phases, it is also relevant to examine the QTL that were identified specifically for traits measured during heat stress i.e. BW28 on GGA26 and BW28-21 on GGA10 and 21. Indeed, the QTL for BW28 on GGA26 was located within the 2-Mb QTL region for BT28 and as discussed above, many QTL for growth have been identified in this region. Similarly, QTL for growth have also been reported in the QTL regions on GGA10 and GGA21 that we detected here.

Although feed represents a large proportion of poultry production costs (51.8 % in 2001 and 68.7 % in 2008) and these costs continue to increase due to the increased demands for grain in other industries such as ethanol for fuel [59], very few QTL related to feed use have been identified (37 of the 4795 QTL listed in http://www.animalgenome.org). In a recent study that evaluated layer hens for feed intake and feed use and performed a GWAS using the 600 K array, eight QTL were identified [60], but none co-localized with the QTL reported here. The three QTL identified for digestibility were all located near previously reported QTL for growth-related traits. Very recently, Mignon-Grasteau et al. [61] identified a QTL for dry excreta weight on GGA19 within the same interval as that reported here for digestibility. The fact that the QTL regions for these two feed-use traits overlap provides evidence that this region on GGA19 is indeed a true QTL and should be further investigated given the economic importance of these traits to the poultry industry.

A strong QTL for breast yield was identified on GGA1. This region contains a large number of QTL related to abdominal fat and growth traits that were detected across diverse chicken populations (http://www.animalgenome.org). Surprisingly, no QTL for breast muscle has been reported in this region even in studies on previous generations of the same AIL. Thus, we suggest that this QTL may be specific to breast muscle growth during heat stress conditions. Furthermore, the QTL for breast yield that we detected on GGA15 and 22 overlap with previously reported QTL for breast muscle on GGA15 [40] and GGA22 [46]. The region on GGA1 warrants further investigation as a QTL specific to heat stress.

Candidate genes

All positional candidate genes were identified for each trait within 200 kb of the most significant SNP. Cellular response to heat stress has been extensively reviewed and involves a range of biological mechanisms, i.e. inhibition of DNA synthesis, transcription, translation, cell cycle arrest, denaturation of proteins, enhanced degradation of proteins by ubiquitin and lysosomal pathways, disruption of the integrity of the cytoskeleton [62] and increased apoptosis [63]. In addition, heat stress induces metabolic changes and increased intracellular ion concentrations. Previously, Coble et al. [64] observed that, in broiler chickens, heat stress induced transcriptional changes and Morimoto [65] reported an increased expression of heat shock protein genes in response to heat stress. Heat shock proteins form an evolutionarily conserved family across all multicellular organisms [66].

Genes involved in the disruption of DNA synthesis, transcription, RNA processing, and translation were identified near the QTL for all traits analyzed in this study except for BW28 and BW28-21. We identified three genes that code for separate subunits of the mediator complex (MED 9, 15, and 31). The mediator complex is required for the regulation of eukaryotic RNA polymerase II transcripts [67]. In yeast, the mediator complex interacts directly with heat shock proteins and serves as a bridge between heat stress and transcriptional regulation of heat shock related genes [68].

Disruption of progression through the cell cycle and resulting apoptosis occur during cellular stress. We identified five genes that are involved in cell cycle progression that were located near some of the QTL detected in our study for all BT measurements and digestibility. Among the cell cycle checkpoints, two crucial checkpoints, between G1/S and G2/M transitions, are arrested in response to heat stress [69]. Apoptosis is induced during extreme stress conditions. In this study, we identified genes involved in apoptosis near QTL for both digestibility and breast yield.

Other categories of functions that were associated with the candidate genes detected in our study include glucose regulation, disruption of the cytoskeleton, free radical damage, and blood vessel development. The ACE gene, involved in blood vessel development, was located near the QTL for BT28-20, thus we hypothesize that it may play a role in reducing BT during periods of heat stress by enhancing blood flow to the body surface [43].

Genes involved in membrane permeability and changes in cellular ion concentrations were located near the QTL for BT, BW, and digestibility. In chickens, Ait-Boulahsen et al. [70] showed that Na+, K+, and Cl− plasma levels increase in response to heat stress, which can have an effect on the endocrine system, for example as a secondary messenger, and subsequently on stress response. Two genes of this functional category were identified near the QTL for digestibility i.e. SLC13A5 that encodes a citrate transporter and PITPNM3 that encodes a calcium ion binding protein; both these genes are involved in ion movement, which is impacted by heat stress.

Genes related to the immune system were identified near the QTL for BT20 and BT28, and breast yield. It has been shown that, compared to animals not exposed to high ambient temperatures, laying hens that are exposed to cyclic heat stress have decreased T-cell and B-cell proliferations and antibody titre to sheep red blood cells and an increased total white blood cell count [71], which supports the hypothesis that the immune function is disrupted during heat stress. We identified the MIF and DDT genes near the QTL for BT22 and BT28, respectively. These genes function as proinflammatory cytokines involved in the immune response [72]. One of the hallmarks of inflammation is to increase BT. If inflammation can be suppressed in birds subjected to heat stress, this might decrease the negative impact of high ambient temperatures. The SOCS2 gene that encodes a suppressor of cytokine signalling was identified near the QTL for breast yield and may be a good candidate gene for future studies on the mechanisms that influence breast muscle yield in chickens.

Candidate genes with a role in cell signalling were identified near the QTL for BT and BW. Cell signalling increases during response to stress. One gene of particular interest in cell signalling is the MAP3K3 gene that was found near the QTL for BT28-20 since that the MAPK signalling pathway is known to be involved in the cellular response to stress [73].

Breast muscle yield is an extremely important trait because of its economic impact in the broiler industry. The 1-Mb window with the largest effect that we identified here was for breast yield and explained 15.4 % of the total genetic variation. The favourable allele of the SNP with the largest effect within the window that explained the largest proportion of genetic variation was fixed in the Fayoumi line but was also segregating in the broiler line; thus, it was not possible to determine which line contributed the favourable allele. The best candidate gene in this region is SOCS2, which has a role in suppressing cytokine signalling. We hypothesize that the effect of this QTL on muscle accretion may be heat-specific because no other QTL for breast yield was identified in this region.

Conclusions

SNPs were identified for BT, BW, digestibility, and breast yield in a unique chicken AIL measured under heat stress. A major QTL for breast yield under heat stress explained more than 24 % of the genetic variation. Exploiting this information for genomic selection to breed heat-tolerant chickens is feasible. The QTL regions that we identified contain many genes with functions that suggest a role in response to heat stress and, thus, these genes are both positional and functional candidates.

References

  1. St-Pierre N, Cobanov B, Schnitkey G. Economic losses from heat stress by US livestock industries. J DairySci. 2003;86:E52–77.

    Article  Google Scholar 

  2. Cershunov A, Cayan D. Recent increase in California heat waves: July 2006 and the last six decades. California Energy Commission, PIER Energy-Related Environmental Research. PIER Report to be released by September; 2008. p. 2008.

  3. Sohail MU, Hume ME, Byrd JA, Nisbet DJ, Ijaz A, Sohail A, et al. Effect of supplementation of prebiotic mannan-oligosaccharides and probiotic mixture on growth performance of broilers subjected to chronic heat stress. Poult Sci. 2012;91:2235–40.

    Article  CAS  PubMed  Google Scholar 

  4. Geraert PA, Padilha JC, Guillaumin S. Metabolic and endocrine changes induced by chronic heatexposure in broiler chickens: growth performance, body composition and energy retention. Br J Nutr. 1996;75:195–204.

    CAS  PubMed  Google Scholar 

  5. Geraert PA, Guillaumin S, Leclercq B. Are genetically lean broilers more resistant to hot climate? Br Poult Sci. 1993;34:643–53.

    Article  CAS  PubMed  Google Scholar 

  6. Wilson HR, Wilcox CJ, Voitle RA, Baird CD, Dorminey RW. Characteristics of White Leghorn chickens selected for heat tolerance. Poult Sci. 1975;54:126–30.

    Article  CAS  PubMed  Google Scholar 

  7. Lu Q, Wen J, Zhang H. Effect of chronic heat exposure on fat deposition and meat quality in two genetic types of chicken. Poult Sci. 2007;86:1059–64.

    Article  CAS  PubMed  Google Scholar 

  8. Minvielle F, Kayang BB, Inoue-Murayama M, Miwa M, Vignal A, Gourichon D, et al. Microsatellite mapping of QTL affecting growth, feed consumption, egg production, tonic immobility and body temperature of Japanese quail. BMC Genomics. 2005;6:87.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Goddard M. Genomic selection: prediction of accuracy and maximisation of long term response. Genetica. 2009;136:245–57.

    Article  PubMed  Google Scholar 

  10. Deeb N, Lamont SJ. Genetic architecture of growth and body composition in unique chicken populations. J Hered. 2002;93:107–18.

    Article  CAS  PubMed  Google Scholar 

  11. National Research Council. Nutrient requirements of poultry. 9th ed. Washington: National Academy Press; 1994.

    Google Scholar 

  12. Bolek KJ. Mitigation strategies to ameliorate acute and chronic heat stress utilizing supplemental methionine or embryonic thermal conditioning in chickens. MSc Thesis; 2013.

  13. Leone JL. Collaborative study of the quantitative determination of titanium dioxide in cheese. J Assoc Off Anal Chem. 1973;56:535–7.

    CAS  PubMed  Google Scholar 

  14. Kranis A, Gheyas AA, Boschiero C, Turner F, Yu L, Smith S, et al. Development of a high density 600 K SNP genotyping array for chicken. BMC Genomics. 2013;14:59.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. JMP®, Version 11.2.1. SAS Institute Inc., Cary, NC, 1989–2007.

  16. Gilmour AR, Gogel B, Cullis B, Thompson R. ASReml user guide release 3.0. Hemel Hempstead: VSN International Ltd; 2009.

    Google Scholar 

  17. Fernando R, Garrick DJ. GenSel: user manual for a portfolio of genomic selection related analyses. Ames: Iowa State University: Animal Breeding and Genetics; 2008.

    Google Scholar 

  18. Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Res. 2015;43:D662–9.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Andreescu C, Avendano S, Brown SR, Hassen A, Lamont SJ, Dekkers JC. Linkage disequilibrium in related breeding lines of chickens. Genetics. 2007;177:2161–9.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Cheng HH, Lamont SJ. Genetics of disease resistance. Diseases of poultry, vol. 12th ed. Oxford: Blackwell Publishing Ltd; 2011.

    Google Scholar 

  21. Deeb N, Lamont S. Use of a novel outbred by inbred F1 cross to detect genetic markers for growth. Anim Genet. 2003;34(3):205–12.

    Article  CAS  PubMed  Google Scholar 

  22. Li H, Deeb N, Zhou H, Mitchell AD, Ashwell CM, Lamont SJ. Chicken quantitative trait loci for growth and body composition associated with transforming growth factor-beta genes. Poult Sci. 2003;82:347–56.

    Article  CAS  PubMed  Google Scholar 

  23. Zhou H, Deeb N, Evock-Clover C, Ashwell C, Lamont S. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. I. Growth and average daily gain. Poult Sci. 2006;85(10):1700–11.

    Article  CAS  PubMed  Google Scholar 

  24. Zhou H, Deeb N, Evock-Clover CM, Ashwell CM, Lamont SJ. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. I. Growth and average daily gain. Poult Sci. 2006;85:1700–11.

    Article  CAS  PubMed  Google Scholar 

  25. Abasht B, Lamont SJ. Genome-wide association analysis reveals cryptic alleles as an important factor in heterosis for fatness in chicken F2 population. Anim Genet. 2007;38:491–8.

    Article  CAS  PubMed  Google Scholar 

  26. Zhou H, Deeb N, Evock-Clover CM, Mitchell AD, Ashwell CM, Lamont SJ. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. III. Skeletal integrity. Poult Sci. 2007;86:255–66.

    Article  CAS  PubMed  Google Scholar 

  27. Zhou H, Evock-Clover CM, McMurtry JP, Ashwell CM, Lamont SJ. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. IV. Metabolic traits. Poult Sci. 2007;86:267–76.

    Article  CAS  PubMed  Google Scholar 

  28. Kaiser MG, Lamont SJ. Microsatellites linked to Salmonella enterica serovar enteritidis burden in spleen and cecal content of young F1 broiler-cross chicks. Poult Sci. 2002;81:657–63.

    Article  CAS  PubMed  Google Scholar 

  29. Kaiser MG, Deeb N, Lamont SJ. Microsatellite markers linked to Salmonella enterica serovar enteritidis vaccine response in young F1 broiler-cross chicks. Poult Sci. 2002;81:193–201.

    Article  CAS  PubMed  Google Scholar 

  30. Liu W, Lamont SJ. Candidate gene approach: potentional association of caspase-1, inhibitor of apoptosis protein-1, and prosaposin gene polymorphisms with response to Salmonella enteritidis challenge or vaccination in young chicks. Anim Biotechnol. 2003;14:61–76.

    Article  CAS  PubMed  Google Scholar 

  31. Chen XY, Wei PP, Xu SY, Geng ZY, Jiang RS. Rectal temperature as an indicator for heat tolerance in chickens. Anim Sci J. 2013;84:737–9.

    Article  PubMed  Google Scholar 

  32. el-Gendy E, Washburn KW. Genetic variation in body temperature and its response to short-term acute heat stress in broilers. Poult Sci. 1995;74:225–30.

    Article  CAS  PubMed  Google Scholar 

  33. Chambers JR. Genetics of growth and meat production in chickens. In: Crawford RD, editor. Poultry breeding and genetics. Amsterdam: Elsevier; 1990. p. 599–643.

    Google Scholar 

  34. Le Bihan-Duval E, Mignon-Grasteau S, Millet N, Beaumont C. Genetic analysis of a selection experiment on increased body weight and breast muscle weight as well as on limited abdominal fat weight. Br Poult Sci. 1998;39:346–53.

    Article  PubMed  Google Scholar 

  35. Mignon-Grasteau S, Muley N, Bastianelli D, Gomez J, Péron A, Sellier N, et al. Heritability of digestibilities and divergent selection for digestion ability in growing chicks fed a wheat diet. Poult Sci. 2004;83:860–7.

    Article  CAS  PubMed  Google Scholar 

  36. Aggrey SE, Karnuah AB, Sebastian B, Anthony NB. Genetic properties of feed efficiency parameters in meat-type chickens. Genet Sel Evol. 2010;42:25.

    Article  PubMed Central  PubMed  Google Scholar 

  37. van Kaam JB, Groenen MA, Bovenhuis H, Veenendaal A, Vereijken AL, Van Arendonk JA. Whole genome scan in chickens for quantitative trait loci affecting growth and feed efficiency. Poult Sci. 1999;78:15–23.

    Article  PubMed  Google Scholar 

  38. Thomas CH, Blow WL, Cockerham CC, Glazener EW. The heritability of body weight, gain, feed consumption, and feed conversion in broilers. Poult Sci. 1958;37:862–9.

    Article  Google Scholar 

  39. Pinard-van der Laan MH, Bed’Hom B, Coville JL, Pitel F, Feve K, Leroux S, et al. Microsatellite mapping of QTLs affecting resistance to coccidiosis (Eimeria tenella) in a Fayoumi × White Leghorn cross. BMC Genomics. 2009;10:31.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Nadaf J, Pitel F, Gilbert H, Duclos MJ, Vignoles F, Beaumont C, et al. QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high-and low-growth chicken lines. Physiol Genomics. 2009;38:241–9.

    Article  CAS  PubMed  Google Scholar 

  41. Demeure O, Duclos MJ, Bacciu N, Le Mignon G, Filangi O, Pitel F, et al. Genome-wide interval mapping using SNPs identifies new QTL for growth, body composition and several physiological variables in an F2 intercross between fat and lean chicken lines. Genet Sel Evol. 2013;45:36.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Navarro P, Visscher PM, Knott SA, Burt DW, Hocking PM, Haley CS. Mapping of quantitative trait loci affecting organ weights and blood variables in a broiler layer cross. Br Poult Sci. 2005;46:430–42.

    Article  CAS  PubMed  Google Scholar 

  43. Wolfenson D, Frei YF, Snapir N, Berman A. Heat stress effects on capillary blood flow and its redistribution in the laying hen. Pflügers Arch. 1981;390:86–93.

    Article  CAS  PubMed  Google Scholar 

  44. Sandercock DA, Hunter RR, Nute GR, Mitchell MA, Hocking PM. Acute heat stress-induced alterations in blood acid-base status and skeletal muscle membrane integrity in broiler chickens at two ages: implications for meat quality. Poult Sci. 2001;80:418–25.

    Article  CAS  PubMed  Google Scholar 

  45. Buitenhuis AJ, Rodenburg TB, Van Hierden YM, Siwek M, Cornelissen SJ, Nieuwland MG, et al. Mapping quantitative trait loci affecting feather pecking behavior and stress response in laying hens. Poult Sci. 2003;82:1215–22.

    Article  CAS  PubMed  Google Scholar 

  46. Lu Y, Chen SR, Liu WB, Hou ZC, Xu GY, Yang N. Polymorphisms in Wnt signaling pathway genes are significantly associated with chicken carcass traits. Poult Sci. 2012;91:1299–307.

    Article  CAS  PubMed  Google Scholar 

  47. Ankra-Badu GA, Le Bihan-Duval E, Mignon-Grasteau S, Pitel F, Beaumont C, Duclos M, et al. Mapping QTL for growth and shank traits in chickens divergently selected for high or low body weight. Anim Genet. 2010;41:400–5.

    CAS  PubMed  Google Scholar 

  48. Ambo M, Moura AS, Ledur MC, Pinto FL, Baron EE, Ruy DC, et al. Quantitative trait loci for performance traits in a broiler × layer cross. Anim Genet. 2009;40:200–8.

    Article  CAS  PubMed  Google Scholar 

  49. Podisi BK, Knott SA, Burt DW, Hocking PM. Comparative analysis of quantitative trait loci for body weight, growth rate and growth curve parameters from 3 to 72 weeks of age in female chickens of a broiler–layer cross. BMC Genet. 2013;14:22.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Cooper MA, Washburn KW. The relationships of body temperature to weight gain, feed consumption, and feed utilization in broilers under heat stress. Poult Sci. 1998;77:237–42.

    Article  CAS  PubMed  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. Sewalem A, Morrice DM, Law A, Windsor D, Haley CS, Ikeobi CO, et al. Mapping of quantitative trait loci for body weight at 3, 6, and 9 weeks of age in a broiler layer cross. Poult Sci. 2002;81:1775–81.

    Article  CAS  PubMed  Google Scholar 

  53. Nassar MK, Goraga ZS, Brockmann GA. Quantitative trait loci segregating in crosses between New Hampshire and White Leghorn chicken lines: II. Muscle weight and carcass composition. Anim Genet. 2012;43:739–45.

    Article  CAS  PubMed  Google Scholar 

  54. McElroy JP, Kim JJ, Harry DE, Dekkers JC, Lamont SJ. Identification of trait loci affecting white meat percentage and other growth and carcass traits in commercial broiler chickens. Poult Sci. 2006;85:593–605.

    Article  CAS  PubMed  Google Scholar 

  55. Le Rouzic A, Álvarez-Castro JM, Carlborg Ö. Dissection of the genetic architecture of body weight in chicken reveals the impact of epistasis on domestication traits. Genetics. 2008;179:1591–9.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Xie L, Luo C, Zhang C, Zhang R, Tang J, Nie Q, et al. Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits. PLoS One. 2012;7:e30910.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Uemoto Y, Sato S, Odawara S, Nokata H, Oyamada Y, Taguchi Y, et al. Genetic mapping of quantitative trait loci affecting growth and carcass traits in F2 intercross chickens. Poult Sci. 2009;88:477–82.

    Article  CAS  PubMed  Google Scholar 

  58. Gu X, Feng C, Ma L, Song C, Wang Y, Da Y, et al. Genome-wide association study of body weight in chicken F2 resource population. PLoS One. 2011;6:e21872.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Donohue M, Cunningham DL. Effects of grain and oilseed prices on the costs of US poultry production. J Appl Poult Res. 2009;18:325–37.

    Article  Google Scholar 

  60. Yuan J, Wang K, Yi G, Ma M, Dou T, Sun C, et al. Genome-wide association studies for feed intake and efficiency in two laying periods of chickens. Genet Sel Evol. 2015;47:82.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  61. Mignon-Grasteau S, Rideau N, Gabriel I, Chantry-Darmon C, Boscher MY, Sellier N, et al. Detection of QTL controlling feed efficiency and excretion in chickens fed a wheat-based diet. Genet Sel Evol. 2015;47:74.

    Article  PubMed Central  PubMed  Google Scholar 

  62. Fuquay JW. Heat stress as it affects animal production. J Anim Sci. 1981;52:164–74.

    CAS  PubMed  Google Scholar 

  63. Matsuki S, Iuchi Y, Ikeda Y, Sasagawa I, Tomita Y, Fujii J. Suppression of cytochrome c release and apoptosis in testes with heat stress by minocycline. Biochem Biophys Res Commun. 2003;312:843–9.

    Article  CAS  PubMed  Google Scholar 

  64. Coble DJ, Fleming D, Persia ME, Ashwell CM, Rothschild MF, Schmidt CJ, et al. RNA-seq analysis of broiler liver transcriptome reveals novel responses to high ambient temperature. BMC Genomics. 2014;15:1084.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Morimoto RI. Regulation of the heat shock transcriptional response: cross talk between a family of heat shock factors, molecular chaperones, and negative regulators. Genes Dev. 1998;12:3788–96.

    Article  CAS  PubMed  Google Scholar 

  66. Waters ER, Lee GJ, Vierling E. Evolution, structure and function of the small heat shock proteins in plants. J Experimental Botany. 1996;47:325–38.

    Article  CAS  Google Scholar 

  67. Malik S, Roeder RG. Dynamic regulation of pol II transcription by the mammalian Mediator complex. Trends Biochem Sci. 2005;30:256–63.

    Article  CAS  PubMed  Google Scholar 

  68. Park JM, Werner J, Kim JM, Lis JT, Kim Y-J. Mediator, not holoenzyme, is directly recruited to the heat shock promoter by HSF upon heat shock. Mol Cell. 2001;8:9–19.

    Article  CAS  PubMed  Google Scholar 

  69. Kühl NM, Rensing L. Heat shock effects on cell cycle progression. Cell Mol Life Sci. 2000;57:450–63.

    Article  PubMed  Google Scholar 

  70. Ait-Boulahsen A, Garlich JD, Edens FW. Effect of fasting and acute heat stress on body temperature, blood acid-base and electrolyte status in chickens. Comp Biochem Physiol A Comp Physiol. 1989;94:683–7.

    Article  CAS  PubMed  Google Scholar 

  71. Mashaly MM, Hendricks GL 3rd, Kalama MA, Gehad AE, Abbas AO, Patterson PH. Effect of heat stress on production parameters and immune responses of commercial laying hens. Poult Sci. 2004;83:889–94.

    Article  CAS  PubMed  Google Scholar 

  72. Calandra T, Roger T. Macrophage migration inhibitory factor: a regulator of innate immunity. Nat Rev Immunol. 2003;3:791–800.

    Article  CAS  PubMed  Google Scholar 

  73. Cowan KJ, Storey KB. Mitogen-activated protein kinases: new signaling pathways functioning in cellular responses to environmental stress. J Exp Biol. 2003;206:1107–15.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

AVG participated in sample collection, carried out all statistical analyses and drafted the manuscript. KJB participated in sample collection and contributed digestibility data for analysis. CMA and MEP participated in the design of the study, sample collection and acquisition of funding. MFR and CJS participated in the design of the study and acquisition of funding. SJL participated in the design of the study, participated in sample collection, participated in acquisition of funding and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors thank the Iowa State University Poultry Research Center staff for animal care; the Lamont, Persia, Rothschild and Ashwell lab personnel for assistance; Jack Dekkers, Dorian Garrick, Rohan Fernando and their research groups for GWAS advice; Melanie Hess and Andrew Hess for their statistical and coding advice. This research was supported by USDA-NIFA-AFRI Climate Change Award #2011-67003-30228; the USDA National Institute of Food and Agriculture, Hatch project #5358; and AVG was supported by a USDA National Needs Fellowship, Award #2011-38420-20050.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

  1. Department of Animal Science, Iowa State University, Ames, IA, USA

    Angelica Van Goor, Max F. Rothschild & Susan J. Lamont

  2. Department of Animal Science, University of California, Davis, CA, USA

    Kevin J. Bolek

  3. Department Poultry Science, North Carolina State University, Raleigh, NC, USA

    Chris M. Ashwell

  4. Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA

    Mike E. Persia

  5. Department of Animal and Food Sciences, University of Delaware, Newark, DE, USA

    Carl J. Schmidt

Authors
  1. Angelica Van Goor
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Kevin J. Bolek
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chris M. Ashwell
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Mike E. Persia
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Max F. Rothschild
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Carl J. Schmidt
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Susan J. Lamont
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Susan J. Lamont.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Van Goor, A., Bolek, K.J., Ashwell, C.M. et al. Identification of quantitative trait loci for body temperature, body weight, breast yield, and digestibility in an advanced intercross line of chickens under heat stress. Genet Sel Evol 47, 96 (2015). https://doi.org/10.1186/s12711-015-0176-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-015-0176-7

Keywords

Genetics Selection Evolution

ISSN: 1297-9686