- Research article
Quantitative trait loci identified for blood chemistry components of an advanced intercross line of chickens under heat stress
BMC Genomicsvolume 17, Article number: 287 (2016)
Heat stress in poultry results in considerable economic losses and is a concern for both animal health and welfare. Physiological changes occur during periods of heat stress, including changes in blood chemistry components. A highly advanced intercross line, created from a broiler (heat susceptible) by Fayoumi (heat resistant) cross, was exposed to daily heat cycles for seven days starting at 22 days of age. Blood components measured pre-heat treatment and on the seventh day of heat treatment included pH, pCO2, pO2, base excess, HCO3, TCO2, K, Na, ionized Ca, hematocrit, hemoglobin, sO2, and glucose. A genome-wide association study (GWAS) for these traits and their calculated changes was conducted to identify quantitative trait loci (QTL) using a 600 K SNP panel.
There were significant increases in pH, base excess, HCO3, TCO2, ionized Ca, hematocrit, hemoglobin, and sO2, and significant decreases in pCO2 and glucose after 7 days of heat treatment. Heritabilities ranged from 0.01-0.21 for pre-heat measurements, 0.01-0.23 for measurements taken during heat, and 0.00-0.10 for the calculated change due to heat treatment. All blood components were highly correlated within measurement days, but not correlated between measurement days. The GWAS revealed 61 QTL for all traits, located on GGA (Gallus gallus chromosome) 1, 3, 6, 9, 10, 12–14, 17, 18, 21–28, and Z. A functional analysis of the genes in these QTL regions identified the Angiopoietin pathway as significant. The QTL that co-localized for three or more traits were on GGA10, 22, 26, 28, and Z and revealed candidate genes for birds’ response to heat stress.
The results of this study contribute to our knowledge of levels and heritabilities of several blood components of chickens under thermoneutral and heat stress conditions. Most components responded to heat treatment. Mapped QTL may serve as markers for genomic selection to enhance heat tolerance in poultry. The Angiopoietin pathway is likely involved in the response to heat stress in chickens. Several candidate genes were identified, giving additional insight into potential mechanisms of physiologic response to high ambient temperatures.
Climate change has increased the frequency of severe heat waves and the global temperature is projected to become increasingly warmer . Heat stress in poultry negatively impacts animal production and welfare resulting in economic losses estimated to be between $125-165 million for the U.S. broiler poultry industry . During a severe heat wave in Iowa, over 1.5 million layer hens died .
To reduce core body temperature during periods of heat stress, blood flow to internal organs decreases and blood flow to the combs and other surface tissues increases in chickens . During periods of heat stress, blood volume and oxygen carrying capacity are altered  and dehydration, caused by increased respiration, can increase hematocrit . Energy availability, as determined by plasma glucose level, is increased in chickens exposed to heat stress .
During high ambient temperatures, chickens reduce feed intake by as much as 17 %, which reduces growth . However, metabolic and endocrine changes during heat stress also contribute to reduction in growth in broilers, as demonstrated by a pair-feeding study .
A major change in blood components is caused by heat-induced increased respiration, which results in respiratory alkalosis, a disturbance in the acid base balance characterized by an increase in blood pH accompanied by a decrease in pCO2. Respiratory alkalosis occurs in broilers during heat stress and is associated with reduced growth rate . Metabolic alkalosis is an additional measure of disturbances in acid base balance and is defined by a decrease in the fixed acid concentrations and an increase in fixed base concentrations within the extracellular fluid .
Electrolyte balance is essential for acid base balance, maintenance of cellular homeostasis, synthesis of tissue protein, electrical potential of cell membranes, enzymatic reactions, and maintaining osmotic pressure . Altering electrolyte amounts in feed partially ameliorates the negative impacts of heat stress in broiler chickens .
The goal of the current study was to identify the physiological changes and genomic regions associated with response to heat stress in chickens as characterized by the blood chemistry components, including pH, pCO2, pO2, base excess (BE), HCO3, TCO2, K, Na, ionized Ca (iCa), hematocrit (Hct), hemoglobin (Hb), sO2, and glucose (Glu). In a commercial egg laying population, developmental measures have been established with hopes of using measures of blood chemistry components for selection . To date, few studies have identified quantitative trait loci (QTL) for blood components in chicken [15–18]. We used a 600 K SNP panel to identify QTL regions associated with levels of blood components of chickens under thermoneutral and heat stress conditions, and changes induced by heat.
Blood component measurements and heritabilities
Phenotypic means and heritabilities are given in Table 1 for blood components measured pre-heat (day 20 of age), after 7 days of heat treatment (day 28 of age), and the calculated change due to heat treatment (day 28–20). After 7 days of heat treatment, pH, BE, HCO3, TCO2, iCa, Hct, Hb, and sO2 significantly increased while pCO2 and glucose significantly decreased. There were no significant changes in pO2, K, and Na due to heat treatment.
Heritabilities ranged from 0.01-0.21 for pre-heat measurements, 0.01-0.23 for measurements taken during heat, and 0.00-0.10 for the calculated change due to heat treatment.
Correlations between blood components at each measurement phase are given in Fig. 1 as a heat map. Almost all blood components were positively correlated with all other variables measured on the same day. Very few significant correlations, however, occurred between variables measured on different days.
Of the 480 genotyped birds, 458 Advanced Intercross Line (AIL) and all 12 parental line birds passed the whole animal DishQC criterion. Of the 580,961 SNPs on the array, filtering based on SNP call rate ≥ 95 % removed a small proportion (59,789 SNPs), whereas filtering based on MAF removed a much larger proportion (311,055 SNPs), yielding 210,117 SNPs for subsequent analyses.
The results from the GWAS for each trait are depicted in Fig. 2. A wide range of genetic variation (0.5-9.8 %) was explained by each significant window and detailed information is found in Table 2. Adjacent windows that were significant for a single trait are discussed below as a single QTL region.
Six QTL for pH phenotypes were identified: three for pH20 with two on GGA18 and one on GGA28, one for pH28 on GGA12, and two for pH28-20 with one each on GGA6 and GGA10.
Nine QTL for pCO2 measurements were identified: one for pCO220 on GGA28, four for pCO228 located on GGA1, 9, 10, and 27, and four for pCO228-20 on GGA3, 10, 23, and 28. No QTL were identified for pO220 or for pO228-20. One QTL was identified for pO228 on GGA13.
A total of five QTL were identified for BE traits: two for BE20 on GGA18, three for BE28 with one each on GGA1, 21, and 27, and none for BE28-20. Nine QTL were identified for TCO2 traits: none for TCO220; eight for TCO228 one each on GGA6 and GGA26, and six on GGAZ, and one for TCO228-20 on GGA10. No QTL were identified for HCO320 or HCO328-20, while seven were revealed for HCO328 with one on GGA6 and six on GGAZ.
Five QTL for K traits were identified: four for K20 with two on GGA10, one on GGA12, and one on GGA26, none for K28 and one for K28-20 located on GGA12. No QTL were identified for the Na phenotypes. A single QTL was identified for ionized Ca phenotypes: for iCa28 on GGA26.
We identified five QTL for Hct measurements: none for Hct20 or Hct28-20, and five for Hct28 located one each on GGA1, 10, 14, 22, and two on GGA28. Seven QTL were identified for Hb: none for Hb20, six for Hb28 located one each on GGA1, 10, 14, 22, and two on 28, and one for Hb28-20 on GGA22. There were three QTL for sO2 phenotypes: none for sO220, two for sO228 located on GGA24 and GGA25, and one for sO228-20 on GGA17.
Four QTL were identified for Glu: one for Glu20 on GGA10, and three for Glu28 with one on GGA22 and two on GGAZ.
The pathway analysis of all annotated genes within significant QTL regions across all measured traits, and separately for genes in the regions of QTL co-localization, and the top 20 significant (P ≤ 0.05) canonical pathways for each group are listed in Table 3. Of the 999 genes identified within all significant QTL regions, 682 genes were annotated within IPA and used for the pathway analysis. Two canonical pathways of interest for all identified QTL include the AMPK signalling and Angiopoietin signalling pathways. Of the 226 genes in regions of QTL co-localization, 185 were annotated within IPA and used for pathway analysis. A pathway of particular interest that was revealed was the Cardiac Hypertrophy signalling pathway.
We explored regions of QTL co-localization in detail to identify candidate genes that may give insight into the complex biological mechanisms that control blood component response to heat stress. Candidate genes were identified using Ensemble Biomart within the 1 Mb windows that were significant for 3 or more traits (Additional file 1: Table S1).
The aim of this study was to identify and estimate the effect of QTL, and to perform a functional analysis using positional candidate genes, for blood components (pH, pCO2, pO2, base excess, HCO3, TCO2, K, Na, ionized Ca, Hct, Hb, sO2, and Glu) using a novel AIL of chickens under heat stress and a 600 K SNP panel for genotyping. The blood components measured were within the accepted range reported for chicken . Blood chemistry components are grouped into functional categories (i.e., respiratory alkalosis, metabolic alkalosis, blood volume and oxygen carrying capacity, electrolytes, and glucose) for discussion.
Previous generations of this AIL were used for several QTL mapping studies and allowed the identification of many QTL including 257 for growth and body composition [20–24], 93 for skeletal integrity , 51 for metabolic traits , 12 for response to Salmonella enteritidis challenge [26–28], and 35 for response to heat stress . Therefore, collectively, a wide range of traits have been associated with a large number of QTL in this AIL. The continued erosion of Linkage Disequilibrium (LD) in this population over subsequent generations, combined with the availability of larger SNP panels, creates a unique opportunity to more finely map the location of QTL that are in LD with a causal mutation.
During periods of intense heat, chickens increase the depth and frequency of respiration to decrease core body temperature . Broilers that are heat stressed increase panting and display signs of respiratory alkalosis , which is caused by an increase in the amount of CO2 expelled from the lungs, and a consequent increase in pH within the blood, and an increase in pO2 within the blood. We investigated blood pH, pCO2, and pO2 to characterize respiratory alkalosis induced by heat stress.
Occurrence of respiratory alkalosis was clearly demonstrated in the current study by a significant increase in blood pH and significant decrease in pCO2 due to heat treatment, in agreement with previous studies. Heat stress for two hours at 32 °C in broilers at 35 days of age significantly increases blood pH and decreases pCO2  and, in another study using broilers, heat stressed at 32 °C for 2 weeks at 28 days of age in birds that were panting . We found pO2 increased in response to heat treatment, although not significantly. In a study using 35 day old broilers, blood pO2 significantly increased after cyclical heat stress for 10 days at 35 °C .
Only one other published study has estimated heritabilities of blood components in chickens under thermal stress . The current study, therefore, adds substantially to the body of information on response of birds to thermal stress by estimating heritabilities of blood component levels and changes under heat stress and thermoneutral conditions. In broiler chickens at 22 days of age reared under cold stress conditions, heritabilities for blood pH, pCO2, and pO2 were estimated at 0.15, 0.15, and 0.03, respectively , in agreement with the current study’s estimates for thermoneutral and heat conditions. Our estimates for the changes in these blood components due to heat treatment was much lower, suggesting that the ability to select for the response to heat stress may be difficult.
To our knowledge, QTL for blood pH, pCO2, and pO2 in chickens have not been previously reported. Identification of QTL for blood pH on different chromosomes across measurement phases, indicates that genetic control of these traits exists and is partly dependent on the environment. Co-localized QTL for pCO220 and pCO228-20 on GGA28, and for pCO228 and pCO228-20 on GGA10, suggest that the same genetic regions contribute to control of pCO2 level independent of environmental temperature. The presence of co-localized QTL between measurement phases was not expected, based on the lack of phenotypic correlations (r = 0.00).
Metabolic alkalosis occurs when there is a disturbance in the fixed acids and bases in the extracellular fluid . Imbalance of dietary Na, K, or Ca can result in metabolic alkalosis , which is characterized by an increase in blood pH, HCO3, and base excess, and can be induced in growing layers by high levels of calcium in feed .
Base excess is considered a comprehensive measure of the metabolic components of bases, which reflects the nonrespiratory contribution to changes in acid–base disturbances . Base excess can be altered by changing the cation:anion ratio in the diet of broiler chickens and is associated with body weight and bone density . In the current study, base excess significantly increased after heat treatment, which is consistent with the hypothesis that chickens experience metabolic alkalosis under heat stress.
HCO3 is the most abundant buffer in the blood, is primarily regulated by the kidneys, and is a metabolic component of acid–base balance . We observed a significant increase in HCO3 due to heat treatment. These results contrasted with a previous study using broilers at 28 days of age in which blood HCO3 significantly decreased in panting birds under acute heat stress , and another study using male broilers that reported a decrease in HCO3 after a heat stress at 32 °C for 10 h . TCO2 also increased in response to heat treatment. It was unexpected to observe a decrease in base excess, consistent with metabolic alkalosis, while HCO3 and TCO2 increased, because the traits are highly positively correlated within all treatment phases (r ≥ 0.95).
We estimated heritability of base excess between 0.00-0.10, of HCO3 between 0.03-0.23, and of TCO2 between 0.01-0.13. In broiler chickens at 22 days of age reared under cold stress conditions, blood HCO3 and TCO2 heritability were both estimated at 0.19 .
We are the first to report QTL in chickens for blood base excess, HCO3, and TCO2, which are related to metabolic alkalosis. QTL for base excess are located on separate chromosomes for all measurement phases, indicating a strong genetics by environmental (G x E) temperature interaction. The phenotypic correlations for base excess between measurement phases were both very low (r = 0.03). The QTL for base excess on GGA18 overlap with pH measured at thermoneutrality and were highly correlated (r = 0.78). Surprisingly, QTL for HCO3 were only identified during heat treatment and were on GGA6 and GGAZ. Ten of the eleven QTL for TCO2 measured during heat co-localized with QTL for HCO3 and these co-localized regions were located on GGA6, 26, and Z.
Blood K and Na levels numerically increased and iCa statistically increased in response to heat treatment. This is in disagreement with previous reports of decreasing levels of both K and Na in response to heat stress, likely due to increased water intake which results in decreased concentrations of electrolytes within the blood [6, 13, 38].
Heritability of K and Na blood levels in humans has been estimated to be very low, 0.03 and 0.04, respectively , in agreement with our low heritability estimates during heat and for the calculated differential. In contrast, our estimates for heritability under thermoneutral conditions for K and Na were higher, 0.20 and 0.08, respectively. Estimated heritability was 0.02 for ionized Ca measured during heat stress, lower than the 0.19 of mice in thermoneutral conditions . The estimated heritability was low, for both thermoneutral (0.04) and the differential due to heat (0.01), indicating the genetic component for ionized Ca is dependent upon environmental conditions at the time of measurement. The low heritabilities of these traits during heat and for the calculated differential due to heat treatment suggest it may be difficult to select for these traits.
This research is the first to describe QTL for the electrolyte-balance traits of blood K, Na, and ionized Ca in the chicken. In swine, QTL have been identified for these traits . QTL for blood K were located on GGA10, 12, and 26. QTL were identified for K across the thermoneutral and differential due to heat measurement phases, indicating genetic control of this component in this region on GGA12 despite environmental temperature. The correlation between thermoneutral and the differential was moderate (r = 0.10). No significant QTL for Na were identified in the current study and a single QTL for ionized Ca was located on GGA26 for the measurement taken during heat.
Blood volume and oxygen saturation
Changes in blood volume and oxygen carrying capacity occur in chickens during periods of heat stress . Both hematocrit and hemoglobin significantly increased due to heat treatment, which may be the result of dehydration. This result contrasts with a previous study using male broilers in which both decreased after an acute heat stress at 32 °C for 10 h . Blood sO2 is a measure of oxyhemoglobin in relation to total hemoglobin that is able to bind oxygen , and this significantly increased during heat treatment.
The heritability of Hct was estimated as very low at 0.01 and 0.02 for pre-heat and the differential, respectively, while during heat was moderately heritable at 0.21. Heritability has been estimated for hematocrit at 0.39 in domestic fowl . The increase in heritability when measured during heat stress indicates that this trait may be useful for selection. Heritability estimates of sO2 were very low (0.01-0.03), which is in general agreement with a previously reported value of 0.07 in cold-stressed broiler chickens at 22 days of age .
Seven QTL for haematocrit have been identified in chickens (www.animalgenome.org). In a broiler by layer F2 intercross, QTL for hematocrit were located on GGA1, 2, 6, and 14 ; in a Fayoumi by Leghorn F2 intercross on GGA1 and GGA15 , and in a broiler by layer cross on GGA1 . Our current work confirmed previously identified QTL for Hct28 on GGA1 and GGA14. Novel QTL for Hct were on GGA10, 22, and 28. Most of the QTL identified in the current study for Hb co-localized with those identified for Hct, with the addition of a relatively large QTL for Hb28-20 on GGA22, explaining 1.7 % of the genetic variation. The co-localization of QTL among Hct and Hb is expected because they have very high positive phenotypic correlations across all measurement phases (r ≥ 0.99). We identified novel QTL for sO2 on GGA17, 24, and 25, none of which overlapped between measurement phases, indicating separate genetic control of this trait dependent upon environmental temperature. A previous study using a commercial broiler line identified one on GGA16 . Thus, QTL for sO2 appear to be population specific.
Glucose is the body’s primary source of energy, and blood Glu significantly decreased due to heat treatment in the current study. In contrast, male broilers had a significant increase in Glu after heat stress at 32 °C for 10 h , and in broiler chicks of 5 weeks of age at 35-40 °C . In chicken lines divergently selected for blood glucose concentration, the low glucose line was less efficient at food utilization compared to the high glucose line , which may indicate that the decrease in glucose we see during heat stress may contribute to inefficiency in food utilization.
The current study estimated heritabilities for glucose ranging between 0.02-0.19. In a study using chickens divergently selected for blood glucose concentration, heritability was estimated at 0.25 .
We identified QTL for Glu20 and Glu28 on GGA10, 22, and Z, while QTL were mapped to GGA2, 7, and Z in the F2 generation of the same chicken population under thermoneutral conditions . The two studies may have detected the same QTL on chromosome Z and, due to the breakdown of LD over the generations, the current study may have mapped the QTL more accurately. In an F2 intercross between fat and lean broilers, QTL were identified for blood glucose on GGA3 and GGA18 , and for fasting plasma glucose on GGA5, 6, 13, and 26 . A study using an F2 of broilers divergently selected for growth, identified QTL for plasma glucose on GGA20 and GGA27 . Thus, QTL location for blood glucose level appears to be heat and/or population specific.
Considering all measured traits, we identified a total of 32 unique QTL. All annotated genes within the QTL regions were used for pathway analysis using IPA and many significantly associated canonical pathways were identified including AMPK signalling and Angiopoietin signalling were identified. The top 20 pathways are found in Table 3. AMPK is a master metabolic regulator involved in metabolism  and, thus, may be a pathway which warrants further investigation for involvement in production traits during heat stress. During high ambient temperatures chickens redirect blood flow to the body surface to decrease body temperature , and the angiopoietin signalling pathway functions in blood vessel development which may help alleviate temperature stress.
The co-localized regions resulted in many significant canonical pathways and the top 20 pathways are found Table 3. Of particular interest is the Cardiac Hypertrophy signalling pathway (P = 4.35E-02). QTL for hemoglobin and hematocrit represent 3 (7 total) regions of co-localization and there is a positive linear relationship between hematocrit and heart weight in chickens under heat stress ; therefore, this pathway likely contributes to the response to heat stress in chickens.
Candidate genes for co-localized QTL
The QTL regions that co-localized for three or more traits were further investigated for positional, functional candidate genes to give further insight into the biological mechanisms involved in the response of blood components to heat stress. The identified genes are located in Additional file 1: Table S1.
There are 51 genes in the region on GGA10 between 3–6 Mb that contained QTL for Glu20, pCO228, and TCO228-20. With 2 of these 3 traits associated with CO2 concentration, CA12 (carbonic anhydrase) is a likely candidate gene involved in the CO2 response to heat stress. Carbonic anhydrases catalyse the reaction of CO2 and H2O to form HCO3 and H+, and thus may stabilize blood acid base balance during heat stress. Another strong functional candidate in this region is HSP40, a member of the heat shock protein family that functions as a molecular chaperone to prevent cellular damage during heat stress . A candidate gene in this region for glucose level is GCNT3, a glucosamine acetyl transferase which is associated with glucose metabolism in humans .
Fourteen genes were identified on GGA10 between 16–17 Mb, where QTLs co-localized for pH28-20, Hct28, Hb28, and K20. Many QTL in chicken have been identified in this region including those related to growth [22, 53–55], abdominal fat [23, 49, 56], and the stress-associated trait of fear response . A strong candidate gene is ALDH6 (aldehyde dehydrogenase) which functions to convert aldehydes to carboxylic acids. This gene may function to maintain blood acid base balance during heat stress. Another gene in this region is IGF1 (insulin like growth factor 1), which has many roles and is a biomarker for growth .
Four genes were identified on GGA22 between 3–4 Mbs, where QTL were co-localized for Hct28, Hb28, Hb28-20, and Glu28. To our knowledge, no QTL have been reported in this region. Because all traits were measured during heat treatment or as the differential, we propose these to be heat specific QTL. Candidate genes TGFA (pretransforming growth factor) and ADRA1A (adrenergic receptor) both regulate cell growth. It is known that metabolic changes occur during periods of heat in chickens that contribute to reduction in growth, independent upon feed intake .
There are 48 genes in the 1 Mb region on GGA26 between 3–4 Mbs, where QTL co-localized for TCO228, K20, and iCa28. Notably, a QTL for tibia bone mineral density identified in a commercial broiler and layer cross is located within this region . This co-localization suggests that this locus might be involved in both blood calcium and bone density, and therefore, may be an ideal candidate for further investigation to understand the physiological response to heat stress on bone mineral density.
There are 86 genes in the 2 Mb region on GGA28 between 3–5 Mb where QTLs co-localize for pH20, Hb28, Hct28, pCO220, and pCO228-20. A QTL for heart weight, relating to susceptibility of pulmonary hypertension  co-localizes with those identified here. Many of these genes are related to membrane transport of solutes and DNA transcription. The solute carriers SLC39A3, SLC25A42 and SLC35E1 were identified, as well as CHERP and CIB3, involved in calcium homeostasis. Transcription-related genes include SUGP1, which is involved in RNA splicing; RFXANK, a DNA-binding protein; NR2C2AP, a nuclear receptor protein; DDX49, an RNA helicase; ELL, an RNA polymerase II elongation factor; and SIN3B a transcriptional regulator.
On GGAZ, 2 genes were identified between 5–7 Mbs, where QTL co-localize for Glu28, HCO328, and TCO228. The only reported QTL near this region is for antibody response to KLH antigen . Heat stress is known to reduce antibody titre in chickens , and this locus may be involved in the complex interaction of heat and antibody titre. Although, antibody levels were not measured in the current study. During periods of heat stress, DNA transcription, RNA translation, and cellular proliferation are altered  and we observed several genes in this region related to these particular responses including: KIAA1328, involved in chromosomal integrity during mitosis; and TPGS2, involved in tubulin formation.
On GGAZ, 21 genes were identified between 69–71 Mbs, where QTL co-localize for Glu28, HCO328, and TCO228. The one QTL that is near this region was identified in a previous generation of the same AIL as the current study, and is for bone mineral density . A recent study found that heat stress in broilers results in decreased bone mineral density . In humans, low serum bicarbonate levels are associated with decreased bone mineral density . Although this relationship has yet to be elucidated in the chicken, further studies should investigate the association between blood chemistry variables and bone mineral density. The genes identified in the current study that are primarily involved in DNA transcription include XPA, which is a DNA repair protein, FOXE3 which is part of the forkhead box, and SNORA66 which is small nuclear RNA. Additionally, microRNAs gga-mir-2954, gga-mir-2131, and gga-mir-1583 were identified in this region. An additional gene of interest identified was DNAJA1, which is part of the heat shock family of proteins.
QTL for blood components reveal orthologous genes between chicken and swine
QTL for blood pCO2 in the current study were located on GGA1, 3, 9, 10, 23, 27, and 28. In swine, QTL for blood pCO2 are on chromosomes 6, 7, 8, 9, and X . We identified a region of synteny between chicken GGA1, 110–111 Mb, and pig chromosome X, 43–44 Mb (Fig. 3a), which contains a pCO2 QTL and several orthologous genes including FUNDC1, EFHC2, NDP, and MAOA. Another region of synteny exists between chicken chromosome 10, 1–4 Mb, and pig chromosome 7, 53–65 Mb (Fig. 3b/c), which contains several orthologous genes including, but not limited to, UBE2Q2, DNAJ, GRAMD2, ADPGK, NEO1, CLK3, SCAMP5, CSK, and MPI. This region contains the carboxylic anhydrase gene (CA12) in chicken, which is involved in calcium metabolism, but this gene maps on pig chromosome 1, a chromosome on which no QTL have been reported for blood chemistry measurements. The region on GGA10, 1–4 Mb, contains QTL for Glu20, pCO228, pCO228-20, and TCO228-20. The syntenic region in swine contains co-localized QTL for pCO2, HCO3, TCO2, and base excess .
A QTL for blood K level mapped to syntenic regions in chicken GGA10, 16–17 Mb, in our line and swine chromosome 1, 63–226 Mb (Fig. 3d) in a previous study . An orthologous gene of interest in this region is IGF-1.
The results of this study contribute to the currently sparse knowledge of levels and heritabilities of several blood components under thermoneutral and heat stress conditions in chickens. Most blood components changed in response to heat treatment. Mapped QTL may serve as markers for genomic selection to enhance heat tolerance in poultry and several candidate genes were identified which may give additional insight into mechanisms of physiologic response to high ambient temperatures.
Animal experiments were approved by the Institutional Animal Care and Use Committee of Iowa State University: Log #4-11-7128-G.
We used the F18 and F19 generations of an AIL between chicken lines divergent for thermotolerence created by crossing a single broiler sire to six highly inbred Fayoumi dams . Birds were reared in floor pens with wood shavings bedding and had ad libitum access to water and feed that met all NRC requirements .
Heat stress experimental design
A total of 631 birds from four hatches (two hatches in each of the two generations) were used for independent heat stress experiments (four replicates). At 17 days of age, birds were transferred to environmentally controlled chambers and acclimated for five days. Multiple chambers, each containing 6 pens, were used per replicate. Ten to 12 birds were placed in each pen. From day 22 to 28 of age, the chambers heated to 35 °C for 7 h per day and remained at 25 °C at all other times.
Blood variable measurements
Blood was collected from the wing vein on day 20 (pre-heat) and day 28 (during heat) using a heparinized syringe and needle, and analysed immediately using an iSTAT Portable Clinical Analyser . The iSTAT CG8+ cartridge was utilized to measure thirteen blood variables including; pH, pCO2, pO2, base excess, HCO3, TCO2, K, Na, ionized Ca, hematocrit, hemoglobin, sO2, and glucose.
DNA isolation and genotyping
Blood was collected from the wing vein by using an EDTA-coated syringe and needle, and 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 while the supernatant remained. The supernatant was combined with 70 % ethanol to precipitate out DNA. The DNA isolated from 468 AIL, 6 broiler, and 6 Fayoumi chickens was genotyped on the Affymetrix 600 K chicken SNP axiom array  by GeneSeek Inc., Lincoln, NE. SNP chromosomal locations were based on the Gallus_gallus_4.0 assembly through Ensembl.
Calculations of means and standard errors, fixed effects and covariates for the GWAS were calculated based on ANOVA (analysis of variance), and significant terms were fit as fixed effects with a P value ≤ 0.05 using JMP statistical software . Heritabilities were estimated with an animal model using ASReml software .
Parameters for inclusion of SNP genotypes included SNP call rate ≥ 95 % and minor allele frequency ≥ 5 %. Genotyping console (Affymetrix) software was used to create genotyping calls and quality control based on whole animal DishQC score ≥ 0.7. The SNPolisher (Affymetrix) R package was used for quality control of individual SNP in all animals with passing DishQC scores.
The GWAS of phenotypic traits with SNP genotypes was done using GenSel software . Bayes B, which fits all SNPs simultaneously as random effects, was used for the analysis. The mixed model used for the GWAS:
Where y = vector form of phenotypes, X = incidence matrix to account for fixed effects on phenotypes, b = vector of fixed effects, z j = 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 = allele substitution effect for SNP j, δ j = whether SNP j was included in the Markov chain Monte Carlo (MCMC) chain, and ε is the error associated with the analysis.
The genomic markers were split into 1001 non-overlapping 1 Mb windows across the genome. A total of 41,000 MCMC iterations were run for each analysis and the first 1000 iterations were discarded (burn in). The δj was set so that π = 0.9978 to avoid fitting more SNPs than number of animals in a given iteration. In 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 ≥ 0.5 % of the total genetic variation, corresponding to 5 times more observed than expected.
To further investigate QTL regions, we conducted a pathway analysis using Ingenuity Pathway Analysis (IPA) software. All annotated genes within significant (explaining ≥ 0.05 % of the genetic variation) 1 Mb windows for any measured trait were identified using Ensemble biomart. This gene list was used as input into IPA and a core analysis was completed using default parameters to identify significant (P ≤ 0.05) canonical pathways and the top 20 significant pathways were reported. Additionally, a gene list was created using the regions of QTL co-localization (3 or more traits) and analysed as described for all QTL regions.
Candidate genes were identified for regions of QTL co-localization (3 or more traits). All genes within the region were identified using ENSEMBL biomart .
Syntenic regions between chicken and swine
To identify syntenic regions for reported QTL for the same blood chemistry component measurements between chicken and pig, the Comparative Genomics option was used in Ensembl .
Availability of Data and Materials
The dataset supporting the conclusions of this article is available in the Animal QTLdb (animalgenome.org) repository and can be found at http://www.animalgenome.org/cgi-bin/QTLdb/GG/pubtails?PUBMED_ID=ISU0082. The phenotypic dataset supporting the conclusions of this article is included within the article as an additional file (Additional file 2: Table S2).
advanced intercross line
genome wide association study
quantitative trait loci
USDA. Climate Change Science White Paper. 2012.
St-Pierre N, Cobanov B, Schnitkey G. Economic losses from heat stress by US livestock industries. J Dairy Sci. 2003;86:E52–77.
Administration NOaA. Summer of 1995 heat wave. Online. 1995.
Wolfenson D, Frei YF, Snapir N, Berman A. Heat stress effects on capillary blood flow and its redistribution in the laying hen. Pflugers Arch. 1981;390(1):86–93.
Yahav S, Straschnow A, Plavnik I, Hurwitz S. Blood system response of chickens to changes in environmental temperature. Poult Sci. 1997;76(4):627–33.
Borges S, Da Silva AF, Majorka A, Hooge D, Cummings K. Physiological responses of broiler chickens to heat stress and dietary electrolyte balance (sodium plus potassium minus chloride, milliequivalents per kilogram). Poult Sci. 2004;83(9):1551–8.
Garriga C, Hunter RR, Amat C, Planas JM, Mitchell MA, Moretó M. Heat stress increases apical glucose transport in the chicken jejunum. Am J Physiol Regul Integr Comp Physiol. 2006;290(1):R195–201.
Austic R. Feeding poultry in hot and cold climates. 1985.
Lin H, Malheiros R, Moraes V, Careghi C, Decuypere E, Buyse J. Acclimation of broiler chickens to chronic high environmental temperature. Archiv Fur Geflugelkunde. 2004;39:46.
Teeter R, Smith M, Owens F, Arp S, Sangiah S, Breazile J. Chronic heat stress and respiratory alkalosis: occurrence and treatment in broiler chicks. Poult Sci. 1985;64(6):1060–4.
Galla JH. Metabolic alkalosis. J Am Soc Neph. 2000;11(2):369–75.
Borges S, Fischer Da Silva A, Maiorka A. Acid–base balance in broilers. World Poult Sci J. 2007;63(01):73–81.
Borges S, Da Silva AF, Ariki J, Hooge D, Cummings K. Dietary electrolyte balance for broiler chickens exposed to thermoneutral or heat-stress environments. Poult Sci. 2003;82(3):428–35.
Schaal T, Arango J, Wolc A, Brady J, Fulton J, Rubinoff I, et al. Commercial Hy-Line W-36 pullet and laying hen venous blood gas and chemistry profiles utilizing the portable i-STAT® 1 analyzer. Poult Sci. 2015;95:pev350.
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 Gen. 2009;38(3):241–9.
Park H-B, Jacobsson L, Wahlberg P, Siegel PB, Andersson L. QTL analysis of body composition and metabolic traits in an intercross between chicken lines divergently selected for growth. Physiol Gen. 2006;25(2):216–23.
Jacobsson L, Park H-B, Wahlberg P, Fredriksson R, Perez-Enciso M, Siegel PB, et al. Many QTLs with minor additive effects are associated with a large difference in growth between two selection lines in chickens. Genet Res. 2005;86(02):115–25.
Zhou H, Evock-Clover C, McMurtry J, Ashwell C, Lamont S. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. IV. Metabolic traits. Poult Sci. 2007;86(2):267–76.
Martin MP, Wineland M, Barnes HJ. Selected blood chemistry and gas reference ranges for broiler breeders using the i-STAT® handheld clinical analyzer. Avian Dis. 2010;54(3):1016–20.
Deeb N, Lamont S. Use of a novel outbred by inbred F1 cross to detect genetic markers for growth. An Genet. 2003;34(3):205–12.
Li H, Deeb N, Zhou H, Mitchell A, Ashwell C, Lamont S. Chicken quantitative trait loci for growth and body composition associated with transforming growth factor-beta genes. Poult Sci. 2003;82(3):347–56.
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.
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. II. Body composition. Poult Sci. 2006;85(10):1712–21.
Abasht B, Lamont S. Genome-wide association analysis reveals cryptic alleles as an important factor in heterosis for fatness in chicken F2 population. An Genet. 2007;38(5):491–8.
Zhou H, Deeb N, Evock-Clover C, Mitchell A, Ashwell C, Lamont S. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. III. Skeletal integrity. Poult Sci. 2007;86(2):255–66.
Kaiser M, Lamont S. Microsatellites linked to Salmonella enterica Serovar Enteritidis burden in spleen and cecal content of young F1 broiler-cross chicks. Poult Sci. 2002;81(5):657–63.
Kaiser M, Deeb N, Lamont S. Microsatellite markers linked to Salmonella enterica serovar enteritidis vaccine response in young F1 broiler-cross chicks. Poult Sci. 2002;81(2):193–201.
Liu W, Lamont S. 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. An Biotech. 2003;14(1):61–76.
Van Goor A, Bolek KJ, Ashwell CM, Persia ME, Rothschild MF, Schmidt CJ, 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. 2015;47(1):1–13.
Randall W, Hiestand W. Panting and temperature regulation in the chicken. Am J Physiol. 1939;127(4):761–7.
Sandercock D, Hunter R, Nute G, Mitchell M, Hocking P. 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(4):418–25.
Deyhim F, Teeter R. Research note: sodium and potassium chloride drinking water supplementation effects on acid–base balance and plasma corticosterone in broilers reared in thermoneutral and heat-distressed environments. Poult Sci. 1991;70(12):2551–3.
Closter A, Van As P, Groenen M, Vereijken A, Van Arendonk J, Bovenhuis H. Genetic and phenotypic relationships between blood gas parameters and ascites-related traits in broilers. Poult Sci. 2009;88(3):483–90.
Daghir N. Nutrient requirements of poultry at high temperature. Poultry production in hot climate. 2008;133:60.
Guoa X, Huanga K, Chenga H, Luo J, Pana C. High dietary calcium causes metabolic alkalosis in growing layers. Poult Sci. 2008;87(7):1353–7.
Chin CD, Linder V, Sia SK. Commercialization of microfluidic point-of-care diagnostic devices. Lab Chip. 2012;12(12):2118–34.
Halley J, Nelson T, Kirby L, Johnson Z. Effect of altering dietary mineral balance on growth, leg abnormalities, and blood base excess in broiler chicks. Poult Sci. 1987;66(10):1684–92.
Ait-Boulahsen A, Garlich J, Edens F. Effect of fasting and acute heat stress on body temperature, blood acid–base and electrolyte status in chickens. Comp Biochem Physiol. 1989;94(4):683–7.
Meyer TE, Verwoert GC, Hwang S-J, Glazer NL, Smith AV, Van Rooij F, et al. Genome-wide association studies of serum magnesium, potassium, and sodium concentrations identify six loci influencing serum magnesium levels. PLoS Genet. 2010;6(8), e1001045.
Tordoff MG, Bachmanov AA, Reed DR. Forty mouse strain survey of voluntary calcium intake, blood calcium, and bone mineral content. Physiol Behav. 2007;91(5):632–43.
Reiner G, Fischer R, Köhler F, Berge T, Hepp S, Willems H. Heritabilities and quantitative trait loci for blood gases and blood pH in swine. An Genet. 2009;40(2):142–8.
Washburn K. Heritability of packed red blood cell volume in the domestic fowl. Poult Sci. 1967;46(4):1025–6.
Boschiero C, Jorge EC, Ninov K, Nones K, Do Rosário MF, Coutinho LL, et al. Association of IGF1 and KDM5A polymorphisms with performance, fatness and carcass traits in chickens. J Appl Genet. 2013;54(1):103–12.
Der Laan M-H P-v, Bed’Hom B, Coville J-L, 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(1):31.
Navarro P, Visscher P, Knott S, Burt D, Hocking P, Haley C. Mapping of quantitative trait loci affecting organ weights and blood variables in a broiler layer cross. Brit Poult Sci. 2005;46(4):430–42.
Ewald S, Ye X, Avendano S, McLeod S, Lamont S, Dekkers J. Associations of BF2 alleles with antibody titres and production traits in commercial pure line broiler chickens. An Genet. 2007;38(2):174–6.
Khan W, Khan A, Anjuman A, Rehman Z. Effects of induced heat stress on some biochemical values in broiler chicks. Int J Agric Biol. 2002;4:74–5.
Leclercq B, Simon J, Ricard F. Effects of selection for high and low plasma glucose concentration in chickens. Brit Poult Sci. 1987;28(4):557–65.
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.
Mihaylova MM, Shaw RJ. The AMPK signalling pathway coordinates cell growth, autophagy and metabolism. Nat Cell Biol. 2011;13(9):1016–23.
Glover JR, Lindquist S. Hsp104, Hsp70, and Hsp40: a novel chaperone system that rescues previously aggregated proteins. Cell. 1998;94(1):73–82.
Dostrovsky N, Towheed T, Hudson R, Anastassiades T. The effect of glucosamine on glucose metabolism in humans: a systematic review of the literature. Osteoarth Cart. 2011;19(4):375–80.
Nassar M, Goraga Z, Brockmann G. Quantitative trait loci segregating in crosses between New Hampshire and White Leghorn chicken lines: IV. Growth performance. An Genet. 2015.
Ambo M, Moura A, Ledur M, Pinto L, Baron E, Ruy D, et al. Quantitative trait loci for performance traits in a broiler × layer cross. An Genet. 2009;40(2):200–8.
Zhou Y, Liu Y, Jiang X, Du H, Li X, Zhu Q. Polymorphism of chicken myocyte-specific enhancer-binding factor 2A gene and its association with chicken carcass traits. Molec Biol Rep. 2010;37(1):587–94.
Campos R, Nones K, Ledur M, Moura A, Pinto L, Ambo M, et al. Quantitative trait loci associated with fatness in a broiler–layer cross. An Genet. 2009;40(5):729–36.
Buitenhuis A, Rodenburg T, Siwek M, Cornelissen S, Nieuwland M, Crooijmans R, et al. Identification of QTLs involved in open-field behavior in young and adult laying hens. Behav Genet. 2004;34(3):325–33.
Renehan AG, Zwahlen M, Minder CT, O’Dwyer S, Shalet SM, Egger M. Insulin-like growth factor (IGF)-I, IGF binding protein-3, and cancer risk: systematic review and meta-regression analysis. Lancet. 2004;363(9418):1346–53.
Schreiweis M, Hester P, Settar P, Moody D. Identification of quantitative trait loci associated with egg quality, egg production, and body weight in an F2 resource population of chickens1. An Genet. 2006;37(2):106–12.
Rabie T, Crooijmans R, Bovenhuis H, Vereijken A, Veenendaal T, Poel JD, et al. Genetic mapping of quantitative trait loci affecting susceptibility in chicken to develop pulmonary hypertension syndrome. An Genet. 2005;36(6):468–76.
Sun Y, Biscarini F, Bovenhuis H, Parmentier H, Poel VDJ. Genetic parameters and across‐line SNP associations differ for natural antibody isotypes IgM and IgG in laying hens. An Genet. 2013;44(4):413–24.
Mashaly M, Hendricks G, Kalama M, Gehad A, Abbas A, Patterson P. Effect of heat stress on production parameters and immune responses of commercial laying hens. Poult Sci. 2004;83(6):889–94.
Kültz D. Molecular and evolutionary basis of the cellular stress response. Annu Rev Physiol. 2005;67:225–57.
Hosseini-Vashan S, Golian A, Yaghobfar A. Growth, immune, antioxidant, and bone responses of heat stress-exposed broilers fed diets supplemented with tomato pomace. Int J Biomet. 2015;1:10.
Chen W, Melamed ML, Abramowitz MK. Serum bicarbonate and bone mineral density in US adults. Am J Kid Dis. 2015;65(2):240–8.
Reiner G, Clemens N, Fischer R, Köhler F, Berge T, Hepp S, et al. Mapping of quantitative trait loci for clinical–chemical traits in swine. An Genet. 2009;40(1):57–64.
Deeb N, Lamont S. Genetic architecture of growth and body composition in unique chicken populations. J Hered. 2002;93(2):107–18.
NRC. Nut Requ Poult. National Research Council. USA: National Academy Press Washington; 1994.
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(1):59.
Institute S. JMP statistical discovery software. 2000.
Gilmour AR, Gogel B, Cullis B, Thompson R, Butler D. ASReml user guide release 3.0. VSN Int Ltd. 2009.
Fernando R, Garrick D. GenSel user manual for a portfolio of genomic selection related analyses. Ames: Animal Breeding and Genetics, Iowa State University; 2008.
Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Res. 2015;43(D1):D662–9.
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 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 partly supported by a USDA National Needs Fellowship, Award #2011-38420-20050.
The authors have no competing interests.
AVG participated in sample collection, carried out all statistical analyses and drafted the manuscript. 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.