- Research article
- Open Access
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 Genetics volume 14, Article number: 22 (2013)
Comparisons of quantitative trait loci (QTL) for growth and parameters of growth curves assist in understanding the genetics and ultimately the physiology of growth. Records of body weight at 3, 6, 12, 24, 48 and 72 weeks of age and growth rate between successive age intervals of about 500 F2 female chickens of the Roslin broiler-layer cross were available for analysis. These data were analysed to detect and compare QTL for body weight, growth rate and parameters of the Gompertz growth function.
Over 50 QTL were identified for body weight at specific ages and most were also detected in the nearest preceding and/or subsequent growth stage. The sum of the significant and suggestive additive effects for bodyweight at specific ages accounted for 23-43% of the phenotypic variation. A single QTL for body weight on chromosome 4 at 48 weeks of age had the largest additive effect (550.4 ± 68.0 g, 11.5% of the phenotypic variation) and a QTL at a similar position accounted 14.5% of the phenotypic variation at 12 weeks of age. Age specific QTL for growth rate were detected suggesting that there are specific genes that affect developmental processes during the different stages of growth. Relatively few QTL influencing Gompertz growth curve parameters were detected and overlapped with loci affecting growth rate. Dominance effects were generally not significant but from 12 weeks of age they exceeded the additive effect in a few cases. No evidence for epistatic QTL pairs was found.
The results confirm the location for body weight and body weight gain during growth that were identified in previous studies and were consistent with QTL for the parameters of the Gompertz growth function. Chromosome 4 explained a relatively large proportion of the observed growth variation across the different ages, and also harboured most of the detected QTL for Gompertz parameters, confirming its importance in controlling growth. Very few QTL were detected for body weight or gain at 48 and 72 weeks of age, probably reflecting the effect of differences in reproduction and random environmental effects.
Combining growth models with QTL mapping facilitates the understanding of the genetics underlying physiological aspects of quantitative traits . Whereas there are many reports of quantitative trait loci (QTL) for body weight and several for rate of growth [2–8] before sexual maturity, there is relatively little information on QTL for body weight at maturity or growth curve parameters. This is partly due to the fact that most meat chickens are slaughtered long before they reach maturity and growth curves cannot usefully be fitted to the data. Furthermore, several weights on each bird are required for reliable model fitting.
Detection of QTL influencing the parameters of the growth curve has been reported for mice , sheep , dairy cattle [11, 12] and pigs . QTL based on the analysis of growth curve models detects similar QTL as age-specific body weight provided the chosen growth functions fit the data satisfactorily . Understanding the biology of the model parameters and their relationships can assist in developing a breeding strategy to modify the shape of the growth curve that would not be feasible with age specific QTL. It has been demonstrated that parameters of the growth curve are heritable and that the curve can be modified through selection on bodyweight at different ages [14–16]. This might be a useful strategy to minimise weight-related problems associated with rapid early growth, such as skeletal disease, for example.
A study was conducted to identify QTL for body weight, growth rate and parameters of the growth curve for chickens from 3 to 72 weeks of age of an F2 cross of a broiler male line and a White Leghorn layer. A Gompertz growth curve was fitted to the body weights for each individual and the estimated parameters of the growth curve for each bird were analysed. Age-specific body weights and body weight gains were also analysed for comparison with QTL for growth curve parameters.
Trait means, variation and phenotypic correlations
Overall means, standard deviations and range of body weight at different ages and estimates of parameters of the Gompertz equation are presented in Additional file 1: Table S1. Mature body weight averaged 3.9 kg (range 2.0 – 5.8 kg); estimated age at the point of inflection of the growth curve was 64 d with a 95% confidence interval of 50–78 d.
Phenotypic correlations between body weights, growth rates and parameters of the Gompertz model are presented in Additional file 2: Table S2. Phenotypic correlations between weights at successive ages were relatively high, whereas correlations between growth rates at different ages were low, generally accounting for less than 10% of the variation in growth rate at later ages. There was a wide range of correlations among the parameters of the Gompertz equation. Surprisingly, correlations with mature weight (Wa) were low; high negative correlations existed between Ti and K and between W0 and both K and L whereas there were large positive correlations between K and L.
Analyses of epistasis
There was no evidence of epistatic QTL pairs for any of the growth-related traits except that a suggestive QTL pair involving chromosomes 2 and 3 was detected for growth rate from 24 to 48 weeks of age (results not shown).
Body weight QTL
QTL for body weight at different ages are presented in Table 1. A total of 31 genome significant QTL, including those that affected different traits, were identified on chromosomes 1, 2, 3, 4, 6, 8, 11, 27, 28 and Z. Two significant QTL for body weight were detected on chromosomes 1 (at 12 weeks), 4 (at 6 and 12 weeks), and 3 (at 48 weeks). A further 25 suggestive QTL were identified for body weight at different ages. More body weight QTL were detected for the early growth stages (3 – 12 weeks of age) than for mature growth (after 24 weeks of age) (Table 1).
The largest proportion of the phenotypic variation explained by a QTL was 14.5% for 12 week body weight at 177 cM on chromosome 4 (Table 1). The contribution of most of the QTL varied across ages e.g. the chromosome 4 QTL contribution peaked just before sexual maturity that occurred at an average of 19 weeks . Several QTL affected body weight at successive ages, as expected from the phenotypic correlations (Additional file 2: Table S2). The proportion of phenotypic variation explained by significant and suggestive QTL for body weight at 3, 6, 12, 24, 48 and 72 weeks respectively was 23, 43, 39, 28, 35 and 24%.
Most of the significant QTL had positive additive effects implying that the alleles from the broiler line increased body weight. Dominance effects were generally not significant and in cases where they were significant, had negative values. A QTL for body weight at 48 weeks of age segregating on chromosome 4 at 153 cM had the largest additive effect (550.4 ± 68.0 g) and explained 11.5% of the phenotypic variation. The largest dominance effect (-194.7 ± 87.1 g) was for a QTL on chromosome 9 at 81 cM and the QTL effect accounted for 1.7% of the phenotypic variation for body weight at 48 weeks of age.
Growth rate QTL
A total of 12 significant QTL for growth rate at different ages were identified on chromosomes 1, 2, 3, 4, 8 and 27 (Table 2). Two highly significant QTL for Gr3-6 were detected on chromosome 4 and two significant QTL on chromosome 3 for Gr24-48. More QTL were detected for the early growth rates and only suggestive QTL were identified for Gr24-48. One of the QTL for Gr3-6 on chromosome 4 explained the highest proportion of the phenotypic variation (7.2%) among growth rate QTL and jointly with the second QTL on this chromosome explained 10.9% of the phenotypic variation for this trait. Most QTL acted additively but some significant dominance effects were found for QTL on chromosomes 1, 8 and 3 respectively for Gr6-12, Gr12-24 and Gr24-48.
The proportion of phenotypic variation explained by significant and suggestive QTL for Gr3-6, Gr6-12, Gr12-24, Gr24-48 and Gr48-72 respectively was 33, 15, 5, 14 and 10%.
Gompertz curve parameter QTL
The QTL results for the lnGompertz curve parameter estimates are presented in Table 3. Except for mature weight, most QTL for growth curve parameter estimates were merely suggestive. Several QTL detected were for asymptotic body weight (WA) on chromosomes 2, 4, 8 and 27, and one QTL was identified for age at maximum growth, Ti, on chromosome 11. Relatively few suggestive QTL were segregating for for K, L and W0, The proportion of phenotypic variation explained by significant and suggestive QTL for WA, Ti, K, L and W0 was 30, 8, 7, 3 and 5%.
Comparisons of three QTL detection approaches
A comparison of the locations of QTL for body weight, rate of body weight gain and the Gompertz parameters is presented in Table 4. Growth rate QTL co-locate with most of the bodyweight QTL at the respective ages, as expected from the part-whole relationship of these traits. However, the number of detected QTL for growth rate is far fewer than those detected for body weight (Tables 1 and 2).
Different QTL on chromosome 8 at approximately 60 cM affected body weight at 3, 6 and 12 weeks, mature weight and growth rate between 12 and 24 weeks of age. A QTL at 12–15 cM affected body weight at 24 weeks, growth from 24–48 weeks and WA. This QTL region was also the only significant QTL for Gr12-24 which spans the age at the onset of sexual maturity. QTL were detected for growth rate across the different growth phases on chromosome 4, and to a lesser extent chromosome 3, except at 12–24 weeks of age corresponding to the age of onset of sexual maturity.
Most of the QTL for growth curve parameter estimates were for asymptotic body weight (WA) and were similar to QTL on chromosomes 2, 3, 4, 8, 15 and 27 for body weight at different ages. The suggestive QTL segregating on chromosome 4 at 144 cM for the rate parameter (L) is close to the location of the QTL detected for WA at 166 cM and is similar to the QTL segregating across different ages at 142–183 cM on chromosome 4. The QTL for L at 541 cM on chromosome 1 is located in the same region flanked by markers ROS0081 and LEI0079 as body weight QTL at 3, 6 and 12 weeks of age (Table 1).
No evidence for epistasis
Failure to detect significant epistasis for the evaluated traits was unexpected given that other studies working on similar traits detected epistatic QTL pairs for body weight at early growth stages and for Gompertz curve parameters [4, 7, 18, 19]. Most approaches to detect epistasis are subject to a high level of false positive results  and the detection of only one suggestive epistatic QTL pair could be due to the stricter thresholds enforced in this analysis compared to those reported by other studies, the low density of marker coverage or the relatively small size of the population.
Bodyweight QTL at Specific Ages
Body weight QTL at 3 weeks on chromosomes 1, 4, 13 and Z (Table 1) were similar to those detected earlier in a similar population raised as broilers for maximum growth . In a previous report, chromosomes 1, 2, 3, 6 and 11 were shown to harbour QTL for body weight at 46 days  and also appear in the list of significant and suggestive QTL for weight at 6 weeks of age in this study (Table 1). The positions of the QTL for 6 week weight on chromosome 3 at 235 cM and 161 cM on chromosome 4 are similar to the QTL at 252 and 149 cM respectively on chromosomes 3 and 4 reported by Jacobsson et al. . Most of the QTL for this trait are similar to those identified on chromosomes, 1, 2, 4, 6, 8, 13 and Z in a parallel broiler study .
QTL detected in the latter study for body weight at 9 weeks on chromosomes 1, 2, 4, 6, 8, and 13 are similar to those at 12 weeks of age. Positions for the QTL on chromosomes 4 and 13 respectively were at 177 and 15 cM in the earlier report and at 177 and 7 cM in this study. These are probably the same QTL because of the high correlations expected between body weight at 9 and 12 weeks of age. The significant 12 week body weight QTL on chromosomes 1, 4 and 27 are at similar locations as carcass weight QTL in a related broiler experiment .
QTL for body weight at 6 weeks on chromosomes 2, 3 (two locations), 4 and 6 were also identified at 48 weeks of age. Whereas body weight at two ages includes a part-whole component and the phenotypic correlation was only 0.52 (Additional file 2: Table S2), the results emphasise the importance of early growth in determining adult body weight.
Growth rate QTL
QTL for Gr3-6 on chromosomes 1, 2, 3, 4 and 13 in this study were similar to those reported for growth from 2 to 4 weeks on chromosomes 1 and 2 [4, 23]; and the QTL for Gr6-12 on chromosomes 1 and 3 are similar to those previously published for growth from 6 to 12 weeks .
QTL for Gompertz parameters
The Laird form of the Gompertz equation is a function of initial growth rate, relative growth rate at the point of inflection (the rate of exponential decay of the relative growth rate) and body weight at time t0 compared to the original Gompertz model which is a function of mature body weight . It was chosen to maximise the number of data points available for analysis (requiring the estimation of only 3 parameters) and provided an excellent fit to the data (results not shown).
The Gompertz parameter estimates (Additional file 1: Table S1) fall within the range reported in the literature: WA: 2483 – 5698 g, L: 0.0908 – 0.141 g/d, K: 0.0224 – 0.031, Ti: 42.2 – 63 d, W0: 39.8 – 64 g [14, 25–27]. The similar locations of significant and suggestive QTL for body weight and Gompertz parameters (Table 4) confirm assertions made earlier that these parameters are genetically determined and can be exploited to improve traits through selection [1, 8].
Whereas sexual maturity occurs earlier in the male-line broilers than the White Leghorn layers used in this experiment , there was evidence of only one significant QTL for the point of inflection of the growth curve. The broiler allele of this QTL decreased the point of inflection by 1.8 d or 3.6 d in the homozygous state and if the three significant and suggestive effects from Table 3 are summed the effect is substantial (almost 7 d or 14 d in the homozygous state).
We chose to use a constant 14 h light (L): 10 h dark (D) photoschedule rather than use a typical 8 L:16D during rearing gradually increasing to, for example, 16 L:8D during lay starting when the birds are between 16 and 22 weeks of age. This was done in order to avoid compromising the growth of the birds by artificially photostimulating them at the same age when individuals would be at different physiological states caused at least in part by genetic differences. This practice therefore ensured that all the birds received the same photostimulation and that growth per se was not affected by differences in the reproductive responses to increasing photoperiods. It also ensured that growth was not constrained by limiting opportunity for ad libitum feed consumption to only 8 h per day. Significant QTL for mature body weight (WA) were detected on chromosomes 2, 4, 8 and 27 while suggestive QTL were detected on chromosomes 3, 7, 9 and 15. Asymptotic or mature body weight QTL have been reported by Le Rouzic et al. (2008) on chromosomes 2 and 27 and additional QTL that were not detected in our study on chromosome 1, 6 and 11 . A similar QTL for age at the point of inflection was detected on chromosome 11 in both studies whereas suggestive QTL for this trait on chromosomes 2 and 4 contrast with significant QTL on chromosomes 1, 12 and 20 . No genome wide significant QTL were detected for W0, K or L. W0 may be determined more by egg size and therefore maternal QTL than by QTL inherited by the chick. Genetic variation for K and L may exist but be explained more by the other parameters, particularly mature weight. Changing growth rate while maintaining acceptable hatch and mature weights may be a desirable goal in commercial meat production systems and growth curves have been modified by differential selection on early growth in chickens and quail [28, 29]. Selection based on growth curve parameters is therefore likely to be effective but difficult to implement because of the length of time required to obtain body weights at older ages on which to base parameter estimates. The results suggest that relatively few areas of the genome may contain the genes largely responsible for controlling growth curves suggesting that whole genome selection based on high density SNP chips could ameliorate this problem. However, in contrast to the successful selection experiments noted above, Ibanez and Blasco  suggest that genomic selection to change the growth curve will be difficult and require constant re-evaluation of the associations between the SNPs and the genes determining curve parameters.
The architecture of growth QTL
Examination of Table 4 suggests that relatively few chromosomal locations affect growth to an extent that they can be detected in this population. Taking the large confidence intervals into consideration (Tables 1, 2, 3), the results suggest that there were at least 15 such locations: 2 on each of chromosomes 1, 2, 3 and 4, one each on chromosomes 6, 8, 9, 11 and Z with an additional 2 on microchromosomes 11 and 27. These are generally consistent with other experiments and provide a smaller list of QTL in the search for causative mutations compared with the list of QTL for body weight identified at different ages such as in Table 1.
Evidence of a QTL at 166 cM on chromosome 4 for mature body weight had the largest additive effect and explained the highest proportion (9.5%) of the phenotypic variation. The results generally confirm earlier observations about the critical role of chromosome 4 in controlling growth and other traits of economic importance [3, 31, 32].
QTL detected for growth rate intervals (Table 2) were generally similar to those for body weight at the corresponding ages (Table 1). Both the body weight and growth rate approaches identified more significant QTL than analysis of the Gompertz curve parameters but all methods identified a significant QTL for adult body weight (WA) on chromosomes 2, 4, and 8.
More QTL were detected for growth before sexual maturity than for later growth and suggest that genetic variation for growth is more important during early life (Tables 1 and 2). This may in part be associated with differential development of the reproductive organs and fat deposition as the birds approach sexual maturity and by changes in reproductive status with increasing age. These effects are likely to decrease the apparent importance of growth QTL whereas those QTL associated with reproductive senescence or fatness may become more important and overshadow the growth QTL identified during the rearing of the birds. The latter are more likely to be associated with growth of the lean tissue mass (muscle, bone and nervous tissues) and after peak rates of lay by QTL for fatness. Alternatively the environmental variance may become larger with age and time in the cages. These observations are consistent with the lack of correlation between WA and growth from 48–72 weeks (Additional file 2: Table S2). Further research will be necessary to determine the relative role of non-growth QTL on variation in adult weight including differences in fatness.
Age specific body weight QTL were detected for each growth stage which supports earlier observations that there are different genes and gene actions involved in growth at different developmental stages [1, 4, 8, 9, 33]. Chromosomes 3, 4 and 8 had QTL involved with body weight throughout the lifetime of the birds. Chromosome 11 harboured QTL involved mainly in very early growth and those on chromosome 1, 13 and Z from 3 to 24 weeks age. QTL on chromosomes 15, 27, and 28 mainly affected growth from 24 weeks. The few QTL that were detected at older ages (48 and 72 weeks) compared to earlier periods may reflect large environmental effects, particularly with respect to individual differences in the age that egg laying started to decline or cease altogether, and the accompanying changes in fat deposition, as discussed above.
Moderately high negative correlations between early body weight and age at inflection of the growth curve are consistent with previous research on the age at the onset of puberty which showed that large body weight QTL were associated with early onset of egg laying . Thus large, faster growing birds will tend to mature earlier and therefore the age at inflection of the growth curve will also be earlier than in slower growing, smaller birds.
A large number of significant QTL for body weight at specific ages was detected and most of the identified QTL affected growth rate in the preceding and/or subsequent growth stages but overall a smaller number of about 15 chromosomal locations have a substantial effect on growth in this cross. Most of the detected QTL at ages before sexual maturity have been reported in other studies and the results confirm many age specific QTL. QTL influencing Gompertz parameters were detected and these QTL also overlapped with loci affecting growth and carcass traits reported in other studies. Age specific growth QTL may be associated with specific genes and gene actions that orchestrate developmental processes during different stages of growth. Some loci featured predominantly in early growth to maturity. However, there was little or no evidence for epistasis for any growth or body weight trait. Chromosome 4 explained a relatively large proportion of the observed growth variation across the different ages, and also harboured most of the detected QTL for Gompertz parameters, confirming its importance in controlling growth. Very few QTL were detected for body weight or gain at 48 and 72 weeks of age, probably reflecting the effect of differences in reproductive decline with age and the accumulation of random environmental effects.
Animals and husbandry
The F2 population was created by crossing two males and two females from the broiler male line with two females and two males from the White Leghorn line to produce an F1 generation that was intercrossed as described . Briefly, 8 males and 32 females of the F1 generation were selected and mated in a balanced mating scheme to produce over 2000 F2 birds. The female chicks from 9 hatches were reared in floor pens and moved to individual cages (40 cm wide × 45 cm deep × 80 cm high) at 12 weeks of age. The birds were fed ad libitum on a conventional poultry diet and had permanent access to water. The birds were exposed to a constant photoperiod of 14 hours per day from hatch to the end of the experiment.
The body weight of each bird was recorded at 3, 6, 12, 24, 48 and 72 weeks of age. Body weights at specific ages, growth rates between successive age intervals and growth curve parameter estimates for each individual were analysed as different traits. Growth rate at each age interval was derived by dividing the body weight gained in the interval by the difference in age (d) and were denoted as Gr3-6, Gr6-12, Gr12-24, Gr24-48 and Gr48-72 respectively for the daily weight gained from 3–6, 6–12, 12–24, 24–48 and 48–72 weeks of age.
All animal husbandry and procedures were conducted under legislation to protect the welfare of animals and were licensed after ethical review.
Growth curve modelling
The Laird form of the Gompertz curve equation as described by Aggrey , was fitted to data from 453 individuals that had a minimum of 4 data points after removing extreme outliers (>3 SD from the mean). The equations for deriving the parameters of the Laird form of the Gompertz curve were defined following  as
where Wt is the body weight of a bird at time t, W0 the estimated initial hatching weight, L the rate of growth (g/d) at time t = 0 and K the rate of exponential decay of the relative growth rate. The following derived parameters were computed: age at the point of inflection, T i where Ti = (1/K)log(L/K) and asymptotic or mature body weight, WA; where WA = W0.exp(L/K). The growth curves were fitted using the Genstat program (http://www.vsni.co.uk/software/genstat) and the parameter estimates of the growth curves were extracted for each bird. The model converged for all the data points analysed.The estimated parameters were not normally distributed and to approximate normality the natural logarithm (ln) of the Laird form of the Gompertz curve was adopted.
DNA analysis and map construction
DNA was extracted from blood samples as described previously . Genotyping was performed using a total of 106 microsatellite markers covering 25 autosomal linkage groups and the Z chromosome (Additional file 3: Table S3). Genetic linkage maps were constructed using the prepare, flips and fixed options of the CRIMAP programme .
Mapping and significance testing were conducted by the interval mapping method for QTL analysis adapted for epistasis detection in GridQTL  as described earlier . Significance thresholds for detection of single QTL with significant marginal effects were determined through 5000 permutations  and 1000 bootstraps were used to generate 95% confidence intervals for the QTL positions . An F value greater than the P > 0.05 and P > 0.01 experiment-wide threshold values were used to identify a significant and highly significant QTL  respectively. QTL that achieved an F ratio exceeding the P > 0.05 chromosome-wide threshold were considered to be suggestive. Significance testing for epistatic pairs used F ratio tests for model comparisons in a nested test framework following .
Different models with additive, dominance and parent-of-origin genetic effects with family and pen as fixed effects (hatch was confounded with pen) were evaluated in preliminary analyses. There was no evidence for a parent-of-origin effect  and parent of origin effects were ignored in subsequent analyses. An additive plus dominance model was evaluated for all chromosomes simultaneously except for the Z chromosome which was analysed with an additive genetic effects model. The epistasis analysis did not include the Z chromosome.
Wu W, Zhou Y, Li W, Mao D, Chen Q: Mapping of quantitative trait loci based on growth models. Theor Appl Genet. 2002, 105 (6–7): 1043-1049.
Ambo M, Moura ASAMT, Ledur MC, Pinto LFB, Baron EE, Ruy DC, Nones K, Campos RLR, Boschiero C, Burt DW: Quantitative trait loci for performance traits in a broiler × layer cross. Anim Genet. 2009, 40 (2): 200-208. 10.1111/j.1365-2052.2008.01824.x.
Sewalem A, Morrice DM, Law A, Windsor D, Haley CS, Ikeobi CO, Burt DW, Hocking PM: Mapping of quantitative trait loci for body weight at three, six, and nine weeks of age in a broiler layer cross. Poult Sci. 2002, 81 (12): 1775-1781.
Carlborg O, Hocking PM, Burt DW, Haley CS: Simultaneous mapping of epistatic QTL in chickens reveals clusters of QTL pairs with similar genetic effects on growth. Genet Res. 2004, 83 (3): 197-209. 10.1017/S0016672304006779.
Terčič D, Holcman A, Dovč P, Morrice DR, Burt DW, Hocking PM, Horvat S: Identification of chromosomal regions associated with growth and carcass traits in an F3 full sib intercross line originating from a cross of chicken lines divergently selected on body weight. Anim Genet. 2009, 40 (5): 743-748. 10.1111/j.1365-2052.2009.01917.x.
Uemoto Y, Sato S, Odawara S, Nokata H, Oyamada Y, Taguchi Y, Yanai S, Sasaki O, Takahashi H, Nirasawa K: Genetic mapping of quantitative trait loci affecting growth and carcass traits in F2 intercross chickens. Poult Sci. 2009, 88 (3): 477-482. 10.3382/ps.2008-00296.
Le Rouzic A, Alvarez-Castro JM, Carlborg O: Dissection of the genetic architecture of body weight in chicken reveals the impact of epistasis on domestication traits. Genetics. 2008, 179 (3): 1591-1599. 10.1534/genetics.108.089300.
Gao Y, Hu XX, Du ZQ, Deng XM, Huang YH, Fei J, Feng JD, Liu ZL, Da Y, Li N: A genome scan for quantitative trait loci associated with body weight at different developmental stages in chickens. Anim Genet. 2006, 37 (3): 276-278. 10.1111/j.1365-2052.2006.01428.x.
Long FEI, Qing Chen Y, Cheverud JM, Wu R: Genetic mapping of allometric scaling laws. Genet Res. 2006, 87 (03): 207-216. 10.1017/S0016672306008172.
Hadjipavlou G, Bishop SC: Age-dependent quantitative trait loci affecting growth traits in Scottish Blackface sheep. Anim Genet. 2009, 40 (2): 165-175. 10.1111/j.1365-2052.2008.01814.x.
Lund MS, Sorensen P, Madsen P, Jaffrézic F: Detection and modelling of time-dependent QTL in animal populations. Genet Sel Evol. 2008, 40 (2): 177-194. 10.1186/1297-9686-40-2-177.
Rodriguez-Zas SL, Southey BR, Heyen DW, Lewin HA: Detection of quantitative trait Lloci influencing dairy traits using a model for longitudinal data. J Dairy Sci. 2002, 85 (10): 2681-2691. 10.3168/jds.S0022-0302(02)74354-3.
Varona L, Gómez-Raya L, Rauw WM, Ovilo C, Clop A, Noguera JL: The value of prior information for detection of QTL affecting longitudinal traits: an example using Von Bertalanffy growth function. J Anim Breed Genet. 2005, 122 (1): 37-48. 10.1111/j.1439-0388.2004.00477.x.
Mignon-Grasteau S: Genetic parameters of growth curve parameters in male and female chickens. Br Poult Sci. 1999, 40 (1): 44-51. 10.1080/00071669987827.
Grossman M, Bohren BB: Logistic growth curve of chickens: heritability of parameters. J Hered. 1985, 76 (6): 459-462.
Mignon-Grasteau S, Beaumont C, Ricard FH: Genetic analysis of a selection experiment on the growth curve of chickens. Poult Sci. 2001, 80 (7): 849-854.
Carlborg O, Kerje S, Schutz K, Jacobsson L, Jensen P, Andersson L: A global search reveals epistatic interaction between QTL for early growth in the chicken. Genome Res. 2003, 13 (3): 413-421. 10.1101/gr.528003.
Ankra-Badu G, Shriner D, Le Bihan-Duval E, Mignon-Grasteau S, Pitel F, Beaumont C, Duclos M, Simon J, Porter T, Vignal A: Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate. BMC Genom. 2010, 11 (1): 107-10.1186/1471-2164-11-107.
Wei WH, Knott S, Haley CS, de Koning DJ: Controlling false positives in the mapping of epistatic QTL. Heredity. 2010, 104 (4): 401-409. 10.1038/hdy.2009.129.
Jacobsson L, Park H-B, Wahlberg P, Fredriksson R, Perez-Encison M, Siegel PB, Andersson L: 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-125. 10.1017/S0016672305007767.
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-442. 10.1080/00071660500158055.
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 (10): 1700-1711.
Koncagul S, Cadirci S: Analysis of growth of broilers with restricting and unrestricting initial body weight in Gompertz-Laird model in different environments. Ital J Anim Sci. 2010, 9 (1): 20-25.
Darmani Kuhi H, Kebreab E, Lopez S, France J: An evaluation of different growth functions for describing the profile of live weight with time (age) in meat and egg strains of chicken. Poult Sci. 2003, 82 (10): 1536-1543.
Podisi BK, Knott SA, Dunn IC, Law AS, Burt DW, Hocking PM: Overlap of quantitative trait loci for early growth rate, and for body weight and age at onset of sexual maturity in chickens. Reproduction. 2011, 141 (3): 381-389. 10.1530/REP-10-0276.
Aggrey SE: Comparison of three nonlinear and spline regression models for describing chicken growth curves. Poult Sci. 2002, 81 (12): 1782-1788.
Mignon-Grasteau S, Beaumont C, Ricard FH: Genetic analysis of a selection experiment on the growth curve of chickens. Poult Sci. 2001, 80 (7): 849-854.
Hyankova L, Knizetova H, Dedkova L, Hort J: Divergent selection for shape of growth curve in Japanese quail. 1. Responses in growth parameters and food conversion. Br Poult Sci. 2001, 42 (5): 583-589. 10.1080/00071660120088371.
Ibanez-Escriche N, Blasco A: Modifying growth curve parameters by multitrait genomic selection. J Anim Sci. 2011, 89 (3): 661-668. 10.2527/jas.2010-2984.
Schreiweis MA, Hester PY, Moody DE: Identification of quantitative trait loci associated with bone traits and body weight in an F2 resource population of chickens. Genet Sel Evol. 2005, 37 (6): 677-698. 10.1186/1297-9686-37-7-677.
Tuiskula-Haavisto M, Honkatukia M, Vilkki J, de Koning DJ, Schulman NF, Maki-Tanila A: Mapping of quantitative trait loci affecting quality and production traits in egg layers. Poult Sci. 2002, 81 (7): 919-927.
Wu R, Ma C-X, Lin M, Casella G: A general framework for analyzing the genetic architecture of developmental characteristics. Genetics. 2004, 166 (3): 1541-1551. 10.1534/genetics.166.3.1541.
Aggrey SE: Comparison of three nonlinear and spline regression models for describing chicken growth curves. Poult Sci. 2002, 81 (12): 1782-1788.
Green P, Falls K, Crooks S: Cri-map (Version 2.4). 1990, St. Louis, Mo. USA: Washington University, School of Medicine, 24
Churchill GA, Doerge RW: Empirical Threshold Values For Quantitative Trait Mapping. Genetics. 1994, 138 (3): 963-971.
Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143: 1013-1020.
Kruglyak L, Lander ES: A nonparametric approach for mapping quantitative trait loci. Genetics. 1995, 139 (3): 1421-1428.
Knott SA, Marklund L, Haley CS, Andersson K, Davies W, Ellegren H, Fredholm M, Hansson I, Hoyheim B, Lundstrom K: Multiple marker mapping of quantitative trait loci in a cross between outbred wild boar and Large White pigs. Genetics. 1998, 149 (2): 1069-1080.
The authors gratefully acknowledge the technical assistance Robert Bernard and Graeme Robertson, and the use of UK National Grid Service in carrying out this work. Aviagen Ltd, Newbridge, Midlothian kindly provided the broiler chicks used in these experiments. The research was based on experiments funded by Defra and BBSRC. The Roslin Institute is funded by a core strategic grant from the BBSRC and B.K. Podisi was sponsored by The Government of Botswana.
The authors declare that they have no competing interests.
PMH, SK and DWB contributed to the planning of the experiment. PMH contributed in the data collection. BKP and PMH performed the data analysis, discussed the results and drafted the manuscript. All authors read, provided comments and agreed on the contents of the manuscript.