- Open Access
Linkage disequilibrium in crossbred and pure line chickens
© Fu et al.; licensee BioMed Central. 2015
- Received: 14 June 2014
- Accepted: 26 January 2015
- Published: 26 February 2015
Both genome-wide association (GWA) studies and genomic selection depend on the level of non-random association of alleles at different loci, i.e. linkage disequilibrium (LD), across the genome. Therefore, characterizing LD is of fundamental importance to implement both approaches. In this study, using a 60K single nucleotide polymorphism (SNP) panel, we estimated LD and haplotype structure in crossbred broiler chickens and their component pure lines (one male and two female lines) and calculated the consistency of LD between these populations.
The average level of LD (measured by r 2) between adjacent SNPs across the chicken autosomes studied here ranged from 0.34 to 0.40 in the pure lines but was only 0.24 in the crossbred populations, with 28.4% of adjacent SNP pairs having an r 2 higher than 0.3. Compared with the pure lines, the crossbred populations consistently showed a lower level of LD, smaller haploblock sizes and lower haplotype homozygosity on macro-, intermediate and micro-chromosomes. Furthermore, correlations of LD between markers at short distances (0 to 10 kb) were high between crossbred and pure lines (0.83 to 0.94).
Our results suggest that using crossbred populations instead of pure lines can be advantageous for high-resolution QTL (quantitative trait loci) mapping in GWA studies and to achieve good persistence of accuracy of genomic breeding values over generations in genomic selection. These results also provide useful information for the design and implementation of GWA studies and genomic selection using crossbred populations.
- Linkage Disequilibrium
- Genomic Selection
- Pure Line
- Single Nucleotide Polymorphism Array
- Female Line
Genome-wide association (GWA) studies and genomic selection are currently widely applied in animal genetics research and animal breeding programs. Using thousands of genetic markers, the ultimate goal of most GWA studies is to find causal polymorphisms that affect a phenotype , whereas a reliable prediction of the total genetic value of selection candidates is the major goal of genomic selection . Despite their different aims, the success of these two approaches depends primarily on the level of linkage disequilibrium (LD) between markers and causal polymorphisms [3,4]. For example, long-distance LD is useful in GWA studies when using a relatively low-density set of markers [5-7]. However, mapping resolution is expected to be lower when LD extends over long distances because multiple markers across a wide chromosomal region may be in LD with the causal polymorphism and all show significant associations with the trait. Conversely, for the opposite reason, a low level and extent of LD can be useful for high-resolution association mapping.
Progress in next-generation sequencing and high-density single nucleotide polymorphism (SNP) genotyping technologies offer unprecedented opportunities for detecting causal polymorphisms or achieving high accuracies of prediction in genomic selection [8-12]. However, taking full advantage of these new technologies may be limited in livestock populations due to the high level and extent of LD. Although LD extends over long distances in most livestock populations, comparisons of LD patterns between populations show that shared haplotype segments are much shorter when the population consists of multiple purebred populations [13-15]; this indicates that LD decays more quickly in multi-breed or crossbred populations than in purebred populations. Therefore, in cases where LD does not extend over long distances, multi-breed and crossbred populations can be potentially useful for fine mapping and identification of causal polymorphisms. Furthermore, as explained in the Discussion section, using a crossbred population as reference population in genomic selection can also be advantageous, particularly in livestock species with crossbreeding programs, such as poultry, pigs and beef cattle.
In this study, our aim was to characterize the consistency of LD and differences in LD between crossbred and their component purebred populations. As for previous studies on LD in layer  and village chickens , we used the Illumina 60K chicken SNP panel  which contains over 10 times more SNPs than that of most other studies on LD in chickens [19-23].
All 2341 genotyped crossbred chickens were produced by a three-way cross of B × [C × D], in which males of line B were mated with CD crossbred females, which is a crossbred female product produced by a two-way C × D cross, in which males of line C were mated with females of line D. Crossbred individuals were sampled from two broiler chicken populations: (i) broiler chickens from the field (end product meat birds), which will be referred to as field crossbred chickens (or fBCD, n = 1093) and (ii) broiler chickens from a pedigree house (produced for genetic evaluation of the pedigree B1 sires), which will be referred to as pedigree crossbred chickens (or pBCD, n = 1248) (Figure 1).
To assess the extent to which the LD pattern in crossbred populations can be predicted based on the genotypes of their component pure lines, we created a combined population by simply combining the genotype data of a random sample of 200 chickens of line B, 100 of line C and 100 of line D. The 2:1:1 proportion of lines B:C:D was used to mimic the expected genetic contribution of these lines to the autosomes of crossbred individuals. This combined population will be referred to as the combined BCD (or cBCD) population.
The Illumina 60K SNP chip contains 57636 SNPs . In this study, we used only SNPs with assigned positions on autosomes (based on the latest reference genome, Gallus gallus 4.0 UCSC, May 2012). Within each pure line and crossbred population, we discarded SNPs with a call rate less than 90%, Mendelian inconsistency greater than 0.001 and minor allele frequency (MAF) less than 0.05. Also, SNPs that strongly deviated from Hardy–Weinberg equilibrium (p value < 0.001) in the pure lines were discarded, as well as SNPs on chromosome 16 and two linkage groups because there were too few SNPs in the 60K SNP panel for these chromosomes.
Haplotype and haploblock analyses
Determining haplotype phase and frequency is necessary to estimate LD and can provide useful information about breed-specific haplotypes and the history of artificial selection. We used BEAGLE (Version 3.3.2)  to infer haplotype phase for each genotyped individual in each population. As in Badke et al. , we set BEAGLE to run 100 iterations of the phasing algorithm and to sample 100 haplotype pairs for each individual per iteration.
In theory, haplotype homozygosity is defined as the likelihood of randomly sampling two identical haplotypes from a population, which is calculated as the sum of squares of haplotype frequencies . Based on the results of the haplotype phases obtained with BEAGLE, haplotype homozygosity was estimated using haplotype frequencies for 250-kb sliding windows, with a step size of 25 kb. For each population, Haploview (Version 4.2)  was used to define haplotype blocks (haploblocks) with the built-in algorithm suggested by Gabriel et al. . In this model, the confidence interval of observed values of LD measured by D' was estimated to determine the upper and lower confidence bounds of D’ (5% tails of the overall probability distribution of D'), and the blocking structure was determined by defining SNP pairs to be in “strong LD” if the upper confidence bound was above 0.98 and the lower bound was above 0.7 .
Estimation of linkage disequilibrium
where r ij is the correlation between alleles at SNP i (alleles M and m) and alleles at SNP j (alleles N and n); f(MN) is the observed frequency of haplotype MN, which can be simply obtained from the phasing results; and f(M), f(m), f(N) and f(n) are observed frequencies of alleles M, m, N and n, respectively .
Previous studies on LD in chickens showed that the extent of LD over physical distances varies greatly among the different categories of avian chromosomes: macro-chromosomes (GGA1 to 5, GGA for Gallus gallus chromosome), intermediate chromosomes (GGA6 to 10) and micro-chromosomes (GGA11 to 38) . Thus, we estimated LD separately for each category of chromosomes within each population. To visualize the LD pattern for each category of chromosomes in different populations, r 2 values were ordered in ascending order based on the physical distance between the corresponding SNP pairs, and then a rolling average of LD was calculated as the arithmetic mean of all r 2 values for SNP pairs in 25-kb intervals and plotted against physical distance between SNPs.
Estimation of consistency of LD
Consistency of LD between two populations was calculated as the correlation of r between SNP pairs. We used the SNPs that were common to the populations to estimate the consistency of LD as the correlation r ij between the same pairs of SNPs within a given interval in two populations. For simplicity, this will be referred to as the correlation of r. To visualize and compare the correlation of r among different pairs of populations, the pairwise correlation of r was estimated separately for each category of chromosomes in 50-kb non-overlapping intervals and plotted against physical distance.
Quality control criteria and number of SNPs discarded in each population
Chromosomes not included
SNPs not called2
SNPs with a call rate < 0.9
SNPs with a MAF < 0.05
HWE (p value < 0.001)
SNPs in use
The mean r 2 between adjacent SNPs across all autosomes studied here ranged from 0.34 to 0.40 in the pure lines but was on average only equal to 0.24 in the crossbred populations. Due to the different densities of SNPs on each chromosome in the 60K SNP chip, mean r 2 values were similar in the three categories of chromosomes. Furthermore, in the pure lines, at least 53.2% of adjacent SNP pairs had an r 2 greater than 0.2 and 42.2% had an r 2 greater than 0.3, but in the crossbred populations, only 39.5% and 28.4% showed an r 2 greater than 0.2 and 0.3, respectively.
Haploblock and haplotype homozygosity
Summary statistics of haploblock structure across different populations
Median block size (kb)
Maximum block size (kb)
Genome coverage (Mb)2
Mean ± SD nBSNPs5
3.5 ± 3.4
4.2 ± 4.6
3.5 ± 3.9
2.6 ± 2.1
Consistency of LD
The correlation of r measures the degree of agreement of the direction and level of LD for SNP pairs between two populations [8,9]. In principle, these correlations can range from −1 to 1: a high positive value indicates high LD and the same haplotype phase in the two populations, and a high negative value results from a high LD between two markers but with the opposite phase in the two populations .
It should be noted that in this study correlations of r were high and positive (>0.99) between field and pedigree crossbred populations and also between either of these real crossbred populations and the combined BCD population at distances between markers less than 50 kb (see Figure 4). The correlation gradually decayed as the distance between markers increased but still remained high (0.87 to 0.93), even if SNPs were about 5 Mb apart.
Our results show that correlations of r between the pure line C and field crossbred population (0.86 to 0.92; < 10 kb) and between the pure line D and field crossbred population (0.83 to 0.88; < 10 kb) were generally similar but lower than those between the pure line B and field crossbred population (0.91 to 0.94; < 10 kb). Among the pure lines, the correlation of r between the two female lines was slightly higher than that between the male line and either of the female lines.
As with the decay of LD, correlations of r decreased with physical distance and this decrease was greater for the smaller chromosomes. For example, correlations of r between the pure line B and field crossbred population were 0.94, 0.94 and 0.91 in the interval of 0 to 10 kb on macro-, intermediate and micro-chromosomes, respectively, but decreased to 0.39, 0.21 and 0.13, respectively, at an interval of approximately 5 Mb.
Studies of LD in farm animals have been mostly limited to purebred populations and there is limited information about the extent of LD in crossbred populations and the consistency of LD phase between crossbred populations and their parental pure lines . In this study, we characterized the consistency of phase and level of LD between crossbred broiler chickens and their parental pure lines. The crossbred chickens in our study were a three-way cross of B × [C × D] that was produced using three pure lines B, C and D, which contributed 50, 25, and 25% of genetic material, respectively, to the autosomes of the crossbred animals. Our analyses used SNP genotypes on 27 chicken autosomes (GGA1 to 28, excluding GGA16 and other micro-chromosomes because the 60K SNP chip did not include enough markers on these chromosomes). To our knowledge, there is no published report that characterizes LD in crossbred chicken populations and that compares the consistency of phase and level of LD between crossbred and their parental pure lines.
Rapid decay of LD in crossbred populations and micro-chromosomes
The level and pattern of LD differed between the populations used in our study. The two field and pedigree crossbred populations had very similar levels of LD for all physical distances. However, compared with the pure lines, crossbred populations showed a small extent and rapid decay of LD by distance for all three chromosome categories. For example, for the macro-chromosomes, the mean r 2 for SNPs that were 0 to 10 kb apart was 0.32 in the field crossbred population but greater than 0.44 in the pure lines, and the mean r 2 dropped to less than 0.2 at distances between SNPs of ~50 kb in the field crossbred population, whereas in the pure lines, this drop was observed for much greater distances (~225 kb). Similarly, the extent of LD was smaller in crossbred beef cattle than in purebred Angus and Charolais cattle . The rapid decay of LD by distance in crossbred populations can be useful for high-resolution mapping of causal polymorphisms. Indeed, if the extent of LD is small, it is less likely that SNPs far away from a causal polymorphism will be in LD with the polymorphism, which confines associations to SNPs that are in close vicinity of the causal polymorphism, thereby increasing map resolution. Moreover, when using a higher SNP density, the small extent of LD in crossbred populations may be an advantage for genomic selection because the tight linkage between causal polymorphisms and adjacent SNPs is less likely to be broken down by recombination, and therefore the accuracy of genomic predictions will persist over more generations .
An average LD (r 2) greater than 0.2 [3,33,34] or 0.3 [11,35] between adjacent SNPs has been recommended to detect SNPs associated with causal polymorphisms or to achieve a reasonable accuracy of prediction in genomic selection. Although in the pure lines at least 53.2% of adjacent SNP pairs of the 60K SNP panel had an r 2 greater than 0.2 and 42.2% had an r 2 greater than 0.3, in the field crossbred population, only 28.4% of adjacent SNP pairs showed an r 2 greater than 0.3, which suggests that a higher density SNP panel would be an advantage for GWA studies or to implement genomic selection in commercial crossbred populations.
In our study, the extent of LD varied greatly between chromosome categories in different populations and decreased as the distance between SNPs increased (Figure 3). Consistent with previous studies [16,20,23], our results showed that small-size autosomes had less LD than large-size autosomes. The differences in LD between small- and large-size autosomes have been attributed mainly to differences in recombination rates per unit of physical distance, with micro-chromosomes showing the largest recombination rate per Mb (6.4 cM/Mb), followed by the intermediate chromosomes (3.9 cM/Mb), and then macro-chromosomes (2.8 cM/Mb) .
Small haploblock size in crossbred populations
Consistent with the results of the LD analysis, results of the haploblock analyses showed that the crossbred populations had the lowest genomic coverage (229.0 Mb) in haploblocks and the smallest average number of SNPs forming haploblocks. Given the size of the autosomes (~918.9 Mb; GGA1 to 28 without GGA16; Gallus gallus 4.0, November 2011), this means that haploblocks covered only 24.9% of the genome of the crossbred populations, which is nearly half that of the pure lines. In a study using commercial chickens, the genome coverage in haploblocks ranged from 337.1 to 599.4 Mb, with an F1 cross between two layer populations showing a lower genome coverage (337.1 Mb) than the layer pure line . Although the pure lines in this study were not the parental lines of this F1 cross, the low genome coverage in haploblocks in this two-way layer cross is consistent with what we observed in the three-way broiler cross.
The percentage of SNPs forming haploblocks differed between populations (Table 2). Only 31.4% of SNPs formed haploblocks in the field crossbred population, compared to more than 60% in the three pure lines. This finding shows that, because of the small extent of LD, most markers did not form haploblocks in the crossbred populations. The small percentage of SNPs forming haploblocks and the small haploblock size in the crossbred populations also indicate that the 60K SNP chip used in our study does not have an adequate SNP density for high-resolution characterization of the haplotype structure in crossbred chickens. Megens et al.  investigated the LD and haplotype diversity on four ~1-cM regions on macro- and micro-chromosomes and suggested that whole-genome marker assays would need to contain at least 100 K informative SNPs to exploit haplotype information in commercial chicken populations. Consistent with this, a study on the haplotype structure of traditional and village chickens suggested using at least 90 to 110 K SNPs to construct a whole-genome haplotype map for these populations . Therefore, the recently available 600K SNP chicken panel  is likely to provide a higher resolution haplotype map in crossbred chickens.
Haplotype homozygosity is a measure of haplotype diversity in a population. In a previous study, a relatively small number (1 to 7) of haplotypes accounted for most of the haplotype diversity (>90%) found on the macro-chromosomes, but not on the micro-chromosomes . These results are consistent with our findings across all populations. Within 250-kb windows, the pure-line chickens showed limited haplotype sharing, which is consistent with a moderate correlation of r (<0.5) at the same distance between these populations (Figure 4). In some chromosomal regions, high levels of HH extended over longer distances in one pure line but not in the others. Because one of the key characteristics of positive selection is the presence of unusually long-range HH , these differences in HH patterns could be considered as evidence of recent positive selections in these pure lines. Therefore, it would be interesting to conduct an analysis of selection sweeps in these pure lines.
Consistency of LD from pure lines to crossbreds
Because of the pyramidal structure of conventional chicken breeding programs, there is approximately four years of lag time from pedigree pure line birds to end-product crossbred birds . To assess if LD persists between pedigree pure lines and commercial crossbred birds, we estimated the correlation r between pure lines in the pedigree program (top of the pyramid) and crossbred chickens that were sampled from the field (bottom of the pyramid). Correlations of r were high (0.83 to 0.94) between these populations for closely spaced SNPs (0 to 10 kb) but these correlations decreased as the distance between SNPs increased; correlations dropped by 4 to 9% from 0 to 10 kb to 10 to 50 kb distances between markers across the three chromosome categories. In our study, the 60K SNP panel provided an average marker spacing of one informative SNP per ~25 kb; therefore, it is expected that using a higher density SNP panel (such as the recently available 600K SNP array) will increase the accuracy of genomic selection of pure lines for crossbred performance. This conclusion is consistent with results from two simulation studies [32,41] that showed that training in crossbred populations led to slightly less accurate predictions of true breeding values of selection candidates in purebred populations compared with training only in the selected purebred population. However, by increasing the density of the SNP panel, differences in accuracies became much smaller.
The correlation of r between crossbred populations and pure lines differed in the three chromosome categories, with macro-chromosomes showing the highest levels of correlation and micro-chromosomes showing the lowest. These results indicated that, in GWA studies or in genomic selection programs, micro- and intermediate chromosomes would require a higher SNP density per kb than macro-chromosomes. In our study, each pure line showed a different level of correlation of r with the crossbred populations and, of all comparisons, this value was the highest between line B and the crossbred populations. As mentioned before, line B was the terminal male line for crossbred chickens B × [C × D], which means that this line contributed 50% of genetic material to the autosomes of crossbred individuals; whereas the female lines C and D were expected to contribute each 25%. As for female lines, the correlation of r was greater between line D and the crossbred populations than between line C and the crossbred populations, which could be due to a greater correlation of r between lines D and B than between lines C and B.
Consistency of LD phase among pure lines
The correlation of r can reflect the relative degree of similarities and divergences between purebred animals. In a large-scale genome-wide survey of SNP variation in cattle breeds, the correlation of r declined as the divergence between breeds increased . In our study, the correlation of r between female lines was slightly but consistently higher than that between female and male lines, which suggests that the two female lines may share a more similar genetic background than each of these lines does with the male line. This is consistent with the fact that male and female lines originated from different breeds, i.e. the male line from Cornish, a meat type breed, and the female lines from dual-purpose breeds. Furthermore, since selection goals of male and female lines are different, this may have contributed to the similarity between female lines and to the divergence between the female lines and the male line. Unlike the male line, which was selected primarily for growth-related traits, the female lines were selected for both reproductive and growth traits.
Using markers on chromosomes 1 and 4, Andreescu et al.  estimated correlations of LD within 500 kb among nine purebred chicken lines from a commercial broiler breeding company; the correlations over all pairs of lines ranged from 0.21 to 0.94, with an average correlation of 0.52. Badke et al.  reported a correlation of LD for distances between markers less than 10 kb that was equal to 0.92 between Landrace and Yorkshire breeds and 0.87 between these breeds and the Duroc breed; these values decreased to 0.41 to 0.57 for distances between markers around 1 Mb. Moreover, a study in cattle found that the correlation of LD for distances between markers less than 10 kb was 0.97 between Dutch black-and-white Holstein-Friesian vs. Dutch red-and-white Holstein-Friesian and New Zealand Friesian vs. Zealand Jersey . In our study, none of the chromosome categories reached this high level of consistency of LD between pure lines for distances between markers less than 10 kb, and correlations of r ranged from 0.73 to 0.82 for distances between markers less than 10 kb, even on macro-chromosomes. It is likely that the genetic diversity of the chicken lines in our study was greater (average Fst > 0.20) than that of both the cattle and pig breeds used in the aforementioned studies. Another possible explanation for this difference is the overall higher recombination rates per unit of physical distance on chicken chromosomes compared with the average ratio of 1 cM/1 Mb in mammalian livestock animals.
Consistency of LD between crossbred populations and the combined BCD population
To assess the extent to which the LD pattern in crossbred populations can be predicted using the genotypic information of their component pure lines, we created a combined BCD population (see Methods section) and studied the differences in LD between this hypothetical population and the actual crossbred populations. Across all three chromosome categories, levels of LD were almost the same in these crossbred and the combined BCD populations (Figure 3). This was also reflected by the consistency of LD between these populations, since the correlations of r were very high (>0.99) for distances between markers less than 50 kb between the crossbred and the combined BCD populations, and decayed gradually as the distance between markers increased, but still remained high (0.87 to 0.93) for markers that were about 5 Mb apart (see Figure 4). These results indicate that by using only genotype information on the pure lines, one can predict the LD in crossbred populations with very high accuracy, as well as the correlation r between crossbred populations and their component pure lines.
In conclusion, our study characterized the extent and consistency of LD in commercial broiler populations from different angles and showed that, between crossbred populations and their component pure lines, the consistency of the level and phase of LD for short distances between markers (0 to 10 kb) is remarkably high. Compared with the pure lines, the crossbred populations showed a considerably lower level of LD and a smaller haploblock size, which suggests that using crossbred animals as a reference population can be an advantage for high-resolution mapping of causal polymorphisms in GWA studies and to achieve better persistence of the accuracy of genomic estimated breeding values over generations in genomic selection programs. However, our results also suggest that a higher SNP density, particularly on micro-chromosomes, is necessary to take full advantage of crossbred populations in GWA studies or in genomic selection programs. Finally, our results prove that LD for short and long distances between markers and haplotype phase for short distances between markers in a crossbred population can be predicted with very high accuracy using genotype information of its parental pure lines.
We would like to thank Heritage Breeders for in-kind contribution (60K SNP genotype data) and financial support; and USDA Chicken GWMAS Consortium, Cobb Vantress, and Hendrix Genetics for access to the developed 60K SNP Illumina iSelect chicken array.
- Risch N, Merikangas K. The future of genetic studies of complex human diseases. Science. 1996;273:1516–7.View ArticlePubMedGoogle Scholar
- Goddard ME, Hayes BJ. Genomic selection. J Anim Breed Genet. 2007;124:323–30.View ArticlePubMedGoogle Scholar
- Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157:1819–29.PubMed CentralPubMedGoogle Scholar
- Zhao H, Nettleton D, Dekkers JCM. Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between single nucleotide polymorphisms. Genet Res. 2007;89:1–6.View ArticlePubMedGoogle Scholar
- Khatkar MS, Thomson PC, Tammen I, Raadsma HW. Quantitative trait loci mapping in dairy cattle: review and meta-analysis. Genet Sel Evol. 2004;36:163–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Rothschild MF, Hu ZL, Jiang Z. Advances in QTL mapping in pigs. Int J Biol Sci. 2007;3:192–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Abasht B, Dekkers JCM, Lamont SJ. Review of quantitative trait loci identified in the chicken. Poult Sci. 2006;85:2079–96.View ArticlePubMedGoogle Scholar
- Badke YM, Bates RO, Ernst CW, Schwab C, Steibel JP. Estimation of linkage disequilibrium in four US pig breeds. BMC Genomics. 2012;13:24.View ArticlePubMed CentralPubMedGoogle Scholar
- De Roos APW, Hayes BJ, Spelman RJ, Goddard ME. Linkage disequilibrium and persistence of phase in Holstein-Friesian, Jersey and Angus cattle. Genetics. 2008;179:1503–12.View ArticlePubMed CentralPubMedGoogle Scholar
- Lu D, Sargolzaei M, Kelly M, Li C, Vander Voort G, Wang Z, et al. Linkage disequilibrium in Angus, Charolais, and Crossbred beef cattle. Front Genet. 2012;3:152.View ArticlePubMed CentralPubMedGoogle Scholar
- Sargolzaei M, Schenkel FS, Jansen GB, Schaeffer LR. Extent of linkage disequilibrium in Holstein cattle in north America. J Dairy Sci. 2008;91:2106–17.View ArticlePubMedGoogle Scholar
- Uimari P, Tapio M. Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds. J Anim Sci. 2011;89:609–14.View ArticlePubMedGoogle Scholar
- Karlsson EK, Baranowska I, Wade CM, Salmon Hillbertz NHC, Zody MC, Anderson N, et al. Efficient mapping of mendelian traits in dogs through genome-wide association. Nat Genet. 2007;39:1321–8.View ArticlePubMedGoogle Scholar
- Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005;438:803–19.View ArticlePubMedGoogle Scholar
- Setoguchi K, Furuta M, Hirano T, Nagao T, Watanabe T, Sugimoto Y, et al. Cross-breed comparisons identified a critical 591-kb region for bovine carcass weight QTL (CW-2) on chromosome 6 and the Ile-442-Met substitution in NCAPG as a positional candidate. BMC Genet. 2009;10:43.View ArticlePubMed CentralPubMedGoogle Scholar
- Qanbari S, Hansen M, Weigend S, Preisinger R, Simianer H. Linkage disequilibrium reveals different demographic history in egg laying chickens. BMC Genet. 2010;11:103.View ArticlePubMed CentralPubMedGoogle Scholar
- Wragg D, Mwacharo JM, Alcalde JA, Hocking PM, Hanotte O. Analysis of genome-wide structure, diversity and fine mapping of Mendelian traits in traditional and village chickens. Heredity (Edinb). 2012;109:6–18.View ArticlePubMed CentralGoogle Scholar
- Groenen MAM, Megens H-J, Zare Y, Warren WC, Hillier LW, Crooijmans RPMA, et al. The development and characterization of a 60K SNP chip for chicken. BMC Genomics. 2011;12:274.View ArticlePubMed CentralPubMedGoogle Scholar
- Heifetz EM, Fulton JE, O’Sullivan N, Zhao H, Dekkers JCM, Soller M. Extent and consistency across generations of linkage disequilibrium in commercial layer chicken breeding populations. Genetics. 2005;171:1173–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Aerts J, Megens HJ, Veenendaal T, Ovcharenko I, Crooijmans R, Gordon L, et al. Extent of linkage disequilibrium in chicken. Cytogenet Genome Res. 2007;117:338–45.View ArticlePubMedGoogle Scholar
- Andreescu C, Avendano S, Brown SR, Hassen A, Lamont SJ, Dekkers JCM. Linkage disequilibrium in related breeding lines of chickens. Genetics. 2007;177:2161–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Abasht B, Sandford E, Arango J, Settar P, Fulton JE, O’Sullivan NP, et al. Extent and consistency of linkage disequilibrium and identification of DNA markers for production and egg quality traits in commercial layer chicken populations. BMC Genomics. 2009;10:S2.View ArticlePubMed CentralPubMedGoogle Scholar
- Megens H-J, Crooijmans RPMA, Bastiaansen JWM, Kerstens HHD, Coster A, Jalving R, et al. Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken. BMC Genet. 2009;10:86.View ArticlePubMed CentralPubMedGoogle Scholar
- Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009;84:210–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Sabatti C, Risch N. Homozygosity and linkage disequilibrium. Genetics. 2002;160:1707–19.PubMed CentralPubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.View ArticlePubMedGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–9.View ArticlePubMedGoogle Scholar
- Devlin B, Risch N. A comparison of linkage disequilibrium measures for fine-scale mapping. Genomics. 1995;29:311–22.View ArticlePubMedGoogle Scholar
- International Chicken Genome Sequencing Consortium. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004;432:695–716.View ArticleGoogle Scholar
- Wright S. The genetical structure of populations. Ann Eugen. 1951;15:323–54.View ArticlePubMedGoogle Scholar
- Hartl DL, Clark AG. Principles of Population Genetics, vol. 116. Sunderland: Sinauer associates; 1997.Google Scholar
- Toosi A, Fernando RL, Dekkers JCM. Genomic selection in admixed and crossbred populations. J Anim Sci. 2010;88:32–46.View ArticlePubMedGoogle Scholar
- Gray MM, Granka JM, Bustamante CD, Sutter NB, Boyko AR, Zhu L, et al. Linkage disequilibrium and demographic history of wild and domestic canids. Genetics. 2009;181:1493–505.View ArticlePubMed CentralPubMedGoogle Scholar
- Shea J, Agarwala V, Philippakis AA, Maguire J, Banks E, Depristo M, et al. Comparing strategies to fine-map the association of common SNPs at chromosome 9p21 with type 2 diabetes and myocardial infarction. Nat Genet. 2011;43:801–5.View ArticlePubMed CentralPubMedGoogle Scholar
- Ardlie KG, Kruglyak L, Seielstad M. Patterns of linkage disequilibrium in the human genome. Nat Rev Genet. 2002;3:299–309.View ArticlePubMedGoogle Scholar
- Nielsen R, Signorovitch J. Correcting for ascertainment biases when analyzing SNP data: applications to the estimation of linkage disequilibrium. Theor Popul Biol. 2003;63:245–55.View ArticlePubMedGoogle Scholar
- Muir WM, Wong GK-S, Zhang Y, Wang J, Groenen MAM, Crooijmans RPMA, et al. Genome-wide assessment of worldwide chicken SNP genetic diversity indicates significant absence of rare alleles in commercial breeds. Proc Natl Acad Sci U S A. 2008;105:17312–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Kranis A, Gheyas AA, Boschiero C, Turner F, Yu L, Smith S, et al. Development of a high density 600K SNP genotyping array for chicken. BMC Genomics. 2013;14:59.View ArticlePubMed CentralPubMedGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.View ArticlePubMedGoogle Scholar
- Pollock DL. A geneticist’s perspective from within a broiler primary breeder company. Poult Sci. 1999;78:414–8.View ArticlePubMedGoogle Scholar
- Ibánez-Escriche N, Fernando RL, Toosi A, Dekkers JCM. Genomic selection of purebreds for crossbred performance. Genet Sel Evol. 2009;41:12.View ArticlePubMed CentralPubMedGoogle Scholar
- Bovine HapMap Consortium. Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009;324:528–32.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.