- Research article
- Open Access
Search for QTL affecting the shape of the egg laying curve of the Japanese quail
BMC Genetics volume 7, Article number: 26 (2006)
Egg production is of critical importance in birds not only for their reproduction but also for human consumption as the egg is a highly nutritive and balanced food. Consequently, laying in poultry has been improved through selection to increase the total number of eggs laid per hen. This number is the cumulative result of the oviposition, a cyclic and repeated process which leads to a pattern over time (the egg laying curve) which can be modelled and described individually. Unlike the total egg number which compounds all variations, the shape of the curve gives information on the different phases of egg laying, and its genetic analysis using molecular markers might contribute to understand better the underlying mechanisms. The purpose of this study was to perform the first QTL search for traits involved in shaping the egg laying curve, in an F2 experiment with 359 female Japanese quail.
Eight QTL were found on five autosomes, and six of them could be directly associated with egg production traits, although none was significant at the genome-wide level. One of them (on CJA13) had an effect on the first part of the laying curve, before the production peak. Another one (on CJA06) was related to the central part of the curve when laying is maintained at a high level, and the four others (on CJA05, CJA10 and CJA14) acted on the last part of the curve where persistency is determinant. The QTL for the central part of the curve was mapped at the same position on CJA06 than a genome-wide significant QTL for total egg number detected previously in the same F2.
Despite its limited scope (number of microsatellites, size of the phenotypic data set), this work has shown that it was possible to use the individual egg laying data collected daily to find new QTL which affect the shape of the egg laying curve. Beyond the present results, this new approach could also be applied to longitudinal traits in other species, like growth and lactation in ruminants, for which good marker coverage of the genome and theoretical models with a biological significance are available.
From the same natural oviposition cyclic process, wild birds generally lay only one or a few clutches of up to two dozen eggs per year  whereas quail  and chicken  from highly selected commercial or experimental lines may produce several hundred eggs over the same period. Yet, their egg production changes with time, increasing from the onset of egg laying to peak production, then plateauing and finally declining steadily. It may be described by drawing an egg laying curve which can be modelled by theoretical equations [e.g.: [4–7]] with biologically meaningful parameters. Another, more global, measure of egg production is the total number of laid eggs which has been extensively studied from the standpoint of quantitative genetics [8, 9], and for which QTL were recently detected in chicken [10, 11] and Japanese quail . This trait, however, conveys little information for explaining the shape of the egg laying curve and its variation. On the contrary, estimations of the parameters of the individual egg laying curves fitted to the observed egg production might be used as phenotypes, and searched for underlying QTL, but this approach has not been tested until now.
The objective of the present work was to produce the first set of QTL detected in birds for the parameters which determine the shape of the egg laying curve. For that purpose, in an F2 experiment with two quail lines selected for early egg production  and for high duration of tonic immobility , which had led already to the detection of QTL for total egg number, age at first egg and clutch length , the individual egg laying curves were fitted to egg production using three theoretical models [4, 6, 7]. The Gifu microsatellite marker panel  and the microsatellite map  developed from the same F2 reference families were used to detect QTL for the parameters of the curves.
Egg laying curves
Table 1 shows the equations of the three models, the estimations of their parameters, and their overall fit to the F2 data from the 359 female quail. It also lists the estimations of the same parameters for the two F0 lines (DD and LTI) used to produce the F1 generation, parent of the F2. Observed and predicted egg laying curves of the F2 population are in Figure 1 (Model 1 for egg laying rate) and Figure 2 (Models 2 and 3 for 3-week egg production).
Table 2 lists the parameters of the curves for which the detection of QTL was performed, and gives their biological meaning. It also shows the corresponding statistics (number, mean, and range) from the set of the individual curves adjusted to the data. Five, 7 and 51 curves could not be obtained by nonlinear regression under Models 1, 2 and 3, respectively.
Table 3 shows the location of the chromosome-wide significant QTL, their position on the chromosome, the maximum F value obtained at this position, their genetic effects, the reduction of the residual variance obtained by fitting a QTL at this location, and the corresponding genome-wide significance.
Out of the 8 QTL found to be chromosome-wide significant (p c <0.05) on 5 of the 12 autosomes scanned in this study, one was genome-wide suggestive (p g <0.20) and none was genome-wide significant (p g <0.05). One chromosome-wide significant QTL for the parameter t 2 related to the ascending phase of the egg laying curve in Model 2 was found on CJA13. Another QTL was close to chromosome-wide significance on CJA06 for the parameter y P which describes the plateaued phase of the curve in the same model. Chromosome-wide significant QTL for the parameters associated with persistency in Models 1 and 3 were found on CJA10 and CJA14 for b, and on CJA05 and CJA10 for k2. QTL were also found on CJA10 for the two parameters, c and d, of Model 1 which are related to sexual maturity.
In most cases, the genetic effect of the QTL was dominant, but the QTL was strictly additive in two instances, for y P and t 2 . The reduction of the residual variance due to fitting a QTL was small and did not exceed 2.2%.
Egg laying curves
For all three populations, the F2 and the two F0 lines, the three models provided a good overall fit (R2) to the egg laying data as they accounted for around 90% of the total variation. Similar figures have been obtained on quail with Model 1 , but it seems that Model 3 underestimated the peak production of the F2 (Figure 2). Moreover, it was more difficult to adjust this model to individual curves as indicated by the high proportion (14%) of instances for which the numerical estimation of the parameters did not converge to an acceptable set of values. Of course, the coefficient of determination R2 of individual curves fitted for each F2 quail was variable and often lower (in the range 99–58%) than for the global curve. Correspondingly, the variation of the biological parameters of the individual curves was much larger because they were estimated from 19 "observations" or less (only one value by bird and 3-week period). Because of the individual between-bird variations in the actual shape of the laying curves, it might be difficult to achieve more uniform goodness of fit for the within-bird curve parameters than that obtained in the present study, and as observed also with chicken [6, 7]. This means that the precision on the phenotypes (the estimated parameters of the curve) used for the QTL analysis varied between birds, possibly making the detection of QTL more difficult to achieve than for traits measured directly.
The majority of the chromosome-wide significant QTL corresponding to different parts of the laying curve obtained for the three models showed some dominance, which could be expected since there is consistent heterosis for egg production in poultry . The significant genome-wide additive QTL for total egg number found in the previous analysis of the same set of quail  was mapped on CJA06 at the same position as the additive QTL found for y p (the number of eggs laid in 3 weeks at the time of maximum egg laying) in the present study. It is unlikely this convergence was coincidental, and it is an indication that the QTL which were only detected at the chromosome-wide level in the present work might truly correspond to regions with genes with some effect on the shape of the laying curve. Also, it is interesting to notice that the QTL analysis of the curves obtained under Model 1(for egg laying rate) and Model 3 (for 3-week egg number) led to the detection of a QTL for the decreasing phase of the curve (related to persistency) on CJA10. To our knowledge, there is no published report on similar QTL analyses of curves fitted to repeated observations of the same phenotype. Yet, this approach might be a good alternative to splitting the egg production into consecutive periods and searching for QTL for each period .
For the first time, a genome scan was carried out to detect QTL for the parameters of curves modelling a biological process estimated by nonlinear regression and used as phenotypes for the analysis. The potential interest of this approach is in the integration of all the information to discover genes which affect the expression of the trait through time, as illustrated in the present work by the egg laying curve. Only chromosome-wide QTL were obtained in this study, which might be related to the limited value of the available models for describing accurately the individual egg laying curve and to the "imprecise" nature of the statistics used as a phenotype, and could possibly be overcome by using improved models and a larger data set. In any case, the approach is interesting, as a QTL for the maximum level of production reached for one of the curves was located at the same position as a genome-wide significant QTL previously reported for the total egg number. Of course, the results should be further tested on chicken and quail for egg laying. Interestingly, similar work might be carried out on longitudinal data in other species, like the growth curve or the lactation curve in ruminants.
Quail and husbandry
Two lines, LTI and DD, with different origins [19, 20] and selected for early egg production  or for high duration of tonic immobility  were crossed reciprocally (12 single-pair matings) to produce the F1 generation. Ten males and 30 egg-laying females were selected at random across the F1 families to produce the F2 generation. Each F1 sire was mated to three full sisters from another F1 family to obtain 434 F2 female chicks (39 to 45 birds per sire family and 12 to 19 quail per full-sib family), in three consecutive hatches. The 57–week egg laying test analyzed in this work was completed by 359 birds.
The F2 quail were given successively standard starter (2901 kcal ME/kg, 11.5% moisture, 7% ash, 27% total protein, 8% fat and 4% crude fiber) and commercial layer (2709 kcal ME/kg, 11.5% moisture, 12% ash, 20% total protein, 4% fat and 4% crude fiber) diets. At 5 weeks of age they were assigned at random to individual cages of a 4-tier battery maintained at 22°C, in which they remained until the end of the experiment with free access to feed and drinking water.
Egg production was recorded daily and individually from 5 weeks of age onwards. Global egg laying curves for all the F2 (and F0) quail and individual curves fitted to the 19 consecutive 3-week egg production of each F2 quail were obtained under three different nonlinear theoretical models [4, 6, 7] of egg laying. The parameters of the models with a biological meaning related to egg production were chosen as traits for which QTL analyses were carried out. For each model, the phenotypes were the estimates of the parameters of the individual laying curves.
Model 1 was: elr(t) = a exp(-bt)/(1+exp(-c(t-d))), where t is the number of weeks in test, elr is the egg laying rate (%), a is a scale parameter, b is the rate of decrease of egg laying, 1/c is an indicator of the variance in sexual maturity, and d is the mean number of weeks in test until sexual maturity . The traits for which QTL were sought under this model were b, c, and d. QTL analysis was not carried out on a because this parameter was highly correlated to b (r = 0.89) and has no biological interpretation.
Model 2 was: y(t) = (yP/t2)t - r (yP/t2) Ln [(exp(t/r) + exp(t2/r))/(1+exp(t2/r))] + r b4 Ln [(exp(t/r) + exp((t2+P)/r))/(1+ exp((t2+P)/r))], where t is the number of weeks in test, y is the number of eggs laid over the last 3 weeks, r is the duration of transition and was taken equal to 0.3 as in , y P is the level of constant production (that is, the 3-week production at the egg laying plateau), t 2 is the time at transition from rapid increase to constant egg production (that is, the number of weeks in test until the egg laying plateau), b4 is the rate of decline in production (that is, the weekly rate of change in egg laying after the plateau) and P is the persistency of constant production (that is, the duration of the plateau, in weeks) . The traits for which QTL were sought under this model were y P , t2, and b4. QTL analysis was not carried out on P because it was highly correlated to y P (r = 0.71).
Model 3 was: y(t) = [21k1 (1 - exp(-t))/(1+exp(-t))] - [21(k1-k2) (1 - exp(-t))/(1+exp(-(t-c2)))], where t is the number of weeks in test, y is the number of eggs laid over the last 3 weeks, k1 is the proportion of maximum production for the increasing phase of the curve, k2 is the proportion of maximum production for the decreasing phase of the curve, and c2 is the point of inflexion from the upper level of the increasing phase to the lower level of the decreasing phase of the curve (that is, a measure of persistency of egg production) . The traits for which QTL were sought under this model were k1, k2, and c2.
The estimations of the parameters were obtained using the NLIN procedure of the SAS software . Bounds were imposed on estimates of b4 (≤ 0), k1 (≤ 1) and k2 (≥ 0) to insure the coherence between these estimates and the biological meaning of the corresponding parameters of the curve. Under Models 1 and 2, individual egg laying curves could not be estimated for 5 and 7 quail, respectively, because convergence was not reached after 400 iterations, despite trying multiple starting values and using the Marquardt method. This problem was more severe under Model 3, however, for which 51 individual curves could not be obtained.
Genotyping was carried out at Gifu University, and all quail were typed for the 72 microsatellites listed in the Gifu augmented panel . The genotypic data were validated and stored in the MAPGENA database. All 24 F0, 40 F1 and 434 F2 quail in the present experiment belonged to the resource family set up to establish the first microsatellite quail map  built from 72 loci resolved into 13 linkage groups, and to carry out a previous QTL study of the same population . Since then, additional work [22–24] made it possible to assign all 13 linkage groups (total map distance of 576 cM with a 10 cM average spacing between loci) to chromosomes. Consequently, the unassigned linkage groups Q03, Q04, Q08, Q09, Q10 and Q11 of the first map have been renamed respectively, CJA03, CJA13, CJA09, CJA04, CJA20 and CJA10 in the present work. The sex chromosome was not included in the present study.
Standard QTL detection was performed using the line cross method  implemented in the QTL Express software , and it took into account a hatch effect (3 classes). Five thousand permutations  were carried out to set significance levels (p c ) for the most likely chromosome-wide QTL. Genome-wide significance (p g ) for a QTL detected on a given chromosome was obtained from p c by: p g = 1-(1-p c )1/r, where r was the ratio of the length of this chromosome over the total chromosome length spanned by the present study. Chromosome-wide significance was set up at p c = 0.05. Genome-wide significant and suggestive thresholds were set up respectively at p g = 0.05 and 0.20.
Klomp H: The determination of clutch-size in birds: a review. Ardea. 1970, 58: 1-124.
Minvielle F: The future of Japanese quail for research and production. Worlds Poult Sci J. 2004, 60: 500-507. 10.1079/WPS200433.
Albers GAA: Future trends in poultry breeding. Proceedings of the 10th European Poultry Conference: 21–26. Edited by: Organizing committee. 1998, WPSA Israel Branch. Jerusalem, 16-20. June ; Jerusalem
Yang N, Wu C, McMillan I: New mathematical model of poultry egg production. Poult Sci. 1989, 68: 476-481.
Fialho FB, Ledur MC: Segmented polynomial model forestimation of egg production curves in laying hens. Br Poult Sci. 1997, 38: 66-73.
Grossman M, Gossman TN, Koops WJ: A model for persistency of egg production. Poult Sci. 2000, 79: 1715-1724.
Grossman M, Koops WJ: A model for individual egg production in chickens. Poult Sci. 2001, 80: 859-867.
Fairfull RW, Gowe RS: Genetics of egg production inchickens. Poultry Breeding and Genetics. Edited by: CrawfordRD. Amsterdam: Elsevier. 1990, 705-759.
Marks HL: Genetics of egg production in other galliforms. Poultry Breeding and Genetics. Edited by: Edited by Crawford RD. Amsterdam: Elsevier. 1990, 761-770.
Hocking PM: Review of QTL mapping results in chickens. Worlds Poult Sci J. 2005, 61: 215-223. 10.1079/WPS200461.
Schreiweis MA, Hester PY, Settar PME: Identification of quantitative trait loci associated with egg quality, egg production, and body weight in an F2 resource population of chickens. Anim Genet. 2006, 37: 106-112. 10.1111/j.1365-2052.2005.01394.x.
Minvielle F, Kayang BB, Inoue-Murayama M, Miwa M, Vignal A, Gourichon D, Neau A, Monvoisin JL, Ito S: Microsatellite mapping of QTL affecting growth, feed consumption, egg production, tonic immobility and body temperature of Japanese quail. BMC Genomics. 2005, 6: 87-10.1186/1471-2164-6-87.
Minvielle F, Monvoisin JL, Costa J, Frénot A, Maeda Y: Changes in heterosis under within-line selection or reciprocal recurrent selection: an experiment on early egg production in Japanese quail. J Anim Breed Genet. 1999, 116: 363-377. 10.1046/j.1439-0388.1999.00218.x.
Mills AD, Faure JM: Divergent selection for duration of tonic immobility and social reinstatement behavior in Japanese quail (Coturnix coturnix japonica) chicks. J Comp Psychol. 1991, 105: 25-38. 10.1037/0735-7036.105.1.25.
Kayang BB, Inoue-Murayama M, Hoshi T, Matsuo K, Takahashi H, Minezawa M, Mizutani M, Ito S: Microsatellite loci in Japanese quail and cross-species amplification in chicken and guinea fowl. Genet Sel Evol. 2002, 34: 233-253. 10.1051/gse:2002006.
Kayang BB, Vignal A, Inoue-Murayama M, Miwa M, Monvoisin JL, Ito S, Minvielle F: A first-generation microsatellite linkage map of the Japanese quail. Anim Genet. 2004, 35: 195-200. 10.1111/j.1365-2052.2004.01135.x.
Minvielle F, Monvoisin JL, Costa J, Maeda Y: Long-term egg production and heterosis in quail lines after within-line or reciprocal recurrent selection for high early egg production. Br Poult Sci. 2000, 41: 150-157.
Fairfull RW: Heterosis. Poultry Breeding and Genetics. Edited by: Crawford RD. Amsterdam: Elsevier. 1990, 913-933.
Minvielle F, Coville JL, Krupa A, Monvoisin JL, Maeda Y, Okamoto S: Genetic similarity and relationships of DNA fingerprints with performance and with heterosis in Japanese quail lines from two origins and under reciprocal recurrent or within-line selection for early egg production. Genet Sel Evol. 2000, 32: 289-302. 10.1051/gse:2000119.
Minvielle F, Mills AD, Faure JM, Monvoisin JL, Gourichon D: Fearfulness and performance related traits in selected lines of Japanese quail (Coturnix japonica). Poult Sci. 2002, 81: 321-326.
SAS Institute: SAS®User's Guide: Statistics. 1998, Cary, SAS Institute Inc
Miwa M, Inoue-Murayama M, Kayang BB, Vignal A, Minvielle F, Monvoisin JL, Takahashi H, Ito S: Mapping of plumage colour and blood protein loci on the microsatellite linkage map of the Japanese quail. Anim Genet. 2005, 36: 396-400. 10.1111/j.1365-2052.2005.01335.x.
Miwa M, Inoue-Murayama M, Kobayashi N, Kayang BB, Mizutani M, Takahashi H, Ito S: Mapping of panda plumage color locus on the microsatellite linkage map of the Japanese quail. BMC Genet. 2006, 7: 2-10.1186/1471-2156-7-2.
Kayang BB, Fillon V, Inoue-Murayama M, Miwa M, Leroux S, Fève K, Monvoisin JL, Pitel F, Vignoles M, Mouilhayrat C, Beaumont C, Ito S, Minvielle F, Vignal A: Integrated maps in quail (Coturnix japonica) confirm the high degree of synteny conservation with chicken (Gallus gallus) despite 35 million years of divergence. BMC Genomics. 2006, 7: 101-10.1186/1471-2164-7-101.
Haley CS, Knott SA, Elsen JM: Mapping quantitative traitloci in crosses between outbred lines using least squares. Genetics. 1994, 136: 1195-1207.
Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTLExpress: mapping quantitative trait loci in simple and complexpedigrees. Bioinformatics. 2002, 18: 339-340. 10.1093/bioinformatics/18.2.339.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.
The authors wish to thank Jean-Michel Faure (INRA-SRA, Nouzilly) for kindly providing F0 birds of Line LTI. Expert care of experimental F0, F1 and F2 quail, and data collection by Chantal Moussu and the team at the INRA Unité Expérimentale de Génétique Avicole (Nouzilly, France) are gratefully acknowledged.
FM coordinated the study, analyzed the laying curves, and wrote the paper. BBK, SI, MIM and MM designed, organized, and conducted all the microsatellite work. AV led the mapping part. AN was responsible for the data bank, and DG and JLM supervised and carried out the data collection.
About this article
Cite this article
Minvielle, F., Kayang, B.B., Inoue-Murayama, M. et al. Search for QTL affecting the shape of the egg laying curve of the Japanese quail. BMC Genet 7, 26 (2006) doi:10.1186/1471-2156-7-26
- Japanese Quail
- Tonic Immobility
- Individual Curf
- Total Chromosome Length
- Female Quail