- Research article
- Open Access
Genetic analysis of an F2 intercross between two chicken lines divergently selected for body-weight
© Wahlberg et al; licensee BioMed Central Ltd. 2009
- Received: 11 November 2008
- Accepted: 27 May 2009
- Published: 27 May 2009
We have performed Quantitative Trait Loci (QTL) analysis of an F2 intercross between two chicken lines divergently selected for juvenile body-weight. In a previous study 13 identified loci with effects on body-weight, only explained a small proportion of the large variation in the F2 population. Epistatic interaction analysis however, indicated that a network of interacting loci with large effect contributed to the difference in body-weight of the parental lines. This previous analysis was, however, based on a sparse microsatellite linkage map and the limited coverage could have affected the main conclusions. Here we present a revised QTL analysis based on a high-density linkage map that provided a more complete coverage of the chicken genome. Furthermore, we utilized genotype data from ~13,000 SNPs to search the genome for potential selective sweeps that have occurred in the selected lines.
We constructed a linkage map comprising 434 genetic markers, covering 31 chromosomes but leaving seven microchromosomes uncovered. The analysis showed that seven regions harbor QTL that influence growth. The pair-wise interaction analysis identified 15 unique QTL pairs and notable is that nine of those involved interactions with a locus on chromosome 7, forming a network of interacting loci. The analysis of ~13,000 SNPs showed that a substantial proportion of the genetic variation present in the founder population has been lost in either of the two selected lines since ~60% of the SNPs polymorphic among lines showed fixation in one of the lines. With the current marker coverage and QTL map resolution we did not observe clear signs of selective sweeps within QTL intervals.
The results from the QTL analysis using the new improved linkage map are to a large extent in concordance with our previous analysis of this pedigree. The difference in body-weight between the parental chicken lines is caused by many QTL each with a small individual effect. Although the increased chromosomal marker coverage did not lead to the identification of additional QTL, we were able to refine the localization of QTL. The importance of epistatic interaction as a mechanism contributing significantly to the remarkable selection response was further strengthened because additional pairs of interacting loci were detected with the improved map.
- Quantitative Trait Locus
- Quantitative Trait Locus Analysis
- Epistatic Interaction
- Quantitative Trait Locus Region
- Major Quantitative Trait Locus
This study is based on two unique chicken lines that have been developed by bi-directional selection for body-weight at 56 days of age . The two lines originated from the same base population and 40 generations of divergent selection resulted in more than an eight-fold difference in body-weight at the age of selection between the High-Weight (HWS) and the Low-Weight selected (LWS) lines. We have bred a large F2 intercross between these lines as a resource population for QTL analysis [2, 3]. Despite a powerful experimental design (~850 F2) a previous QTL analysis could only explain about 13% of the residual phenotypic variation for body-weight at selection age in the F2 population and at most 50% of the difference between the parental lines . These results indicated that the phenotypic difference between the two chicken lines may be attributed to many QTL each with small individual effects in concordance with the infinitesimal model for polygenic inheritance. We identified 13 QTL (denoted Growth1 – Growth13) and at each locus the allele from the HWS line was associated with increased growth consistent with the line difference .
To further study the genetic architecture causing the phenotypic difference between the HWS and LWS lines Carlborg et al.  analyzed the QTL data using a genetic model testing for pair-wise interaction between loci (epistasis). The epistatic analysis revealed a network of six interacting loci, all interactions involving the Growth9 QTL on chromosome 7. Growth9 was also detected as the only QTL reaching genome-wide significance in the standard QTL analysis of body-weight at 56 days. The growth-promoting effect of the interacting loci was highly dependent on the genotype at Growth9 as HWS homozygosity at this locus was required for the other interacting loci to have a significant effect on growth. Moreover, LWS line homozygosity at either Growth9 or all three of the loci interacting with Growth9 effectively led to the network having no effect on growth. Hence, the significant marginal effect observed in the one-locus QTL analysis for Growth9 and the suggestive evidence for marginal effects of the three other major loci in the network is manifested solely by the interactions between Growth9 and the other loci in the network. The genetic network was estimated to account for about 45% of parental line difference at 56 days of age .
The previous analyses [2–4] were carried out using a rather sparse microsatellite-based linkage map. The map lacked markers on 13 microchromosomes and there were also several regions with poor coverage. Therefore, it was possible that important QTL were missed due to lack of markers. A high-density linkage map would increase the power to detect QTL and enhance the precision in the QTL analysis.
Here we report results from a QTL analysis, based on a comprehensive linkage map, for body-weight traits using our intercross between the HWS and LWS lines. First, we performed a standard genome-wide one-dimensional QTL analysis to fit a one-locus model to detect QTL with main effects. Secondly, we performed a genome-wide pair-wise search to detect epistatic interactions between loci. Finally, we used data on about 13,000 SNPs to compare allele frequency distributions between the HWS and LWS lines to investigate if we could detect any molecular footprints in the form of selective sweeps caused by selection for body-weight.
Chicken lines and intercross
A large F2 intercross pedigree was set-up by crossing individuals from the HWS and LWS lines [2, 3]. The two parental lines originated from a common base population formed by intercrossing seven partially inbred White Plymouth Rock lines of chicken. The selected lines have been developed and maintained at Virginia Polytechnic Institute and State University, Blacksburg, VA. The sex-average mean body-weight, at the age of selection, for the HWS line was at generation forty 1522 g (SD: ± 36), compared to the LWS line with a sex-average mean of 181 g (SD: ± 5.4). Comprehensive information regarding the development, maintenance and characteristics of the selection lines have been previously reviewed by Dunnington and Siegel .
The development of the intercross and a detailed description of the husbandry have previously been described . Briefly, to produce the reciprocal intercross (i.e., HL and LH F2 progeny) we mated 10 HWS line males with 22 LWS line females and eight LWS line males to 19 HWS line females. From the F1 generation four HL males were mated to 37 LH females and, four LH males were mated to 38 HL females to produce an F2 population consisting of 874 chickens. The F1 and F2 animals were housed in the same facilities and fed the same dietary as the parental selection lines. Chickens were provided ad libitum access to water and feed. The F2 generation was produced in a single hatch and all animals were sacrificed at 70 days of age.
Phenotypic traits analyzed in the F2 intercross.
Mean ± SD
75.2 ± 14.9
38 – 144
179.1 ± 56.8
53 – 385
365.5 ± 113.1
76 – 729
621.6 ± 186.9
134 – 1179
943.3 ± 262.1
182 – 1627
47.4 ± 14.7
12 – 116
103.8 ± 47.3
-10 – 241
179.5 ± 68.1
0 – 372
251.7 ± 88.6
-103 – 466
320.7 ± 94.9
-82 – 778
DNA isolation, marker selection and genotyping
DNA was extracted from blood by AGOWA (Berlin, Germany) from chickens included in the intercross pedigree (P, F1 and F2). Methods for SNP selection and the genotyping protocol used to generate the genotypes for ~13,000 SNP markers for 15 individuals from each parental strain have been published elsewhere . A set of 384 highly informative SNPs were selected from the 13,000 SNPs, to enable the construction of a linkage map comprising about one marker every ten cM throughout the chicken genome. SNP genotyping was conducted with the GoldenGate assay (Illumina, CA) at a DNA concentration of 50 ng/ul. Animals in the pedigree were genotyped at the SNP technology platform in Uppsala (Sweden) while the ~13,000 SNP panel was analyzed at CNG in Paris (France).
Construction of linkage maps
We used the CRI-MAP software package to construct linkage maps . We have previously reported a microsatellite-based linkage map constructed using this intercross pedigree . To build the current autosomal linkage map, genotype data from 138 microsatellites were merged with data from 351 SNPs. The dataset was checked for non-Mendelian inheritance errors using the prepare option in CRI-MAP. Markers were then ordered according to the May 2006 genome assembly and the fixed option in CRI-MAP was used to estimate sex-average map distances, in Kosambi cM, between loci. Markers that were not included in the genome assembly were incorporated at the most likely position in the linkage map based on results from a two-point linkage analysis. The order of markers was evaluated using the functions flips and chrompic in CRI-MAP.
Factors affecting phenotypic trait variation in the F2 population were identified using analysis of variance (ANOVA). Sex was included as a fixed effect in the QTL model for all body-weight traits. In contrary to the analysis performed by Park et al.  and Jacobsson et al.  using the same F2 intercross, but consistent with our previous analysis of epistatic interaction  we did not include family (F2 full sibling family) as a fixed effect in our QTL model because the number of full-sib families was large (75 in total) and because maternal effects in chicken are expected to be small. Furthermore, we did not include reciprocal mating in the QTL model because our previous analysis showed that mating type did not have a significant effect on growth .
The QTL analysis of autosomes was performed using interval mapping as first described by Lander and Botstein  and later modified for analyzing crosses between outbred lines by Haley and co-workers [9, 10]. We used QTLExpress  to calculate QTL genotype probabilities conditional on marker genotypes. These were then used in a least square regression-based QTL analysis conducted using an in-house QTL mapping software . Empirical significance thresholds levels were established for each trait by permutation  and the same significance thresholds were used for all traits as the variation between different traits were minimal. The 5% genome-wide significance level was set to an F-value > 9 and as genome-wide suggestive level we used the chromosome-wide significance level for chromosome 4 (F-value > 6), because this chromosome represents approximately 5% of the chicken linkage map. With this suggestive level of significance we expect to find one false positive QTL (type-1 error) per trait and genome scan. Epistatic interaction analysis was conducted as described by Carlborg et al. [4, 12]. In short, a statistical model including sex as a fixed effect, the additive, dominance and all pair-wise epistatic effects (additive-by-additive, additive-by-dominance, dominance-by-additive and dominance-by-dominance) of QTL pairs was used. This epistatic analysis increases the power to identify loci whose effect is dependent on the genotype at another locus. QTL pairs that reached the 5% genome-wide significance threshold in a randomization test for the joint effect of the epistatic pair (no QTL vs. two interacting QTL or one QTL with significant individual effect vs. two interacting QTL) and a 1% significance threshold in a model-selection randomization test for the joint effect of the epistatic parameters (two non-interacting QTL vs. two interacting QTL) are reported as epistatic pairs. The randomization testing procedure was previously described in detail . QTL peaks located within 25 cM of each other were assumed to represent the same locus.
To visualize the effect of individual QTL pairs and also to estimate the total contribution of the epistatic network, we produced a dataset including only F2 individuals where the genotype for all QTL in the main four locus network affecting body-weight at 56 days of age could be determined with high confidence (n = 727). Only individuals where the genotype probabilities exceeded 0.70 for all four loci in the network were selected. The individuals were then classified according to the most probable genotype (i.e. HH, HL or LL) for each locus. We calculated sex-averaged phenotypic means for all nine two-locus genotype classes and plotted the data. The combined effect of the epistatic loci was visualized as follows. First, the dataset was stratified according to the genotype at the main locus Growth9, one group to include the HWS homozygotes, a second group to contain the heterozygotes and a third group including LWS homozygotes individuals (HH; n = 176, HL; n = 365, LL; n = 186). For each of the three groups we then created four more strata with individuals that were HWS homozygous for 0, 1, 2 or 3 of the loci Growth4, Growth6 and Growth12. A one-sided two-sample t-test was used to test for differences in mean body-weight at 56 days of age between the groups.
Fst values  were calculated for each of the segregating SNP markers using the GenePop software (version 4.0)  to estimate the degree of differentiation between lines because the two lines are expected to show marked allele frequency differences at those loci that have responded to the divergent selection. Observed homozygosity was measured as the number of homozygous individuals divided by the total number of individuals. Confidence intervals for QTL regions for the body-weight trait with the highest statistical significance were defined with the one-LOD drop method .
A linkage map including 434 genetic markers
Summary of the linkage map.
Map distance (cM)
Physical coverage (Mb)1
Chr. length (Mb)2
Twenty out of the 138 microsatellites were excluded from the current map either because of low information content or because they showed inheritance errors that became apparent with the high-density marker map. Further, 35 SNPs were excluded, as they did not give reliable genotype calls, and additionally 33 SNPs out of the 384, were not considered in the current study as they are located on the Z chromosome. This updated sex-average autosomal linkage map now includes 434 genetic markers and spans totally 2954 cM, with an average distance between adjacent markers of 6.8 cM. The total map length is considerably shorter than expected from the old consensus map  but in good agreement with a new high-resolution, SNP-based consensus map for chicken . It is clear that the inflated map length in the previous consensus map is because it was primarily based on microsatellite markers which have a higher rate of genotyping errors than Illumina Golden gate SNP assays as used in the present study and the new consensus map.
The current linkage map includes markers on chromosome 1–28, and in addition we also have markers on unassigned linkage group E22C19W28E50C23. Two genetic markers (GCT004, rs16748775) did not show linkage to any other marker in the dataset. It is possible that they are located on microchromosomes lacking other markers, or alternatively the genetic distance to neighboring markers was too far to be detected in the linkage analysis.
By calculating the distance between the first and last marker for each linkage group, we have physical coverage over 890.3 Mb corresponding to approximately 93% of the autosomal chicken genome assembly (May, 2006).
Standard QTL analysis
Descriptive statistics for phenotypic traits analyzed in this study are presented in Table 1. Below we present results from standard genome-wide scans for marginal QTL effects on body-weight traits in the F2 population.
QTL with significant marginal effects on body-weight and growth traits.
a ± SE3
d ± SE3
25.4 ± 4.2
1.8 ± 6.5
59.0 ± 12.2
5.0 ± 8.2
14.5 ± 4.2
5.0 ± 6.2
34.2 ± 9.2
7.8 ± 13.6
2.6 ± 0.8
-2.3 ± 1.3
2.4 ± 0.8
-2.7 ± 1.3
2.5 ± 0.7
-3.0 ± 1.1
2.4 ± 0.7
-3.1 ± 1.1
16.3 ± 3.9
-2.7 ± 5.9
15.9 ± 3.7
1.9 ± 5.7
51.4 ± 11.1
-7.2 ± 16.9
36.3 ± 8.3
-10.6 ± 12.6
8.3 ± 2.1
-3.7 ± 3.1
10.0 ± 2.5
-3.9 ± 3.8
18.3 ± 5.1
-7.1 ± 7.7
12.8 ± 3.5
-1.5 ± 5.1
6.1 ± 2.3
-9.0 ± 3.4
8.2 ± 2.7
-10.4 ± 4.1
12.2 ± 4.0
-10.9 ± 6.0
23.1 ± 5.5
-8.8 ± 8.0
13.1 ± 3.3
-4.7 ± 4.9
43.2 ± 9.0
-12.0 ± 13.0
53.6 ± 11.9
-11.6 ± 17.3
17.7 ± 4.1
0.2 ± 5.8
10.5 ± 2.9
1.3 ± 4.4
22.4 ± 6.1
-4.9 ± 9.8
3.8 ± 0.8
0.9 ± 1.3
3.7 ± 0.8
0.6 ± 1.3
Three QTL regions influenced early growth (0–14 days) whereas four regions appeared to primarily affect later growth (14–70 days). For all QTL the allele originating from the HWS line increased growth. Each individual QTL explained only a small portion of the residual phenotypic variation in the F2 population, ranging from 1.5 – 4.4% (Table 3).
We detected a QTL influencing growth between 0 – 14 days and body-weight at 14 days of age with a QTL peak at 470 cM at chromosome 1. This QTL was not detected in our previous analysis and may represent a novel QTL. However, because it occurs in the vicinity of Growth1 we cannot exclude the possibility that it represents the same locus although Growth1 affects primarily later growth (Table 3).
Detection of epistatic interactions
Detection of selective sweeps using a dense SNP screen
We searched for potential selective sweeps by screening 15 birds from each line using a SNP panel comprising about 13,000 loci evenly distributed throughout the genome. Two birds from the LWS line were excluded from the analysis as results from a phylogenetic analysis indicated that there had been a sample mix up. We then calculated Fst and homozygosity values to identify the regions showing the most pronounced divergence between lines.
Mean Fst and observed homozygosity across the genome and within QTL confidence intervals.
(Mean ± SD)
Fst ± SD
0.33 ± 0.29
0.74 ± 0.24
0.75 ± 0.22
0.38 ± 0.33
0.75 ± 0.24
0.81 ± 0.22
171.9 – 176.0
0.51 ± 0.38
0.76 ± 0.21
0.98 ± 0.05
5.1 – 35.0
0.40 ± 0.35
0.78 ± 0.24
0.82 ± 0.21
21.8 – 33.3
0.45 ± 0.29
0.79 ± 0.24
0.71 ± 0.28
6.0 – 12.0
0.29 ± 0.30
0.67 ± 0.23
0.86 ± 0.16
We calculated Fst-values for each marker located in QTL intervals for body-weight at 56 days of age and compared the mean Fst to the average genome-wide Fst-values. The QTL regions were here defined by their confidence intervals calculated by the one-LOD drop method . The results are presented in Table 4 together with observed homozygosity estimates. Only minor differences were found between the average Fst-values in QTL intervals and the genome average. This suggests that the size of the haplotype blocks associated with a causative QTL mutation are small compared to the size of the QTL confidence intervals obtained in our F2 segregation analysis. However, the confidence interval for Growth1 on chromosome 1 showed a markedly higher divergence between lines and near homozygosity within the LWS line (Table 4). This result suggests that a selective sweep has occurred in the LWS line. The interval showing near fixation within the LWS line is 3.5 Mb in size and spans from 173.4 Mb to 176.9 Mb.
The results of the present study are consistent with the major conclusions obtained in our previous QTL analysis of this large intercross pedigree based on a less complete linkage map [2–4]. Despite the large phenotypic differences between the High-Weight and Low-Weight selected lines we were not able to detect any major QTL with a large marginal effect. Only five loci, Growth1, 4, 6, 9 and 10, reached genome-wide significance for at least one body-weight or growth trait and none explained more than 4.4% of the residual variance for a growth-related trait (Table 3). However, the importance of epistatic interaction as a mechanism contributing significantly to the remarkable response to selection was further strengthened because additional pairs of interacting loci were detected in the present study. We are currently using an Advance Intercross Line (AIL) to replicate and further explore interactions among these loci.
The HWS and LWS lines originated from the same base population that was formed by crossing seven partially inbred lines of White Plymouth Rock broiler chickens . That the results of the QTL analysis infer that no QTL with large individual effects segregated in the base population is consistent with the steady selection response in both directions that have been obtained during the 40 generations of divergent selection [1, 2]. This pattern is also consistent with a rather slow change in allele frequencies at QTL during the course of selection or a gradual release of selectable additive genetic variation from the epistatic loci .
The remarkable response to the divergent selection for juvenile body-weight in the HWS and LWS lines must reflect allele frequency changes between lines at those loci that have responded to selection. The chromosomal regions harboring such loci are expected to show reduced variation within lines and also a higher divergence between lines. Such selective sweeps are caused by hitchhiking of closely linked loci during selection . The size of a selective sweep depends on the number of generations that have passed since the QTL mutation occurred and the local recombination rate in the region. Thus, a recently derived mutation will be associated with a large haplotype block whereas an old mutation that well predates the initiation of the selection experiment is expected to be associated with a smaller haplotype block.
The outcome of the analysis of Fst and homozygosities based on the screening of 13,000 SNPs in the HWS and LWS lines are consistent with the results of the QTL analysis and the observed selection response. We found only minor differences in the average homozygosity as well as average Fst values within QTL intervals compared with the genome average (Table 4). The result suggests that selection has primarily been acting on standing genetic variation that existed well before the initiation of the selection experiment, which is supported by results of a recent simulation study exploring the potential role of epistatic interactions in response to directional selection . The causal mutation will occur in many haplotype combinations if it has been transmitted through a large number of meiotic events before it reaches a high allele frequency in the selected population and as a consequence the selective sweep will be short and not detected by an Fst analysis based on the rather sparse marker set used in this study (~1 polymorphic SNP/200 kb). The results suggest that the release of genetic variance due to epistatic interaction  is a more likely explanation for the long-term selection response in these lines rather than the occurrence of new mutations during the course of selection .
An indication of a selective sweep was, however, observed for Growth1 on chromosome 1 where there was a higher Fst value between lines than the genome average. The confidence interval for Growth1 was smaller than for the other growth QTL which facilitated the detection of footprints of selection. The homozygosity in this QTL interval for the HWS line was close to the genome average whereas the region showed almost complete fixation in the LWS line. The pattern suggests that while this locus has contributed to the selection for low growth in the LWS line it may have been selectively neutral in HWS line. The causal mutation(s) for this QTL is expected to be located within this region of high homozygosity. The region with near homozygosity in the LWS line spans from position 173.4 to 176.9 Mb on chicken chromosome 1, a region still too large and containing too many genes to pin-point strong positional candidate genes.
To be an exhaustive search for QTLs, a genome scan should be based on a marker set that covers the entire genome. The map used in the current study covered 93% of the assembled chicken genome a clear improvement to the ~80% coverage achieved with our previous microsatellite-based map. In particular, we have added markers on six microchromosomes that previously lacked markers. There are, however, still seven microchromosomes lacking markers in the current map. This is because they are missing in the genome assembly and no markers from these chromosomes have yet been reported. These missing microchromosomes are all small, on the order of five Mb or less but are expected to have a higher gene density than the macrochromosomes . It is apparent that parts of the chicken genome, including some microchromosomes, are difficult to clone in bacterial vectors [24, 25] and this is the reason why they are missing from the current genome assembly. Hopefully, the use of new sequencing technologies that do not require vector-based cloning will allow us to fill the holes in the current chicken genome assembly to eventually making a complete genome scan feasible.
If we compare our current QTL analysis based on an improved linkage map with our previous QTL analysis  there are as many as six previously reported loci (Growth2, Growth3, Growth8, Growth10, Growth11 and Growth13) that did not reach statistical significance in the present study. Our previous analysis included family as a fixed effect in the regression model whereas we decided to not include family as a fixed effect in the statistical model used in the present study. To test whether the new information provided by the improved linkage map or the change in regression model caused this discrepancy we reran the analysis including family as a fixed effect for those growth loci that we could not replicate. With a QTL model including Family as a fixed effect, the F-value associated with Growth2 increased from 5.4 to 6.3 and therefore reached the suggestive significance threshold. For Growth3, the F-value increased from 4.5 to 6.0 with family effect in the model. A more dramatic increase in significance was observed for Growth8 since the F-value increased from not reaching even the suggestive threshold to become genome-wide significant with an F-value of 9.3 when the Family effect was included in the model. It is very likely that these three loci (Growth2, 3 and 8) represents true QTL effects and at each locus the allele from the HWS line is increasing growth as expected from the line differences in growth. In contrast, for Growth10, 11 and 13 there was only a minor difference in statistical significance between the two models and in these cases it is more likely that the previous suggestive evidence for QTL at these positions in fact were false positives.
The eight-fold difference in body-weight at selection age between two divergently selected chicken lines is determined by many QTL each with a small individual effect. Although no major QTL explaining a large proportion of the residual phenotypic variance was found, a network of interacting loci with a large effect was identified using the current extended linkage map. Furthermore, results from the QTL study was consistent with our analysis of Fst and homozygosity using a marker set including ~13,000 SNP markers. The results from that analysis failed to detect any major selective sweeps within QTL region suggesting that selection has been acting on standing genetic variation present in the founder population rather than recently derived mutations.
ÖC was supported by grants from the Swedish Foundation for Strategic Research, the Swedish Research Council, the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning and the European Science Foundation (EURYI). LA was supported by grants from the Swedish Foundation for Strategic Research and The Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning. The work at CEA/IG-Centre National de Génotypage was supported by the French Ministry of Research (Ministère délégué à la Recherche). The SNP technology platform in Uppsala was supported by the Knut and Alice Wallenberg Foundation (via Wallenberg Consortium North). We thank Tomas Axelsson and Kristina Larsson for assistance with genotyping in Uppsala.
- Dunnington EA, Siegel PB: Long-term divergent selection for eight-week body weight in white Plymouth rock chickens. Poult Sci. 1996, 75 (10): 1168-1179.View ArticlePubMedGoogle Scholar
- Jacobsson L, Park HB, Wahlberg P, Fredriksson R, Perez-Enciso 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 (2): 115-125. 10.1017/S0016672305007767.View ArticlePubMedGoogle Scholar
- Park HB, 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 Genomics. 2006, 25 (2): 216-223. 10.1152/physiolgenomics.00113.2005.View ArticlePubMedGoogle Scholar
- Carlborg O, Jacobsson L, Ahgren P, Siegel P, Andersson L: Epistasis and the release of genetic variation during long-term selection. Nat Genet. 2006, 38 (4): 418-420. 10.1038/ng1761.View ArticlePubMedGoogle Scholar
- Wahlberg P, Stromstedt L, Tordoir X, Foglio M, Heath S, Lechner D, Hellstrom AR, Tixier-Boichard M, Lathrop M, Gut IG, et al: A high-resolution linkage map for the Z chromosome in chicken reveals hot spots for recombination. Cytogenet Genome Res. 2007, 117 (1–4): 22-29. 10.1159/000103161.View ArticlePubMedGoogle Scholar
- Green P, Falls K, Crooks S: Documentation for CRI-MAP, version 2.4. 1990, Washington School of Medicine, St Louis, MOGoogle Scholar
- Jacobsson L, Park HB, Wahlberg P, Jiang S, Siegel PB, Andersson L: Assignment of fourteen microsatellite markers to the chicken linkage map. Poult Sci. 2004, 83 (11): 1825-1831.View ArticlePubMedGoogle Scholar
- Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121 (1): 185-199.PubMed CentralPubMedGoogle Scholar
- Haley CS, Knott SA: A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity. 1992, 69 (4): 315-324.View ArticlePubMedGoogle Scholar
- Haley CS, Knott SA, Elsen JM: Mapping quantitative trait loci in crosses between outbred lines using least squares. Genetics. 1994, 136 (3): 1195-1207.PubMed CentralPubMedGoogle Scholar
- Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. Bioinformatics. 2002, 18 (2): 339-340. 10.1093/bioinformatics/18.2.339.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.PubMed CentralPubMedGoogle Scholar
- Carlborg O, Andersson L: Use of randomization testing to detect multiple epistatic QTLs. Genet Res. 2002, 79 (2): 175-184. 10.1017/S001667230200558X.View ArticlePubMedGoogle Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, Columbia University Press, New YorkGoogle Scholar
- Raymond M, Rousset F: GENEPOP Version 1.2: population genetics software for exact tests and ecumenicism. J Hered. 1995, 86 (3): 248-249.Google Scholar
- Paterson AH, Damon S, Hewitt JD, Zamir D, Rabinowitch HD, Lincoln SE, Lander ES, Tanksley SD: Mendelian factors underlying quantitative traits in tomato: comparison across species, generations, and environments. Genetics. 1991, 127 (1): 181-197.PubMed CentralPubMedGoogle Scholar
- Groenen MA, Cheng HH, Bumstead N, Benkel BF, Briles WE, Burke T, Burt DW, Crittenden LB, Dodgson J, Hillel J, et al: A consensus linkage map of the chicken genome. Genome Res. 2000, 10 (1): 137-147.PubMed CentralPubMedGoogle Scholar
- Groenen MA, Wahlberg P, Foglio M, Cheng HH, Megens HJ, Crooijmans RP, Besnier F, Lathrop M, Muir WM, Wong GK, et al: A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate. Genome Res. 2009, 19 (3): 510-519. 10.1101/gr.086538.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlborg Ö, Andersson L, Kinghorn B: The use of a genetic algorithm for simultaneous mapping of multiple interacting quantitative trait loci. Genetics. 2000, 155 (4): 2003-2010.PubMed CentralPubMedGoogle Scholar
- Le Rouzic A, Siegel PB, Carlborg Ö: Phenotypic evolution from genetic polymorphisms in a radial network architecture. BMC Biol. 2007, 5: 50-10.1186/1741-7007-5-50.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith JM, Haigh J: The hitch-hiking effect of a favourable gene. Genet Res. 1974, 23 (1): 23-35.View ArticlePubMedGoogle Scholar
- Keightley PD, Hill WG: Quantitative genetic variation in body size of mice from new mutations. Genetics. 1992, 131 (3): 693-700.PubMed CentralPubMedGoogle Scholar
- International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432 (7018): 695-716. 10.1038/nature03154.View ArticleGoogle Scholar
- Gordon L, Yang S, Tran-Gyamfi M, Baggott D, Christensen M, Hamilton A, Crooijmans R, Groenen M, Lucas S, Ovcharenko I, et al: Comparative analysis of chicken chromosome 28 provides new clues to the evolutionary fragility of gene-rich vertebrate regions. Genome Res. 2007, 17 (11): 1603-1613. 10.1101/gr.6775107.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.